Skip to main content
Download PDF

Analysis of the Brucella melitensis epidemic in Xinjiang: genotyping, polymorphism, antibiotic resistance and tracing

Annals of Clinical Microbiology and Antimicrobials volume 23, Article number: 71 (2024) Cite this article

Abstract

Brucella spp. are facultative intracellular pathogens that cause zoonosis- brucellosis worldwide. There has been a trend of the re-emergence of brucellosis worldwide in recent years. The epidemic situation of brucellosis is serious in Xinjiang. To analyze the epidemic situation of Brucella spp. in Xinjiang among humans and animals, this study identified 144 Brucella isolates from Xinjiang using classical identification and 16 S rRNA sequencing. MLVA, drug resistance testing, and wgSNP detection were also performed. At the same time, analysis was conducted based on the published data of Brucella isolates worldwide. The results showed that the dominant species was B. melitensis biovar 3, which belonged to GT42 (MLVA-8 typing) and the East Mediterranean lineage. The correlation among isolates was high both in humans or animals. The isolates in Xinjiang exhibited higher polymorphism compared to other locations in China, with polymorphism increasing each year since 2010. No amikacin/kanamycin-resistant strains were detected, but six rifampicin-intermediate isolates were identified without rpoB gene variation. The NJ tree of the wgSNP results indicated that there were three main complexes of the B. melitensis epidemic in Xinjiang. Based on the results of this study, the prevention and control of brucellosis in Xinjiang should focus on B. melitensis, particularly strains belonging to B. melitensis bv.3 GT42 (MLVA-8 typing) and East Mediterranean lineage. Additionally, the rifampicin- and trimethoprim-sulfamethoxazole- resistance of isolates in Xinjiang should be closely monitored to avoid compromising the therapeutic efficacy and causing greater losses. These results provide essential data for the prevention and control of brucellosis in Xinjiang and China. Although the isolates from Xinjiang have significant characteristics among Chinese isolates and can reflect the epidemiological situation of brucellosis in China to some extent, this study cannot represent the characteristics of isolates from other regions.

Introduction

Brucella spp. are facultative intracellular pathogens that cause the zoonosis- brucellosis. Brucellosis not only affects public health and safety but also influences animal husbandry and economic development, leading fever, fatigue and reduced labor in humans, and reproductive disorders in animals [1]. Brucellosis is epidemic in South America, the Middle East, and parts of Asia [2], with more than half a million cases annually and incidence rates exceeding 100 per million in endemic countries [3]. However, the actual incidence is estimated to be 10~25 times higher than reported, according to the data from World Health Organization (WHO) [4]. Six classical species and six new (non-classical) species have been identified based on biochemical characteristics and host preferences [5]. Among them, B. melitensis, B. abortus and B. suis are the primary causes of human brucellosis [5]. In China, existing results showed that B. melitensis biovar 1 was the dominant strain in previous years, but B. melitensis biovar 3 was the main isolate in recent years [6, 7].

