Skip to main content

Molecular characterization of fluoroquinolone and/or cephalosporin resistance in Shigella sonnei isolates from yaks



Members of the genus Shigella are intestinal pathogens and a major cause of seasonal outbreaks of bacterial diarrhea worldwide. Although humans are the conventional hosts of Shigella species, expansion of the Shigella host range to certain animals was recently reported. To investigate the prevalence of Shigella sonnei (S. sonnei) in yaks and perform molecular characterization, we analyzed 1132 fresh yak diarrheal stool samples and collected a total of 44 S. sonnei isolates.


We performed multilocus sequence typing (MLST) and pulsed-field gel electrophoresis (PFGE) with XbaI-digested DNA to study genetic relatedness among the 44 isolates, which were differentiated into 4 sequence types (STs) and 32 PFGE types (PTs). All isolates harbored virulence genes, and 87.36% tested positive for invasion plasmid antigen H (ipaH), invasion associated locus (ial) and the Shigella enterotoxin gene sen. According to the results of antimicrobial susceptibility tests, 45.45% (20/44) were resistant to fluoroquinolones and/or cephalosporin. By sequencing the quinolone resistance determining region (QRDR) genes, we identified double mutations in gyrA (Ser83-Leu and Asp87-Asn) and a single mutation in parC (Ser80-Ile). All 12 fluoroquinolone-resistant S. sonnei isolates tested positive for the aac(6′)-Ib-cr gene but negative for qepA. Three isolates harbored qnr genes, including two with qnrS and one with qnrB. In addition, three types of β-lactamase genes, bla TEM-1 , bla OXA-1 and bla CTX-M-14/79 , were detected in cephalosporin-resistant isolates.


The findings of this study have enriched our knowledge of fluoroquinolone- and/or cephalosporin-resistant S. sonnei isolates from yaks, which has important public health significance.


Bacteria of the genus Shigella are important members of the Enterobacteriaceae family and cause acute diarrhea in humans and animals [1,2,3]. This pathogen has been historically divided into 4 subgroups based on biochemical properties and group-specific O antigens in the outer membrane, specifically subgroups A (S. dysenteriae), B (S. flexneri), C (S. boydii) and D (S. sonnei). All four Shigella species cause shigellosis and are phylogenetically distinct from Escherichia coli [4, 5], with S. flexneri and S. sonnei being the most prevalent. Notably, S. sonnei has become the dominant subgroup in Asian countries in recent years [6, 7].

The specific pathogenicity of a given Shigella isolate is derived from the expression of diverse virulence genes that are associated with colonization, invasion/penetration and toxin-mediated disease [8]. The pathogenicity of these virulence genes is often multifactorial, and they are coordinately regulated [9]. An investigation of the virulence determinant genes in Shigella would help us better understand its pathogenicity. However, knowledge regarding the distribution of these genes in Shigella isolates from animals is limited.

For Shigella infections, prompt therapy with effective antimicrobial agents reduces the duration and severity of illness [10]. However, the progressive increase in resistance to commonly used antimicrobials, especially fluoroquinolones and extended-spectrum cephalosporins, among Shigella spp. poses a major therapeutic challenge to the control of diseases caused by these species [10, 11]. In addition, the emergence and global dissemination of multi-drug resistant (MDR) pathogens is increasing rapidly due to the unique ability of bacteria to acquire resistance factors (transmissible genes) from the environment or from other bacteria [12].

Foodborne pathogens are usually carried by the host animal and contaminate the food during the slaughtering process [13]. And food is the main route to human infection with Shigella, including antibiotic-resistant strains [14]. Foods implicated in human cases of shigellosis include fresh fruit and vegetables, raw oysters, deli meats, and unpasteurized milk [15]. Multidrug-resistant Shigella has been found in bovine-derived foods such as beef, milk and cheese, which predict potential threats to human [16].

The primary goal of this study was to investigate the incidence of S. sonnei in yaks with diarrheal disease in a plateau area. We analyzed the biochemical and serological characteristics, virulence gene profiles, antimicrobial resistance profiles, antimicrobial resistance genes (ARGs) and genotypes and genetic relatedness diversity of S. sonnei isolates.


Identification of S. sonnei isolates

A total of 44 S. sonnei isolates were collected from 1132 fresh diarrheal stool samples from yaks in three Chinese provinces, Gansu, Qinghai and Tibet, from 2014 to 2016. Among these isolates, 24 (54.55%) were collected from Gansu, 15 (34.09%) were collected from Qinghai, and 5 (11.36%) were collected from Tibet. Based on the results of biochemical characterization assays, we observed that all 44 S. sonnei isolates possessed the typical biochemical features of Shigella species and demonstrated the ability to ferment I-D-galactopyranoside (ONPG), glucose (GLU), mannitol (MAN), melibiose (MEL) and arabinose (ARA). Although both serotypes were isolated in each region, phase I (33/44, 75%) was the predominant serotype, and the detection rate of S. sonnei phase I was three times higher than that of phase II (11/44, 25%).

