Skip to main content

Brucella spp. distribution, hosting ruminants from Greece, applying various molecular identification techniques



Brucellosis still remains an endemic disease for both livestock and human in Greece, influencing the primary sector and national economy in general. Although farm animals and particularly ruminants constitute the natural hosts of the disease, transmission to humans is not uncommon, thus representing a serious occupational disease as well. Under this prism, knowledge concerning Brucella species distribution in ruminants is considered a high priority. There are various molecular methodologies for Brucella detection with however differential discriminant capacity. Hence, the aim of this survey was to achieve nationally Brucella epidemiology baseline genotyping data at species and subtype level, as well as to evaluate the pros and cons of different molecular techniques utilized for detection of Brucella species. Thirty-nine tissue samples from 30 domestic ruminants, which were found positive applying a screening PCR, were tested by four different molecular techniques i.e. sequencing of the 16S rRNA, the BP26 and the OMP31 regions, and the MLVA typing panel 1 assay of minisatellite markers.


Only one haplotype was revealed from the 16S rRNA sequencing analysis, indicating that molecular identification of Brucella bacteria based on this marker might be feasible solely up to genus level. BP26 sequencing analysis and MLVA were in complete agreement detecting both B. melitensis and B. abortus. An interesting exception was observed in 11 samples, of lower quality extracted DNA, in which not all expected MLVA amplicons were produced and identification was based on the remaining ones as well as on BP26. On the contrary OMP31 failed to provide a clear band in any of the examined samples.


The present study reveals the constant circulation of Brucella bacteria in ruminants throughout Greece. Further, according to our results, BP26 gene represents a very good alternative to MLVA minisatellite assay, particularly in lower quality DNA samples.

Peer Review reports


Brucellosis is globally one of the most severe and debilitating zoonoses affecting livestock and humans [1, 2]. Ruminants are also considered the major natural hosts, maintaining the causative agent of the disease in the environment [3,4,5]. The disease is transmitted to humans either by direct contact with infected animals or by consumption of unpasteurized dairy products [6]. In this context, high prevalence of the disease is accompanied by economic collapse for the stakeholders owing to the produced milk reduction, abortions and forced slaughter. Brucellosis is endemic in countries around the Mediterranean basin, western Asia and parts of Africa as well as Latin America [7] while it constitutes an important occupational disease for breeders, employees in slaughterhouses and veterinarians as well.

The etiological agents of brucellosis are the proteobacteria belonging in the genus Brucella (family: Brucellaceae). Twelve species have been described, with B. abortus mostly infecting cattle while B. melitensis mostly infecting small ruminants. It should be also mentioned that Brucella melitensis is responsible for the vast majority of human cases worldwide [8]. In a relevant study in Greece, more than 90% of Brucella spp. strains from clinical specimens was identified as B. melitensis [9]. It should be noted that in Greece, an ongoing vaccination and monitoring program takes place during the last two decades, in an effort to eradicate the disease by slaughtering infected animals and protecting young ones [10,11,12,13].

For the effective monitoring of both brucellosis control programs and human disease, it is important to have reliable tests to differentiate vaccine and field strains. Many molecular approaches have been developed to detect vaccine strains [14,15,16,17].

More specifically, many PCR based attempts have been conducted to develop an efficient protocol, able to identify Brucella genus and distinguish species. At genus specific level, Brucella detection molecular tools have been developed targeting various conserved genomic regions, such as the 16S–23S intergenic transcribed [18, 19], the BCSP31 gene [20, 21], the 16S rRNA [22], the perosamine synthetase (per) gene [23], the gene encoding the Omp2a protein antigen [24], the outer membrane proteins (omp2b, omp2a and omp31) [25] and the proteins of the omp25/omp31 family of Brucella spp. [26]. Nevertheless, the above assays vary greatly in sensitivity and specificity [27].

PCR based assays targeting species specific level identification are also numerous [28,29,30,31,32,33], whereas the multi-locus variable-number tandem-repeat analysis (MLVA) assay developed by Le Flèche et al. [30] is considered the gold standard molecular methodology for Brucella typing [34]. An appropriate PCR based detection protocol should be able to discriminate vaccine Brucella strains (RB51 and Rev1) from pathogenic ones, with many available assays towards this direction [16, 29, 35,36,37].

