The genome sequence of Clostridium tetani, the causative agent of tetanus disease

Dissection of the caffeate respiratory chain in the acetogen Acetobacterium woodii: identification of an Rnf-type NADH dehydrogenase as a potential coupling site

The anaerobic acetogenic bacterium Acetobacterium woodii couples caffeate reduction with electrons derived from hydrogen to the synthesis of ATP by a chemiosmotic mechanism with sodium ions as coupling ions, a process referred to as caffeate respiration. We addressed the nature of the hitherto unknown enzymatic activities involved in this process and their cellular localization. Cell extract of A. woodii catalyzes H(2)-dependent caffeate reduction. This reaction is strictly ATP dependent but can be activated also by acetyl coenzyme A (CoA), indicating that there is formation of caffeyl-CoA prior to reduction. Two-dimensional gel electrophoresis revealed proteins present only in caffeate-grown cells. Two proteins were identified by electrospray ionization-mass spectrometry/mass spectrometry, and the encoding genes were cloned. These proteins are very similar to subunits alpha (EtfA) and beta (EtfB) of electron transfer flavoproteins present in various anaerobic bacteria. Western blot analysis demonstrated that they are induced by caffeate and localized in the cytoplasm. Etf proteins are known electron carriers that shuttle electrons from NADH to different acceptors. Indeed, NADH was used as an electron donor for cytosolic caffeate reduction. Since the hydrogenase was soluble and used ferredoxin as an electron acceptor, the missing link was a ferredoxin:NAD(+) oxidoreductase. This activity could be determined and, interestingly, was membrane bound. A search for genes that could encode this activity revealed DNA fragments encoding subunits C and D of a membrane-bound Rnf-type NADH dehydrogenase that is a potential Na(+) pump. These data suggest the following electron transport chain: H(2) --> ferredoxin --> NAD(+) --> Etf --> caffeyl-CoA reductase. They also imply that the sodium motive step in the chain is the ferredoxin-dependent NAD(+) reduction catalyzed by Rnf.

Identification and characterization of the surface-layer protein ofClostridium tetani

Many bacterial species produce a paracrystalline layer, the surface layer, which completely surrounds the exterior of the cell. In some bacteria, the surface layer is implicated in pathogenesis. Two proteins present in cell wall extracts from Clostridium tetani have been investigated and identified one of these has been unambiguously as the surface-layer protein (SLP). The gene, slpA, has been located in the genome of C. tetani E88 that encodes the SLP. The molecular mass of the protein as determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis is considerably larger than that predicted from the gene; however the protein does not appear to be glycosylated. Furthermore, analysis of five C. tetani strains, including three recent clinical isolates, shows considerable variation in the sizes of the SLP.

Genome sequence of Fusobacterium nucleatum subspecies polymorphum - a genetically tractable fusobacterium

Fusobacterium nucleatum is a prominent member of the oral microbiota and is a common cause of human infection. F. nucleatum includes five subspecies: polymorphum, nucleatum, vincentii, fusiforme, and animalis. F. nucleatum subsp. polymorphum ATCC 10953 has been well characterized phenotypically and, in contrast to previously sequenced strains, is amenable to gene transfer. We sequenced and annotated the 2,429,698 bp genome of F. nucleatum subsp. polymorphum ATCC 10953. Plasmid pFN3 from the strain was also sequenced and analyzed. When compared to the other two available fusobacterial genomes (F. nucleatum subsp. nucleatum, and F. nucleatum subsp. vincentii) 627 open reading frames unique to F. nucleatum subsp. polymorphum ATCC 10953 were identified. A large percentage of these mapped within one of 28 regions or islands containing five or more genes. Seventeen percent of the clustered proteins that demonstrated similarity were most similar to proteins from the clostridia, with others being most similar to proteins from other gram-positive organisms such as Bacillus and Streptococcus. A ten kilobase region homologous to the Salmonella typhimurium propanediol utilization locus was identified, as was a prophage and integrated conjugal plasmid. The genome contains five composite ribozyme/transposons, similar to the CdISt IStrons described in Clostridium difficile. IStrons are not present in the other fusobacterial genomes. These findings indicate that F. nucleatum subsp. polymorphum is proficient at horizontal gene transfer and that exchange with the Firmicutes, particularly the Clostridia, is common.

Replication-associated purine asymmetry may contribute to strand-biased gene distribution

Among prokaryotic genomes, the distribution of genes on the leading and lagging strands of the replication fork is known to be biased. Several hypotheses explaining this strand-biased gene distribution (SGD) have been proposed, but none have been tested or supported by sufficient data analyses. In this work we have analyzed 211 prokaryotic genomes in terms of compositional strand asymmetries and the presence or absence of polC and have found that SGD correlates not only with polC, but also with purine asymmetry (PAS). Furthermore, SGD, PAS, and polC are all features associated with a group of low-GC, gram-positive bacteria (Firmicutes). We conclude that PAS is a characteristic of organisms with a heterodimeric DNA polymerase III alpha-subunit constituted by polC and dnaE, which may play a direct role in the maintenance of SGD.

Regulation of sugar catabolism in Lactococcus lactis

The increasing number of genomic and post-genomic studies on Gram-positive organisms and especially on lactic acid bacteria brings a lot of information on sugar catabolism in these bacteria. Like for many other bacteria, glucose is the most preferred source of carbon and energy for Lactococcus lactis. Other carbon sources can induce their own utilization in the absence of well-metabolized sugar. These processes engage numbers of genes and undergo complex mechanisms of regulation. In this review, we discuss various biochemical and genetic control mechanisms involved in sugar catabolism, like regulation by repressors, activators, antiterminators or carbon catabolite repression control.

Regulation of neurotoxin complex expression in Clostridium botulinum strains 62A, Hall A-hyper, and NCTC 2916

