Highly precise diagnostics and forensic assays can be developed through a combination of evolutionary analysis and the exhaustive examination of genomic sequences. In Bacillus anthracis, whole-genome sequencing efforts revealed ca. 3,500 single-nucleotide polymorphisms (SNPs) among eight different strains and evolutionary analysis provides the identification of canonical SNPs. We have previously shown that SNPs are highly evolutionarily stable, and the clonal nature of B. anthracis makes them ideal signatures for subtyping this pathogen. Here we identified SNPs that define the lineage of B. anthracis that contains the Ames strain, the strain used in the 2001 bioterrorist attacks in the United States. Sequencing and real-time PCR were used to validate these SNPs across B. anthracis strains, including (i) 88 globally and genetically diverse isolates; (ii) isolates that were shown to be genetic relatives of the Ames strain by multiple-locus variable number tandem repeat analysis (MLVA); and (iii) several different lab stocks of the Ames strain, including a clinical isolate from the 2001 letter attack. Six SNPs were found to be highly specific for the Ames strain; four on the chromosome, one on the pX01 plasmid, and one on the pX02 plasmid. All six SNPs differentiated the B. anthracis Ames strain from the 88 unique B. anthracis strains, while five of the six separated Ames from its close genetic relatives. The use of these SNPs coupled with real-time PCR allows specific and sensitive (<100 fg of template DNA) identification of the Ames strain. This evolutionary and genomics-based approach provides an effective means for the discovery of strain-specific SNPs in B. anthracis.The 2001 anthrax letter attack illustrated the "real-world" efficacy of Bacillus anthracis as a bioterror agent. Forensic and epidemiological analysis of clinical samples and weaponized spores from the letter attack included the identification of the B. anthracis strain as Ames (3,8,17). This was initially accomplished using multiple-locus variable number tandem repeat analysis (MLVA) (3,6,8), which was one of the few technical approaches possible for this highly monomorphic pathogen (5). Identifying the strain and establishing its identity in attack locations were important in linking the dispersed anthrax cases and suggesting a possible source for the weaponized material. During the anthrax attack crisis, diagnostic speed, specificity, and sensitivity were often limiting factors, necessitating extraordinary efforts by public health officials and forensic labs (3). Further advancement of molecular diagnostics allows for more-efficient responses in disease outbreaks, whether natural or bioterrorist mediated.Single-nucleotide polymorphisms (SNPs) are increasingly recognized as important markers for detecting and subtyping bacterial pathogens, including B. anthracis (1,2,4,8,16,18,19). Recent comparative full-genome sequencing allowed the discovery of about 3,500 SNPs among eight strains of B. anthracis (15, 17; J. Ravel, unpublished data) and repre...

BackgroundThe global pattern of distribution of 1033 B. anthracis isolates has previously been defined by a set of 12 conserved canonical single nucleotide polymorphisms (canSNP). These studies reinforced the presence of three major lineages and 12 sub-lineages and sub-groups of this anthrax-causing pathogen. Isolates that form the A lineage (unlike the B and C lineages) have become widely dispersed throughout the world and form the basis for the geographical disposition of "modern" anthrax. An archival collection of 191 different B. anthracis isolates from China provides a glimpse into the possible role of Chinese trade and commerce in the spread of certain sub-lineages of this pathogen. Canonical single nucleotide polymorphism (canSNP) and multiple locus VNTR analysis (MLVA) typing has been used to examine this archival collection of isolates.ResultsThe canSNP study indicates that there are 5 different sub-lineages/sub-groups in China out of 12 previously described world-wide canSNP genotypes. Three of these canSNP genotypes were only found in the western-most province of China, Xinjiang. These genotypes were A.Br.008/009, a sub-group that is spread across most of Europe and Asia; A.Br.Aust 94, a sub-lineage that is present in Europe and India, and A.Br.Vollum, a lineage that is also present in Europe. The remaining two canSNP genotypes are spread across the whole of China and belong to sub-group A.Br.001/002 and the A.Br.Ames sub-lineage, two closely related genotypes. MLVA typing adds resolution to the isolates in each canSNP genotype and diversity indices for the A.Br.008/009 and A.Br.001/002 sub-groups suggest that these represent older and established clades in China.ConclusionB. anthracis isolates were recovered from three canSNP sub-groups (A.Br.008/009, A.Br.Aust94, and A.Br.Vollum) in the western most portion of the large Chinese province of Xinjiang. The city of Kashi in this province appears to have served as a crossroads for not only trade but the movement of diseases such as anthrax along the ancient "silk road". Phylogenetic inference also suggests that the A.Br.Ames sub-lineage, first identified in the original Ames strain isolated from Jim Hogg County, TX, is descended from the A.Br.001/002 sub-group that has a major presence in most of China. These results suggest a genetic discontinuity between the younger Ames sub-lineage in Texas and the large Western North American sub-lineage spread across central Canada and the Dakotas.

Disease introduction into the New World during colonial expansion is well documented and had a major impact on indigenous populations; however, few diseases have been associated with early human migrations into North America. During the late Pleistocene epoch, Asia and North America were joined by the Beringian Steppe ecosystem which allowed animals and humans to freely cross what would become a water barrier in the Holocene. Anthrax has clearly been shown to be dispersed by human commerce and trade in animal products contaminated with Bacillus anthracis spores. Humans appear to have brought B. anthracis to this area from Asia and then moved it further south as an ice-free corridor opened in central Canada ∼13,000 ybp. In this study, we have defined the evolutionary history of Western North American (WNA) anthrax using 2,850 single nucleotide polymorphisms (SNPs) and 285 geographically diverse B. anthracis isolates. Phylogeography of the major WNA B. anthracis clone reveals ancestral populations in northern Canada with progressively derived populations to the south; the most recent ancestor of this clonal lineage is in Eurasia. Our phylogeographic patterns are consistent with B. anthracis arriving with humans via the Bering Land Bridge. This northern-origin hypothesis is highly consistent with our phylogeographic patterns and rates of SNP accumulation observed in current day B. anthracis isolates. Continent-wide dispersal of WNA B. anthracis likely required movement by later European colonizers, but the continent's first inhabitants may have seeded the initial North American populations.

To determine potential links between the clinical isolate to animal products and their geographic origin, we genotyped (MLVA-8, MVLA-15, and canSNP analysis) 80 environmental and 12 clinical isolates and 2 clinical specimens from five cases of anthrax (California in 1976 [n = 1], New York in 2006 [n = 1], Connecticut in 2007 [n = 2], and New Hampshire in 2009[n = 1]) resulting from recreational handling of animal products. For the California case, four clinical isolates were identified as MLVA-8 genotype (GT) 76 and in the canSNP A.Br.Vollum lineage, which is consistent with the Pakistani origin of the yarn. Twenty eight of the California isolates were in the A.Br.Vollum canSNP lineage and one isolate was in the A.Br. 003/004 canSNP sub-group. All 52 isolates and both clinical specimens related to the New York and Connecticut cases were MLVA-8 GT 1. The animal products associated with the NY and CT cases were believed to originate from West Africa, but no isolates from this region are available to be genotyped for comparison. All isolates associated with the New Hampshire case were identical and had a new genotype (GT 149). Isolates from the NY, CT and NH cases diverge from the established canSNP phylogeny near the base of the A.Br.011/009. This report illustrates the power of the current genotyping methods and the dramatically different epidemiological conditions that can lead to infections (i.e., contamination by a single genotype versus widespread contamination of numerous genotypes). These cases illustrate the need to acquire and genotype global isolates so that accurate assignments can be made about isolate origins.

No abstract