Molecular Subtyping of Bacillus anthracis and the 2001 Bioterrorism-Associated Anthrax Outbreak, United States. (Bioterrorism-Related Anthrax).

Molecular subtyping of Bacillus anthracis played an important role in differentiating and identifying strains during the 2001 bioterrorism-associated outbreak. Because B. anthracis has a low level of genetic variability, only a few subtyping methods, with varying reliability, exist. We initially used multiple-locus variable-number tandem repeat analysis (MLVA) to subtype 135 B. anthracis isolates associated with the outbreak. All isolates were determined to be of genotype 62, the same as the Ames strain used in laboratories. We sequenced the protective antigen gene (pagA) from 42 representative outbreak isolates and determined they all had a pagA sequence indistinguishable from the Ames strain (PA genotype I). MLVA and pagA sequencing were also used on DNA from clinical specimens, making subtyping B. anthracis possible without an isolate. Use of high-resolution molecular subtyping determined that all outbreak isolates were indistinguishable by the methods used and probably originated from a single source. In addition, subtyping rapidly identified laboratory contaminants and nonoutbreak-related isolates.

**********

The recent bioterrorism-associated anthrax outbreak demonstrated the need for rapid molecular subtyping of Bacillus anthracis isolates. Numerous methods, including multiple-locus enzyme electrophoresis (MEE) and multiple-locus sequence typing (MLST), have shown the lack of genetic diversity of B. anthracis (1-4, unpub, data). Despite this low diversity, methods have been developed that can detect differences between B. anthracis isolates. Amplified fragment length polymorphism (AFLP) analysis has been used to detect differences between B. anthracis isolates and to examine phylogenetic relationships between B. anthracis and its close relatives, B. cereus and B. thuringiensis (4,5). Keim et al. (6) reported on multiple-locus variable-number tandem repeat analysis (MLVA) for subtyping B. anthracis, which unlike AFLP is designed to subtype B. anthracis specifically and cannot be used to address phylogenetic relationships between Bacillus species. MLVA determines the copy number of variable-number tandem repeats (VNTR) at eight genetic loci (six chromosomal and one on each of the two plasmids). Recently, MLVA has been used to differentiate 426 B. anthracis isolates into 89 distinct genotypes and to study the ecology of anthrax (6,7). MLVA is relatively simple, has excellent reproducibility, can subtype multiple strains on a single gel, and gives results in <8 hours.

Protective antigen (PA) is one of the three anthrax toxin proteins and is key to developing immunity to anthrax. Sequencing the gene that encodes PA (pagA) has been used to subtype 26 diverse B. anthracis isolates into six PA genotypes (8). Although sequencing of pagA results in limited numbers of subtypes, it does have the added benefit of determining if the pagA gene has been altered or engineered.

During the 2001 bioterrorism-associated anthrax outbreak, we used MLVA to subtype isolates from patients, the environment, and powders. Subtyping of B. anthracis allowed anthrax cases to be linked to environmental specimens and powders and provided information about potential sources. Sequencing of pagA was also performed on a subset of these B. anthracis isolates, and we confirmed that the pagA sequence was not altered. In addition, we used these methods on DNA extracted from select clinical specimens to detect and subtype B. anthracis directly from clinical specimens. During the outbreak, laboratories throughout the United States and around the world received an increased number of specimens to be tested for B. anthracis. With such large numbers, occasional contamination or detection of non-outbreak strains was inevitable, and molecular subtyping was used to clarify these situations on several occasions. Overall, the recent anthrax outbreak has dramatically illustrated the importance of rapid molecular subtyping during a bioterrorism event.

