Molecular epidemiology of anthrax cases associated with recreational use of animal hides and yarn in the United States.

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.

Cost-effective interrogation of single nucleotide polymorphisms using the mismatch amplification mutation assay and capillary electrophoresis.

The ability to characterize SNPs is an important aspect of many clinical diagnostic, genetic and evolutionary studies. Here, we designed a multiplexed SNP genotyping method to survey a large number of phylogenetically informative SNPs within the genome of the bacterium Bacillus anthracis. This novel method, CE universal tail mismatch amplification mutation assay (CUMA), allows for PCR multiplexing and automatic scoring of SNP genotypes, thus providing a rapid, economical and higher throughput alternative to more expensive SNP genotyping techniques. CUMA delivered accurate B. anthracis SNP genotyping results and, when multiplexed, saved reagent costs by more than 80% compared with TaqMan real-time PCR. When real-time PCR technology and instrumentation is unavailable or the reagents are cost-prohibitive, CUMA is a powerful alternative for SNP genotyping.

Accurate, rapid and high-throughput detection of strain-specific polymorphisms in Bacillus anthracis and Yersinia pestis by next-generation sequencing.

BACKGROUND: In the event of biocrimes or infectious disease outbreaks, high-resolution genetic characterization for identifying the agent and attributing it to a specific source can be crucial for an effective response. Until recently, in-depth genetic characterization required expensive and time-consuming Sanger sequencing of a few strains, followed by genotyping of a small number of marker loci in a panel of isolates at or by gel-based approaches such as pulsed field gel electrophoresis, which by necessity ignores most of the genome. Next-generation, massively parallel sequencing (MPS) technology (specifically the Applied Biosystems sequencing by oligonucleotide ligation and detection (SOLiD™) system) is a powerful investigative tool for rapid, cost-effective and parallel microbial whole-genome characterization. RESULTS: To demonstrate the utility of MPS for whole-genome typing of monomorphic pathogens, four Bacillus anthracis and four Yersinia pestis strains were sequenced in parallel. Reads were aligned to complete reference genomes, and genomic variations were identified. Resequencing of the B. anthracis Ames ancestor strain detected no false-positive single-nucleotide polymorphisms (SNPs), and mapping of reads to the Sterne strain correctly identified 98% of the 133 SNPs that are not clustered or associated with repeats. Three geographically distinct B. anthracis strains from the A branch lineage were found to have between 352 and 471 SNPs each, relative to the Ames genome, and one strain harbored a genomic amplification. Sequencing of four Y. pestis strains from the Orientalis lineage identified between 20 and 54 SNPs per strain relative to the CO92 genome, with the single Bolivian isolate having approximately twice as many SNPs as the three more closely related North American strains. Coverage plotting also revealed a common deletion in two strains and an amplification in the Bolivian strain that appear to be due to insertion element-mediated recombination events. Most private SNPs (that is, a, variant found in only one strain in this set) selected for validation by Sanger sequencing were confirmed, although rare false-positive SNPs were associated with variable nucleotide tandem repeats. CONCLUSIONS: The high-throughput, multiplexing capability, and accuracy of this system make it suitable for rapid whole-genome typing of microbial pathogens during a forensic or epidemiological investigation. By interrogating nearly every base of the genome, rare polymorphisms can be reliably discovered, thus facilitating high-resolution strain tracking and strengthening forensic attribution.

Bacillus anthracis in China and its relationship to worldwide lineages.

BACKGROUND: The 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. RESULTS: The 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. CONCLUSION: B. 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.

Pre-Columbian origins for North American anthrax.

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 approximately 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.

Single nucleotide polymorphism typing of Bacillus anthracis from Sverdlovsk tissue.

A small number of conserved canonical single nucleotide polymorphisms (canSNP) that define major phylogenetic branches for Bacillus anthracis were used to place a Sverdlovsk patient's B. anthracis genotype into 1 of 12 subgroups. Reconstruction of the pagA gene also showed a unique SNP that defines a new lineage for B. anthracis.

Role of law enforcement response and microbial forensics in investigation of bioterrorism.

The risk and threat of bioterrorism and biocrime have become a large concern and challenge for governments and society to enhance biosecurity. Law enforcement plays an important role in assessing and investigating activities involved in an event of bioterrorism or biocrime. Key to a successful biosecurity program is increased awareness and early detection of threats facilitated by an integrated network of responsibilities and capabilities from government, academic, private, and public assets. To support an investigation, microbial forensic sciences are employed to analyze and characterize forensic evidence with the goal of attribution or crime scene reconstruction. Two different molecular biology-based assays--real time polymerase chain reaction (PCR) and repetitive element PCR--are described and demonstrate how molecular biology tools may be utilized to aid in the investigative process. Technologies relied on by microbial forensic scientists need to be properly validated so that the methods used are understood and so that interpretation of results is carried out within the limitations of the assays. The three types of validation are preliminary, developmental, and internal. The first is necessary for rapid response when a threat is imminent or an attack has recently occurred. The latter two apply to implementation of routinely used procedures.

Strain-specific single-nucleotide polymorphism assays for the Bacillus anthracis Ames strain.

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.