Comparison of assembled Clostridium botulinum A1 genomes revealed their evolutionary relationship

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Discovery of novel virulence mechanisms in Clostridium botulinum type A3 using genome-wide analysis

OBJECTIVE

Clostridium botulinum type A is a neurotoxin-producing, spore-forming anaerobic bacterium that causes botulism in humans. The evolutionary genomic context of this organism is not yet known to understand its molecular virulence mechanisms in the human intestinal tract. Hence, this study aimed to investigate the mechanisms underlying virulence and pathogenesis by comparing the genomic contexts across species, serotypes, and subtypes.

METHODS

A comparative genomic approach was used to analyze evolutionary genomic relationships, intergenomic distances, syntenic blocks, replication origins, and gene abundance with phylogenomic neighbors.

RESULTS

Type A strains have shown genomic proximity to group I strains with distinct accessory genes and vary even within subtypes. Phylogenomic data showed that type C and D strains were distantly related to a group I and group II strains. Synthetic plots indicated that orthologous genes might have evolved from Clostridial ancestry to subtype A3 strains, whereas syntonic out-paralogs might have emerged between subtypes A3 and A1 through α-events. Gene abundance analysis revealed the key roles of genes involved in biofilm formation, cell-cell communication, human diseases, and drug resistance compared to the pathogenic Clostridia. Moreover, we identified 43 unique genes in the type A3 genome, of which 29 were involved in the pathophysiological processes and other genes contributed to amino acid metabolism. The C. botulinum type A3 genome contains 14 new virulence proteins that can provide the ability to confer antibiotic resistance, virulence exertion and adherence to host cells, the host immune system, and mobility of extrachromosomal genetic elements.

CONCLUSION

The results of our study provide insight into the understanding of new virulence mechanisms to discover new therapeutics for the treatment of human diseases caused by type A3 strains.

How Botulinum Neurotoxin Light Chain A1 Maintains Stable Association with the Intracellular Neuronal Plasma Membrane

Botulinum neurotoxin serotype A (BoNT/A) is the most potent protein toxin for humans and is utilized as a therapy for numerous neurologic diseases. BoNT/A comprises a catalytic Light Chain (LC/A) and a Heavy Chain (HC/A) and includes eight subtypes (BoNT/A1-/A8). Previously we showed BoNT/A potency positively correlated with stable localization on the intracellular plasma membrane and identified a low homology domain (amino acids 268-357) responsible for LC/A1 stable co-localization with SNAP-25 on the plasma membrane, while LC/A3 was present in the cytosol of Neuro2A cells. In the present study, steady-state- and live-imaging of a cytosolic LC/A3 derivative (LC/A3V) engineered to contain individual structural elements of the A1 LDH showed that a 59 amino acid region (275-334) termed the MLD was sufficient to direct LC/A3V from the cytosol to the plasma membrane co-localized with SNAP-25. Informatics and experimental validation of the MLD-predicted R1 region (an α-helix, residues 275-300) and R2 region (a loop, α-helix, loop, residues 302-334) both contribute independent steps to the stable co-localization of LC/A1 with SNAP-25 on the plasma membrane of Neuro-2A cells. Understanding how these structural elements contribute to the overall association of LC/A1 on the plasma membrane may identify the molecular basis for the LC contribution of BoNT/A1 to high potency.

Identification of foodborne pathogenic bacteria using confocal Raman microspectroscopy and chemometrics

Rapid and accurate identification of foodborne pathogenic bacteria is of great importance because they are often responsible for the majority of serious foodborne illnesses. The confocal Raman microspectroscopy (CRM) is a fast and easy-to-use method known for its effectiveness in detecting and identifying microorganisms. This study demonstrates that CRM combined with chemometrics can serve as a rapid, reliable, and efficient method for the detection and identification of foodborne pathogenic bacteria without any laborious pre-treatments. Six important foodborne pathogenic bacteria including S. flexneri, L. monocytogenes, V. cholerae, S. aureus, S. typhimurium, and C. botulinum were investigated with CRM. These pathogenic bacteria can be differentiated based on several characteristic peaks and peak intensity ratio. Principal component analysis (PCA) was used for investigating the difference of various samples and reducing the dimensionality of the dataset. Performances of some classical classifiers were compared for bacterial detection and identification including decision tree (DT), artificial neural network (ANN), and Fisher's discriminant analysis (FDA). Correct recognition ratio (CRR), area under the receiver operating characteristic curve (ROC), cumulative gains, and lift charts were used to evaluate the performance of models. The impact of different pretreatment methods on the models was explored, and pretreatment methods include Savitzky-Golay algorithm smoothing (SG), standard normal variate (SNV), multivariate scatter correction (MSC), and Savitzky-Golay algorithm 1st Derivative (SG 1st Der). In the DT, ANN, and FDA model, FDA is more robust for overfitting problem and offers the highest accuracy. Most pretreatment methods raised the performance of the models except SNV. The results revealed that CRM coupled with chemometrics offers a powerful tool for the discrimination of foodborne pathogenic bacteria.