The spread of Brucella spp. is facilitated by livestock trade and tourism, leading to a pattern of re-emergence and constant expansion of brucellosis [4, 8]. According to the data collected by WOAH in 2021 (https://wahis.woah.org), France, Cyprus, Croatia, Israel, Australia, and other countries reported brucellosis outbreaks from 2019 to 2022. The last reported brucellosis outbreaks in these countries occurred 10 years ago or even longer. There also has been a trend of the re-emergence of brucellosis in China in the last decade [9]. According to the data released by the Chinese Center for Disease Control and Prevention (CCDC), the numbers of newly reported brucellosis cases in China from 2020 to 2022 were 48,455, 69,767 and 65,384 (https://www.phsciencedata.cn/Share/). Previous studies showed that the main epidemic strains of B. melitensis in China belonged to 9 complexes, of which one isolate from Xinjiang Uygur Autonomous Region (referred to as Xinjiang in this study) was associated with all complexes [10], and an amikacin-resistant B. abortus strain was first isolated from Xinjiang [11]. The location of Xinjiang is very important. Xinjiang is located in northwestern China and in the hinterland of the Eurasian continent. Xinjiang borders the provinces of Qinghai, Tibet and Gansu in China, and borders Mongolia, Russia, Kazakhstan, Afghanistan, Pakistan and India. Brucellosis prevalence in these areas is relatively serious. These factors haveurged us to consider the brucellosis epidemic trend in Xinjiang.

In this study, the prevalence of Brucella spp. in Xinjiang was comprehensively analyzed, including genotyping, polymorphism, antibiotic resistance and tracing. Identification typing and polymorphism of Brucella spp. were based on conventional and molecular methods. Compared to the conventional methods (such as culturing, agglutination testing, phage typing, and biochemical analyses), molecular methods (mainly based on multiple PCR molecular typing) were shown to be rapid, relevant, and efficient for typing and clustering Brucella strains [12]. Multiple Locus Variable Number of Tandem Repeats Analysis (MLVA) was wildly used in typing and polymorphism of Brucella spp [13]. In recent years, SNP (single nucleotide polymorphism), especially wgSNP (whole genomic SNP), has been widely used for tracing of Brucella isolates [14, 15]. The aims of this study were as follows: (1) MLVA was used to analyze the genotype, polymorphism, and individual differences of Brucella isolates in Xinjiang to better understand the changes in the population; (2) amikacin/kanamycin-resistance detection was used to analyze whether the amikacin/kanamycin -resistant isolates epidemic were in ruminants and humans; (3) the minimum inhibitory concentration (MIC) values of the clinic brucellosis treatment drug for isolates were detected to analyze resistance; and (4) a phylogenetic tree was constructed with the wgSNP to analyze the relationship between the isolates in Xinjiang and those Brucella isolates in China and to trace these strains.

Materials and methods

The strains and genomes of Brucella spp. used in this study

The genome sequences of the public B. melitensis with information of collection date and location were downloaded from NCBI GenBank (www.ncbi.nlm.nih.gov/genome/?term=brucella) (Table S 1). The wgSNP of Chinese Brucella isolates involved in this study had been published previously [10]. Cultured and inactivated Brucella isolations used in this study were carried out in a BSL-3 laboratory of the National Reference Laboratory for Animal Brucellosis. The information of all the isolates used in this study was listed in Table S1.

Identification

Classical identification of Brucella strains may be performed by using a combination of the following detections [16]: CO2 requirement test, production of H2S (detected by lead acetate papers) test, urease and oxidase tests, the slide agglutination test with anti-A, -M or -R monospecific sera, and/ or phage Tbilisi (Tb) lysis test. In this study, positive control using B. melitensis bv.1 strain 16 M and negative control using saline water were established for each test.

MLVA typing

After the Brucella isolates from Xinjiang were identified as Brucella species by 16 S rRNA sequencing in our laboratory, the biochemical characteristic detection was used to identify the biovar. The genomic DNA of these isolates was extracted using the Wizard Genomic DNA Purification kit (Promega, USA) according to the manufacturer’s instructions. The primers for MLVA-16 typing were referred to the published study [17], and all primers used in this study were listed in Table S1. A phylogenetic tree using the UPGMA method was generated with BioNumerics software (http://www.applied-maths.com/bionumerics) using default software settings [6]. Meanwhile, a minimum spanning tree (MST) was also generated with BioNumerics software [10, 18].

This study also downloaded and collected published MLVA typing data from MLVA bank databases (http://mlva.i2bc.paris-saclay.fr/brucella/) and published studies [19,20,21,22,23,24,25,26,27] (Table S2). If multiple data points were available from the same province/site/country, the data for MLVA-16 analysis were chosen in this study.

Polymorphism analysis

The polymorphism analysis was based on the results of MLVA typing. All the results of MLVA sites were uploaded to the Comparing Partitions Website (http://www.comparingpartitions.info/index.php?link=Tut4), and polymorphism analysis was performed.

Amikacin/kanamycin-resistance detection

Based on a previous study [11], our lab established a real time PCR (qPCR) method to detect the aph(3’)-IIa gene. The upstream primer sequence was 5’- AACTGTTCGCCAGGCTCAAG − 3’, and the downstream primer sequence was 5’- AAGCGGCCATTTTCCACCAT − 3’. The probe sequence was 5’-VIC- TCGTGACCCATGGCGATGCCTGCTT-BHQ-3’. The specific primers for Brucella species were described in a published study [28]. The amplification system was 20 µL / sample, including 10 µL of premix ex Taq (probe qPCR) (Takara, Japan), 5.6 µL of ddH2O, 0.4 µL of primers and probes (4 primers and 2 probes, the final concentration was 0.2 µM), and 2 µL of samples. For two-step amplification detection, the conditions were pre denaturation at 95 °C for 30 s, reaction for 45 cycles, 95 °C for 5 s, 60 °C for 30 s, and the corresponding fluorescence signal was collected at the end of each reaction. This method used double probes to detect amikacin/kanamycin-resistant Brucella strains with good specificity. The minimum Brucella content detected was 102 CFU/µL. All the primers and probes were synthesized by Sangon Biotech (Shanghai) Co., Ltd (Shanghai, China).

MIC value

According to the recommendation of the WHO and the clinical medication for brucellosis in China, rifampicin (REP), doxycycline (DOX), streptomycin (SM), trimethoprim-sulfamethoxazole (TMP/SMZ, trimethoprim: sulfamethoxazole = 1:19), ceftriaxone sodium (CRO), gentamicin (GM), levofloxacin (LfX), ampicillin (AMP) and ciprofloxacin (CIP) were selected to detect the MIC value. The susceptibility of isolates to these drugs (National Institutes for Food and Drug Control, China) was tested by the double dilution method in 96-well plates according to CLSI guidelines. Isolates were cultured on Brucella agar (BD, USA) media and incubated at 35 ± 2 ℃ for 72~96 h. Isolates were suspended in saline water to 0.5 McF turbidity and suspended in Brucella broth (BD, USA). Drugs were diluted in a 96-well plate ranging from 0.125 to 256 µg/mL. Then 100 µL dilute solution (approximately 1 × 105 CFU/well) was added. The plate was cultured at 35 ± 2 °C for 48 h. The quality control strain was Streptococcus pneumoniae ATCC 49619 (refer to CISL_M45 (2016)) (to test whether the drug dilution is qualified), and the positive control strain was B. melitensis bv.3 strain Ether.

Whole genome sequencing and SNPs identification

Library preparation of isolates was performed using the MGIEasy Universal DNA Library Kit (MGI, China) according to the manufacturer’s manual. Genomic DNA was sheared randomly to construct three read-libraries with lengths of 350 bp physico-chemical methods by Beijing Genomics Institute (BGI) (Shenzhen, China). The genomes were sequenced above 100-fold (100×) genome coverage. SNPs and indels identification and annotation were performed as described in a previous study [10]. Some SNPs (not all) was identified by PCR amplification and sequencing.

Clustering analyses and phylogenetic tree construction

The genome of B. melitensis bv.2 str. ATCC 23457 (GCA0000_22625.1) was selected as the reference. SNP calling and annotation were performed as described in a previous study [29]. PhyloSNP [30] was used to generate phylogenetic trees. The neighbor-joining (NJ) method with 1000 bootstraps was used.

Statistical analysis

Statistical analyses were performed using Excel and SPSS. A P value < 0.05 was considered significant when using one-way analysis of variance (ANOVA). The figure of MST was generated with BioNumerics software, and the figure of polymorphism of isolates in China was generated with ArcMap (version 10.2). The figure of the NJ tree based on wgSNP was exported by iTOL (https://itol.embl.de). Other figures were drawn by Excel software.

Results

All isolates belonged to B. melitensis biovar 3 in both humans and animals

In this study, 144 Brucella strains isolated from cattle and sheep in 2017~2019 were obtained. Among them, 43 strains were isolated from Tarbagatay. The numbers of isolates from Changji, Altay, Ili, Aksu, Bayingol, Kashgar and Tulufan were 16, 27, 32, 16, 6, 2 and 2, respectively. Additionally, 30 Brucella isolates were from cattle, and the rest were from sheep.

The results of classical identification showed that all the isolates were B. melitensis biovar 3. Compared with published data of isolates from human in Xinjiang [18], the dominant species was B. melitensis biovar 3 both in humans and animals.

The isolates exhibited high correlation in Xinjiang by MLVA typing

MLVA-8 typing. The results of MLVA were uploaded to the MLVA bank for MLVA-8 typing. The results showed that the isolates were divided into 6 genotypes (GT), which included two new genotypes named NEW1 and NEW2 (Fig. 1A). The MST showed that GT 42 was centrally located and related to other isolates. GT 63, GT 45, NEW1 and NEW2 had only one VNTR site difference when compared to GT42. There was also one VNTR site difference between GT 112 and GT 63. GT 45, NEW1 and NEW2 may vary from GT 42.

MLVA-16 typing. The results showed that the isolates were divided into 58 genotypes, including 40 genotypes that contained only one isolate. All the genotypes belonged to the East Mediterranean lineage (Fig. 1B).

Fig. 1
figure 1

MST based on the data of MLVA typing.  A was the MST based on the data of MLVA-8 typing. B was the MST based on the data of MLVA-16 typing. Different colors represented different locations. The area of the circle indicated the number of strains. And the larger the area of the circle, the more strains of the genotype. Numbers and/or letters outside the circle represented genotypes. For example, in Fig. 1A, the number outside the largest circle was 42, indicating that GT 42 had the most strains, and different colors in the circle mean the locations and number of strains that composed this genotype. To distinguish from the results of MLVA-8 typing, the pre-genotype of MLVA-16 (Fig. 1B) was marked as “M”

Full size image

The polymorphism of isolates in Xinjiang was higher compared to other locations increased year by year

The Simpson index showed that 16 VNTR sites had different degrees of polymorphism. Only two VNTR sites exhibited high polymorphism. Compared with Panel 1 (Bruce06, Bruce08, Bruce11, Bruce12, Bruce42, Bruce43, Bruce45 and Bruce55) and Panel 2 A (Bruce18, Bruce19 and Bruce21), the polymorphism of Panle 2B (Bruce04, Bruce07, Bruce09, Bruce16, Bruce30) was higher. Among them, the polymorphism of Bruce16 was highest, and the Simpson index was 0.826. The isolates from Xinjiang in this study exhibited high genetic relatedness (Table 1). Statistical polymorphism by location showed that the isolates from Altay had low polymorphism, while those from Aksu had high polymorphism.

Table 1 Simpson index of VNTR sites in Xinjiang
Full size table

Published data about Brucella isolates from Xinjiang (323 strains, including the isolates in this study) were also analyzed for polymorphism. The results showed that the VNTR sites of Bruce21, Bruce55, Bruce45 and Bruce08 had no polymorphism in isolates from Xinjiang. Compared with the polymorphism results of isolates from Xinjiang among different years, the polymorphism of Brucella isolates in Xinjiang increased year by year.

This study also analyzed polymorphisms in published data from Chinese Brucella isolates based on the results of MLVA-16 typing (1754 strains, including the isolates in this study). The results showed that Inner Mongolia, Xinjiang, Qinghai, Ningxia and Hainan had higher polymorphism, while Guangdong had lower polymorphism among the locations with more than 50 isolates (Fig. 2).

Fig. 2
figure 2

The polymorphism of isolates in China based on the data of VNTR sites.  Different colors represented different VNTR sites. Each circle was composed of different colors. The more colors indicated the higher polymorphism of isolates in the location, and the less colors indicated the lower polymorphism of isolates. The larger the color area in the circle, the higher polymorphism of these VNTR sites

Full size image

Six rifampicin intermediate Brucella isolates were identified

All the isolates were tested for amikacin-resistant with qPCR based on the methods established by our laboratory, and no amikacin-resistant isolates were detected.

The drugs used in this study were double diluted from 0.125 to 256 µg/mL in 96-well plates (except for TMP/SMZ). The MIC values of B. melitensis bv. 3 strain Ether was used as a reference. The MIC value of isolates of REP, DOX, SM, TMP/SMZ, CRO, GM, LfX, AMP and CIP were 1~2 µg/mL, 0.06~0.12 µg/mL, 0.5~4 µg/mL, 0.03/0.06~0.24/4.8 µg/mL, 0.12~2 µg/mL, 0.015~0.06 µg/mL, 0.25~2 µg/mL, 0.5~4 µg/mL and 0.5~1 µg/mL, respectively. The MIC value of reference strain for these drugs were 0.5 µg/mL, 0.12 µg/mL, 2 µg/mL, 0.015/0.3 µg/mL, 0.5 µg/mL, 0.03 µg/mL, 0.5 µg/mL, 1 µg/mL and 0.5 µg/mL, respectively. Brucella isolates were sensitive to GM, SM, DOX and TMP-SMZ based on CLSI M45. According to the standard of other slow-growing bacteria of CLSI for rifampin, any MIC value ≤ 1 µg/mL was susceptible, 2 µg/mL was intermediate, and ≥ 4 µg/mL was resistant. Six rifampicin-intermediate strains were found in this study. Therefore, the highly variable sites of the rpoB gene were detected according to a published study [31], and no rpoB gene mutation was found (Fig. 3).

Fig. 3
figure 3

Number of isolations distribution of fold change of MIC values for different drugs (Compared with reference strain of B. melitensis bv.3 strain Ether).  Different colors represented different drugs. The area of the circle indicated the number of strains. For REF, the MIC value of reference was 0.5 µg/mL, the MIC value of isolates was 1 µg/mL, then the fold change of MIC value was 2

Full size image

Three complex Brucella isolates were epidemic in Xinjiang

Based on the results of MLVA-16 typing, five regions with more than 5 isolates were selected for further analysis. Two isolates were randomly selected from each region. The selected strains were N2-17, N2-55, N8-7, N8-121, N9-1, N9-21, N11-12, N11-65, Y13-103 and Y13-113. The wgSNPs of 10 sequenced strains in this study, B. melitensis isolates around the world and standard strains of different species were used to construct the phylogenetic tree by the NJ method (Figure S1). The results showed that the eight strains sequenced in this study were clustered with the Chinese isolates, and the two strains isolated from Tarbagatay, were in the same branch as the strains isolated from Turkey, Russia, India, etc.

The 10 isolates were combined with the published Chinese isolates and SNPs involved in previous studies to construct a genome-wide NJ phylogenetic tree (Fig. 4). The results showed that the two isolates from Tarbagatay, which were closely related to foreign isolates, were not in the same branch as those from other regions in China, indicating that the isolate was only prevalent in Xinjiang. The strains isolated from Altay, Changji and Ili were in the same branch. The strains of this branch were all isolated from Xinjiang and were not in the same branch as the strains isolated from other regions. The other isolates from Bayingol and Changji were in the same branch as those from Guangdong. Another isolate from the Ili region was in the same branch as the isolates from Inner Mongolia and Gansu and was also the closest isolate to the isolates from other regions in China. This was also the main epidemic branch in China at present.

Fig. 4
figure 4

NJ tree based on wgSNP of Brucella isolates.  The red letters were the isolates sequenced in this study.The colored ranges were the complex including sequenced isolates. The length of the tree scale represented one SNP constructed the phylogenetic tree

Full size image

Discussion

Based on the results of this study, the prevention and control of brucellosis in Xinjiang should focus on B. melitensis, particularly strains belonging to B. melitensis bv.3 GT42 (MLVA-8 typing) and East Mediterranean lineage. The results showed that there were several genotypes of Brucella isolates in China; however, GT42 was the most prevalent (MLVA-8) [7, 23]. The results of MST based on MLVA typing showed that the isolates from China were highly correlated with those from Kazakhstan and Turkey [32]. MLVA typing could be used to distinguish species, biovar or some isolates and analyze the relationships between the isolates [33]. However, it could not be used for tracing due to the risk of homology of these repeated sequences [12]. The MST results constructed by MLVA typing in this study were similar to the published studies in China [32, 34], and the isolates in Xinjiang exhibited a close relationship with those in other regions in China.

The isolates from regions that are major points of brucellosis spread in China need more attention due to their high polymorphism. The polymorphism index of isolates showed the differences in polymorphism in different regions. The higher the Simpson index, the greater the polymorphism. It was found that strains from Inner Mongolia, Xinjiang, Qinghai, Ningxia and Gansu with high polymorphism were developed in animal husbandry and endemic regions of brucellosis. Strains from Hainan, Fujian and Zhejiang, which had low polymorphism, belonged to cattle and sheep product or dairy product inflow regions (https://data.stats.gov.cn/index.htm). All of these regions belonged to the major regions or important points of the spread of brucellosis in China [6]. The regions with high polymorphism exhibited higher genetic diversity and genomic variation of Brucella spp. Therefore, there was a possibility of changes or alterations in the dominant epidemic Brucella strains. The isolates of these regions should be taken seriously, and the frequency of surveillance needs to be increased to avoid the emergence of new variations in these regions. Further research could focus on the correlation between high polymorphism and specific factors to better understand the evolution of B. melitensis.

The rifampicin- and trimethoprim-sulfamethoxazole- resistance of isolates in Xinjiang should be closely monitored to avoid compromising the therapeutic efficacy and causing greater losses. This study found no transmission of amikacin/kanamycin-resistant strains was found in cattle and sheep in this study. The results of drug resistance analysis indicated that B. melitensis in Xinjiang was sensitive to common brucellosis treatment drugs, but there was a trend of resistance to rifampicin and trimethoprim-sulfamethoxazole (with higher MIC values). In recent years, drug resistance detection of Brucella spp. has been carried out worldwide, and some resistant strains have been found. Current studies have shown that rifampicin-resistant Brucella strains were isolated in Turkey [35], Egypt [36], Norway [37], Brazil [38], Kazakhstan [39], Qatar [40], Malaysia [41], South Korea [42], China [43], Iran [44], Italy [45] and other countries. Trimethoprim-resistant strains were detected in Brazil [38], Turkey [46], China [43] and Iran [44]. In addition, ciprofloxacin- and streptomycin-resistant strains were observed in Brazil [38], Turkey [46] and South Korea [42]. Studies of drug resistance for Brucella spp. were useful to analyze the variation and epidemic trends of Brucella under antibiotic selection pressure. At the same time, those results could guide clinical medication and promote the prevention, control, and eradication of brucellosis.

The results of genomic analysis show that isolates in China were in the same branch but not in the same branch as those in other regions of the world [15]. Chinese isolates were closely related to those in the Middle East, such as in Turkey. A study used the published genome of Brucella to construct an NJ tree based on wgSNP and showed that the isolates named BL from Xinjiang were in the same branch as the isolates from Turkey, Russia and other regions [19]. The phylogenetic tree based on wgSNP of two Tarbagatay isolates obtained in this study showed that these were in the same branch as the isolates from Turkey, and this branch also included the isolates from Russia, India and other regions. However, these isolates were not in the same branch as the isolates from other regions of China.

The results of the phylogenetic tree constructed in this study based on wgSNP indicated that there were mainly three complexes of the B. melitensis epidemic in Xinjiang. The first was the Xinjiang-specific complex, mainly epidemic in different locations in Xinjiang. The second complex was mainly related to or similar to the strains in most regions of China. The third complex was mainly similar to the isolates from Turkey and other regions. The epidemic area of this third complex in Xinjiang was limited, and it was not epidemic in other regions in China. Due to the enforcement of brucellosis detection in cross-border trade, it was speculated that foreign strains were transmitted through wild animals. A study which conducted serum detection for 258 wild ruminants in Xinjiang found that the Brucella spp. serum positive rate was 2.3% (95 confidence interval 0.5~4.2%) [47], which to some extent supports the results of this study. Meanwhile, the transmission of B. melitensis also deserves further investigation.

Conclusion

In conclusion, the dominant Brucella species exhibits higher polymorphism in Xinjiang than other locations in China, and there is an emerging trend of rifampicin-resistant trend in Xinjiang. This study identified 144 isolates of Brucella and their biovars. MLVA typing, drug resistance testing, and wgSNP detection were also performed. Additionally, analysis based on the published data of Brucella isolates worldwide was conducted. The results showed that the dominant species was B. melitensis biovar 3, and the correlation among isolates was high. Since 2010, the polymorphism of isolates in Xinjiang has been increasing each year. No reported amikacin/kanamycin-resistant strains were found, but 6 rifampicin-intermediate isolates were found without rpoB gene variation. The NJ tree of the wgSNP results indicated that there were three main complexes of the B. melitensis epidemic in Xinjiang. This study analyzes the prevalence, genotyping, diversity, drug resistance, and tracing of B. melitensis in Xinjiang. Based on the results of this study, the prevention and control of brucellosis in Xinjiang should focus on B. melitensis in both humans and animals, particularly strains belonging to B. melitensis bv.3 GT42 (by MLVA-8 typing) and East Mediterranean lineage. Additionally, the rifampicin- and trimethoprim-sulfamethoxazole- resistance of isolates in Xinjiang should be closely monitored to avoid compromising the therapeutic efficacy and causing greater losses.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

References

  1. Nelson-Jones A. Brucellosis. Postgrad Med J. 1952;28(324):529–34.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Liu Z, Gao L, Wang M, Yuan M, Li Z. Long ignored but making a comeback: a worldwide epidemiological evolution of human brucellosis. Emerg Microbes Infect. 2024;13(1):2290839.

    Article  PubMed  Google Scholar 

  3. Bukhari EE. Pediatric brucellosis. An update review for the new millennium. Saudi Med J. 2018;39(4):336–41.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos EV. The new global map of human brucellosis. Lancet Infect Dis. 2006;6(2):91–9.

    Article  PubMed  Google Scholar 

  5. Kurmanov B, Zincke D, Su W, Hadfield TL, Aikimbayev A, Karibayev T, Berdikulov M, Orynbayev M, Nikolich MP, Blackburn JK. Assays for identification and differentiation of Brucella species: a review. Microorganisms. 2022;10(8):1584.

  6. Tian GZ, Cui BY, Piao DR, Zhao HY, Li LY, Liu X, Xiao P, Zhao ZZ, Xu LQ, Jiang H, et al. Multi-locus variable-number tandem repeat analysis of Chinese Brucella strains isolated from 1953 to 2013. Infect Dis Poverty. 2017;6(1):89.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Jiang H, Fan M, Chen J, Mi J, Yu R, Zhao H, Piao D, Ke C, Deng X, Tian G, et al. MLVA genotyping of Chinese human Brucella melitensis biovar 1, 2 and 3 isolates. BMC Microbiol. 2011;11:256.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Rohr JR, Barrett CB, Civitello DJ, Craft ME, Delius B, DeLeo GA, Hudson PJ, Jouanard N, Nguyen KH, Ostfeld RS, et al. Emerging human infectious diseases and the links to global food production. Nat Sustain. 2019;2(6):445–56.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Jiang H, O’Callaghan D, Ding JB. Brucellosis in China: history, progress and challenge. Infect Dis Poverty. 2020;9(1):55.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Yang X, Piao D, Mao L, Pang B, Zhao H, Tian G, Jiang H, Kan B. Whole-genome sequencing of rough Brucella melitensis in China provides insights into its genetic features. Emerg Microbes Infect. 2020;9(1):2147–56.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Yang X, Wang Y, Li J, Chen J, Liu J, Tian G, Zhao H, Piao D, Fan Y, Jiang H. Genetic characteristics of an amikacin-resistant Brucella abortus strain first isolated from Marmota himalayana. Microb Pathog. 2022;164:105402.

    Article  PubMed  CAS  Google Scholar 

  12. Whatmore AM, Foster JT. Emerging diversity and ongoing expansion of the genus Brucella. Infect Genet Evol. 2021;92:104865.

    Article  PubMed  CAS  Google Scholar 

  13. Whatmore AM. Current understanding of the genetic diversity of Brucella, an expanding genus of zoonotic pathogens. Infect Genet Evol. 2009;9(6):1168–84.

    Article  PubMed  CAS  Google Scholar 

  14. Suarez-Esquivel M, Chaves-Olarte E, Moreno E, Guzman-Verri C. Brucella Genomics: Macro and Micro Evolution. Int J Mol Sci. 2020;21(20):7749.

  15. Tan KK, Tan YC, Chang LY, Lee KW, Nore SS, Yee WY, Mat Isa MN, Jafar FL, Hoh CC, AbuBakar S. Full genome SNP-based phylogenetic analysis reveals the origin and global spread of Brucella melitensis. BMC Genomics. 2015;16(1):93.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Alton GG, Jones LM, Pietz DE. Laboratory techniques in brucellosis. Monogr Ser World Health Organ. 1975;55:1–163.

    Google Scholar 

  17. Le Fleche P, Jacques I, Grayon M, Al Dahouk S, 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.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Liu Z, Wang C, Wei K, Zhao Z, Wang M, Li D, Wang H, Wei Q, Li Z. Investigation of genetic relatedness of Brucella Strains in Countries along the Silk Road. Front Vet Sci. 2020;7:539444.

    Article  PubMed  Google Scholar 

  19. Zhu X, Zhao Z, Ma S, Guo Z, Wang M, Li Z, Liu Z. Brucella melitensis, a latent travel bacterium, continual spread and expansion from Northern to Southern China and its relationship to worldwide lineages. Emerg Microbes Infect. 2020;9(1):1618–27.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. An CH, Nie SM, Sun YX, Fan SP, Luo BY, Li Z, Liu ZG, Chang WH. Seroprevalence trend of human brucellosis and MLVA genotyping characteristics of Brucella melitensis in Shaanxi Province, China, during 2008–2020. Transbound Emerg Dis. 2022;69(4):e423–34.

    Article  PubMed  CAS  Google Scholar 

  21. Li Z, Wang XM, Zhu X, Wang M, Cheng H, Li D, Liu ZG. Molecular characteristics of Brucella isolates collected from humans in Hainan Province, China. Front Microbiol. 2020;11:452.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Yuan HT, Wang CL, Liu LN, Wang D, Li D, Li ZJ, Liu ZG. Epidemiologically characteristics of human brucellosis and antimicrobial susceptibility pattern of Brucella melitensis in Hinggan League of the Inner Mongolia Autonomous Region, China. Infect Dis Poverty. 2020;9(1):79.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Sun M, Jing Z, Di D, Yan H, Zhang Z, Xu Q, Zhang X, Wang X, Ni B, Sun X, et al. Multiple locus variable-number Tandem-repeat and single-nucleotide polymorphism-based Brucella typing reveals multiple lineages in Brucella melitensis currently endemic in China. Front Vet Sci. 2017;4:215.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Zhao ZJ, Li JQ, Ma L, Xue HM, Yang XX, Zhao YB, Qin YM, Yang XW, Piao DR, Zhao HY, et al. Molecular characteristics of Brucella melitensis isolates from humans in Qinghai Province, China. Infect Dis Poverty. 2021;10(1):42.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Wang H, Xu WM, Zhu KJ, Zhu SJ, Zhang HF, Wang J, Yang Y, Shao FY, Jiang NM, Tao ZY, et al. Molecular investigation of infection sources and transmission chains of brucellosis in Zhejiang, China. Emerg Microbes Infect. 2020;9(1):889–99.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Liu ZG, Wang M, Zhao HY, Piao DR, Jiang H, Li ZJ. Investigation of the molecular characteristics of Brucella isolates from Guangxi Province, China. BMC Microbiol. 2019;19(1):292.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Cao X, Li S, Li Z, Liu Z, Ma J, Lou Z, Zhou J, Liu Y, Jing Z, Fu B. Enzootic situation and molecular epidemiology of Brucella in livestock from 2011 to 2015 in Qingyang, China. Emerg Microbes Infect. 2018;7(1):58.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Bounaadja L, Albert D, Chenais B, Henault S, Zygmunt MS, Poliak S, Garin-Bastuji B. Real-time PCR for identification of Brucella spp.: a comparative study of IS711, bcsp31 and per target genes. Vet Microbiol. 2009;137(1–2):156–64.

    Article  PubMed  CAS  Google Scholar 

  29. Yang X, Wu T, Liu W, Tian G, Zhao H, Piao D, Jiang H, Wu Q. Cell membrane components of Brucella melitensis play important roles in the resistance of low-level rifampicin. PLoS Negl Trop Dis. 2020;14(12):e0008888.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Faison WJ, Rostovtsev A, Castro-Nallar E, Crandall KA, Chumakov K, Simonyan V, Mazumder R. Whole genome single-nucleotide variation profile-based phylogenetic tree building methods for analysis of viral, bacterial and human genomes. Genomics. 2014;104(1):1–7.

    Article  PubMed  CAS  Google Scholar 

  31. Valdezate S, Navarro A, Medina-Pascual MJ, Carrasco G, Saez-Nieto JA. Molecular screening for rifampicin and fluoroquinolone resistance in a clinical population of Brucella melitensis. J Antimicrob Chemother. 2010;65(1):51–3.

    Article  PubMed  CAS  Google Scholar 

  32. Sun MJ, Di DD, Li Y, Zhang ZC, Yan H, Tian LL, Jing ZG, Li JP, Jiang H, Fan WX. Genotyping of Brucella melitensis and Brucella abortus strains currently circulating in Xinjiang, China. Infect Genet Evol. 2016;44:522–9.

    Article  PubMed  CAS  Google Scholar 

  33. Scholz HC, Vergnaud G. Molecular characterisation of Brucella species. Rev Sci Tech. 2013;32(1):149–62.

    Article  PubMed  CAS  Google Scholar 

  34. Zhang F, Li Z, La X, Ma X, Zhang Y, Ji P, Jiang M, Hu J, Zhang Z, Lu X, et al. Multiple-locus variable-number tandem-repeat analysis of Brucella isolates from patients in Xinjiang China. Int J Clin Exp Med. 2015;8(9):15716–23.

    PubMed  PubMed Central  CAS  Google Scholar 

  35. Sayan M, Kilic S, Uyanik MH. Epidemiological survey of rifampicin resistance in clinic isolates of Brucella melitensis obtained from all regions of Turkey. J Infect Chemother. 2012;18(1):41–6.

    Article  PubMed  CAS  Google Scholar 

  36. Abdel-Maksoud M, House B, Wasfy M, Abdel-Rahman B, Pimentel G, Roushdy G, Dueger E. In vitro antibiotic susceptibility testing of Brucella isolates from Egypt between 1999 and 2007 and evidence of probable rifampin resistance. Ann Clin Microbiol Antimicrob. 2012;11:24.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Johansen TB, Scheffer L, Jensen VK, Bohlin J, Feruglio SL. Whole-genome sequencing and antimicrobial resistance in Brucella melitensis from a Norwegian perspective. Sci Rep. 2018;8(1):8538.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Barbosa Pauletti R, Reinato Stynen AP, Pinto da Silva Mol J, Seles Dorneles EM, Alves TM, de Sousa Moura Souto M, Minharro S, Heinemann MB, Lage AP. Reduced susceptibility to Rifampicin and Resistance to multiple Antimicrobial agents among Brucella abortus isolates from cattle in Brazil. PLoS ONE. 2015;10(7):e0132532.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Shevtsov A, Syzdykov M, Kuznetsov A, Shustov A, Shevtsova E, Berdimuratova K, Mukanov K, Ramankulov Y. Antimicrobial susceptibility of Brucella melitensis in Kazakhstan. Antimicrob Resist Infect Control. 2017;6:130.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Deshmukh A, Hagen F, Sharabasi OA, Abraham M, Wilson G, Doiphode S, Maslamani MA, Meis JF. In vitro antimicrobial susceptibility testing of human Brucella melitensis isolates from Qatar between 2014–2015. BMC Microbiol. 2015;15:121.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Hashim R, Ahmad N, Mohamed Zahidi J, Tay BY, Mohd Noor A, Zainal S, Hamzah H, Hamzah SH, Chow TS, Wong PS, et al. Identification and in vitro antimicrobial susceptibility of Brucella species isolated from human brucellosis. Int J Microbiol. 2014;2014:596245.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Heo EJ, Kang SI, Kim JW, Her M, Cho D, Cho YS, Hwang IY, Moon JS, Wee SH, Jung SC, et al. In vitro activities of antimicrobials against Brucella abortus isolates from cattle in Korea during 1998–2006. J Microbiol Biotechnol. 2012;22(4):567–70.

    Article  PubMed  CAS  Google Scholar 

  43. Liu ZG, Di DD, Wang M, Liu RH, Zhao HY, Piao DR, Zhao ZZ, Hao YQ, Du YN, Jiang H, et al. In vitro antimicrobial susceptibility testing of human Brucella melitensis isolates from Ulanqab of Inner Mongolia, China. BMC Infect Dis. 2018;18(1):43.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Torkaman Asadi F, Hashemi SH, Alikhani MY, Moghimbeigi A, Naseri Z. Clinical and diagnostic aspects of brucellosis and Antimicrobial susceptibility of Brucella isolates in Hamedan, Iran. Jpn J Infect Dis. 2017;70(3):235–8.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Ilhan Z, Solmaz H, Ekin IH. In vitro antimicrobial susceptibility of Brucella melitensis isolates from sheep in an area endemic for human brucellosis in Turkey. J Vet Med Sci. 2013;75(8):1035–40.

    Article  PubMed  CAS  Google Scholar 

  47. Wu JY, Li JJ, Wang DF, Wei YR, Meng XX, Tuerxun G, Bolati H, Liu KK, Muhan M, Shahan A, et al. Seroprevalence of five zoonotic pathogens in Wild ruminants in Xinjiang, Northwest China. Vector Borne Zoonotic Dis. 2020;20(12):882–7.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to express our deep gratitude to the staff for their effort and assistance with the strain collection at different levels in CCDC and CADC. Furthermore, we are grateful to Haijian Zhou (National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention), Qingmin Wu (College of Veterinary Medicine, China Agricultural University) and Hai Jiang (State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention) for their advices in setting up this study.

Funding

This work was supported by the Innovation Program of Chinese Academy of Agricultural Sciences (No. CAAS-CSLPDCP-202403), Beijing Municipal Natural Science Foundation (No. 7222120), The Agricultural Science and Technology Innovation Program (ASTIP) (No. ASTIP-IAS-15), Central Public-interest Scientific Institution Basal Research Fund (No. 2023-YWF-ZYSQ-09), National Natural Science Foundation of China (No. 81902105), and grant from the State Key Laboratory of Infectious Disease Prevention and Control, China CDC) (No. 2021SKLID508).

Author information

Authors and Affiliations

  1. Key Laboratory of Animal Biosafe Risk Prevention and Control (North), Ministry of Agriculture and Rural Affairs, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China

    Xiaowen Yang, Lin Liang, Hui Jiang & Jiabo Ding

  2. College of Veterinary Medicine, Xinjiang Agricultural University, Urumuqi, 830052, China

    Yan Liu

  3. Institute of Laboratory Animals Science, CAMS &PUMC, Beijing, 100021, China

    Na Li

  4. National Reference Laboratory for Animal Brucellosis, China Institute of Veterinary Drug Control, Beijing, 102600, China

    Xiaowei Peng, Yinghui Zhang, Xiaoqian Zhang & Zengjie Bian

Authors
  1. Xiaowen Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yan Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Na Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xiaowei Peng
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yinghui Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xiaoqian Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Lin Liang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Zengjie Bian
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Hui Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Jiabo Ding
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization and study design: XWY, HJ, JBD; strains collection and identification: YL, ZJB; data analysis: XWY, LL; laboratory work: XWY, YL, NL, XWP, YHZ, XQZ, LL, ZJB; writing-review & editing: HJ, JBD. All the authors read and approved the final version of the manuscript.

Corresponding authors

Correspondence to Hui Jiang or Jiabo Ding.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, X., Liu, Y., Li, N. et al. Analysis of the Brucella melitensis epidemic in Xinjiang: genotyping, polymorphism, antibiotic resistance and tracing. Ann Clin Microbiol Antimicrob 23, 71 (2024). https://doi.org/10.1186/s12941-024-00724-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12941-024-00724-0

Keywords

Annals of Clinical Microbiology and Antimicrobials

ISSN: 1476-0711

Contact us