MLVA and MLST typing of Brucella from Qinghai, China
- Jun-Ying Ma†1,
- Hu Wang†1,
- Xue-Fei Zhang1,
- Li-Qing Xu1,
- Gui-Ying Hu1,
- Hai Jiang2,
- Fang Zhao3,
- Hong-Yan Zhao2,
- Dong-Ri Piao2,
- Yu-Min Qin1,
- Bu-Yun Cui2Email author and
- Gong-Hua Lin3Email author
© Ma et al. 2016
Received: 21 October 2015
Accepted: 23 March 2016
Published: 13 April 2016
The Qinghai-Tibet Plateau (QTP) of China is an extensive pastoral and semi-pastoral area, and because of poverty and bad hygiene conditions, Brucella is highly prevalent in this region. In order to adequately prevent this disease in the QTP region it is important to determine the identity of Brucella species that caused the infection.
A total of 65 Brucella isolates were obtained from human, livestock and wild animals in Qinghai, a Chinese province in east of the QTP. Two molecular typing methods, MLVA (multi-locus variable-number tandem-repeat analysis) and MLST (multi locus sequence typing) were used to identify the species and genotypes of these isolates.
Both MLVA and MLST typing methods classified the 65 isolates into three species, B. melitensis, B. abortus and B. suis, which included 60, 4 and 1 isolates respectively. The MLVA method uniquely detected 34 (Bm01 ~ Bm34), 3 (Ba01 ~ Ba03), and 1 (Bs01) MLVA-16 genotypes for B. melitensis, B. abortus and B. suis, respectively. However, none of these genotypes exactly matched any of the genotypes in the Brucella2012 MLVA database. The MLST method identified five known ST types: ST7 and ST8 (B. melitensis), ST2 and ST5 (B. abortus), and ST14 (B. suis). We also detected a strain with a mutant type (3-2-3-2-?-5-3-8-2) of ST8 (3-2-3-2-1-5-3-8-2). Extensive genotype-sharing events could be observed among isolates from different host species.
There were at least three Brucella (B. melitensis, B. abortus and B. suis) species in Qinghai, of which B. melitensis was the predominant species in the area examined. The Brucella population in Qinghai was very different from other regions of the world, possibly owing to the unique geographical characteristics such as extremely high altitude in QTP. There were extensive genotype-sharing events between isolates obtained from humans and other animals. Yaks, sheep and blue sheep were important zoonotic reservoirs of brucellosis causing species found in humans.
KeywordsBrucella Molecular identification Genotype Zoonotic host
Please see Additional file 1 for translations of the abstract into the six official working languages of the United Nations.
Brucellosis is one of the most widespread and severe zoonotic diseases across the world. This disease affects mammals and is transmitted to humans by animals through direct contact with contaminated animal fluids or indirectly through consumption of unpasteurized milk products . Brucellosis in domestic mammals and humans is the cause of huge economic burden and medical challenges globally, especially in poor regions [2–5]. The causal agent of brucellosis, Brucella genus, currently contains 10 species: B. abortus, B. suis, B. melitensis, B. ovis, B. canis, B. neotomae, B. ceti, B. pinnipedialis, B. microti, and B. inopinata. Some species are further subdivided into biovars, i.e., B. melitensis bv.1 to 3, B. abortus bv. 1 to 6 and 9, and B. suis bv. 1 to 5 . Most species with a preferential host and three other species (B. abortus, B. melitensis, and B. suis) remain the principal causes of zoonotic potential .
Since different species or biovars within a species of Brucella have different epidemiological significances, in order to effectively prevent this disease, it is important to identify the species or strain of Brucella that causes infection. Most previous molecular subtyping tools and “classical biotyping” methods lack sufficient discriminatory power for epidemiological investigations. Recently, two molecular typing methods called MLVA (multi-locus variable-number tandem-repeat analysis) and MLST (multi locus sequence typing) have emerged as useful tools for identifying and genotyping Brucella isolates [3, 7–9].
Geographic distribution of Brucella samples in Qinghai, China
Ma37, Ma51, Ma52, Ma53, Ma54, Ma55, Ma56, Ma57, Ma58
Ma13, Ma16, Ma18, Ma19, Ma20, Ma21, Ma26,
Ma27, Ma28, Ma29, Ma31, Ma32, Ma33
Ma03, Ma12, Ma22, Ma24, Ma25, Ma30, Ma34, Ma35
Ma39, Ma40, Ma41, Ma42, Ma43, Ma44, Ma45, Ma46
Ma36, Ma59, Ma60, Ma61, Ma62, Ma63, Ma64
Ma01, Ma02, Ma05, Ma06, Ma08, Ma09, Ma15
Ma07, Ma14, Ma65
Ma47, Ma48, Ma49
This study is a retrospective investigation of our historical collection with modern typing methods. Therefore, the study does not involve the collection or reporting of patient data. No animal work was carried out with the obtained results.
Bacterial strains and DNA preparation
A total of 65 isolates were examined: 31 recovered from human, 15 from sheep (Ovis aries), 11 from blue sheep (Pseudois nayaur), 6 from yaks (Bos mutus grunniens), and 1 each from cattle (Bos primigenius taurus) and Tibetan gazelle (Procapra picticaudata). Bacterial strains were cultured on Trypticase soy agar containing 5 % sheep blood (BD Diagnostic Systems, China Ltd., China) at 37 °C for 48 h. Total genomic DNA was extracted using the DNeasy Blood & Tissue Kit (Qiagen, Germany) as per the manufacturer’s recommended protocol.
MLVA was performed as described earlier . Briefly, 16 primer pairs were designed and classified as belonging to one of the three panels: panel 1 (bruce06, bruce08, bruce11, bruce12, bruce42, bruce43, bruce45, and bruce55), panel 2A (bruce18, bruce19, and bruce21), and panel 2B (bruce04, bruce07, bruce09, bruce16, and bruce30). PCR amplifications were performed in 40 μL reaction volumes, 5 μL of these were loaded in to 2 % (panel 1) or 3 % (panels 2A and 2B) agarose gels containing 0.5 μg/ml of ethidium bromide, visualized under UV light, and photographed. Band intensities were estimated using BioNumerics version 5.1 (Applied Maths, Belgium) and were then converted to repeat units by following the published allele numbering system .
In order to improve the genotyping accuracy, PCR products were also directly sequenced after purification. The sequences were aligned and the numbers of repeat units were checked in MEGA version 5  using B. melitensis 16 M genome sequence as a reference (GenBank accession number NC_003317.1 and NC_003318.1). Clustering analysis was based on the categorical coefficient and unweighted pair group method using arithmetic averages (UPGMA) method provided in BioNumerics. The resulting genotypes were compared using the web-based Brucella2012 MLVA database (http://mlva.u-psud.fr/).
Nine distinct genomic loci, including seven housekeeping genes (gap, aroA, glk, dnaK, gyrB, trpE, and cobQ), one outer membrane protein (omp25), and one intergenic fragment (int-hyp), were selected for MLST genotyping . Similar to the MLVA process described, PCR amplifications were performed in 40 μL reaction volumes using primer sequences as previously described . Sequences obtained from purified PCR products were aligned using MEGA program and verified by visualizing. MLST sequences (GenBank accession number AM694191 ~ AM695630) of the strains described by Whatmore et al.  were downloaded and a local BLAST database was built using makeblastdb program of BLAST+ program package version 2.2.31 . The sequences were matched with the local database using blastn program of the BLAST+ program package.
The BLAST results were parsed using BLAST Parser program version 1.2.6 . Distinct alleles identified at the nine selected loci were each given a numerical designation according to sequence of defined alleles. If the sequence was different from those defined previously, it was designated as a new allele. Each sequence type over all loci (ST) was predicted using web-based MLST service (BrucellaBase, http://188.8.131.52/brucellabase/mlst.html). Phylogenetic relationships of combined sequences were inferred using the Neighbor-Joining method in BrucellaBase.
Based on a previous study , the B. melitensis strains were sub-grouped into “East Mediterranean” group (genotype 42 and 43) and “American” group (genotype 47). Considering all three MLVA panels, 34 (Bm01 ~ Bm34), 3 (Ba01 ~ Ba03), and 1 (Bs01) MLVA-16 genotypes were identified for B. melitensis, B. abortus and B. suis isolates respectively (Fig. 2). No genotypes completely matched any of the genotypes in the Brucella2012 MLVA database.
A total of five known MLST genotypes were identified: ST7 (3-5-3-2-1-5-2-10-2; N = 2), ST8 (3-2-3-2-1-5-3-8-2; N = 57), ST2 (2-1-2-2-1-3-1-1-1; N = 2), ST5 (2-1-1-2-1-4-1-1-1; N = 2), ST14 (1-6-4-1-4-3-5-2-1; N = 1). The Neighbor-Joining clustering showed that the 65 isolates formed four main clusters – a, d, e and f. Of these, clusters a, e and f correspond to genotypes ST8, ST7 and ST14 respectively while cluster d corresponds to two genotypes, ST2 and ST5. According to Whatmore et al. , genotypes ST7 and ST8 belong to B. melitensis, ST2 and ST5 belong to B. abortus, and ST14 belongs to B. suis.
The gyrB sequence of strain Ma06 was different from any of known alleles (allele 1 ~ 6). Its sequence closely matched that of allele 1 (identical except for a di-nucleotide change of GC to AA on site 405 ~ 406 relative to the gyrB allele 1 sequence or on 2504 ~ 2505 sites relative to the ST8 sequence). Considering all MLST alleles, the strain Ma06 could be viewed as a mutant variant (3-2-3-2-?-5-3-8-2) of ST8 (3-2-3-2-1-5-3-8-2).
Based on MLVA-16, ten genotypes were shared by two or more host species. To elaborate further, four genotypes (Bm09, Bm19, Bm20, and Bm25) were shared between humans and sheep, and two (Bm08 and Ba03) were shared between humans and yaks (Fig. 2). Based on MLST, the ST8 type was also shared between humans and four other animals (sheep, yaks, blue sheep and Tibetan gazelle); the ST5 type was shared between humans and yaks (Fig. 3). Sheep and yaks are the main livestock on QTP and it is therefore not surprising that these animals pose a heavy zoonotic disease threat to humans. It should be noted that, all sheep and four out of six yaks were infected by B. melitensis, indicating that this bacterial species should be paid special attention to in both host species. It is well known that different bacterial species need different vaccine strains. For instance, there are three vaccine strains M5, S19 and S2, which are used to immunize animals in China for B. melitensis, B. abortus, and B. suis, respectively . Because yaks could be infected by either B. melitensis or B. abortus, we suggest based on our study that both the B. melitensis type (M5) and B. abortus type vaccines (S19) should be tried in this animal.
The blue sheep were infected by all three Brucella species including B. suis, which is generally carried by pigs [1, 22]. More importantly, this animal shared Brucella genotypes with human (Bm07, Bm15 and Bm25), yaks (Bm06), and Tibetan gazelle (Bm01), indicating this animal is an important natural repertory for Brucella. Since there is obvious spatial and diet overlap between blue sheep and domestic livestock , we suggest that this animal is also an important infection source of Brucella in humans.
Additionally, due to high genetic diversity of B. melitensis, we also discuss the geographic distribution features of the strains. Figure 1 showed that 7 counties (locality ID: 2, 3, 5, 8, 9, 10, and 12) were dominated by A branch strains; 3 counties (locality ID: 1, 4, and 6) were dominated by B branch strains; 3 counties (locality ID: 7, 11, and 13) contained only C branch strains; and the E branch genotypes appeared only in 1 county (locality ID: 2). These results showed that the B. melitensis prevalence in the QTP region might have a considerable regional characteristic i.e., there might be genetic variations among different regions.
Our results for the first time elucidated the composition of species and genotypes of Brucella samples in the QTP region. We also presented the host- as well as geographic distributions of the species and genotypes. These results may have some implications for the future disease control programmes in QTP. First, since B. melitensis was the predominant species in the study area, special attention should be paid on this species in brucellosis control programs such as livestock vaccination. Second, the genotypes of Brucella in the QTP region were very different from other regions, but whether the genetic difference may result in different vaccination efficiencies should be elucidated in the further studies. Third, since there were extensive genotype-sharing events between humans and animals, brucellosis surveillance should be regularly executed on both livestock (sheep and yaks) and wild ungulates (blue sheep).
It should be noted that, as mentioned above, none of the genotypes found in our study were identical to any of the genotypes in the Brucella2012 database, in order to prevent over-interpreting our results, we did not compare our strains with those from other studies. This might limit our understanding of genetic relationships between the Brucella populations from QTP and from other regions. Moreover, because of imbalance of sample sizes among different counties and among different hosts, in this study, we were not able to do a detail epidemiological investigation. We suggest that the accumulation of more samples will enable us to further elucidate the genetic characteristics of Brucella species in QTP.
Our study drew three main conclusions: (i) there were at least three Brucella (B. melitensis, B. abortus and B. suis) species in the east of QTP and B. melitensis was the predominant species in the area studied; (ii) the Brucella population in the QTP region was very different from other regions probably due to the unique geographical characteristics e.g., extremely high altitude in QTP; and (iii) there were extensive genotype-sharing events between humans and animals; sheep, yaks and blue sheep were important zoonotic hosts of brucellosis to humans in the area studied.
multi locus sequence typing
multi-locus variable-number tandem-repeat analysis
This study was supported by the Science and Technology Plan of Qinghai Province (NO. 2013-Z-749) and the Youth Innovation Promotion Association of Chinese Academy of Sciences (NO. 2015352).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- 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;65:523–41.View ArticlePubMedGoogle Scholar
- Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos EV. The new global map of human Brucellosis. Lancet Infect Dis. 2006;6:91–9.View ArticlePubMedGoogle Scholar
- Whatmore AM, Perrett LL, MacMillan AP. Characterisation of the genetic diversity of Brucella by multilocus sequencing. BMC Microbiol. 2007;7:34.View ArticlePubMedPubMed CentralGoogle Scholar
- Kansiime C, Rutebemberwa E, Asiimwe BB, Makumbi F, Bazira J, Mugisha A. Annual trends of human brucellosis in pastoralist communities of south-western Uganda: a retrospective ten-year study. Infec Dis Poverty. 2015;4:39.View ArticleGoogle Scholar
- Ducrotoy MJ, Ammary K, Lbacha HA, Zouagui Z, Mick V, Prevost L, Bryssinckx W, Welburn SC, Benkirane A. Narrative overview of animal and human brucellosis in Morocco: intensification of livestock production as a driver for emergence? Infec Dis Poverty. 2015;4:57.View ArticleGoogle Scholar
- Garin-Bastuji B, Mick V, Le Carrou G, Allix S, Perrett LL, Dawson CE, et al. Examination of taxonomic uncertainties surrounding Brucella abortus bv. 7 by phenotypic and molecular approaches. Appl Environ Microbiol. 2014;80(5):1570–9.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:9.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Y, Zhen Q, Wang Y, Xu J, Sun Y, Li T, et al. Development of an extended multilocus sequence typing for genotyping of Brucella isolates. J Microbiol Methods. 2011;86(2):252–4.View ArticlePubMedGoogle Scholar
- Di D, Cui B, Wang H, Zhao H, Piao D, Tian L, et al. Genetic polymorphism characteristics of Brucella canis isolated in China. PLoS One. 2014;9(1):e84862.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu T. The Qinghai-Tibetan plateau: how high do Tibetans live? High Alt Med Biol. 2001;2:489–99.View ArticlePubMedGoogle Scholar
- Wiener G, Han J, Long R. The yak. 2nd ed. Bangkok: RAP publication; 2003.Google Scholar
- Lang X, Qu H, Bu Z, Yang Y, Sun C, Li X, et al. Seroprevalence of Brucella infection in yaks (Bos grunniens) on the Qinghai-Tibet plateau of China. Trop Anim Health Prod. 2011;43:305–6.View ArticleGoogle Scholar
- Xue H. Analysis on infection and the related factor of brucellosis in high-risk population in Qinghai Province of 2006 to 2010. Taiyuan, Shanxi Province, China: Master thesis of Shanxi Medical University; 2012.
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Y, Ke Y, Wang Y, Yuan X, Zhou X, Jiang H, et al. Changes of predominant species/biovars and sequence types of Brucella isolates, Inner Mongolia. China BMC Infect Dis. 2013;13:514.View ArticlePubMedGoogle Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.View ArticlePubMedPubMed CentralGoogle Scholar
- Marco C, Roberta B. BLAST Parser. BITS Conference. Societa di Bioinformatica Italiana. Bologna. Societa di Bioinformatica Italiana Bologna Italy. 2006;5:28–9.Google Scholar
- Jiang H, Fan M, Chen J, Mi J, Yu R, Zhao H, et al. MLVA genotyping of Chinese human Brucella melitensis biovar 1, 2 and 3 isolates. BMC Microbiol. 2011;11:256.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Al Dahouk S, Le Flèche P, Nöckler P, Jacques I, Grayon M, Scholz HC, et al. Evaluation of Brucella MLVA typing for human brucellosis. J Microbiol Methods. 2007;69:137–45.View ArticlePubMedGoogle Scholar
- Shang D, Xiao D, Yin J. Epidemiology and control of brucellosis in China. Vet Microbiol. 2002;90:165–82.View ArticleGoogle Scholar
- Burniston S, Okello AL, Khamlome B, Inthavong P, Gilbert J, Blacksell SD, Allen J, Welburn SC. Cultural drivers and health-seeking behaviours that impact on the transmission of pig-associated zoonoses in Lao People’s Democratic Republic. Infec Dis Poverty. 2015;4:11.View ArticleGoogle Scholar
- Mishra C, Van Wieren SE, Ketner P, Heitkonig IMA, Prins HHT. Competition between domestic livestock and wild bharal Pseudois nayaur in the Indian Trans-Himalaya. J Appl Ecol. 2004;41(2):344–54.View ArticleGoogle Scholar