Journal Information
Vol. 25. Issue 1.
(January - February 2021)
Download PDF
More article options
Vol. 25. Issue 1.
(January - February 2021)
Original article
Open Access
Determination of investigation of the link between human and animal Brucella isolates in Iran using multiple-locus variable number tandem repeat method comprising 16 loci (MLVA-16)
Shiva Mirkalantari
Corresponding author

Corresponding author.
, Faramarz Masjedian
Corresponding author

Corresponding author.
, Amraee Fateme
Iran University of Medical Sciences, Faculty of Medicine, Microbiology Department, Tehran, Iran
This item has received

Under a Creative Commons license
Article information
Full Text
Download PDF
Figures (2)

Epidemiological studies are important tools to assess the diversity of Brucella isolates and to estimate their epidemiological relationship among isolates from different geographical origins. In this study the MLVA16 (multiple-locus variable number tandem repeat analysis based on 16 loci) was employed to investigate the diversity of Brucella spp. Isolated from humans and animals for epidemiological purposes and to determine the most common Brucella genotypes in Iran.


We designed a molecular-based study to evaluate the potential reservoirs of human brucellosis. After isolation and identification of 54 Brucella spp human and animal specimens from three regions of Iran, bacterial genomic DNA was extracted MLVA with three panel was used for the genotyping of isolates. The size of PCR products were analyzed and converted to repeat unit numbers using a published allele numbering system and data set was imported into Bionumerics.


Three isolates (5.55%) were identified as Brucella abortus and 51 (94.44%) as Brucella melitensis. Two isolates of Brucella abortus were from humans and one from an animal. Thirty-four Brucella melitensis isolates were from humans and 17 from animals. Using MLVA16-genotyping, 54 isolates with genetic similarity coefficient of 80% were divided into 46 genotypes and 22 genotypes were represented by a single isolate, while 4, 2, 1 and 2 genotypes were represented by 2, 3, 4 and 7 isolates, respectively. The most prevalent genotype was represented by 14 isolates. There were two other frequent genotypes each represented by seven isolates, among which only one was restricted to a geographic region. Discriminatory power for each locus was determined in this study and panel 2B shows the high discretionary power [Bruce04 (0.837), Bruce30 (0.806), Bruce 09 (0.787), Bruce 07 (0.772), Bruce16 (0.766)].


MLVA16 analysis of 54 Brucella isolates showed high level polymorphism in their genotypes. Only two genotypes, each observed in seven isolates, were related to one another and only one of these genotypes were found in to two separate regions.

Full Text

Brucellosis is a zoonotic bacterial disease caused by the genus Brucella, which is transmitted through direct or indirect contact to animals.1 Based on differences in pathogenicity, phenotypic characteristics and hosts, the genus Brucella is classified into 11 species.2,3 Four species, including Brucella abortus, Brucella melitensis, Brucella canis and Brucella suis are known to infect humans. Enormous economic losses and public health problems results from abortion, infertility in livestock, weak offspring, decreased milk production, and morbidity in humans. Although brucellosis has been eradicated in the USA, Canada, North Europe, and Australia, this infection is still highly prevalent in central Asia, Middle East, Mediterranean region, Africa, and Latin America.4 Iran is an endemic area for brucellosis and there is a major risk of Brucella transmission from eastern and western neighboring countries such as Iraq, Pakistan, and Afghanistan due to lack of high quality veterinary services for controlling animal diseases.5,6 Although, 500,000 human cases of brucellosis are annually reported worldwide, some cases remain undetected or neglected.7 Key approaches for preventing brucellosis in humans include controlling the disease in animals using epidemiological studies to assess diversity among strains, and estimating the epidemiological relationship between isolates from different geographical origins.8,9

Classical phenotyping methods, such as serotyping, phage typing, metabolic profiling and sensitivity to dyes have been employed for Brucella subtyping; however, these techniques are only available in reference laboratories, have a limited discriminatory power, require manipulation, and lack any standardized interpretation of this method cause difficulties.10

