Previous Article
Applied and Environmental Microbiology, May 2007, p. 3446-3449, Vol. 73, No. 10
0099-2240/07/$08.00+0 doi:10.1128/AEM.02478-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Differentiation of Bacillus anthracis, B. cereus, and B. thuringiensis by Using Pulsed-Field Gel Electrophoresis
Wenwan Zhong,2,
Yulin Shou,1
Thomas M. Yoshida,2 and
Babetta L. Marrone1*
Bioscience,1
Chemistry Divisions, Los Alamos National Laboratory, Los Alamos, New Mexico 875452
Received 23 October 2006/
Accepted 18 March 2007

ABSTRACT
A pulsed-field gel electrophoresis (PFGE) method was developed
for discriminating
Bacillus anthracis from
B. cereus and
B. thuringiensis. A worldwide collection of 25
B. anthracis isolates
showed high-profile homology, and these isolates were unambiguously
distinguished from
B. cereus and
B. thuringiensis isolates by
cluster analysis of the whole-genome macrorestriction enzyme
digestion patterns generated by NotI.

INTRODUCTION
Bacillus anthracis,
B. cereus, and
B. thuringiensis are closely
related species but have different virulence characteristics
(
2,
5,
16). While
B. anthracis is the etiological agent of anthrax
and has recognized potential for use as a biological weapon,
B. cereus is often associated with human food poisoning and
B. thuringiensis produces insect toxins (
13). Therefore, reliable
methods for the discrimination of
B. anthracis from
B. cereus and
B. thuringiensis are important for controlling disease outbreaks
and for determining their sources, which could be naturally
occurring or of malicious intent. However, the high degree of
morphological and biochemical similarities shared by these three
species poses great difficulties for differentiation. Pulsed-field
gel electrophoresis (PFGE) has been adopted by the Centers for
Disease Control and Prevention and all 50 state public health
labs to establish PulseNet, a public health laboratory network
for food-borne disease surveillance (
15). Because
B. anthracis and
B. cereus infections may be food borne, a reliable PFGE
method for distinguishing
B. anthracis from its
B. cereus and
B. thuringiensis near neighbors would be useful for enhancing
the public health response to food-borne outbreaks and to suspected
bioterrorism. The objective of this study was to develop a PFGE
method for reliably genotyping
B. anthracis and distinguishing
B. anthracis isolates from those of
B. cereus and
B. thuringiensis.
Table 1 lists the bacterial strains used in this study, which were obtained either from Dugway Proving Ground, Utah, or from in-house collaborators. A diverse panel of B. anthracis strains was chosen to represent different phylogenetic groups (9, 10). The B. cereus strains D17, 3A, S2-8, F1-15, and ATCC 4342 were isolated from food-borne sources. B. thuringiensis 97-27 has apparent pathogenic properties that are unusual for B. thuringiensis. The other two B. thuringiensis strains, HD 571 and Al Hakam, are nonpathogenic (12). The remaining B. cereus and B. thuringiensis isolates either were demonstrated to be closely related to B. anthracis by amplified fragment length polymorphism analysis and multiple-locus sequence typing or caused anthraxlike symptoms in animals or humans (B. cereus E33L, G924, and TX strains [6, 7, 8, 14]). Two B. atrophaeus samples were also included in our study to expand the scope of comparison.
The PFGE protocol was tailored from previous methods (
3) specifically
to facilitate lysis and macrorestriction digestion of spore-forming
Bacillus species. Cultures were enriched for vegetative cells
by diluting 5 µl of an overnight culture 1,000-fold in
5 ml of Trypticase soy broth medium and incubating for another
4 h at 37°C. Lysozyme (Sigma-Aldrich Co., St. Louis, MO;
250 mg/ml) was added to the diluted cells and molten agarose
during plug casting to boost the standard lysis step, and a
heat shock step (70°C for 30 min) was added after proteinase
K treatment. The
Bacillus genomic DNA was then digested with
the NotI restriction enzyme (New England Biolabs, Inc., Beverly,
MA), and the resulting macrorestriction DNA fragments were separated
by PFGE. Each gel was loaded with 22 sample lanes and 8 reference
lanes (Fig.
1). Thiourea (100 µM; Sigma-Aldrich Co.) was
added to the PFGE running buffer, which improved the PFGE image
quality of samples that contained residual nuclease activity
(
11). Three replicates were run on separate gels for each strain
analyzed.
BioNumerics (Applied Maths, Inc., Austin, TX) was used to analyze
the PFGE patterns. PFGE gel images were normalized using lambda
ladder markers (New England Biolabs) with the "snap-to-peak"
algorithm. Fragment sizes of the three sequenced strains (
4),
B. anthracis Ames (A2084),
B. cereus E33L, and
B. thuringiensis 97-27, were highly correlated (
R2 > 0.98) with those obtained
from virtual enzyme digestion using Kodon (Applied Maths, Inc.,
Austin, TX), validating the high quality of fragment separation
and normalization.
A dendrogram (Fig. 2) was produced from the PFGE profiles using band-based cluster analysis with the Dice coefficient and unweighted-pair group method with arithmetic mean algorithm. By using the "tolerance and optimization analysis" function, the optimal position tolerance was found at 0.72% and a 64.3% cluster cutoff value for the dendrogram was obtained. All B. anthracis isolates fell into the same cluster, with an overall similarity of 71.1%, and were well separated from B. cereus, B. thuringiensis, and B. atrophaeus isolates. Visual examination of the dendrogram showed that most of the B. anthracis strains shared the same profile, with similarities of over 94%, confirming previous reports of high homology (9).
Further improvements in the application of BioNumerics, along
with the adaptation of a more standard PulseNet experimental
design, e.g., the use of the PulseNet universal reference standard
(
Salmonella serotype Braenderup strain H9812) and fewer lanes
per gel, will enhance the quality of the analysis for use in
public health laboratories. However, the conditions and settings
used in this study resulted in the accurate grouping of all
B. anthracis isolates together and the clear differentiation
of
B. anthracis strains from
B. cereus,
B. thuringiensis, and
B. atrophaeus strains.
B. cereus and B. thuringiensis species were not segregated by this method, confirming reports using other genotyping methods (1, 5, 6), and individual isolates showed greater molecular diversity than did B. anthracis. The two B. atrophaeus strains were clustered with 100% similarity and clearly differentiated from those of the group 1 Bacilli. Principal component analysis was performed on the band-matching matrix, and the score plot verified the grouping of B. anthracis strains and their separation from other samples by the first principal component. Other commonly used macrorestriction enzymes, SmaI, XbaI, and XmaI, were also tested in this study; however, none of these performed as well as NotI did in the differentiation of B. anthracis from B. cereus and B. thuringiensis. Only two to three fragments were generated by SmaI digestion, and more than 50 bands were created by XmaI and XbaI digestions, indicating greater difficulties in using PFGE profiles for differentiation. None of the enzymes tested resulted in PGFE profiles that could be used to differentiate the nine phylogenetic groups within B. anthracis as found by the multiple-locus variable-number tandem-repeat analysis system (10); moreover, neither could the composite data sets formed by profiles from any combinations of the four tested enzymes.
Our results show that B. anthracis can be distinguished from its closest relatives by using the PFGE profiles from NotI digestion, after a few modifications to the PFGE procedure to facilitate complete lysis. PFGE may be a useful technique for the initial genotyping of Bacillus species in a research or public health laboratory, where the equipment and methods used in this technique are familiar, readily standardized, and widely available. PFGE is a whole-genome characterization method that is unbiased relative to other genotyping methods that are more rapid but are based on specific target sequences. Such a capability added to the public health laboratory toolbox of diagnostic methods would enhance the response to food-borne outbreaks and suspected acts of bioterrorism.