The kinetics of botulinum toxin gene expression have been investigated in Clostridium botulinum type A strains 62A, Hall A-hyper, and type A(B) strain NCTC 2916 during the growth cycle. The analyses were performed in TPGY and type A Toxin Production Media (TPM). The mRNA transcript levels encoding the proteins of the neurotoxin complex were determined using Northern analyses. Neurotoxin concentrations in culture supernatants and lysed cell pellets were assayed using ELISA, Western blots, and mouse bioassay. Proteolytic activation of botulinum neurotoxin during the growth cycle was evaluated by Western blots. For all three strains, mRNA transcripts for the toxin complex genes were initially detected in early log phase, reached peak levels in early stationary phase, and rapidly decreased in mid-to-late stationary phase and during lysis. Toxin expression varied depending on the strain and growth medium. Toxin production was highest in strain Hall A-hyper, followed by NCTC 2916 and 62A. For C. botulinum strain Hall A-hyper, cell lysis and toxin release into the supernatant occurred rapidly for cells grown in TPM, while cells grown in TPGY remained in stationary phase with minimal lysis and toxin release through 96 h of growth. In contrast, strains 62A and NCTC 2916 lysed more extensively than Hall A-hyper in TPGY. TPM supported higher toxin production and activation than TPGY in strains 62A and Hall A-hyper. These data support that the genes of the botulinum neurotoxin complex are temporally expressed during late-log and early stationary phase and that toxin complex formation depends on the strain and growth medium. Botulinum toxin synthesis and activation appears to be a complex process that is highly regulated by nutritional and environmental conditions. Further research is needed to elucidate the sensing mechanisms and genetic regulatory factors controlling these processes.

Organization and regulation of the neurotoxin genes in Clostridium botulinum and Clostridium tetani

Botulinum and tetanus neurotoxins are structurally and functionally related 150 kDa proteins that are potent inhibitors of neuroexocytosis. Botulinum neurotoxin associates with non-toxic proteins to form complexes of various sizes. The botulinum neurotoxin and non-toxic protein genes are clustered in a DNA segment called the botulinum locus. This locus is probably located on a mobile or degenerate mobile element, which accounts for the various genomic localizations (chromosome, plasmid, phage) in different Clostridium botulinum types. The botulinum neurotoxin and non-toxic protein genes are organized in two polycistronic operons (ntnh-bont and ha operons) transcribed in opposite orientations. The gene that separates the two operons of the botulinum locus in C. botulinum A encodes a 21 kDa protein BotR/A, which is a positive regulator of the expression of the botulinum locus genes. Similarly, in Clostridium tetani, the gene located immediately upstream of the tetanus toxin gene, encodes a positive regulatory protein, TetR. BotR and TetR are possibly alternative sigma factors related to TxeR and UviA, which regulate C. difficile toxin and C. perfringens bacteriocin production, respectively. TxeR and UviA define a new sub-group of the sigma(70) family of RNA polymerase initiation factors. In addition, the C. botulinum genome contains predicted two-component system genes, some of which are possibly involved in regulation of toxinogenesis.

Genome sequence of a proteolytic (Group I) Clostridium botulinum strain Hall A and comparative analysis of the clostridial genomes

Clostridium botulinum is a heterogeneous Gram-positive species that comprises four genetically and physiologically distinct groups of bacteria that share the ability to produce botulinum neurotoxin, the most poisonous toxin known to man, and the causative agent of botulism, a severe disease of humans and animals. We report here the complete genome sequence of a representative of Group I (proteolytic) C. botulinum (strain Hall A, ATCC 3502). The genome consists of a chromosome (3,886,916 bp) and a plasmid (16,344 bp), which carry 3650 and 19 predicted genes, respectively. Consistent with the proteolytic phenotype of this strain, the genome harbors a large number of genes encoding secreted proteases and enzymes involved in uptake and metabolism of amino acids. The genome also reveals a hitherto unknown ability of C. botulinum to degrade chitin. There is a significant lack of recently acquired DNA, indicating a stable genomic content, in strong contrast to the fluid genome of Clostridium difficile, which can form longer-term relationships with its host. Overall, the genome indicates that C. botulinum is adapted to a saprophytic lifestyle both in soil and aquatic environments. This pathogen relies on its toxin to rapidly kill a wide range of prey species, and to gain access to nutrient sources, it releases a large number of extracellular enzymes to soften and destroy rotting or decayed tissues.

Polyketide synthases and nonribosomal peptide synthetases: the emerging view from bacterial genomics

A total of 223 complete bacterial genomes are analyzed, with 281 citations, for the presence of genes encoding modular polyketide synthases (PKS) and nonribosomal peptide synthetases (NRPS). We report on the distribution of these systems in different bacterial taxa and, whenever known, the metabolites they synthesize. We also highlight, in the different bacterial lineages, the PKS and NRPS genes and, whenever known, the corresponding products.

The genome of Syntrophus aciditrophicus: life at the thermodynamic limit of microbial growth

Biochemically, the syntrophic bacteria constitute the missing link in our understanding of anaerobic flow of carbon in the biosphere. The completed genome sequence of Syntrophus aciditrophicus SB, a model fatty acid- and aromatic acid-degrading syntrophic bacterium, provides a glimpse of the composition and architecture of the electron transfer and energy-transducing systems needed to exist on marginal energy economies of a syntrophic lifestyle. The genome contains 3,179,300 base pairs and 3,169 genes where 1,618 genes were assigned putative functions. Metabolic reconstruction of the gene inventory revealed that most biosynthetic pathways of a typical Gram-negative microbe were present. A distinctive feature of syntrophic metabolism is the need for reverse electron transport; the presence of a unique Rnf-type ion-translocating electron transfer complex, menaquinone, and membrane-bound Fe-S proteins with associated heterodisulfide reductase domains suggests mechanisms to accomplish this task. Previously undescribed approaches to degrade fatty and aromatic acids, including multiple AMP-forming CoA ligases and acyl-CoA synthetases seem to be present as ways to form and dissipate ion gradients by using a sodium-based energy strategy. Thus, S. aciditrophicus, although nutritionally self-sufficient, seems to be a syntrophic specialist with limited fermentative and respiratory metabolism. Genomic analysis confirms the S. aciditrophicus metabolic and regulatory commitment to a nonconventional mode of life compared with our prevailing understanding of microbiology.

Mutations in fusA associated with posttherapy fusidic acid resistance in Clostridium difficile

In silico, we identified fusA (2,067 bp) in Clostridium difficile 630. Sequencing of fusA in posttherapy fusidic acid-resistant C. difficile isolates from 12 patients with C. difficile-associated diarrhea (CDAD) identified fusA mutations, one or two nonsynonymous substitutions, or in one case a deletion of one codon associated with resistance. Five of these mutations have previously been described in fusA of fusidic acid-resistant Staphylococcus aureus, but seven were novel fusA mutations. Fusidic acid monotherapy for CDAD seemed to rapidly select conserved resistant mutants.