Altogether, the main scope of the present study was the investigation of Brucella species in ruminants from Greece, at both species and subtype level. Particularly, our objectives were to identify Brucella species and subtypes of all the strains, which are very important for epidemiologic surveillance and investigation of outbreaks in brucellosis endemic regions such as Greece [38, 39]. Eventually, in an effort to evaluate the advantages and drawbacks of different molecular techniques utilized for detection and discrimination of Brucella species, an additional goal of the study was to evaluate the efficacy and reliability of different Brucella identification methodologies.


16S rRNA sequencing and phylogenetic analysis

Initially all 39 examined samples were subjected to amplification, sequencing and neighbor joining phylogenetic analysis of 16S rRNA. Sequencing results revealed only one identical haplotype in all derived sequences, confirming that all samples belonged to the genus Brucella with 100% similarity with other species of Brucella (Fig. 1). This haplotype was deposited in the GenBank database and given the accession number OM570553. Those inferences demonstrate that 16S rRNA based molecular identification of Brucella bacteria may be feasible solely up to genus level.

Fig. 1
figure 1

Phylogenetic dendrogram of the examined Brucella samples in comparison with haplotypes retrieved from GenBank, based on the 16S rRNA (a) and on the BP26 gene (b). Haplotypes found in the present study are indicated with a square

MLVA assay

The panel 1 of MLVA minisatellite markers was applied afterwards in the 39 samples with each PCR being performed at least twice. Firstly, three pairs of primers (Bruce06, Bruce11, Bruce42) were used since they provide an easy to observe differentiation method between the species. Bruce06 amplifies an expected amplicon of 408 bp in case of Brucella melitensis and 542 bp in Brucella abortus, Bruce11 a product of 257 bp and 383 bp, and Bruce42 539 bp and 289 bp respectively. Afterwards, the pairs of Bruce08 and Bruce12 were used, followed by the remaining three markers i.e. Bruce42, Bruce 43 and Bruce55, which were applied only for samples in which a positive result was observed with at least one of the first five markers, as Bruce42, Bruce 43 and Bruce55 result in an identical amplicon for Brucella melitensis and Brucella abortus.

The majority of the samples (28 out of the 39) was successfully identified with all the 8 primer pairs of the MLVA assay Panel 1 (Supplementary Material 1). The remaining 11 samples gave the expected amplicon with two or three primer pairs (from the first five markers) and failed to give an amplicon or gave non-specific bands with the rest of the first five primer pairs (Table 1). Interestingly, these samples were of lower DNA quality in terms of purity as measured by the 260/280 absorbance ratio (Table 1).

Table 1 Geographic, animal and tissue origin of the examined samples, and Brucella identification based on various molecular techniques

BP26 and OMP31 analysis

Identification based on the BP26 marker was in complete agreement with MLVA results, as only two haplotypes were defined, one assigned to B. melitensis and one to B. abortus, clearly differentiated from each other and identical with conspecific sequences obtained from GenBank (Fig. 1). These two haplotypes were deposited in the GenBank database and given the accession numbers OM628689 and OM628690, for B. melitensis and for B. abortus respectively. Nevertheless, this was not always the case regarding the identification with the screening PCR, where there were two samples in which the two methods did not agree (Table 1). Finally, OMP31 marker failed to provide a clear band or sequence in any of the examined samples.

Brucella infection summarized results

In total, concerning billy goats all of the examined samples were identified as Brucella melitensis, whereas in rams and bulls, only B. melitensis and B. abortus were identified, respectively. Furthermore, in 3 billy goats and in 1 ram which were PCR positive for vaccine strain Rev1, it was confirmed by MLVA. Regarding females, in 2/3 goats B. melitensis was detected, with the remaining on being B. abortus, whereas 5/6 ewes hosted B. melitensis with again the remaining one being B. abortus, and all 17 examined cows hosting B. abortus.


