Genome sequence of Shigella flexneri strain SP1, a diarrheal isolate that encodes an extended-spectrum β-lactamase (ESBL)
- Ping Shen†1,
- Jianzhong Fan†2,
- Lihua Guo1,
- Jiahua Li3,
- Ang Li1,
- Jing Zhang1,
- Chaoqun Ying1,
- Jinru Ji1,
- Hao Xu1,
- Beiwen Zheng1Email author and
- Yonghong Xiao1
© The Author(s) 2017
Received: 9 February 2017
Accepted: 4 May 2017
Published: 12 May 2017
Shigellosis is the most common cause of gastrointestinal infections in developing countries. In China, the species most frequently responsible for shigellosis is Shigella flexneri. S. flexneri remains largely unexplored from a genomic standpoint and is still described using a vocabulary based on biochemical and serological properties. Moreover, increasing numbers of ESBL-producing Shigella strains have been isolated from clinical samples. Despite this, only a few cases of ESBL-producing Shigella have been described in China. Therefore, a better understanding of ESBL-producing Shigella from a genomic standpoint is required. In this study, a S. flexneri type 1a isolate SP1 harboring blaCTX-M-14, which was recovered from the patient with diarrhea, was subjected to whole genome sequencing.
The draft genome assembly of S. flexneri strain SP1 consisted of 4,592,345 bp with a G+C content of 50.46%. RAST analysis revealed the genome contained 4798 coding sequences (CDSs) and 100 RNA-encoding genes. We detected one incomplete prophage and six candidate CRISPR loci in the genome. In vitro antimicrobial susceptibility testing demonstrated that strain SP1 is resistant to ampicillin, amoxicillin/clavulanic acid, cefazolin, ceftriaxone and trimethoprim. In silico analysis detected genes mediating resistance to aminoglycosides, β-lactams, phenicol, tetracycline, sulphonamides, and trimethoprim. The bla CTX-M-14 gene was located on an IncFII2 plasmid. A series of virulence factors were identified in the genome.
In this study, we report the whole genome sequence of a blaCTX-M-14-encoding S. flexneri strain SP1. Dozens of resistance determinants were detected in the genome and may be responsible for the multidrug-resistance of this strain, although further confirmation studies are warranted. Numerous virulence factors identified in the strain suggest that isolate SP1 is potential pathogenic. The availability of the genome sequence and comparative analysis with other S. flexneri strains provides the basis to further address the evolution of drug resistance mechanisms and pathogenicity in S. flexneri.
Shigella species are a major causative cause of gastrointestinal infections throughout the world, especially in developing countries . Globally, there are approximately 164.7 million cases per annum, of which 1.1 million people are estimated to die from Shigella infections . Based on biochemical and serological properties, the genus Shigella comprises four serogroups: Shigella dysenteriae, Shigella flexneri, Shigella boydii and Shigella sonnei . S. flexneri is endemic in many developing countries and causes more deaths than any other Shigella serotypes . In China, shigellosis is the most common gastrointestinal infections  and most frequently isolated species responsible for shigellosis in mainland China is S. flexneri .
Antibiotic treatment is usually recommended for shigellosis as it reduces the duration and severity of symptoms, reduces the excretion of organisms and prevents potentially lethal complications . However, the emergence of extended-spectrum β-lactamase (ESBL)-producing S. flexneri is a major public health problem in China, as these strains are associated with critical infections . Most of the CTX-M, SHV, and TEM-type ESBLs genes are located on conjugative plasmids . Moreover, co-existence with other antimicrobial resistance genes is frequently observed in ESBL-producers , which makes the choice of effective treatment extremely limited. Increasing instances of ESBL-producing Shigella strains isolated from Asia have been reported [1, 2]. So far, only a few cases of ESBL-producing Shigella have been described in China . Therefore, a better understanding of ESBL-producing Shigella using a genomics approach is required.
The first genome sequence of S. flexneri was reported by Jin et al. in 2002 . So far, more than 145 S. flexneri strains have been sequenced and analyzed . To extend our understanding of the resistance mechanisms and pathogenesis of ESBL-producing S. flexneri, we performed sequencing and genomic analysis of the ESBL-harbouring S. flexneri SP1. Comparative genomics analysis of S. flexneri SP1 with other S. flexneri genomes may improve our understanding of the antibiotic resistance and virulence factors present in Shigella.
Shigella flexneri 1a isolate SP1 was isolated from the stool sample of a 76-year-old female diarrhea patient. Standard biochemical tests were performed using the Vitek II system (BioMerieux, France) and species-specific 16S rRNA sequencing was used to confirm the identity of isolate SP1. Serotyping was performed with specific antiserum (Denka-Seiken). Genomic DNA of SP1 was extracted from a single colony of the pure bacterial culture. Possible contamination with other DNA and misassemblies were assessed by performing a BLAST search against the non-redundant database as described previously . The whole genome of SP1 is in the expected size range for a Shigella genome and the coverage of the reads was consistent throughout the genome. The assembled draft genome sequence of SP1 was further verified by comparative analysis with the published complete genome sequences of S. flexneri strains.
Antimicrobial susceptibility testing
Susceptibility testing for ampicillin, amoxicillin/clavulanic acid, amikacin, aztreonam, piperacillin/tazobactam, cefazolin, cefoxitin, ciprofloxacin, ceftriaxone, cefepime, ertapenem, imipenem, gentamycin, tobramycin, levofloxacin, tigecycline, nitrofurantoin, and trimethoprim were performed by the microbroth dilution method and interpreted according to the CLSI guidelines . Escherichia coli strain ATCC 25922 was used as the control strain for susceptibility studies. Late log phase cells were harvested, and genomic DNA was extracted from the pure culture using the DNeasy Blood & Tissue kit (Qiagen, Germany) according to the manufacturer’s instructions. For the purpose of bacterial identification, we amplified the 16S rRNA gene with a 16S rRNA universal primer set and the PCR product was sequenced. Ethical approval was granted by the Ethics Committee of the First Affiliated Hospital of Zhejiang University.
Genome sequencing and assembly
The extracted DNA was visualized by agarose gel electrophoresis and quantitated by Qubit 2.0. Whole-genome sequencing was performed on the Illumina HiSeq 4000-PE150 platform. DNA was tailed, ligated to paired-end adaptors and PCR amplified with a 500 bp insert size and a mate-pair libraries with an insert size of 5 kb were used for library construction at the Beijing Novogene Bioinformatics Technology Co., Ltd. Illumina PCR adapter reads and low quality reads from the paired-end and mate-pair library were filtered by the quality control step using Novogene pipeline. All high quality paired reads were assembled using Velvet 1.2.10  into a number of scaffolds. The filtered reads were then passed handled by the next step of the gap-closing.
Genome annotation included the prediction of coding genes, transfer RNAs, ribosomal RNA, prophage, and clustered regularly interspaced short palindromic repeat sequences (CRISPR). Open reading frames (ORFs) were identified and classified using the Rapid Annotation using Subsystem Technology (RAST) server . Protein classification into functional groups was performed using the Clusters of Orthologous Groups of proteins (COGs) . Transfer RNAs and ribosomal RNA genes rRNAs were detected by tRNAscan-SE  and RNAmmer 1.2 software , respectively. PHASTER  was used to identify prophage and putative phage-like elements and CRISPRFinder  was used to identify CRISPR sites. The plasmid replicon was predicted by the PlasmidFinder Tool . ISfinder  was employed to search for IS sequences in the genome, with an e-value of 1E−3. plasmidSPAdes was used to produce plasmid sequences from the WGS data .
Antibiotic resistance genes prediction and virulence factors analysis
Antibiotic resistance genes were annotated using the comprehensive antibiotic resistance database (CARD)  and Resfinder  with default parameters. We further verified all putative antibiotic resistance genes (ARGs) through a BLAST search with cut-off e-value of 1E−0.5. Virulence factors were predicted by using BLAST to search against the VFDB database  with an e-value threshold of 1E−5 and also with VirulenceFinder 1.5 .
PlasmidFinder 1.3 was used for identify the incompatibility group of the plasmid present in S. flexneri SP1 . The plasmid sequence of bla CTX-M-14-harboring plasmid from isolate SP1 (named pSP1) was assembled with plasmidSPAdes . Assignment of the plasmid to an incompatibility (Inc) group was performed by multiplex PCR. PCRs were performed as described previously .
Phylogenetic analysis and comparative genomic analysis
Comparative genomic analysis was performed by orthology identification method as previously described [11, 28]. Genome sequences of the following representative S. flexneri strains were downloaded from the NCBI genome database: S. flexneri 2a strain 981 (CP012137), S. flexneri S7737 (AMJY00000000), S. flexneri 5a strain M90T (CM001474), S. flexneri CDC 796.83 (AERO00000000), S. flexneri 4343.70 (AFHC00000000), S. flexneri NCTC1 (LM651928), S. flexneri 2a strain 301 (AE005674), S. flexneri G1663 (CP007037), S. flexneri Shi06HN006 (CP004057), S. flexneri 2003036 (CP004056), S. flexneri str 4S BJ10610 (JMRK00000000), S. flexneri 4c strain 1205 (CP012140), S. flexneri 2002017 (CP001383) and S. flexneri 1a strain 0228 (CP012735). Phylogenetic reconstruction and analysis was performed wih the phangorn package, written in the R language . VennDiagram  was used to generate the Venn plots of S. flexneri SP1, S. flexneri str 4S BJ10610, S. flexneri 4c strain 1205 (CP012140), S. flexneri 2002017, and S. flexneri 1a strain 0228.
Results and discussion
Antimicrobial susceptibility profiles and antibiotic resistance genes
The in vitro antimicrobial susceptibility testing demonstrated that the strain SP1 was resistant to ampicillin, amoxicillin/clavulanic acid, cefazolin, ceftriaxone and trimethoprim, but susceptible to piperacillin/tazobactam, cefoxitin, cefepime, aztreonam, imipenem, amikacin, gentamicin, tobramycin, ciprofloxacin, levofloxacin, tigecycline and nitrofurantoin (Additional file 2: Table S2). We then screened the antibiotic resistance genes (ARGs) in the genome to further explore the genetic basis of multidrug resistance in this strain (Additional file 3: Table S3). In silico analysis revealed the presence of some putative ARGs for different drug classes. We detected genes mediating resistance to aminoglycosides (aadA24, strA and strB), β-lactams (bla CTX-M-14 and bla OXA-1), phenicol (catA1), tetracycline (tetD), sulphonamides (sul2), and trimethoprim (dfrA1). CTX-M-14 was the most frequent ESBL variant detected in Shigella isolates in China, followed by CTX-M-15 [4, 33]. Hitherto, only a few studies have reported the presence of ESBLs in S. flexneri . Interestingly, all of the bla CTX-M-14-harbouring S. flexneri strains were isolated from China [4, 7, 34, 35]. A previous study reported a high prevalence of extended-spectrum cephalosporin resistance Shigella mediated mainly by bla CTX-M (mainly bla CTX-M-14, 14.1%) in Hangzhou City, Zhejiang Province, China . Our study further indicates the existence of bla CTX-M-14-harbouring S. flexneri clone may responsible for these Shigella infections in this area. In addition, the insertion sequence (IS) elements were frequently detected in upstream of bla CTX-M genes . ISfinder was thus employed to scan the blaCTX-M-14 flanking sequences in a range of 6-kb for IS sequences and junction associated proteins. ISEcp1 and IS903B were found upstream of blaCTX-M-14. ISEcp1 belongs to the IS1380 family, which may enhance the expression of bla CTX-M-14/-18, bla CTX-M-17 and bla CTX-M-19 β-lactamase genes .
Genetic context of bla CTX-M-14 gene
Pathogenesis and virulence factors
Shigella flexneri remains a public health concern throughout the world and its pathogenesis should be further investigated. We performed a BLASTP search against the VRDB database and found several known virulence factors. These virulence factors include Shigella extracellular protein A (sepA), glutamate decarboxylase (gadA), invasion plasmid antigen (ipaH9.8), long polar fimbriae (lpfA), hexosyltransferase homolog (capU), invasion protein Shigella flexneri (ipaD), serine protease autotransporters of Enterobacteriaceae (pic), VirF transcriptional activator (virF), and Shigella IgA-like protease homologue (sigA). The ability to withstand an acid-challenge of pH 2.5 by S. flexneri is a necessary virulence trait, which requires acid-induction of a functional GdaA in the stationary-growth phase . IpaH9.8 is a member of the IpaH family of Shigella, which are encoded on the 220 kb virulence plasmid or chromosome and has been shown to be secreted into the host cell where it is targeted to the nucleus . Moreover, LpfA has been linked to virulence in enterohemorrhagic E. coli . Earlier reports indicated the SigA is cytopathic for HEp-2 cells and contributes to the intestinal fluid accumulation associated with S. flexneri infections . Translocation of effector proteins into host cells and the surrounding space is a common strategy used by S. flexneri to target signaling pathways in the host cell , and IpaD and VirF are required to facilitate bacterial invasion of host cells . More importantly, Pic is secreted by pathogenic Gram-negative bacteria through the autotransporter pathway and targets a broad range of human leukocyte adhesion proteins, which represent unique immune-modulating bacterial virulence factors . These data suggest that isolate SP1 is potential pathogenic, which is consistent with the isolation of SP1 from a diarrhea patient.
Comparative analysis with other S. flexneri strains
As more S. flexneri genomes have been sequenced in recent years, comparative genomic studies have progressed rapidly. So far, whole genome studies of S. flexneri have exclusively focused on the historical global spread and recent local persistence among these isolates. This work is the first description of the draft genome of a bla CTX-M-14-harbouring S. flexneri isolate and demonstrates the compares the genome of strain SP1 to other S. flexneri isolates. However, the data presented here is a preliminary report on the virulence profile and antibiotic resistance of S. flexneri strain SP1. Future studies involving more ESBL-encoding S. flexneri isolates from China are urgently needed to study the dynamics of the dissemination of ESBL genes, especially the complete sequence of plasmids carrying these ESBL genes.
BZ, PS and JF were involved in the overall experimental design. PS, JF, LG, JL, HX, JJ and CY performed microbiology and molecular biology experiments. BZ, LG and AL generated and analyzed the sequencing data. BZ participated in all discussions of data analysis and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
Whole Genome Shotgun project of Shigella flexneri strain SP1 has been deposited at DDBJ/EMBL/GenBank under the accession MTPK00000000. The version described in this paper is version MTPK00000000.
Ethics approval and consent to participate
This research was approved by the Research Ethics Committee of the First Affiliated Hospital, School of Medicine, Zhejiang University.
This work was supported by the National Natural Science Foundation of China (Nos. 81301460 and 81301461); National Key Research and Development Program of China (No. 2016YFD0501105).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Taneja N, Mewara A, Kumar A, Verma G, Sharma M. Cephalosporin-resistant Shigella flexneri over 9 years (2001–2009) in India. J Antimicrob Chemother. 2012;67:1347–53.View ArticlePubMedGoogle Scholar
- Tajbakhsh M, Garcia 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:1128–33.View ArticlePubMedGoogle Scholar
- Bennish ML, Wojtyniak BJ. Mortality due to shigellosis: community and hospital data. Rev Infect Dis. 1991;13(Suppl 4):S245–51.View ArticlePubMedGoogle Scholar
- Zhang R, Zhou HW, Cai JC, Zhang J, Chen GX, Nasu M, Xie XY. Serotypes and extended-spectrum beta-lactamase types of clinical isolates of Shigella spp. from the Zhejiang province of China. Diagn Microbiol Infect Dis. 2011;69:98–104.View ArticlePubMedGoogle Scholar
- Yang C, Li P, Zhang X, Ma Q, Cui X, Li H, Liu H, Wang J, Xie J, Wu F, et al. Molecular characterization and analysis of high-level multidrug-resistance of Shigella flexneri serotype 4 s strains from China. Sci Rep. 2016;6:29124.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu Y, Zhuang L, Kang H, Ma P, Xu T, Pan S, Gu B. Prevalence, resistance patterns, and characterization of integrons of Shigella flexneri isolated from Jiangsu Province in China, 2001–2011. Eur J Clin Microbiol Infect Dis. 2016;35:1347–53.View ArticlePubMedGoogle Scholar
- Zhang W, Luo Y, Li J, Lin L, Ma Y, Hu C, Jin S, Ran L, Cui S. Wide dissemination of multidrug-resistant Shigella isolates in China. J Antimicrob Chemother. 2011;66:2527–35.View ArticlePubMedGoogle Scholar
- Zheng B, Dong H, Xu H, Lv J, Zhang J, Jiang X, Du Y, Xiao Y, Li L. Coexistence of MCR-1 and NDM-1 in clinical Escherichia coli isolates. Clin Infect Dis. 2016;63:1393–5.View ArticlePubMedGoogle Scholar
- Jin Q, Yuan Z, Xu J, Wang Y, Shen Y, Lu W, Wang J, Liu H, Yang J, Yang F, et al. Genome sequence of Shigella flexneri 2a: insights into pathogenicity through comparison with genomes of Escherichia coli K12 and O157. Nucleic Acids Res. 2002;30:4432–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Ashton PM, Baker KS, Gentle A, Wooldridge DJ, Thomson NR, Dallman TJ, Jenkins C. Draft genome sequences of the type strains of Shigella flexneri held at Public Health England: comparison of classical phenotypic and novel molecular assays with whole genome sequence. Gut Pathog. 2014;6:7.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang T, Jiang X, Feng C, Li A, Dong H, Wu S, Zheng B. Whole genome sequencing uncovers a novel IND-16 metallo-beta-lactamase from an extensively drug-resistant Chryseobacterium indologenes strain J31. Gut Pathog. 2016;8:47.View ArticlePubMedPubMed CentralGoogle Scholar
- CLSI. Performance standards for antimicrobial susceptibility testing; 23rd informational supplement. Wayne: M100-S23: Clinical and Laboratory Standards Institute; 2013.Google Scholar
- Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, et al. The RAST server: rapid annotations using subsystems technology. BMC Genom. 2008;9:75.View ArticleGoogle Scholar
- Natale DA, Shankavaram UT, Galperin MY, Wolf YI, Aravind L, Koonin EV. Towards understanding the first genome sequence of a crenarchaeon by genome annotation using clusters of orthologous groups of proteins (COGs). Genome Biol. 2000;1:RESEARCH0009.View ArticlePubMedPubMed CentralGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:0955–64.View ArticleGoogle Scholar
- Lagesen K, Hallin P, Rødland EA, Stærfeldt H-H, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, Wishart DS. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44:W16–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35:W52–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Carattoli A, Zankari E, Garcia-Fernandez A, Voldby Larsen M, Lund O, Villa L, Moller Aarestrup F, Hasman H. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58:3895–903.View ArticlePubMedPubMed CentralGoogle Scholar
- Siguier P, Varani A, Perochon J, Chandler M. Exploring bacterial insertion sequences with ISfinder: objectives, uses, and future developments. Methods Mol Biol. 2012;859:91–103.View ArticlePubMedGoogle Scholar
- Antipov D, Hartwick N, Shen M, Raiko M, Lapidus A, Pevzner PA. plasmidSPAdes: assembling plasmids from whole genome sequencing data. Bioinformatics. 2016;32:3380–7.View ArticlePubMedGoogle Scholar
- Jia B, Raphenya AR, Alcock B, Waglechner N, Guo P, Tsang KK, Lago BA, Dave BM, Pereira S, Sharma AN, et al. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017;45:D566–73.View ArticlePubMedGoogle Scholar
- Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, Aarestrup FM, Larsen MV. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67:2640–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen L, Zheng D, Liu B, Yang J, Jin Q. VFDB 2016: hierarchical and refined dataset for big data analysis–10 years on. Nucleic Acids Res. 2016;44:D694–7.View ArticlePubMedGoogle Scholar
- Kleinheinz KA, Joensen KG, Larsen MV. Applying the ResFinder and VirulenceFinder web-services for easy identification of acquired antibiotic resistance and E. coli virulence genes in bacteriophage and prophage nucleotide sequences. Bacteriophage. 2014;4:e27943.View ArticlePubMedPubMed CentralGoogle Scholar
- Carattoli A, Bertini A, Villa L, Falbo V, Hopkins KL, Threlfall EJ: Identification of plasmids by PCR-based replicon typing. J Microbiol Methods 2005;63:219–28.View ArticlePubMedGoogle Scholar
- Zheng B, Li A, Jiang X, Hu X, Yao J, Zhao L, Ji J, Ye M, Xiao Y, Li L. Genome sequencing and genomic characterization of a tigecycline-resistant Klebsiella pneumoniae strain isolated from the bile samples of a cholangiocarcinoma patient. Gut Pathog. 2014;6:40.View ArticlePubMedPubMed CentralGoogle Scholar
- Schliep KP. phangorn: phylogenetic analysis in R. Bioinformatics. 2011;27:592–3.View ArticlePubMedGoogle Scholar
- Lam F, Lalansingh CM, Babaran HE, Wang Z, Prokopec SD, Fox NS, Boutros PC. VennDiagramWeb: a web application for the generation of highly customizable Venn and Euler diagrams. BMC Bioinform. 2016;17:401.View ArticleGoogle Scholar
- Rossi O, Baker KS, Phalipon A, Weill FX, Citiulo F, Sansonetti P, Gerke C, Thomson NR. Draft genomes of Shigella strains used by the STOPENTERICS consortium. Gut Pathog. 2015;7:14.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang N, Lan R, Sun Q, Wang J, Wang Y, Zhang J, Yu D, Hu W, Hu S, Dai H, et al. Genomic portrait of the evolution and epidemic spread of a recently emerged multidrug-resistant Shigella flexneri clone in China. J Clin Microbiol. 2014;52:1119–26.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:447–51.View ArticlePubMedGoogle Scholar
- 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:232–42.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Song C, Duan G, Zhu J, Yang H, Xi Y, Fan Q. Transposition of ISEcp1 modulates blaCTX-M-55-mediated Shigella flexneri resistance to cefalothin. Int J Antimicrob Agents. 2013;42:507–12.View ArticlePubMedGoogle Scholar
- Poirel L, Decousser JW, Nordmann P. Insertion sequence ISEcp1B is involved in expression and mobilization of a bla(CTX-M) beta-lactamase gene. Antimicrob Agents Chemother. 2003;47:2938–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Riccobono E, Di Pilato V, Villagran AL, Bartoloni A, Rossolini GM, Pallecchi L. Complete sequence of pV404, a novel IncI1 plasmid harbouring blaCTX-M-14 in an original genetic context. Int J Antimicrob Agents. 2014;44:374–6.View ArticlePubMedGoogle Scholar
- Valat C, Forest K, Billet M, Polizzi C, Saras E, Madec JY, Haenni M. Absence of co-localization between pathovar-associated virulence factors and extended-spectrum beta-lactamase (blaCTX-M) genes on a single plasmid. Vet Microbiol. 2016;192:163–6.View ArticlePubMedGoogle Scholar
- Bhagwat AA, Bhagwat M. Comparative analysis of transcriptional regulatory elements of glutamate-dependent acid-resistance systems of Shigella flexneri and Escherichia coli O157:H7. FEMS Microbiol Lett. 2004;234:139–47.View ArticlePubMedGoogle Scholar
- Edwards DJ, Streich FC Jr, Ronchi VP, Todaro DR, Haas AL. Convergent evolution in the assembly of polyubiquitin degradation signals by the Shigella flexneri IpaH9.8 ligase. J Biol Chem. 2014;289:34114–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Dogan B, Suzuki H, Herlekar D, Sartor RB, Campbell BJ, Roberts CL, Stewart K, Scherl EJ, Araz Y, Bitar PP, et al. Inflammation-associated adherent-invasive Escherichia coli are enriched in pathways for use of propanediol and iron and M-cell translocation. Inflamm Bowel Dis. 2014;20:1919–32.View ArticlePubMedGoogle Scholar
- Al-Hasani K, Navarro-Garcia F, Huerta J, Sakellaris H, Adler B. The immunogenic SigA enterotoxin of Shigella flexneri 2a binds to HEp-2 cells and induces fodrin redistribution in intoxicated epithelial cells. PLoS ONE. 2009;4:e8223.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang F, Jiang Z, Li Y, He X, Zhao J, Yang X, Zhu L, Yin Z, Li X, Wang X, et al. Shigella flexneri T3SS effector IpaH4.5 modulates the host inflammatory response via interaction with NF-kappaB p65 protein. Cell Microbiol. 2013;15:474–85.View ArticlePubMedGoogle Scholar
- Tobe T, Nagai S, Okada N, Adler B, Yoshikawa M, Sasakawa C. Temperature-regulated expression of invasion genes in Shigella flexneri is controlled through the transcriptional activation of the virB gene on the large plasmid. Mol Microbiol. 1991;5:887–93.View ArticlePubMedGoogle Scholar
- Ruiz-Perez F, Wahid R, Faherty CS, Kolappaswamy K, Rodriguez L, Santiago A, Murphy E, Cross A, Sztein MB, Nataro JP. Serine protease autotransporters from Shigella flexneri and pathogenic Escherichia coli target a broad range of leukocyte glycoproteins. Proc Natl Acad Sci USA. 2011;108:12881–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Li P, Yang C, Wang J, Yi S, Li H, Wang Y, Qiu S, Song H. Draft genome sequence of a new Shigella flexneri subserotype, 4S BJ10610. Genome Announc. 2014;2:e00715.PubMedPubMed CentralGoogle Scholar
- Connor TR, Barker CR, Baker KS, Weill FX, Talukder KA, Smith AM, Baker S, Gouali M, Pham Thanh D, Jahan Azmi I, et al. Species-wide whole genome sequencing reveals historical global spread and recent local persistence in Shigella flexneri. Elife. 2015;4:e07335.View ArticlePubMedPubMed CentralGoogle Scholar