Materials and Methods

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[degrees]C. A single colony was transferred and dispersed into 0.22-gm centrifugal filter units (Millipore, Bedford, MA) containing 200 [micro]L 10 mM Tris-HCl (pH 8.0). The suspension was heated at 95[degrees]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[degrees]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.sub.1], CG3, and [vrrA]), reaction 2 ([vrrB.sub.2], pXO1-aat, and pXO2-at), reaction 3 (vrr[C.sub.1]), and reaction 4 (vrr[C.sub.2]). 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 10mM Tris-HCL (pH 8.3), 50 mM KCl, 4.0 mM Mg[Cl.sub.2], 0.4 mM of each forward and reverse primer, 100 gM of each deoxynucleotide, 2.0 U of Taq DNA polymerase (Roche), and 2 [micro]L of bacterial lysate. Reactions were heated at 94[degrees]C for 5 min and then cycled 35 times at 94[degrees]C for 30 s, 51[degrees]C for 30 s, and 72[degrees]C for 1.5 min, with a final extension of 72[degrees]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 manufacturer's 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.

Results

By MLVA, all 135 outbreak-related B. anthracis isolates had the following loci sizes: vrrA = 313, vrr[B.sub.1] = 229, vrr[B.sub.2] = 153, vrr[C.sub.1] = 583, vrr[C.sub.2] = 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).

Discussion

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). This is likely the result of a lower limit of detection for vrrA compared to the other loci; however, limits of detection for each of the loci have 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.

Table 1. Primers used for amplification and sequencing of the
pagA gene

Primer                             Sequence (5'-3')

PA1575F (a)                   CGA ACT GAT ACA CGT ATT TTA G
PA4191R (a)                  AGG ATT ATG ATG ATT TAG ATT ACT
PA1566F (a)                  TTT ATC CGA ACT GAT ACA CGT ATT
PA4205R (a)                  ACA AAC AAT CTC AAA GGA TTA TGA
PA-1F (a)                 ATA TTT ATA AAA GTT CTG TTT AAA AAG CC
PA-1R (a)                    TAA ATC CTG CAG ATA CAC TCC CAC
PA-2F (a)               ATA AGT AAA AAT ACT TCT ACA AGT AGG ACA C
PA-2R (a)               GAT TTA GAT TAC TGT TTA AAA CAT ACT CTC C
PA-3                        TCA TGT AAC AAT GTG GGT AGA TGA C
PA-4                        CTC TAT GAG CCT CCT TAA CTA CTG AC
PA-5F                       ATC CTA GTG ATC CAT TAG AAA CGA C
PA-5R                      CTT CTC TAT GAG CCT CCT TAA CTA CTG
PA-5[F.sub.nest]                AGT GAT CCA TTA GAA ACG AC
PA-5[R.sub.nest]                TAA CTA CTG ACT CAT CCG C
PA-U2121                     TAC ATT TGC TAC TTC CGC TGA TAA
PA-L3892                       TGT TTT TCC ATC TTG CCG TAA
2121R                        TTA TCA GCG GAA GTA GCA AAT GTA
3892F                          TTA CGG CAA GAT GGA AAA ACA
2557R                          AGC CGT GCT CCA TTT TTC AGG
3318R                          TGC GGT AAC ACT BTCA CTC CAG
2560F                          GAA AAA TGG AGC ACG GCT TCT
2924F                          CTT GGG CTG AAA CAA TGG GTT

(a) Primers used for amplification of the pagA gene. All other primers
were used for sequencing.

Table 2. MLVA and pagA genotyping of Bacillus anthracis isolates (a,b)

B. anthracis strain          No. strains    vrrA     vrr[B.sub.1]

Outbreak-associated              135         313         229
Ames                               1         313         229
NH (2000032764)                    1         301         229
NH (2000032760)                    1         313         229
RI (2000032763)                    1         313         229
RI (2000032761)                    1         313         229
MA (2000032762)                    1         313         229
State A (2002017388)               1         313         229
Pasteur                            1         313         229
Country B (2002007581)             1         313         229
Country B (2002007648)             1         313         229
Country B (2002007649)             1         313         229
Country B (2002007650)             1         313         229
Country B (2002007651)             1         313         229

B. anthracis strain           vrr[B.sub2]    vrr[C.sub.1]        vrr
                                                              [C.sub.2]