Characterization of an exceedingly active NADH oxidase from the anaerobic hyperthermophilic bacterium Thermotoga maritima

An NADH oxidase from the anaerobic hyperthermophilic bacterium Thermotoga maritima was purified. The enzyme was very active in catalyzing the reduction of oxygen to hydrogen peroxide with an optimal pH value of 7 at 80 degrees C. The V(max) was 230 +/- 14 mumol/min/mg (k(cat)/K(m) = 548,000 min(-1) mM(-1)), and the K(m) values for NADH and oxygen were 42 +/- 3 and 43 +/- 4 muM, respectively. The NADH oxidase was a heterodimeric flavoprotein with two subunits with molecular masses of 54 kDa and 46 kDa. Its gene sequences were identified, and the enzyme might represent a new type of NADH oxidase in anaerobes. An NADH-dependent peroxidase with a specific activity of 0.1 U/mg was also present in the cell extract of T. maritima.

Organization of butyrate synthetic genes in human colonic bacteria: phylogenetic conservation and horizontal gene transfer

Butyrate producers constitute an important bacterial group in the human large intestine. Butyryl-CoA is formed from two molecules of acetyl-CoA in a process resembling beta-oxidation in reverse. Three different arrangements of the six genes coding for this pathway have been found in low mol% G+C-content gram-positive human colonic bacteria using DNA sequencing and degenerate PCR. Gene arrangements were strongly conserved within phylogenetic groups defined by 16S rRNA gene sequence relationships. In the case of one of the genes, encoding beta-hydroxybutyryl-CoA dehydrogenase, however, sequence relationships were strongly suggestive of horizontal gene transfer between lineages. The newly identified gene for butyryl-CoA CoA-transferase, which performs the final step in butyrate formation in most known human colonic bacteria, was not closely linked to these central pathway genes.

The genome and transcriptomes of the anti-tumor agent Clostridium novyi-NT

Bacteriolytic anti-cancer therapies employ attenuated bacterial strains that selectively proliferate within tumors. Clostridium novyi-NT spores represent one of the most promising of these agents, as they generate potent anti-tumor effects in experimental animals. We have determined the 2.55-Mb genomic sequence of C. novyi-NT, identifying a new type of transposition and 139 genes that do not have homologs in other bacteria. The genomic sequence was used to facilitate the detection of transcripts expressed at various stages of the life cycle of this bacterium in vitro as well as in infections of tumors in vivo. Through this analysis, we found that C. novyi-NT spores contained mRNA and that the spore transcripts were distinct from those in vegetative forms of the bacterium.

Complementation of a Clostridium perfringens spo0A mutant with wild-type spo0A from other Clostridium species

To evaluate whether C. perfringens can be used as a model organism for studying the sporulation process in other clostridia, C. perfringens spo0A mutant IH101 was complemented with wild-type spo0A from four different Clostridium species. Wild-type spo0A from C. acetobutylicum or C. tetani, but not from C. botulinum or C. difficile, restored sporulation and enterotoxin production in IH101. The ability of spo0A from C. botulinum or C. difficile to complement the lack of spore formation in IH101 might be due, at least in part, to the low levels of spo0A transcription and Spo0A production.

The mosaic nature of intergenic 16S-23S rRNA spacer regions suggests rRNA operon copy number variation in Clostridium difficile strains

Clostridium difficile is a major spore-forming environmental pathogen that causes serious health problems in patients undergoing antibiotic therapy. Consequently, reliable and sensitive methods for typing individual strains are required for epidemiological and environmental studies. Ribotyping is generally considered the best method, but it fails to account for sequence diversity which might exist in intergenic 16S-23S rRNA spacer regions (ISRs) within and among strains of this organism. Therefore, this study was undertaken to compare the sequence of each individual ISR in five strains of C. difficile to explore the extent of this diversity and see whether such information might provide the basis for more sensitive and discriminatory strain typing methods. After targeted PCR amplification, cloning, and sequencing, the diversity of the ISRs was used as a measure of rRNA operon copy number. In C. difficile strains 630, ATCC 43593, A, and B, 11, 11, 7, and 8 ISR length variants, respectively, were found (containing different combinations of sequence groups [i to xiii]), suggesting 11, 11, 7, and 8 rrn copies in the respective strains. Many ISRs of the same length differed markedly in their sequences, and some of these were restricted in occurrence to a single strain. Most of these ISRs did not contain any tRNA genes, and only single copies of the tRNA(Ala) gene were found in those that did. The presence of ISR sequence groups (i to xiii) varied between strains, with some found in one, two, three, four, or all five strains. We conclude that the intergenic 16S-23S rRNA spacer regions showed a high degree of diversity, not only among the rrn operons in different strains and different rrn copies in a single strain but also among ISRs of the same length. It appears that C. difficile ISRs vary more at the inter- and intragenic levels than those of other species as determined by empirical comparison of sequences. The precise characterization of these sequences has demonstrated a high level of mosaic sequence block rearrangements that are present or absent in multiple strain-variable rrn copies within and between five different strains of C. difficile.

Structures of apo- and holo-tyrosine phenol-lyase reveal a catalytically critical closed conformation and suggest a mechanism for activation by K+ ions

Tyrosine phenol-lyase, a tetrameric pyridoxal 5'-phosphate dependent enzyme, catalyzes the reversible hydrolytic cleavage of L-tyrosine to phenol and ammonium pyruvate. Here we describe the crystal structure of the Citrobacter freundii holoenzyme at 1.9 A resolution. The structure reveals a network of protein interactions with the cofactor, pyridoxal 5'-phosphate, and details of coordination of the catalytically important K+ ion. We also present the structure of the apoenzyme at 1.85 A resolution. Both structures were determined using crystals grown at pH 8.0, which is close to the pH of the maximal enzymatic activity (8.2). Comparison of the apoenzyme structure with the one previously determined at pH 6.0 reveals significant differences. The data suggest that the decrease of the enzymatic activity at pH 6.0 may be caused by conformational changes in the active site residues Tyr71, Tyr291, and Arg381 and in the monovalent cation binding residue Glu69. Moreover, at pH 8.0 we observe two different active site conformations: open, which was characterized before, and closed, which is observed for the first time in beta-eliminating lyases. In the closed conformation a significant part of the small domain undergoes an extraordinary motion of up to 12 A toward the large domain, closing the active site cleft and bringing the catalytically important Arg381 and Phe448 into the active site. The closed conformation allows rationalization of the results of previous mutational studies and suggests that the observed active site closure is critical for the course of the enzymatic reaction and for the enzyme's specificity toward its physiological substrate. Finally, the closed conformation allows us to model keto(imino)quinonoid, the key transition intermediate.