Based on these limitations, bacterial typing has shifted towards molecular identification by which epidemiological relationships are investigated among isolates. Thus far, several molecular typing methods, such as random amplified polymorphic DNA (RAPD)-PCR, amplified fragment-length polymorphism (AFLP), pulsed field gel electrophoresis (PFGE), polymerase chain reaction restriction fragment length polymorphism (RFLP), multiple-locus sequence typing (MLST) and multiple-locus VNTR (variable number tandem-repeats) analysis (MLVA) have been introduced. MLVA method is considered to be an effective and rapid tool for monitoring variations of copy numbers of tandem repeat units (TRs) with a high discriminatory power.11 This method is not only used in outbreaks and epidemiological and trace-back investigations, but also in confirmatory laboratories or screening of food-borne acquired infections. TR sequences are multiple alleles that can be presented on a single locus and be easily determined through agarose electrophoresis or capillary electrophoresis based on their size differences.12 MLVA has been proven to be a good technique for the assessment of pathogenic bacteria, such as Brucella, that display very little genomic diversity. MLVA schemes with 21, 15 and 16 loci (MLVA-21, MLVA-15 and MLVA-16, respectively) have been presented for discriminating Brucella spp. MLVA-16 has been proposed by Al-Dahank et al. with eight minisatellite markers (panel1: Bruce-06, Bruce-08, Bruce-11, Bruce-12, Bruce-42, Bruce-43, Bruce-45, and Bruce-55) for identification of bacteria spp and eight microsatellite markers (panel 2A: Bruce-18, Bruce-19, Bruce-21 and panel 2B: Bruce-04, Bruce-07, Bruce-09, Bruce-16 and Bruce-30) for screening the most common genotypes and classification of isolates to subspecies.13 The genetic diversity of Brucella isolated from human and animal specimens has not yet been investigated in Iran. The main objectives of this study was to employ the MLVA-16 assay to investigate and assess the diversity of Brucella isolates for epidemiological purposes and to determine the most common genotypes of Brucella isolates in Iran.

MethodsEthical statement

Ethical Committee of Iran University of medical science approved this study under the ethics code R.IUMS.REC 1396.33003.

Study design and sample size

This study is a molecular-based and cross-sectional study to evaluate the possible reservoirs of human brucellosis in Iran.

Isolation and determination of Brucella spp

The clinical samples were collected from humans and animals suffering from brucellosis. Culture isolation of Brucella spp was achieved by agar supplemented with 10% horse serum. The inoculated plates were incubated at 37 °C aerobically for 3–5 days. Single colonies were analyzed microscopically using Gram staining and biochemical (urease, oxidase and catalase) and motility tests were carried out for bacterial identification. Colonies confirmed as Brucella spp were subjected to PCR methods for species-level identification.

Preparation of Brucella DNA and amplification of loci

DNA of Brucella was extracted using the PCR template preparation Kit (Roche Diagnostics, Germany). MLVA genotyping was performed as previously described by Flech et al. and completed by Al Dahouk et al. This method includes three panels named panel1, panel 2A and panel 2B. The 16 primer pairs were divided into three groups. Panel1 consists of eight minisatellite loci (Bruce 06, 08, 11, 12, 42, 43, 45, 55), while panel 2A (Bruce18, 19, 21) and panel 2B (Bruce 04, 07, 09, 16 and 30) consisted of three and five minisatellite loci, respectively. PCR amplification was performed in a total volume of 15 µl containing 1 ng of DNA, 1X PCR reaction buffer, 1U of Taq DNA polymerase, 200 mM of each deoxynucleotide triphosphate, and 0.3 mM of each flanking primer. Amplification was performed in an Eppendrof termocycler as follow: an initial denaturation step at 96 °C for 5 min, followed by 30 cycles of denaturation at 96 °C for 30 s, annealing at 60 °C for 30 s, and elongation at 70 °C for 1 min. The final extension step was performed at 70 °C for 5 min. Five microliters of amplification product was loaded on 2.5% agarose gel and run under a voltage of 8 V/cm. Gel images were then recorded.

Analysis of MLVA data

