Previous Article | Next Article ![]()
Journal of Clinical Microbiology, September 2006, p. 3352-3360, Vol. 44, No. 9
0095-1137/06/$08.00+0 doi:10.1128/JCM.00561-06
Patricia P. Wilkins,1
Swati B. Avashia,3,4,
Rahsaan Drumgoole,4
Charles H. Helma,5
Lawrence O. Ticknor,6
Richard T. Okinaka,5 and
Paul J. Jackson7
Epidemiologic Investigations Laboratory, Meningitis and Special Pathogens Branch, MS G34,1 Office of the Director,2 Epidemic Intelligence Service, MS E92, Centers for Disease Control and Prevention, 1600 Clifton Road, Atlanta, Georgia 30333,3 Infectious Disease Epidemiology and Surveillance Division, T-801, Texas Department of State Health Services, Austin, Texas 78756,4 Bioscience,5 Decision Applications Divisions, Los Alamos National Laboratory, Los Alamos, New Mexico 87545,6 Defense Biology Division, Lawrence Livermore National Laboratory, Livermore, California 945517
Received 15 March 2006/ Returned for modification 19 June 2006/ Accepted 10 July 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
B. cereus G9241, which was associated with severe pneumonia in a welder from Louisiana in 1994, has been well characterized, and its genome has been sequenced and analyzed (13). Genomic analysis and multilocus sequence typing (MLST) of this isolate revealed it to be closely related to B. anthracis. Several methods have shown that B. cereus isolates closely related to B. anthracis tend to be of clinical rather than environmental origin (9, 10, 11). However, it is very uncommon for B. cereus isolates, even those that are closely related to B. anthracis, to carry B. anthracis virulence plasmids (23). B. cereus G9241 carries an almost complete pXO1 plasmid, designated pBCXO1. This isolate also harbors a 218-kb circular plasmid (pBC218) and a cryptic bacteriophage (pBClin29) (13). Another unique feature of G9241 relative to other B. cereus isolates is the presence of a capsule. However, this capsule is not composed of D-glutamyl polypeptides and is not encoded by the B. anthracis cap genes normally located on the B. anthracis pXO2 plasmid. Instead, it has been hypothesized to be a polysaccharide and be encoded by a putative polysaccharide capsule biosynthetic operon located on pBC218 (13). While the presence of these plasmids in an isolate that causes severe disease similar to inhalation anthrax is intriguing, their roles, if any, in the virulence of the isolate or the presentation of disease has not yet been determined.
In October and November 2003, two fatal cases of B. cereus pneumonia occurred in metal workers (a welder and a muller operator) at different locations in Texas. A detailed report of these cases and the epidemiologic investigation is in preparation (S. B. Avashia, submitted for publication). In this report, we describe the initial molecular genetic characterization of two clinical and one environmental B. cereus isolate from the investigation of these two fatal cases. We also describe how these isolates are related to each other, to B. anthracis, and to other previously characterized B. cereus and B. thuringiensis isolates.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
DNA isolation and purification. Five milliliters of nutrient broth was inoculated with bacteria from a single colony of each isolate, and each culture was incubated overnight with shaking at 28°C. Four colonies from each isolate were used to produce four different cultures. Bacterial cells were harvested by centrifugation at 1,000 x g for 15 min. The resulting pellets were subjected to three freeze-thaw cycles. DNA was isolated from the disrupted cells using a QIAamp tissue kit (QIAGEN, Inc., Valencia, CA) following the protocol provided by the manufacturer. The quantity and quality of the isolated DNA were determined by gel electrophoresis of a small amount of the sample through a 1.0% agarose gel dissolved in a solution containing 10 mM Tris borate (pH 8.3) and 1 mM EDTA for 1 h at 80 V (26). Gels were stained for 20 min with a solution containing 1 µg ethidium bromide/ml, destained in distilled water, and visualized and photographed under UV light. Known concentrations of DNA molecular-weight standards were compared to the isolated bacterial DNA to determine the size and approximate concentration of the bacterial DNA.
PCR amplification of DNA. DNA oligomers that functioned as PCR primers to amplify representative B. anthracis gene fragments are listed in Table S1 in the supplemental material. Primers designed to amplify a portion of all bacterial 16S rRNA genes were used as positive PCR controls. Initial PCR contained 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl2, 0.001% (wt/vol) gelatin, 0.2 mM of each dNTP, 20 pmol of each primer, 2.5 U of AmpliTaq DNA polymerase, and approximately 1 ng template DNA in a 100-µl total reaction volume. PerkinElmer reagents (Foster City, CA) were used for all reactions. Template DNA was initially denatured by heating at 94°C for 2 min. This was followed by 35 cycles of denaturation at 94°C for 1 min, annealing at 55°C for 1 min, and primer extension at 72°C for 1 min. Incubation for 5 min at 72°C followed to complete the extension. PCR was conducted in a PerkinElmer GeneAmp PCR system 9600 thermal cycler. Control reactions containing B. anthracis 91-429C-2 DNA demonstrated whether each PCR primer set functioned as expected to amplify the appropriate B. anthracis plasmid or genomic amplicon.
Analysis of PCR amplicons. PCR amplicons were analyzed by electrophoresis through a 3.0% agarose gel dissolved in a solution containing 10 mM Tris borate (pH 8.3) and 1 mM EDTA. Electrophoresis was for 1 h at 80 V. Gels were stained for 20 min with a solution containing 1 µg of ethidium bromide/ml, destained in distilled water, and then visualized and photographed under UV light. Images were captured electronically using a Stratagene Eagle Eye II still video system (Stratagene).
Protective antigen (PA) gene sequencing and multilocus sequence typing. The pagA gene was amplified and sequenced as described previously (12, 24). However, additional primers amplifying a 1,076-nucleotide fragment between pXO1 nucleotide 134,396 and 135,472 (see Table S1 in the supplemental material) were used in an attempt to amplify the 5' end of pagA in B. cereus 03BB102 because this gene failed to amplify using previously described primers. MLST was based on the analysis of seven housekeeping gene partial sequences (glpF, gmk, ilvD, pta, pur, pycA, and tpi) and is described online at http://pubmlst.org/bcereus (25). Phylogenetic and molecular evolutionary analyses were conducted using MEGA, version 2.1 (18).
AFLP analysis of DNA samples. Amplified fragment length polymorphism (AFLP) analysis was accomplished as previously described (11, 14, 15, 26). Briefly, 100 ng of DNA was digested with the restriction endonucleases EcoRI and MseI and the resulting fragments were ligated to double-stranded adapters. The digested and ligated DNA was then amplified by PCR using EcoRI and MseI +0/+0 primers (5'-GTAGACTGCGTACCAATTC-3' and 5'-GACGATGAGTCCTGAGTAA-3', respectively). The + 0/+0 PCR product was analyzed by agarose gel electrophoresis to determine the size range of amplified fragments. Three microliters of each product was used as a template in subsequent selective amplifications using the +1/+1 primer combination of 6-carboxyfluorescein-labeled EcoRI-C (5'-GTAGACTGCGTACCAATTCC-3') and MseI-G (5'-GACGATGAGTCCTGAGTAAG-3'). Selective amplifications were performed in 20-µl reaction mixtures. The resulting products (0.5 to 1.0 µl) were mixed with a solution containing a mixture of DNA size standards (GeneScan-500 [Applied Biosystems Inc., Foster City, CA] and MapMarker 500 [BioVentures, Inc., Murfreesboro, TN]) labeled with N,N,N,N-tetramethyl-6-carboxyrhodamine. Following a 2-min heat denaturation at 90°C, the reactions were loaded onto 5% Long Ranger DNA sequencing gels (Cambrex Bio Science, Rockland, ME) and visualized on an ABI 377 automated fluorescent sequencer. To make sure that the DNA assayed in these experiments was representative of each of the isolates, four individual colonies from each isolate were inoculated into medium and DNA purified. The DNA preparations were then analyzed using AFLP. Forty DNA fragments generated from the AFLP experiments were used to represent each of the DNA preparations. These DNA preparations from individual colonies of the same isolate were indistinguishable from each other. Therefore, only one DNA preparation was used in the other experiments. Each set of reactions also contained an AFLP reaction using B. anthracis Vollum DNA as an experimental control. The inclusion of such a reaction in each analysis set allowed a comparison of results from earlier archived analysis sets run at different times or on different gels. GeneScan analysis software (Applied Biosystems, Inc.) was used to determine the lengths of the sample fragments by comparisons to the DNA fragment length size standards included with each sample.
To minimize gel electrophoresis artifacts, each labeling reaction was run in triplicate. Samples were loaded on three different gels in a random order. AFLP data analysis was performed as previously described by Ticknor et al. (26). Sample fragments of between 100 and 500 bp, with fluorescence above 50 arbitrary units in all three runs on the ABI sequencer, were used in the analysis.
Similarities among samples were determined by the Jaccard coefficient. The 40 tallest peaks for each sample fingerprint were used to calculate the Jaccard coefficient among samples. Dendrograms were produced by using the similarity matrix of Jaccard coefficients and the unweighted pair-group mean average method (F. J. Rohlf, NTSYS-PC numerical taxonomy and multivariate analysis system, version 1.8; Exeter Software, Setauket, NY). All statistical data manipulations were performed using codes developed in S-Plus (Data Analysis Products Division, MathSoft, Seattle, WA).
| RESULTS |
|---|
|
|
|---|
The two-component B. anthracis-specific DFA (6) was performed on all three isolates and produced positive results for the cell wall antigen but was negative for the B. anthracis capsule antigen (poly-
-D-glutamic acid) when cells were incubated under several conditions, including conditions used for B. anthracis capsule expression. B. cereus clinical isolates are sometimes positive for the cell wall antigen. B. cereus G9241, isolated from a Louisiana welder (13), also tested positive for this cell wall antigen. However, unlike the results obtained for B. cereus G9241, India ink staining of the two Texas clinical isolates did not initially reveal any type of capsule when grown on sheep blood agar at 37°C in ambient atmosphere due to a very small percentage of the cells being encapsulated (Table 1). Adjusting the medium and temperature resulted in increased capsule formation, but not to the level exhibited by B. cereus G9241. A capsule was not detected on the environmental isolate 03BB108.
Plasmid PCR analysis. PCR analysis was performed using primers designed to amplify sequences from pXO1, pXO2, and pBC218. Table S1 in the supplemental material shows the sequences and targets for a battery of PCR primers used to assay for the presence of different pXO1, pXO2, and pBC218 open reading frames (ORFs). Table 2 shows the results of PCR assays to detect pXO1 and pXO2 sequences. All primers generated the expected amplicon when the control strain (B. anthracis 91-429C-2) DNA was used as the template.
|
B. cereus 03BB102 and 03BB108, isolated from the Texas welder and his worksite, respectively, did not contain sequences of sufficient homology to generate PCR amplicons from many of the pXO1 ORFs tested (Table 2). The amplification profile of B. cereus 03BB108 was a subset of that for B. cereus 03BB102. This latter isolate contained pXO1 ORFs that were not present in B. cereus 03BB108. A possible explanation of this result is a loss of DNA from B. cereus 03BB102 to produce the profile seen in B. cereus 03BB108, but the pattern of ORFs missing from B. cereus 03BB108 relative to B. cereus 03BB102 is complex relative to the order of these ORFs in the B. anthracis pXO1 plasmid. While B. cereus 03BB102 tested positive for several genes located within the pXO1 pathogenicity island (22), including the three toxin genes (pagA, lef, and cya), the environmental isolate (B. cereus 03BB108) was missing the majority of the pathogenicity island ORFs.
Analyses to detect pXO2-specific ORFs revealed that, while little if any of this plasmid was present in either B. cereus 03BB102 or 03BB108, the capA, capC, and capB genes, which are required for biosynthesis of the B. anthracis capsule, were detected in both isolates (Table 2). Genomic analysis of B. cereus G9241 identified a 218-kb plasmid, pBC218, which includes a polysaccharide capsule operon that was hypothesized to encode the observed capsule for this isolate (13). PCR analysis using primers that are complementary to pBC218 (Table 3) showed that B. cereus 03BB87 has an identical PCR profile to B. cereus G9241 for 32 different pBC218 markers scattered around the plasmid, strongly suggesting the presence of pBC218 or a very similar plasmid in B. cereus 03BB87. In particular, B. cereus 03BB87 possessed the putative polysaccharide polymerase gene (ORF pBC218-0073) thought to participate in the synthesis of a polysaccharide capsule, while B. cereus 03BB102 and 03BB108 did not possess this gene.
|
pagA sequencing. The B. anthracis protective antigen gene pagA was detected by PCR in B. cereus 03BB87 and 03BB102. To determine similarities between the B. anthracis pagA gene, G9241 pagA, and the pagA sequences from these isolates, we attempted to amplify and sequence pagA from each of the Texas clinical isolates. The pagA sequence in B. cereus 03BB87 was identical to the pagA sequence from B. cereus G9241 (13) and has been assigned the designation of PA genotype 9 (Table 1). Attempts to amplify the extreme 5' end of pagA from B. cereus 03BB102 were not successful, despite using additional primers. We were therefore able to generate only a partial pagA sequence that was missing the first 15 bases relative to the published B. anthracis pagA sequence (see reference 24 and other GenBank entries for this gene). A comparison of this partial sequence to the pagA gene sequence amplified from B. cereus 03BB87 and G9241 revealed two additional mutations at positions 1862 and 1898, which further distinguishes them from the other two isolates and previously sequenced B. anthracis isolates. These two mutations result in serine-to-asparagine and lysine-to-threonine changes, respectively, and are located in the N-terminal 20-kDa PA fragment (PA20) that is cleaved from PA63 (63-kDa C-terminal fragment) by a class of furin proteases (21).
Molecular subtyping. All three Texas isolates and the Louisiana isolate (G9241) were analyzed using AFLP to understand their genetic relationship to one another, to B. anthracis, and to other previously characterized B. cereus and B. thuringiensis isolates. AFLP is a multilocus sampling method that relies on the generation of a fingerprint or signature that is created by restriction enzyme digestion of the entire genome and any extrachromosomal material that might be present. After digestion with restriction enzymes and further selective PCR amplification, the AFLP fingerprint for an isolate was represented by approximately 40 fragments between 100 and 500 bp in length that were present in all three replicates of each isolate. These fingerprint fragment sizes were then used to generate a phylogenetic tree that illustrates the genetic relationships among the isolates analyzed. The fingerprints generated from the three Texas isolates were compared to a fingerprint generated from total DNA isolated from the B. anthracis Vollum strain that is included in every AFLP analysis as a control. They were also compared to fingerprints generated for a battery of B. cereus and B. thuringiensis isolates that are known to be closely related to B. anthracis. Figure 1 shows that, relative to the majority of B. cereus and B. thuringiensis isolates studied, all three Texas isolates (B. cereus 03BB87, 03BB102, and 03BB108) and the Louisiana isolate (B. cereus G9241) are very closely related to B. anthracis (11, 26). However, they are clearly distinct from this pathogen. They mapped to a branch of the phylogenetic tree that includes many of the known pathogenic and toxigenic B. cereus and B. thuringiensis isolates so far characterized (11), while most B. cereus and B. thuringiensis isolates map to other branches of the tree. AFLP analysis could not differentiate B. cereus 03BB87, isolated from the Texas muller operator, from B. cereus G9241, isolated from the Louisiana welder in 1994. In contrast, B. cereus 03BB102 is phylogenetically significantly different from these two clinical isolates. Its AFLP fragment profile was similar to, but distinct from, the profile generated for B. cereus 03BB108, the isolate collected from the working environment of the B. cereus 03BB102-infected victim (Texas welder). Figure 2 shows a direct comparison of a major portion of the AFLP DNA fragment profiles for the two clinical isolates and that of G9241. A direct comparison of the AFLP profiles of B. cereus G9241, isolated in 1994 from a Louisiana welder, and B. cereus 03BB87, isolated in 2003 from a Texas muller operator (Fig 2A), showed that, within the resolution of the method, the isolates were virtually identical. In contrast, Fig. 2B shows the B. cereus 03BB102 profile. There were distinct differences between this profile and the two profiles shown in Fig. 2A. The AFLP profile of B. cereus 03BB102 from the Texas welder, was similar to but distinct from the profile for B. cereus 03BB108, which was collected from the welder's work environment (Fig. 1).
|
|
|
| DISCUSSION |
|---|
|
|
|---|
B. cereus 03BB102 and 03BB108, from the Texas welder and his workplace, respectively, are distinctly different from B. cereus 03BB87, collected from the Texas muller operator, based on AFLP and MLST analyses (Fig. 1, 2, and 3). These two isolates are similar to, but distinct from, one another. AFLP analysis places them in a small subbranch of the phylogenetic tree (Fig. 1). MLST results (Fig. 3) are consistent with this finding, showing the two isolates to be similar but distinct. PCR analysis (Table 2) also reveals differences in their pXO1 gene profiles. The isolates also manifest different colony morphologies. It is clear from these results that B. cereus 03BB108 was not derived from B. cereus 03BB102 by a simple loss of plasmid DNA material since the two isolates are genetically distinct from one another beyond this level. A comparison of AFLP and MLST profiles for B. cereus 03BB102 and 03BB108 to those generated from B. anthracis and other B. cereus and B. thuringiensis isolates (Fig. 1 and 3) demonstrates that, while these isolates may contain fewer sequences homologous to B. anthracis pXO1 and pXO2 sequences, they are phylogenetically more similar to B. anthracis than are B. cereus 03BB87 and G9241. They are also closely related to two other pathogenic Bacillus isolates (B. thuringiensis 97-27 and B. cereus E33L) and two toxigenic B. cereus isolates (D17 and 3a) (Fig. 1).
The presence of the capA, B, and C genes in B. cereus 03BB102 and 03BB108 suggests that these genes may play a role in virulence as they do in B. anthracis. However, we did not observe the expression of the poly-
-D-glutamic acid capsule in these isolates under in vitro conditions where it is normally expressed in B. anthracis and we do not yet know whether the entire genes are present or under what conditions they might be expressed. The presence of the cap gene sequences in the absence of any other detectable pXO2 sequences is unique relative to other isolates so far characterized. Most of the differences in the detection of pXO1 genes between these isolates are within the pathogenicity island (22).
Despite the fact that B. anthracis toxin genes are found in several of the B. cereus isolates associated with severe pneumonia in metal workers, we do not know whether these genes are appropriately expressed or whether any expressed toxins are functional or play a role in virulence. At least one such isolate, described by Miller et al. (20), did not contain the pagA gene yet still caused a fatal infection in a previously healthy individual (13). However, relative to the general population, professional exposure to a variety of environmental factors may place welders at additional risk of acquiring such infections. Reports suggest welders experience higher frequencies of pulmonary infections with increased severity and duration (1, 2). The majority of data suggest that while welders may not be at increased risk for infections overall, certain conditions and practices may impact susceptibility to pneumonias (1, 4, 7). It has also been suggested that increased exposure to these types of B. cereus isolates and other respiratory pathogens may play a role in increased pneumonia in foundry workers, including welders (5, 17).
All of the isolates clustered in a branch of the AFLP-based phylogenetic tree occupied by B. anthracis and several other pathogenic B. cereus and B. thuringiensis isolates so far characterized (Fig. 1). It is also intriguing that they were collected in a relatively small geographic area of the United States (Texas and Louisiana) and that this area coincides with regions where anthrax occurs naturally in herbivores. If these isolates acquired these plasmids or plasmid sequences via horizontal transfer from B. anthracis in the environment, it is tempting to speculate that such isolates may be restricted to or at least be more common in areas where anthrax is endemic.
B. cereus is ubiquitous in the environment, which can make it difficult to link clinical cases to their environmental sources. Although the source of the clinical isolate was not identified for either case, the isolation of 03BB108 from the welder's worksite illustrates that closely related B. cereus isolates were present. Increased sampling and analysis of isolates may have allowed for the detection of isolates in the environment that matched the clinical isolate. It is also possible that the source of infection was present in material to which both victims were occupationally exposed because they shared similar work environments. This scenario would be more likely if the same strain infected both Texas victims. However, B. cereus 03BB87, isolated from the Texas muller operator, is very similar or identical to B. cereus G9241, isolated from the Louisiana welder 10 years earlier. Therefore, the possibility of a common source is very intriguing, even though these cases were separated temporally by 10 years.
Taken together with previous reports, these results clearly demonstrate the need for increased awareness of the potential for B. cereus to cause serious systemic infections, even in patients who appear to be otherwise healthy. They also demonstrate the need for a better understanding of environmental and occupational risk factors that may increase the susceptibility of individuals to infection by B. cereus.
| ACKNOWLEDGMENTS |
|---|
Funding for a portion of this work was provided by the United States Department of Homeland Security.
| FOOTNOTES |
|---|
Supplemental material for this article may be found at http://jcm.asm.org/. ![]()
Present address: Department of Microbiology and Immunology, Emory University School of Medicine, 3004 Rollins Research Center, Atlanta, GA 30322. ![]()
Present address: Internal Medicine/Pediatrics, Seton Kozmetsky Clinic, 3706 South First Street Austin, TX 78704. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|