MLST-based genotype analysis

MLST was performed to analyze the genotypic diversity of S. sonnei isolates based on 15 housekeeping genes. The 44 isolates were divided into 4 sequence types (STs): ST76, ST116, ST125 and ST155. All four STs were previously reported and belonged to the same clonal complex (CC; CC29). The allele number for each locus and the ST designation are listed in Additional file 1: Table S1 and shared on the EcMLST website. Among these STs, ST116 (n = 27) and ST155 (n = 13) were the most common, accounting for 90.91% of all STs. ST116 encompassed 16 phase I and 11 phase II serotype isolates; furthermore, all phase II serotype S. sonnei isolates in the study were ST116. Interestingly, all four STs were detected in Gansu province, while only ST116 was observed in Qinghai and Tibet provinces (Table 1).

Table 1 Statistical analysis of the distribution of 4 MLST types among 44 S. sonnei isolates in different provinces

PFGE-based genotype analysis

The genotypes and genetic relatedness diversity of the 44 S. sonnei isolates were assessed by PFGE. At 80% similarity, XbaI-digested S. sonnei DNA generated 32 reproducible unique PFGE banding patterns (PT), yielding three major groups: A (n = 17), B (n = 25), and C (n = 2) (Fig. 1). Most isolates contained in PT group A were collected from the same geographical region (Gansu province) with the exception of one Tibet isolate (SS036). Group B contained the highest number of isolates, accounting for 56.82% of total isolates. In addition, almost all Qinghai isolates were clustered in group B2, with the exception of SS017. Interestingly, we also observed that different serotypes affect the PT cluster. In this study, all 17 group A isolates belonged to the phase I serotype, while all of the phase II isolates clustered into group B. The large number of PT types suggested a capricious genotypic and genetic diversity among S. sonnei strains, which adapt to different environments.

Fig. 1

Dendrogram of 44 XbaI-digested S. sonnei isolates based on the cluster analysis of PFGE patterns. The dendrogram was constructed using the UPGMA clustering method. The corresponding antibiotic resistance profile, PFGE pattern and background information for each strain is listed on the right side of the dendrogram. R = resistance and no bacteriostatic ring; r = resistance and positive bacteriostatic ring; N = sensitive and intermediary bacteriostatic ring

Virulence genes

The virulence gene profiles of the isolates are listed in Additional file 2: Table S2. Only 3 virulence genes (ipaH, ial and sen) were observed in this study. Among them, the invasive plasmid gene ipaH was observed in all isolates, followed by ial (38/44, 86.36%) and sen (37/44, 84.09%). All isolates tested negative for set1A, set1B and stx.

According to their virulence gene distribution, the 44 S. sonnei isolates formed 3 gene profile patterns (Additional file 3: Table S3). VT3, which was positive for ipaH, ial and sen, contained the highest number of isolates (37, 87.36%) and was the major VT present in each province. VT2, which was positive for ipaH and ial, was only detected in one Qinghai isolate (SS005). Six isolates (13.64%), one in Gansu and five in Qinghai, belonged to VT1, which was only positive for ipaH.

Antimicrobial susceptibility profiles

The antimicrobial resistance profiles of the 44 S. sonnei isolates to 21 antimicrobials are listed in Fig. 1 and Table 2. The antimicrobial susceptibility results showed that all S. sonnei isolates were sensitive to ciprofloxacin, ofloxacin, amoxicillin/clavulanic acid, cefepime, meropenem, imipenem, and amikacin but resistant to penicillin G. The resistance rates to other antimicrobials are as follows: ampicillin (31/44, 70.45%), tetracycline (30/44, 68.18%), cephazolin (14/44, 31.82%), streptomycin (14/44, 31.82%), cephalothin (13/44, 29.55%), enrofloxacin (12/44, 27.27%), gentamicin (12/44, 27.27%), norfloxacin (10/44, 22.73%), ceftriaxone (9/44, 20.45%), cefotaxime (5/44, 11.36%), levofloxacin (4/44, 9.09%), cefoxitin (4/44, 9.09%), and chloramphenicol (4/44, 9.09%). Although the observed resistance to most antimicrobials (with the exception of penicillin G, ampicillin and tetracycline) was less than 50%, multidrug-resistant (MDR) S. sonnei isolates, which were defined as being resistant to 3 or more antimicrobial agent subclasses, accounted for 43.18% of isolates (19/44).

Table 2 Statistical analysis of the results of antimicrobial susceptibility testing of 21 antibiotics for 44 S. sonnei isolates