In Greece, brucellosis still remains an endemic zoonosis infecting a wide range of animal species with an influence both in public health and in national economy. Although in 1977, a national program against ruminants’ brucellosis was designed and initiated, the disease is not eradicated yet. The generally accepted pathogenic agent of bovine brucellosis is B. abortus, wherea B. melitensis is only occasionally detected. Respectively, the common pathogenic agent of caprine and ovine brucellosis is B. melitensis, whereas B. abortus is only rarely found in small ruminants. Thus, small ruminants are considered as the main hosts for B. melitensis [5, 40]. Nevertheless, even though Giantzis et al. [41, 42] detected B. melitensis in bovine aborted fetuses, they found no B. abortus in aborted fetuses from ewes and goats during a five-year period of 1976-1981. The clinical, pathological and epidemiological picture of caprine brucellosis due to B. melitensis is similar to B. abortus infection in cattle. The dominant strain for human brucellosis is B. melitensis [43,44,45,46]. On the other hand, Giannakopoulos et al. [47] referred human cases in Western Greece where B. abortus was identified as an equally frequent pathogenic agent. On top of that, it is of high importance to investigate whether interspecies transmission of B. melitensis, B. abortus and the vaccine strain Rev1 may occur naturally and cause clinical disease in domestic ruminants or may perplex the standard laboratory exams (RBT & CFT). It should be emphasized that control and eradication policies may have to be sometimes readapted to the new data.

In our study, the presence of B. abortus was revealed in two such cases, detected from one ewe and one goat. This could be explained by several reasons such as the existence of mixed livestock pastures, the coexistence of sheep and goats in the same shelter with cattle, or occasionally owing to herd movements with infected but non-detected animals. The results suggest cross-species infection of B. abortus from cattle to small ruminants raised in close contact [48].

The fact that in four non-vaccinated male small ruminants, the vaccine strain Rev1 was detected (Table 1), possibly indicates that during the vaccine administration in females, errors may have occurred. As a result of the accidental vaccination, the male animals may be characterized as false seropositive during the control program. Finally, due to the aforementioned vaccine administration errors, male animals end up to the slaughterhouse while the farm comes under quarantine for a period of at least 2 to 6 weeks. For these reasons, it is important to improve molecular methods to detect the exact pathogenic agent from samples from alive seropositive animals or in case of slaughter to detect it in a rapid and effective way.

In particular, development and improvement of efficient DNA-based methods is a consistent demand, as well their comparative efficiency. In this context, here we utilized and evaluated four different molecular tools towards this scope. Particularly, apart from the “gold standard” MLVA assay, the 16S rRNA, the OMP31 and the BP26 genes were analyzed. Assays targeting Brucella non-coding genomic tandem repeats, such as the MLVA minisatellites [30] or the microsatellite loci proposed earlier [15, 49], have been developed based on the principle of the high genetic homogeneity levels of the different strains, which fail to be distinguished by classic conserved genomic region markers such as the 16S rRNA. Nevertheless, microsatellites, although are excellent markers for genotyping and genetic differentiation studies, on account of hypervariability and high mutation rates may be inadequate to assign the biovar correctly or an isolate at species level. Minisatellites (MLVA) [30] on the contrary, possess a better species identification capability and are moderately variable higher, resulting in higher discriminatory power. On the other hand, gene sequences possess lower mutation rates. In line with Gupta et al. [50] who concluded that for genus-based identification of Brucella species, 16S rRNA and 16S-23 rRNA gene are the best target, in our case examination of 16S rRNA produced one identical haplotype in all samples, thus not capable to discriminate the different Brucella species. Our results are in agreement with the same study [50] concerning identification based on BP26 gene targets.

Interestingly, BP26 worked better in lower quality samples, where not all minisatellite markers gave product (Table 1). In PCR based methods the quality and purity of Brucella spp. DNA is a crucial prerequisite before performing these methods, especially for multiplex PCR methods. Any inhibitor in DNA samples from any source can affect the result of a PCR based method. Sensitivity and specificity of most PCR-based methods are not well established and their real capability for use with clinical samples and hence diagnosis has not been validated [27]. Specifically, for investigation of Brucella presence and species identification, samples may originate from slaughterhouses or dead animals, sometimes travelling for days and eventually received at the labs highly degraded. Sequencing of the BP26 gene represents a very good alternative to MLVA minisatellite assay that according to our results constitutes an effective alternative marker in lower quality DNA samples.


The present research provides clear evidence for the continuing circulation of Brucella species in ruminants throughout Greece, a country that still remains therefore endemic in Brucellosis. All but one molecular techniques tested were proved informative and efficient. However our results clearly demonstrate a better performance of BP26 gene marker in lower quality samples, in terms of optical density ratio.


Sample collections

