During the 2001 anthrax outbreak investigation, 135 B. anthracis isolates were subtyped. The identity of all strains was confirmed with standard microbiologic procedures and the Laboratory Response Network (LRN) testing algorithm (9,10). Isolates were obtained from patients with laboratory-confirmed anthrax (n=10), powders (n=4), and environmental specimens (n=121). For comparison purposes, five B. anthracis isolates originating from New England in the 1960s and 1970s, the Ames strain, and the Pasteur strain were included.

DNA extractions of 28 clinical specimens from six patients with confirmed inhalational anthrax were used for molecular subtyping. These specimens included blood, pleural fluid, blood cultures, serum, cerebrospinal fluid (CSF), lung tissue, and lymph node tissue.

DNA from all strains was prepared with a heat lysis method. Isolates were streaked onto trypticase soy agar containing 5% sheep blood (Becton Dickinson Microbiology Systems, Cockeysville, MD) and incubated overnight at 37°C. A single colony was transferred and dispersed into 0.22-µm centrifugal filter units (Millipore, Bedford, MA) containing 200 µL 10 mM Tris-HCl (pH 8.0). The suspension was heated at 95°C for 20 min and cooled to room temperature. The filter units were then centrifuged in a microfuge at 6,000 x g for 2 min and the filter discarded. The resulting lysate was stored at –20°C until use. DNA from clinical specimens was extracted with a Qiagen DNA Mini Kit per manufacturer’s instructions (Qiagen Inc., Valencia, CA).

MLVA typing was done as described by Keim et al. (6). Briefly, the eight loci were amplified in four reactions: reaction 1 (vrrB₁, CG3, and vrrA), reaction 2 (vrrB₂, pXO1-aat, and pXO2-at), reaction 3 (vrrC₁), and reaction 4 (vrrC₂). In some instances CG3 was removed from reaction 1 and amplified as a 5th reaction because of weak amplification. Each amplicon was labeled with one of three different dyes. Products were separated by polyacrylamide gel electrophoresis under denaturing conditions on an ABI 377 automated DNA sequencer (Applied Biosystems, Foster City, CA), and allele sizes were determined using ABI Genescan software (Applied Biosystems).

The amplification and sequencing of pagA were performed on 42 B. anthracis isolates and 22 clinical specimens as described by Price et al. (8), with the following modifications. Initially, synthetic oligonucleotide polymerase chain reaction (PCR) primers PA-1F and PA-1R and PA-2F and PA-2R (Table 1) were used to amplify two overlapping fragments (1,119 bp and 1,449 bp, respectively) together totaling 2,531 bp and containing the pagA open reading frame (ORF) (8). Because of inconsistent amplification with PA-2F and PA-2R and to generate a single template for sequencing, PCR amplification was performed using primers 1566F and 4205R. In some instances, possibly from the method of DNA purification, 1566F and 4205R did not amplify sufficiently and thus nested PCR was performed using 1575F and 4191R. The primers used in this study were a combination of both published primers (8) and primers designed from the published DNA sequence of the virulence plasmid pXO1 (GenBank accession no. AF065404) (Table 1). The pagA sequencing template was amplified by PCR using the Expand High Fidelity PCR system (Roche, Mannheim, Germany). Fifty-microliter PCR mixtures contained 10 mM Tris-HCL (pH 8.3), 50 mM KCl, 4.0 mM MgCl₂, 0.4 mM of each forward and reverse primer, 100 µM of each deoxynucleotide, 2.0 U of Taq DNA polymerase (Roche), and 2 µL of bacterial lysate. Reactions were heated at 94°C for 5 min and then cycled 35 times at 94°C for 30 s, 51°C for 30 s and 72°C for 1.5 min, with a final extension of 72°C for 5 min. PCR amplicons were purified using QIAquick PCR purification kit (Qiagen, Inc.) and the resulting purified amplicons were used in the subsequent sequencing reactions.

Sequencing was performed on an Applied Biosystems 3100 genetic analyzer (Applied Biosystems) using BigDye terminator cycle sequencing ready reaction mix according to manufacturers instructions (Applied Biosystems). All sequence data were analyzed with the Lasergene 99 (DNASTAR, Madison, WI) software, which comprises several different programs: DNA sequences were assembled using the SeqMan program, and MegAlign was used to do sequence comparisons.

By MLVA, all 135 outbreak-related B. anthracis isolates had the following loci sizes: vrrA = 313, vrrB₁ = 229, vrrB₂ = 153, vrrC₁ = 583, vrrC₂ = 532, CG3 = 158, pXO1 = 123, and pXO2 = 141, resulting in genotype 62, as described by Keim et al. (6). In addition, the entire 2,294-bp pagA gene was sequenced from a subset of 42 isolates: including ten patient isolates, all four powder isolates, and 28 select environmental isolates. All isolates had an indistinguishable sequence, PA genotype I (Table 2).

Before B. anthracis was detected in the mail, we subtyped several isolates from cutaneous anthrax cases that occurred in the 1960s and 1970s in the eastern United States to determine if any were similar to the outbreak strain. Two isolates from Rhode Island were MLVA genotype 71, one New Hampshire isolate was genotype 78, while an additional New Hampshire isolate and a Massachusetts isolate each had unique combinations of alleles resulting in new genotypes. The pagA sequence of all five New England isolates was identified as PA genotype VI (Table 2).