Remarkably, the rates of resistance to fluoroquinolones and cephalosporin were 27.275% (12/44) and 31.82% (14/44), respectively, and 13.64% (6/44) of isolates were resistant to both fluoroquinolones and cephalosporin. Moreover, all fluoroquinolone-resistant isolates belonged to the MDR group. Interestingly, the fluoroquinolone- and/or cephalosporin-resistant isolates demonstrated diverse antimicrobial resistance profiles, and no isolates had the same resistance profile (Table 3 and Fig. 1).

Table 3 Analysis of amino acid types in QRDRs and ARGs in S. sonnei isolates with resistance to fluoroquinolones and/or cephalosporin

Molecular analysis of antimicrobial resistance determinants

In the present study, all four QRDR genes (gyrA, gyrB, parC, and parE) were sequenced and compared with reference sequences. According to the results, double mutations in gyrA (Ser83-Leu and Asp87-Asn) and a single mutation in parC (Ser80-Ile) were associated with fluoroquinolone resistance in S. sonnei isolates (Table 3). The encoded Ser83-Leu amino acid substitution in gyrA was observed in every fluoroquinolone-resistant isolate, and 66.67% (8/12) of fluoroquinolone-resistant isolates possessed the Asp87-Asn mutation in gyrA. Only four isolates harbored all three amino acid substitutions. Interestingly, all four of these isolates belonged to the phase I serotype.

In addition to the identification of substitutions in QRDR genes, PMQR genes were amplified and sequenced. According to the results, all fluoroquinolone-resistant isolates harbored the aac(6′)-Ib-cr gene. The qepA gene, which encodes an efflux pump, was not observed in the studied isolates. Only three Qinghai isolates tested positive for the qnr gene, with two having qnrS and one having qnrB (Table 3).

For the 14 cephalosporin-resistant isolates, only three β-lactamase gene types (bla TEM , bla OXA and bla CTX-M ) were observed. All isolates tested positive for the bla TEM gene and showed 100% identity with the bla TEM − 1 subtype. And all the third-cephalosporin resistant isolates were positive for bla OXA (bla OXA-1 ), except SS027. Five out of the fourteen cephalosporin-resistant isolates tested positive for bla CTX-M , 4 of which had the bla CTX-M-14 gene, and 1 had the bla CTX-M-79 gene.


Bacillary dysentery caused by Shigella is endemic throughout the world. Humans are the natural hosts of Shigella, although reports of Shigella infections in animals (fish, calves, piglets, and chickens) have emerged [2, 17,18,19]. We previously confirmed the existence of diverse S. sonnei isolates in yak in local epidemiological studies. The rate of S. sonnei isolation from the diarrhea of yaks was 3.89% (44/1132), and the predominant serotype was phase I, which was isolated three times more frequently than phase II. In addition, 54 S. flexneri isolates with diverse serotypes were reported in our previous study [19].

MLST and PFGE have shown comparable discrimination in terms of their ability to subtype S. sonnei. In this study, 44 S. sonnei isolates from yaks were grouped into 4 STs based on 15 housekeeping genes. Comparison of these genes with the standard gene types on the EcMLST website revealed that the different STs of these isolates were primarily attributable to three genes, aspC, mutS, and ropS. In addition, all isolates belonged to the same clonal complex, suggesting a relatively close genetic relationship among them. In our study, the majority of isolates (n = 27 and n = 13) belonged to ST116 and ST155, respectively. However, ST76, ST116, and ST114 were previously reported in human isolates [20, 21]. PFGE is a broadly applicable typing method that has high resolution for the molecular characterization of several enteric bacteria, including Shigella [20]. Based on our PFGE analysis, the S. sonnei isolates in the present study were heterogeneous and distributed into 32 PTs. The clustering of these diverse PTs allows us to learn more about the epidemiological characteristics of S. sonnei in specific geographical regions. PFGE appeared to have a higher discriminatory power for genotypes and genetic relatedness diversity analysis than MLST.

Virulence genes are responsible for the invasion of virulent Shigella strains into intestinal epithelial cells, resulting in dysentery in hosts [22]. ipaH, which is typically used as a marker of Shigella, was harbored by all isolates. Additionally, in this study, most isolates simultaneously harbored the ial and sen genes. The diversity of the observed virulence genes suggested the isolates were pathogenic. Moreover, virulence genes were frequently located in megaplasmids, which allow Shigella to acquire and disseminate bacterial invasion genes through mobile genetic elements [23].

Recently, antimicrobial-resistant and MDR S. sonnei isolates have been widely isolated [20, 23]. In this study, we isolated MDR S. sonnei isolates from yaks. More than half of the isolates were resistant to conventional and commonly used antimicrobials: penicillin G (100%), ampicillin (70.45%), and tetracycline (68.18%). Fortunately, for most antimicrobials, the rates and levels of resistance of Shigella isolates from yaks were less than those from isolates from livestock [24]. Notably, although only low rates of resistance and MIC values to fluoroquinolones and/or cephalosporin were observed, some S. sonnei had acquired the ability to survive the selective pressures of these antimicrobials, which will limit the control and treatment of shigellosis.