Skewed genomic variability in strains of the toxigenic bacterial pathogen, Clostridium perfringens

Clostridium perfringens is a Gram-positive, anaerobic spore-forming bacterium commonly found in soil, sediments, and the human gastrointestinal tract. C. perfringens is responsible for a wide spectrum of disease, including food poisoning, gas gangrene (clostridial myonecrosis), enteritis necroticans, and non-foodborne gastrointestinal infections. The complete genome sequences of Clostridium perfringens strain ATCC 13124, a gas gangrene isolate and the species type strain, and the enterotoxin-producing food poisoning strain SM101, were determined and compared with the published C. perfringens strain 13 genome. Comparison of the three genomes revealed considerable genomic diversity with >300 unique "genomic islands" identified, with the majority of these islands unusually clustered on one replichore. PCR-based analysis indicated that the large genomic islands are widely variable across a large collection of C. perfringens strains. These islands encode genes that correlate to differences in virulence and phenotypic characteristics of these strains. Significant differences between the strains include numerous novel mobile elements and genes encoding metabolic capabilities, strain-specific extracellular polysaccharide capsule, sporulation factors, toxins, and other secreted enzymes, providing substantial insight into this medically important bacterial pathogen.

Sodium ion pumps and hydrogen production in glutamate fermenting anaerobic bacteria

Anaerobic bacteria ferment glutamate via two different pathways to ammonia, carbon dioxide, acetate, butyrate and molecular hydrogen. The coenzyme B12-dependent pathway in Clostridium tetanomorphum via 3-methylaspartate involves pyruvate:ferredoxin oxidoreductase and a novel enzyme, a membrane-bound NADH:ferredoxin oxidoreductase. The flavin- and iron-sulfur-containing enzyme probably uses the energy difference between reduced ferredoxin and NADH to generate an electrochemical Na+ gradient, which drives transport processes. The other pathway via 2-hydroxyglutarate in Acidaminococcus fermentans and Fusobacterium nucleatum involves glutaconyl-CoA decarboxylase, which uses the free energy of decarboxylation to generate also an electrochemical Na+ gradient. In the latter two organisms, similar membrane-bound NADH:ferredoxin oxidoreductases have been characterized. We propose that in the hydroxyglutarate pathway these oxidoreductases work in the reverse direction, whereby the reduction of ferredoxin by NADH is driven by the Na+ gradient. The reduced ferredoxin is required for hydrogen production and the activation of radical enzymes. Further examples show that reduced ferredoxin is an agent, whose reducing energy is about 1 ATP 'richer' than that of NADH.

The multidrug-resistant human pathogen Clostridium difficile has a highly mobile, mosaic genome

We determined the complete genome sequence of Clostridium difficile strain 630, a virulent and multidrug-resistant strain. Our analysis indicates that a large proportion (11%) of the genome consists of mobile genetic elements, mainly in the form of conjugative transposons. These mobile elements are putatively responsible for the acquisition by C. difficile of an extensive array of genes involved in antimicrobial resistance, virulence, host interaction and the production of surface structures. The metabolic capabilities encoded in the genome show multiple adaptations for survival and growth within the gut environment. The extreme genome variability was confirmed by whole-genome microarray analysis; it may reflect the organism's niche in the gut and should provide information on the evolution of virulence in this organism.

Botulinum neurotoxin types A, B, and E: fragmentations by autoproteolysis and other mechanisms including by O-phenanthroline-dithiothreitol, and association of the dinucleotides NAD(+)/NADH with the heavy chain of the three neurotoxins