State A reported isolating B. anthracis (2002017388) from an envelope. This state was not in the vicinity of the 2001 outbreak. By MLVA, the isolate was shown to have been cured of pXO1 and had the same genotype as the Pasteur strain, used in laboratories as a control strain for various tests (Table 2).

Country B sent an isolate (2002007581) that was reportedly isolated from a letter to a private physician. MLVA identified the strain as genotype 45, which clearly distinguished it from the strain associated with the ongoing outbreak in the United States. In addition, four other isolates from the same facility were assayed by MLVA (2002007648–51), resulting in the identification of two Sterne strains, one Pasteur strain, and one additional strain of genotype 45 (Table 2).

MLVA and pagA sequencing were performed on clinical specimens collected from seven patients with laboratory-confirmed inhalational anthrax during the 2001 bioterrorism-associated anthrax outbreak. These methods have an unproven utility on clinical specimens, and further testing will be necessary for full evaluation. A total of 28 clinical specimens were analyzed by using MLVA, including: blood, CSF, pleural fluid, serum, lung tissue, and lymph node tissue (Table 3). All eight loci were detected in three specimens (two pleural fluids and one lymph node from patient 10). Of the eight loci examined, vrrA was detected in all nine specimens in which any of the MLVA loci were detected and on two occasions was the only locus detected. The pagA gene was successfully amplified and sequenced from 5 of 22 specimens analyzed (Table 3).

During the 2001 anthrax investigation, molecular subtyping of B. anthracis by MLVA and pagA sequencing was important in linking cases to each other and to contaminated sites and in distinguishing isolates that were not related to this event. We used two methods for the molecular subtyping of B. anthracis: pagA sequencing and MLVA. All outbreak-associated isolates were identified as MLVA genotype 62 and PA genotype I. To date, MLVA genotype 62 has only been associated with a few isolates from herbivores in Texas and has not been identified in any B. anthracis strains originating in eastern United States or anywhere else in the world. None of the New England isolates analyzed in this study were MLVA genotype 62 or PA genotype I. All five were of PA genotype VI, while MLVA identified two isolates as genotype 71, one as genotype 78, and two as new genotypes. Genotype 62 is also the genotype of the Ames strain commonly used in research laboratories worldwide and frequently used in animal challenge studies (11–16). The sequence of pagA from the outbreak strain, PA genotype I, was also identical to that of the Ames strain; thus, the outbreak B. anthracis strain is indistinguishable from the Ames strain based on the examination of the eight MLVA loci and the pagA sequence. Recently, comparative genome sequencing detected only four differences between the chromosomes of the outbreak strain (Florida isolate) and Ames (Porton) isolate (17).

Molecular subtyping of isolates immediately upon their arrival to the laboratory allowed for instant confirmation that the cases were caused by the same strain and thus for linking cases to environmental contamination and to the powder-containing envelopes. The speed of the MLVA allowed for genotype identification within 8 hours of receiving the isolates. In addition to linking the cases, molecular subtyping was invaluable in determining if B. anthracis, isolated from around the world during the same time period, were potentially related to the ongoing outbreak in the United States. The level of discrimination provided by MLVA, allowed for non-outbreak isolates to be rapidly and easily distinguished.

While both MLVA and pagA sequencing are primarily used for molecular subtyping of isolates, we were also able to amplify the eight MLVA loci and pagA directly from a limited number of available clinical specimens. Although this event was not a prospective case-control study, amplification was most successful from pleural fluid and lymph node specimens. Similar results were demonstrated with a B. anthracis specific real-time PCR assay (18). Amplification of the MLVA loci and pagA was not very successful from blood cultures even when taken before antibiotic therapy. The lack of success with blood cultures was not because of a complete inhibition of PCR since the B. anthracis LRN PCR assay was positive on these specimens. Of the MLVA loci, vrrA was the most readily amplified (9 of 28), probably because vrrA amplification is more sensitive compared to the other MLVA loci. However, limits of detection for each of the loci has not been evaluated.

For a single patient (patient 10), we were able to amplify all eight MLVA loci and determine the genotype of the B. anthracis strain without having the isolate itself. In this instance, MLVA was used directly on DNA extracted from pleural fluid and genotype 62 was identified. B. anthracis was not successfully cultured from that same pleural fluid sample. The pagA gene was amplified and sequenced from the same DNA specimen and identified as PA genotype I. MLVA and pagA amplification were attempted on DNA extracted from blood drawn from this patient the same day as the pleural fluid but failed to detect any of the loci, suggesting more efficient clearance of the bacilli from the blood or less sensitivity of these molecular approaches on blood compared to pleural fluid. Again, the negative result on blood was not because of complete inhibition of PCR since the B. anthracis LRN PCR assay on this sample was positive. At a later date, when the Centers for Disease Control and Prevention received the isolate originally cultured from this patient at the local medical facility where the patient was treated, the isolate was confirmed to be genotype 62. Despite the fact that B. anthracis was not successfully cultured from any of these clinical specimens taken after the initiation of antimicrobial therapy, we were able to amplify the MLVA loci and pagA from some of these specimens.

The entire chromosomal sequence of the B. anthracis Ames strain (available from: URL: www.tigr.org) is now available and has been compared to the chromosomal sequence of the outbreak (Florida) isolate (17). While sequencing and comparing B. anthracis genomes are not likely to be useful for rapidly identifying isolates during an outbreak investigation, the data generated from such comparisons may identify new loci, which could be targets for methods such as MLVA and can be done rapidly on large numbers of isolates from patients, the environment, and on DNA from clinical specimens.