Fluoroquinolone resistance is primarily mediated by the accumulation of sequential mutations in four QRDR genes: gyrA and gyrB, which encode DNA gyrase, and parC and parE, which encode and topoisomerase IV [25]. The primary mutations are located in gyrA and parC. According to a previous study, the major mutations in gyrA codons 81 (Gly → Ser), 83 (Ser → Leu), 87 (Asp→Asn/Gly/Tyr), and 211 (His→Tyr) and parC codons 80 (Ser → Ile), 83 (Ser → Leu), and 129 (Ser → Pro) are associated with quinolone and/or fluoroquinolone resistance in Shigella [11, 24, 26, 27]. Furthermore, the continuous accumulation of mutations at multiple positions and different mutations at the same position may lead to diverse levels of resistance [23]. In this study, we identified QRDR mutations in all fluoroquinolone-resistant S. sonnei, observing only a few classic amino acid substitutions in gyrA and parC.

The PMQR determinants are typically located on mobile or transposable genetic elements, which may facilitate dissemination among Shigella species and other members of the Enterobacteriaceae family [28, 29]. Recently, the diversity of PMQR determinants has emerged as an important issue worldwide [27]. Among the PMQR genes, aac(6′)-Ib-cr, which encodes aminoglycoside acetyl transferase and is responsible for reduced fluoroquinolone activity, was the most prevalent in our study. The qnr-encoded family of proteins protects DNA gyrase from quinolones and confers low-level resistance [30]. QnrS were majority prevalence in Shigella (S.sonnei and S.flexneri) in China, while qnrB previously in India [23, 26]. In addition, the plasmid-mediated efflux pump-encoding gene qepA was detected at lower frequencies [27]. The presence of PMQR genes may help microorganisms develop resistance by conferring mutations and facilitating the selection of higher-level quinolone resistance [31].

ESBLs continue to be the major defense mechanism against cephalosporins among Gram-negative bacteria [32]. Since first being identified in the 1980s, different types of ESBLs belonging to the TEM, OXA, SHV and CTX-M families have been widely reported in Enterobacteriaceae, including Shigella species [33, 34]. The global spread of ESBL-producing Shigella is difficult to control and eradicate because these bacteria reside in the intestines of humans and animals and harbor epidemic resistance plasmids. Therefore, the prevalence of ESBL-producing S. sonnei in yaks represents a substantial threat.


Antimicrobial-resistant Shigella remains a serious threat to humans and animals; however, knowledge of the molecular epidemiology of this pathogen in yaks is limited. Shigella species typically harbor various plasmids that are associated with invasion into host intestinal epithelial cells and antimicrobial resistance and may be transferred between different bacteria. In the present study the S. sonnei strains isolated from yaks that tested positive for diverse virulence genes and ARGs represent a potential threat to hosts, including humans. Therefore, continuous and extensive surveillance will be necessary to control and reduce the threat of foodborne pathogens, such as Shigella, Salmonella and E. coli.


Bacterial isolates and serological confirmation

A total of 44 S. sonnei isolates were obtained from yaks with diarrhea from October 2014 to December 2016. All yaks have hardly been exposed to any antibiotic except for a few of penicillin to against the disease. All isolates were collected directly from fresh stool samples, plated on Salmonella-Shigella (SS) selective agar and confirmed on MacConkey (MAC) agar at 37 °C for 18 to 24 h. Colorless, semitransparent, smooth, and moist circular plaques were considered presumptive Shigella isolates and were selected for further confirmation [34]. The presumptive positive Shigella isolates were identified using API 20 E kits (bioMerieux, Marcy-l’Etoile, France), and serotypes were tested with a commercially available kit (Denka Seiken; Tokyo, Japan) according to the manufacturer’s recommendations.

Multilocus sequence typing (MLST)

