Genome-wide unique insertion sequences among five Brucella species and demonstration of differential identification of Brucella by multiplex PCR assay

Brucellosis is a neglected zoonotic disease caused by alpha proteobacterial genus Brucella comprising of facultative intracellular pathogenic species that can infect both animals and humans. In this study, we aimed to identify genome-wide unique insertion sequence (IS) elements among Brucella abortus, B. melitensis, B. ovis, B. suis and B. canis for use in species differentiation by conducting an intensive in silico-based comparative genomic analysis. As a result, 25, 27, 37, 86 and 3 unique ISs were identified respectively and they had a striking pattern of distribution among them. To explain, a particular IS would be present in four species with 100% identity whereas completely absent in the fifth species. However, flanking regions of that IS element would be highly identical and conserved in all five species. Species-specific primers designed on these flanking conserved regions resulted in two different amplicons grouping the species into two: one that possesses IS and the other that lacks it. Seeking for species-specific amplicon size for particular species was sufficient to identify it irrespective of biovar. A multiplex PCR developed using these primers resulted in successful differentiation of the five species irrespective of biovars with significant specificity and sensitivity when examined on clinical samples.

other 16 . Thus, the stability of IS711 was not only relevant for Brucella typification: its mobility was implicated in the generation of genetic diversity and speciation, as shown by the distribution of IS711 among Brucella species. However, in 2011, Mancilla and coworkers reported the mobility of IS711 in B. abortus field strains under natural conditions, causing polymorphism within the strains of B. abortus. The study reported B. abortus strains isolated from cattle carrying seven copies of complete IS711 instead of six copies, and this new copy was found to be located in a completely new locus that was not previously reported to harbor IS711. The reason for the presence of this one additional copy of IS was predicted to be replicative transposition within the genome 15 . However, in another study, novel insertion of IS711 within the Omp31 gene of an atypical B. ovis isolate recovered from an infected ram in Hungary was encountered and was not correctly identified in the Bruce-ladder multiplex PCR assay 17 . Based on these new case reports, the fact that the most acclaimed IS element, IS711, could be a transposon element in Brucella species was established, and its reliability in differentiating Brucella species was weakened as it demonstrated polymorphism within the strains of a species. With this, the diagnostic application of IS711 reveals the need for a more improved approach to decipher polymorphisms among Brucella species that is even more unique and reliable in differentiating species, even at the biovar level, if possible.
In the current study, we describe a novel in silico-based comparative genomic analysis approach to decipher species-associated unique and precise DNA block distributions among five species of Brucella, namely, B. abortus, B. melitensis, B. ovis, B. suis and B. canis irrespective of biovars. B. abortus, B. melitensis, B. suis and B. canis were selected as representative zoonotic species of Brucella, and B. ovis was selected as a representative species of non-zoonotic Brucella. Due to the unavailability of the remaining seven species in our repository, they were not included in this study. Preliminary in silico characterization of the conserved flanking regions indicated that the DNA blocks were inserted in transposition sites, indicating that the DNA blocks were actually IS elements. This was achieved by conducting an intensive in silico-based comparative genomic analysis of Brucella species using the Perl language and the vast data available online on whole genome sequences of Brucella species. The distribution of the precise stretch of DNA blocks among the five species was unique, as a stretch of nucleotide sequence was embedded within a conserved region of the genome in four of five species but absent in the fifth species. However, the conserved flanking regions (where the primer pairs are located) demonstrated 100% homology in all five species. Furthermore, these highly specific nucleotide sequences were evaluated for their ability to detect targets, and a novel multiplex PCR was developed wherein the five species were differentiated irrespective of biovars.

Materials and Methods
Bacterial strains. The Brucella strains (Table 1) were clinical isolates preserved in the repository of the Defence Food Research Laboratory (DFRL). The clinical strains recovered from infected animals and patients were identified by biochemical tests and DNA sequencing. Reference strains of Pseudomonas aeruginosa, P. azotogenesis, Aeromonas hydrophila, Proteus mirabilis, Citrobacter freundii and Yersinia enterocolitica were used as negative controls in the specificity study. All bacterial strains were aseptically recovered from glycerol stocks and grown in trypticase soy broth (TSB) overnight at 37 °C. Genomic DNA from harvested cultures was prepared using a QIAamp DNA Mini Kit (Qiagen, Germany) according to the manufacturer's instructions. DNA samples were determined using a Nanodrop TM 1000 spectrophotometer (Thermo Scientific, MA, USA), stored at −20 °C and used in PCR.
Biosafety procedures. All procedures using Brucella strains were performed in a biosafety level 3 laboratory. Procedures for culture of non-Brucella species and DNA extraction were conducted in a class II type A2 biological cabinet (Thermo Scientific, USA).  Table 1). Briefly, the whole genome sequences of the selected Brucella strains were cut into fragments of 1000 bp by implementation of a script written in the Perl programming language. The nucleotide fragments of each species (for instance, B. abortus) were aligned against the genomic sequences of all other biovars within the species using the BLASTN program. Nucleotide fragments that matched all biovars with an E-value less than 10 −200 were identified as highly specific targets for B. abortus irrespective of biovars (Fig. 1). The procedure was repeated with WGSs of B. melitensis, B. suis, B. canis and B. ovis. Thus, the conserved DNA fragments from each species of Brucella were determined and further aligned against the genomic sequences of the other four species to evaluate their specificity.

Genomic analysis of
Analysis of Brucella species-specific DNA fragments. The determined species-specific fragments of Brucella were analyzed for conservation and variation by multiple sequence alignment using the CLUSTALW tool available on the EMBL-EBI website (http://www.ebi.ac.uk/Tools/msa/clustalw2/) and by BLAST analysis against the vast WGS database available in the National Center for Biotechnology Information database (http:// www.ncbi.nlm.nih.gov/) and the Broad Institute (http://www.broadinstitute.org). The respective protein(s) or peptide(s) encoded by each of the species-specific fragments of Brucella were determined using Gene Runner 3.0 software (http://www.generunner.net/). Thus, the determined protein(s) or peptide(s) sequences were analyzed by BLASTP analysis of the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/). Protein sequence features, including the identification of conserved domains, protein superfamily and further protein classification, were performed using the National Center for Biotechnology Information.
Primer design. Primers for specific PCR amplification of Brucella species-specific fragments were designed using Gene Runner 3.0 software (http://www.generunner.net/). The primers were designed to have T m values ranging from 52-68 °C and ΔG values for primer duplex below −10 kcal/mol and not to form hairpin loops or primer dimers. The strategy of designing primers is illustrated in Fig. 2. The primers were designed with the intent of obtaining specific yet easily distinguishable and interpretable amplicons under identical PCR conditions and were thus usable in the development of a multiplex PCR system. Designed primers were analyzed for specificity by using the BLASTN program prior to synthesis. To prevent false negative results due to the interference of inhibitory substances in test samples, competitive internal amplification control (IAC) was generated using the pUC19 plasmid as the target. The IAC primers were designed such that they had 5′ overhanging ends identical to the primer sequences of the B. abortus-specific primer pair and 3′ ends complementary to the pUC19 plasmid, amplifying 808 bp of the pUC19 plasmid and yielding a product of 847 bp (inclusive of the 5′ B. abortus-specific flanking regions of the primer pair). The sequences of the designed primers are presented in Table 2. All primers were synthesized by Eurofins India, India.
Monoplex PCR. For initial standardization of the detection system, monoplex PCR was performed using the primer sets and respective genomic DNA as template. Each 25 µL PCR contained the following: 1 U of Taq polymerase, 250 µM each deoxynucleotide triphosphate, 1x PCR buffer (including 1.5 mM MgCl 2 ), 10 pmol of each oligonucleotide primer, and 25 ng of template DNA. Thermal cycling was carried out in an Eppendorf thermal cycler (Eppendorf, USA) for 35 cycles of 30 s at 94 °C, annealing for 30 s at 56 °C, and extension at 72 °C for 1 min, with a final 8 min extension at 72 °C. Approximately 2 µl of each PCR product was visualized by agarose gel electrophoresis. The template DNA of each species was used with its respective specific primer in each monoplex PCR.
IAC detection limit. The 847 bp IAC amplicon amplified using chimeric oligonucleotide primers was purified using a QIAquick PCR Purification Kit (Qiagen, Germany) and further subjected to 1000-fold serial dilutions from 1:1000 to 1:10000. From each dilution, 1 µl of the amplicon was used as template along with 25 ng of B. abortus gDNA and PCR-amplified in a single reaction using B. abortus-specific primers. The PCR components and amplification conditions were the same as described in the previous methods section. The lowest dilution of www.nature.com/scientificreports www.nature.com/scientificreports/ IAC DNA resulting in an amplicon of the same intensity as B. abortus DNA was selected for standardization of multiplex PCR.
Multiplex PCR assay. Twenty-five nanograms of genomic DNA from each Brucella species, namely, B. abortus, B. melitensis, B. canis, B. suis and B. ovis, along with a 1:9000 dilution of IAC DNA, was subjected to each multiplex PCR in a 50 µl reaction mixture containing 1x PCR buffer, 1.5 mM MgCl 2 , 200 mM each deoxynucleotide triphosphate, 10 pmol of each oligonucleotide primer and 1 U of Taq polymerase (Sigma Aldrich, India). PCR amplification was carried out as follows: initial denaturation at 94 °C for 10 min; 30 cycles of 94 °C for 1 min, 56 °C for 1 min and 72 °C for 1.30 min; and a final elongation step at 72 °C for 8 min. Amplicons were loaded onto a 2% agarose gel containing ethidium bromide that was run at 100 V for 1 h. To evaluate the reproducibility of the developed assay, the PCR was repeated three times in three different thermal cyclers: Mastercycler Gradient (Eppendorf, USA), C1000 Touch TM (Biorad, USA) and Techne TC-3000G (Techne, UK).

Validation of the multiplex PCR assay for specificity.
To evaluate the specificity of the multiplex PCR assay, purified genomic DNA isolated from all bacteria listed in Table 1 was employed. Nucleic acid was isolated using the QIAamp DNA Mini Kit (Qiagen) by following the manufacturer's instructions. Ten microliters of the eluted DNA was used to perform a multiplex PCR assay.
Statistical analysis. Descriptive statistics were performed to evaluate the specificity of the multiplex PCR assay. Reproducibility of the results for multiplex PCR optimization and overall specificity percentage was examined by repeating the blind experiments three times. Kappa statistical analysis was used to measure false negative results. Kappa estimates were defined as follows: values <0.2 imply poor agreement, values from 0.2-0.6 imply fair to moderate, and values >0.6 imply good or very good agreement 18 .   Table 2).
Ethical statement. This article does not contain any studies with human participants or animals performed by any of the authors. Brucella abortus specific fragment. The B. abortus-specific fragment of 1154 bp and its respective rescue fragment of 1478 bp from the four other Brucella species of interest were found to be comprised of an unknown conserved region 299 bp from the forward primer position followed by a conserved hypothetical protein encoding region in an inverted position ( Supplementary Fig. 6). Interestingly, no features of insertional sequences, including direct repeats, terminal inverted repeats, translational frame shifting or translational termination, were observed.

Mining of
Brucella melitensis-specific fragment. The B. melitensis-specific fragment of 745 bp was comprised of a C-terminal portion (including a stop codon) of an unknown gene encoding a hypothetical protein, an unannotated region 223 bp length and a portion (366 bp) of the gene encoding Omp31, including its open reading frame. In the case of its respective rescue fragment of 977 bp from the other four Brucella species, an insertion of a truncated stretch of the ORF of a gene encoding the N-terminal portion of Omp31b within the ORF of the Omp31 gene was observed. Integration of this insertional sequence was accompanied by the duplication of a 14 bp direct target repeat upstream of the insertional sequence ( Supplementary Fig. 7).
Brucella suis-specific fragment. The B. suis-specific fragment of 290 bp and its respective rescue fragment of 383 bp from the four other Brucella species of interest were identified as a portion of the gene encoding flagellum-specific ATP synthase, fliI, belonging to the flagellum-specific ATPase/type III secretory pathway virulence-related protein. ATPases of this group are responsible for the export of flagellum and virulence-related proteins. The bacterial flagellar motor is similar to the F0F1-ATPase, in that they both are proton-driven rotary molecular devices. However, the main function of the bacterial flagellar motor is to rotate the flagellar filament for cell motility. Intracellular pathogens such as Salmonella and Chlamydia also have proteins that are similar to www.nature.com/scientificreports www.nature.com/scientificreports/ the flagellar-specific ATPase but function in the secretion of virulence-related proteins via the type III secretory pathway. The gene was found to be directional within the genome of the Brucella species. Similar to the case of B. abortus-specific and rescue fragments, no features of insertional sequences, including direct repeats, terminal inverted repeats, translational frame shifting or translational termination, were observed ( Supplementary Fig. 8).
Brucella canis-specific fragment. The B. canis-specific fragment of 224 bp and its respective rescue fragment of 521 bp from four other Brucella species of interest were identified to be the structural domain of the acyl-CoA acyltransferase enzyme from GNAT family N-acetyltransferase. This domain has a 3-layer alpha/beta/ alpha structure that contains mixed beta-sheets and can be found in proteins including N-acetyl transferase (NAT) family members, aminoglycoside N-acetyltransferases, diamine acetyltransferase 1, autoinducer synthetases such as protein LasI and acyl-homoserinelactone synthase EsaI, the histone acetyltransferase domain of P300/CBP-associating factor PCAF and the catalytic domain of GCN5 histone acetyltransferase. The insertional sequence had a length of 298 bp, and the structural domain of this gene encoding acyl-CoA acyltransferase was also annotated to be a part of this structural domain inserted between flanking direct "AT" repeats at the 128 th and 426 th nucleotide positions of the rescue fragment ( Supplementary Fig. 9).
Brucella ovis-specific fragment. The B. ovis-specific fragment of 446 bp comprised a C-terminal portion (336 bp) of phosphoglycerate kinase encoded by the pgk gene, a 33 bp DNA block encoding transposase, a 9 bp direct repeat, a 35 bp DNA block encoding DNA translocase and a 33 bp DNA sequence encoding the N-terminus of a peptidase. On the other hand, the rescue fragment of the same conserved stretch in the other four Brucella species contained a 150 bp DNA block comprising an 87 bp DNA fragment encoding an outer membrane autotransporter barrel domain, 24 bp of an unknown nucleotide sequence and a 39 bp truncated stretch of DNA sequence encoding a portion of a protein belonging to the ABC transporter ATP-binding protein group. The insertion of this 150 bp DNA block was found to be after the direct repeat observed in the specific fragment ( Supplementary Fig. 10).

PCR assay optimization.
For initial standardization of the detection system, monoplex PCR was performed using each primer pair and respective genomic DNA as template. The result of each monoplex PCR is represented in Fig. 3. Monoplex PCR using B. abortus, B. melitensis, B. suis, B. ovis and B. canis in the presence of their specific primer pairs yielded amplicons of 1154 bp, 745 bp, 290 bp, 446 bp and 224 bp, respectively. Sequencing of each PCR amplicon revealed the complete identity of the amplified products with their specific genes. When the reaction was performed using genomic DNA from all five Brucella species in a single PCR test with any one pair of primers as mentioned in the methods section, each primer pair amplified a specific amplicon of a particular size with its respective species and a rescue amplicon with the remaining four species that was longer in length when compared to the respective specific amplicons; therefore, the amplicons were grouped into separate group altogether (Fig. 4). The rescue amplicons would provide additional information on the presence of other Brucella species in the test sample and hence rules out the chance of having nonspecific amplification with other contaminating organisms. For instance, the primer specific for B. abortus resulted in a 1154 bp amplicon with B. abortus and a 1478 bp amplicon with four other species, namely, B. melitensis, B. suis, B. ovis and B. canis. The multiplex PCR described in this study includes a total of six pairs of primers: one pair for each of the five Brucella species, namely, B. abortus, B. melitensis, B. canis, B. suis and B. ovis, and one pair for the internal amplification control. In our multiplex PCR assay, each primer pair had the capability to detect a specific Brucella species and differentiate it from the other four species. This species-dependent differential amplification was noticed in the case of all specific primers, wherein the primer-specific species would result in an amplicon different from the amplicon size yielded by the remaining four species.
Interpretation of a multiplex PCR system. The multiplex PCR described here includes a total of six pairs of primers: one pair for amplification of each of the five Brucella species, namely, B. abortus, B. melitensis, B. canis, B. suis and B. ovis, and one for the internal amplification control (Fig. 5). Each pair of Brucella species primers mentioned above can yield PCR amplicons for all five species, resulting in a specific amplicon of a particular size with its respective species and a rescue amplicon longer than the specific amplicon with the other four species, thus differentiating them from the fifth Brucella species. For example, the B. abortus-specific primer pair results in a 1154 bp amplicon with B. abortus and a 1478 bp amplicon with 4 other species, namely, B. melitensis, B. suis, B. ovis and B. canis. The interpretation of the multiplex PCR assay described here is illustrated in Fig. 6. In an instance where multiplex PCR is performed with a test sample contaminated with all five Brucella species, namely, B. abortus, B. melitensis, B. suis, B. canis and B. ovis; the result will be as follows: (1)

Discussion
The characteristic feature of bacterial ISs is their capacity to generate mutations in an organism by transposition within the genome. ISs occur in various copy numbers in a particular genome and can move within the genome or even between genomes horizontally, resulting in clonal expansion of species in a genus 19 . Intensive studies with respect to insertion sequences have unraveled genetic features, namely, direct repeats, terminal inverted repeats, translational frame shifting or translational termination, segmental genome duplication and several other features associated with genetic insertion and/or translocation of DNA blocks in bacterial genomes 20 .
The transposition of ISs is generally identified by the presence of short terminal inverted repeats that flank these DNA blocks. Many ISs also carry one or two open reading frames encoding genes that are essential for their transposition. The mechanism of transposition is primarily dependent on catalysts, including transposase and inverted repeats necessary for the binding of transposase, which in turn results in donor DNA cleavage and strand  www.nature.com/scientificreports www.nature.com/scientificreports/ transfer 20 . The presence and transposition of IS elements within and between genomes are associated with gene inactivation, modulating virulence, the expression of neighboring genes, antibiotic resistance and metabolism 21 . Irrespective of IS families, elevated transposition of ISs between genomes serves to amplify both deleterious  B. abortus, B. melitensis, B. ovis, B. suis   www.nature.com/scientificreports www.nature.com/scientificreports/ and advantageous fitness changes, but this accelerates the extinction of elements, particularly if the mean fitness effect of IS-induced mutations is negative. However, transposition bursts create genetic diversity, which occasionally generates advantageous mutations and thus help organisms adapt to new environments 22 . In this line of IS-mediated mutations among prokaryotes, several bacterial models, including Enterobacteriaceae members, Pseudomonas sp., Acinetobacter baumannii, Lactobacillus sp., Burkholderia sp., Staphylococcus aureus and Vibrio sp. have been successfully explored and studied. These studies have supported the establishment of a knowledge base on the involvement of ISs in the plasticity and adaptability of bacterial genomes 21 . However, in the case of the genus Brucella, very little information is available on the diversity of ISs within the genomes, and only IS711 has been extensively studied in this respect until now. With this background, we designed our study to explore novel genetic polymorphisms among Brucella species that can be exploited to establish the differentiation among them for diagnostic applications.
The distinctively remarkable pattern of insertions and deletions of insertional sequences with precise lengths with the determined target sequences in our study, which was exploited in the development of a novel multiplex PCR assay for use in the diagnosis of brucellosis, encouraged us to understand the mechanism involved in the presence or absence of such precise lengths of IS elements. Our opinion was that one mechanism was possibly involved in the unique pattern of inserting DNA blocks with identical lengths in four different species of Brucella, but not one species. To our surprise, the diversity was not limited to the range of complexity observed among each fragment but to the possible mechanisms involved in the transposition of the DNA blocks.
The typical segmental genome duplication observed in the case of the insertion of the Omp31b gene fragment within the conserved Omp31 gene is a feature usually observed among the IS21 family. Segmental duplications usually occur at more than one site within a given genome with a high level of sequence identity, as confirmed by both in silico and in situ hybridization experiments, and such patterns are usually observed in the human genome 23,24 . Such segmental duplications are observed as an indicative sign of evolutionary rearrangement 25,26 .
In the case of the B. abortus and B. suis rescue fragments in our study, no features of insertional sequences were observed. A lack of such features is an indication of possible homologous intermolecular or intramolecular recombination that would have occurred between two insertional elements each with different direct repeats. Such recombination usually results in the formation of a hybrid element carrying one direct repeat of each parent. In some cases, formation of adjacent deletions from duplicative intramolecular transpositions can also result in such insertional sequence elements. A single copy of the direct repeat could be located on each of the deletion products in such cases. However, extensive studies are required to confirm the exact mechanism involved in such insertion events on a case-by-case basis throughout the genome.
The generation of short directly repeated sequences flanking the inserted DNA sequences is a general feature of insertional sequence elements. Such insertions are facilitated by the attack of each DNA strand at the target direct repeat site by one of the two transposon ends in a staggered way. Usually, the lengths of such direct repeats range between 2 and 14 bp. In the case of the B. canis rescue fragment, the flanking direct repeat was 2 bp in length. Usually, in the case of such direct repeats, the insertion sequence will be generated as a genome duplication of a fixed length, but nevertheless, generation of small variations in the transposition complex is also observed at a low frequency. Latter condition was observed in the rescue fragment of B. canis in our study.
The intention of the present work was to develop a detection system for simultaneous and differential identification of five Brucella species irrespective of biovars and strains based on comparative genomic analysis using the updated database of Brucella genomes. The specificity of the intended detection system is crucial for the diagnosis of brucellosis, which solely depends on the target genes. With the comparative whole genomic analysis approach for mining species-specific nucleotide sequences, a remarkable pattern of insertion or deletion of certain length nucleotides was encountered in the genomes of B. abortus, B. melitensis, B. suis, B. ovis and B. canis. Interestingly, the deciphered sequences were not limited to genes required for virulence or other specialized metabolic functions but also included noncoding sequences.
The distinctively remarkable pattern of insertions and deletions of nucleotide sequences or DNA blocks with particular lengths with the determined target sequences separated them into two groups, one group with deletions and the other four groups with insertions, supporting the development of the exemplary detection tool. The unique sequences were differentially present in four species but absent in one and were intergenic and/or intragenic in nature, inserting within conserved regions of the Brucella genus as detection targets. Species-specific deletions in the case of Brucella were reported earlier by Rajashekara and coworkers (2004) based on microarray analysis, and these deletions were found to be highly conserved. Successful exploitation of those deletions in the design of species-specific primers for the differentiation of B. neotomae and B. canis was described in the previously mentioned study. The primers in that study not only resulted in different amplicon sizes in the cases of B. neotomae and B. canis but also showed one or completely different PCR amplification bands in the case of the other Brucella species screened. Similar differential patterns were derived for specific and differential identification of other Brucella species by the exploitation of different deletions. However, as a limitation in that study, biovars within a species showed different banding patterns, making interpretation a difficult and challenging task 27 .
In the current study, irrespective of biovars, in the unique sequence of B. abortus, which we selected randomly, to include as a detection target, a deletion www.nature.com/scientificreports www.nature.com/scientificreports/ the rescue amplicons of B. canis, its specific amplicon had a 297 bp deletion from the 129 th nucleotide position, making it a 224 bp product.
Five pairs of oligonucleotide primers were designed to exploit the newly identified species-specific nucleotide sequences for the development of a multiplex PCR assay that identifies five Brucella species in one single reaction. Initially, each species-specific primer was exploited in the development of a monoplex PCR system for each species with specific amplicons of 1154 bp, 745 bp, 446 bp, 290 bp and 224 bp for B. abortus, B. melitensis, B. ovis, B. suis and B. canis, respectively. However, when the same five species-specific primer pairs were employed in a multiplex PCR for the detection of each species, a unique multiplex PCR system resulted. This was because all the primers would anneal to the conserved regions flanking all five nucleotide sequences of the Brucella genome, resulting in four additional amplicons termed "rescue amplicons" along with the species-specific amplicon. Nevertheless, in a multiplex PCR in the presence of all five Brucella species, an integrative format comprising all the species-specific amplicons was obtained for the identification and differentiation of all five species together in a single lane by gel electrophoresis. Furthermore, these primers were employed in the stabilization of conventional monoplex and multiplex PCR systems. To avoid false negative results, competitive internal amplification control (IAC) was also incorporated in the multiplex PCR system with pUC19 plasmid DNA.
The stabilized monoplex and multiplex PCRs were evaluated for specificity with reference strains of the said five Brucella species as well as non-Brucella species. As an outcome, irrespective of biovars, a highly specific PCR-based detection format was developed that could successfully and differentially identify B. abortus, B. melitensis, B. ovis, B. suis and B. canis. The major advantage of our assay is that the use of any one primer pair was also sufficient to identify Brucella at the genus level and species level by seeking specific and/or rescue amplicons irrespective of biovars. Each pair of primers was useful not only as a species-specific marker but also as a genus-specific marker.
The developed multiplex PCR was used to evaluate its applicability in detecting 80 strains of Brucella species. The assay successfully identified the five target species and differentiated them from each other. This result reflected the unambiguous specificity of the assay and thus its applicability in the diagnosis of brucellosis. However, to examine the robustness of this tool in clinical cases, multiplex PCR assays require extensive validation in several reference laboratories with a large number of clinical samples.