ACKNOWLEDGMENTS
We gratefully acknowledge Karen Hill for providing the strains
of
B. cereus and
B. thuringiensis used in this study.
Funding was provided by the Department of Homeland Security, National Bioforensics Analysis Center, and the Federal Bureau of Investigation. W. Zhong was supported by the National Flow Cytometry Resource (NCRR grant RR-01315).
The statements and conclusions herein are those of the authors and do not necessarily represent the views of the FBI.

FOOTNOTES
* Corresponding author. Mailing address: Bioscience Division, MS M888, Los Alamos National Laboratory, Los Alamos, NM 87545. Phone: (505) 667-3279. Fax: (505) 665-3024. E-mail:
blm{at}lanl.gov 
Published ahead of print on 30 March 2007. 
Present address: Department of Chemistry, University of California, Riverside, Riverside, CA 92521. 

REFERENCES
1 - Carlson, C. R., D. A. Caugant, and A.-B. Kolstø. 1994. Genotypic diversity among Bacillus cereus and Bacillus thuringiensis strains. Appl. Environ. Microbiol. 60:1719-1725.[Abstract/Free Full Text]
2 - Daffonchio, D., A. Cherif, and S. Borin. 2000. Homoduplex and heteroduplex polymorphisms of the amplified ribosomal 16S-23S internal transcribed spacers describe genetic relationships in the "Bacillus cereus group." Appl. Environ. Microbiol. 66:5460-5468.[Abstract/Free Full Text]
3 - Ferris, M. M., X. Yan, R. C. Habbersett, Y. Shou, C. L. Lemanski, J. H. Jett, T. M. Yoshida, and B. L. Marrone. 2004. Performance assessment of DNA fragment sizing by high-sensitivity flow cytometry and pulsed-field gel electrophoresis. J. Clin. Microbiol. 42:1965-1976.[Abstract/Free Full Text]
4 - Han, C. S., G. Xie, J. F. Challacombe, M. R. Altherr, S. S. Bhotika, D. Bruce, C. S. Campbell, M. L. Campbell, J. Chen, O. Chertkov, C. Cleland, M. Dimitrijevic, N. A. Doggett, J. J. Fawcett, T. Glavina, L. A. Goodwin, K. K. Hill, P. Hitchcock, P. J. Jackson, P. Keim, A. R. Kewalramani, J. Longmire, S. Lucas, S. Malfatti, K. McMurry, L. J. Meincke, M. Misra, B. L. Moseman, M. Mundt, A. C. Muck, R. T. Okinaka, B. Parson-Quintana, L. P. Reilly, P. Richardson, D. L. Robinson, E. Rubin, E. Saunders, R. Tapia, J. G. Tesmer, N. Thayer, L. S. Thompson, H. Tice, L. O. Ticknor, P. L. Wills, T. S. Brettin, and P. Gilna. 2006. Pathogenomic sequence analysis of Bacillus cereus and Bacillus thuringiensis isolates closely related to Bacillus anthracis. J. Bacteriol. 188:3382-3390.[Abstract/Free Full Text]
5 - Helgason, E., O. A. Okstad, D. A. Caugant, H. A. Johansen, A. Fouet, M. Mock, I. Hegna, and A.-B. Kolstø. 2000. Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensisone species on the basis of genetic evidence. Appl. Environ. Microbiol. 66:2627-2630.[Abstract/Free Full Text]
6 - Hill, K. K., L. O. Ticknor, R. T. Okinaka, M. Asay, H. Blair, K. A. Bliss, M. Laker, P. E. Pardington, A. P. Richardson, M. Tonks, D. J. Beecher, J. D. Kemp, A. B. Kolstø, A. C. L. Wong, P. Keim, and P. J. Jackson. 2004. Fluorescent amplified fragment length polymorphism analysis of Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis isolates. Appl. Environ. Microbiol. 70:1068-1080.[Abstract/Free Full Text]
7 - Hoffmaster, A. R., J. Ravel, D. A. Rasko, G. D. Chapman, M. D. Chute, C. K. Marston, B. K. De, C. T. Sacchi, C. Fitzgerald, L. W. Mayer, M. C. J. Maiden, F. G. Priest, M. Barker, L. Jiang, R. Z. Cer, J. Rilsone, S. N. Peterson, R. S. Weyant, D. R. Galloway, T. D. Read, T. Popovic, and C. M. Fraser. 2004. Identification of anthrax toxin genes in a Bacillus cereus associated with an illness resembling inhalation anthrax. Proc. Natl. Acad. Sci. USA 101:8449-8454.[Abstract/Free Full Text]
8 - Jernigan, J. A., D. S. Stephens, D. A. Ashford, C. Omenaca, M. S. Topiel, M. Galbraith, M. Tapper, T. L. Fisk, S. Zaki, T. Popovic, R. F. Meyer, C. P. Quinn, S. A. Harper, S. K. Fridkin, J. J. Sejvar, C. W. Shepard, M. McConnell, J. Guarner, W. J. Shieh, J. M. Malecki, J. L. Gerberding, J. M. Hughes, and B. A. Perkins. 2001. Bioterrorism-related inhalational anthrax: the first 10 cases reported in the United States. Emerg. Infect. Dis. 7:933-944.[Medline]
9 - Keim, P., A. Kalif, J. M. Schupp, K. K. Hill, S. E. Travis, K. Richmond, D. M. Adair, M. E. Hugh-Jones, C. R. Kuske, and P. Jackson. 1997. Molecular evolution and diversity in Bacillus anthracis as detected by amplified fragment length polymorphism markers. J. Bacteriol. 179:818-824.[Abstract/Free Full Text]
10 - Keim, P., L. B. Price, A. M. Klevytska, K. L. Smith, J. M. Schupp, R. Okinaka, P. J. Jackson, and M. E. Hugh-Jones. 2000. Multiple-locus variable-number tandem repeat analysis reveals genetic relationships within Bacillus anthracis. J. Bacteriol. 182:2928-2936.[Abstract/Free Full Text]
11 - Klaassen, C. H. W., H. A. van Hanneke, and A. M. Horrevorts. 2002. Molecular fingerprinting of Clostridium difficile isolates: pulsed-field gel electrophoresis versus amplified fragment length polymorphism. J. Clin. Microbiol. 40:101-104.[Abstract/Free Full Text]
12 - Radnedge, L., P. G. Agron, K. K. Hill, P. J. Jackson, L. O. Ticknor, P. Keim, and G. L. Andersen. 2003. Genome differences that distinguish Bacillus anthracis from Bacillus cereus and Bacillus thuringiensis. Appl. Environ. Microbiol. 69:2755-2764.[Abstract/Free Full Text]
13 - Rasko, D. A., M. R. Altherr, C. S. Han, and J. Ravel. 2005. Genomics of the Bacillus cereus group of organisms. FEMS Microbiol. Rev. 29:303-329.[CrossRef][Medline]
14 - Read, T. D., S. L. Salzberg, M. Pop, M. Shumway, L. Umayam, L. Jiang, E. Holtzapple, J. D. Busch, K. L. Smith, J. M. Schupp, D. Solomon, P. Keim, and C. M. Fraser. 2002. Comparative genome sequencing for discovery of novel polymorphisms in Bacillus anthracis. Science 296:2028-2032.[Abstract/Free Full Text]
15 - Swaminathan, B., T. J. Barrett, S. B. Hunter, and R. V. Tauxe. 2001. PulseNet: the molecular subtyping network for foodborne bacterial disease surveillance, United States. Emerg. Infect. Dis. 7:382-389.[Medline]
16 - Vilas-Boas, G., V. Sanchis, D. Lereclus, and M. V. Lemos. 2002. Genetic differentiation between sympatric populations of Bacillus cereus and Bacillus thuringiensis. Appl. Environ. Microbiol. 68:1414-1424.[Abstract/Free Full Text]
Applied and Environmental Microbiology, May 2007, p. 3446-3449, Vol. 73, No. 10
0099-2240/07/$08.00+0 doi:10.1128/AEM.02478-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.