|
|
||||||||
Diagnostic Systems Division1 and Bacteriology Division2 , United States Army Medical Research Institute of Infectious Diseases, Fort Detrick, Frederick, MD 21702-5011, USA
Correspondence
Ricky L. Ulrich
Ricky.Ulrich{at}AMEDD.ARMY.MIL
Received 30 September 2005
Accepted 13 January 2006
Abbreviations: CT value, crossing threshold value; EPF, end-point fluorescence; MGB probe, minor groove binding probe; p.e., post exposure.
| INTRODUCTION |
|---|
|
|
|---|
Acute human glanders is normally fatal without immediate antibiotic treatment. Infected individuals may present clinically with swollen lymph nodes, ulcerating nodules in the alimentary and respiratory tracts, and weight loss, as well as numerous subcutaneous abscesses (DeShazer & Waag, 2004; Srinivasan et al., 2001). The route of exposure for human glanders is normally via contact with infectious material (i.e. nasal discharge from infected horses) through skin abrasions or cuts, which usually results in a chronic form of disease. In contrast, and likely the route of exposure in the deliberate release of B. mallei, contraction from aerosols causes acute glanders, and if not treated immediately, death occurs in a matter of days (Howe & Miller, 1947; Hunter, 1936).
With the potential of B. mallei weaponization, rapid methods for the identification of this infectious Burkholderia species are essential to ensure that adequate prophylactic therapies can be initiated following exposure. The conventional laboratory methods for detecting B. mallei rely on biochemical testing, colony morphology, and motility assays, all of which can take up to 7 days to obtain preliminary results. In addition, it has recently been demonstrated that the API 20NE and RapID NF commercial kits are 060 % successful in positively identifying B. mallei (Glass & Popovic, 2005).
With the physiological and genetic relatedness between B. mallei, Burkholderia pseudomallei and Burkholderia thailandensis (Brett et al., 1997; Godoy et al., 2003; Holden et al., 2004; Nierman et al., 2004), PCR-based and biochemical techniques for identifying and differentiating these Burkholderia species are difficult. For B. pseudomallei, the 16S and 23S rRNA, 16S23S intergenic region, fliC, heat-shock protein 70 (Hsp70), and a type three secretion (TTS) system have been targeted for the development of molecular-based assays (Antonov et al., 2005; Bauernfeind et al., 1998; Gee et al., 2003; Hagen et al., 2002; Lee et al., 2005; Sprague et al., 2002; Tanpiboonsak et al., 2004; Thibault et al., 2004; Tomaso et al., 2004, 2005; Tyler et al., 1995; U'Ren et al., 2005). One study required the use of two separate assays and a negative PCR result to discriminate B. mallei from B. pseudomallei. Several other investigations could differentiate B. pseudomallei from B. mallei, but these reports were based on a single base pair difference between the two species. The lack of a target region with substantial nucleotide variation between B. pseudomallei and B. mallei further demonstrates the genetic similarity between these Burkholderia species.
In this investigation, we have designed a real-time PCR assay for the definitive identification of B. mallei and separation from B. pseudomallei. The assay was tested for specificity against a panel of diverse bacterial species, and for the ability to detect B. mallei within the lungs, liver and spleen of aerosol-challenged BALB/c mice.
| METHODS |
|---|
|
|
|---|
|
Bacterial and tissue DNA extraction. Burkholderia species were cultured as described above, and genomic DNA from cultures as well as collected tissues was prepared using the Qiagen QIAamp DNA Mini kit with 180 µl of an overnight culture or homogenized tissue, with the following modifications to the manufacturer's protocol. The sample (180 µl) was mixed with 20 µl proteinase K (17·8 mg ml1) and 200 µl buffer AL, and incubated at 55 °C for 1 h. Ethanol (96100 %, 210 µl) was added and mixed, and the sample was loaded onto the column and the column washed using buffers AW1 and AW2. Columns were centrifuged at 16 000 g for 3 min to dry residual ethanol from the column. Buffer AE (100 µl), preheated to 70 °C, was added to each column, incubated at 70 °C for 5 min, and DNA was eluted.
The DNA used in this study for cross-reactivity testing was from reference material obtained from recognized culture collections, commercial vendors, clinics, or unique entries from previous United States Army Medical Research Institute of Infectious Diseases (USAMRIID) collections. Pure cultures of each agent were grown on appropriate media at optimal temperatures and harvested for DNA isolation. The DNA was extracted using either Bactozol Kits (Molecular Research Center) according to the manufacturer's recommendations or Qiagen kits as previously described (Coyne et al., 2004). Extracted DNA was quantified and checked for purity with a spectrophotometer by measuring absorbance at 260, 280 and 320 nm. The ability of the DNA to be amplified was established with a universal primer set to the 16S rRNA gene (Bricker & Halling, 1994). The presence of high-molecular-mass DNA was confirmed by running 400 ng ml1 of each genomic DNA on 0·7 % agarose gels. Identity was confirmed for each DNA by sequence analysis of the entire 16S rRNA gene and comparison with published data.
Assay design and optimization.
The real-time PCR assay primer and probe sequences are listed in Table 1
. The B. mallei bimAma (Burkholderia intracellular motility A) homologue encoded by gene locus BMAA0749 was used for assay development. It has recently been demonstrated that the N-terminal nucleotide sequence of bimAma contains a unique B. mallei region not present in the B. pseudomallei and B. thailandensis bimA genes (Stevens et al., 2005a). The specific primer and TaqMan minor groove binding (MGB) probe sequences were designed using Visual OMP version 4.1.0.0 for Windows (DNA Software). All primer and probe sequences were analysed using BLASTN to ensure specificity for their respective targets. All primers were purchased from Invitrogen. The probes were synthesized by Applied Biosystems, and contained 6-carboxyfluorescein (FAM) at the 5' end and a non-fluorescent quencher with the MGB protein at the 3' end. Primer candidates were used to amplify 10 pg target template in the presence of SYBR green dye. Melt curves and agarose gels were analysed to eliminate inefficient and/or dimer-producing primer pairs. MgCl2 and primer concentrations were optimized sequentially as follows: MgCl2 was optimized in 1 mM increments from 3 to 7 mM, and primers were optimized symmetrically in 0·1 µM increments from 0·1 to 1·0 µM. The combinations that exhibited the earliest crossing threshold (CT) value and generated the highest end-point fluorescence (EPF) value were chosen as the optimal conditions for each assay (Table 1
).
|
Assay linearity and estimated detection limit. The linearity of the assay was assessed by standard curve analysis in order to determine the amplification efficacy, efficiency, and quantitative utility. Additionally, serial dilutions were used to estimate the limit of detection. Serial dilutions in Tris/EDTA buffer (Eppendorf) were prepared at the following per-reaction concentrations: 1 ng, 500 pg, 100 pg, 50 pg, 10 pg, 5 pg, 1 pg, 500 fg, 100 fg. Triplicate assays were performed at each dilution, the standard curve was generated by the LightCycler 4.0 software, and slopes were used to calculate amplification efficiency and standard error. An estimate of the detection limit was determined by the lowest concentration of template per reaction that produced three positive results.
Inclusivity and exclusivity experiments.
A general panel of 60 organisms was analysed using the optimized conditions to establish inclusivity and exclusivity (Table 2
). In addition, a panel of 30 genetic near neighbours was also analysed (Table 3
). Panels were constructed to include: other threat organisms; nearest genetic neighbours to threat organisms; organisms sharing a clinical niche with B. mallei, especially respiratory pathogens, opportunists, and normal respiratory flora; organisms observed repeatedly in clinical and environmental samples. In all cases, 100 pg genomic DNA per reaction was used to determine if the assays cross-reacted with nucleic acid from other organisms. Inclusivity and exclusivity testing was performed on the LightCycler 2.0 instrument and utilized the automated calling capability of LightCycler software version 4.0.
|
| RESULTS |
|---|
|
|
|---|
Assay optimization
Primers were designed from a portion of the 5' region of the B. mallei bimAma sequence and initially tested by rapid PCR amplification in the presence of SYBR Green. All primer combinations resulting in PCR products smaller than 160 bp were tested for amplification efficiency. Initial testing was performed with five forward and five reverse primers in a matrix that produced 13 combinations. All combinations amplified B. mallei DNA from strain ATCC 23344. Two combinations, F49 with R144 and F52 with R146, generated the largest amount of amplified DNA. In addition, no template controls with these two primer sets produced minimal to no primer dimers. After primers were tested and selected using standard PCR, five compatible probes were tested for performance with both of these primer sets. For both primer sets, probe P104R consistently produced the earliest CT and the highest EPF values. Assays were further optimized for primer and MgCl2 concentrations to produce maximum sensitivity, earliest CT and highest EPF. The final primer and probe pairs and optimized assay conditions are listed in Table 1
.
Sensitivity/efficiency
Standard curve analysis was utilized to assess the amplification efficiency and approximate detection limit of the assays. Standard curves were designed to contain triplicate samples at a minimum of six detectable concentrations of DNA. The assays had amplification efficiencies of 1·97 (F49/R144) and 1·92 (F52/R146) with errors of 0·018 and 0·057, respectively. An amplification efficiency of 2·00 is considered ideal, and corresponds to a doubling of copy number for every PCR cycle. Both assays detected all three replicates at 1 pg DNA template per reaction, which corresponds to 424 genome copies, based on a genome size of 2·3 Mbp (this corresponds to B. mallei ATCC 23344 chromosome 2).
Inclusivity/exclusivity
Using standard PCR (Epicentre FailSafe kit with buffer J, Epicentre Technologies), with primers targeting the unique 5' DNA sequence of the B. mallei bimAma (AT5, 5'-TTCGATCGATTCCTGCTATC-3', and AT4, 5'-GCGTTAAACGCCGTACTTTC-3'), we have previously tested a collection of 29 B. mallei and 34 B. pseudomallei isolates to determine if this unique B. mallei nucleotide sequence is conserved (Ulrich et al., 2006). Surprisingly, given the extensive genetic similarity between B. mallei and B. pseudomallei (Godoy et al., 2003), only B. mallei genomic DNA produced a detectable amplicon with PCR primers (listed above) specifically designed for bmaAma. Likewise, both real-time PCR assays (Table 1
) were analysed for their ability to specifically detect strains of B. mallei in the absence of non-specific detection of both genetic near neighbours as well as a broader collection of micro-organisms. The near-neighbour analysis included a panel of 30 DNAs, including eight strains of B. thailandensis, eight strains of B. mallei, eight strains of B. pseudomallei and six strains of other Burkholderia species, including B. cepacia, B. multivorans, B. stabilis and B. vietnamiensis. As observed with primer pairs AT4 and AT5 (standard PCR primers targeting bimAma), both assays detected only B. mallei, and consistent CT values were observed across all strains and assays tested (Table 3
). Both assays were negative for other genetic near neighbours, as well as for a panel of 60 micro-organisms from a general microbial DNA panel (Tables 2 and 3![]()
).
In vivo analysis
To determine if our real-time PCR assays were capable of detecting B. mallei within infected tissues and blood, we challenged (10 LD50 or 1000 c.f.u.) 20 female BALB/c mice via the aerosol route and followed organ colonization for 48 h post exposure (p.e.). At 1 h p.e., our real-time PCR primer and probe sets definitively identified B. mallei within the lungs of mice, indicated by mean CT values of 30·06 for primer pair 49F/144R and 30·15 for set 52F/146R (Table 4
). At 24 h p.e., dissemination of B. mallei to the liver was detected in two out of three animals with primer pairs 49F/144R (mean CT of 34·56) and 52F/146R (mean CT of 35·19) (Table 4
). Within the spleen, primers 52F/146R detected B. mallei in two out of three animals (mean cycle value of 35·57), whereas primers 49F/144R only detected B. mallei in one out of the three animals (CT of 33·42) at 24 h p.e. (Table 4
). At 48 h post challenge, all organ extracts (lungs, livers and spleens) tested positive for B. mallei (Table 4
). The mean CT values for the triplicate lungs were 23·50 for primer pair 49F/144R and 23·59 for primer pair 52F/146R. For extracted livers, the mean CT values were 28·50 and 28·82 for 49F/144R and 52F/146R, respectively, and for spleen tissue, 31·56 with pair 49F/144R and 31·60 with set 52F/146R. It should be noted that BALB/c mice challenged by the aerosol route with B. mallei develop acute glanders and normally succumb to infection 4872 h post challenge; therefore, additional time-points were not analysed. All control organ extracts, from non-infected animals housed within the same biosafety level 3 containment laboratory, were negative for the presence of B. mallei by analysis with the real-time PCR primers and probes described in this study (Table 4
). Surprisingly, and likely a result of the low bacterial burden in the blood, our real-time PCR assays were not capable of detecting B. mallei in the blood of infected mice (Table 4
). The presence of PCR inhibitors in all sample eluates was tested by using an inhibition assay previously described (Hartman et al., 2005). The extracted DNA material from whole blood required a dilution of 1 : 8 to overcome the effects of PCR inhibitors present in the sample (Table 4
).
|
| DISCUSSION |
|---|
|
|
|---|
As previously mentioned, the current laboratory methods for identifying B. mallei in clinical samples rely on biochemical and phenotypical testing. However, it has been shown that commercially available biochemical testing kits (API 20NE and RapID NF) cross-react with other non-virulent bacterial species, leading to false-positive reactions (Glass & Popovic, 2005; Inglis et al., 1998). An ELISA has been developed using irradiation-killed B. mallei whole cells; however, cross-reactivity with irradiation-killed B. pseudomallei was observed, suggesting that this immunological assay is not B. mallei specific (unpublished results).
The development of PCR-based methodologies for the identification and differentiation of B. mallei from B. pseudomallei and B. thailandensis has been troublesome due to the level of genetic similarity between these Burkholderia species (Brett et al., 1998; Godoy et al., 2003; Holden et al., 2004; Nierman et al., 2004). Various chromosomal markers have been used to detect B. pseudomallei, and include the 16S and 23S rRNA, 16S23S intergenic region, the filament-forming flagellin fliC gene, heat-shock protein Hsp70, and a putative type three secretion operon (Antonov et al., 2005; Bauernfeind et al., 1998; Gee et al., 2003; Hagen et al., 2002; Lee et al., 2005; Sprague et al., 2002; Thibault et al., 2004; Tomaso et al., 2004, 2005; Tyler et al., 1995; U'Ren et al., 2005). Several of the assays that have been described are targeted toward single-nucleotide differences in the 23S rDNA region, the 16S region, and a putative antibiotic resistance gene. However, targeting a single base pair could be problematic for the identification of B. mallei, as it has been shown that the B. mallei genome is unstable upon passage. B. mallei ATCC 23344 has been passaged in culture, in a mouse and in a horse, and two isolates have been obtained from the accidental infection of a scientist (Srinivasan et al., 2001). The genomes of these isolates were sequenced and compared to the published B. mallei sequence. Single-nucleotide polymorphisms were noted after passage, mainly located within coding regions of the genome. The two human isolates obtained from the accidental infection each possessed a distinct set of sequence variations as well as two sequence variations in common (unpublished data). These data suggest that the B. mallei population in a host is not identical and that accumulated genome-sequence variation may occur.
The B. mallei-specific real-time PCR assay described in this study was designed using a unique 5' nucleotide region encoded within bimAma (BMAA0749) (Stevens et al., 2005a). In B. pseudomallei, BimA has been shown to be involved in host-cell actin polymerization (Stevens et al., 2005b). Surprisingly, the unique DNA sequence encoding the N-terminal region of BimAma, the portion that is exposed at the bacterial cell surface, is conserved among numerous B. mallei isolates (29 strains obtained from various clinical sources over a 70 year period; Ulrich et al., 2006), suggesting that this sequence divergence occurred early in the evolution of B. mallei from B. pseudomallei. These findings also suggest that this chromosomally encoded gene is genetically stable and will provide an ideal target for the rapid identification of B. mallei. Despite these deviations in nucleotide sequence, heterologous expression of the B. mallei bimAma in a B. pseudomallei bimAps mutant fully complements the defective actin motility phenotype observed in this strain, suggesting that bimAma encodes a functional protein (Stevens et al., 2005a).
The specificity of our B. mallei real-time PCR assays was tested against an extensive cross-reactivity panel comprising bacterial (Gram-positive and -negative species), viral, fungal and human DNAs (Table 2
). In addition to the cross-reactivity panel, a set of 30 closely related Burkholderia species was tested, including B. mallei, B. pseudomallei, B. thailandensis, B. stabilis, B. multivorans, B. cepacia and B. vietnamiensis (Table 3
). With our real-time PCR primer and probe sets, B. mallei was definitively identified and differentiated from B. pseudomallei and other closely related species (Table 3
).
To further test the applicability of our B. mallei-specific real-time PCR assay, serial animal sacrifices were performed on BALB/c mice infected by the aerosol route with B. mallei ATCC 23344. As anticipated, B. mallei was positively detected in organ extracts processed at 1 h (lungs), 24 h (lungs, spleen and liver) and 48 h p.e. (lungs, spleen and liver) (Table 4
). The data produced in this study indicate the proliferation of bacterial cells in the organs of infected mice based on the decreasing CT values for each of the time-points tested. These findings suggest that the real-time PCR primer and probe sets described in this work will provide a rapid method for the identification of B. mallei in clinical samples. Unlike previously described techniques for detecting B. mallei in human clinical samples (i.e. culturing abscess biopsies, GLC of cellular fatty acids, and 16S rRNA gene-sequence analysis), real-time PCR provides an alternative method for identifying B. mallei and eliminates the need for bacterial pre-enrichment (i.e. total DNA can be isolated from organ abscesses and real-time PCR performed). Unfortunately, the sensitivity (see below) of our B. mallei-specific real-time PCR assays was not sufficient to detect B. mallei in the blood of infected animals. These findings are likely the result of either the low bacterial load found within the blood of acutely infected BALB/c mice or the presence of PCR inhibitors (Fritz et al., 2000). Even when incorporating spiked BALB/c control blood with B. mallei at concentrations of 1000, 500 and 100 c.f.u. ml1, we failed to identify this pathogen using the primer and probe sets described in this work (data not shown). However, we were able to detect B. mallei in spiked human blood (human use protocol FY04-14) at a concentration of 500 c.f.u. ml1 with primer set 49F/144R (three out of three replicates) (data not shown). These findings suggest that, despite using a commonly employed commercial DNA extraction kit (Qiagen QIAamp DNA Mini kit), mouse blood contains PCR inhibitors not present in human blood. Further studies to address these extraction limitations will be needed. Similar reductions in the ability to detect B. mallei and B. pseudomallei in spiked blood have been reported by Tomaso et al. (2005), who targeted the B. pseudomallei and B. mallei 16S rDNA, fliC and rpsU genes.
For serological detection of B. mallei, the Russian scientists Gelman and Kalning developed the mallein assay in the late 1800s (DeShazer & Waag, 2004). This immunological test utilized antigens prepared from 48 month culture filtrates of B. mallei and was utilized by the USA and Canada to identify horses infected with B. mallei (Steele, 1979). Despite the rapid (48 h) identification of B. mallei using the mallein test in horses, serological conversion does not occur in humans for 34 weeks p.e., which suggests that this assay is ineffective for detecting human glanders. It should be noted that this serological test also cross-reacts with B. pseudomallei (DeShazer & Waag, 2004). To date, the only B. mallei-specific detection assay reported in the literature incorporates the bacterial phage phiE125, which has been shown to be specific for B. mallei capable of biosynthesizing LPS (Woods et al., 2002).
The findings of this investigation provide two real-time PCR assays for the specific identification of B. mallei and its differentiation from B. pseudomallei and closely related Burkholderia species. The described assays did not cross-react with any of the bacterial species included in our culture collection, and were capable of detecting B. mallei in the lungs, spleen and liver of infected mice beginning at 1 h p.e. This B. mallei-specific real-time PCR assay will undoubtedly complement other methodologies used in the clinical laboratory for the rapid identification of this obligate mammalian pathogen.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | J MED MICROBIOL | MICROBIOLOGY | J GEN VIROL | ALL SGM JOURNALS |