Previous Article | Next Article ![]()
Journal of Clinical Microbiology, December 2008, p. 3935-3940, Vol. 46, No. 12
0095-1137/08/$08.00+0 doi:10.1128/JCM.00464-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Department of Pathology and Laboratory Medicine,1 Department of Microbiology and Immunology, American University of Beirut, Beirut, Lebanon,3 World Health Organization Regional Collaborating Center on Human Brucellosis, Beirut, Lebanon,2 Université Paris-Sud 11, CNRS, UMR8621, Institut de Génétique et Microbiologie, Orsay 91405, France,4 Department of Analytical Microbiology, Centre d'Etudes du Bouchet, Vert le Petit 91710, France,5 DGA/D4S-Mission pour la Recherche et l'Innovation Scientifique, Bagneux 92220, France6
Received 9 March 2008/ Returned for modification 6 July 2008/ Accepted 30 September 2008
|
|
|---|
|
|
|---|
Species identification and subtyping of Brucella culture isolates is important for epidemiologic surveillance; investigation of outbreaks in regions of both endemicity and nonendemicity; and distinguishing cases of human reinfection from relapse, thereby influencing clinical therapeutic decisions (3).
Terrestrial Brucella spp. are homogeneous and harbor >80% interspecies homology by DNA-DNA hybridization studies (31, 38), identical 16S rRNA sequences (15), and >98% sequence similarity by comparative genomics (18, 29). Molecular methods commonly used for the subtyping of isolates of other bacterial species, such as multilocus enzyme electrophoresis (12), pulsed-field gel electrophoresis (4), random amplified polymorphic DNA analysis (34), enterobacterial repetitive intergenic consensus sequence PCR (26, 35), repetitive intergenic palindromic sequence PCR (25), amplified fragment length polymorphism analysis (42), and monolocus (such as omp2a and omp2b) sequence analysis (9) or multilocus sequence typing (40), are able to segregate Brucella isolates into the recognized species and certain species biovars at best. However, these methods are not sufficiently discriminatory for the routine subtyping of isolates for epidemiological trace-back purposes.
Tandemly repeated sequences used in forensic investigations have also been found in bacteria (for a review, see reference 39). Minisatellites have repeat unit sizes of 9 bp or greater, and microsatellites have repeat unit sizes of up to 8 bp (39). Combinations of minisatellite and microsatellite repeats in multilocus variable-number tandem-repeat analysis (MLVA) have proven highly discriminatory in the epidemiological subtyping of isolates belonging to monomorphic bacterial species, such as Bacillus anthracis, Yersinia pestis (8, 23), and more recently, Brucella spp. In Brucella, MLVA schemes with 21 loci (MLVA-21) and MLVA-16 (1, 22, 41) that use a combination of repeat markers distributed across the Brucella genome were able to distinguish isolates of Brucella spp. of widespread temporal and geographical origins or of very close origins (13). Furthermore, with some rare exceptions, which might be due to a lack of resolution and incorrect clustering produced by the biotyping methods, the isolates formed MLVA groups that corresponded to the known Brucella species. A study of 128 human isolates of B. melitensis by MLVA-16 identified the Americas, West Mediterranean, and East Mediterranean groups and detected remarkable degrees of diversity within each group (1). The relative stability of MLVA loci has been estimated directly by repetitive subculturing of a few strains in the laboratory (41) and by the typing of B. melitensis isolates of the Rev.1 vaccine strain obtained from different sources (13), so that the most recent MLVA data analyses have started to assign different weights to different markers, taking into account the mutation rate and the associated levels of homoplasy (1). The MLVA schemes that have been devised are very promising, but the strength of a typing method depends as much on the technique itself as on the associated genotyping database. These typing assays still need to be applied to larger numbers of isolates from as many countries as possible to validate their reliability in investigating strain relatedness both in countries of endemicity, especially where efforts to eradicate human disease may be attempted, and in brucellosis-free countries where a cluster of infections may be suspected to be due to a bioterrorism event. The primary purpose of high-resolution typing in a format compatible with the production of large-scale international databases is the long-term surveillance of infectious diseases and the early detection of abnormal events, including temporal or geographical shifts in the population structure, e.g., the emergence and spread of new clones or lineages. The current Brucella2007 MLVA database, hosted at http://mlva.u-psud.fr, contains data derived from more than 500 animal and human Brucella isolates. It is expected that in the future, similar and compatible databases with data for isolates from wider geographical origins will be developed and made jointly and freely accessible to querying via the use of web-based analysis tools (16).
In this study, we applied the MLVA-16 scheme (1, 22) to a series of clinical isolates of Brucella accrued over a 5-year period at a single tertiary-care center in Lebanon, a small country where human brucellosis is endemic. The genotypes of these isolates were compared to the genotypes previously determined by MLVA-16 typing and are included in the Brucella2007 MLVA database.
(This study was presented in part at the 107th General Meeting of the American Society for Microbiology, Toronto, Ontario, Canada, 2007.)
|
|
|---|
Brucella MLVA-16 genotyping scheme. MLVA was performed with all 42 isolates according to the scheme initially proposed by Le Flèche et al. (22), which includes 15 tandem-repeat loci (MLVA-15), and modified by Al-Dahouk et al. to include 1 additional locus, bruce19 (MLVA-16) (1). PCR amplification of eight minisatellite loci in panel 1 (the bruce06, bruce08, bruce11, bruce12, bruce42, bruce43, bruce45, and bruce55 loci), three microsatellite loci in panel 2A (the bruce18, bruce19, and bruce21 loci), and five microsatellite loci in panel 2B (the bruce04, bruce07, bruce09, bruce16, and bruce30 loci) was carried out by a protocol similar to the one described previously (1).
DNA preparation. DNA was extracted from one loopful of bacterial cells grown for 48 h on chocolate agar, and single colonies were isolated by using the tissue protocol of the QIAamp DNA minikit (Qiagen, Hilden, Germany). DNA concentrations were measured by UV spectrophotometry (Shimadzu, Japan).
PCR amplification. PCR amplification was performed in a total volume of 15 µl containing 5 ng of DNA, 1x PCR buffer, 1 U of Taq DNA polymerase, 200 µM of each deoxynucleoside triphosphate, 0.3 mM each flanking primer, and 1 M of betaine (catalog no. B2629; Sigma-Aldrich). The PCR cycling parameters were as follows: initial denaturation at 96°C for 5 min, followed by 30 cycles of denaturation at 96°C for 30 s, primer annealing at 60°C for 30 s, and extension at 70°C for 1 min, with a final extension step at 70°C for 5 min. For detection of the PCR products, 5 µl of each of the PCR amplification products was subjected to gel electrophoresis with a 2% 3:1 HRB agarose gel (catalog no. E776; Amresco) for panel 1 and a 3% HRB agarose gel for panel 2 in 0.5x Tris-borate-EDTA buffer. Electrophoresis was carried until the bromophenol blue had run at least 20 cm. The size markers used were a 100-bp ladder (EZ Load 100-bp PCR molecular ruler; Bio-Rad) for panel 1 and a 20-bp ladder (EZ Load 20-bp molecular ruler; Bio-Rad) for panel 2. Gels that had been stained with ethidium bromide (0.5 µg/ml) were digitally acquired with a gel documentation system (DigiDoc-It; UVP, Upland, CA).
Analysis of MLVA data. Minisatellite repeat numbers (panel 1 [loci bruce06, bruce08, bruce11, bruce12, bruce42, bruce43, bruce45, and bruce55]) and microsatellite repeat numbers (panel 2A [loci bruce18, bruce19, and bruce21] and panel 2B [loci bruce04, bruce07, bruce09, bruce16, and bruce30) were calculated on the basis of the band sizes after the gel images were normalized by use of the Bionumerics (version 5.1) software package (Applied Maths, Belgium) and the previously published allele numbering convention (22). The reproducibilities of the MLVA profiles were assessed by comparing the typing results for experiments performed independently at two laboratory testing sites in Lebanon and France. The similarity matrix for each panel was calculated with Bionumerics software by using the repeat number at each locus as a character type with the categorical similarity coefficient. This categorical approach considers all alleles at each locus and within each panel to be equally distant (37). For MLVA-16 clustering analysis, an aggregated similarity matrix was produced by using a composite data set and giving weights of 20, 5, and 1 to panels 1, 2A, and 2B, respectively, according to the principles reported previously (1). Weights were assigned to the individual loci in each panel because all loci do not evolve at the same speed but have different mutation rates, as determined empirically. Clustering was done by use of the unweighted paired group method for arithmetic averages algorithm. Genetic diversity was determined by calculating the Hunter-Gaston (or Simpson's) diversity index (HGDI) of each MLVA marker and for the combination of markers for each of panel 1, panel 2A, panel 2B, panel 2 (panel 2A plus panel 2B), and panels 1 and 2 (20). The 95% confidence intervals of the HGDIs for indices of >0.85 were calculated by the method of Grundmann et al. (17). The MLVA-16 genotypes, defined by the combination of alleles, of the 42 isolates were compared to the corresponding data obtained for the reference strains and field isolates investigated previously (1, 13, 14, 22, 31, 32, 33). The present and previous data have been incorporated into a single database designated Brucella2007, which can be queried at http://mlva.u-psud.fr by the use of web-based analysis tools.
|
|
|---|
The typeability of the various MLVA markers in this series approached 1 (37). Only the bruce04 locus in strain S29 and the bruce09 locus in strain S19 could not be amplified, despite multiple attempts. All other amplification reactions were successful. The reproducibility of MLVA-16 was excellent, as experiments performed in two different laboratories in Lebanon and France by use of the same protocol and with the results interpreted separately in a blinded fashion yielded 100% concordant results. This is not surprising for a sequence-based approach and was recently confirmed by a ring trial for the panel 1 loci among 15 Brucella reference laboratories (http://mlva.u-psud.fr/BRUCELLA/). In that trial, the few discrepancies in typing results observed were due to errors in allele calling (24a). The 42 isolates in the series used in this study formed a homogeneous group for which the differences did not exceed 15% (Fig. 1). Limited diversity was observed with the panel 1 and panel 2A loci (Table 1), as would be expected for strains originating from a very restricted geographical area and by the use of markers that are the most phylogenetically informative. MLVA-16 yielded a total of 29 genotypes with seven clusters and 22 singleton genotypes. Previous work identified a total of 68, 52, and more than 300 genotypes for panel 1, panel 2A, and panel 2B, respectively (1, 13, 14, 22, 31, 32, 33) (see Fig. S1 in the supplemental material). Twenty-six panel 1 genotypes, numbered 41 to 66 (1, 22), corresponded to B. melitensis. No new panel 1 genotypes were identified in the present study. Thirty-seven isolates shared panel 1 genotype 43. Panel 1 genotypes 42, 44, and 57 were observed in a single isolate each; and genotype 60 was observed in isolates S18 and S24, which were derived from the same patient. Therefore, only 5 different panel 1 genotypes were observed among the 26 panel 1 B. melitensis genotypes observed so far. Accordingly, the bruce06, bruce11, bruce12, and bruce55 loci from panel 1 and the bruce21 locus from panel 2 showed no diversity among the 42 isolates (HGDI = 0). All five panel 1 genotypes represented in the present study correspond to the East Mediterranean group proposed previously (1) (see Fig. S1 in the supplemental material). The panel 2B markers displayed the most diversity in the present collection. This was as expected, given the high mutation rate of some panel 2B markers (13), with bruce04 and bruce16 exhibiting the highest HGDIs, followed by bruce30 (HGDIs > 0.7). Bruce07 and bruce09 in panel 2B had HGDIs of <0.5. Panel 2B and panels 2A and 2B yielded 27 and 28 genotypes, respectively, and had similar HGDIs. The combined panel 1 and 2 markers yielded 29 genotypes and an HGDI that was only slightly higher than that for panel 2B alone. The HGDIs of panel 2B, panels 2A and 2B, and panels 1 and 2 were all >0.95, which is considered optimal for a typing system (37). These data might indicate that in the setting of a local outbreak investigation, panel 1 and panel 2A could be omitted and the highly polymorphic panel 2B would be sufficient. However, typing with panel 1 and 2A is necessary before typing with panel 2B is used in order to check that the isolate is indeed of the expected local type, as it is important to keep in mind that the HGDI value is not a sufficient indicator for the selection of tandem repeat loci and validation of an MLVA assay. Three of the five panel 2B loci (bruce04, bruce09, bruce30) are part of the hypervariable octameric oligonucleotide fingerprint variable-number tandem repeats (HOOF-Prints) MLVA panel of loci that have been described by Bricker et al. (6, 7) and that have proven to be highly mutable in terrestrial Brucella spp. (6, 7, 36, 41). The most frequent genotype (largest cluster) comprised seven isolates (isolates S3 to S9) obtained over 7 weeks in the late spring and early summer of 2002, demonstrating an outbreak from a common source. Another genotype was shared by isolates S18 and S24, obtained 5 months apart from one patient who had no at-risk occupational history and who experienced a recurrence of his febrile illness. The two isolates from this patient differed from the rest of the isolates at the bruce45 minisatellite locus from panel 1 (panel 1 genotype 60, two versus three repeats), indicating relapsing disease in this patient. Two isolates obtained 2 months apart (isolates S37 and S41) and four pairs obtained over 1 year apart (isolates S25 and S30, S12 and S39, S2 and S23, and S28 and S42) had identical genotypes. These five seemingly unrelated pairs may represent either epidemiologically unrelated isolates with homoplasy at MLVA-16 loci (most likely panel 2B) or persistent circulating strains causing sporadic infections. All other isolates had distinct genotypes reflecting sporadic cases. In a previous study (1), two pairs of epidemiologically related human B. melitensis isolates in Germany displayed identical MLVA-16 profiles. One pair of isolates was obtained from one patient: one isolate was obtained before treatment for brucellosis and one isolate was obtained 3 months later when he experienced a relapse related to inadequate therapy. Another pair of isolates came from a married couple who visited Turkey, where they contracted their infection from the consumption of unpasteurized cheese. Three isolates from a small outbreak of B. suis infections in pigs in Spain also shared the same MLVA-16 genotype (14). The clustering of epidemiologically related isolates identified in the current and previous studies support the use of MLVA-16 as a valuable tool for investigations of outbreaks of both human and animal brucellosis.
![]() View larger version (42K): [in a new window] |
FIG. 1. Dendrogram based on the MLVA-16 genotyping assay showing the relationships of the 42 human isolates of Brucella. Under panels 1, 2A, and 2B are shown the individual MLVA-16 loci and the number of tandem repeats for each isolate. Key, the serial number for the isolate in the Brucella2007 MLVA database (http://mlva-u-psud.fr/) (see Fig. S1 in the supplemental material); strain, the local strain designation; panels 1, 2A, and 2B, genotypes corresponding to each isolate in the database for each set of loci; date isolation, the dates (day/month/year) that samples were received in the laboratory for culturing. The information under comments provides the specimen source-patient gender (M, male; F, female)/patient age (in years).
|
|
View this table: [in a new window] |
TABLE 1. Summary of MLVA-16 results for 42 human isolates of B. melitensis
|
Although phenotypic biovar typing was not performed in this study, previous studies have shown that the results of the MLVA-16 or the MLVA-21 typing schemes that have been devised (1, 41) show no correlation with those of biovar typing for B. melitensis. Our isolates were closely related to B. melitensis isolates from various countries that belonged to all three biovars (biovars 1, 2, and 3). Neither MLVA nor multilocus sequence typing distinguished the B. melitensis biovars, which are mere serotypes (1, 40, 41). This indicates that variable-number tandem-repeat loci and the single-nucleotide polymorphisms which provide congruent data evolve independently of the putative genetic determinants for these biovars. This is also the case for B. abortus, but to a lesser extent. In contrast, both MLVA-16 and MLVA-21 separate B. suis biovar 2 from the other B. suis biovars, biovars 1, 3, and 4 (see Fig. S1 in the supplemental material) (14, 41).
Both MLVA-16 and MLVA-21 segregated Brucella strains into groups corresponding to the recognized phylogenetic lineages in the genus Brucella (22, 41). When it was applied to more than 500 strains of terrestrial Brucella, panel 1 in MLVA-16 subdivided these isolates into 65 genotypes that were, with one exception, exclusive to each of the Brucella species (see Fig. S1 in the supplemental material). Similarly, by MLVA-21, six minisatellites were found to be sufficient for species designation (41). In both schemes, B. canis clustered with B. suis biovar 4. These data suggest that MLVA is a generally useful method for primary species assignment of terrestrial Brucella spp., especially in regions of endemicity where B. melitensis, B. abortus, and B. suis are the predominant causes of both human and animal brucellosis.
The diversity of the panel 1 loci and most panel 2 loci in our series parallels that found in a study of 24 human B. melitensis isolates from Sicily in Italy (24). The two studies are comparable in that they focused on a restricted geographic area. Two panel 1 genotypes (genotypes 49 and 51) were observed in the Sicilian study, and both were from the West Mediterranean group (see Fig. S1 in the supplemental material). We found higher levels of diversity for bruce30 (panel 2B) and bruce19 (panel 2A) and a lower level of diversity for bruce08 (panel 1). The Sicilian isolates formed a homogeneous group that could be distinguished from the East Mediterranean group at three of the panel 1 minisatellite loci. These data illustrate the usefulness of panel 1 not only for the species identification of Brucella isolates but also for inference of their geographical origins. In contrast, panel 2 markers afford a higher discriminatory power for investigation of strain relatedness in regions of endemicity. An abbreviated panel consisting of five microsatellite loci, panel 2B, can be used as a practical short-term epidemiological tool in this setting, but the species cannot be deduced from the data produced by this panel only. Conversely, in regions of nonendemicity or whenever the introduction of an exogenous strain might be suspected, all loci should be used to pinpoint the region where the isolates may have originated. For that reason, there is a need to build expanded searchable and shared international databases of MLVA fingerprints to enable such comparisons. The currently existing Brucella2007 MLVA database (http://mlva.u-psud.fr) provides one constituent within the framework of such a concerted effort (16).
In summary, MLVA-16 displayed a good discriminatory power for the typing of B. melitensis isolates from a small region of endemicity, and either panel 2 or panels 1 and 2 may be used as epidemiological tools for the resolution of strains and for distinguishing relapses from reinfections in patients with brucellosis.
Published ahead of print on 15 October 2008. ![]()
Supplemental material for this article may be found at http://jcm.asm.org/. ![]()
|
|
|---|
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»