MLST was performed for each isolate according to the protocols described on the EcMLST website ( The following PCR amplification conditions were used for 15 housekeeping genes: 94 °C for 5 min, followed by 30 cycles of 94 °C for 30 s, 55 °C for 30 s, and 72 °C for 1 min, with incubation for 10 min at 72 °C with ExTaq DNA polymerase (Takara; Dalian, China). The PCR products were sequenced bi-directionally, and sequences of the 15 housekeeping genes were edited using SeqMan 7.0. Finally, the sequences were uploaded to the EcMLST website for comparison, which allowed us to determine the gene and ST type [35].

Pulsed-field gel electrophoresis (PFGE)

All S. sonnei isolates were analyzed for genetic relatedness by PFGE after digestion with the restriction enzyme XbaI (TaKaRa; Japan) according to the recommendations of a previous study [36]. Salmonella enterica serotype Braenderup strain H9812 was used as a molecular size maker. Electrophoresis was performed with a Bio-Rad CHEF Mapper XA system (Bio-Rad; USA) in a 1% SeaKem Gold agarose gel (Lonza; USA) in 0.5X TBE buffer with a size range of 30–700 kb. The PFGE run conditions were 6 V/cm with a switch from 2.16 to 54.17 s for 21 h at 14 °C.

After gel electrophoresis, gels were stained with ethidium bromide, rinsed, and photographed with a Universal Hood II (Bio-Rad; USA). The PFGE patterns were analyzed with BioNumerics using the Dice similarity coefficient, unweighted pair-group method with arithmetic mean (UPGMA) and 1.0% tolerance.

Detection of virulence genes by PCR

All S. sonnei isolates were examined for the presence of six virulence genes, invasion plasmid antigen H (ipaH), invasion associated locus (ial), Shiga toxin gene (stx), and Shigella enterotoxin genes (set1A, set1B, and sen), by PCR according to published procedures [23, 37]. The primers for these virulence genes are listed in Additional file 4: Table S4. Amplification products were separated by 1% agarose gel electrophoresis and stained with ethidium bromide.

Antimicrobial susceptibility

Antimicrobial susceptibility testing was performed by the Kirby-Bauer disc diffusion method on Muller-Hinton agar (MHA) plates following the guidelines of the Clinical and Laboratory Standards Institute (CLSI) using commercially available antimicrobial discs (Oxoid, UK). Twenty-one antimicrobial discs were used in this study: norfloxacin (NOR, 10 μg), enrofloxacin (ENR, 5 μg), levofloxacin (LEV, 5 μg), ciprofloxacin (CIP, 5 μg), ofloxacin (OFX, 5 μg), penicillin G (P, 10 μg), ampicillin (AMP, 10 μg), amoxicillin/clavulanic acid (AMC, 30 μg), cephalothin (KF, 30 μg), cephazolin (KZ, 30 μg), cefoxitin (FOX, 30 μg), ceftriaxone (CRO, 30 μg), cefotaxime (CTX, 30 μg), cefepime (FEP, 30 μg), imipenem (IPM, 10 μg), meropenem (MEM, 10 μg), chloramphenicol (C, 30 μg), tetracycline (TE, 30 μg), streptomycin (S, 10 μg), gentamicin (CN, 10 μg), and amikacin (AK, 30 μg). E. coli ATCC 25922 was used as a control strain. The standard of antibiotic susceptibility for CLSI were list in Additional file 5: Table S5.

Detection of β-lactamase and quinolone resistance genes

Cephalosporin- and/or fluoroquinolone-resistant S. sonnei isolates were assayed by PCR using 24 primer panels to determine the underlying mechanism conferring resistance (Additional file 6: Table S6). Specifically, 6 plasmid-mediated quinolone resistance (PMQRs) determinant genes (qnrA, qnrB, qnrD, qnrS, qepA and aac(6′)-Ib-cr) and 4 quinolone resistance determining region (QRDR) genes (gyrA, gyrB, parC and parE) were amplified to clarify the underlying mechanism conferring resistance to quinolones [38,39,40]. Extended-spectrum β-lactamase (ESBL) genes (bla CTX-M , bla SHV , bla TEM , and bla OXA ) were detected in cephalosporin-resistant isolates [39, 41, 42].



Extended-spectrum β-lactamase


Multilocus sequence typing


Pulsed-field gel electrophoresis


Plasmid-mediated quinolone resistance


Quinolone resistance determining region


Virulence gene profile types


  1. 1.

    Stirling J, Griffith M, Blair I, Cormican M, Dooley JS, Goldsmith CE, Glover SG, Loughrey A, Lowery CJ, Matsuda M, McClurg R, McCorry K, McDowell D, McMahon A, Cherie Millar B, Nagano Y, Rao JR, Rooney PJ, Smyth M, Snelling WJ, Xu J, Moore JE. Prevalence of gastrointestinal bacterial pathogens in a population of zoo animals. Zoonoses Public Health. 2008;55(3):166–72.

    Article  PubMed  CAS  Google Scholar 

  2. 2.

    Shi R, Yang X, Chen L, Chang HT, Liu HY, Zhao J, Wang XW, Wang CQ. Pathogenicity of Shigella in chickens. PLoS One. 2014;9(6):e100264.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. 3.

    Parajuli P, Adamski M, Verma NK. Bacteriophages are the major drivers of Shigella flexneri serotype 1c genome plasticity: a complete genome analysis. BMC Genomics. 2017;18(1):722.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Pupo GM, Lan R, Reeves PR. Multiple independent origins of Shigella clones of Escherichia coli and convergent evolution of many of their characteristics. Proc Natl Acad Sci U S. 2000;7(19):10567–72.

    Article  Google Scholar 

  5. 5.

    Xia S, Xu B, Huang L, Zhao JY, Ran L, Zhang J, Chen H, Pulsrikarn C, Pornruangwong S, Aarestrup FM, Hendriksen RS. Prevalence and characterization of human Shigella infections in Henan Province, China, in 2006. J Clin Microbiol. 2011;49(1):232–42.

    Article  PubMed  Google Scholar 

  6. 6.

    Bangtrakulnonth A, Vieira AR, Lo Fo Wong DM, Pornreongwong S, Pulsrikarn C, Sawanpanyalert P, Hendriksen RS, Aarestrup FM. Shigella from humans in Thailand during 1993 to 2006: spatial-time trends in species and serotype distribution. Foodborne Pathog Dis. 2008;5(6):773–84.

    Article  PubMed  Google Scholar 

  7. 7.

    Qu F, Bao C, Chen S, Cui E, Guo T, Wang H, Zhang J, Wang H, Tang YW, Mao Y. Genotypes and antimicrobial profiles of Shigella sonnei isolates from diarrheal patients circulating in Beijing between 2002 and 2007. Diagn Microbiol Infect Dis. 2012;74(2):166–70.

    Article  PubMed  CAS  Google Scholar 

  8. 8.

    Zhang CL, Liu QZ, Wang J, Chu X, Shen LM, Guo YY. Epidemic and virulence characteristic of Shigella spp. with extended-spectrum cephalosporin resistance in Xiaoshan District, Hangzhou, China. BMC Infect Dis. 2014;14:260.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. 9.

    Qu M, Zhang X, Liu G, Huang Y, Jia L, Liang W, Li X, Wu X, Li J, Yan H, Kan B, Wang Q. An eight-year study of Shigella species in Beijing, China: serodiversity, virulence genes, and antimicrobial resistance. J Infect Dev Ctries. 2014;8(7):904–8.

    Article  PubMed  CAS  Google Scholar 

  10. 10.

    Azmi IJ, Khajanchi BK, Akter F, Hasan TN, Shahnaij M, Akter M, Banik A, Sultana H, Hossain MA, Ahmed MK, Faruque SM, Talukder KA. Fluoroquinolone resistance mechanisms of Shigella flexneri isolated in Bangladesh. PLoS One. 2014;9(7):e102533.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. 11.

    Sabra AH, Araj GF, Kattar MM, Abi-Rached RY, Khairallah MT, Klena JD, Matar GM. Molecular characterization of ESBL-producing Shigella sonnei isolates from patients with bacilliary dysentery in Lebanon. J Infect Dev Ctries. 2009;3(4):300–5.

    PubMed  CAS  Google Scholar 

  12. 12.

    Ud-Din AI, Wahid SU, Latif HA, Shahnaij M, Akter M, Azmi IJ, Hasan TN, Ahmed D, Hossain MA, Faruque AS, Faruque SM, Talukder KA. Changing trends in the prevalence of Shigella species: emergence of multi-drug resistant Shigella sonnei biotype g in Bangladesh. PLoS One. 2013;8(12):e82601.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. 13.

    Schlundt J. Emerging food-borne pathogens. Biomed Environ Sci. 2001;14(1–2):44–52.

    PubMed  CAS  Google Scholar 

  14. 14.

    Newell DG, Koopmans M, Verhoef L, Duizer E, Aidara-Kane A, Sprong H, Opsteegh M, Langelaar M, Threfall J, Scheutz F, van der Giessen J, Kruse H. Food-borne diseases - the challenges of 20 years ago still persist while new ones continue to emerge. Int J Food Microbiol. 2010;139(Suppl 1):S3–15.

    Article  PubMed  Google Scholar 

  15. 15.

    EFSA (European Food Safety Authority). Scientific opinion of the panel on biological hazards on a request from the European food safety authority on foodborne antimicrobial resistance as a biological hazard. EFSA J. 2008;765:1–87.

    Google Scholar 

  16. 16.

    Ahmed AM, Shimamoto T. Molecular characterization of multidrug-resistant Shigella spp. of food origin. Int J Food Microbiol. 2015;194:78–82.

    Article  PubMed  CAS  Google Scholar 

  17. 17.

    Onyango DM, Wandili S, Kakai R, Waindi EN. Isolation of Salmonella and Shigella from fish harvested from the Winam gulf of Lake Victoria. Kenya J Infect Dev Ctries. 2009;3(2):99–104.

    PubMed  Google Scholar 

  18. 18.

    Priamukhina NS, Kilesso VA, Tikhomirov ED. Animal carriers of Shigella and their possible epidemiological importance. Zh Mikrobiol Epidemiol Immunobiol. 1984;11:20–4.

    Google Scholar 

  19. 19.

    Maurelli AT, Routh PR, Dillman RC, Ficken MD, Weinstock DM, Almond GW, Orndorff PE. Shigella infection as observed in the experimentally inoculated domestic pig, Sus scrofa domestica. Microb Pathog. 1998;25(4):189–96.

    Article  PubMed  CAS  Google Scholar 

  20. 20.

    Li S, Wang J, Wei X, Liu Y, You L, Luo X, Tang G, Sun Q, Ye C, Xu J, Wang D. Molecular characterization of Shigella sonnei: an increasingly prevalent etiologic agent of shigellosis in Guizhou Province, Southwest of China. PLoS One. 2016;11(5):e0156020.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. 21.

    Cheng T, Shi X, Yong W, Wang J, Xie G, Ding J. Molecular typing of Shigella sonnei isolates circulating in Nanjing, China, 2007-2011. J Infect Dev Ctries. 2014;8(12):1525–32.

    Article  PubMed  Google Scholar 

  22. 22.

    Ashida H, Ogawa M, Mimuro H, Kobayashi T, Sanada T, Sasakawa C. Shigella are versatile mucosal pathogens that circumvent the host innate immune system. Curr Opin Immunol. 2011;23(4):448–55.

    Article  PubMed  CAS  Google Scholar 

  23. 23.

    Das A, Natarajan M, Mandal J. The emergence of quinolone resistant Shigella sonnei, Pondicherry, India. PLoS One. 2016;11(8):e0160290.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. 24.

    Zhu Z, Cao M, Zhou X, Li B, Zhang J. Epidemic characterization and molecular genotyping of Shigella flexneri isolated from calves with diarrhea in Northwest China. Antimicrob Resist Infect Control. 2017;6:92.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Taneja N, Kumar A, Appannanavar S, Verma G, Sharma M. Plasmid-mediated quinolone resistance in Shigella isolates over a decade in India. J Glob Antimicrob Resist. 2014;2(1):59–60.

    Article  PubMed  Google Scholar 

  26. 26.

    Cui X, Wang J, Yang C, Liang B, Ma Q, Yi S, Li H, Liu H, Li P, Wu Z, Xie J, Jia L, Hao R, Wang L, Hua Y, Qiu S, Song H. Prevalence and antimicrobial resistance of Shigella flexneri serotype 2 variant in China. Front Microbiol. 2015;6:435.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Gu B, Qin TT, Fan WT, Bi RR, Chen Y, Li Y, Ma P. Novel mutations in gyrA and parC among Shigella sonnei strains from Jiangsu Province of China, 2002-2011. Int J Infect Dis. 2017;59:44–9.

    Article  PubMed  CAS  Google Scholar 

  28. 28.

    Folster JP, Pecic G, Bowen A, Rickert R, Carattoli A, Whichard JM. Decreased susceptibility to ciprofloxacin among Shigella isolates in the United States, 2006 to 2009. Antimicrob Agents Chemother. 2011;55(4):1758–60.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. 29.

    Rodríguez-Martínez JM, Machuca J, Cano ME, Calvo J, Martínez-Martínez L, Pascual A. Plasmid-mediated quinolone resistance: two decades on. Drug Resist Updat. 2016;29:13–29.

    Article  PubMed  Google Scholar 

  30. 30.

    Jacoby GA, Gacharna N, Black TA, Miller GH, Hooper DC. Temporal appearance of plasmid-mediated quinolone resistance genes. Antimicrob Agents Chemother. 2009;53(4):1665–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. 31.

    Yang H, Duan G, Zhu J, Zhang W, Xi Y, Fan Q. Prevalence and characterisation of plasmid-mediated quinolone resistance and mutations in the gyrase and topoisomerase IV genes among Shigella isolates from Henan, China, between 2001 and 2008. Int J Antimicrob Agents. 2013;42(2):173–7.

    Article  PubMed  CAS  Google Scholar 

  32. 32.

    Li J, Li B, Ni Y, Sun J. Molecular characterization of the extended-spectrum beta-lactamase (ESBL)-producing Shigella spp. in Shanghai. Eur J Clin Microbiol Infect Dis. 2015;34(3):447–51.

    Article  PubMed  CAS  Google Scholar 

  33. 33.

    Robin F, Beyrouthy R, Bonacorsi S, Aissa N, Bret L, Brieu N, Cattoir V, Chapuis A, Chardon H, Degand N, Doucet-Populaire F, Dubois V, Fortineau N, Grillon A, Lanotte P, Leyssene D, Patry I, Podglajen I, Recule C, Ros A, Colomb-Cotinat M, Ponties V, Ploy MC, Bonnet R. Inventory of Extended-Spectrum-β-Lactamase-Producing Enterobacteriaceae in France as assessed by a multicenter study. Antimicrob Agents Chemother. 2017;61(3):e01911–16.

  34. 34.

    Tajbakhsh M, García Migura L, Rahbar M, Svendsen CA, Mohammadzadeh M, Zali MR, Aarestrup FM, Hendriksen RS. Antimicrobial-resistant Shigella infections from Iran: an overlooked problem? J Antimicrob Chemother. 2012;67(5):1128–33.

    Article  PubMed  CAS  Google Scholar 

  35. 35.

    Yang C, Li P, Zhang X, Ma Q, Cui X, Li H, Liu H, Wang J, Xie J, Wu F, Sheng C, Du X, Qi L, Su W, Jia L, Xu X, Zhao J, Xia S, Zhou N, Ma H, Qiu S, Song H. Molecular characterization and analysis of high-level multidrug-resistance of Shigella flexneri serotype 4s strains from China. Sci Rep. 2016;6:29124.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. 36.

    Qu F, Cui E, Guo T, Li H, Chen S, Liu L, Han W, Bao C, Mao Y, Tang YW. Nasal colonization of and clonal transmission of methicillin-susceptible Staphylococcus aureus among Chinese military volunteers. J Clin Microbiol. 2010;48(1):64–9.

    Article  PubMed  CAS  Google Scholar 

  37. 37.

    Talukder KA, Mondol AS, Islam MA, Islam Z, Dutta DK, Khajanchi BK, Azmi IJ, Hossain MA, Rahman M, Cheasty T, Cravioto A, Nair GB, Sack DA. A novel serovar of Shigella dysenteriae from patients with diarrhoea in Bangladesh. J Med Microbiol. 2007;56(Pt 5):654–8.

    Article  PubMed  CAS  Google Scholar 

  38. 38.

    Hu LF, Li JB, Ye Y, Li X. Mutations in the GyrA subunit of DNA gyrase and the ParC subunit of topoisomerase IV in clinical strains of fluoroquinolone-resistant Shigella in Anhui. China J Microbiol. 2007;45(2):168–70.

    PubMed  CAS  Google Scholar 

  39. 39.

    Tariq A, Haque A, Ali A, Bashir S, Habeeb MA, Salman M, Sarwar Y. Molecular profiling of antimicrobial resistance and integron association of multidrug-resistant clinical isolates of Shigella species from Faisalabad, Pakistan. Can J Microbiol. 2012;58(9):1047–54.

    Article  PubMed  CAS  Google Scholar 

  40. 40.

    Colobatiu L, Tabaran A, Flonta M, Oniga O, Mirel S, Mihaiu M. First description of plasmid-mediated quinolone resistance determinants and β-lactamase encoding genes in non-typhoidal Salmonella isolated from humans, one companion animal and food in Romania. Gut Pathog. 2015;7:16.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. 41.

    Zong Z, Partridge SR, Thomas L, Iredell JR. Dominance of blaCTX-M within an Australian extended-spectrum beta-lactamase gene pool. Antimicrob Agents Chemother. 2008;52(11):4198–202.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. 42.

    Liu J, Zeng T, Su G, Lin LY, Zhao Y, Yang WQ, Xie WX, Zhao ZG, Li GM. The dissemination mode of drug-resistant genes in Enterobacter cloacae. Indian J Med Microbiol. 2015;33(Suppl):87–92.

    PubMed  Google Scholar 

Download references


The authors warmly thank to all veterinarians and personnel of the farms that provided the samples for this study.


This work was financed by grants from the National Natural Science Foundation of China (31272603, 31101836).

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its Additional files]. If the reader requires further information to help interpret the results, please contact the corresponding author.

