Previous Article | Next Article ![]()
Journal of Clinical Microbiology, December 2004, p. 5825-5831, Vol. 42, No. 12
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.12.5825-5831.2004
and
David Norwood1*
United States Army Medical Research Institute of Infectious Diseases, Fort Detrick, Maryland,1 Clinical Research Management, Hinckley, Ohio2
Received 27 April 2004/ Returned for modification 4 June 2004/ Accepted 9 August 2004
|
|
|---|
|
|
|---|
Past attempts to develop a chromosomal real-time PCR assay have failed due to the close genetic relationship of Bacillus species. B. anthracis, Bacillus cereus, and Bacillus thuringiensis have very little variability and are genetically indistinguishable by multilocus enzyme electrophoresis (10). Recent work by repetitive PCR has shown that the previously listed species of Bacillus, as well as Bacillus mycoides, Bacillus pseudomycoides, and Bacillus weihenstephanensis, do have some genetic differences (3). Real-time PCR assays based on the chromosomal rpoB and gyrA genes of B. anthracis have been developed (5, 13, 25). However, these assays are based on single-nucleotide differences between B. anthracis and other Bacillus species. In our hands, the performance of assays based on point mutations has been notoriously dependent on assay conditions and cycling parameters, and small alterations in these conditions can result in the loss of specificity, especially with TaqMan-based probes.
Thus, a real-time PCR assay based on a unique chromosomal nucleotide sequence in B. anthracis is needed to complement other available assays for identification of the organism. Plasmids can be transmitted among bacterial species; therefore, a chromosomal marker to identify B. anthracis as the organism harboring the plasmids is important in a time of bioengineering and terrorist threats. While the presence of both pXO1 and pXO2 is needed to give B. anthracis its virulence, it is conceivable that these plasmids could be passed to its genetic neighbors, with unknown implications.
Bacterial subtraction has previously been used to find unique sequences in bacterial species (6, 8, 24, 26). In this study, nucleic acid isolated from B. cereus ATCC 21769 and genomic DNA from plasmid-cured B. anthracis strain Delta ANR-SWS (dANR) was used in a genomic subtraction-hybridization to identify possible chromosomal sequences unique to B. anthracis that could be used to develop a real-time PCR assay.
|
|
|---|
Bacterial subtraction. A PCR-Select Bacterial Genome Subtraction kit (Clontech, Palo Alto, Calif.) was used to identify sequences unique to B. anthracis. Strain BACI177 was chosen for use as the driver DNA, and strain dANR was chosen for use as the tester. The genomic Escherichia coli provided in the kit was used as a control. The majority of the subtraction-hybridization was performed according to the manufacturer's specifications, but the modifications described below were used. Restriction enzyme digestion was performed according to the instructions provided with the kit, and the digest was analyzed for efficiency by gel electrophoresis, as described above. The digested nucleic acid was purified, and adaptor ligations were performed with the tester and the control DNA. Ligation efficiency analysis was performed by PCR with the following components: 35.5 µl of molecular biology-grade (MBG) H2O, 2 µl of a 10 mM deoxynucleoside triphosphate (dNTP) mix (Roche, Indianapolis, Ind.), 5 µl of 10x PCR buffer, 4 µl of 25 mM MgCl2, 0.5 µl of AmpliTaq Gold (Applied Biosystems, Foster City, Calif.), 1 µl of each primer, and 1 µl of ligated DNA per 50-µl reaction mixture. The primers used with the tester DNA were 704F (5'-AGATTTTCCGACGGCAGGTT-3') at a concentration of 50 pmol/µl and 829R (5'-TTTCAATCAATCGCGCCTTATT-3') at a concentration of 50 pmol/µl; these primer sequences span a region of the gyrA gene in B. anthracis. One tester DNA reaction mixture incorporated primer 704F with PCR primer 1 at a concentration of 0.2 µM/µl, while primer 704F and primer 829R were used in the other tester DNA reactions. The PCR conditions were 72°C for 2 min and 95°C for 10 min, followed by 35 cycles of 95°C for 2 min, 60°C for 1 min, and 72°C for 1 min, with a final extension at 72°C for 10 min. The PCRs were performed with an MJ Research (Waltham, Mass.) PTC-100 thermal cycler. The DNA products were analyzed by gel electrophoresis on a 2% gel run at 150 V for 30 min. Following the ligation efficiency analysis, the two hybridization steps were performed at an incubation temperature of 62°C for both the tester DNA and the driver DNA due to their lower G+C contents. The PCRs used to amplify tester nucleic acid followed the manufacturer's instructions, with the exception that the nested PCR was run for 17 cycles, and were performed on an MJ Research PTC-100 thermal cycler. The DNA products were analyzed by gel electrophoresis, as described directly above.
Bacterial cloning. The nested PCR products subtracted from strain dANR were cloned into the pCR2.1-TOPO vector by using the TOPO TA Cloning kit (Invitrogen, Carlsbad, Calif.), according to the manufacturer's instructions for chemical transformation. White colonies were selected and boiled in 50 µl of H2O for 10 min to break up the cells and release the nucleic acid. PCR was performed with the following components per 100-µl reaction mixture: 71 µl of MBG H2O, 4 µl of a 10 mM dNTP mixture, 10 µl of 10x PCR buffer, 8 µl of 25 mM MgCl2, 1 µl of AmpliTaq Gold, 2 µl of M13 forward primer (50 pmol/µl; Invitrogen), 2 µl of M13 reverse primer (50 pmol/µl; Invitrogen), and 2 µl of supernatant from the boiled preparation. PCR conditions were 72°C for 2 min and 95°C for 10 min, followed by 30 cycles of 95°C for 2 min, 55°C for 1 min, and 72°C for 1 min, with a final extension at 72°C for 10 min. The PCRs were performed with an MJ Research PTC-100 thermal cycler. PCR products were analyzed by gel electrophoresis, as described above. Clones with positive inserts were selected for sequencing and further analysis.
Bacterial sequencing. PCR products from the selected clones were purified with the Montage PCR96 cleanup kit (Millipore, Billerica, Mass.). Sequencing reactions were carried out with 8 µl of Big Dye (Applied Biosystems), 1 µl of M13 forward primer (50 pmol/µl) or 1 µl of M13 reverse primer (50 pmol/µl; Invitrogen), 1 µl of purified PCR product, and 10 µl of MBG H2O. The sequencing reactions were cycled on an MJ Research PTC-100 thermal cycler at 85°C for 30 s, followed by 25 cycles of 96°C for 10 s, 50°C for 5 s, and 60°C for 4 min, with a final extension at 60°C for 10 min. Sequencing was performed on an ABI 3100 sequencer (Applied Biosystems), and the sequences were aligned by using DNASTAR software (DNASTAR, Inc., Madison, Wis.). Searches with the BLAST program were conducted for all sequences to eliminate regions not unique to B. anthracis and those with possible cross-reactivities with other organisms. A single clone (clone B26) was picked for further analysis.
Primer and probe design and testing. Prospective primers and probes were generated with Primer Express software (version 2.0; Applied Biosystems) and the imported sequence from clone B26. Primers were further analyzed for palindromes, hairpins, dimer formation, and annealing temperature with NetPrimer software (BioSoft International, Palo Alto, Calif.).
Primer testing by standard PCR. Two forward primers and two reverse primers were assessed by standard PCR on an MJ Research PTC-100 thermal cycler with the following components per 20-µl reaction mixture: 10.44 µl of MBG H2O, 2 µl of 10x 2 mM dNTPs (Idaho Technology, Salt Lake City, Utah), 2 µl of 10x PCR buffer containing 50 mM MgCl2 (Idaho Technology), 0.2 µl of forward primer (50 pmol/µl), 0.2 µl of reverse primer (50 pmol/µl), 0.16 µl of of Platinum Taq DNA polymerase (5 U/µl; Invitrogen), and 5 µl of Bacillus species DNA at either 1 ng/µl or 100 pg/µl. The reaction mixtures were cycled at 94°C for 2 min, followed by 34 cycles of 94°C for 30 s, 60°C for 30 min, and 72°C for 1 min, with a final extension at 72°C for 10 min, on an MJ Research PTC-100 thermal cycler. Reactions were analyzed by gel electrophoresis as described above. Primer F41 (5'-TGGCGGAAAAGCTAATATAGTAAAGTA-3') and primer R146 (5'-CCACATATCGAATCTCCTGTCTAAAA-3') were selected for further testing. Forty-two B. anthracis strains and 53 Bacillus species were tested by PCR under the conditions defined above. The reaction products were analyzed by gel electrophoresis, as described above.
Probe testing by real-time PCR. Two Taqman minor groove binder (MGB) probes were tested by real-time PCR on the Ruggedized Advanced Pathogen Identification Device (RAPID; Idaho Technology) with primers F41 and R146. The background fluorescence of the probes was obtained by diluting the 100 µM MGB Taqman probes to 1 µM so that the fluorescence was between 10 and 20 with a gain setting of 16. The probes were tested against 16 B. anthracis isolates and 13 B. cereus isolates. The master mixture contained 8.44 µl of MBG H2O, 2 µl of 10x 2 mM dNTPs, 2 µl of 10x PCR buffer containing 50 mM MgCl2, 0.2 µl of primer F41 (50 pmol/µl), 0.2 µl of primer R146 (50 pmol/µl), 2 µl of 1 µM probe, 0.16 µl of Platinum Taq (5 U/µl), and 5 µl of DNA (1 ng/µl) per 20-µl reaction mixture. The PCR cycling conditions for all runs on the RAPID were 94°C for 2 min, followed by 45 cycles of 94°C for 20 s and 60°C for 20 s. Probe BALB26588-MGB (P88; 5'-ACTTCTAAAAAGCAGATAGAAAT-3') was chosen for optimization of the RAPID. Five concentrations of MgCl2 (3, 4, 5, 6, and 7 mM) and 10 concentrations of primers F41 and R146 (0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0 µM) were tested in triplicate with 1 pg of B. anthracis Ames DNA. The combination that exhibited the earliest crossing threshold (CT) and highest end fluorescence (Ef) was chosen as the optimal condition, which consisted of 5 mM MgCl2 and 0.7 µM each primer. The optimized master mixture contained 8.28 µl of H2O, 2 µl of 10x 2 mM dNTPs, 2 µl of 10x PCR buffer containing 50 mM MgCl2, 0.28 µl of primer F41 (50 pmol/µl), 0.28 µl of primer R146 (50 pmol/µl), 2 µl of probe P88 (1 mM), and 0.16 µl of Platinum Taq (5 U/µl) per 20-µl reaction mixture. Tenfold dilutions of Ames strain DNA from 4 pg/µl to 0.4 fg/µl were made so that 2.5 µl of each dilution along with 2.5 µl of human DNA (5 ng) or 2.5 µl of H2O could be added to 15 µl of the optimized master mixture.
Cross-reactivity and other tests under real-time optimized conditions. A panel of 77 organisms was run by real-time PCR with the optimized conditions to test for cross-reactivity. The master mixture contained 5 µl of DNA (20 pg/µl) isolated from each organism. Sensitivity and specificity testing was also accomplished as described above by using 45 B. anthracis strains and 62 additional Bacillus species.
Limit of detection and efficacy of the assay. The limit of detection of the assay was established by running 60 samples of strain Ames DNA at a single concentration with two negative controls. The lowest concentration that was positive for at least 57 of the 60 samples (95%) was determined to be the limit of detection. The efficacy, or efficiency, of the PCR assay was determined by testing 10-fold serial dilutions from 100 pg to 10 fg in triplicate. The crossing point of each dilution was determined by the second derivative method in the LightCycler software (version 3.3; Roche). Amplification efficacy (E) was calculated on the basis of the slope of the crossing point versus the input DNA concentration by the equation E = 10(1/slope) (23).
|
|
|---|
![]() View larger version (132K): [in a new window] |
FIG. 1. DNA products following nested PCR. Lanes 1 and 8, DNA PCR ladder (Gene Choice, Frederick, Md.); lane 2, subtracted B. anthracis dANR; lane 3, unsubtracted B. anthracis dANR; lane 4, subtracted E. coli; lane 5, unsubtracted E. coli; lane 6, subtracted control E. coli; lane 7, empty. The values at the sides of the gel are in base pairs.
|
Probe testing by real-time PCR. After the primers were tested and selected by standard PCR, two probes previously designed with the primers were tested against 16 B. anthracis and 13 B. cereus strains with the RAPID. The specificity and sensitivity were each 100% for each probe tested. Both probes generated similar CT values for the samples tested, but P88 generated a higher Ef for the samples tested, so it was chosen for optimization. After optimization of the MgCl2 and primer concentrations, serial 10-fold dilutions made with strain Ames DNA gave positive results to 1 fg in triplicate assays (data not shown). When the reactions were spiked with 5 ng of human DNA, positive results to 1 fg in triplicate assays were still obtained, with no changes in CT or Ef observed.
Cross-reactivity panel. The optimized assay was tested with a cross-reactivity panel of 76 organisms (Table 1), including 11 B. anthracis isolates. All 11 isolates were detected, regardless of the plasmid profile, by the chromosomal marker test, while the remaining 66 organisms failed to produce a signal. In addition, sensitivity and specificity were evaluated by running 45 B. anthracis isolates and 62 Bacillus species by the optimized assay (these isolates included the original 16 B. anthracis isolates and 13 B. cereus isolates tested), and the results are shown in Table 2. All 45 B. anthracis strains were found to be positive, and all of the other 62 Bacillus species were negative, including BACI179, providing 100% sensitivity and specificity.
|
View this table: [in a new window] |
TABLE 1. Organisms used to test cross-reactivity by the real-time chromosomal PCR assay
|
|
View this table: [in a new window] |
TABLE 2. Bacillus species used for sensitivity and specificity testing of the real-time chromosomal PCR assay
|
![]() View larger version (50K): [in a new window] |
FIG. 2. Real-time PCR results for determination of efficacy. NTC, no template control; PGAMES, picograms of strain Ames; FGAMES, femtograms of strain Ames.
|
|
|
|---|
The importance of designing a chromosomal assay for B. anthracis has arisen due to the discovery of isolates cured of plasmids pX01 and pX02 (27, 31). Identification of avirulent B. anthracis strains could be a potential indicator of virulent strains in the environment (31). Real-time PCR based on chromosomal markers could provide a powerful tool for identifying B. anthracis strains either containing or lacking virulence plasmids. To date, the majority of B. anthracis targets for real-time PCR assays have been developed from genes located on the pXO1 and pXO2 plasmids (2, 7, 20). These assays are specific because the targets code for the toxins and capsule proteins essential for virulence. However, complete identification of B. anthracis should involve the detection of both plasmids and a chromosomal marker due to the possibility of the transmission, transfer, or loss of one or both plasmids between Bacillus species seen both in vivo and in vitro (18, 27, 29, 31). Bioengineering may also allow the transfer of virulence plasmids to other species, including bioengineered organisms (12), exemplifying the need to detect these agents made by humans.
The development of chromosomal real-time PCR assays has been slow due to the genetic similarities between B. anthracis and its neighbors (10). To date there are very few chromosomal markers to distinguish B. anthracis from its close relatives. Analysis of the 16S rRNA genes from B. anthracis, B. cereus, B. mycoides, and B. thuringiensis revealed that they have almost identical sequences (9); and some argue that due to the similarity of the B. anthracis, B. thuringiensis, and B. cereus genomes, they should be considered to belong to the same species (10). Some chromosomal markers are available, such as the Ba813 fragment (21) and the vrrA variable number tandem repeat (1, 15, 17); however, these markers have failed because they are not specific to B. anthracis (4, 12), and assays incorporating these targets require extra steps after the PCR, making them impractical for the rapid identification of B. anthracis (30). Recently, the rpoB and gyrA genes, located on the B. anthracis chromosome, have been targeted for the development of a real-time PCR assay (13, 25). These assays are based on theoretical differences in melting temperatures due to single-nucleotide differences, but in practice, false-positive results with other Bacillus species have occurred (7, 25). In fact, Bacillus sp. strain Ba813+ #11, which was found to be positive by the published rpoB assay (25), was negative by the chromosomal assay developed here. The present assay has the potential for increased specificity because it is not dependent on probable point mutations or theoretical differences in melting temperatures, and all primers and probes are based on a sequence unique to B. anthracis.
The assay developed in this study is based on a chromosomal sequence unique to B. anthracis. This region was found by subtractive hybridization incorporating a plasmid-cured B. anthracis (tester) strain and a B. cereus (driver) isolate. Previous reports have shown success with the use of this method to find sequences in the tester strain's genome that are absent in the driver (8, 24), and recently, it has been used to find unique targets in B. anthracis (26). The forward and reverse primers as well as the reverse complements from the previous report (26) were checked for alignment with our sequence, without success. By the use of subtractive hybridization, a sequence unique to B. anthracis was found and was used to develop a chromosomal real-time PCR assay that had 100% sensitivity in tests with all B. anthracis strains studied and 100% specificity when it was tested with other Bacillus species. In addition, a cross-reactivity panel of 76 organisms that included 27 Bacillus strains was used to test the sensitivity and the specificity of the assay. The organisms in this panel were chosen because of their genetic proximity to B. anthracis and their possible similar environmental and clinical niches with B. anthracis.
The recent release of the B. anthracis genome allowed the comparison of the B26 clone sequence with this new sequence information by use of the BLAST program. Forward primer F41 and probe P88 produced significant alignments with B. anthracis strain Ames and strain Sterne, having the highest scores compared to any other alignment and E values of 7 x 104. Reverse primer R146 also produced the highest scores with B. anthracis strain Ames and strain Sterne, with E values of 2 x 105; however, an alignment to B. cereus ATCC 10987 with an E value of 0.005 also occurred. This B. cereus strain has been tested and was not detected by the assay. A search of the assay sequence and surrounding nucleotides (527 nucleotides in total) against all protein sequences with the BLAST program produced a high score of 191, with an E value of 5 x 1048 for the homology of B. anthracis (GenBank accession number NP_656566.1) with an abhydrolase (alpha/beta-hydrolase) fold. The 106-bp product produced in this investigation is just outside the open reading frame of this protein. Other significant alignments included a prolyl aminopeptidase of an Oceanobacillus iheyensis strain with a score of 115 and an E value of 2 x 1029 and a proline iminopeptidase of B. cereus ATCC 14579 with a score of 94 and an E value of 4 x 1021.
Further testing with other B. anthracis strains from other sources must be performed to ensure that this sequence is conserved between additional B. anthracis strains. In addition, the phenotypic necessity of this genetic region has yet to be determined. While amplified fragment length polymorphism analysis (16) has shown little variability among B. anthracis isolates, despite the diversity that exists, previous work with multiple-locus variable number tandem repeat analysis has classified B. anthracis into 89 genotypes (17). In our study, genotype information was available for 17 of the 45 B. anthracis strains studied. These 17 B. anthracis strains represented 8 of the 89 genotypes. Recent sequencing of a B. anthracis isolate from a victim of the 2001 anthrax attack has revealed single-nucleotide polymorphisms, insertions, deletions, and changes in tandem repeats when the strain's sequence is compared to those of reference B. anthracis strains (28). It is conceivable that the region targeted in this study may not be conserved. Genetic variation and molecular diversity in sequences could potentially decrease the sensitivity and the specificity of this particular assay. However, molecular identification of B. anthracis is a polyphasic process, requiring assays for pXO1, pXO2, and the chromosome. Furthermore, an integrated diagnostic approach that uses multiple genetic markers decreases the risk of technological surprise. The assay developed in this study could be an important component of an integrated diagnostic system for the detection of B. anthracis.
The present findings demonstrate that the assay developed is sensitive and specific for the identification of the B. anthracis chromosome. The assay did not cross-react with any genetic neighbors available for the study, and all strains of B. anthracis were identified. The assay was independent of the plasmid profile of the strains tested and had a DNA detection limit of 50 fg. In summary, the assay developed from the sequence identified in this study is a rapid, sensitive, and specific tool for the identification of B. anthracis. The assay can be used to complement assays for virulence markers by providing the genetic background of the organism containing plasmid pXO1 or pXO2, or both.
The assay described in this report has a provisional patent pending as of the date of publication.
Present address: 300B West Knight St., Dugway, UT 84022. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»