Previous Article | Next Article ![]()
Journal of Clinical Microbiology, August 2004, p. 3649-3654, Vol. 42, No. 8
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.8.3649-3654.2004
Anne M. Whitney, Ryan T. Novak, and Tanja Popovic
Meningitis and Special Pathogens Branch, Division of Bacterial and Mycotic Diseases, National Center for Infectious Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia 30333
Received 2 January 2004/ Returned for modification 22 February 2004/ Accepted 24 April 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Brucellosis impacts public health and agricultural economies worldwide because of its high infectivity rate (13). Brucella spp. have also long been considered a potential biological weapon (20), and currently, with the renewed threat of biological warfare and agricultural terrorism, B. melitensis, B. suis, and B. abortus are listed as category B biothreat agents by the Centers for Disease Control and Prevention Strategic Planning Group (34).
Infections in humans generally result from (i) transmission via the gastrointestinal route by the consumption of unpasteurized diary products (18, 32) and contaminated meat (25), (ii) airborne transmission in animal husbandry by inhaling dust contaminated by aborted tissues (23), and (iii) transmission caused by laboratory-associated exposure to aerosols (26). In the event of a bioterror attack, the preferred method of dissemination would most likely be via aerosol (3, 20, 34).
Because of the limited availability of animal vaccines, cost of animal inoculations, a lack of vaccines for human use (13, 37), and its low infectious dose for humans, human brucellosis is endemic in many parts of the world, including the Mediterranean region, Latin America, Asia, and Africa. The reported incidence varies from <0.01 to >200 per 100,000 population (13). Human brucellosis is rare in the United States, with approximately 100 human cases reported per year, mostly caused by consumption of unpasteurized diary products and, to a lesser degree, occupational exposure (10, 38).
The rapidity with which a laboratory diagnosis can be obtained is an important component of every outbreak investigation, because knowledge about the causative agent plays a pivotal role in implementing appropriate public health decisions in a timely manner (21, 34). At present, diagnosis of human brucellosis depends primarily on isolation of Brucella spp. from blood (43), followed by performing a set of bacteriological tests that allows for reliable identification to the subspecies (biotype) level but that requires up to 7 days for completion (1, 42). Another diagnostic option is the use of serologic assays; however, the specificity of these tests is inconsistent because of cross-reactivity with bacteria closely related to the genus Brucella (29).
Extensive efforts have been expended on the development of molecular diagnostic assays based on amplification of different genomic targets by the PCR for the identification of Brucella spp., as recently reviewed by Bricker (6). Although these molecular approaches appeared to be faster and more sensitive than traditional bacteriological tests, they are not considered to be confirmatory tests due to limited specificity (8, 9, 14). Recently a multiplex system has been developed that is sensitive for Brucella spp. and is able to differentiate between B. melitensis and B. abortus (33). However, discrepant results were observed with some B. abortus isolates. So far, none of these assays has been accepted for common use in diagnostic laboratories.
In collaboration with the Association of Public Health Laboratories, the Centers for Disease Control and Prevention Laboratory Response Network has developed standard algorithms based on traditional bacteriological tests (1) for the identification of Brucella spp. by Laboratory Response Network laboratories. Our current efforts are aimed at developing and validating more rapid methods that will allow confirmatory identification of Brucella spp. from cultures and, subsequently, directly in clinical specimen samples.
We chose to evaluate 16S rRNA gene sequencing to determine its utility for confirmatory identification of Brucella spp. The advantage of this method is that results can be obtained within 1 day as compared to 7 days by traditional microbiological testing. Previous work on other bacteria has indicated that differences in 16S rRNA gene sequences may be useful for subtyping or for the differentiation of virulent subtypes from nonvirulent subtypes (30, 36). Low variability in the 16S rRNA locus has been noted as an impediment in using 16S rRNA gene sequencing to discriminate at the species level (41); recent studies of other biothreat select agents have indicated that even subtle differences in the 16S rRNA gene sequence may be used for differentiating and identifying closely related species, which are often cross-reactive in biochemical identification systems commonly used in diagnostic laboratories (19, 35). In this study, we have primarily focused on the applicability of 16S rRNA gene sequencing as a rapid confirmatory identification tool and consequently have sequenced the near-full-length 16S rRNA genes of a reference set of Brucella isolates that represent the species commonly encountered in animal and human infections. In addition, we sequenced 17 closely related diagnostically challenging strains.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
16S rRNA gene sequencing. Sequencing primers were chosen from a panel of eubacterial primers: F8, R1492 (described above), F357 (5' TACGGGAGGCAGCAG 3'), R357 (5' CTGCTGCCTCCCGTA 3'), F530 (5' CAGCAGCCGCGGTAATAC 3'), R530 (5' GTATTACCGCGGCTGCTG 3'), F790 (5' ATTAGATACCCTGGTAG 3'), R790 (5' CTACCAGGGTATCTAAT 3'), F1068 (5' GTCGTCAGCTCGTGTCGTGAG 3'), F1083 (5' CGTGACATGTTGGGTTAAGTC 3'), F981 (5' CCCGCAACGAGCGCAACCC 3'), and R981 (5' GGGTTGCGCTCGTTGCGGG 3') (17, 28, 35). For non-Brucella strains sequenced for comparison, four additional eubacterial primers were used: R1333 (5' CTAGCGATTCCGACTTCATGC 3'), R180 (5' TCTCTCAAGACGTATGCGGTA 3'), R591 (5' CATCCTGCTTAAGTAACCGTC 3'), and F1127 (5' ATTAGTTGCCATCATTCAGTT 3'). Sequencing was performed with an Applied Biosystems (ABI, Foster City, Calif.) BigDye terminator cycle sequencing kit (version 2.0, 1.1, or 1.0) per the manufacturer's instructions, with the exception of using 6 µl of BigDye instead of 8 µl. Sequencing products were purified by using Centri-Sep spin columns (Princeton Separations, Adelphia, N.J.), and these sequencing products were resolved with an ABI model 3100 automated DNA sequencing system (ABI).
Computer analysis of 16S rRNA gene sequences. (i) Evaluation of the 16S rRNA gene sequences of Brucella spp. The 16S rRNA gene sequencing of the amplified products obtained with PCR primers F8 and R1492 generated 1,412-bp-long sequences that closely matched the size of previously published Brucella 16S rRNA gene sequences. The length of these sequences was close to the full length of the gene, which is 1,485 bp (5, 7, 16), and thus allowed direct comparisons.
The software package used for all data analysis was the Accelrys, Inc. (San Diego, Calif.), GCG Wisconsin package, version 10.2. The raw trace files from the ABI 3100 sequencer were visually examined, aligned, and edited (SEQMERGE).
We compared the Brucella sequences by using DISTANCES to determine the level of similarity among them. The sequences generated in this study were compared with previously published Brucella spp. sequences by using BESTFIT to determine levels of similarity between sequences.
(ii) Comparisons of 16S rRNA gene sequences of Brucella spp. to sequences of closely related strains. The eubacterial primers F8 and R1492 were also used for the panel of related strains likely to present a differential diagnostic challenge. The 16S rRNA gene sequences from the amplification products that resulted from F8 and R1492 varied in length (1,377 to 1,491 bp). For the Agrobacterium, Ochrobactrum anthropi, Oligella urethralis, and Vibrio cholerae strains, the sequences were organized into sets by genus (Table 2) and aligned (PILEUP). A core segment of 1.4 to 1.5 kb was selected for each that was common to all the sequences in a given set by trimming sequences at the 5' and 3' ends. The level of similarity of the sequences in a set was then determined by using Jukes-Cantor correction (DISTANCES). The sequences within the Agrobacterium set matched each other 100%, as did the sequences within the O. anthropi and V. cholerae sets (Table 2). The most divergent sequences within the O. urethralis set matched at 99.9% (Table 2). A representative sequence from each set, as well as the 16S rRNA gene sequences from Escherichia coli O157:H7, Bartonella henselae, Afipia broomeae, and Afipia felis, were then compared to the Brucella consensus sequence to determine the level of similarity by BESTFIT (Table 3).
|
| RESULTS |
|---|
|
|
|---|
A BLAST search (10 October 2003) on GenBank indicated that there were only 17 Brucella 16S rRNA gene sequences (from 10 strains) near 1.4 kb available in this public database (2). Eleven of these gene sequences were a 100% match to our Brucella consensus sequence. Six of these 11 sequences were from B. melitensis strain 16 M and B. suis strain 1330, for which full genome sequences were published in 2002 (15, 31). Both strains have two chromosomes, and each strain has three copies of the 16S rRNA gene with two copies on chromosome I and one copy on chromosome II (15, 31). A BLAST analysis (8 October 2003) of both genome sequences on The Institute for Genomic Research website (http://tigrblast.tigr.org/cmr-BLAST/) (2) revealed that all three copies from each strain were a 100% match to our Brucella 16S rRNA gene consensus sequence. Two additional 16S rRNA sequences that were a 100% match to our Brucella consensus sequence are those of the recently published marine Brucella strain M2357/93 and B. abortus strain 544 (GenBank accession no. AF091353 and AF091354) (7). Finally, the last three 16S rRNA gene sequences of the 11 are from a study published in 2000 that also used the B. melitensis strain 16 M (accession no. AF220149, AF220148, and AF220147) (5).
There were six near-full-length Brucella 16S rRNA gene sequences that were not 100% matches to the consensus sequence. These sequences were entered into GenBank between 1989 and 1995. The first one is that of B. ovis strain ATCC 25840 (accession no. L26168) (K. H. Wilson, unpublished data), which differs from the consensus sequence by 3 bp. Compared to the consensus sequence, this B. ovis sequence has an A instead of a G at position 600, as well as Ns at positions 426 and 700. We sequenced three B. ovis isolates, including the ATCC 25840 strain, and all three sequences were a 100% match to the Brucella consensus sequence. The Brucella consensus sequence also differs substantially from the 16S rRNA gene sequence generated in 1989 for B. abortus strain 11/19 (accession no. X13695) (16). A BESTFIT analysis indicates that there is a 19-bp difference, including an N as well as a 3-bp deletion compared to our Brucella consensus sequence. A BLAST query on GenBank indicates that the sequence of accession no. X13695 is unique and does not perfectly match any other 16S rRNA gene sequence in the public database. The remaining four Brucella 16S rRNA gene sequences in GenBank are from B. suis strain ATCC 23444, B. melitensis strain ATCC 23456, B. neotomae strain ATCC 23459 (accession no. L26169, L26166, and L26167, respectively) (Wilson, unpublished), and B. canis strain ATCC 23365 (accession no. L37584) (K. H. Wilson and H. G. Hills, unpublished data), which differed from the consensus sequence by 1 to 6 bp. Three of these four sequences had Ns in the sequences, indicating an inability to identify bases.
B. abortus strain 11/19 was not available for resequencing: however, this strain is the same as U.S. strain 19 (E. Moreno, personal communication). We sequenced the 16S rRNA genes from B. abortus U.S. strain 19, B. melitensis strain ATCC 23456, B. neotomae strain ATCC 23459, and B. canis strain ATCC 23365. All four sequences were a 100% match to the Brucella consensus sequence.
Comparisons with the panel of 17 related strains likely to present a differential diagnostic challenge indicated that the sequence most similar to the Brucella consensus sequence was that from an O. anthropi strain at 98.8% similarity (Table 3), which is consistent with the relationship seen in previous phylogenetic studies (39). There is a difference of 16 bp, as well as a deletion and an insertion, which easily differentiates the Brucella consensus sequence from this sequence. The other 16S rRNA gene sequences from this panel of strains had even more pronounced differences (Table 3).
| DISCUSSION |
|---|
|
|
|---|
The critical advantage of 16S rRNA gene sequencing is based on its ability to be used for identification of unknown isolates. By using the broad-range eubacterial primers, as demonstrated in this study, a 16S rRNA gene sequence can be readily obtained, allowing for rapid confirmatory identification of an unknown isolate as a Brucella species by simple comparison to the consensus sequence. Confirmatory identification can then serve as a guide to continue all subsequent work with biosafety level 3 practices, while the rapidity with which it is obtained may be the critical factor needed for quickly implementing control and preventative measures in an emergency biothreat situation (6, 34). The PCR tests previously proposed limit the laboratorian to either positively identifying gene regions associated with a Brucella species or ruling out Brucella spp., which would then require subsequent diagnostic tests to identify the etiologic agent. Sequencing the 16S rRNA gene enables querying a database such as GenBank to yield an alternative identification without the need for any additional laboratory work. Consequently, given that the clinical presentation of brucellosis can resemble diseases caused by other select agents, this global diagnostic approach will be invaluable in a potential biothreat setting. We have added to the robustness of this approach by submitting to GenBank 82 high-quality 16S rRNA gene sequences from Brucella spp. and close relatives, which allows for a more accurate comparison of sequences.
One additional advantage of using 16S rRNA gene sequencing is that this method may be able to identify the causative organism from clinical specimens of patients who have received antimicrobial agents prior to the collection of specimens, as has been demonstrated recently during the 2001 bioterrorism-associated anthrax investigation (35). Because antimicrobial and other supportive treatment will not differ regardless of the Brucella species causing the illness (4), detection of the Brucella-specific 16S rRNA gene sequence will be advantageous for timely treatment of patients even without species identification.
Our result of 100% identity of the 16S rRNA gene sequence among the large number of Brucella strains tested is consistent with previous studies using a limited number of strains (5, 7) and could also be used to support proposals that Brucella may be a monospecific genus, based on the DNA hybridization studies that indicate a >77% level of homogeneity for the Brucella spp. (40). Vizcaíno et al. observed that there was some variability in the 16S rRNA locus, but not enough to separate the Brucella species (41). Similarly, there are some Brucella 16S rRNA gene sequences in GenBank that differ from the consensus sequence. We resequenced all but one of these discrepant GenBank submissions. Each was a 100% match to the Brucella consensus sequence. We do not believe the differences in the sequences from the original GenBank entries represent microheterogeneity in the 16S rRNA genes of these Brucella strains. Rather, these discrepant sequences were submitted from 1989 to 1995 before the advent of automated sequencing technologies with higher fidelity. This is supported by the presence of unresolved bases in some of these discrepant sequences.
A comparison of the Brucella 16S rRNA consensus sequence with those of close relatives indicates that sequencing of the first 500 bp is sufficient to differentiate them. The closest match, O. anthropi, was different by 2 bp in this segment, which is a sufficient level of divergence to be detected by current automated sequencing technology. For laboratories that may wish to conserve resources, this may be a viable option compared to sequencing 1.4 kb.
This work demonstrates that 16S rRNA gene sequencing can provide rapid confirmatory identification of an isolate as a Brucella species, allowing prompt and appropriate public health responses to be implemented. This approach will be of even greater value if shown to be effective when used directly on clinical specimens; studies are under way in our laboratory to evaluate its sensitivity and specificity on clinical specimens of patients with culture-confirmed brucellosis.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Present address: Saskatchewan Health, Provincial Laboratory, 3211 Albert St., Regina, Saskatchewan, Canada S4S 5W6. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|