MLVA and MLST typing of Brucella from Qinghai, China

Similar documents
Genetic polymorphisms identify in species/ biovars of Brucella isolated in China between 1953 and 2013 by MLST

Microbiología, Instituto de Salud Carlos III, Majadahonda, Madrid, Spain.

MOLECULAR EPIDEMIOLOGY OF BRUCELLA MELITENSIS STRAINS CAUSING OUTBREAKS IN CROATIA AND BOSNIA AND HERZEGOVINA

Overview of animal and human brucellosis in EU: a controlled disease?

Recent Topics of Brucellosis

A rapid minor groove binder PCR method for distinguishing the vaccine strain Brucella abortus 104M

Link between geographical origin and occurrence of Brucella abortus biovars in cattle

Multiple-Locus Variable-Number Tandem-Repeat Analysis Genotyping of Human Brucella Isolates from Turkey

Occurrence, molecular characterization and predominant genotypes of Enterocytozoon bieneusi in dairy cattle in Henan and Ningxia, China

Trends of reported human brucellosis cases in mainland China from 2007 to 2017: an exponential smoothing time series analysis

OIE RL for Rabies in China: Activities and Challenges

of Conferences of OIE Regional Commissions organised since 1 June 2013 endorsed by the Assembly of the OIE on 29 May 2014

ERG on multidrug-resistant P. falciparum in the GMS

Promoting rational antibiotic prophylaxis in clean surgeries in China

Food safety related to camelids products: Brucellosis and its impact on Public Health and the consumers as an example

Drd. OBADĂ MIHAI DORU. PhD THESIS ABSTRACT

Sylvia Valdezate,* Ana Navarro, Pilar Villalón, Gema Carrasco, and Juan A. Saéz-Nieto

Rapid molecular testing to detect Staphylococcus aureus in positive blood cultures improves patient management. Martin McHugh Clinical Scientist

Comparing DNA Sequences to Understand Evolutionary Relationships with BLAST

The melanocortin 1 receptor (mc1r) is a gene that has been implicated in the wide

EFSA s activities on Antimicrobial Resistance

Impact of neglected diseases on animal productivity and public health in Africa

IMPACT OF NEGLECTED DISEASES ON ANIMAL PRODUCTIVITY AND PUBLIC HEALTH IN AFRICA

COMPARING DNA SEQUENCES TO UNDERSTAND EVOLUTIONARY RELATIONSHIPS WITH BLAST

Classificatie: intern

FAO-APHCA/OIE/USDA Regional Workshop on Prevention and Control of Neglected Zoonoses in Asia July, 2015, Obihiro, Japan.

EPIDEMIOLOGY OF CAMPYLOBACTER IN IRELAND

Co-transfer of bla NDM-5 and mcr-1 by an IncX3 X4 hybrid plasmid in Escherichia coli 4

Index. Note: Page numbers of article titles are in boldface type.

Downloaded from irje.tums.ac.ir at 8:43 IRST on Sunday February 17th 2019

PCR detection of Leptospira in. stray cat and

The OIE Manual of Diagnostic Tests and Vaccines for Terrestrial & Aquatic Animals

Jefferson County High School Course Syllabus

MULTILOCUS SEQUENCE TYPING OF BRUCELLA ISOLATES FROM THAILAND