The first evidence of autoproteolytic activity of the approximately 50-kDa light chain of the clostridial neurotoxins (NT) is traceable to the observations that the light chains of botulinum NT serotypes A and E, separated from their approximately 100-kDa heavy chain conjugate, were found cleaved at the amino side of Tyr250 and Arg244, respectively [DasGupta and Foley (1989). Biochimie 71: 1183-1200]. Specific cleavages of the recombinant light chain of NT type A, including at Tyr249-Tyr250, firmly established that the cleavages reported earlier were due to autoproteolysis [Ahmed et al. (2001). J. Protein Chem. 20: 221-231; Ahmed et al. (2003). Biochemistry 42:12539-12549] and not by contaminating proteases or non-enzymatic. We now report many cleavages in the NT types A, B and E and also in their separated light and heavy chains, and identification of several of the peptide bonds cleaved. None of the identified cleaved bonds (-P1-P1' -) in one serotype (except Asp-Pro) was found common in other serotypes or cleaved within itself at a second site. After separation from the heavy chain self-cleavages of the light chains of type A, B and E at Tyr249-Tyr250, Gln258-Ser259 and Ile243-Arg244, respectively indicate an intriguing feature (in the aligned sequences these bonds of type A and B are 2 and type A and E are 4 peptide bonds apart) that may have some role in the NT's structure-function relationship yet to be understood. We point out that autoproteolysis of a single peptide bond (Phe418-Thr419 or Phe422-Glu423) in NT type A reported by Ahmed et al. (2001) can potentially generate proteolytically active light chain freed of the heavy chain; this is an efficient pathway, that by-passes nicking by a trypsin-like protease(s) inside the intrachain disulfide bridge and its reductive cleavage. We offer probable explanations for the observed cleavages such as acid- and metal-mediated (non-catalytic and non-stoichiometric) reactions in addition to autoproteolysis but cannot predict which mechanism(s) of cleavage occur or prevail following NT's entry in the body as poison or therapeutic agent. The metal chelator O-phenanthroline (above critical miceller concentration) in the presence of dithiothreitol cleaved type E NT at limited sites generating discrete 114-, 87-, 49-, 42-, and 31-kDa fragments but degraded NTs type A and B extensively. The limited cleavage of type E NT was dependent on the presence of metal ion(s) bound to the protein and its native (urea sensitive) conformation. The self-cleavage of the NTs at specific sites prompted us to search for specific binding sites on the NTs analogous to SNARE-motifs-the 9-residuelong motifs present on the NT's natural substrates (SNAP-25, syntaxin, VAMP/synaptobrevin); such putative binding motifs (sites) noted on all clostridial NTs are reported here. Their relationship to the observed autoproteolysis remains to be determined experimentally. The dinucleotide NAD(+)/NADH associated with the NTs type A, B and E (2-3 NADH per protein molecule) via their H-chains, and a portion of the H-chain (toward the C-terminus) appears to exhibit limited amino acid sequence homology with lactate dehydrogenase-a representative NAD(+)/NADH binding protein.

Predicting Shine-Dalgarno sequence locations exposes genome annotation errors

In prokaryotes, Shine-Dalgarno (SD) sequences, nucleotides upstream from start codons on messenger RNAs (mRNAs) that are complementary to ribosomal RNA (rRNA), facilitate the initiation of protein synthesis. The location of SD sequences relative to start codons and the stability of the hybridization between the mRNA and the rRNA correlate with the rate of synthesis. Thus, accurate characterization of SD sequences enhances our understanding of how an organism's transcriptome relates to its cellular proteome. We implemented the Individual Nearest Neighbor Hydrogen Bond model for oligo-oligo hybridization and created a new metric, relative spacing (RS), to identify both the location and the hybridization potential of SD sequences by simulating the binding between mRNAs and single-stranded 16S rRNA 3' tails. In 18 prokaryote genomes, we identified 2,420 genes out of 58,550 where the strongest binding in the translation initiation region included the start codon, deviating from the expected location for the SD sequence of five to ten bases upstream. We designated these as RS+1 genes. Additional analysis uncovered an unusual bias of the start codon in that the majority of the RS+1 genes used GUG, not AUG. Furthermore, of the 624 RS+1 genes whose SD sequence was associated with a free energy release of less than -8.4 kcal/mol (strong RS+1 genes), 384 were within 12 nucleotides upstream of in-frame initiation codons. The most likely explanation for the unexpected location of the SD sequence for these 384 genes is mis-annotation of the start codon. In this way, the new RS metric provides an improved method for gene sequence annotation. The remaining strong RS+1 genes appear to have their SD sequences in an unexpected location that includes the start codon. Thus, our RS metric provides a new way to explore the role of rRNA-mRNA nucleotide hybridization in translation initiation.

Electron transport in the pathway of acetate conversion to methane in the marine archaeon Methanosarcina acetivorans

A liquid chromatography-hybrid linear ion trap-Fourier transform ion cyclotron resonance mass spectrometry approach was used to determine the differential abundance of proteins in acetate-grown cells compared to that of proteins in methanol-grown cells of the marine isolate Methanosarcina acetivorans metabolically labeled with 14N versus 15N. The 246 differentially abundant proteins in M. acetivorans were compared with the previously reported 240 differentially expressed genes of the freshwater isolate Methanosarcina mazei determined by transcriptional profiling of acetate-grown cells compared to methanol-grown cells. Profound differences were revealed for proteins involved in electron transport and energy conservation. Compared to methanol-grown cells, acetate-grown M. acetivorans synthesized greater amounts of subunits encoded in an eight-gene transcriptional unit homologous to operons encoding the ion-translocating Rnf electron transport complex previously characterized from the Bacteria domain. Combined with sequence and physiological analyses, these results suggest that M. acetivorans replaces the H2-evolving Ech hydrogenase complex of freshwater Methanosarcina species with the Rnf complex, which generates a transmembrane ion gradient for ATP synthesis. Compared to methanol-grown cells, acetate-grown M. acetivorans synthesized a greater abundance of proteins encoded in a seven-gene transcriptional unit annotated for the Mrp complex previously reported to function as a sodium/proton antiporter in the Bacteria domain. The differences reported here between M. acetivorans and M. mazei can be attributed to an adaptation of M. acetivorans to the marine environment.

Characterization of a unique type IA topoisomerase in Bacillus cereus

Bacillus cereus topoisomerase IIIbeta (bcTopo IIIbeta) has been cloned, overexpressed and biochemically characterized. This enzyme exhibits 64% and 33% sequence identity to Bacillus subtilis topoisomerase III (bsTopo III) and Escherichia coli topoisomerase III (ecTopo III) respectively. The enzymatic properties of bcTopo IIIbeta differ substantially from other bacterial type IA topoisomerases, including E. coli type IA topoisomerases and B. cereus topoisomerase I (bcTopo I) and IIIalpha (bcTopo IIIalpha). bcTopo IIIbeta only partially relaxes negatively supercoiled DNA and appears incapable of generating fully relaxed topoisomers. In contrast to ecTopo III and bcTopo IIIalpha, bcTopo IIIbeta is not a decatenase. bcTopo IIIbeta is unable to compensate the loss of ecTopo III in vivo. Therefore, bcTopo IIIbeta is a unique prokaryotic type IA topoisomerase that is different from previously characterized topoisomerases.

Origin of replication in circular prokaryotic chromosomes

To predict origins of replication in prokaryotic chromosomes, we analyse the leading and lagging strands of 200 chromosomes for differences in oligomer composition and show that these correlate strongly with taxonomic grouping, lifestyle and molecular details of the replication process. While all bacteria have a preference for Gs over Cs on the leading strand, we discover that the direction of the A/T skew is determined by the polymerase-alpha subunit that replicates the leading strand. The strength of the strand bias varies greatly between both phyla and environments and appears to correlate with growth rate. Finally we observe much greater diversity of skew among archaea than among bacteria. We have developed a program that accurately locates the origins of replication by measuring the differences between leading and lagging strand of all oligonucleotides up to 8 bp in length. The program and results for all publicly available genomes are available from http://www.cbs.dtu.dk/services/GenomeAtlas/suppl/origin.

A comparative genomic view of clostridial sporulation and physiology

Clostridia are anaerobic, endospore-forming prokaryotes that include strains of importance to human and animal health and physiology, cellulose degradation, solvent production and bioremediation. Their differentiation and related developmental programmes are not well understood at the molecular level. Recent genome sequencing and transcriptional-profiling studies have offered a glimpse of their inner workings and indicate that a better understanding of the orchestration of the molecular events that underlie their unique physiology, capabilities and diversity will pay major dividends.

The impact of genomics on vaccine design

After 200 years of practice, vaccinology has gained new perspectives for preventing infectious diseases. Sequencing of complete bacterial genomes led to the development of new large-scale technologies, such as bioinformatics, proteomics and DNA microarrays. By examining genetic content, as well as transcription and expression profiles, a more detailed understanding of bacterial pathogenesis can be reached. Moreover, the whole-genome perspective is expected to provide an instrumental contribution to vaccine development, particularly to target those pathogens for which the traditional approaches have failed so far. In this review, we describe how genomic approaches can be used to identify novel vaccine candidates or create safer live-attenuated vaccines.

Evolutionary change - patterns and processes

The present review considered: (a) the factors that conditioned the early transition from non-life to life; (b) genome structure and complexity in prokaryotes, eukaryotes, and organelles; (c) comparative human chromosome genomics; and (d) the Brazilian contribution to some of these studies. Understanding the dialectical conflict between freedom and organization is fundamental to give meaning to the patterns and processes of organic evolution.

Demonstration of cotranscription and 1-methyl-3-nitroso-nitroguanidine induction of a 30-gene operon of Borrelia burgdorferi: evidence that the 32-kilobase circular plasmids are prophages

The Borrelia genome is comprised of linear and circular elements, including a group of 32-kb circular plasmids (cp32s). Earlier analyses identified a bacteriophage, varphiBB-1, that may package cp32s, suggesting that these plasmids are prophages. cp32-8, cp32-9, and cp32-1 (plasmids L, N, and P, respectively) encode virulence factors such as the factor H binding, OspE proteins (BBL39, BBN38, and BBP38). Here the expression patterns of cp32-8 open reading frames (ORFs) in in vitro-cultivated 1-methyl-3-nitroso-nitroguanidine (MNNG)-treated and untreated spirochetes and during infection were assessed. ORFs BBL42 through BBL28, which encode several bacteriophage protein homologs, were found to be cotranscribed and expression was upregulated by MNNG. Immunoblotting revealed that MNNG-induced transcription led to increased protein production. The expression of several genes that reside outside of the BBL42-BBL28 operon was not affected by MNNG. Some of these genes, including OspE (BBL39), appear to represent morons. Real-time reverse transcription-PCR of spirochetes in mouse tissue revealed that although the phage operon was not induced during infection, transcription of BBL23 (previously designated BlyA), a putative holin, was upregulated. This observation indicates that some genes within the operon can be independently transcribed from internal promoters. Additional transcriptional analyses of the operon identified multiple transcriptional start sites and provided evidence for the expression of a homologous operon from other cp32s. The data support the hypothesis put forth by C. Eggers and D. S. Samuels (J. Bacteriol. 181:7308-7313, 1999) that the cp32s are prophages, a finding with broad implications for our understanding of Borrelia pathogenesis and Borrelia genome evolution.

The genome sequence of Clostridium botulinum type C neurotoxin-converting phage and the molecular mechanisms of unstable lysogeny

Botulinum neurotoxins (BoNTXs) produced by Clostridium botulinum are among the most poisonous substances known. Of the seven types of BoNTXs, genes for type C1 and D toxins (BoNTX/C1 and D) are carried by bacteriophages. The gene for exoenzyme C3 also resides on these phages. Here, we present the complete genome sequence of c-st, a representative of BoNTX/C1-converting phages. The genome is a linear double-stranded DNA of 185,682 bp with 404-bp terminal direct repeats, the largest known temperate phage genome. We identified 198 potential protein-coding regions, including the genes for production of BoNTX/C1 and exoenzyme C3. Very exceptionally, as a viable bacteriophage, a number of insertion sequences were found on the c-st genome. By analyzing the molecular structure of the c-st genome in lysogens, we also found that it exists as a circular plasmid prophage. These features account for the unstable lysogeny of BoNTX phages, which has historically been called "pseudolysogeny." The PCR scanning analysis of other BoNTX/C1 and D phages based on the c-st sequence further revealed that BoNTX phages comprise a divergent phage family, probably generated by exchanging genomic segments among BoNTX phages and their relatives.

Molecular characterization of Clostridium tetani strains by pulsed-field gel electrophoresis and colony PCR

Pulsed-field gel electrophoresis and PCR were applied for the first time to the molecular characterization of Clostridium tetani. Among five strains tested, one (CN1339) turned out to contain a mixture of two genetically different clones and two (D11 and G761) to contain bacteria differing by the presence or absence of the 74-kb plasmid harboring the tetX gene.

Improvement in laboratory diagnosis of wound botulism and tetanus among injecting illicit-drug users by use of real-time PCR assays for neurotoxin gene fragments

An upsurge in wound infections due to Clostridium botulinum and Clostridium tetani among users of illegal injected drugs (IDUs) occurred in the United Kingdom during 2003 and 2004. A real-time PCR assay was developed to detect a fragment of the neurotoxin gene of C. tetani (TeNT) and was used in conjunction with previously described assays for C. botulinum neurotoxin types A, B, and E (BoNTA, -B, and -E). The assays were sensitive, specific, rapid to perform, and applicable to investigating infections among IDUs using DNA extracted directly from wound tissue, as well as bacteria growing among mixed microflora in enrichment cultures and in pure culture on solid media. A combination of bioassay and PCR test results confirmed the clinical diagnosis in 10 of 25 cases of suspected botulism and two of five suspected cases of tetanus among IDUs. The PCR assays were in almost complete agreement with the conventional bioassays when considering results from different samples collected from the same patient. The replacement of bioassays by real-time PCR for the isolation and identification of both C. botulinum and C. tetani demonstrates a sensitivity and specificity similar to those of conventional approaches. However, the real-time PCR assays substantially improves the diagnostic process in terms of the speed of results and by the replacement of experimental animals. Recommendations are given for an improved strategy for the laboratory investigation of suspected wound botulism and tetanus among IDUs.

The Clostridium botulinum GerAB germination protein is located in the inner membrane of spores

Clostridium botulinum dormant spores germinate in presence of l-alanine via a specific receptor composed of GerAA, GerAB and GerAC proteins. In Bacillus subtilis spores, GerAA and GerAC proteins were located in the inner membrane of the spore. We studied the location of the GerAB protein in C. botulinum spore fractions by Western-blot analysis, using an antipeptidic antibody. The protein GerAB was in vitro translated and used to confirm the specificity of the antibodies. GerAB was not present in a coat and spore outer membrane fraction but was present in a fraction of decoated spores containing inner membrane. These results strongly suggest that the protein GerAB is located in the inner membrane of the spore.

Genomics of clostridial pathogens: implication of extrachromosomal elements in pathogenicity

The recently decoded genomes of the major clostridial toxin-producing pathogens Clostridium perfringens, Clostridium tetani, Clostridium botulinum and Clostridium difficile have provided a huge amount of new sequence data. Recent studies have focused on the identification and investigation of pathogenic determinants and the regulatory events governing their expression. The sequence data revealed also the genomic background of virulence genes, as well as the contribution of extrachromosomal elements to a pathogenic phenotype. This has generated new insights in clostridial pathogenesis - and will continue to do so in the future - and has deepened our understanding of the anaerobic lifestyle of clostridial species.

Radical-mediated dehydration reactions in anaerobic bacteria

Most dehydratases catalyse the elimination of water from β-hydroxy ketones, β-hydroxy carboxylic acids or β-hydroxyacyl-CoA. The electron-withdrawing carbonyl functionalities acidify the α-hydrogens to enable their removal by basic amino acid side chains. Anaerobic bacteria, however, ferment amino acids via α- or γ-hydroxyacyl-CoA, dehydrations of which involve the abstraction of a β-hydrogen, which is ostensibly non-acidic (pKca. 40). Evidence is accumulating that β-hydrogens are acidified via transient conversion of the CoA derivatives to enoxy radicals by one-electron transfers, which decrease the pKto 14. The dehydrations of (R)-2-hydroxyacyl-CoA to (E)-2-enoyl-CoA are catalysed by heterodimeric [4Fe-4S]-containing dehydratases, which require reductive activation by an ATP-dependent one-electron transfer mediated by a homodimeric protein with a [4Fe-4S] cluster between the two subunits. The electron is further transferred to the substrate, yielding a ketyl radical anion, which expels the hydroxyl group and forms an enoxy radical. The dehydration of 4-hydroxybutyryl-CoA to crotonyl-CoA involves a similar mechanism, in which the ketyl radical anion is generated by one-electron oxidation. The structure of the FAD- and [4Fe-4S]-containing homotetrameric dehydratase is related to that of acyl-CoA dehydrogenases, suggesting a radical-based mechanism for both flavoproteins.

Hypothesis: a motor neuron toxin produced by a clostridial species residing in gut causes ALS

We hypothesize that a yet-to-be-identified motor neuron toxin produced by a clostridial species causes sporadic amyotrophic lateral sclerosis (ALS) in susceptible individuals. This clostridial species would reside undetected in the gut and chronically produce a toxin that targets the motor system, like the tetanus and botulinum toxins. After gaining access to the lower motor neuron, the toxin would be transported back to the cell body, as occurs with the tetanus toxin, and destroy the lower motor neuron - the essential feature of ALS. Again like the tetanus toxin, some of the toxin would cross to neighboring cells and to the upper motor neuron and similarly destroy these motor neurons. Weakness would relentlessly progress until not enough motor neurons remained to sustain life. If this hypothesis were correct, treatment with appropriate antibiotics or antitoxins might slow or halt progression of disease, and immunization might prevent disease.

Bacillus cereus DNA topoisomerase I and IIIalpha: purification, characterization and complementation of Escherichia coli TopoIII activity

The Bacillus cereus genome possesses three type IA topoisomerase genes. These genes, encoding DNA topoisomerase I and IIIα (bcTopo I, bcTopo IIIα), have been cloned into T7 RNA polymerase-regulated plasmid expression vectors and the enzymes have been overexpressed, purified and characterized. The proteins exhibit similar biochemical activity to their Escherichia coli counterparts, DNA topoisomerase I and III (ecTopo I, ecTopo III). bcTopo I is capable of efficiently relaxing negatively supercoiled DNA in the presence of Mg²⁺ but does not possess an efficient DNA decatenation activity. bcTopo IIIα is an active topoisomerase that is capable of relaxing supercoiled DNA at a broad range of Mg²⁺ concentrations; however, its DNA relaxation activity is not as efficient as that of bcTopo I. In addition, bcTopo III is a potent DNA decatenase that resolves oriC-based plasmid replication intermediates in vitro. Interestingly, bcTopo I and bcTopo IIIα are both able to compensate for the loss of ecTopo III in E.coli cells that lack ecTopo I. In contrast, ecTopo I cannot substitute for ecTopo III under these conditions.

Essential cysteine residues in Bacillus subtilis spore photoproduct lyase identified by alanine scanning mutagenesis

Endospore-forming bacteria (Bacillus and Clostridium spp.) are highly ultraviolet (UV) resistant and repair UV-induced DNA damage in part using the spore-specific DNA repair enzyme spore photoproduct (SP) lyase. SP lyase in all known sporeformers contains four conserved cysteine residues; three absolutely conserved residues are located at the "Radical SAM" consensus (C91xxxC95xxC98), which presumably participates in [4Fe-4S] cluster formation. A fourth conserved cysteine, the function of which is unknown, is located at C141 in SP lyase from all Bacillus spp. sequenced to date. To probe the function of the fourth cysteine, each conserved cysteine in the B. subtilis SP lyase was systematically altered to alanine by site-directed mutagenesis. UV-visible spectroscopy of wild-type and mutant SP lyases indicated that C141 does not participate in [4Fe-4S] formation and redox chemistry; however, in vivo SP lyase activity was abolished in all mutants, indicating an essential role for C141 in SP lyase activity.

Genomic analysis of the protein secretion systems in Clostridium acetobutylicum ATCC 824

Consistent information about protein secretion in Gram-positive bacteria is essentially restricted to the model organism Bacillus subtilis. Among genome-sequenced clostridia, Clostridium acetobutylicum has been the most extensively studied from a physiological point of view and is the organism for which the largest variety of molecular biology tools have been developed. Following in silico analyses, both secreted proteins and protein secretion systems were identified. The Tat (Twin arginine translocation; TC #2.A.64) pathway and ABC (ATP binding cassette) protein exporters (TC #3.A.1.) could not be identified, but the Sec (secretion) pathway (TC #3.A.5) appears to be used prevalently. Similarly, a flagella export apparatus (FEA; TC #3.A.6.), holins (TC #1.E.), and an ESAT-6/WXG100 (early secreted antigen target of 6 kDa/proteins with a WXG motif of approximately 100 residues) secretion system were identified. Here, we report for the first time the identification of a fimbrilin protein exporter (FPE; TC #3.A.14) and a Tad (tight adherence) export apparatus in C. acetobutylicum. This investigation highlights the potential use of this saprophytic bacterium in biotechnological and biomedical applications as well as a model organism for studying protein secretion in pathogenic Gram-positive bacteria.

Mosaic structure of a multiple-drug-resistant, conjugative plasmid from Campylobacter jejuni

Partial sequence analysis of a tet(O) plasmid from a multiple-drug-resistant clinical isolate of Campylobacter jejuni revealed 10 genes or pseudogenes encoding different aminoglycoside inactivating enzymes, transposase-like genes, and multiple unknown genes from a variety of pathogenic and commensal bacteria. The plasmid could be mobilized by a P incompatibility group plasmid into Escherichia coli, where it apparently integrated into the chromosome and expressed high-level resistance to multiple aminoglycoside antibiotics. This work provides new information about both the nature of drug resistance in C. jejuni and the ability of C. jejuni to exchange genes with other bacterial species.

Distribution of F- and A/V-type ATPases in Thermus scotoductus and other closely related species

The presence of an A/V-type ATPase in different Thermus species and in the deeper branching species Meiothermus ruber and Deinococcus radiodurans suggests that the presence of the archaeal-type ATPase is a primitive character of the Deinococci that was acquired through horizontal gene transfer (HGT). However, the presence of a bacterial type F-ATPases was reported in two newly identified Thermus species (Thermus scotoductus DSM 8553 and Thermus filiformis DSM 4687). Two different scenarios can explain this finding, either the recent replacement of the ancestral A/V-type ATPase in Thermus scotoductus and Thermus filiformis with a newly acquired F-type ATPase or a long-term persistence of both F and A type ATPase in the Deinococci, which would imply several independent losses of the F-type ATPase in the Deinococci. Using PCR with redundant primers, sequencing and Southern blot analyses, we tried to confirm the presence of an F-type ATPase in the genome of Thermus scotoductus and Thermus filiformis, and determine its phylogenetic affinities. Initial experiments appeared to confirm the presence of an F-type ATPase in Thermus scotoductus that was similar to the F-ATPases found in Bacillus. However, further experiments revealed that the detection of an F-ATPase was due to a culture contamination. For all the Thermus and Deinococcus species surveyed, including Thermus scotoductus, cultures that were free of contamination only contained an A/V-type ATP synthases.

Characterization of theEscherichia coliO59 and O155 O-antigen gene clusters: The atypicalwzxgenes are evolutionary related

O-antigens are highly polymorphic. The genes specifically involved in O-antigen synthesis are generally grouped together on the chromosome as a gene cluster. In Escherichia coli, the O-antigen gene clusters are characteristically located between the housekeeping genes galF and gnd. In this study, the O-antigen gene clusters of E. coli O59 and E. coli O155 were sequenced. The former was found to contain genes for GDP-mannose synthesis, glycosyltransferase genes and the O-antigen polymerase gene (wzy), while the latter contained only glycosyltransferase genes and wzy. O unit flippase genes (wzx) were found immediately downstream of the gnd gene, in the region between the gnd and hisI genes in these two strains. This atypical location of wzx has not been reported before, and furthermore these two genes complemented in trans despite the fact that different O-antigen structures are present in E. coli O59 and O155. A putative acetyltransferase gene was found downstream of wzx in both strains. Comparison of the region between gnd and hisI revealed that the wzx and acetyltransferase genes are closely related between E. coli O59 and O155, indicating that the two gene clusters arose recently from a common ancestor. This work provides further evidence for the O-antigen gene cluster having formed gradually, and selection pressure will eventually bring O-antigen genes into a single cluster. Genes specific for E. coli O59 and O155, respectively, were also identified.

Prokaryotic K(+) channels: from crystal structures to diversity

The deep roots and wide branches of the K(+)-channel family are evident from genome surveys and laboratory experimentation. K(+)-channel genes are widespread and found in nearly all the free-living bacteria, archaea and eukarya. The conservation of basic structures and mechanisms such as the K(+) filter, the gate, and some of the gate's regulatory domains have allowed general insights on animal K(+) channels to be gained from crystal structures of prokaryotic channels. Since microbes are the great majority of life's diversity, it is not surprising that microbial genomes reveal structural motifs beyond those found in animals. There are open-reading frames that encode K(+)-channel subunits with unconventional filter sequences, or regulatory domains of different sizes and numbers not previously known. Parasitic or symbiotic bacteria tend not to have K(+) channels, while those showing lifestyle versatility often have more than one K(+)-channel gene. It is speculated that prokaryotic K(+) channels function to allow adaptation to environmental and metabolic changes, although the actual roles of these channels in prokaryotes are not yet known. Unlike enzymes in basic metabolism, K(+) channel, though evolved early, appear to play more diverse roles than revealed by animal research. Finding and sorting out these roles will be the goal and challenge of the near future.

Novel RepA-MCM proteins encoded in plasmids pTAU4, pORA1 and pTIK4 from Sulfolobus neozealandicus

Three plasmids isolated from the crenarchaeal thermoacidophile Sulfolobus neozealandicus were characterized. Plasmids pTAU4 (7,192 bp), pORA1 (9,689 bp) and pTIK4 (13,638 bp) show unusual properties that distinguish them from previously characterized cryptic plasmids of the genus Sulfolobus. Plasmids pORA1 and pTIK4 encode RepA proteins, only the former of which carries the novel polymerase-primase domain of other known Sulfolobus plasmids. Plasmid pTAU4 encodes a mini-chromosome maintenance protein homolog and no RepA protein; the implications for DNA replication are considered. Plasmid pORA1 is the first Sulfolobus plasmid to be characterized that does not encode the otherwise highly conserved DNA-binding PlrA protein. Another encoded protein appears to be specific for the New Zealand plasmids. The three plasmids should provide useful model systems for functional studies of these important crenarchaeal proteins.

Transcription activation of a UV-inducible Clostridium perfringens bacteriocin gene by a novel sigma factor

Expression of the plasmid-encoded Clostridium perfringens gene for bacteriocin BCN5 was shown to depend in vivo and in vitro on the activity of UviA protein. UviA, also plasmid-encoded, proved to be an RNA polymerase sigma factor and was also partly autoregulatory. The uviA gene has two promoters; one provided a UviA-independent, basal level of gene expression while the stronger, UviA-dependent promoter was only utilized after the cell experienced DNA damage. As a result, BCN5 synthesis is induced by treatment with UV light or mitomycin C. UviA is related to a special class of sigma factors found to date only in Clostridium species and responsible for activating transcription of toxin genes in Clostridium difficile, Clostridium tetani, and Clostridium botulinum.