Anthrax, caused by the bacterium Bacillus anthracis, is a disease of historical and current importance that is found throughout the world. The basis of its historical transmission is anecdotal and its true global population structure has remained largely cryptic. Seven diverse B. anthracis strains were whole-genome sequenced to identify rare single nucleotide polymorphisms (SNPs), followed by phylogenetic reconstruction of these characters onto an evolutionary model. This analysis identified SNPs that define the major clonal lineages within the species. These SNPs, in concert with 15 variable number tandem repeat (VNTR) markers, were used to subtype a collection of 1,033 B. anthracis isolates from 42 countries to create an extensive genotype data set. These analyses subdivided the isolates into three previously recognized major lineages (A, B, and C), with further subdivision into 12 clonal sub-lineages or sub-groups and, finally, 221 unique MLVA15 genotypes. This rare genomic variation was used to document the evolutionary progression of B. anthracis and to establish global patterns of diversity. Isolates in the A lineage are widely dispersed globally, whereas the B and C lineages occur on more restricted spatial scales. Molecular clock models based upon genome-wide synonymous substitutions indicate there was a massive radiation of the A lineage that occurred in the mid-Holocene (3,064–6,127 ybp). On more recent temporal scales, the global population structure of B. anthracis reflects colonial-era importation of specific genotypes from the Old World into the New World, as well as the repeated industrial importation of diverse genotypes into developed countries via spore-contaminated animal products. These findings indicate humans have played an important role in the evolution of anthrax by increasing the proliferation and dispersal of this now global disease. Finally, the value of global genotypic analysis for investigating bioterrorist-mediated outbreaks of anthrax is demonstrated.

Bacillus anthracis is one of the most genetically homogeneous pathogens described, making strain discrimination particularly difficult. In this paper, we present a novel molecular typing system based on rapidly evolving variable-number tandem repeat (VNTR) loci. Multiple-locus VNTR analysis (MLVA) uses the combined power of multiple alleles at several marker loci. In our system, fluorescently labeled PCR primers are used to produce PCR amplification products from eight VNTR regions in the B. anthracis genome. These are detected and their sizes are determined using an ABI377 automated DNA sequencer. Five of these eight loci were discovered by sequence characterization of molecular markers (vrrC 1 , vrrC 2 , vrrB 1 , vrrB 2 , and CG3), two were discovered by searching complete plasmid nucleotide sequences (pXO1-aat and pXO2-at), and one was known previously (vrrA). MLVA characterization of 426 B. anthracis isolates identified 89 distinct genotypes. VNTR markers frequently identified multiple alleles (from two to nine), with Nei's diversity values between 0.3 and 0.8. Unweighted pair-group method arithmetic average cluster analysis identified six genetically distinct groups that appear to be derived from clones. Some of these clones show worldwide distribution, while others are restricted to particular geographic regions. Human commerce doubtlessly has contributed to the dispersal of particular clones in ancient and modern times.

Clostridium botulinum is a taxonomic designation for many diverse anaerobic spore-forming rod-shaped bacteria that have the common property of producing botulinum neurotoxins (BoNTs). The BoNTs are exoneurotoxins that can cause severe paralysis and death in humans and other animal species. A collection of 174 C. botulinum strains was examined by amplified fragment length polymorphism (AFLP) analysis and by sequencing of the 16S rRNA gene and BoNT genes to examine the genetic diversity within this species. This collection contained representatives of each of the seven different serotypes of botulinum neurotoxins (BoNT/A to BoNT/G). Analysis of the16S rRNA gene sequences confirmed previous identifications of at least four distinct genomic backgrounds (groups I to IV), each of which has independently acquired one or more BoNT genes through horizontal gene transfer. AFLP analysis provided higher resolution and could be used to further subdivide the four groups into subgroups. Sequencing of the BoNT genes from multiple strains of serotypes A, B, and E confirmed significant sequence variation within each serotype. Four distinct lineages within each of the BoNT A and B serotypes and five distinct lineages of serotype E strains were identified. The nucleotide sequences of the seven toxin genes of the serotypes were compared and showed various degrees of interrelatedness and recombination, as was previously noted for the nontoxic nonhemagglutinin gene, which is linked to the BoNT gene. These analyses contribute to the understanding of the evolution and phylogeny within this species and assist in the development of improved diagnostics and therapeutics for the treatment of botulism.Clostridium botulinum is a taxonomic collection of several distinct species of anaerobic gram-positive spore-forming bacteria that produce the most poisonous substance known, botulinum neurotoxin (BoNT) (1,8). These organisms, along with related neurotoxin-producing species that, for a variety of reasons, were not included under the C. botulinum taxon, pose global health problems that affect both infant and adult humans and can also affect wildlife, waterfowl, and domestic animals. They cause intoxication through ingestion of the neurotoxin in contaminated foods. Toxicoinfections can also occur after contact with bacteria or bacterial spores (6, 17). These pathogens are ubiquitous and can be found in soils and sedi-

No abstract