- Methodology article
- Open Access
Lab on a chip genotyping for Brucella spp. based on 15-loci multi locus VNTR analysis
BMC Microbiology volume 9, Article number: 66 (2009)
Brucellosis is an important zoonosis caused by the genus Brucella. In addition Brucella represents potential biological warfare agents due to the high contagious rates for humans and animals. Therefore, the strain typing epidemiological tool may be crucial for tracing back source of infection in outbreaks and discriminating naturally occurring outbreaks versus bioterroristic event. A Multiple Locus Variable-number tandem repeats (VNTR) Analysis (MLVA) assay based on 15 polymorphic markers was previously described. The obtained MLVA band profiles may be resolved by techniques ranging from low cost manual agarose gels to the more expensive capillary electrophoresis sequencing. In this paper a rapid, accurate and reproducible system, based on the Lab on a chip technology was set up for Brucella spp. genotyping.
Seventeen DNA samples of Brucella strains isolated in Sicily, previously genotyped, and twelve DNA samples, provided by MLVA Brucella VNTR ring trial, were analyzed by MLVA-15 on Agilent 2100. The DNA fragment sizes produced by Agilent, compared with those expected, showed discrepancies; therefore, in order to assign the correct alleles to the Agilent DNA fragment sizes, a conversion table was produced. In order to validate the system twelve unknown DNA samples were analyzed by this method obtaining a full concordance with the VNTR ring trial results.
In this paper we described a rapid and specific detection method for the characterization of Brucella isolates. The comparison of the MLVA typing data produced by Agilent system with the data obtained by standard sequencing or ethidium bromide slab gel electrophoresis showed a general concordance of the results. Therefore this platform represents a fair compromise among costs, speed and specificity compared to any conventional molecular typing technique.
The genus Brucella comprises Gram-negative bacteria that are the etiological agents of brucellosis, a zoonosis that represents one of the most important worldwide disease affecting animals and humans, especially in the Mediterranean areas [1, 2]. The disease can be transmitted to humans directly by contact with infected animals or indirectly by contaminated dairy products. Brucella infections can cause chronic debilitating diseases in humans with the involvement of multiple organs and low mortality while in the domesticated animals can lead to reproductive failure with subsequent economic loss. Six classically species are recognized within the genus Brucella: Brucella abortus (7 biovars), Brucella melitensis (3 biovars), Brucella suis (5 biovars), Brucella ovis, Brucella canis, and Brucella neotomae . Isolation and characterization of novel different Brucella strains from a wide variety wildlife species from terrestrial and marine mammals, recently classified as three additional species, B. ceti (cetaceans) B. pinnipedialis and B. microti underline the importance of wildlife as reservoir for zoonotic brucellosis [4, 5]. In addition, Brucella spp. represent potential biological warfare agents due to the high contagious rates for humans and animals [6, 7]. Therefore it is very important to have a strain typing epidemiological tool for source trace back in outbreaks of infection. Characterization of Brucella at species and biovar level can be performed by differential microbiological approaches used for phenotyping . However, these biotyping methods are time consuming and potentially hazardous for laboratory personnel. In addition, the limited variation among some species and biovars may result in conflicting data and complicate interpretation [9, 1]. Biotyping methods are not therefore suitable for epidemiological tracking where a more accurate identification is necessary. Thus, genetic characterization using molecular DNA technology has been investigated and several molecular techniques for subtyping have been proposed. In the last years PCR-based assays have been employed in an attempt to find suitable DNA markers for the molecular typing of Brucella. Examples are repetitive intergenic palindromic sequence-PCR (REP-PCR) , random amplified polymorphic DNA-PCR (RAPD-PCR) , arbitrary primed-PCR (AP-PCR) , amplified fragment length polymorphism (AFLP) , single nucleotide polymorphism (SNP) [14, 15], and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) . Recently, a highly discriminatory typing system called MLVA, based on multilocus characterization of variable number of tandem repeats (VNTRs), has been proposed for several pathogens such as Haemophilus influenzae , Mycobacterium tuberculosis , Bacillus anthracis and Yersinia pestis [19–21]. The multilocus VNTR analysis methods (MLVA-15 and -16) were also described for Brucella genotyping [22, 23, 6]. The MLVA-15 is based on 15 polymorphic markers subdivided in two panels: the panel 1, consisting of 8 minisatellite markers with a good species identification capability, and the panel 2, consisting of 7 microsatellite markers with higher discriminatory power. MLVA-15 assay generates multiple band profiles or fingerprint patterns and represents a useful molecular approach for Brucella typing in spite of the high DNA homology among Brucella species (90%) . The obtained MLVA band profiles may be resolved by different techniques ranging from low cost manual agarose gels to the more expensive capillary electrophoresis sequencing systems. Recently, a rapid and inexpensive method based on the Lab on a chip technology (Agilent Technologies) has been proposed. This miniaturized platform for electrophoresis applications is able to size and quantitate PCR fragments, and was previously used for studying the genetic variability of bcIA gene of Bacillus cereus  and for genotyping of Bacillus anthracis and Yersinia pestis . In this paper we evaluate the possibility to use the Agilent 2100 Bioanalyzer equipment for MLVA typing of Brucella strains using the selected subset of 15 loci proposed by La Flèche . On that account, we analyzed by Agilent technology the collection of seventeen human Brucella isolates whose MLVA fingerprinting profiles were previously resolved by capillary electrophoresis sequencing system  and results compared, while twelve DNA samples, provided in 2007 for a MLVA VNTR ring trial, were de novo genotyped.
In order to set up this method, we analyzed, by MLVA-15  on Agilent 2100 Bioanalyzer, seventeen Brucella strains previously identified as B. melitensis biovar 3 by direct sequencing . The 15 VNTR markers amplified for bacterial genotyping were arranged into three duplex and nine singleplex PCR reactions. The arrangement of different loci in the same multiplex was such to avoid overlapping of VNTR markers size ranges; in this way it was possible to use a single 12-wells chip (DNA 1000 LabChip Kit) for strain analysis. After PCR amplification 1 μl of each reaction was loaded into the chip and the amplification product size estimates were obtained by the Agilent Software. These values showed discrepancies compared with those expected , making difficult the allele assignment directly from rough data. In order to assign the correct alleles to the Agilent DNA fragment sizes, a conversion table containing for each locus the expected size, the range of observed sizes, including arithmetical average ± standard deviation, and the corresponding allele was produced (Additional File 1). We could establish experimentally the variability range for each allele. Even if we didn't conduct an extensive study on migration, these ranges were determined considering interchip/intrachip variability from the same amplification product or different amplification of the same strain allele or, for the same alleles, different strain amplification (data not shown). The data are considered valuable only if standard deviation is lower than the 50% of the repeat unit length. In this way, we could measure a variation proportionated to the relative number of nucleotides in each repeat unit. All the allele measurements satisfied this criterion, allowing the unambiguous assignment of the correct allele to each observed value (Additional File 1). In order to validate this platform, we analyzed twelve unknown samples provided by Dr. Falk Melzer for MLVA Brucella 2007 ring trial . The Agilent 2100 Bioanalyzer MLVA-15 products were separated and DNA fragment sizes were correlated to the alleles by the conversion table. The resulting fingerprint for each strain was matched against the MLVA database Brucella test , allowing identification of samples and their genetic relationship with the other database strains (Table 1). The identified species were compared with the VNTR ring trial results , obtaining a full concordance.
The renewed threat of biological weapons and the appearance of zoonotic infections caused by Brucella spp. have reinforced the need to provide a rapid and specific detection method for the characterization of Brucella isolates. In fact, the discrimination between natural outbreaks and/or intentional release of micro-organism agents may be of crucial importance in the context of the bioterrorism. Brucella strains differentiation in the last years was obtained by assays based on the analysis of the polymorphism of tandemly repeated DNA sequences. This strategy has proven to be appropriate for pathogenic bacterial species typing with a high genetic homogeneity. Recently a new multilocus VNTR analysis (MLVA-15) method, based on a subset of fifteen tandem repeat loci, was described demonstrating high discriminatory power . The different alleles, amplified by standard PCR techniques, can be analyzed by several electrophoretic techniques that can be time-consuming as agarose gel, or require high costs for reagents and instruments as capillary electrophoresis sequencing. Our attention was addressed on the Agilent 2100 "Lab on a Chip" platform, in order to obtain a more rapid and inexpensive method of discrimination. The loading of the fifteen amplification products in a single chip provides an increased time-reduction (chip run time is approximately 30 minutes) as compared to assay run on agarose gels for the same analysed markers. After data collection and analysis, we observed that the size estimated by the Agilent 2100 Bioanalyzer were shifted by a variable value (offset) in respect to the real size. Indeed, the estimated sizes were shifted of a variable value (offset) respect to the real size estimated by sequencing and each offset showed a variable range. These differences could be due to the different nature of the gel matrix or to a slightly biased flanking sequence or repeat unit specific mobility pattern . These discrepancies between observed and expected sizes have been overcome creating a correspondence table which allows for each observed value to assign the expected size and corresponding allele (Additional File 1). The comparison of the results obtained by the multilocus VNTR analysis (MLVA-15) method  on the Agilent 2100 "Lab on a Chip" platform showed a good size correlation with nucleotide sequence analysis, confirming the rapidity and efficiency of this platform respect to standard sequencing or ethidium bromide slab gel electrophoresis. Furthermore, the possibility to compare different chip data in different times made the system suitable for epidemiological purposes. We consider this system one the most promising platforms for MLVA-15 typing because offers a fair compromise among costs, speed and specificity compared to any of the conventional molecular typing techniques.
In this paper is described a rapid, accurate and reproducible system for genotyping of Brucella. The method is based on Multiple Locus Variable Number Tandem Repeats (VNTR) Analysis (MLVA) assay with 15 markers (MLVA-15), previously described, and the Agilent 2100 Bioanalyzer, for the separation of nucleic acid molecules. The validation of the system was performed reanalyzing seventeen DNA samples, previously genotyped by traditional techniques, and twelve DNA provided for a VNTR ring trials, obtaining a general concordance of the results. Therefore MLVA typing data produced by Agilent system represents an alternative to standard sequencing or ethidium bromide slab gel electrophoresis.
Brucella strains and DNA extraction
In this study, seventeen Brucella strains isolate from Sicilian hospitalized patients with acute brucellosis , and twelve DNA samples, provided by Dr. Falk Melzer for the Ring trial Brucella 2007, were analysed. DNA was extracted using proteinase K and sodium dodecyl sulfate method. Pellets were resuspended in 50 μl of nuclease-free water. Twenty nanograms of DNA template were used for PCR amplifications.
VNTR amplifications were performed according to Le Flèche et al . Fifteen sets of primers previously proposed were used: Bruce06, Bruce08, Bruce11, Bruce12, Bruce42, Bruce43, Bruce45, Bruce55 (panel 1), and Bruce04, Bruce07, Bruce09, Bruce16, Bruce18, Bruce21 and Bruce30 (panel 2). The 15 markers were arranged into 3 duplex, indicated as multiplex 1a, 2b and 3c respectively for the loci bruce 43 and bruce 08, bruce 12 and bruce 18 and bruce 11 and bruce 21 and 9 singleplex. Amplification reaction mixtures were prepared in 15 μl volumes using 1U FastStart polymerase Taq (Roche) and containing 1 ng of DNA, 1× PCR Roche reaction buffer (10 mM Tris-HCl, 2,5 mM MgCl2, 50 mM KCl pH 8.3), 0.2 mM dNTPs (Roche) and 0,3 μM of each flanking primer. Thermal cycling, conducted on a Peltier Thermal Cycler DNA Engine DYAD (MJ Research), was performed as follows: The optimized protocol was, after an initial heating at 95°C for 5 min, 35 cycles denaturation at 95°C for 30 sec, annealing at 60°C for 30 sec and extension at 70°C for 60 sec. A final extension was performed at 70°C for 5 min.
The amplification products were loaded into chip wells prepared according to manufacturer recommendations (DNA 1000 LabChip Kit). Each chip contains 16 wells: 12 for the samples, 3 for gel mix. After gel preparation, each sample well was loaded with 1 μl of PCR reaction and 5 μl of internal marker (containing two MW size standards of 15 and 1500 bp). One microliter of DNA ladder was loaded in the ladder well. Finally, the chip was vortexed for 60 sec and inserted into Agilent 2100 Bioanalyzer. During the separation of the fragments, the samples were analyzed sequentially and electropherograms, virtual gel images and table data were shown. Amplification product size estimates were obtained by using the Agilent 2100 Expert Software version B.02.03.SI307 firmware C.01.055 (Agilent Technologies).
Corbel MJ: Brucellosis: an overview. Emerg Infect Dis. 1997, 3: 213-21.
Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos EV: The new global map of human brucellosis. Lancet Infect Dis. 2006, 6: 91-99. 10.1016/S1473-3099(06)70382-6.
Corbel MJ, Brinley-Morgan WJ: Genus Brucella Meyer and Shaw 1920, 173AL. Bergey's Manual of Systematic Bacteriology 1984. Edited by: Krieg NR, Holt JG. 1984, Baltimore: Williams and Wilkins, 1: 377-390.
Foster G, Osterman BS, Godfroid J, Jacques I, Cloeckaert A: Brucella ceti sp. nov. and Brucella pinnipedialis sp. nov. for Brucella strains with cetaceans and seals as their preferred hosts. Int J Syst Evol Microbiol. 2007, 57: 2688-2693. 10.1099/ijs.0.65269-0.
Scholz HC, Hubalek Z, Sedlácek I, Vergnaud G, Tomaso H, Al Dahouk S, Melzer F, Kämpfer P, Neubauer H, Cloeckaert A, Maquart M, Zygmunt MS, Whatmore AM, Falsen E, Bahn P, Göllner C, Pfeffer M, Huber B, Busse HJ, Nöckler K: Brucella microti sp. nov., isolated from the common vole Microtus arvalis. Int J Syst Evol Microbiol. 2008, 58: 375-382. 10.1099/ijs.0.65356-0.
Al Dahouk S, Le Fleche P, Nockler K, Jacques I, Grayon M, Scholz HC, Tomaso H, Vergnaud G, Neubauer H: Evaluation of Brucella MLVA typing for human brucellosis. J Microbiol Methods. 2007, 69: 137-145. 10.1016/j.mimet.2006.12.015.
Whatmore AM, Perrett LL, MacMillan AP: Characterization of the genetic diversity of Brucella by multilocus sequencing. BMC Microbiol. 2007, 7: 34-10.1186/1471-2180-7-34.
Alton GG, Jones LM, Angus RD, Verger JM: Techniques for the brucellosis laboratory. 1988, Institut National de la Recherche Agronomique, Paris, France
Banai M, Mayer I, Cohen A: Isolation, identification, and characterization in Israel of Brucella melitensis biovar 1 atypical strains susceptible to dyes and penicillin, indicating the evolution of a new variant. J Clin Microbiol. 1990, 28: 1057-1059.
Tscherneva E, Rijpens N, Naydensky C, Herman LMF: Repetitive element sequence based polymerase chain reaction for typing of Brucella strains. Vet Microbiol. 1996, 51: 169-178. 10.1016/0378-1135(96)00036-3.
Tscherneva E, Rijpens N, Jersek B, Herman LMF: Differentiation of Brucella species by random amplified polymorphic DNA analysis. J Appl Microbiol. 2000, 88: 69-80. 10.1046/j.1365-2672.2000.00945.x.
AlMomin S, Saleem M, Al-Mutawa Q: The use of an arbitrarily primed PCR product for the specific detection of Brucella. World Journal of Microbiology & Biotechnology. 1999, 15: 381-385. 10.1023/A:1008975616874.
Whatmore AM, Murphy TJ, Shankster S, Young E, Cutler S, Macmillan AP: Use of amplified fragment length polymorphism to identify and type Brucella isolates of medical and veterinary interest. J Clin Microbiol. 2005, 43: 761-769. 10.1128/JCM.43.2.761-769.2005.
Marianelli C, Ciuchini F, Tarantino M, Pasquali P, Adone R: Molecular characterization of the rpoB gene in Brucella species: new potential molecular markers for genotyping. Microbes Infect. 2006, 8: 860-865. 10.1016/j.micinf.2005.10.008.
Scott JC, Koylass MS, Stubberfield MR, Whatmore AM: Multiplex Assay based on single-nucleotide polymorphisms for rapid identification of Brucella isolates at the species level. Appl Environ Microbiol. 2007, 73: 7331-7337. 10.1128/AEM.00976-07.
Al Dahouk S, Tomaso H, Prenger-Berninghoff E, Splettstoesser WD, Scholz HC, Neubauer H: Identification of Brucella species and biotypes using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Crit Rev Microbiol. 2005, 31: 191-196. 10.1080/10408410500304041.
Schouls LM, Ende van der A, Pol van de I, Schot C, Spanjaard L, Vauterin P, Wilderbeek D, Witteveen S: Increase in genetic diversity of Haemophilus influenzae serotype b (Hib) strains after introduction of Hib vaccination in The Netherlands. J Clin Microbiol. 2005, 43: 2741-2749. 10.1128/JCM.43.6.2741-2749.2005.
Le Flèche P, Fabre M, Denoeud F, Koeck JL, Vergnaud G: High resolution, on-line identification of strains from the Mycobacterium tuberculosis complex based on tandem repeat typing. BMC Microbiol. 2002, 2: 37-10.1186/1471-2180-2-37.
Keim P, Price LB, Klevytska AM, Smith KL, Schupp JM, Okinaka R, Jackson PJ, Hugh-Jones ME: Multiple-Locus Variable-Number Tandem Repeat Analysis reveals genetic relationships within Bacillus anthracis. J Bacteriol. 2000, 182: 2928-2936. 10.1128/JB.182.10.2928-2936.2000.
Le Flèche P, Hauck Y, Onteniente L, Prieur A, Denoeud F, Ramisse V, Sylvestre P, Benson G, Ramisse F, Vergnaud G: A tandem repeatsdatabase for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis. BMC Microbiol. 2001, 1: 2-10.1186/1471-2180-1-2.
Lista F, Faggioni G, Samina Valjevac S, Ciammaruconi A, Vaissaire J, le Doujet C, Gorgé O, De Santis R, Carattoli A, Ciervo A, Fasanella A, Orsini F, D'Amelio R, Pource C, Cassone A, Vergnaud G: Genotyping of Bacillus anthracis strains based on automated capillary 25-loci Multiple Locus Variable-Number Tandem Repeats Analysis. BMC Microbiology. 2006, 6: 33-10.1186/1471-2180-6-33.
Kattar MM, Jaafar RF, Araj GF, Le Flèche P, Matar MG, Rached RA, Khalife S, Gilles Vergnaud G: Evaluation of a Multilocus Variable-Number Tandem-Repeat Analysis Scheme for Typing Human Brucella Isolates in a Region of Brucellosis Endemicity. J Clin Microbiol. 2008, 45: 3935-3940. 10.1128/JCM.00464-08.
Le Flèche P, Jacques I, Grayon M, Al Dahouk S, Bouchon P, Denoeud F, Nöckler K, Neubauer H, Guilloteau LA, Vergnaud G: Evaluation and selection of tandem repeat loci for a Brucella MLVA typing assay. BMC Microbiology. 2006, 6: 9-10.1186/1471-2180-6-9.
Whatmore AM, Shankster SJ, Perrett LL, Murphy TJ, Brew SD, Thirlwall RE, Cutler SJ, MacMillan AP: Identification and characterization of Variable-Number Tandem-Repeat Markers for typing of Brucella spp. J Clin Microbiol. 2006, 44: 1982-1993. 10.1128/JCM.02039-05.
Castanha ER, Swiger RR, Senior B, Fox A, Waller LN, Fox K: Strain discrimination among B. anthracis and related organisms by characterization of bclA polymorphisms using PCR coupled with agarose gel or microchannel fluidics electrophoresis. J Microbiol Methods. 2006, 64: 27-45. 10.1016/j.mimet.2005.04.032.
Ciammaruconi A, Grassi S, De Santis R, Faggioni G, Pittiglio V, D'Amelio R, Carattoli A, Cassone A, Vergnaud G, Lista F: Fieldable genotyping of Bacillus anthracis and Yersinia pestis based on 25-loci Multi Locus VNTR Analysis. BMC Microbiology. 2008, 8: 21-10.1186/1471-2180-8-21.
Marianelli C, Graziani C, Santangelo C, Xibilia MT, Imbriani A, Amato R, Neri D, Cuccia M, Rinnone S, Di Marco V, Ciuchini F: Molecular epidemiological and antibiotic susceptibility characterization of Brucella isolates from humans in Sicily, Italy. J Clin Microbiol. 2007, 45: 2923-2928. 10.1128/JCM.00822-07.
Brucella MLVA genotyping. [http://mlva.u-psud.fr/BRUCELLA/spip.php?article74&var_recherche=ring%20trial]
MLVAbank for Bacterial Genotyping. [http://mlva.u-psud.fr/]
This work was part of the European Biodefence project CEPA13.14 involving biodefence institutions from Sweden, Norway, the Nederlands, Germany, France and Italy. We especially thank Dr Gilles Vergnaud and Dr Falk Melzer for the organization of the Ring trial Brucella 2007 very helpful for this work.
RDS and AC did the set up of the Brucella 15-MLVA assay. RDS, AC and CM participated to typing work. FL, RD'A and CM did the error checking analysis. RDS and GF did various sequence analysis. FL and RDS were in charge of the database and clustering analyses. FL, AC, RD'A and RDS conceived the study. FL and RDS wrote the report. All authors read and approved the final manuscript
Electronic supplementary material
Additional file 1: Comparison between Brucella product sizes inferred by Agilent 2100. Bioanalyzer software – Observed size and their arithmetic average (x) ± standard deviation (σ) – and actual sizes obtained by direct sequencing of the PCR product or data available in Genbank (Expected size). Unit Length size (UL bps). (DOC 258 KB)
About this article
Cite this article
De Santis, R., Ciammaruconi, A., Faggioni, G. et al. Lab on a chip genotyping for Brucella spp. based on 15-loci multi locus VNTR analysis. BMC Microbiol 9, 66 (2009). https://doi.org/10.1186/1471-2180-9-66
- Variable Number Tandem Repeat
- Yersinia Pestis
- Brucella Abortus
- Brucella Melitensis