Outbreak-associated               153             583             532
Ames                              153             583             532
NH (2000032764)                   153             538             604
NH (2000032760)                   153             538             604
RI (2000032763)                   162             538             604
RI (2000032761)                   162             538             604
MA (2000032762)                   153             538             604
State A (2002017388)              162             613             604
Pasteur                           162             613             604
Country B (2002007581)            162             613             532
Country B (2002007648)            162             613             532
Country B (2002007649)            162             613             604
Country B (2002007650)            162             583             532
Country B (2002007651)            162             583             532

B. anthracis strain               CG3    pXO1    pXO2      MLVA type

Outbreak-associated               158     123     141         62
Ames                              158     123     141         62
NH (2000032764)                   158     132     139         78
NH (2000032760)                   158     123     139       New (e)
RI (2000032763)                   158     132     139         71
RI (2000032761)                   158     132     139         71
MA (2000032762)                   158     132     143       New (e)
State A (2002017388)              153      -      137     Pasteur (c)
Pasteur                           153      -      137     Pasteur (c)
Country B (2002007581)            158     129     141         45
Country B (2002007648)            158     129     141         45
Country B (2002007649)            153      -      137     Pasteur (c)
Country B (2002007650)            158     129              Sterne (c)
Country B (2002007651)            158     129              Sterne (c)

B. anthracis strain          PA genotype

Outbreak-associated               I
Ames                              I
NH (2000032764)                  VI
NH (2000032760)                  VI
RI (2000032763)                  VI
RI (2000032761)                  VI
MA (2000032762)                  VI
State A (2002017388)           NA (d)
Pasteur                        NA (d)
Country B (2002007581)            I
Country B (2002007648)            I
Country B (2002007649)         NA (d)
Country B (2002007650)            I
Country B (2002007651)            I

(a) MLVA, multiple-locus variable-number tandem repeat analysis; pagA,
protective antigen gene; PA, protective antigen; --, loci not detected;
NA, not applicable.

(b) Allele size for each VNTR locus is shown in addition to the MLVA
and PA genotypes.

(c) No MLVA genotype assigned due to the lack of one of the virulence
plasmids (pXO1 or pXO2).

(d) pagA not present in pXO1-cured strains and thus could not be
assigned a PA genotype.

(e) New combination of alleles resulting in a new genotype. Genotype
no. to be assigned at a later date.

Table 3. Molecular subtyping by MLVA and pagA sequencing performed on
28 clinical specimens from seven patients with inhalational anthrax
(a,b)

                                                 Interval after ant
Patient no. (b)          Specimen type            microbial therapy
                                                      (days) (d)

1                      Pleural fluid (c)                     4
                       Pleural fluid (c)                     4
                           Blood (c)                         4
                            Lung (c)                         4
                            Lung (c)                         4
                        Heart blood (c)                      4
                     Pericardial blood (c)                   4
2                     Thoracentesis fluid                    4
                             Serum                          10
                        Respiratory wash                     4
                         Pleural fluid                       4
3                        Blood culture                       0
5                        Blood culture                       0
                         Blood culture                       0
6                        Blood culture                       0
10                       Pleural fluid                       1
                         Pleural fluid                       1
                             Blood                           1
                            CSF (c)                          3
                           Lung ( c)                         3
                         Lymph node (c)                      3
11                       Pleural fluid                       2
                             Blood                           2
                         Blood culture                      -1
                         Blood culture                      -1
                         Blood culture                      -1
                         Blood culture                      -1
                         Lymph node (c)                      4

Patient no. (b)          MLVA loci detected              pagA

1                           All negative                Negative
                            All negative                Negative
                            All negative                Negative
                            All negative                   ND
                                vrrA                       ND
                         vrrA, vrr[B.sub.1]             Negative
                            All negative                Negative
2                               vrrA                    Negative
                            All negative                Negative
                            All negative                Negative
                            All negative                Negative