Author information




ZZ and YXS contributed equally to first author in this work. ZZ conceived, designed, performed the experiments and drafted the manuscript the study. YXS collected the samples from Qinghai, performed PFGE and MLST work. XZZ collected the samples from Gansu and isolated Shigella sonneri. BL participated in PCR and sequencing experiments. JYZ designed the study, helped with and supervised the execution of laboratory work, as well as the statistical interpretation of the data and reviewed the final paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jiyu Zhang.

Ethics declarations

Ethics approval

This study has gained ethical approval for this study. Our study was conducted according to the Ethics Committee of Animal Experiments at the Institute of Husbandry and Pharmaceutical Sciences of CAAS in Lanzhou, China. And we gained consent from the owners of the animals for them to be used in the study.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. MLST allelic profiles and ST designations of 44 S. sonnei isolates from this study. (DOCX 45 kb)

Additional file 2:

Table S2. Statistical analysis of the presence of virulence genes in each S. sonnei isolate. (DOCX 39 kb)

Additional file 3:

Table S3. Statistical analysis of the occurrence of each virulence gene profile in different provinces. (DOCX 36 kb)

Additional file 4:

Table S4. Primers used to detect virulence genes. (DOCX 39 kb)

Additional file 5:

Table S5. The standard of antibiotic susceptibility for K-B disc-diffusion method. (DOCX 59 kb)

Additional file 6:

Table S6. Primers used to detect antibiotic resistance determinant genes. (DOCX 81 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhu, Z., Shi, Y., Zhou, X. et al. Molecular characterization of fluoroquinolone and/or cephalosporin resistance in Shigella sonnei isolates from yaks. BMC Vet Res 14, 177 (2018).

Download citation


  • Shigella sonnei
  • MLST
  • PFGE
  • Virulence gene
  • Antimicrobial resistance genes