The sizes of PCR products were analyzed and converted to repeat unit numbers using a published allele numbering system and the data set was imported to Bionumerics. For quantification of polymorphism at each locus, Hunted and Gaston diversity index (HGDI), available on the human protein HPA website ( was used. Categorical coefficient and unweighted pair group methods were applied for clustering analysis. Minimum spanning trees were also constructed.

Statistical analysis

Statistical analyses were conducted using the SPSS version 18, (Inc, Chicago, IL, USA).

ResultsCharacteristics of patients and Brucella isolates

A total of 54 Brucella spp. were respectively isolated from blood (62.96%) and cerebrospinal fluid (3.70%) of humans, as well as blood (16.66%), spleen (5.55%) and liver (11.11%) of animals. Human patients included 26 men (48.14%) and 10 women (18.51%) and animals included 15 sheep (27.77%) and three cattle (5.55%). Isolates were collected from three states of Iran, including Hamedan (59.25%), Arak (7.40%), and Tehran (33.33%). All isolates were identified either as Brucella melitensis or Brucella abortus. Three isolates (5.55%) were identified as Brucella abortus and 51 isolates (94.44%) were identified Brucella melitensis. The mean age of the 36 patients was 45, ranging from 6 to 80 years-old and male to female ratio was 2.6.

MLVA-16 genotyping results

The complete MLVA16 assay was performed targeting panel 1, 2A and 2B loci. The repeat unit size of panel 1 loci was ≥ 9 bp, while that of panel 2A and 2B was up to 8 bp. The PCR products for 16 loci were converted to copy number of TRs. Based on the published data, the polymorphism of each 16 VNTR loci was analyzed by the MLVA-16 method. Based on this method 54 isolates with genetic similarity coefficient of 80% were divided into 46 genotypes. Of those, 22 genotypes were singletons, while 4, 2, 1, and 2 genotypes were represented by 2, 3, 4, and 7 isolates, respectively (Fig. 1). The most prevalent genotype was revealed in 14 strains. There were two more frequent genotypes, each observed in 7 isolates, one genotype was restricted to one geographic region (Hamadan), whereas the other was present in two regions (Hamadan and Tehran). Genotypes present in two or three isolates were from Tehran and Hamadan provinces (Fig. 2). Discriminatory power was determined for each studied locus in this study and according to the results, panel 2B showed the highest discriminatory power [ Bruce-04 (0.837), Bruce-30 (0.806), Bruce-09 (0.787), Bruce-07 (0.772) and Bruce16-(0.766)] while panel 2A displayed a moderate variability [Bruce-18 (0.805), Bruce-19 (0.568), and Bruce-21 (0.358)], and panel 1 demonstrated a limited diversity [Bruce-06 (0.204), Bruce-08 (0.654), Bruce-11 (0.673), Bruce-12 (0.598), Bruce-42 (0.408), Bruce-43 (0.541), Bruce-45 (0.713), Bruce-55 (0.730)].

Fig. 1.

Dendrogram based on MLVA-16 assays for 54 isolates of Brucella melitensis and Brucella abortus. The columns represent the number of strains, geographical origin and bacterial sources.

Fig. 2.

Minimum spanning tree for Brucella isolates using MLVA-16 data.


In the current study a total of 54 Brucella isolates were collected from human and animal samples from three regions of Iran and MLVA-16 was used to assess the genetic diversity among the isolates. Thirty-six (66.66%) isolates were from humans and 18 (33.33%) from animals. Three isolates (5.55%) were identified as Brucella abortus and the remaing 51 were Brucella melitensis. Brucellosis is an atypical zoonotic disease and constant monitoring of animals will largely contribute to the control human brucellosis.14,15 Human brucellosis is most often linked with animal husbandry or consuming unpasteurized milk.16–18 MLVA-16 is a suitable typing method for establishing an epidemiological relationship between isolates in an outbreak to avoid excessive classical epidemiological investigations.19,20 MLVA-16 with a genetic similarity coefficient of 0.8, yielded a total of 46 genotypes, among which 22 were singletons. The frequency of different MLVA genotypes varied in the three studied regions. Three B. melitensis isolates collected from two sheep samples (blood and spleen), one isolates from cattle (liver) and one from human (CSF) in Arak were distinctive in genotypes. Two B. abortus isolated from human blood and one from cattle blood was also distinctive in genotypes. The results of this study showed no correlation between the genotypes of B. abortus isolated from cattle and those isolated from men.

Fifty B. melitensis isolates from human blood in Hamadan shared the same genotype as B. melitensis isolates from sheep blood and B. abortus isolates from cattle spleen in Hamadan, suggesting a possible epidemiological connection among them. These results were consistent to the study of Foster et al. who reported the spread of Brucella from animals to other animals or humans, and the epidemiologic relationship among them.21 Four other genotypes isolated from human blood in Hamedan were very closely related and differed by a single repeat unit at one or two of the most variant loci. Variety in single or double loci may reflect the microevolution of B. melitensis isolates. Additionally, three B. melitensis strains isolated from the blood of three human patients in Hamadan and Tehran, two B. melitensis strains isolated from the blood of animals, two isolates from blood of a human patient and one B. melitensis strain isolated from the liver of sheep were completely matched. These results suggest transmission of Brucella from animals to humans. There are some control measurements and financial compensation for the slaughtered seropositive cattle. However, no such measurement exists for sheep or goat, and therefore, sheep infected with Brucella are one of the main sources for human and animal brucellosis in Iran. This finding is similar to that reported in study from Hai Jiang, China22,23. Similar to previous studies, the highest diversity discriminatory power was found for panel 2B (0.7936) followed by panel 2A (0.557) and panel 1 (0.4760). According to these results, the loci of panel 1 including minisatellite loci with repeat unit length above 9 bp are more conserved than those in panel 2A and panel 2B with more heterogeneous microsatellite loci. The highest discriminatory power was indicated for Bruce30 (0.806) of panel 2B followed by Bruce18 (0.805) of panel 2A.


The MLVA16 analysis on 54 Brucella isolates showed high polymorphism of genotypes. Only two genotypes, each observed in seven isolates, were related to one another and only one of these genotypes belonged to two separate regions.

Availability of data and materials

Data are available upon request from corresponding author.

Authors' contributions

SM initiated and supervised the study. SM and FM designed the experiments. FM and AN conducted the experiments. AN conducted the analysis. SM, FM and AN has written the manuscript. All authors have read and approved the manuscript.


There is no funding for the present study.

Ethics approval and consent to participate

This study was carried according to approved protocol by Ethics committees of Iran University of medical sciences.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.


Not applicable.

G. Vergnaud, Y. Hauck, D. Christiany, B. Daoud, C. Pourcel, I. Jacques, et al.
Genotypic Expansion Within the Population Structure of Classical Brucella Species Revealed by MLVA16 Typing of 1404 Brucella Isolates From Different Animal and Geographic Origins, 1974-2006.
Front Microbiol., 9 (2018), pp. 1545
E. Shevtsova, G. Vergnaud, A. Shevtsov, A. Shustov, K. Berdimuratova, K. Mukanov, et al.
Genetic Diversity of Brucella melitensis in Kazakhstan in Relation to World-Wide Diversity.
Front Microbiol., 10 (2019), pp. 1897
K. Mohseni, R. Mirnejad, V. Piranfar, S. Mirkalantari.
A Comparative Evaluation of ELISA, PCR, and Serum Agglutination Tests For Diagnosis of Brucella Using Human Serum.
Iran J Pathol., 12 (2017), pp. 371-376
S. Mirkalantari, A.H. Zarnani, M. Nazari, G.R. Irajian, N. Amirmozafari.
Brucella melitensis VirB12 recombinant protein is a potential marker for serodiagnosis of human brucellosis.
Ann Clin Microbiol Antimicrob., 16 (2017), pp. 8
S. Mirkalantari, N. Amirmozafari, B. Kazemi, G. Irajian.
Molecular cloning of virB12 gene of Brucella melitensis 16M strain in pET28a vector.
Asian Pac J Trop Med., 5 (2012), pp. 511-513
M. Vatankhah, N. Beheshti, S. Mirkalantari, N. Khoramabadi, H. Aghababa, M. Mahdavi.
Recombinant Omp2b antigen-based ELISA is an efficient tool for specific serodiagnosis of animal brucellosis.
Braz J Microbiol., 50 (2019), pp. 979-984
M.N. Xavier, ÉA. Costa, T.A. Paixão, R.L. Santos.
The genus Brucella and clinical manifestations of brucellosis.
Ciência Rural., 39 (2009), pp. 2252-2260
Organization WH.
Brucellosis in humans and animals.
World Health Organization, (2006),
E. Thorne.
In Infectious Diseases of Wild Mammals, 3rd ed, pp. 372-395
S. Al Dahouk, P.L. Fleche, K. Nockler, I. Jacques, M. Grayon, H.C. Scholz, et al.
Evaluation of Brucella MLVA typing for human brucellosis.
J Microbiol Methods., 69 (2007), pp. 137-145
E. Tcherneva, N. Rijpens, B. Jersek, L.M. Herman.
Differentiation of Brucella species by random amplified polymorphic DNA analysis.
J Appl Microbiol., 88 (2000), pp. 69-80
B.J. Adams, N.P. Herrmann.
Estimation of living stature from selected anthropometric (soft tissue) measurements: applications for forensic anthropology.
J Forensic Sci., 54 (2009), pp. 753-760
M.M. Kattar, R.F. Jaafar, G.F. Araj, P. Le Fleche, G.M. Matar, R. Abi Rached, et al.
Evaluation of a multilocus variable-number tandem-repeat analysis scheme for typing human Brucella isolates in a region of brucellosis endemicity.
J Clin Microbiol., 46 (2008), pp. 3935-3940
J. Godfroid.
Brucellosis in livestock and wildlife: zoonotic diseases without pandemic potential in need of innovative one health approaches.
Arch Public Health., 75 (2017), pp. 34
M.Z. Khan, M. Zahoor.
An Overview of Brucellosis in Cattle and Humans, and its Serological and Molecular Diagnosis in Control Strategies.
Trop Med Infect Dis., 3 (2018),
K.A. Franc, R.C. Krecek, B.N. Hasler, A.M. Arenas-Gamboa.
Brucellosis remains a neglected disease in the developing world: a call for interdisciplinary action.
BMC Public Health., 18 (2018), pp. 125
N.C. Hull, B.A. Schumaker.
Comparisons of brucellosis between human and veterinary medicine.
Infect Ecol Epidemiol., 8 (2018),
M. Dadar, Y. Shahali, A.M. Whatmore.
Human brucellosis caused by raw dairy products: A review on the occurrence, major risk factors and prevention.
Int J Food Microbiol., 292 (2019), pp. 39-47
G. Garofolo, E. Di Giannatale, F. De Massis, K. Zilli, M. Ancora, C. Camma, et al.
Investigating genetic diversity of Brucella abortus and Brucella melitensis in Italy with MLVA-16.
Infect Genet Evol., 19 (2013), pp. 59-70
Z.G. Liu, Wang M. Di DD, R.H. Liu, H.Y. Zhao, D.R. Piao, et al.
MLVA Genotyping Characteristics of Human Brucella melitensis Isolated from Ulanqab of Inner Mongolia.
China. Front Microbiol., 8 (2017), pp. 6
J. Foster, F. Walker, B. Rannals, M. Hussain, K. Drees, et al.
African Lineage Brucella melitensis Isolates from Omani Livestock.
Front. Microbiol., (2018),
H. Jiang, M. Fan, J. Chen, et al.
MLVA genotyping of Chinese human Brucella melitensisbiovar 1, 2 and 3 isolates.
G. Garofolo, A. Fasanella, E. Di Giannatale, et al.
Cases of human brucellosis in Sweden linked to Middle East and Africa.
The Brazilian Journal of Infectious Diseases
Article options