- Open Access
Low distribution of genes encoding virulence factors in Shigella flexneri serotypes 1b clinical isolates from eastern Chinese populations
Gut Pathogensvolume 9, Article number: 76 (2017)
The ability of Shigella to invade, colonize, and eventually kill host cells is influenced by many virulence factors. However, there is no analysis of related genes in Jiangsu Province of China so far. Shigella flexneri was collected from 13 cities of Jiangsu Province through the provincial Centers for Disease Control (CDC) for analysis of distribution of major virulence genes (ipaH, ipaBCD, ial, virF, virB, sigA, set1A, sepA, sat, pic, set1B and sen) detected by PCR technology.
A total of 545 isolates received were confirmed as S. flexneri which belongs to 11 serotypes of S. flexneri, among which serotype 2a was the most predominant (n = 223, 40.9%). All isolates were positive for ipaH gene, followed by sat (94.1%), sigA (78.9%), set1B (78.0%), pic (77.6%), set1A (74.5%), virF (64.8%), sepA (63.5%), sen (56.9%), ipaBCD (50.5%), ial (47.0%) and virB (47.0%). The presence of virulence genes in different serotypes was distinct. The existence of virulence genes of serotype 1b was generally lower than other serotype-the positive rate for virulence genes was between 0.0 and 14.1% except for ipaH and sat. In addition, virulence genes also fluctuated in different regions and at different times in Jiangsu province. The result of analysis on the relationship between virulence genes of S. flexneri showed that the existence of virulence genes of Shigella could be well represented by multiplex PCR combination ipaH + ial + set1A, which had a high clinical value.
The present study was designed to explore the prevalence of 12 S. flexneri-associated virulence genes. The data showed high diversity of virulence genes with regard to periods, regions and serotypes in Jiangsu Province of China.
Shigellosis is major health burden in many parts of the world. It is an acute invasive enteric infection caused by four members of Shigella species (S. dysenteriae, S. flexneri, S. boydii, and S. sonnei). Different serotypes for these species exist including more than 30 serotypes of S. flexneri which are categorized based on their O antigens, . Although the role of shigellosis in contributing to childhood mortality has been decreased significantly over the past few years, there are still about 28,000 children younger than 5 years of age who died of shigellosis every year . In a systematic review,  it was reported that due to low economic conditions and large population density in Asian countries, over 125 million Shigella-related infections led to 14,000 deaths per year. There are some factors contributing to the high prevalence of human Shigella infection. One reason for the high infection rate in some developing countries is the low sanitary conditions, knowing that Shigella spp. is transmitted via the fecal–oral route. Another important factor is that S. flexneri possess protective mechanisms that help it to survive even at high levels of acid in the stomach, which makes it highly infectious with only 10–100 microorganisms required to cause a disease .
In children, main symptoms of shigellosis vary from mild to severe which include: diarrhoea characterized by presence of blood in stool, abdominal cramping, fever, among other gastrointestinal complications. Its clinical phenotypes are determined by different virulence genes and the activity of immune system of the host. Among the many Shigella spp.—associated virulence factors, invasion plasmid antigen (ipa) B, C, D, and H as well as invasion-associated locus (ial) facilitates its penetration into intestinal cells . As with gram-negative bacteria, these genes are important for S. flexneri because they are components of the type III secretion system (T3SS) which is important for S. flexneri and other gram-negative pathogenic or symbiotic bacteria in manipulating the host cell processes and establish a successful infection . Shigella enterotoxin 1 (ShET-1), Shigella enterotoxin 2 (ShET-2) and shiga toxin (stx) are among virulence genes encoding Shigella enterotoxins. A group of genes mostly found in S. flexneri serotype 2 clinical samples encode ShET1, a 55 kDa protein complex [6, 7]. ShET2 has been reported in different species of Shigella . The stx is produced exclusively by S. dysenteriae 1, but this species is rare in China . The transcription of invasion-related genes is controlled by two proteins, virF and virB (InvE) which are derived from plasmids . Finally, Shigella spp. harbors toxic factors like serine protease autotransporters of enterobacteriaceae (SPATE) of which there are two phylogenetic classes . Shigella IgA-like protease homologue (sigA) and secreted autotransporter toxin (sat) belong to class 1 which are toxic to epithelial cells, while non-toxic SPATE class 2 toxins includes sepA, which facilitates intestinal inflammation and pic, a mucinase associated with colonization.
Although some studies reported the prevalence and distribution of S. flexneri virulence genes in some regions in China, investigations to dozen virulence genes of Shigella spp. mentioned above are still rare throughout the world, and to the best of our knowledge there is no report in China. To develop effective control strategies, it is important to conduct an epidemic study about Shigella in terms of its drug resistance and genetic features . For this reason, we sought to explore the distribution profile and prevalence of 12 Shigella-related virulence genes obtained from patients with diarrhea in Jiangsu Province of China, and discussed the genetic diversity and clinical applications of these genes.
Collection of bacteria isolates
CDC-based real-time surveillance program in 13 cities of Jiangsu Province from 2010 to 2015 (Fig. 1) were conducted by collecting suspecting Shigella spp. isolates from different patients with either diarrhea or dysentery in different hospitals in 13 cities by using routine biochemical techniques. Shigella is a class B infectious disease in China. The bacteria detected in any local hospital must be reported to the provincial CDC by the city’s CDC. The study was conducted in collaboration with the Provincial CDC, so the collection of Shigella was the most comprehensive.
Bacteria identification and serotyping
By use of Rapid ID32E strips (bioMérieux Corp., Singapore) and automated biochemical analyzer (Hitachi 917; Boehringer Mannheim, Japan), the collected samples were processed and screened. O and H antigens were examined using hyperimmune sera through slide agglutination test (Ningbo Tianrong Bio-pharmaceutical Company Limited), and thereafter, the serotypes were grouped according to the Kauffmann–White scheme.
Polymerase chain reaction (PCR) assay for virulence genes
Qiagen DNA mini kit was used for the extraction of DNA in line with the manufacturer’s instructions. The polymerase chain reaction (PCR) assays were performed targeting virulence genes using previously reported primers listed in Table 1. A reaction mixture, Green Taq Mix (Vazyme, Nanjing, China) was prepared as per manufacturer’s guidelines. Amplification was done in a thermocycler programmed with the following sequence: a 5-min initial denaturation at 95 °C, then 30 cycles including a 50 s denaturation at 95 °C, annealing for 45 s (annealing temperature is shown in Table 1), and 72 °Cfor 1 min and a single final extension at 72 °C for 7 min. For each virulence gene detected, a representative amplicon was sequenced to confirm that the gene was amplified by its specific primer.
Statistical analyses were performed using SPSS 16.0 database software. Distribution of different virulence genes in serotypes, periods and regions were analyzed by the Chi square test, and statistical differences between groups were considered to be significant for p < 0.05.
S. flexneri isolated during the study period
Altogether 545 isolates were confirmed as S. flexneri during the 6-year period of investigation. These 545 strains were distributed in 13 cities of Jiangsu province. The prevalence ranked as follows: Taizhou (n = 68, 12.5%), Suzhou (n = 57, 10.5%), and Lianyungang (n = 54, 9.9%). What’s more, only 20 (3.7%) strains were from Nanjing, the capital of Jiangsu Province (Fig. 1). In the year variation of S. flexneri, it was found that the most of the S. flexneri were isolated in 2012, and there was a trend of decrease in the following years (Table 2).
Serotypes of S. flexneri
All the 545 isolates of S. flexneri belonged to 11 serotypes of S. flexneri. Of these, S. flexneri serotype 2a was the most predominant (n = 223, 40.9%) compared with the other serotypes. 2a, 2b, 1a, 1b, x and 4c were the six most frequently isolated serotypes, accounting for 96.5% of all S. flexneri, while other serotypes accounted for less than 1.5%. Six major serotypes had an obvious fluctuation over time (Table 2). The other infrequently observed serotypes, including Y, 4, 4a, 3b and 6, were only found in a small amount within a certain period of one or 2 years.
Prevalence of virulence genes
The detection of invasion-associated genes in 545 S. flexneri showed that ipaH had the highest frequency (100%) followed by followed by ipaBCD (50.5%) and ial (47.0%) (Table 3). There were 253 (46.4%) strains with three invasive genes. The invasion-associated genes had serotype differences. The ial gene was detected in 66.7% of 18 S. flexneri 4c isolates, 60.1% of 223 S. flexneri 2a isolates, 50.9% of 106 S. flexneri 2b isolates, 32.9% of 85 S. flexneri 1a isolates and 56.7% of 30 S. flexneri x isolates. However, the ial gene was just found in 1.6% of 64 S. flexneri 1b. The ipaBCD gene was detected in 66.7% of S. flexneri 4c, 62.8% of S. flexneri 2a, 55.7% of S. flexneri 2b, 36.5% of S. flexneri 1a, 66.7% of S. flexneri x and none of the S. flexneri 1b.
A total of 353 (64.8%) isolates tested positive for virF genes while 256 (47.0%) isolates were positive for virB genes (Table 3). 254 (45.0%) bacteria simultaneously had two genes. In addition, virB gene existed in 66.7% of S. flexneri 4c and 66.7% of S. flexneri x, followed by S. flexneri 2a (60.1%), S. flexneri 2b (46.2%), S. flexneri 1a (34.1%) and S. flexneri 1b (0.0%). virF gene was found in 93.3% of S. flexneri x, followed by S. flexneri 4c (77.8%), S. flexneri 2a (75.8%), S. flexneri 2b (73.6%), S. flexneri 1a (52.9%), S. flexneri 1b (1.6%).
Of the 545 S. flexneri strains tested, 99.3% contained genes that encode SPATE proteins, such as class II (SepA, Pic) and/or class I (SigA, Sat) (Table 3). A total of 35.8% of the strains had only one or two class I SPATE genes and no class II SPATE. No strains had only class II SPATE, and 54.5% had both two class I SPATE genes and two class II SPATE genes. The most common SPATE among S. flexneri strains was sat (94.1% of strains), following by sigA (78.9%), pic (77.6%), sepA (63.5%). Regarding differences between serotypes, a total of 86.7% S. flexneri x had two class I SPATE genes and two class II SPATE genes, following by 72.2% S. flexneri 4c, 69.1% S. flexneri 2a, 57.5% S. flexneri 2b, 38.8% S. flexneri 1a, and 0.0% S. flexneri 1b.
The set1A gene was present in 406(74.5%) S. flexneri isolates (198 serotype 2a, 78 serotype 2b, 74 serotype 1a, 10 serotype 1b, 25 serotype x, 16 serotype 4c), and set1B was present in 425 (78.0%) S. flexneri isolates (206 serotype 2a, 86 serotype 2b, 72 serotype 1a, 6 serotype 1b, 28 serotype x, 17 serotype 4c) (Table 3). Both set1A and set1B were detected in 375 (68.8%) of S. flexneri isolates. The sen was present in 310 (56.9%) S. flexneri isolates (153 serotype 2a, 61 serotype 2b, 38 serotype 1a, 9 serotype 1b, 22 serotype x, 14 serotype 4c).
Fluctuation in time and place
The existence of these genes had a fluctuation over time and place. In general, there were two epidemic peaks in virulence genes in 2011 and 2014 (Fig. 2). Except for ShET-1, the positivity of virulence genes in S. flexneri was the lowest in 2012. There was no regular change in virulence genes between regions, such as the positive rate of invasion-associated genes was the highest in Yangzhou, followed by Lianyungang, but the highest positive rate of ShET-1 existed in Lianyungang, followed by Xuzhou (Fig. 2). When taking into account the different serotypes of year and regional changes, some interesting phenomena were noticed. The positive change in virulence genes of serotype 2a was consistent with the overall change in virulence genes, and serotype 2b had the highest existence of virulence gene in 2013. In addition, the number of virulence genes of serotype 1a isolated in 2012 was obviously smaller than that in other years. It should be noted, however, that the ShET-1 was generally independent of these changes (Additional file 1). Serotype 1a in Zhenjiang was a low virulence gene carrying type. In general, the variation of virulence genes among different serotypes was general not particularly obvious (Additional file 2).
Due to inadequate supply of quality water and low hygienic conditions in less developed countries, Shigella—a cause of inflammatory diarrhea and dysentery, poses major challenges to public health sectors. S. flexneri was the most common of the four species in many developing countries [19, 20]. However, in developed countries, S. sonnei is the commonest Shigella species isolated [21, 22]. The reason for this difference is unclear, however, it is apparent that efforts to boost sanitation and local hygiene have greatly decreased the prevalence of shigellosis and even changed the pattern in which Shigella species are most distributed. Jiangsu Province is located in the eastern part of China, with a population about 80 million. Epidemiological analysis of Shigella will be beneficial to the prevention and control of the infectious diseases in the region. The results of analysis of the distribution characteristics of S. flexneri in Jiangsu Province in the present study showed that S. flexneri 2a was the most common of the eleven serotype, which is different from the study conducted in Beijing in China reporting that S. flexneri 4c was the most prevalent serotype among 19 serotypes . In Jiangsu Province, serotype 4c accounted for only 3.3%. However, our results matched the findings in developing countries [19, 24] and Zhejiang Province of China . Even in Jiangsu Province, there were also differences between the various cities (Additional file 2). For example, most prevalent serotypes in Nanjing are serotypes 2b, serotypes 1b in Zhenjiang and serotypes 1a in Taizhou. What’s more, some rare serotypes were detected at specific times in specific cities, such as serotype 6 was only separated from Nanjing in 2010. High heterogeneity with regard to temporal distribution was noted in Shigella species and serotypes, which further suggested the need for serotype-level identification to enhance the effectiveness of control strategies.
Since the information on the variety of Shigella virulence genes in China is limited, to fully understand its pathogenicity, further research is required to advance the search for virulence-related genes for Shigella. In the present study, the prevalence and distribution of 12 such genes was examined. In the present study, ipaH gene was highly conserved in various serotypes. Similar findings have been shown in many other studies . The presence of many copies of this gene i.e. seven in chromosomes and five in plasmids may explain why the gene tested positive in all strains. Considering that this gene can be detected even after the loss of plasmid, it is promising target for diagnostic purposes. In Shigella, the ability to enter host cells depends on the availability of type III-secretion-system (T3SS) which are encoded by large virulence plasmids [26, 27]. ial gene has been identified in invasion processes and on inv plasmid . Many proteins form part of the T3SS complex which includes a needle-shaped oligomer that connects the inner and outer membrane of the bacteria. The oligomer contains invasive plasmid antigens ipaB, ipaC, and ipaD at its tip end [26,27,28,29], which can be identified using upstream region of ipaB, acting as marker. The effects of deleting ial and ipaBCD on invasiveness of S. flexneri are not known. Numerous studies have shown that there is a link between the ability of the Shigella spp. strains to cause diarrhea and the presence of invasive genes in the bacteria. Mokhtari et al.  showed that, unlike in asymptomatic patients, isolates from stools of patients with diarrhea contained invasive genes, ial and/or ipaBCD. A study by Phantouamath et al. , showed that ial gene was found only in isolates from cases. In our study, 47.0% S. flexneri’ isolates were positive for ial gene, and 50.5% S. flexneri’ isolates were positive for ipaBCD gene. Comparison with other similar studies, 78.9% S. flexneri’ isolates were positive for ial gene in Iran , and even 100% in Zhejiang of China . For the ipaBCD gene, our result is similar to that of a study in Peru (49%) , but lower than that of a study in Brazil (100%) . in this sense, the invasive ability of S. flexneri in Jiangsu Province is not strong compared with other areas. Moreover, prevalence of virulence genes showed obvious serotype characteristics, such as none S. flexneri 1b expressed both ial and ipaBCD strains. But it should be noted that the pathogenicity of S. flexneri is also related to both the number of infected bacteria and the immunity of infected people.
Expression of Shigella virulence genes is regulated by heat-stable nucleoid structural protein (H-NS) which downregulates their transcription during unfavorable conditions for invasion. In response to favorable environmental signals, transcription of a series of genes is activated starting from AraC-like protein gene virF, which subsequently turns on transcription of virB regulatory genes. Thereafter, virB protein reverses the H-NS-induced inhibition on transcription which eventually turns on the virulence genes on the plasmid [9, 34].In the present study, both virF and virB were found in 45.0% S. flexneri isolates, indicating that there might be other pathways for regulating gene expression. In addition, virF but not virB was found in 19.8% S. flexneri isolates, suggesting that virF regulated virulence genes not only through virB pathway. Interestingly, of the 545 S. flexneri, 11 strains had only virB, which may be due to loss of the virF gene. On the other hand, because of the importance of virF in regulating virulence genes, potential novel antibiotics targeting virF have gained increasing attention [35, 36]. However, only 64.8% of the positive rate of this gene might limit this antibiotics application.
Two new enterotoxins have recently been described in S. flexneri. One is called Shigella enterotoxin 1 (ShET-1), which is encoded in the set1 chromosomal gene. It has been suggested that in its active form, the ShET-1 toxin is composed of a subunit A (encoded by set1A) and five B subunits (encoded by set1B) . Other is plasmid-encoded ShET-2 (encoded by sen). ShET-1 and ShET-2 could alter electrolyte and water transport in the small intestine , which is closely related to the symptoms of dehydration in the shigellosis. Prior studies reported that set1 genes were only detected in S. flexneri serotype 2 (2a and 2b) isolates and less so in other serotypes. In contrast, in the current study, many S. flexneri serotypes tested positive for set1 genes [7, 12, 38]. In some serotypes, however, the prevalence of set1 (set1A and/or set1B) was significantly lower than in other serotypes, such as S. flexneri 1b, S. flexneri 3b (Table 3). And interestingly, 14.9% of S. flexneri had only one subunit of ShET-1, the question about whether a single subunit would affect the pathogenicity of ShET-1 remains to be answered, but which needs further study for verification. The association remains to be further studied. sen gene was found in 11 serotypes, with a majority between 40 and 80%, but the serotype 1b positive rate was only 14%. The low positive rate of ShET-1 and ShET-2 in S. flexneri 1b means that this serotype has a low ability to cause dehydration.
Another factor that possess virulence activities is the Serine protease autotransporters of Enterobacteriaceae (SPATEs), which are toxins secreted from gram-negative bacteria. Nevertheless, only a few studies have searched for the presence of their encoding genes in large Shigella collections. A similar study in Iran found that the sat gene was present in all S. flexneri isolates, and the presence of sigA, pic and sepA genes simultaneously were existed in 35.5% of S. flexneri . Comparing the similar study, unsurprising, the most common SPATEs among Shigella was sat in our study, but the positive rate of the other three genes of SPATEs was significantly higher than that of Iran. Interestingly, sat is now recognized as a pathogenic E. coli, although it was initially studied in uropathogenic E. coli strains. In comparison with previous studies on the frequencies of sat gene in E. coli [39, 40], however, the presence of sat gene in Shigella was found to be higher. It should be noted that except for sat gene, SPATEs of serotype 1b was significantly less than that of the other serotypes.
The virulence gene can be used to identify Shigella, which had been confirmed by previous studies. Some studies [41, 42] reported that the positive rate of detecting Shigella by a PCR assay targeting the ipaH gene was higher than that by the traditional culture method. The disadvantage of this method is that it can only identify one virulence gene at a time, though this disadvantage could probably be overcome by multiple PCR techniques by screening the amplified genes in view of the difficulty of multiple PCR and the restriction of the number of amplified genes. IpaH can be used as a marker gene of Shigella to detect the Shigella. Four genes (pic, set1A, set1B and sigA) are located on the chromosome SHI-1 Island, and the pic gene overlaps with set1A and set1B. When Shigella flexneri set1A gene was positive for Shigella flexneri, 94.1% Shigella set1B was positive, and 92.4% Shigella isolates were positive for pic and sigA. set1A positive Shigella had a stronger representation of the integrity of this segment of the gene. Because of the high expression of sat in Shigella, the clinical value of its amplification is not significant. Other virulence genes include ial, ipaBCD, virF, virB, sen and sepA, all of which are located on the large virulence plasmid (140 MDa). To reflect these virulence genes of Shigella, we chose the lowest existent ial gene as a marker and found that the positive rate of ial positive S. flexneri, ipaBCD was 98.8%, the positive rate of virF was 96.1%, the positive rate of virB was 92.6%, and the positive rate of sen and sepA was 94.5%. To sum up, multiplex PCR combination ipaH + set1A + ial can comprehensively reflect the virulence of Shigella.
In the present study, we provided some baseline information about the distribution of some virulence genes in clinical strains of S. flexneri in Jiangsu Province in China. It was found that the prevalence of these virulence genes varied greatly, leading to different severities of the disease. The profile of these virulence genes correlated with serotype, period and region. We found a low pathogenicity serotype (1b) and combination between those genes. These findings may help better control and identify Shigella strains.
- S. flexneri :
invasion plasmid antigen H
invasion associated locus
invasion plasmid antigen B
Shigella enterotoxin 1
Shigella enterotoxin 2
secreted autotransporter toxin
Shigella IgA-like protease homologue
Center for Disease Prevention and Control
Muthuirulandi Sethuvel DP, Devanga Ragupathi NK, Anandan S. Update on: Shigella new serogroups/serotypes and their antimicrobial resistance. Lett Appl Microbiol. 2017;64:8–18.
Lanata CF, Fischer-Walker CL, Olascoaga AC, Torres CX, Aryee MJ, Black RE. Global causes of diarrheal disease mortality in children < 5 years of age: a systematic review. PLoS ONE. 2013;8:e72788.
Bardhan P, Faruque AS, Naheed A. Decrease in shigellosis-related deaths without Shigella spp.—specific interventions, Asia. Emerg Infect Dis. 2010;16:1718–23.
Schroeder GN. Molecular pathogenesis of Shigella spp.: controlling host cell signaling, invasion, and death by type III secretion. Clin Microbiol Rev. 2008;21:134–56.
Lima IF, Havt A. Update on molecular epidemiology of Shigella infection. Curr Opin Gastroenterol. 2015;31:30–7.
Roy S, Thanasekaran K, Dutta Roy AR. Distribution of Shigella enterotoxin genes and secreted autotransporter toxin gene among diverse species and serotypes of shigella isolated from Andaman Islands, India. Trop Med Int Health. 2006;11:1694–8.
Niyogi SK, Vargas M. Prevalence of the sat, set and sen genes among diverse serotypes of Shigella flexneri strains isolated from patients with acute diarrhoea. Clin Microbiol Infect. 2004;10:574–6.
Wang XY, Tao F, Xiao D, Lee H, Deen J, Gong J, Zhao Y, Zhou W, Li W, Shen B, Song Y, Ma J, Li ZM, Wang Z, Su PY, Chang N. Trend and disease burden of bacillary dysentery in China (1991–2000). Bull World Health Organ. 2006;84:561–8.
Jost BH. Site of transcriptional activation of virB on the large plasmid of Shigella flexneri 2a by VirF, a member of the AraC family of transcriptional activators. Microb Pathog. 1993;14:481–8.
Parham NJ, Pollard SJ, Desvaux M, Scott-Tucker A, Liu C, Fivian A, Henderson IR. Distribution of the serine protease autotransporters of the Enterobacteriaceae among extraintestinal clinical isolates of Escherichia coli. J Clin Microbiol. 2005;43:4076–82.
Qin T, Bi R, Fan W, Kang H, Ma P, Gu B. Novel mutations in quinolone resistance-determining regions of gyrA, gyrB, parC and parE in Shigella flexneri clinical isolates from eastern Chinese populations between 2001 and 2011. Eur J Clin Microbiol Infect Dis. 2016;35:2037–45.
Sangeetha AV, Parija SC, Mandal I. Clinical and microbiological profiles of shigellosis in children. J Health Popul Nutr. 2014;32:580–6.
Talukder KA, Mondol AS, Islam MA, Islam Z, Dutta DK, Khajanchi BK, Azmi IJ, Hossain MA, Rahman M, Cheasty T, Cravioto A. A novel serovar of Shigella dysenteriae from patients with diarrhoea in Bangladesh. J Med Microbiol. 2007;56:654–8.
Müller D, Greune L, Heusipp G, Karch H, Fruth A, Tschäpe H, Schmidt MA. Identification of unconventional intestinal pathogenic Escherichia coli isolates expressing intermediate virulence factor profiles by using a novel single-step multiplex PCR. Appl Environ Microbiol. 2007;73:3380–90.
Vidal M, Kruger E, Durán C, Lagos R, Levine M, Prado V, Toro C, Vidal R. Single multiplex PCR assay to identify simultaneously the six categories of diarrheagenic Escherichia coli associated with enteric infections. J Clin Microbiol. 2005;43:5362–5.
Boisen N, Ruiz-Perez F, Scheutz F, Krogfelt KA, Nataro JP. Short report: high prevalence of serine protease autotransporter cytotoxins among strains of enteroaggregative Escherichia coli. Am J Trop Med Hyg. 2009;80:294–301.
Restieri C, Garriss G, Locas MC. Autotransporter-encoding sequences are phylogenetically distributed among Escherichia coli clinical isolates and reference strains. Appl Environ Microbiol. 2007;73:1553–62.
Faruque SM, Khan R, Kamruzzaman M, Yamasaki S, Ahmad QS, Azim T, Nair GB, Takeda Y, Sack DA. Isolation of Shigella dysenteriae type 1 and S. flexneri strains from surface waters in Bangladesh: comparative molecular analysis of environmental Shigella isolates versus clinical strains. Appl Environ Microbiol. 2002;68:3908–13.
Lluque A, Mosquito S, Gomes C, Riveros M, Durand D, Tilley DH, Bernal M, Prada A, Ochoa TJ, Ruiz J. Virulence factors and mechanisms of antimicrobial resistance in Shigella strains from periurban areas of Lima (Peru). Int J Med Microbiol. 2015;305:480–90.
Casabonne C, González A, Aquili V. Prevalence and Virulence Genes of Shigella spp. Isolated from Patients with diarrhea in Rosario, Argentina. Jpn J Infect Dis. 2016;69:477–81.
Gupta A, Polyak CS, Bishop RD, Sobel J, Mintz ED. Laboratory-confirmed shigellosis in the United States, 1989–2002: epidemiologic trends and patterns. Clin Infect Dis. 2004;38:1372–7.
Nógrády N, Király M, Borbás K, Tóth Á, Pászti J, Tóth I. Antimicrobial resistance and genetic characteristics of integron-carrier Shigella isolated in Hungary (1998–2008). J Med Microbiol. 2013;62:1545–51.
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:904–8.
Bangtrakulnonth A, Vieira AR, Lo Fo Wong DM, Pornreongwong S, Pulsrikarn C, Sawanpanyalert P, Hendriksen RS. Shigella from humans in Thailand during, to 2006: spatial-time trends in species and serotype distribution. Foodborne Pathog Dis. 1993;2008(5):773–84.
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.
Coster TS, Hoge CW, VanDeVerg LL, et al. Vaccination against shigellosis with attenuated Shigella flexneri 2a strain SC602. Infect Immun. 1999;67:3437–43.
Sansonetti PJ. Rupture, invasion and inflammatory destruction of the intestinal barrier by Shigella, making sense of prokaryote-eukaryote cross-talks. FEMS Microbiol Rev. 2001;25:3–14.
Nataro JP, Seriwatana J, Fasano A, Maneval DR, Guers LD, Noriega F, Dubovsky F, Levine MM, Morris JG. Identification and cloning of a novel plasmid-encoded enterotoxin of enteroinvasive Escherichia coli and Shigella strains. Infect Immun. 1995;63:4721–8.
Marteyn B, Gazi A. Shigella: a model of virulence regulation in vivo. Gut Microbes. 2012;3:104–20.
Mokhtari W, Nsaibia S, Majouri D, Ben Hassen A, Gharbi A, Aouni M. Detection and characterization of Shigella species isolated from food and human stool samples in Nabeul, Tunisia, by molecular methods and culture techniques. J Appl Microbiol. 2012;113:209–22.
Phantouamath B, Sithivong N, Insisiengmay S, Ichinose Y, Higa N, Song T, Iwanaga M. Pathogenicity of Shigella in healthy carriers: a study in Vientiane, Lao People’s Democratic Republic. Jpn J Infect Dis. 2005;58:232–4.
Hosseini Nave H, Mansouri S, Emaneini M. Distribution of genes encoding virulence factors and molecular analysis of Shigella spp. isolated from patients with diarrhea in Kerman, Iran. Microb Pathog. 2016;92:68–71.
da Cruz CB, de Souza MC, Serra PT, Santos I, Balieiro A, Pieri FA, Nogueira PA, Orlandi PP. Virulence factors associated with pediatric shigellosis in Brazilian Amazon. Biomed Res Int. 2014. https://doi.org/10.1155/2014/539697.
Gao X, Zou T, Mu Z, Qin B, Yang J, Waltersperger S, Wang M, Cui S, Jin Q. Structural insights into VirB-DNA complexes reveal mechanism of transcriptional activation of virulence genes. Nucleic Acids Res. 2013;41:10529–41.
Emanuele AA, Adams NE, Chen YC, Maurelli AT, Garcia GA. Potential novel antibiotics from HTS targeting the virulence-regulating transcription factor, VirF, from Shigella flexneri. J Antibiot (Tokyo). 2014;67:379–86.
Emanuele AA. Mechanism of action and initial, in vitro SAR of an inhibitor of the Shigella flexneri virulence regulator VirF. PLoS ONE. 2015;10:e0137410.
Ghosh S, Pazhani GP, Niyogi SK, Nataro JP, Ramamurthy T. Genetic characterization of Shigella spp. isolated from diarrhoeal and asymptomatic children. J Med Microbiol. 2014;63:903–10.
Vargas M, Gascon J, Jimenez De Anta MT. Prevalence of Shigella enterotoxins 1 and 2 among Shigella strains isolated from patients with traveler’s diarrhea. J Clin Microbiol. 1999;37:3608–11.
Andrade FB, Abreu AG, Nunes KO, Gomes TA, Piazza RM, Elias WP. Distribution of serine protease autotransporters of Enterobacteriaceae in typical and atypical enteroaggregative Escherichia coli. Infect Genet Evol. 2017;50:83–6.
Tapader R, Chatterjee S, Singh AK, Dayma P, Haldar S, Pal A, Basu S. The high prevalence of serine protease autotransporters of Enterobacteriaceae (SPATEs) in Escherichia coli causing neonatal septicemia. Eur J Clin Microbiol Infect Dis. 2014;33:2015–24.
Vu DT, Sethabutr O, Von Seidlein L, Tran VT, Do GC, Bui TC, Le HT, Lee H, Houng HS, Hale TL, Clemens JD, Mason C. Detection of Shigella by a PCR assay targeting the ipaH gene suggests increased prevalence of shigellosis in Nha Trang, Vietna. J Clin Microbiol. 2004;42:2031–5.
Wang SM, Ma JC, Hao ZY, Zhang ZY, Mason C, Sethabutr O, von Seidlein L, Wang XY, Xu ZY. Surveillance of shigellosis by real-time PCR suggests underestimation of shigellosis prevalence by culture-based methods in a population of rural China. J Infect. 2010;61:471–5.
WF, HQ, BG and PM conceived and designed the study; WS, XZ, SC and HQ contributed by collecting samples; WF and CY performed the experiments; WF, BG and PM contributed to analysis and interpretation; WF and CY drafted the manuscript. All authors read and approved the final manuscript.
We are very grateful to the provincial CDC and city CDC staff for their efforts to collect strains.
The authors declare that they have no competing interests.
Availability of data and materials
All data and material will be available upon requests to the corresponding author.
Consent for publication
Ethics approval and consent to participate
Ethical approval was obtained from the Ethics Committees of Xuzhou medical university. For this type of study, formal consent was not required.
This research was supported by the National Natural Science Foundation of China (81471994, 81702103), National Natural Science Foundation of China (81702103), Jiangsu Provincial Natural Science Foundation (BK20170252), Jiangsu Provincial Medical Talent (ZDRCA2016053), Six talent peaks project of Jiangsu Province (WSN-135), Advanced health talent of six-one project of Jiangsu Province (LGY2016042).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.