Therapeutic efficacy of new botulinum toxin identified in CCUG 7968 strain

Botulinum neurotoxin type A (BoNT/A) induces muscle atrophy by cleaving synaptosomal-associated protein 25. Thus, BoNT/A has been actively utilized for the treatment of masseter and gastrocnemius hypertrophy. In this study, INI101 toxin was newly identified from the CCUG 7968 strain, and its therapeutic efficacy was evaluated both in vitro and in vivo. The INI101 toxin showed identical genetic sequence, amino acid sequence, and protein subunit composition to BoNT/A produced from strain Hall A. Electromyography (EMG), and immunofluorescence staining demonstrated that INI101 (at 2 ~ 8 U/rat) effectively blocked the neuromuscular junction with no toxicity in a rat model. The EMG results showed INI101 toxin-induced weight loss and volume reduction of the gastrocnemius, similar to the effects of Botox® (BTX). Histological and immunofluorescence staining was consistent with this EMG result, showing that INI101 toxin caused muscle fiber reduction in the gastrocnemius. Notably, INI101 toxin diffused less into adjacent muscle tissue than BTX, indicating that INI101 toxin may reduce potential side effects due to diffusion into normal tissues. INI101 toxin isolated from the novel strain CCUG 7968 is a newly identified meaningful biopharmaceutical comparable to the conventional BoNT/A in the medical field. KEY POINTS: • Botulinum neurotoxin type A (BoNT/A, INI101) was identified from the CCUG 7968 strain. • INI101 toxin showed similar safety and therapeutic efficacy comparable to conventional BoNT/A both in vitro and in vivo. • INI101 toxin is a meaningful biopharmaceutical comparable to the conventional BoNT/A in the medical field.

Comparative in silico genome analysis of Clostridium perfringens unravels stable phylogroups with different genome characteristics and pathogenic potential

Clostridium perfringens causes a plethora of devastating infections, with toxin production being the underlying mechanism of pathogenicity in various hosts. Genomic analyses of 206 public-available C. perfringens strains´ sequence data identified a substantial degree of genomic variability in respect to episome content, chromosome size and mobile elements. However, the position and order of the local collinear blocks on the chromosome showed a considerable degree of preservation. The strains were divided into five stable phylogroups (I–V). Phylogroup I contained human food poisoning strains with chromosomal enterotoxin (cpe) and a Darmbrand strain characterized by a high frequency of mobile elements, a relatively small genome size and a marked loss of chromosomal genes, including loss of genes encoding virulence traits. These features might correspond to the adaptation of these strains to a particular habitat, causing human foodborne illnesses. This contrasts strains that belong to phylogroup II where the genome size points to the acquisition of genetic material. Most strains of phylogroup II have been isolated from enteric lesions in horses and dogs. Phylogroups III, IV and V are heterogeneous groups containing a variety of different strains, with phylogroup III being the most abundant (65.5%). In conclusion, C. perfringens displays five stable phylogroups reflecting different disease involvements, prompting further studies on the evolution of this highly important pathogen.

Botulinum neurotoxin-encoding plasmids can be conjugatively transferred to diverse clostridial strains

Most Group I Clostridium botulinum strains harbor botulinum neurotoxin (bont) genes on their chromosome, while some carry these genes (including bont/a, bont/b, and bont/f) on large plasmids. Prior work in our laboratory demonstrated that Group I BoNT plasmids were mobilized to C. botulinum recipient strains containing the Tn916 transposon. Here, we show that Tn916 is nonessential for plasmid transfer. Relying on an auxotrophic donor phenotype and a plasmid-borne selectable marker, we observed the transfer of pCLJ, a 270 kb plasmid harboring two bont genes, from its host strain to various clostridia. Transfer frequency was greatest to other Group I C. botulinum strains, but the plasmid was also transferred into traditionally nontoxigenic species, namely C. sporogenes and C. butyricum. Expression and toxicity of BoNT/A4 was confirmed in transconjugants by immunoblot and mouse bioassay. These data indicate that conjugation within the genus Clostridium can occur across physiological Groups of C. botulinum, supporting horizontal gene transfer via bont-bearing plasmids. The transfer of plasmids possessing bont genes to resistant Clostridium spp. such as C. sporogenes could impact biological safety for animals and humans. These plasmids may play an environmental role in initiating death in vertebrates, leading to decomposition and nutrient recycling of animal biomass.

The allele combinations of three loci based on, liver, stomach cancers, hematencephalon, COPD and normal population: A preliminary study

