MLVA genotyping of Chinese human Brucella melitensisbiovar 1, 2 and 3 isolates
- Hai Jiang†1,
- Mengguang Fan†2,
- Jingdiao Chen3,
- Jingchuan Mi2,
- Ruiping Yu2,
- Hongyan Zhao1,
- Dongri Piao1,
- Changwen Ke3,
- Xiaoling Deng3,
- Guozhong Tian1 and
- Buyun Cui1Email author
© Jiang et al; licensee BioMed Central Ltd. 2011
Received: 6 May 2011
Accepted: 22 November 2011
Published: 22 November 2011
Since 1950, Brucella melitensis has been the predominant strain associated with human brucellosis in China. In this study we investigated the genotypic characteristics of B. melitensis isolates from China using a multiple-locus variable-number tandem-repeat analysis (MLVA) and evaluated the utility of MLVA with regards to epidemiological trace-back investigation.
A total of 105 B. melitensis strains isolated from throughout China were divided into 69 MLVA types using MLVA-16. Nei's genetic diversity indices for the various loci ranged between 0.00 - 0.84. 12 out 16 loci were the low diversity with values < 0.2 and the most discriminatory markers were bruce16 and bruce30 with a diversity index of > 0.75 and containing 8 and 7 alleles, respectively. Many isolates were single-locus or double-locus variants of closely related B. melitensis isolates from different regions, including the north and south of China. Using panel 1, the majority of strains (84/105) were genotype 42 clustering to the 'East Mediterranean' B. melitensis group. Chinese B. melitensis are classified in limited number of closely related genotypes showing variation mainly at the panel 2B loci.
The MLVA-16 assay can be useful to reveal the predominant genotypes and strain relatedness in endemic or non-endemic regions of brucellosis. However it is not suitable for biovar differentiation of B. melitensis. Genotype 42 is widely distributed throughout China during a long time. Bruce 16 and bruce 30 in panel 2B markers are most useful for typing Chinese isolates.
Brucellosis, recognized as a common zoonotic disease globally, is caused by bacteria of the genus Brucella. B. melitensis, B. abortus, and B. suis remain the principal causes of human brucellosis worldwide and are major public health problems, primarily in Africa, the Middle East and Southeast Asia . Brucellosis is prevalent in China, especially in the northern China, where people are economically dependent on ruminant livestock. Approximately 30,000 human cases are reported annually over the past 5 years . In China, B. melitensis was the predominant strain associated with human brucellosis outbreaks, according to annual report on surveillance of selected infectious disease and vector. Species identification and subtyping of Brucella isolates is very important for epidemiologic surveillance and investigation of outbreaks in Brucella-endemic regions [3, 4].
Recent studies have confirmed that multiple-locus variable-number tandem-repeat analysis (MLVA) is a useful tool for identifying and genotyping Brucella strains and the resultant data can be used for epidemiological trace-back investigations [3, 5–8]. In efforts to better improve surveillance and evaluate the power of epidemiological trace-back in China, the MLVA-16 scheme was used to type a collection of 105 B. melitensis isolates from 18 different regions throughout China.
(This study was presented in part at the 5th Brucellosis International Research Conference of the American Society for Microbiology, Buenos Aires, Argentina, 2011.)
Typing and clustering of B. melitensisisolates by MLVA-16
Main characteristics of 16 VNTR loci in 105 B. melitensis isolates
No. of alleles
No. of repeats
Nei's DI and 95%CI
Analysis of the isolates from Inner Mongolia and Guangdong
Using the complete MLVA-16 assay, the 26 B. melitensis isolates from Inner Mongolia and 39 isolates from Guangdong were clustered in 20 and 27 different genotypes, respectively. The bruce16 loci had 6 and 7 alleles and the bruce30 loci had 6 and 5 alleles in these two population. 7 genotypes of 11 isolates from Inner Mongolia and 9 genotypes of 16 isolates from Guangdong isolates were single-locus or double locus variants of the main subgroups at 90% similarity. Only two of the clusters consisted of epidemiologically related strains exclusively from the same province and the same time. Bru0501 and bru0502 isolates came from a married couple who contacted with sheep fetuses and placenta in May, 1991 from Bai country, Xilinguole, Inner Mongolia. Bru0478, br0479, bru0480 and bru0481 isolates were collected over six weeks from late April to early June, 2010, in Pingsha county, Zhuhai city, Guangdong province. These four patients had not identified at-risk occupational history and experienced fever, debility and joint pains.
Trace-back of a laboratory-acquired Brucellainfection
We report a case of brucellosis affecting a hospital microbiology laboratory technician in Beijing, a non-endemic area of China. To better elucidate the origin of such infection, Brucella strains from both the patient and the laboratory technician were characterized by MLVA-16. The strain BJ06-10 showed the same MLVA type with strain NM06-11 isolated from a patient with acute brucellosis who engaged in fur-making in Inner Mongolia.
Identification of the B. melitensisvaccine strain M5
LB10-01, a B. melitensis biovar 1 strain isolated from Guangdong in 2010 was indistinguishable from the vaccine strain M5 according to the MLVA cluster analysis (MLVA027: 1-5-3-13-2-2-3-2-4-20-8-6-4-3-7-5). This is unexpected since the vaccine strain M5 was not used in Guangdong. Detection of a strain with phenotypic and genotypic properties indistinguishable from the vaccine strain M5 raises the concern of the origin of the wild type strain.
Brucellosis surveillance was started in 1980 in some parts of China. In 2008, 21 surveillance points for animal and human brucellosis were established in the 19 provinces of Heilongjiang, Jilin, Hebei, Henan, Inner Mongolia, Shandong, Guangdong, Guangxi, Sichuan, Tibet, Gansu, Ningxia, Xinjiang, Shanxi, Shan'xi, Zhejiang, Liaoning, Ningxia and Yunnan. Since the established of these surveillance points more than 30 years ago, a huge panel of animals and humans strains have been surveyed. It is significant that the national epidemiological characteristics can be analyzed. It suggests that B. melitensis isolates from different locations and years would reflect the epidemic features of human brucellosis.
Sheep infected with Brucella are one of the main sources for human and animal brucellosis in China . Over the last 20 years, the geographic distribution of brucellosis in China had been changing from pasturing areas to regions of with reduced agricultural interests (or alternatively more industrial concentrations); in these areas the infection rates, reported incidence, and number of outbreaks of brucellosis have increased markedly based on the National Notifiable Disease Surveillance System data. During this period, the cases have mostly been reported from Inner Mongolia, Shanxi, Hebei, Shandong, Henan, Liaoning, Jilin, Heilongjiang, and Shan'xi provinces. It is worth noting that brucellosis is endemic in Guangdong province, one of the wealthiest and industrial provinces in China. This is because of the movement of infected animal to Guangdong, resulting in the change of the geographic distribution of brucellosis.
In the different epidemic regions of China, the predominant strains have been shown to be B. melitensis biovar 1, 2, or 3 . Since 2005, most human cases in China have been caused by B. melitensis biovar 3 . Classical typing systems are unable to subdivide Brucella isolates below the biovar level. Molecular typing methods such as MLVA have been utilized to distinguish between strains of the same biovar in both animal and human isolates [3, 5, 6, 11–13]. In an effort to assess the value of MLVA as a subtyping tool for Brucella strains, genotypic characteristics of 105 B. melitensis isolates were investigated. Cluster analysis of these China strains, based on the eight variable-nucleotide tandem repeat loci included in the MLVA-16 panel 1 grouped them all into the B. melitensis 'East Mediterranean group'  and unique from circulating strains in Northern Africa, Southern Europe ('West Mediterranean group' and 'American group'). For instance, an (panel 1 genotype 42 and 43) clustered separately from most of the other 'West Mediterranean group' (panel 1 genotype 49 and 51) and 'American group'(panel 1 genotype 47). Previous studies have shown that Near Eastern countries frequently report human cases associated with genotypes 42 and 43 [3, 14]. Genotype 42, as we have shown, is widely distributed throughout China, and has previously been reported to be predominant in Turkey, Portugal and Spain . In Spain, human B. melitensis strains clustered into genotypes 42 (Eastern Mediterranean group, 55%), 48 and 53 (Americas group, ~11%) and 51 (Western Mediterranean group, ~8%). Chinese B. melitensis are classified in limited number of closely related genotypes showing variation mainly at the panel 2B loci.
In China, the Inner Mongolia Autonomous Region is the most severe endemic focus of brucellosis, with an annual incidence of the disease varying from 40 to 70/100,000 during 2005-2010 . Inner Mongolia is in close proximity to Heilongjiang, Jilin, Hebei and Shanxi provinces; these provinces are located in the north and east of China, where stocking raising is the most important aspect of the economy. In these regions, B. melitensis genotype 42 strains were predominant, but genotype 42 strains were also common in provinces reporting sporadic cases such as Liaoning, Shandong, Zhejiang, Fujian and Tianjin. These isolates were only single-locus or double-locus variants of B. melitensis from the endemic regions. Of particular note is the apparently stability of genotype 42 in China; genotype 42 strains were isolated from Inner Mongolia in1957 as well as 53 years later.
Guangdong province, which is now considered to be an endemic region for brucellosis, is located in the southern coastal region of China, where the incidence of human brucellosis has increased gradually since 2000. The prevailing panel 1 type is genotype 42 as well. The genotypes for most of the B. melitensis isolates in this series and their close relatedness by MLVA (single-locus variants and in some cases double-locus variants) compared to the relatedness of B. melitensis isolates in other countries reflects microevolution within the regionally important strain due to a few mutational events . We also observed that strains from the north and east of China (eg., Inner Mongolia and Shanxi) had the same MLVA-16 genotype (010) as those from the south of China (eg., Guangdong). This data indicates that the emergence of brucellosis in the south of China is likely to have its origins from the importation of animals from elsewhere in China. The clustering of epidemiologically-related isolates identified in the current and previous studies support the use of MLVA-16 as a valuable tool for investigations of outbreaks of both human and animal brucellosis. In our study, only 4 of 105 isolates (3.8%) had MLVA-16 genotype 030. It is likely that these cases represented a common-source outbreak or infected the herds of the same genotype. Because consistent epidemiological information for the strains is not routinely available, it is impossible to assess the relationship of the cluster results for these data and outbreaks. An urgent integrated, laboratory-based surveillance is needed to address this important public health gap.
To facilitate outbreak investigation, it has been recommended to use an abbreviated MLVA scheme, omitting testing with panel 1 and 2A since panel 2B is highly polymorphic and potentially more discriminating in determining genetic relationships in regions of endemicity . Some apparently unlinked (epidemiologically or otherwise) isolates had identical MLVA-16 profiles also. This led us to hypothesize that these may represent either epidemiologically unrelated isolates with homoplasy at MLVA-16 loci (most likely panel 2B) or persistent circulating strains causing sporadic infections [3, 14]. More detailed genetic investigations such as whole genome sequence comparison, should clarify these relationships.
Results of genotyping confirmed a laboratory-acquired Brucella infection. Laboratory workers who handle infected specimens are at high risk of acquiring Brucella infection, as suggested by the numerous cases of laboratory-acquired brucellosis reported in the literature . We report a case of brucellosis affecting a hospital microbiology laboratory technician in Beijing, a non-endemic area of China.
Human infection with the vaccine strain M5 in China has not been reported. However, in the previous reports, strains were only biotyped using conventional methods and no direct molecular linkage was shown between the isolated and commercial M5 vaccine strain. In this study, LB 10-01 has the identical genotype with M5. This suggests that LB 10-01 might be that a wild-type biovar 1 evolved with a pattern identical to M5 or that the original strain from which M5 was developed still is transmitted. Results obtained by Garcia-Yoldi et al. confirmed B. melitensis vaccine strain Rev 1 group as assayed by MLVA is genetically very homogeneous .
In conclusion, the MLVA assay is very promising to be used as a simple molecular tool for genotypes distribution of Brucella isolates from an endemic region, and might be useful for trace back investigations in non-endemic areas.
Bacterial strains and DNA preparation
The 105 B. melitensis isolates examined in this study
No. of isolates
42, 43, 63, CN-1
MLVA-16 genotyping scheme
MLVA was performed as previously described . The sixteen primer pairs were divided into three groups as previously described: panel 1 (8 loci including bruce06, bruce08, bruce11, bruce12, bruce42, bruce43, bruce45, and bruce55), panel 2A (3 loci including bruce18, bruce19, and bruce21), and panel 2B (5 loci including bruce04, bruce07, bruce09, bruce16, and bruce30). PCR conditions were as follows: initial denaturation at 94°C for 3 min, and then 30 cycles of 94°C for 30 s, 60°C for 30 s and 72°C for 50 s. Five microliters of the amplification products were loaded in to 2% (panel 1) and 3% (panels 2A and 2B) agarose gels containing ethidium bromide (0.5 μg/ml), visualized under UV light, and photographed. The reference strain B. melitensis 16 M, for which the precise molecular mass is known for each primer pair locus, was used for size comparison. To determine the number of repeats from the sample products, PCR products were purified and directly sequenced using an ABI Prism Big Dye Terminator (v3.1) cycle sequencing ready reaction kit (v5.0). The PCR products of samples were sequenced and the sequences were compared to that of B. melitensis 16 M.
Analysis of MLVA data
All data were analyzed using BioNumerics version 5.1 software (Applied Maths, Belgium). Clustering analysis was based on the categorical coefficient and unweighted pair group method using arithmetic averages (UPGMA) method. Polymorphism at each loci was quantified using Nei's diversity index, available in the website of HPA http://www.hpa-bioinformatics.org.uk/cgi-bin/DICI/DICI.pl. Resultant genotypes were compared using the web-based Brucella2010 MLVA database http://mlva.u-psud.fr/.
We thank John Klena for his assistance in improving this manuscript. We also gratefully thank Haijian Zhou for clustering analysis.
This study was funded by the National Basic Research Program (2010CB530201) and National High Technology Research and Development Program (2007AA02Z410) from Ministry of Science and Technology of the People's Republic of China.
- 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.PubMedView ArticleGoogle Scholar
- Zhang WY, Guo WD, Sun SH, Jiang JF, Sun HL, Li SL, Liu W, Cao WC: Human brucellosis, Inner Mongolia, China. Emerg Infect Dis. 2010, 16: 2001-2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Al DS, Fleche PL, 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.View ArticleGoogle Scholar
- Marianelli C, Graziani C, Santangelo C, Xibilia MT, Imbriani A, Amato R, Neri D, Cuccia M, Rinnone S, Di MV, 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Her M, Kang SI, Cho DH, Cho YS, Hwang IY, Heo YR, Jung SC, Yoo HS: Application and evaluation of the MLVA typing assay for the Brucella abortus strains isolated in Korea. BMC Microbiol. 2009, 9: 230-10.1186/1471-2180-9-230.PubMedPubMed CentralView ArticleGoogle Scholar
- Her M, Kang SI, Kim JW, Kim JY, Hwang IY, Jung SC, Park SH, Park MY, Yoo HS: A genetic comparison of Brucella abortus isolates from animals and humans by using an MLVA assay. J Microbiol Biotechnol. 2010, 20: 1750-1755.PubMedGoogle Scholar
- Kang SI, Heo EJ, Cho D, Kim JW, Kim JY, Jung SC, Her M: Genetic Comparison of Brucella canis Isolates by the MLVA Assay in South Korea. J Vet Med Sci. 2011Google Scholar
- Smits HL, Espinosa B, Castillo R, Hall E, Guillen A, Zevaleta M, Gilman RH, Melendez P, Guerra C, Draeger A, Broglia A, Nockler K: MLVA genotyping of human Brucella isolates from Peru. Trans R Soc Trop Med Hyg. 2009, 103: 399-402. 10.1016/j.trstmh.2008.10.004.PubMedView ArticleGoogle Scholar
- Shang DQ, Xiao DL, Yin JM: Epidemiology and control of brucellosis in China. Vet Microbiol. 2002, 90: 165-182. 10.1016/S0378-1135(02)00252-3.View ArticleGoogle Scholar
- Cui BY: Endemic surveillance and control of Brucellosis in China. Disease Surveillance. 2007, 22: 649-651.Google Scholar
- Le FP, Jacques I, Grayon M, Al DS, Bouchon P, Denoeud F, Nockler K, Neubauer H, Guilloteau LA, Vergnaud G: Evaluation and selection of tandem repeat loci for a Brucella MLVA typing assay. BMC Microbiol. 2006, 6: 9-10.1186/1471-2180-6-9.View ArticleGoogle Scholar
- Tiller RV, De BK, Boshra M, Huynh LY, Van Ert MN, Wagner DM, Klena J, Mohsen TS, El-Shafie SS, Keim P, Hoffmaster AR, Wilkins PP, Pimentel G: Comparison of two multiple-locus variable-number tandem-repeat analysis methods for molecular strain typing of human Brucella melitensis isolates from the Middle East. J Clin Microbiol. 2009, 47: 2226-2231. 10.1128/JCM.02362-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Valdezate S, Navarro A, Villalon P, Carrasco G, Saez-Nieto JA: Epidemiological and phylogenetic analysis of Spanish human Brucella melitensis strains by multiple-locus variable-number tandem-repeat typing, hypervariable octameric oligonucleotide fingerprinting, and rpoB typing. J Clin Microbiol. 2010, 48: 2734-2740. 10.1128/JCM.00533-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Kattar MM, Jaafar RF, Araj GF, Le FP, Matar GM, Abi RR, Khalife S, 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, 46: 3935-3940. 10.1128/JCM.00464-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Marianelli C, Petrucca A, Pasquali P, Ciuchini F, Papadopoulou S, Cipriani P: Use of MLVA-16 typing to trace the source of a laboratory-acquired Brucella infection. J Hosp Infect. 2008, 68: 274-276. 10.1016/j.jhin.2008.01.003.PubMedView ArticleGoogle Scholar
- Garcia-Yoldi D, Le FP, Marin CM, de Miguel MJ, Munoz PM, Vergnaud G, Lopea-Goni I: Assessment of genetic stability of Brucella melitensis Rev 1 vaccine strain by multiple-locus variable-number tandem repeat analysis. Vaccine. 2007, 25: 2858-2862. 10.1016/j.vaccine.2006.09.063.PubMedView ArticleGoogle Scholar
- Alton GG, Jones LM, Pietz DE: Laboratory techniques in brucellosis. Monogr Ser World Health Organ. 1975, 1-163.Google Scholar
- 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: 2660-2666.PubMedPubMed CentralGoogle Scholar
- Hunter PR, Gaston MA: Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. J Clin Microbiol. 1988, 26: 2465-2466.PubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.