3                           All negative                Negative
5                           All negative                Negative
                            All negative                Negative
6                vrrA, vrr[B.sub.1], vrr[B.sub.2],      Negative
                            vrr[C.sub.2]
10               vrrA, vrr[B.sub.1], vrr[B.sub.2],      Positive
                  vrr[C.sub.1], vrr[C.sub.2], CG3,
                            pXO 1, pXO2
                 vrrA, vrr[B.sub.1], vrr[B.sub.2],      Positive
                  vrr[C.sub.1], vrr[C.sub.2], CG3,
                             pXO1, pXO2
                            All negative                Negative
                            All negative                Negative
                             vrrA, CG3                  Negative
                 vrrA, vrr[B.sub.1], vrr[B.sub.2],      Positive
                  vrr[C.sub.1], vrr[C.sub.2], CG3,
                             pXO1, pXO2
11                          All negative                Positive
                            All negative                Negative
                            All negative                   ND
                            All negative                   ND
                            All negative                   ND
                            All negative                   ND
                 vrrA, vrr[B.sub.1], vrr[B.sub.2],      Positive
                  vrr[C.sub.1], vrr[C.sub.2], CG3

                 Bacillus anthracis
Patient no. (b)     LRN PCR (e)

1                   Positive
                    Positive
                    Negative
                       ND
                    Negative
                    Negative
                    Positive
2                   Positive
                    Positive
                    Positive
                    Positive
3                   Positive
5                   Positive
                    Positive
6                   Positive
10                  Positive
                    Positive
                    Positive
                    Positive
                    Positive
                    Positive
11                  Positive
                    Negative
                    ND
                    Positive
                    ND
                    Positive
                    Positive

(a) MLVA, multiple-locus variable-number tandem repeat analysis; pagA,
protective antigen gene; LRN, Laboratory Response Network; PCR,
polymerase chain reaction.

(b) patients 1-10 described in Jernigan et al. (19) and patient 11 in
Barakat et al (20)

(c) Specimens collected postmortem.

(d) Number of days the specimen was collected following or before the
initiation of antimicrobial therapy. Specimens collected the same day
as the initiation of therapy were designated as day 0 but were
collected before antibiotic therapy.

(e) Results using the Laboratory Response Network PCR assay for
detection of B. anthracis during the outbreak (18).

Acknowledgments

We are grateful to Paul Keim, Kimothy Smith, and Robbin Weyant for their input and contributions.

Dr. Hoffmaster is a microbiologist in the Meningitis and Special Pathogens Branch, Centers for Disease Control and Prevention. His interests include research on methods for rapid detection and molecular subtyping of Brucella spp., Burkholderia spp., and Bacillus anthracis.

References

(1.) Harrell L J, Andersen GL, Wilson KH. Genetic variability of Bacillus anthracis and related species. J Clin Microbiol 1995;33:1847-50.