During a period of 3 years (2016-2018), 264 samples were collected from a total of 191 farmed ruminants, originating from farming units located in the mainland of Greece throughout the country. Those samples derived either from positive males to Rose Bengal Test (RBT) and/or Component Fixation Test (CFT) either from aborted fetuses of sheep, goats, and cattle. Thirty-nine (39) out of the 264 tissue samples, which derived from 30 animals, were found PCR positive using the screening Multiplex PCR detection molecular methodology of Garcia-Yoldi et al. [29] for the needs of the Greek Ministry of Rural Development [51] and were utilized for the needs of the present study. The 39 samples originated from various organs and tissues from 30 male and female domestic ruminants, i.e. 7 sheep, 7 goats and 16 bovines, including 7 testicles, 1 spleen, 5 lymph nodes, 11 embryonic rumens, 11 embryonic livers, 4 cotyledonary placentas. The above 39 samples were found positive by the Multiplex PCR Assay for Brucella spp. [29] in our laboratory [51]. All samples are displayed in detail in Table 1.

Extraction of genomic DNA from tissue sample

All procedures were performed under the Biosafety Level three (BSL3) guidelines [52]. Tissue samples were initially homogenized within 200μl sterile Phosphate Buffered Saline (PBS) Sigma- Aldrich using a tissue grinder. For homogenization, aseptic processing of all samples was performed by removal of extraneous material and further maceration and chopping into small pieces in PBS. DNA isolation was carried out using the High Pure PCR template preparation kit (Roche, Basel, Switzerland) following the manufacturer instruction with a final elution volume of 80 μl in each sample. The concentration and the quality of the extracted DNA were evaluated in a micro-volume Q5000 UV-Vis Spectrophotometer (Quawell, USA).

Molecular identification of Brucella species

Since our aim was also to evaluate the applicability of different molecular tools for Brucella detection, four different molecular techniques were applied for identification of bacterial species and assignment of Brucella positive samples to species level.

Initially the 16S rRNA gene was amplified using the universal for bacteria species primer pair 27F-1492R (Table 2) that may identify the greatest majority of bacteria species at least at genus level. This pair of primers amplifies nearly the complete 16S rRNA. The 27F and 1492R [53] primers corresponded to positions 8–27 and to positions 1492–1513 of Escherichia coli 16S rRNA, respectively [56].

Table 2 Primers used for Brucella identification

Positive Brucella samples were then subjected to PCRs targeting the BP26 and the OMP31 genes, using the primer pairs BP26F-BP26R [54] and OMP31F-OMP31R [55] (Table 2). BP26 is a conserved gene capable of distinguishing the different Brucella species [57]. On the other hand, OMP31 gene is only present in B. melitensis. Cloning and sequencing of B. melitensis 16M OMP31 [55], a gene coding for a major Brucella outer membrane protein, verified that this gene is missing from B. abortus strains [58].

Finally, the MLVA typing panel 1 assay of minisatellite markers was applied [30, 59]. For MLVA analysis we worked with the set of primers (Bruce06, 08, 11, 12, 42, 43, 45, 55, Table 2), which amplifies 8 minisatellite markers with a good species identification capability [30]. All PCRs were performed in 20 μl final volumes, containing 0.6 pmol of each forward or reverse primer, 10 μl KAPA 2G Fast Hot Start readymix (Merck, Germany), approximately 50 ng extracted DNA and nuclease free water up to the final volume. Reactions were performed in a FastGene ULTRA Cycler (Nippon Genetics, Japan) and amplification program was as follows: after an initial denaturation at 95C for 3 min, 35 cycles were performed, of denaturation at 95C for 30 sec, annealing at 48-60C (Table 2) for 40 sec, and extension at 72C for 40-60 sec, depending on the product length (40 sec for products smaller than 1000 bp; 60 sec for products larger than 1000 bp), and final extension at 72C for 10 min. The amplified products were examined by electrophoresis in a 1.5% agarose gel stained with ethidium bromide (0.5 mg/ml) and photographed by photo documentation system, or in polyacrylamide gel electrophoresis stained with silver nitrate in cases of very small size products differences separation. To ensure reproducibility, each PCR was performed at least twice. Particularly for 16S rRNA, BP26 and OMP31 regions, successfully amplified products were purified using the NucleoSpin Gel and PCR Clean-up kit (Macherey-Nagel, Germany) following the manufacturer’s recommended protocol and sequenced bidirectionally in an ABI 3730xl automatic sequencer. Sequences produced, were aligned using the software MEGA 7 [60] and compared with conspecific and congeneric ones obtained from the GenBank database in Neighbor Joining phylogenetic trees that were created in the same software applying a bootstrap value of 1000 iterations.