Background 1 st, 2 nd and 3 rd FAO-APHCA/OIE Regional Workshop on Brucellosis Diagnosis and Control with an Emphasis on Brucella melitensis (in

Required and Recommended Supporting Information for IUCN Red List Assessments

Garin-Bastuji. In terms of research and development, the work of the Unit concerns:

Antibiotic resistance of bacteria along the food chain: A global challenge for food safety

European Regional Verification Commission for Measles and Rubella Elimination (RVC) TERMS OF REFERENCE. 6 December 2011

Comparing DNA Sequence to Understand

The Identification of the Cryptosporidium ubiquitum in Pre weaned Ovines from Aba Tibetan and Qiang Autonomous Prefecture in China*

OIE Reference Laboratory Reports Activities

Global Coordination of Animal Disease Research. Alex Morrow

Brucellosis situation in Mongolia and Result of Bovine Brucellosis Proficiency Test

ANTIBIOTICS AND ANTIMICROBIAL RESISTANCE: CAUSES AND POSSIBLE SOLUTIONS

Structure of the OIE Manual of Diagnostic Tests and Vaccines for Terrestrial Animals

Development and improvement of diagnostics to improve use of antibiotics and alternatives to antibiotics

Introduction to Biorisk and the OIE Standard

TRYPANOSOMIASIS IN TANZANIA

Molecular Characterization of Mycoplasma agalactiae. Reveals the Presence of an Endemic Clone in Spain

Brucella abortus continues to be an important veterinary pathogen

Current status of the evaluation of genetic diversity in livestock breeds

The Rufford Foundation Final Report

Bovine Brucellosis Control of indirect ELISA kits

Received 9 December 2008/Returned for modification 26 April 2009/Accepted 6 May 2009

Comparing DNA Sequences Cladogram Practice

11-ID-10. Committee: Infectious Disease. Title: Creation of a National Campylobacteriosis Case Definition

FAO Initiatives and Protocols on Brucellosis and Tuberculosis Prevention and Control in Animals

Breeding Icelandic Sheepdog article for ISIC 2012 Wilma Roem

Project Summary. Emerging Pathogens in US Cattle

RICKETTSIA SPECIES AMONG TICKS IN AN AREA OF JAPAN ENDEMIC FOR JAPANESE SPOTTED FEVER

SNP genotypes of olfactory receptor genes associated with olfactory ability in German Shepherd dogs

EVOLUTIONARY GENETICS (Genome 453) Midterm Exam Name KEY

Dog vaccination with EgM proteins against Echinococcus granulosus

Antimicrobial resistance (EARS-Net)

Development and Characterization of Mouse Models of Infection with Aerosolized Brucella melitensis and Brucella suis

o VETERINARY IMMUNODIAGNOSTICS MARKET- GLOBAL OPPORTUNITY ANALYSIS AND INDUSTRY FORECASTS TO 2022 Report ID: MRAM Publishing Date: July, 2017

Disease Ecology: The role of global change on emerging infectious diseases

Within-host evolution of Brucella canis during a canine brucellosis outbreak in a kennel

OIE Reference Laboratory Reports Activities

Curriculum Vitae. : AlBaha University, faculty of Science.

Antibiotic Resistance

Monitoring gonococcal antimicrobial susceptibility

Annex III : Programme for the control and eradication of Transmissible Spongiform Encephalopathies submitted for obtaining EU cofinancing

Understanding and prevention of transmission of antibiotic resistance between bacterial populations and One Health reservoirs

2015 Work Programme of the

Presence of extended spectrum β-lactamase producing Escherichia coli in

QUANTITATIVE AND QUALITATIVE IMPROVEMENT OF THE SHEEP MUTTON PRODUCTION WITH THE HELP OF MOLECULAR MARKER AND GENOME EDITING TECHNOLOGY : A REVIEW

Relative effectiveness of Irish factories in the surveillance of slaughtered cattle for visible lesions of tuberculosis,

1 In 1958, scientists made a breakthrough in artificial reproductive cloning by successfully cloning a

Epidemic characteristics, high-risk townships and space-time clusters of human brucellosis in Shanxi Province of China,

A LABORATORY NETWORK FOR DIAGNOSTIC OF CAMELIDS DISEASES

OIE Reference Laboratory Reports Activities

Strengthening Epidemiology Capacity Using a One Health Framework in South Asia

OIE international standards on Rabies:

Received 7 December 1998/Returned for modification 5 April 1999/Accepted 22 June 1999

Centre for Public Health Research Laboratories

Campylobacter infections in EU/EEA and related AMR

UPDATE ON DEMONSTRATED RISKS IN HUMAN MEDICINE FROM RESISTANT PATHOGENS OF ANIMAL ORIGINS

EUROPEAN COMMISSION HEALTH & CONSUMERS DIRECTORATE-GENERAL. Unit G5 - Veterinary Programmes

Ticks and tick-borne pathogens Jordi Tarrés-Call, Scientific Officer of the AHAW unit

Report by the Director-General

Boosting Biomedical Entity Extraction by Using Syntactic Patterns for Semantic Relation Discovery

A systematic review of zoonoses transmission and livestock/wildlife interactionspreliminary

AKC Canine Health Foundation Grant Updates: Research Currently Being Sponsored By The Vizsla Club of America Welfare Foundation

Genes What are they good for? STUDENT HANDOUT. Module 4

Bioinformatics: Investigating Molecular/Biochemical Evidence for Evolution

Outcome of the Conference Towards the elimination of rabies in Eurasia Joint OIE/WHO/EU Conference

Lecture 11 Wednesday, September 19, 2012

Inheritance of Livershunt in Irish Wolfhounds By Maura Lyons PhD

Transcription:

Ma et al. Infectious Diseases of Poverty (0) : DOI 0./s0-0-0-z SHORT REPORT Open Access MLVA and MLST typing of Brucella from Qinghai, China Jun-Ying Ma, Hu Wang, Xue-Fei Zhang, Li-Qing Xu, Gui-Ying Hu, Hai Jiang, Fang Zhao, Hong-Yan Zhao, Dong-Ri Piao, Yu-Min Qin, Bu-Yun Cui * and Gong-Hua Lin * Abstract Background: 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. Methods: A total of Brucella isolates were obtained from, 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. Findings: Both MLVA and MLST typing methods classified the isolates into three species, B. melitensis, B. abortus and B. suis, which included 0, and isolates respectively. The MLVA method uniquely detected (Bm0 ~ Bm), (Ba0 ~ Ba0), and (Bs0) MLVA- genotypes for B. melitensis, B. abortus and B. suis, respectively. However, none of these genotypes exactly matched any of the genotypes in the Brucella0 MLVA database. The MLST method identified five known ST types: ST and ST (B. melitensis), ST and ST (B. abortus), and ST (B. suis). We also detected a strain with a mutant type (----?----) of ST (--------). Extensive genotype-sharing events could be observed among isolates from different host species. Conclusions: 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 s and other animals. Yaks, and blue were important zoonotic reservoirs of brucellosis causing species found in s. Keywords: Brucella, Molecular identification, Genotype, Zoonotic host Multilingual abstracts Please see Additional file for translations of the abstract into the six official working languages of the United Nations. Introduction Brucellosis is one of the most widespread and severe zoonotic diseases across the world. This disease affects mammals and is transmitted to s by animals * Correspondence: cuibuyun@icdc.cn; lingonghua@nwipb.cas.cn Equal contributors State Key Laboratory for Infectious Disease Prevention and Control, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing 00, China Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 000, Qinghai, China Full list of author information is available at the end of the article through direct contact with contaminated animal fluids or indirectly through consumption of unpasteurized milk products []. Brucellosis in domestic mammals and s is the cause of huge economic burden and medical challenges globally, especially in poor regions [ ]. The causal agent of brucellosis, Brucella genus, currently contains 0 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. to, B. abortus bv. to and, and B. suis bv. to []. Most species with a preferential host and three other species (B. abortus, B. melitensis, and B. suis) remain the principal causes of zoonotic potential []. 0 Ma et al. Open Access This article is distributed under the terms of the Creative Commons Attribution.0 International License (http://creativecommons.org/licenses/by/.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/.0/) applies to the data made available in this article, unless otherwise stated.

Ma et al. Infectious Diseases of Poverty (0) : Page of 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 [, ]. The Qinghai-Tibetan Plateau (QTP) is the largest and highest plateau on Earth, with an area of. 0 km and an average elevation of,000 m above sea level [0]. The QTP is an extensive pastoral and semi-pastoral area, and owing to poverty and bad hygiene, Brucella is highly endemic to this region among s (Bos grunniens) [, ] and s []. It should be mentioned that, probably due to technical difficulties, most previous studies on Brucella in the QTP merely examined infection surveillance, without paying adequate attention to the study of taxonomic identification of Brucella species. In the present study we collected a considerable number of Brucella samples from Qinghai Province, (Table ) which is located in the northeastern corner of the QTP and constitutes > % of the area of the plateau Table Geographic distribution of Brucella samples in Qinghai, China Locality ID County Sample ID Sample size Datong Ma0, Ma0 Dulan Ma, Ma, Ma, Ma, Ma, Ma, Ma, Ma, Ma Gangcha Ma Golmud Ma, Ma, Ma, Ma, Ma0, Ma, Ma,Ma, Ma, Ma, Ma, Ma, Ma Gonghe Ma0, Ma, Ma, Ma, Ma, Ma0, Ma, Ma Qilian Ma, Ma0, Ma, Ma, Ma, Ma, Ma, Ma Qumalai Ma Tianjun Ma0, Ma Tongde Ma, Ma, Ma0, Ma, Ma, Ma, Ma 0 Xinghai Ma0, Ma0, Ma0, Ma0, Ma0, Ma0, Ma Xining Ma0, Ma, Ma Zeku Ma Zhiduo Ma, Ma, Ma (Fig. ). We used both MLVA as well as MLST methods to type these samples and determine the composition characteristics of Brucella in the plateau area. Methods Ethics statement 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 isolates were examined: recovered from, from (Ovis aries), from blue (Pseudois nayaur), from s (Bos mutus grunniens), and each from cattle (Bos primigenius taurus) and Tibetan gazelle (Procapra picticaudata). Bacterial strains were cultured on Trypticase soy agar containing % blood (BD Diagnostic Systems, China Ltd., China) at C for h. Total genomic DNA was extracted using the DNeasy Blood & Tissue Kit (Qiagen, Germany) as per the manufacturer s recommended protocol. MLVA genotyping MLVA was performed as described earlier []. Briefly, primer pairs were designed and classified as belonging to one of the three panels: panel (bruce0, bruce0, bruce, bruce, bruce, bruce, bruce, and bruce), panel A (bruce, bruce, and bruce), and panel B (bruce0, bruce0, bruce0, bruce, and bruce0). PCR amplifications were performed in 0 μl reaction volumes, μl ofthesewere loaded in to % (panel ) or % (panels A and B) agarose gels containing 0. μg/ml of ethidium bromide, visualized under UV light, and photographed. Band intensities were estimated using BioNumerics version. (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 [] using B. melitensis M genome sequence as a reference (GenBank accession number NC_00. and NC_00.). 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 Brucella0 MLVA database (http://mlva.u-psud.fr/).

Ma et al. Infectious Diseases of Poverty (0) : Page of Fig. Geographic distribution of Brucella samples in Qinghai, China. QTP, Qinghai-Tibet Plateau; the dotted counties (NO. ~ ) correspond to the counties (Locality ID ~ ) in Table. The pie charts showed the composition (the relative pie size corresponds to the number of strains) of the four branches (A, B, C, and E) of B. melitensis MLST genotyping Nine distinct genomic loci, including seven housekeeping genes (gap, aroa, glk, dnak, gyrb, trpe, and cobq), one outer membrane protein (omp), and one intergenic fragment (int-hyp), were selected for MLST genotyping []. Similar to the MLVA process described, PCR amplifications were performed in 0 μ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 AM ~ AM0) 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.. []. 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.. []. 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://...0/brucellabase/mlst.html). Phylogenetic relationships of combined sequences were inferred using the Neighbor-Joining method in BrucellaBase. We could not identify biovars within strains, because it was difficult to distinguish among them using either MLVA or MLST markers [,, ]. Findings MLVA results Using panel markers, the present population clustered into eight known genotypes: (-------; N = ), (-------; N = ), (-------; N = ), (-------; N = ), (-------; N = ), (-------; N = ), (-------; N = ), and (---0----; N = ). The Clustering analysis showed that the isolates formed six main clusters (A ~ F). Cluster A had two genotypes ( and ); cluster B, C, E and F had a single genotype,, and, respectively; cluster D had three genotypes (,, and ) (Fig. ). According to Brucella0 MLVA database and based on panel markers we identified our samples as containing three species: B. melitensis (genotype, and ; cluster A, B, C and E), B. abortus (genotype, and ; cluster D) and B. suis (genotype ; cluster F). Based on a previous study [0], the B. melitensis strains were sub-grouped into East Mediterranean group (genotype and ) and American group

00 0 0 0 0 0 0 0 0 ID Ma Ma Ma Ma Ma Ma Ma Ma Ma Ma Ma Ma Ma Ma0 Ma Ma0 Ma0 Ma Ma0 Ma Ma Ma0 Ma0 Ma Ma0 Ma0 Ma Ma Ma Ma Ma Ma Ma0 Ma Ma Ma Ma Ma Ma Ma Ma Ma0 Ma Ma Ma Ma0 Ma Ma Ma0 Ma Ma Ma Ma Ma Ma Ma0 Ma Ma Ma0 Ma Ma Ma0 Ma Ma Ma bruce0 bruce0 bruce bruce 0 bruce bruce bruce bruce bruce bruce 0 bruce bruce0 bruce0 bruce0 bruce 0 0 bruce0 panel Host blue Tibetan gazelle blue blue blue blue blue blue blue catle blue blue blue Isolate 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 00 00 Year 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 A B C D E F MLVA- Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm0 Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm0 Bm0 Bm0 Bm0 Bm0 Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm Bm0 Bm0 Bm Bm Ba0 Ba0 Ba0 Ba0 Bm Bm Bs0 B. m. B. m. B. m. B. s. B. m. B. a. Fig. Dendrogram based on the MLVA genotyping assay showing relationships of the Brucella isolates. ID: serial number for the isolates; MLVA type: panel and MLVA- genotypes; Host: the hosts from which the bacteria was isolated; Isolate: isolate name in the laboratory in which the DNA extraction was done; Year: the years when the strains were collected Ma et al. Infectious Diseases of Poverty (0) : Page of

Ma et al. Infectious Diseases of Poverty (0) : Page of (genotype ). Considering all three MLVA panels, (Bm0 ~ Bm), (Ba0 ~ Ba0), and (Bs0) MLVA- genotypes were identified for B. melitensis, B. abortus and B. suis isolates respectively (Fig. ). No genotypes completely matched any of the genotypes in the Brucella0 MLVA database. MLST results A total of five known MLST genotypes were identified: ST (-------0-; N = ), ST (--------; N = ), ST (--------; N = ), ST (------ --; N = ), ST (--------; N = ). The Neighbor-Joining clustering showed that the isolates formed four main clusters a, d, e and f. Of these, clusters a, e and f correspond to genotypes ST, ST and ST respectively while cluster d corresponds to two genotypes, ST and ST. According to Whatmore et al. [], genotypes ST and ST belong to B. melitensis, ST and ST belong to B. abortus, and ST belongs to B. suis. The gyrb sequence of strain Ma0 was different from any of known alleles (allele ~ ). Its sequence closely matched that of allele (identical except for a dinucleotide change of GC to AA on site 0 ~ 0 relative to the gyrb allele sequence or on 0 ~ 0 sites relative to the ST sequence). Considering all MLST alleles, the strain Ma0 could be viewed as a mutant variant (----?----) of ST (--------). Discussion In this study we used both MLVA and MLST methods to identify Brucella species in the east of the Qinghai- Tibet Plateau. Although there were some minor incongruences for e.g., in clustering results (Fig. and Fig. ), both methods showed a consistent conclusion that there were at least three Brucella species in the sample analyzed. Of the isolates 0,, and were respectively identified to B. melitensis, B. abortus, and B. suis, indicating that B. melitensis was a predominant species on the plateau. Interestingly however, when searching in the Brucella0 MLVA database none of the genotypes found in our study were identical to any of the genotypes in the database. Based on MLST the ST seems to be the main ST type in the QTP region. Moreover, there was also a unique genotype (strain Ma0) which has not been reported anywhere before. These results demonstrated that the Brucella population found in the QTP region was very different from that in other regions. We attribute this to relatively isolated and special environment of the plateau. Extremely high altitudes make it difficult for lowland livestock breeds or wild animals to survive, thereby favoring only the endemic ones. Hence, there was a very limited breeds exchange between QTP and other regions, which may have consequently formed a unique local Brucella population. Based on MLVA-, ten genotypes were shared by two or more host species. To elaborate further, four genotypes (Bm0, Bm, Bm0, and Bm) were shared between s and, and two (Bm0 and Ba0) were shared between s and s (Fig. ). Based on MLST, the ST type was also shared between s and four other animals (, s, blue and Tibetan gazelle); the ST type was shared between s and s (Fig. ). Sheep and s are the main livestock on QTP and it is therefore not surprising that these animals pose a heavy zoonotic disease threat to s. It should be noted that, all and four out of six s 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 M, S and S, which are used to immunize animals in China for B. melitensis, B. abortus, and B. suis, respectively []. Because s could be infected by either B. melitensis or B. abortus, we suggest based on our study that both the B. melitensis type (M) and B. abortus type vaccines (S) should be tried in this animal. The blue were infected by all three Brucella species including B. suis, which is generally carried by pigs [, ]. More importantly, this animal shared Brucella genotypes with (Bm0, Bm and Bm), s (Bm0), and Tibetan gazelle (Bm0), indicating this animal is an important natural repertory for Brucella. Since there is obvious spatial and diet overlap between blue and domestic livestock [], we suggest that this animal is also an important infection source of Brucella in s. Additionally, due to high genetic diversity of B. melitensis, we also discuss the geographic distribution features of the strains. Figure showed that counties (locality ID:,,,,, 0, and ) were dominated by A branch strains; counties (locality ID:,, and ) were dominated by B branch strains; counties (locality ID:,, and ) contained only C branch strains; and the E branch genotypes appeared only in county (locality ID: ). 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.

Ma et al. Infectious Diseases of Poverty (0) : Page of 0.000 d B.a. a B.m. e B.m. f B.s. gap aroa glk dnak gyrb trpe cobq omp int-hyp ST ID Host M a 0 s h e e p M a 0 s h e e p M a 0 h u m a n M a 0 h u m a n M a 0 h u m a n M a 0 h u m a n M a 0 h u m a n M a s h e e p M a s h e e p M a h u m a n M a h u m a n M a h u m a n M a h u m a n M a h u m a n M a h u m a n Ma M a 0 h u m a n M a h u m a n M a h u m a n M a h u m a n M a h u m a n M a h u m a n M a h u m a n Ma M a h u m a n M a h u m a n M a 0 h u m a n M a h u m a n M a h u m a n M a h u m a n M a h u m a n Ma M a h u m a n M a b l u e s h e e p Ma Tibetan gazelle M a s h e e p Ma0 Ma M a s h e e p M a s h e e p M a s h e e p Ma Ma Ma M a s h e e p M a s h e e p M a 0 y a k M a y a k M a y a k M a y a k M a b l u e s h e e p M a b l u e s h e e p M a b l u e s h e e p M a 0 b l u e s h e e p M a b l u e s h e e p M a b l u e s h e e p M a b l u e s h e e p M a 0?? h u m a n M a 0 b l u e s h e e p Ma 0 blue M a 0 h u m a n M a y a k M a 0 y a k M a c a t t l e M a b l u e s h e e p Fig. Dendrogram based on the MLST genotyping assay showing relationships of the Brucella isolates. ID: serial number for the isolates; ST: MLST genotypes; Host: the hosts from which the bacteria was isolated. (Also see Fig. for additional information)

Ma et al. Infectious Diseases of Poverty (0) : Page of 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 s and animals, brucellosis surveillance should be regularly executed on both livestock ( and s) and wild ungulates (blue ). 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 Brucella0 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. Conclusions 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 s and animals;, s and blue were important zoonotic hosts of brucellosis to s in the area studied. Additional file Additional file : Multilingual abstracts in the six official working languages of the United Nations. (PDF 0 kb) Abbreviations MLST: multi locus sequence typing; MLVA: multi-locus variable-number tandem-repeat analysis; QTP: Qinghai-Tibet Plateau. Competing interests The authors declare that they have no competing interests. Authors contributions JYM, HW, BYC and GHL contributed to the design of the study. JYM, BYC and GHL collected the data and wrote the paper. JYM, XFZ, LQX, GYH, HJ, FZ, HYZ, DRP, YMQ and GHL analyzed the data. All authors read and approved the final paper. Acknowledgments This study was supported by the Science and Technology Plan of Qinghai Province (NO. 0-Z-) and the Youth Innovation Promotion Association of Chinese Academy of Sciences (NO. 0). Author details Qinghai Institute for Endemic Disease Prevention and Control, Xining 0, Qinghai, China. State Key Laboratory for Infectious Disease Prevention and Control, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing 00, China. Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 000, Qinghai, China. Received: October 0 Accepted: March 0 Published: April 0 References. Atluri VL, Xavier MN, de Jong MF, den Hartigh AB, Tsolis RM. Interactions of the pathogenic Brucella species with their hosts. Annu Rev Microbiol. 0;:.. Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos EV. The new global map of Brucellosis. Lancet Infect Dis. 00;:.. Whatmore AM, Perrett LL, MacMillan AP. Characterisation of the genetic diversity of Brucella by multilocus sequencing. BMC Microbiol. 00;:.. Kansiime C, Rutebemberwa E, Asiimwe BB, Makumbi F, Bazira J, Mugisha A. Annual trends of brucellosis in pastoralist communities of southwestern Uganda: a retrospective ten-year study. Infec Dis Poverty. 0;:.. Ducrotoy MJ, Ammary K, Lbacha HA, Zouagui Z, Mick V, Prevost L, Bryssinckx W, Welburn SC, Benkirane A. Narrative overview of animal and brucellosis in Morocco: intensification of livestock production as a driver for emergence? Infec Dis Poverty. 0;:.. Garin-Bastuji B, Mick V, Le Carrou G, Allix S, Perrett LL, Dawson CE, et al. Examination of taxonomic uncertainties surrounding Brucella abortus bv. by phenotypic and molecular approaches. Appl Environ Microbiol. 0; 0():0.. 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. 00;:.. 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. 0;():.. 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. 0;():e. 0. Wu T. The Qinghai-Tibetan plateau: how high do Tibetans live? High Alt Med Biol. 00;:.. Wiener G, Han J, Long R. The. nd ed. Bangkok: RAP publication; 00.. Lang X, Qu H, Bu Z, Yang Y, Sun C, Li X, et al. Seroprevalence of Brucella infection in s (Bos grunniens) on the Qinghai-Tibet plateau of China. Trop Anim Health Prod. 0;:0.. Xue H. Analysis on infection and the related factor of brucellosis in high-risk population in Qinghai Province of 00 to 00. Taiyuan, Shanxi Province, China: Master thesis of Shanxi Medical University; 0.. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 0;(0):.. 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. 0;:.. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 00;0:.. Marco C, Roberta B. BLAST Parser. BITS Conference. Societa di Bioinformatica Italiana. Bologna. Societa di Bioinformatica Italiana Bologna Italy. 00;:.. Jiang H, Fan M, Chen J, Mi J, Yu R, Zhao H, et al. MLVA genotyping of Chinese Brucella melitensis biovar, and isolates. BMC Microbiol. 0;:.. 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 s and identification of Brucella suis vaccine strain s from cattle in China. PLoS One. 0;(0):e.

Ma et al. Infectious Diseases of Poverty (0) : Page of 0. Al Dahouk S, Le Flèche P, Nöckler P, Jacques I, Grayon M, Scholz HC, et al. Evaluation of Brucella MLVA typing for brucellosis. J Microbiol Methods. 00;:.. Shang D, Xiao D, Yin J. Epidemiology and control of brucellosis in China. Vet Microbiol. 00;0:.. 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. 0;:.. 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. 00;():. doi:0./s0-0-0-z Cite this article as: Ma et al.: MLVA and MLST typing of Brucella from Qinghai, China. Infectious Diseases of Poverty 0 :. Submit your next manuscript to BioMed Central and we will help you at every step: We accept pre-submission inquiries Our selector tool helps you to find the most relevant journal We provide round the clock customer support Convenient online submission Thorough peer review Inclusion in PubMed and all major indexing services Maximum visibility for your research Submit your manuscript at www.biomedcentral.com/submit