(2.) Henderson ID, Yu D, Tumbull PC. Differentiation of Bacillus anthracis and other `Bacillus cereus group' bacteria using IS231-derived sequences. FEMS Microbiol Lett 1995;128:113-8.

(3.) Helgason E, Okstad OA, Caugant DA, Johansen HA, Fouet A, Mock M, et al. Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis--one species on the basis of genetic evidence. Appl Environ Microbiol 2000;66:2627-30.

(4.) Keim P, Kalif A, Schupp J, Hill K, Travis SE, Richmond K, et al. Molecular evolution and diversity in Bacillus anthracis as detected by amplified fragment length polymorphism markers. J Bacteriol 1997;179:818-24.

(5.) Ticknor LO, Kolsto AB, Hill KK, Keim P, Laker MT, Tonks M, et al. Fluorescent amplified fragment length polymorphism analysis of Norwegian Bacillus cereus and Bacillus thuringiensis soil isolates. Appl Environ Microbiol 2001;67:4863-73.

(6.) Keim P, Price LB, Klevytska AM, Smith KL, Schupp JM, Okinaka R, et al. Multiple-locus variable-number tandem repeat analysis reveals genetic relationships within Bacillus anthracis. J Bacteriol 2000;182:2928-36.

(7.) Smith KL, DeVos V, Bryden H, Price LB, Hugh-Jones ME, Keim P. Bacillus anthracis diversity in Kruger National Park. J Clin Microbiol 2000;38:3780-4.

(8.) Price LB, Hugh-Jones M, Jackson P J, Keim P. Genetic diversity in the protective antigen gene of Bacillus anthracis. J Bacteriol 1999;181:2358-62.

(9.) Khan AS, Morse S, Lillibridge SR. Public-health preparedness for biological terrorism in the USA. Lancet 2000;356:1179-82.

(10.) Logan NA, Tumbull PC. Bacillus and recently derived genera. In: Murray P, editor. Manual of clinical microbiology. Washington: ASM Press; 1991. p. 357-69.

(11.) McBride BW, Mogg A, Telfer JL, Lever MS, Miller J, Tumbull PC, et al. Protective efficacy of a recombinant protective antigen against Bacillus anthracis challenge and assessment of immunological markers. Vaccine 1998;16:810-7.

(12.) Barnard JP, Friedlander AM. Vaccination against anthrax with attenuated recombinant strains of Bacillus anthracis that produce protective antigen. Infect Immun 1999;67:562-7.

(13.) Little SF, Ivins BE, Fellows PF, Friedlander AM. Passive protection by polyclonal antibodies against Bacillus anthracis infection in guinea pigs. Infect Immun 1997;65:5171-5.

(14.) Welkos S, Little S, Friedlander A, Fritz D, Fellows P. The role of antibodies to Bacillus anthracis and anthrax toxin components in inhibiting the early stages of infection by anthrax spores. Microbiol 2001;147:1677-85.

(15.) Ivins BE, Pitt ML, Fellows PF, Farchaus JW, Benner GE, Waag DM, et al. Comparative efficacy of experimental anthrax vaccine candidates against inhalation anthrax in Rhesus macaques. Vaccine 1998;16:1141-8.

(16.) Ivins B, Fellows P, Pitt L, Estep J, Farchaus J, Friedlander A, et al. Experimental anthrax vaccines: efficacy of adjuvants combined with protective antigen against an aerosol Bacillus anthracis spore challenge in guinea pigs. Vaccine 1995;13:1779-84.

(17.) Read TD, Salzberg SL, Pop M, Shumway M, Umayam L, Jiang L, et al. Comparative genome sequenincg for discovery of novel polymorphisms in Bacillus anthracis. Science 2002;296:2028-33.

(18.) Hoffmaster AR, Meyer RF, Bowen M, Marston CK, Weyant RS, Barnett GA, et al. Evaluation and validation of a real-time polymerase chain reaction assay for rapid identification of Bacillus anthracis. Emerg Infect Dis 2002;8;1178-82.

(19.) Jemigan JA, Stephens DS, Ashford DA, Omenaca C, Topiel MS, Galbraith M, et al. Bioterrorism-related inhalational anthrax: the first 10 cases reported in the United States. Emerg Infect Dis 2001;7:933-44.

(20.) Barakat LA, Quentzel HL, Jemigan JA, Kirschke DL, Griffith K, Spear SM, et al. Fatal inhalational anthrax in a 94-year-old Connecticut woman. JAMA 2002;287:863-68.

Address for correspondence: Alex R. Hoffmaster; Meningitis and Special Pathogens Branch, Division of Bacterial and Mycotic Diseases, National Center for Infectious Diseases, Centers for Disease Control and Prevention, Mailstop G34, 1600 Clifton Road, N.E., Atlanta, GA 30333, USA; fax: [phone]; e-mail: [e-mail]

Alex R. Hoffmaster, * Collette C. Fitzgerald, * Efrain Ribot,* Leonard W. Mayer, * and Tanja Popovic *

* Centers for Disease Control and Prevention, Atlanta, Georgia, USA