Availability of data and materials

The datasets generated and analysed during the current study are available in the NCBI GenBank repository, accession numbers OM628689, OM628690 and OM570553” (,,



variable number tandem repeat analysis


Rose Bengal Test


Component Fixation Test


  1. Gwida M, Al Dahouk S, Melzer F, Rösler U, Neubauer H, Tomaso H. Brucellosis–regionally emerging zoonotic disease? Croat Med J. 2010;51(4):289–95.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Dean AS, Crump L, Greter H, Schelling E, Zinsstag J. Global burden of human brucellosis: a systematic review of disease frequency. PLoS Negl Trop Dis. 2012;6(10):e1865.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Blasco JM, Molina-Flores B. Control and eradication of Brucella melitensis infection in sheep and goats. Vet Clin Food Anim Pract. 2011;27(1):95–104.

    Article  Google Scholar 

  4. Elzer PH, Hagius SD, Davis DS, DelVecchio VG, Enright FM. Characterization of the caprine model for ruminant brucellosis. Vet Microbiol. 2002;90(1-4):425–31.

    Article  CAS  PubMed  Google Scholar 

  5. Díaz A. Epidemiology of brucellosis in domestic animals caused by Brucella melitensis, Brucella suis and Brucella abortus. Revue scientifique et technique-Office international des epizooties. 2013;32(1).

  6. Atluri VL, Xavier MN, De Jong MF, Den Hartigh AB, Tsolis RM. Interactions of the human pathogenic Brucella species with their hosts. Annu Rev Microbiol. 2011;13(65):523–41.

    Article  CAS  Google Scholar 

  7. Moreno E. Retrospective and prospective perspectives on zoonotic brucellosis. Front Microbiol. 2014;13(5):213.

    Google Scholar 

  8. Corbel MJ. Food and Agriculture Organization of the United Nations, World Health Organization & World Organisation for Animal Health. Brucellosis in humans and animals. World Health Organization; 2006.

  9. Dougas G, Katsiolis A, Linou M, Kostoulas P, Billinis C. Modelling Human Brucellosis Based on Infection Rate and Vaccination Coverage of Sheep and Goats. Pathogens. 2022;11(2):167.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Tzani M, Katsiolis A, Program for the control and eradication of brucellosis in sheep and goats in Greece. HCDCP e-bulletin. 2012;2(4):19-20 ISSN 1792-9016. accessed on 13 March 2021.

  11. European Commission 2017. SANTE/10769/2017 Report of the brucellosis Task Force sub-group. Meeting held in Athens Greece 29-31 March 2017. accessed online on 10 Jan 2022.

  12. Katsiolis A, Thanou O, Tzani M, Dile C, Korou M, Stournara A, et al. Investigation of the human resources needs for the effective and efficient implementation of the sheep and goat brucellosis program in Greece. Vet J Repub Srpska (Banja Luka). 2018;18(2):270–96.

    Google Scholar 

  13. Strausbaugh LJ, Berkelman RL. Human illness associated with use of veterinary vaccines. Clin Infect Dis. 2003;37(3):407–14.

    Article  Google Scholar 

  14. Cloeckaert A, Grayon M, Grépinet O. Identification of Brucella melitensis vaccine strain Rev. 1 by PCR-RFLP based on a mutation in the rpsL gene. Vaccine. 2002;20(19-20):2546–50.

    Article  CAS  PubMed  Google Scholar 

  15. Bricker BJ, Ewalt DR, Halling SM. Brucella'HOOF-Prints': strain typing by multi-locus analysis of variable number tandem repeats (VNTRs). BMC Microbiol. 2003;3(1):1–3.

    Article  Google Scholar 

  16. Lopez-Goñi I, Garcia-Yoldi D, Marin CM, De Miguel MJ, Munoz PM, Blasco JM, et al. Evaluation of a multiplex PCR assay (Bruce-ladder) for molecular typing of all Brucella species, including the vaccine strains. J Clin Microbiol. 2008;46(10):3484–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Gopaul KK, Sells J, Bricker BJ, Crasta OR, Whatmore AM. Rapid and reliable single nucleotide polymorphism-based differentiation of Brucella live vaccine strains from field strains. J Clin Microbiol. 2010;48(4):1461–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Rijpens NP, Jannes G, Van Asbroeck MA, Rossau R, Herman LM. Direct detection of Brucella spp. in raw milk by PCR and reverse hybridization with 16S-23S rRNA spacer probes. Appl Environ Microbiol. 1996;62(5):1683–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Bricker BJ, Ewalt DR, MacMillan AP, Foster G, Brew S. Molecular characterization of Brucella strains isolated from marine mammals. J Clin Microbiol. 2000;38(3):1258–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Baily GG, Krahn JB, Drasar BS, Stoker NG. Detection of Brucella melitensis and Brucella abortus by DNA amplification. J Trop Med Hyg. 1992;95(4):271–5.

    CAS  PubMed  Google Scholar 

  21. Luna L, Mejía G, Barragán V, Trueba G. Molecular Detection of Brucella Species in Ecuador. Intern J Appl Res Vet Med. 2016;14(2):185–9.

    CAS  Google Scholar 

  22. Romero C, Gamazo C, Pardo M, Lopez-Goñi I. Specific detection of Brucella DNA by PCR. J Clin Microbiol. 1995;33(3):615–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bogdanovich T, Skurnik M, Lubeck PS, Ahrens P, Hoorfar J. Validated 5′ nuclease PCR assay for rapid identification of the genus Brucella. J Clin Microbiol. 2004;42(5):2261–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Leal-Klevezas DS, Martínez-Vázquez IO, Lopez-Merino A, Martínez-Soriano JP. Single-step PCR for detection of Brucella spp. from blood and milk of infected animals. J Clin Microbiol. 1995;33(12):3087–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Imaoka K, Kimura M, Suzuki M, Kamiyama T, Yamada A. Simultaneous detection of the genus Brucella by combinatorial PCR. Jpn J Infect Dis. 2007;60(2/3):137.

    CAS  PubMed  Google Scholar 

  26. Vizcaíno N, Caro-Hernández P, Cloeckaert A, Fernández-Lago L. DNA polymorphism in the omp25/omp31 family of Brucella spp.: identification of a 1.7-kb inversion in Brucella cetaceae and of a 15.1-kb genomic island, absent from Brucella ovis, related to the synthesis of smooth lipopolysaccharide. Microbes Infect. 2004;6(9):821.

    Article  PubMed  CAS  Google Scholar 

  27. Yu WL, Nielsen K. Review of detection of Brucella spp. by polymerase chain reaction. Croat Med J. 2010;51(4):306–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Hinić V, Brodard I, Thomann A, Cvetnić Ž, Makaya PV, Frey J, et al. Novel identification and differentiation of Brucella melitensis, B. abortus, B. suis, B. ovis, B. canis, and B. neotomae suitable for both conventional and real-time PCR systems. J Microbiol Methods. 2008;75(2):375–8.

    Article  PubMed  CAS  Google Scholar 

  29. Garcia-Yoldi D, Marín CM, de Miguel MJ, Munoz PM, Vizmanos JL, López-Goñi I. Multiplex PCR assay for the identification and differentiation of all Brucella species and the vaccine strains Brucella abortus S19 and RB51 and Brucella melitensis Rev1. Clin Chem. 2006;52(4):779–81.

    Article  CAS  PubMed  Google Scholar 

  30. Le Flèche P, Jacques I, Grayon M, Al Dahouk S, Bouchon P, Denoeud F, et al. Evaluation and selection of tandem repeat loci for a Brucella MLVA typing assay. BMC Microbiol. 2006;6(1):1–4.

    Article  CAS  Google Scholar 

  31. Bricker BJ, Halling SM. Differentiation of Brucella abortus bv. 1, 2, and 4, Brucella melitensis, Brucella ovis, and Brucella suis bv. 1 by PCR. J Clin Microbiol. 1994;32(11):2660–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Mayer-Scholl A, Draeger A, Göllner C, Scholz HC, Nöckler K. Advancement of a multiplex PCR for the differentiation of all currently described Brucella species. J Microbiol Methods. 2010;80(1):112–4.

    Article  CAS  PubMed  Google Scholar 

  33. Bricker BJ, Ewalt DR, Olsen SC, Jensen AE. Evaluation of the Brucella abortus species–specific polymerase chain reaction assay, an improved version of the Brucella AMOS polymerase chain reaction assay for cattle. J Vet Diagn Invest. 2003;15(4):374–8.

    Article  PubMed  Google Scholar 

  34. Pelerito A, Nunes A, Grilo T, Isidro J, Silva C, Ferreira AC, et al. Genetic Characterization of Brucella spp. Front Microbiol. 2021 Nov;12:12.

    Article  Google Scholar 

  35. Sangari FJ, Agüero J. Identification of Brucella abortus B19 vaccine strain by the detection of DNA polymorphism at the ery locus. Vaccine. 1994;12(5):435–8.

    Article  CAS  PubMed  Google Scholar 

  36. Vemulapalli R, McQuiston JR, Schurig GG, Sriranganathan N, Halling SM, Boyle SM. Identification of an IS 711 element interrupting the wboA gene of Brucella abortus vaccine strain RB51 and a PCR assay to distinguish strain RB51 from other Brucella species and strains. Clin Diagn Lab Immunol. 1999;6(5):760–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Ewalt DR, Bricker BJ. Validation of the abbreviated Brucella AMOS PCR as a rapid screening method for differentiation of Brucella abortus field strain isolates and the vaccine strains, 19 and RB51. J Clin Microbiol. 2000;38(8):3085–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Al Dahouk S, Le Flèche P, Nöckler K, Jacques I, Grayon M, Scholz HC, et al. Evaluation of Brucella MLVA typing for human brucellosis. J Microbiol Methods. 2007;69(1):137–45.

    Article  PubMed  CAS  Google Scholar 

  39. Marianelli C, Graziani C, Santangelo C, Xibilia MT, Imbriani A, Amato R, et al. Molecular epidemiological and antibiotic susceptibility characterization of Brucella isolates from humans in Sicily, Italy. J Clin Microbiol. 2007;45(9):2923–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Karvounaris PA. Aspects cliniques et épizootiologiques de la brucellose bovine, caprine et ovine, en Grèce [Clinical and epizootiological aspects of bovine, caprine, and ovine brucellosis in Greece]. Dev Biol Stand. 1976;31:254–64 French. PMID: 1261740.

    CAS  PubMed  Google Scholar 

  41. Giantzis DG. Brucella eradications program: Course and considerations serological and microbiological tests results. J Hellenic Vet Med Soc. 1984;35(1):19–25.

    Google Scholar 

  42. Giantzis DG, Xenos G, Pashaleri E. Investigation of abortion caused by infective agents in sheep and goats (in Greece). Ellenike Kteniatrike Hellenic Vet Med. 1984;27:132–41 (in Greek).

    Google Scholar 

  43. Galanakis E, Kefallinou A, Kostoula-Tsiara A, Lapatsanis PD. Childhood brucellosis in Epirus (NW Greece) during the decade 1980-89. Paediatriki. 1993;56:162–8.

    Google Scholar 

  44. Galanakis E, Bourantas KL, Leveidiotou S, Lapatsanis PD. Childhood brucellosis in north-western Greece: a retrospective analysis. Eur J Pediatr. 1996;155(1):1–6.

    Article  CAS  PubMed  Google Scholar 

  45. Tsolia M, Drakonaki S, Messaritaki A, Farmakakis T, Kostaki M, Tsapra H, et al. Clinical features, complications and treatment outcome of childhood brucellosis in central Greece. J Infect. 2002;44(4):257–62.

    Article  CAS  PubMed  Google Scholar 

  46. Fouskis I, Sandalakis V, Christidou A, Tsatsaris A, Tzanakis N, Tselentis Y, et al. The epidemiology of Brucellosis in Greece, 2007–2012: a ‘One Health’approach. Trans R Soc Trop Med Hyg. 2018;112(3):124–35.

    Article  PubMed  Google Scholar 

  47. Giannakopoulos I, Nikolakopoulou NM, Eliopoulou M, Ellina A, Kolonitsiou F, Papanastasiou DA. Presentation of childhood brucellosis in Western Greece. Jap J Infect Dis. 2006;59(3):160.

    Google Scholar 

  48. Wareth G, Melzer F, Tomaso H, Roesler U, Neubauer H. Detection of Brucella abortus DNA in aborted goats and sheep in Egypt by real-time PCR. BMC Res Notes. 2015;8(1):1–5.

    Article  Google Scholar 

  49. Bricker BJ, Ewalt DR. Evaluation of the HOOF-Print assay for typing Brucella abortus strains isolated from cattle in the United States: results with four performance criteria. BMC Microbiol. 2005;5(1):1–0.

    Article  CAS  Google Scholar 

  50. Gupta VK, Shivasharanappa N, Kumar V, Kumar A. Diagnostic evaluation of serological assays and different gene based PCR for detection of Brucella melitensis in goat. Small Ruminant Res. 2014;117(1):94–102.

    Article  Google Scholar 

  51. Katsiolis A, Papanikolaou E, Stournara A, Giakkoupi P, Papadogiannakis E, Zdragas A, Giadinis ND, Petridou E. Molecular detection of Brucella spp. in ruminant herds in Greece. Trop Anim Health Prod. 2022;54:173.

  52. Zaki AN. Biosafety and biosecurity measures: management of biosafety level 3 facilities. Int J Antimicrob Agents. 2010;36:S70–4.

    Article  CAS  PubMed  Google Scholar 

  53. Frank JA, Reich CI, Sharma S, Weisbaum JS, Wilson BA, Olsen GJ. Critical evaluation of two primers commonly used for amplification of bacterial 16S rRNA genes. App Environment Microbiol. 2008;74(8):2461–70.

    Article  CAS  Google Scholar 

  54. Gupta VK, Kumari R, Vohra J, Singh SV, Vihan VS. Comparative evaluation of recombinant BP26 protein for serological diagnosis of Brucella melitensis infection in goats. Small Ruminant Res. 2010;93(2-3):119–25.

    Article  Google Scholar 

  55. Vizcaino N, Cloeckaert A, Zygmunt MS, Dubray G. Cloning, nucleotide sequence, and expression of the Brucella melitensis omp31 gene coding for an immunogenic major outer membrane protein. Infect Immun. 1996;64(9):3744–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Brosius J, Palmer ML, Kennedy PJ, Noller HF. Complete nucleotide sequence of a 16S ribosomal RNA gene from Escherichia coli. Proc Natl Acad Sci. 1978;75(10):4801–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Seco-Mediavilla P, Verger JM, Grayon M, Cloeckaert A, Marín CM, Zygmunt MS, et al. Epitope mapping of the Brucella melitensis BP26 immunogenic protein: usefulness for diagnosis of sheep brucellosis. Clin Vaccine Immunol. 2003;10(4):647–51.

    Article  CAS  Google Scholar 

  58. Vizcaino N, Verger JM, Grayon M, Zygmunt MS, Cloeckaert A. DNA polymorphism at the omp-31 locus of Brucella spp.: evidence for a large deletion in Brucella abortus, and other species-specific markers. Microbiology. 1997;143(9):2913–21.

    Article  PubMed  Google Scholar 

  59. Jiang H, Wang H, Xu L, Hu G, Ma J, Xiao P, et al. MLVA genotyping of Brucella melitensis and Brucella abortus isolates from different animal species and humans and identification of Brucella suis vaccine strain S2 from cattle in China. PLoS One. 2013;8(10):e76332.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


The authors are grateful to all the staff and employees that helped for the completion of the samplings


Funding through Research Committee of Aristotle University of Thessaloniki.

Author information

Authors and Affiliations



A.K. collected the samples, performed the molecular analyses and wrote the first draft of the manuscript. D.K.P. participated in the PCR analyses and wrote the first draft of the manuscript. I.G. corrected and wrote the final draft on the manuscript. K.P participated in sample collections and in molecular analyses. A.Z. conceived and corrected the paper. N.D.G coordinated the study and helped to draft the manuscript. E.P. conceived, coordinated the study and supervised the whole experiments. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ioannis A. Giantsis.

Ethics declarations

Ethics approval and consent to participate

All animal manipulations were carried out according to the EU Directive on the protection of animals’ usage for scientific purposes (2010/63/EU). The research protocol was approved by the General Assembly of the Veterinary Faculty of Aristotle University of Thessaloniki decision 55/27-05-2015. Also, permissions were obtained from all participating farm owners for collecting the samples.

Consent for publication

Not applicable

Competing interests

The authors declare no competing interests

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Katsiolis, A., Papadopoulos, D.K., Giantsis, I.A. et al. Brucella spp. distribution, hosting ruminants from Greece, applying various molecular identification techniques. BMC Vet Res 18, 202 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: