Prevalence and antimicrobial resistance of Shigella flexneri serotype 2 variant in China

Shigella flexneri serotype 2 variant (II:3,4,7,8) was isolated in 2008 and first reported in China in 2013. In the present study, epidemiological surveillance from 2003 to 2013 in China suggested that this serotype first appeared in Guangxi in 2003; it then emerged in Shanghai and Xinjiang in 2004 and in Henan in 2008. Of the 1813 S. flexneri isolates, 58 S. flexneri serotype 2 variant strains were identified. Serotype 2 variant has emerged as a prominent serotype in recent years, with 2a (32.6%), X variant (25.2%), 1a (9.4%), X (6.3%), 2b (5.4%), and 1b (3.6%). According to phenotypic and genotypic analysis, the serotype 2 variant originated from 2a to 2b. A higher antibiotic resistance rate was observed between 2009 and 2013 than that between 2003 and 2008. Among 22 cephalosporin-resistant isolates, blaTEM-1, blaOXA-1, blaCTX-3, blaCTX-14, and blaCTX-79 were detected. Among 22 fluoroquinolone-resistant isolates, a Ser80Ile mutation in parC was present in all of the isolates. Moreover, 21 isolates had three gyrA point mutations (Ser83Leu, His211Tyr, Asp87Asn, or Gly) and one isolate had two gyrA point mutations (Ser83Leu and His211Tyr). The prevalence of His211Tyr in the fluoroquinolone-resistant isolates is concerning, and the mutation was first reported in China. Besides, 22 isolates harbored the aac(6′)-Ib-cr gene, and two isolates harbored qnrS1. In view of the increased epidemic frequency and multidrug-resistant strain emergence, continuous surveillance will be needed to understand the actual disease burden and provide guidance for shigellosis.


Introduction
Shigellosis, caused by Shigella species, is recognized as a major public health burden and continues to be an important cause of diarrheal diseases. Worldwide, Shigella episodes were estimated to be 164.7 million annually, of which, 163.2 million were in developing countries, resulting in 1.1 million deaths (Kotloff et al., 1999). Epidemics generally occur in underdeveloped and developing countries with poor sanitary conditions where transmission from person to person is common, or when food or water is contaminated by the organism (Qu et al., 2012). Despite economic and public health improvements, outbreaks of shigellosis are still reported regularly (Huang et al., 2005;Gaynor et al., 2009;Kuo et al., 2009;Qiu et al., 2011). In China, shigellosis has been ranked third in morbidity, and it has become the number one cause of disease-related death in children (Mathers et al., 2009).
Based on biochemical and serological properties, the genus Shigella is divided into four species or subgroups: S. flexneri, S. dysenteriae, S. boydii, and S. sonnei. Among these subgroups, the epidemic subgroup S. sonnei caused diarrhea in industrialized countries and S. flexneri in developing countries (Kotloff et al., 1999). A previous study showed that the annual shigellosis morbidity rate was 20.28 cases per 100,000 people in mainland China from the national surveillance data of 2009, and S. flexneri (67.3%) and S. sonnei (32.7%) were two major causative species (Sui et al., 2010). S. flexneri is divided into at least 20 serotypes (serotypes 1a,1b,1c,1d,2a,2b,2 variant,3a,3b,4a,4av,4b,5a,5b,X,Xv,Y,Yv,F6,and 7b) based on the combinations of antigenic determinants present on the O antigen of the cell envelope lipopolysaccharide (Qiu et al., 2013;Sun et al., 2013). Serotypes 1c, 4av, 7b, 1d, Yv, 2 variant and Xv are newly reported serotypes in recent years, and some have caused epidemic-level disease (Talukder et al., 2001;Stagg et al., 2008;Ye et al., 2010). Serotype 1c, first appeared in Bangladesh in 1989, was subsequently found to be prevalent in Bangladesh, Egypt, and Vietnam (Carlin et al., 1989;El-Gendy et al., 1999;Stagg et al., 2008). Serotype Xv was first identified in 2000 and later become the predominant serotype in Henan Province between 2002 and 2006. Serotype Xv was also the most prevalent serotype in Gansu and Anhui in 2007 (Ye et al., 2010). Therefore, the rapid expansion and spread of novel serotypes pose a severe threat to public health in areas where shigellosis is endemic.
The S. flexneri serotype 2 variant (II:3,4,7,8) was first reported by Qiu et al. in 2013(Qiu et al., 2013. During our routine surveillance of bacillary dysentery from 2003 to 2013, a total of 58 serotype 2 variant strains were identified. This novel serotype first appeared in Guangxi in 2003, and then emerged in other provinces. However, no study has yet extensively characterized the newly emerging strains of serotype 2 variant; therefore, an extensive study is needed to determine the emergence, antimicrobial resistance pattern, and epidemic trends of the serotype 2 variant.

Materials and Methods
Bacterial Isolates, Serotyping, and Biochemical Characterization Fresh stool samples from diarrhea patients with clinically suspected dysentery were collected in sentinel hospitals. Samples were cultured for Shigella by streaking directly onto Salmonella-Shigella agar and incubated at 37 • C for 18 h. Resultant Shigella colonies (colorless, semitransparent, smooth, and moist circular colonies) were routinely grown in a 37 • C incubator in Luria-Bertani agar plates (Qu et al., 2012). Then, the strains were submitted to our laboratory for further confirmation. This study was approved by the ethics committee of the Academy of Military Medical Sciences (China), and written approval was also obtained from the patients involved in this study. The isolates were confirmed using API 20E test strips (bioMerieux Vitek, Marcyl'Etoile, France) following the manufacturer's recommendations.
Serotypes of Shigella isolates were further determined with two serotyping kits: a commercially available kit (Denka Seiken, Tokyo, Japan) and monoclonal antibody reagents (Reagensia AB, Stockholm, Sweden).

Multilocus Sequence Typing (MLST)
MLST of these isolates was carried out using the protocols described at http://www.shigatox.net/ecmlst/cgi-bin/index. The PCR amplification conditions were as follows: 94 • C for 5 min; 33 cycles of 94 • C for 30 s, 55 • C for 30 s, and 72 • C for 1 min; and 72 • C for 10 min with Ex Taq DNA polymerase (Takara, Dalian, China). The PCR amplicons of 15 housekeeping genes were sequenced, and sequences were edited using SeqMan 7.0. Then, the sequences were uploaded to the EcMLST website for comparison, which allowed us to determine the sequence type.

Pulsed-Field Gel Electrophoresis (PFGE)
All isolates of S. flexneri were analyzed by PFGE according to the standard protocol for S. flexneri developed by the Centers for Disease Control and Prevention of the USA. Salmonella enterica serotype Braenderup H9812 was digested with XbaI and used as molecular weight standard. Slices of agarose plugs were digested with NotI (Takara, Dalian, China) at 37 • C for 3 h. Electrophoresis was carried out in 1% agarose SeaKem Gold gel (Lonza, Rockland, ME, USE) with the CHEF Mapper system (Bio-Rad) with the following run parameters: 6 V/cm and a linear increase in switching times from 2.16 to 54.17 s over a period of 20 h. The interpretation of the PFGE patterns was performed with BioNumerics software version 6.0 (Applied Maths, Sint-Martens-Latem, Belgium). A tree indicating relative genetic similarity was constructed based on the unweighted pair group method of averages and a position tolerance of 1.2%.

PCR Amplification of the Antibiotic-Resistant Determinants and Integrons
To understand the underlying mechanism conferring resistance to third-generation cephalosporins, PCR assays (Ahmed et al., 2006;Pan et al., 2006;Matar et al., 2007;Galani et al., 2010;Tariq et al., 2012) were used to detect the presence of β-lactamase genes such as bla SHV , bla TEM , bla OXA , and bla CTX−M . The variable regions of class 1 integrons and class 2 integrons were also amplified with primers listed in Table 1. The primers hep58 and hep59, hep74 and hep51 were used to detect the variable regions of class 1 integrons and class 2 integrons, respectively. The reverse primers aadA1, aadA2, aadA5, and cmlA1, together with the forward primer hep58, were used to amplify the gene cassettes of class 1 integron-positive strains. The forward primer hep74 and reverse primer hep51 were used to amplify the gene cassettes of class 2 integron-positive strains. To understand the underlying mechanism conferring resistance to quinolones, plasmid-mediated quinolone resistance determinants (PMQRs) such as qnrA, qnrB, qnrD, qnrS, and aac(6 ′ )-Ib-cr and quinolone resistance-determining regions (QRDRs) of the DNA gyrase (gyrA) and topoisomerase IV (parC) genes were amplified (Robicsek et al., 2006;Hu et al., 2007;Pu et al., 2009;Tariq et al., 2012). Purified PCR fragments were sequenced and analyzed by comparison with sequences in GenBank.

MLST and PFGE Analysis
To determine the phylogenetic relationship of these serotype 2 variant strains, an extended MLST scheme of 15 genes was performed. We found that all the isolates in this study belonged    (Figure 3), suggesting a high genetic diversity among the serotype 2 variant strains.

Molecular Analysis of Antibiotic-Resistance Determinants and Integrons
Twenty-two of 58 isolates showing resistance to cephalosporin were selected to test for antibiotic-resistance determinants and integrons, including the bla VIM , bla NDM , bla SHV , bla TEM , bla OXA , bla CTX−M , intI1, and intI2 gene regions ( Table 5). PCR screenings showed that all strains were negative for bla SHV , bla VIM , and bla NDM , but positive for bla TEM , bla OXA , and bla CTX−M . Sequencing results using primers amplifying the whole bla TEM gene showed 100% identity with the bla TEM−1 gene. All except two isolates from Shanghai and one isolate from Xinjiang harbored bla OXA−1 . Besides, 18 isolates contained the bla CTX−M gene; nine isolates from Shanghai, two from Xinjiang, and one from Henan harbored bla CTX−M−14 ; three isolates from Shanghai and one isolate from Henan harbored bla CTX−M−3 ; and one isolate from Xinjiang and one isolate from Henan  harbored bla CTX−M−79 . All isolates harbored class 1 integrons with bla OXA−30 and aadA1 gene cassettes, and one isolate from Shanghai harbored aacA4 and cmlA1 gene cassettes. All of the isolates harbored class 2 integrons following dfrA1, sat1, and aadA1 gene cassettes. In 22 quinolone-resistant isolates, no point mutations in QRDRs of gyrB and parE were found, but point mutations in QRDRs of gyrA and parC were identified in each resistant isolate. All of the norfloxacin-resistant isolates had the gyrA mutation of Ser83Leu and the parC mutations of Ser80Ile and His211Tyr in our study. Interestingly, the strain 2006040 (the only isolate resistant to levofloxacin) from Xinjiang and strain HN08096 from Henan had the gyrA mutation of Asp87Asn, while the rest of the strains except strain HN09103 had the gyrA mutation of Asp87Gly. Only two species of PMQR determinants [qnrS and aac(6 ′ )-Ib-cr] were identified. All of the isolates harbored the aac(6 ′ )-Ib-cr gene, and only two isolates from Shanghai harbored qnrS. Sequencing analysis determined that the two qnrS genes were qnrS1 (Table 6).

Discussion
The emergence of MDR strains has become a serious problem and has complicated the selection of empirical treatment for shigellosis. In the present study, the serotype 2 variant isolates were highly resistant to ampicillin, chloramphenicol, tetracycline, and ticarcillin, and 98.3% of these isolates showed MDR profiles.
One of the reasons for the rapid accumulation of resistance is excessive or inappropriate use of antibiotics in outpatients in China (Pickering, 2004;Zhang et al., 2011).
Our study further indicates that current resistance patterns have changed and that empirical therapy should keep pace with these changes. According to previous results, fluoroquinolones and third-generation cephalosporins are probably the best choice for empiric treatment of severe gastrointestinal infections caused by pathogenic bacteria (Gendrel et al., 2008). However, the resistance rates of these isolates to fluoroquinolones and thirdgeneration cephalosporins in the serotype 2 variant strains have increased at an alarming rate, and continuous surveillance of the resistance pattern will be essential to choose the appropriate antimicrobial therapy.
The increasing antibiotic resistance rate led us to study the genetics and mechanisms of antibiotic resistance. The horizontal transfer of integrons may account for the dissemination of resistance genes, and resistance to some antibiotics is associated with the presence of class 1 and class 2 integrons containing resistance gene cassettes (Ke et al., 2011). In the current study, bla OXA−30 and aadA1 were detected in the gene cassettes of class 1 integrons, and dfrA1+sat1+aadA1 was detected in the gene cassettes of class 2 integrons. In addition, of the 22 cephalosporinresistant isolates, 100% harbored the bla TEM−1 resistance gene, 81.8% harbored the bla CTX resistance gene, and 90.9% harbored the bla OXA resistance gene. Additionally, the strong ability to acquire resistance genes facilitated the survival of bacteria under the pressure of various antibiotics. The bla TEM−1 gene was detected in all of the ESBL-producing isolates; this gene exists at high frequencies in antibiotic-resistant bacteria across the globe and confers resistance to penicillin and other β-lactamase antibiotics (Barlow and Hall, 2003). OXA-type β-lactamases confer resistance to ampicillin and cephalothin, and they are characterized by high hydrolytic activity against oxacillin and cloxacillin (Bradford, 2001). Sequencing analysis showed that all of the bla OXA genes were bla OXA−1 , indicating the host specificity for this subtype in S. flexneri, which is consistent with the result of a previous study (Siu et al., 2000). Fluoroquinolone resistance in Enterobacteriaceae has been mostly attributed to modifications in the QRDRs of the gyrA and parC genes and/or PMQR determinants such as qnr and [aac(6 ′ )-Ib-cr] (Folster et al., 2011). In the current study, 95% of the quinolone-resistant Shigella isolates had at least three point mutations in gyrA (Ser83Leu, Asp87Gly/Asn, and His211Tyr) and parC (Ser80Ile). Furthermore, different amino acid substitutions at the same position may result in different quinolone susceptibility levels (Ruiz, 2003). The His211Tyr mutation in gyrA was detected in resistant S. flexneri serotype 2a strains from Bangladesh (Azmi et al., 2014). Notably, all of the quinolone-resistant isolates in our study had the gyrA mutation of His211Tyr. PMQR determinants are located on mobile genetic elements (i.e., plasmids), which may allow for dissemination among Shigella (Folster et al., 2011). The presence of the qnr gene can facilitate the selection of chromosomal mutations that cause quinolone resistance Martinez-Martinez et al., 2003). Worldwide, S. flexneri serotypes 1a, 2a, 2b, and 4c carrying the qnrS gene were reported with low incidence (Hata et al., 2005;Pu et al., 2009). In our study, only two strains from Shanghai carried the qnrS1 gene with a percentage of 3.4%, and one of the qnrS-positive strains was the only isolate that was resistant to levofloxacin. The aac(6 ′ )-Ib-cr gene is identified in many Enterobacteriaceae and is responsible for low-level resistance to fluoroquinolones (Frasson et al., 2011). The aac(6 ′ )-Ib-cr-positive S. flexneri 2a strain was first isolated in 1998 (Pu et al., 2009), and all of the 17 quinolone-resistant isolates were aac(6 ′ )-Ib-cr-positive, suggesting that the aac(6 ′ )-Ib-cr gene had been present in China for many years.
Serotype conversion is a major mechanism for S. flexneri to escape protective host immune responses, which may drive the emergence of novel and atypical serotypes under pressure (Allison and Verma, 2000). Based on phenotypic and genotypic analysis, the serotype 2 variant is likely to originate from serotypes 2a and 2b, which have experienced changes in their genome to adapt to altered environmental conditions. The serotype 2 variant was circulated as a major serotype during our 11-year study period, and the widespread dissemination of this serotype poses a serious threat to public health not only in China but also in the world.
In summary, studies on the serotype 2 variant have become very important owing to the increased epidemic frequency and the prevalence of MDR strains. In view of the fact that this serotype was prevalent in several provinces in China, it is likely to provoke a major crisis because MDR clones can and may have already spread to other countries. Therefore, continuous surveillance will be needed to determine the distribution and resistance development of this serotype so as to understand the actual disease burden and provide guidance for the clinical treatment of shigellosis.

Funding
This research was supported by the National Major Scientific and Technological Special Projects for Infectious Diseases during the Twelfth Five-year Plan Period (Nos. 2012ZX10004215 and 2013ZX10004607); the National Nature Science Foundation of China (Nos. 81473024, 81202252, 81373053 and 81371854); and the Beijing Science and Technology Nova program (No. xx2013061).