The purpose of this study was to examine the specific allele combinations of three loci connected with the liver cancers, stomach cancers, hematencephalon and patients with chronic obstructive pulmonary disease (COPD) and to explore the feasibility of the research methods. We explored different mathematical methods for statistical analyses to assess the association between the genotype and phenotype. At the same time we still analyses the statistical results of allele combinations of three loci by difference value method and ratio method. All the DNA blood samples were collected from patients with 50 liver cancers, 75 stomach cancers, 50 hematencephalon, 72 COPD and 200 normal populations. All the samples were from Chinese. Alleles from short tandem repeat (STR) loci were determined using the STR Profiler plus PCR amplification kit (15 STR loci). Previous research was based on combinations of single-locus alleles, and combinations of cross-loci (two loci) alleles. Allele combinations of three loci were obtained by computer counting and stronger genetic signal was obtained. The methods of allele combinations of three loci can help to identify the statistically significant differences of allele combinations between liver cancers, stomach cancers, patients with hematencephalon, COPD and the normal population. The probability of illness followed different rules and had apparent specificity. This method can be extended to other diseases and provide reference for early clinical diagnosis.

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.

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. In this paper we present modelling based on the integration of diverse information 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. Here, we integrate these regulatory elements into a combined model of population dynamics and gene regulation to build the 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 that could be used to prevent botulism, and potentially contribute to improved methods for the production of toxin used for therapeutics.

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.

Insights into the Origin of Clostridium botulinum Strains: Evolution of Distinct Restriction Endonuclease Sites in rrs (16S rRNA gene)

Diversity analysis of Clostridium botulinum strains is complicated by high microheterogeneity caused by the presence of 9-22 copies of rrs (16S rRNA gene). The need is to mine genetic markers to identify very closely related strains. Multiple alignments of the nucleotide sequences of the 212 rrs of 13 C. botulinum strains revealed intra- and inter-genomic heterogeneity. Low intragenomic heterogeneity in rrs was evident in strains 230613, Alaska E43, Okra, Eklund 17B, Langeland, 657, Kyoto, BKT015925, and Loch Maree. The most heterogenous rrs sequences were those of C. botulinum strains ATCC 19397, Hall, H04402065, and ATCC 3502. In silico restriction mapping of these rrs sequences was observable with 137 type II Restriction endonucleases (REs). Nucleotide changes (NC) at these RE sites resulted in appearance of distinct and additional sites, and loss in certain others. De novo appearances of RE sites due to NC were recorded at different positions in rrs gene. A nucleotide transition A>G in rrs of C. botulinum Loch Maree and 657 resulted in the generation of 4 and 10 distinct RE sites, respectively. Transitions A>G, G>A, and T>C led to the loss of RE sites. A perusal of the entire NC and in silico RE mapping of rrs of C. botulinum strains provided insights into their evolution. Segregation of strains on the basis of RE digestion patterns of rrs was validated by the cladistic analysis involving six house keeping genes: dnaN, gyrB, metG, prfA, pyrG, and Rho.

Identification and characterization of a serious multidrug resistant Stenotrophomonas maltophilia strain in China

An S. maltophilia strain named WJ66 was isolated from a patient; WJ66 showed resistance to more antibiotics than the other S. maltophilia strains. This bacteraemia is resistant to sulphonamides, or fluoroquinolones, while the representative strain of S. maltophilia, K279a, is sensitive to both. To explore drug resistance determinants of this strain, the draft genome sequence of WJ66 was determined and compared to other S. maltophilia sequences. Genome sequencing and genome-wide evolutionary analysis revealed that WJ66 was highly homologous with the strain K279a, but strain WJ66 contained additional antibiotic resistance genes. Further analysis confirmed that strain WJ66 contained an amino acid substitution (Q83L) in fluoroquinolone target GyrA and carried a class 1 integron, with an aadA2 gene in the resistance gene cassette. Homology analysis from the pathogen-host interaction database showed that strain WJ66 lacks raxST and raxA, which is consistent with K279a. Comparative genomic analyses revealed that subtle nucleotide differences contribute to various significant phenotypes in close genetic relationship strains.

Distributions of allele combination in single and cross loci among patients with several kinds of chronic diseases and the normal population

Genetic research has progressed along with scientific and technological developments. However, it is difficult to identify frequency differences in a particular allele distribution at a single locus. Such differences can be identified by examining the allele combination distribution. We explored different mathematical methods for statistical analyses to assess the association between the genotype and phenotype. We investigated the frequency distributions of alleles, combinations of single-locus genes, and combinations of cross-loci genes at 15 loci using 447 blood samples of 200 normal subjects, 72 patients with chronic obstructive pulmonary resistance, 50 liver cancers, 75 stomach cancers and 50 hematencephalon and identified each population as having a unique gene distribution and that the distribution followed certain rules. The probability of illness followed different rules and had apparent specificity. Differences obtained using statistics of combinations of cross-loci genes are superior to single-locus gene statistics, and combinations of single-locus gene statistics are better than allelic statistics.