PLoS OnePLoS ONEplosplosonePLoS ONE1932-6203Public Library of ScienceSan Francisco, CA USA25617838430529610.1371/journal.pone.0116708PONE-D-14-42480Research ArticleGenetic Characterization of Shigella flexneri Isolates in Guizhou Province, ChinaGenetic Characterization of ShigellaLiShijun1*SunQiangzheng2WeiXiaoyu1KlenaJohn D.34WangJianping2LiuYing1TianKecheng1LuoXia2YeChangyun2XuJianguo2WangDingming1TangGuangpeng1Institute of Communicable Disease Control and Prevention, Guizhou Provincial Center for Disease Control and Prevention, 101 Bageyan Road, Guiyang, 550004, Guizhou, People’s Republic of ChinaState Key Laboratory for Infectious Disease Prevention and Control, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, National Institute for Communicable Disease Control and Prevention, China CDC, P.O. Box 5, Changping, Beijing, ChinaInternational Emerging Infections Program, US Centers for Disease Control and Prevention, Beijing, ChinaGlobal Disease Detection Branch, Division of Global Health Protection, Center for Global Health, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of AmericaChangYung-FuAcademic EditorCornell University, UNITED STATES

Competing Interests: The authors have declared that no competing interests exist.

Conceived and designed the experiments: SL QS JX CY. Performed the experiments: SL YL JW XW. Analyzed the data: SL XL. Contributed reagents/materials/analysis tools: QS JX CY KT GT DW. Wrote the paper: SL QS. Critically revised the manuscript: JK.

‡ These authors contributed equally to this work.

* E-mail: zjumedjun@163.com
24120152015101e01167082192014121220142015Li et alThis is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Shigella flexneri is one of the major etiologic causes of shigellosis in Guizhou Province, China. However, the genetic characteristics of circulating isolates are unknown. Phenotypic and molecular profiles of 60 S. flexneri isolates recovered in Guizhou between 1972 to 1982 and 2008 to 2010 were determined. Nine serotypes (1a, 2a, 3a, 1b, 2b, X, Y, 4av and Yv) were identified. Multi-locus sequence typing differentiated the isolates into 20 sequence types (STs); 18 were novel. Four STs, ST 129, ST 100, ST 126 and ST 18, were most abundant, accounting for 65% of the isolates. Thirty-nine NotI-pulsed field gel electrophoresis patterns (pulsotypes, PTs) were observed; eight PTs were represented by more than one isolate with six isolates sharing the PT 13 profile. Multi-locus variable-nucleotide tandem-repeat analysis recognized 44 different types (MTs); seven MTs were represented by more than one isolate and MT 1 was most commonly encountered. Correlation between genetic relationships and serotypes was observed among the isolates studied; the majority of isolates belonging to the same serotype from different years clustered together based on the molecular data. These clustered isolates were also from similar geographical origins. These results enhance our understanding of genetic relationships between S. flexneri in Guizhou Province and can be used to help understand the changing etiology of shigellosis in China.

This work was supported by the grant of national science and technology major project of China (No. 2009ZX10602-12), Talent Base Funds for Infectious Disease Control and Prevention in Guizhou Provincial Government (No. Qian Ren Ling Fa [2013] 15), National Natural Science Foundation of China (No. 81271788 and 81290345), National Basic Research Priorities Program (No. 2011CB504901) and National Key Program for Infectious Diseases of China (No. 2013ZX10004221, 013ZX10004216-001-002 and 2ZX10004215. We acknowledge the use of the EcMLST database which is operated by the Microbial Evolution Laboratory at Michigan State University. The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Data AvailabilityAll relevant data are within the paper.
Introduction

Shigellosis, gastroenteritis caused by Shigella spp., is a major public health problem in both developing and developed countries [1, 2]. Approximately 164.7 million cases of shigellosis occur annually worldwide, resulting in 1.1 million deaths, primarily among children aged < 5 years [3]. Shigellosis occurs mainly in developing countries due to poor hygiene and limited access to clean drinking water; in industrialized countries the disease mainly results from travel to developing countries and exposure to contaminated foods and/or food handlers [1]. In China, Shigella spp. is one of the most frequently isolated gastrointestinal pathogens [4], accounting for up to 1.7 million episodes of shigellosis annually, with up to 200,000 patients admitted to hospitals [5, 6].

Guizhou Province, with nearly 50 million people, is an under-developed province in the southwest of China. Shigellosis has been one of the primary bacterial diseases in Guizhou in past decades, and from 2007 to 2010, 48,222 cases of shigellosis were reported [7]. Four species of Shigella can cause shigellosis but Shigella flexneri is the predominant species in China. Although there has been an observed shift in prevalence from S. flexneri to Shigella sonnei in recent years, S. flexneri is still one of the major etiologic cause of shigellosis in Guizhou Province [7].

Although Shigella has been a major source of disease over the past decades, information on the genetic relationships of circulating S. flexneri isolates from Guizhou Province is lacking. In this study, the molecular techniques multi-locus sequence typing (MLST), pulsed field gel electrophoresis (PFGE) and multi-locus variable-nucleotide tandem-repeat analysis (MLVA) were used to analyze the relationships between S. flexneri isolates recovered from Guizhou during the periods 1972 to 1982 and 2008 to 2010.

Material and MethodsBacterial isolates and serotyping

Sixty isolates of S. flexneri, including 30 isolates recovered from 1972 to 1982 and 30 recovered from 2008 to 2010 in Guizhou Province, were analyzed (Table 1). The isolates were from Guiyang, Anshun, Qianxinan, Qiandongnan, Qiannan, Tongren and Zunyi Prefectures (seven of the nine prefectures making up Guizhou Province). All S. flexneri isolates were serotyped by slide agglutination using a commercially available monovalent antisera kit (Denka Seiken, Tokyo, Japan) and monoclonal antibody reagents (Reagensia AB, Sweden) per the manufacturer’s instructions [8]. S. flexneri isolates were routinly cultured in a 37°C incubator on Luria-Burtani (LB) agar plates or in an orbital shaker in LB broth.

10.1371/journal.pone.0116708.t001Isolation location, year and serotyping results of 60 <italic>S</italic>. <italic>flexneri</italic> isolates, Guizhou, 1972 to 1982 and 2008 to 2010.
No.Isolate No.PrefectureCountyYearSerotype
F011972GZ01TongrenSongtao19722a
F021972GZ02TongrenSongtao19722a
F031973GZ01GuiyangGuiyang1973y
F041973GZ03GuiyangGuiyang19734av
F051973GZ02GuiyangGuiyang19732a
F061973GZ04GuiyangGuiyang19733a
F071973GZ05GuiyangGuiyang19733a
F081973GZ06GuiyangGuiyang19733a
F091973GZ07GuiyangGuiyang19733a
F101973GZ08GuiyangGuiyang19732a
F111978GZ01QiandongnanKaili1978Yv
F121981GZ01ZunyiQiannan19811a
F131981GZ02ZunyiQiannan19811a
F141982GZ01QianxinanQinglong19823a
F151982GZ02QianxinanPuan19821a
F161982GZ03QianxinanCeheng19822a
F171982GZ04QianxinanZengfeng19821a
F181982GZ05GuiyangGuiyang19821a
F191982GZ06GuiyangGuiyang19821a
F201982GZ07GuiyangGuiyang19821a
F221982GZ09GuiyangGuiyang1982y
F231982GZ10ZunyiZunyi19821a
F241982GZ11GuiyangGuiyang19821a
F251982GZ12GuiyangGuiyang19821a
F261982GZ13ZunyiZunyi19821a
F271982GZ14ZunyiZunyi19821a
F281982GZ15ZunyiZunyi19821a
F291982GZ16ZunyiZunyi19821a
F301982GZ17ZunyiZunyi19821b
F311982GZ18ZunyiZunyi19821a
F322008GZ01GuiyangKaiyang20082a
F332008GZ02GuiyangKaiyang20082a
F342008GZ03AnshunZiyun20083a
F352008GZ07AnshunPingba20082a
F362008GZ08AnshunPingba20082a
F372008GZ11GuiyangKaiyang20083a
F382008GZ12GuiyangKaiyang20083a
F392008GZ13GuiyangKaiyang20083a
F402008GZ14GuiyangKaiyang20083a
F412008GZ15GuiyangKaiyang20083a
F422008GZ16GuiyangKaiyang20083a
F432008GZ17GuiyangKaiyang20083a
F442009GZ01GuiyangKaiyang20092a
F452009GZ02GuiyangKaiyang20092a
F462009GZ03AnshunPingba20092a
F472009GZ04AnshunPingba20092b
F482009GZ05AnshunPingba20092a
F492009GZ06AnshunPingba20092a
F502009GZ23AnshunZiyun20093a
F512009GZ25AnshunZiyun20092a
F522009GZ28GuiyangKaiyang20092a
F532009GZ29GuiyangKaiyang20092a
F542009GZ60GuiyangKaiyang20092a
F552009GZ68AnshunPingba20091a
F562010GZ01AnshunZiyun20102a
F572010GZ02AnshunZiyun20102a
F582010GZ03AnshunZiyun2010x
F592010GZ04AnshunZiyun2010x
F602010GZ05AnshunZiyun20102a
F612010GZ06AnshunZiyun20102a
Preparation of DNA

Genomic DNA for PCR was prepared directly from bacterial colonies by the lysis by boiling method [8]. Briefly, a single colony from an overnight culture at 37°C on LB agar was suspended in 30 μl of distilled water and boiled at 100°C for 10 min. The sample was immediately cooled on ice for 5 min and centrifuged at 13,000 × g at 4°C for 10 min. The supernatant, containing DNA, was used as the template for PCR amplification.

MLST

MLST analysis of 15 housekeeping genes was performed as described on the EcMLST website (http://www.shigatox.net/ecmlst). PCR products were sequenced bi-directionally.

Each unique allele was assigned a different number and the allelic profile (string of fifteen allelic loci) was used to define each isolate’s sequence type (ST). New allele numbers and STs were submitted to the EcMLST curator for confirmation and allocation of a unique identifier. Clustering and minimum spanning tree (MST) analysis was used to infer relationships among the isolates using the fingerprint analysis software BioNumerics version 4.5 (Applied Maths, Kortrijk, Belgium) [9].

PFGE

PFGE analysis was performed using the method described by Ye et al. [9] using the restriction enzyme NotI. PFGE images were analyzed using BioNumerics. A PFGE pulsotype (PT) was defined as a pattern with one or more DNA bands different from other patterns. A dendrogram constructed using PFGE patterns was generated using the UPGMA algorithm with Dice-predicted similarity value of two patterns set at 1.0% pattern optimization and 0.8% band position tolerance.

MLVA

MLVA was performed as previously described [10]. Eight VNTR loci (SF3, SF4, SF6, SF7, SF8, SF9, SF10 and SF25) were selected. The forward primer for each primer set was labeled at its 5′ end with a compatible HEX, FAM, TAMRA, and ROX dye, respectively. The loci were amplified individually, with each 20 μl PCR mixture containing 1 μl each primer, 1 μl DNA template, 10 μl 2 × Taq MasterMix (Cowin Biotech, Beijing, China) and deionized water used to make up volume differences to 20 μl. PCR products were analyzed by capillary electrophoresis on an ABI 3730XL sequencer with GeneScan 500 LIZ Size Standard (Applied Biosystems Incorporated, Carlsbad, CA, USA) as described [11]. The copy number of each VNTR locus was incorporated into BioNumerics software and analyzed as described previously [12]. Each unique allelic string was designated a unique MLVA type (MT). A dendrogram was constructed by UPGMA clustering based on categorical coefficient analysis.

ResultsDistribution of serotypes

The 60 S. flexneri isolates were grouped into nine serotypes (1a, 2a, 1b, 2b, 3a, X, Y, 4av and Yv) (Table 1). Three serotypes, 1a, 2a and 3a, were predominant. Serotype 1a was the most frequently identified serotype (50%, 15/30) among isolates from 1972 to 1982, however 2a was dominated from 2008 to 2010 (56.7%, 17/30) of isolates. Further, 93.8% (15/16) of the 1a isolates were from 1972 to1982, and 72.3% (17/22) of the 2a isolates were isolated during 2008 to 2010; 3a isolates were almost equally distributed across both time periods. Isolate expressing serotype 4av (1973GZ03) and Yv (1978GZ01), recently described by Sun et al [1315], were recovered from diarrheal cases in 1973 and 1978, respectively.

MLST based genotypes

The 60 isolates were divided into 20 STs, among which 2 STs (ST 18 and ST 100) have been previously reported; the remaining 18 STs (ST 120, ST 121, ST 124—ST 139) were unique (Table 2). The most common STs identified were ST 129 (22%), including isolates of serotype 1a; ST 100 (18%) including isolates of serotype 1a, 2a and X; ST 126 (13%) including all 3a isolates and ST 18 (12%) including isolates of serotype 1a, 2a, and Y. Among the most common STs, ST 18 (except for isolate 2008GZ02) and ST 129 isolates were recovered from 1972 to 1982; ST 100 isolates were only recovered from 2009 to 2010; and ST 126 isolates were recovered from both time periods, respectively. Eleven STs (18%) were singletons (Fig. 1). The predominant ST from 1972 to 1982 was ST 129 (43.3%, 13/30), while ST 100 was the predominant ST (36.7%, 11/30) from 2008 to 2010. A MLST cluster tree of the isolates showed they were divided into two clusters, designated A and B, with an overall coefficient of similarity of 50% (Fig. 1). ST 124 (4av, isolate 1973GZ03 isolated in 1973) was the only isolate within cluster B, while the remaining 19 STs formed cluster A. Cluster A was further divided in to subclusters A1 (15 isolates) and A2 (44 isolates); all the isolates in cluster A1 belonged to serotype 3a with the exception of one serotype 2a (2008GZ01) isolate. The STs in cluster A1 included ST 126, ST 127, ST 133 and ST 136–139. Cluster A2 was further divided into two distinct branches A2a (28 isolates) and A2b (16 isolates); A2a included 20 of the 22 S. flexneri 2a isolates, 2 1a, 2 X, and one each of 1b, 2b, Y and Yv. Branch A2b contained 14 of the 16 S. flexneri 1a isolates and a single 2a and Y isolate. The cluster tree indicated that isolates belonging to the same serotype closely clustered based on the time of isolation. A minimum spanning tree (MST), based on the 20 STs indicated that 20 STs were divided into 2 clonal complexes (CCs) (CC 18 and CC 126) and four singletons; CC18 contains isolates expressing serotypes 1a, 2a, 2b, 4b, X and Y, and included 10 STs. In contrast all isolates in CC 126 were serotype 3a and included five STs; ST 124, ST 132, ST 133 and ST 139 were singletons (Fig. 2).

10.1371/journal.pone.0116708.g001MLST clustering tree of <italic>S</italic>. <italic>flexneri</italic> isolates, Guizhou, from 1972 to 1982 and 2008 to 2010.

The 60 isolates from Guizhou province were analyzed using a 15 allele MLST as described in the Materials and Methods.

10.1371/journal.pone.0116708.g002Genetic relationships of <italic>S</italic>. <italic>flexneri</italic> isolates recovered from Guizhou Province based on MLST.

The minimum spanning tree was constructed using the 20 identified STs obtained from the 60 Guizhou Province isolates. Each circle corresponds to a single ST. The shadow zones in different color correspond to different clonal complexes. The size of the circle is proportional to the number of the isolates, and the color within the cycles represents the serotypes of the isolates. The corresponding color, serotype, number of isolates and back ground information are shown on the right of the minimum spanning tree.

10.1371/journal.pone.0116708.t002<italic>S</italic>. <italic>flexneri</italic> MLST allelic profiles and ST designation, Guizhou, 1972 to 1982 and 2008 to 2010.
IsolatesNameAlle profileST
arcAaroEaspCclpXcyaAdNaGfadDgrpEicdAlysPmdhmtlDmutSrpoSuidA
1972GZ018101316101114101611191514151418
1972GZ028101316101114101611191514151418
1973GZ018101316101114101611191514151418
1973GZ0283315916331114101611179153715236124
1973GZ03810131610112301016111915141514125
1973GZ04835151610111430301119315147314126
1973GZ058351516101114302071119315147314127
1973GZ06835151610111430301119315147314126
1973GZ07835151610111430301119315147314126
1973GZ088101316101114101611191514151418
1978GZ-018101516101122910161271915141514128
1981GZ012810131610112301016111915141514129
1981GZ022810131610112301016111915141514129
1982GZ01835151610111430301119315147314126
1982GZ022810131610112301016111915141514129
1982GZ038101316101114101611191514151418
1982GZ042810131610112301016111915141514129
1982GZ052810131610112301016111915141514129
1982GZ062810131610112301016111915141514129
1982GZ072810131610112301016111915141514129
1982GZ0928101318910112301016111915141514130
1982GZ102810131610112301016111915141514129
1982GZ112810131610112301016111915147214131
1982GZ128101316101114101611191514151418
1982GZ132810131610112301016111915141514129
1982GZ142810131610112301016111915141514129
1982GZ152810131610112301016111915141514129
1982GZ162810131610112301016111915141514129
1982GZ17834131891011141016111915141514132
1982GZ182810131610112301016111915141514129
2008GZ01810151610111410207111915147314133
2008GZ028101316101114101611191514151418
2008GZ038351516101114302071119315147314127
2008GZ07810171161011141016111915141514134
2008GZ08810171161011141016231915141514135
2008GZ11835151610111430302319315147314136
2008GZ12835151610111430301119315147314126
2008GZ13835131610111430301119315147314137
2008GZ14835151610111430301119315177314138
2008GZ15835151610111430301119315147314126
2008GZ16835151610111430301119315147314126
2008GZ1783513161011143030111915141514139
2009GZ0181015161011141016231915141514120
2009GZ0281013161011141016231915141514100
2009GZ03810171161011141016231915141514135
2009GZ0481015161011141016231915141514120
2009GZ05810171161011141016231915141514135
2009GZ06810171161011141016231915141514135
2009GZ23835151610111430301119315147314126
2009GZ2581013161011141016231915141514100
2009GZ2881013161011141016231915141514100
2009GZ2981013161011141016231915141514100
2009GZ6081013161011141016231915141514100
2009GZ6881013161011141016231915141514100
2010GZ0181013161011141016231915141514100
2010GZ0281013161011141016231915141514100
2010GZ0381013161011141016231915141514100
2010GZ0481013161011141016231915141714121
2010GZ0581013161011141016231915141514100
2010GZ0681013161011141016231915141514100
PFGE based Genotypes

The genotypes and genetic relatedness of the Guizhou S. flexneri isolates were also determined using PFGE. NotI-digested S. flexneri DNA generated 39 reproducible unique PTs, each with 12–17 bands. Eight patterns were represented by more than one isolate with PT 13 containing the greatest number of isolates, followed by PT 16. Among the isolates from 1972 to 1982, the predominant PT was PT 13 or PT 16, each representing 20.0% (6/30) of the total, while the predominant PT among isolates from 2008 to 2010 was PT 06, PT 19 or PT 31, each representing 13.3% (4/30) of the total. All 60 isolates were related at a coefficient of similarity of 60%, but two main clusters could be distinguished at a 62% similarity value (cluster A and B; Fig. 3). Cluster A was split into two additional broad subgroups, A1 (15 isolates) and A2 (1 isolate). The majority (12 of 14) serotype 3a isolates, with the exception of isolates 1973GZ04 and 1982GZ01, grouped together in A1; single isolates of serotype 2a, 4av and Y complete this group. Subgroup A2 contained a single Y serotype isolate. Cluster B split into two subgroups as well, subgroup B1 (17 isolates) and B2 (27 isolates). The majority of serotype 1a isolates (14 of 16) were found in subgroup B1; the three remaining isolates expressed 1b, 3a and Y. Subgroup 2B contained 21 of the 22 serotype 2a isolates, 2 isolates of 1a and X and single isolates of 3a and Yv. Similar to MLST, the majority of isolates expressing the same serotype were closely clustered and this was also related to the year of isolation. Additionally, the isolate (isolate 1973GZ03 recovered in 1973) of serotype 4av, was formed cluster A3, while the serotype Yv isolate, shared a unique but similar PT and was included in subgroup B2.

10.1371/journal.pone.0116708.g003Relationship of <italic>S</italic>. <italic>flexneri</italic> isolates recovered from Guizhou based on <italic>Not</italic>I-PFGE analysis.

The dendrogram were constructed using UPGMA. The corresponding PFGE pattern, serotype and background information are shown to the right of the dendrogram.

MLVA typing

Using MLVA, the 60 S. flexneri isolates grouped to 44 different MTs (Fig. 4). Seven MTs were represented by more than one isolate with MT 1 occurring most frequently (n = 7) followed by MT 42 (n = 5). The predominant MT among isolates recovered from 2008 to 2010 was MT 1 (23.3%; 7/30), however the predominant MT among isolates from 1972 to 1982 was not obvious. MLVA analysis showed the greatest diversity among the 60 S. flexneri isolates resulting in an extensively branched tree. Like MLST and PFGE, two main clusters, A and B, were observed (Fig. 4). The majority of serotype 1a and 2a isolates were assigned in cluster A, in which the 1a and 2a isolates were further grouped in cluster A1 and A2, respectively, whereas most of the serotype 3a isolates were grouped in cluster B. Isolates of serotype X assigned to cluster A showed relatively close relationship to serotype 2a, while the serotype Y isolates were closely related to serotype 1a or 3a isolates. Isolates belonging to the same serotype but recovered from different years showed clear relatedness, indicated by grouping in the same clusters. For example, serotype 3a isolates recovered from 1972 to 1982 and 2008 to 2010 clustered together, and similar characteristics were also observed for isolates belonging to serotype 1a and 2a (Fig. 4). Serotype 4av and Yv isolates were closely grouped in cluster B1. Additionally, isolates belonging to the same serotype and within a close time span clustered together based on geographical origin. For example, 2a isolates recovered in Guiyang, Ziyun and Pingba shared very similar patterns.

10.1371/journal.pone.0116708.g004Relationship of <italic>S</italic>. <italic>flexneri</italic> isolates recovered from Guizhou based on MLVA.

Isolates were analyzed using an eight VNTR loci MLVA scheme. The dendrogram was constructed using UPGMA. The corresponding MLVA type with copy numbers for the eight VNTRs, serotype, and background information are shown to the right of the dendrogram.

Discussion

In the present study, the genetic characteristics of 60 isolates of S. flexneri recovered from Guizhou Province between 1972 to 1982 and 2008 to 2010 were systematically studied. Thirty of the S. flexneri isolates were from shigellosis cases from 1972 to 1982 and the remaining isolates were recovered from patients with shigellosis from 2008 to 2010. The serotypes of the isolates used in this study included 1a, 2a, 3a 1b, 2b, X, Y, 4av and Yv; isolates were also from seven of the nine prefectures in Guizhou Province. Serotyping results indicated that S. flexneri 1a and 2a were the predominant serotypes from 1972 to 1982 (94%) and 2008 to 2010 (72%), respectively, while 3a isolates were almost equally recovered in both periods. The predominant serotype (2a) in Guizhou Province recovered from 2008 to 2010 are consistent with isolates from Suzhou of Jiangsu, Henan and Shanxi Province [1618], but are different from isolates of other provinces such as Beijing (4a and 4b) and Jiading of Shanghai Province [19, 20]. It is noteworthy that isolates of serotype 4av and Yv were recovered as early as 1973 and 1978, respectively, indicating the early emergence of these serotypes in China.

Recently a number of genotyping methods with higher discriminatory power than serotyping or biochemical testing such as MLST [21, 22], PFGE [5, 23] and MLVA [10, 11, 24] were introduced to characterize Shigella isolates. These methods are based primarily based on changes in isolate genotype, permitting analysis of phylogenetic relationships. Analysis of the isolates can be helpful for clinical diagnosis, treatment, prevention and control of shigellosis. Choi et al. [22] showed that S. flexneri serotypes 1–5, X and Y clustered together in a group containing many allelic variants while serotype 6 formed a distinct group, as previously established [25, 26]. Wang et al.[10] reported that phylogenetic groupings of 242 S. flexneri isolates recovered from shigellosis cases in Taiwan between 1995 to 2008, based on PFGE and MLVA profiles, correlated with serotype and isolate origin. Two distinct clusters for isolates of serotype 3 were shown but only one distinct cluster for each of the serotype groups 1a/1b/NT, 2a/2b/X/NT, 4a/Y, and 6 were revealed. Serologically different isolates including serotype Y and subserotype 4a; serotype X and subserotype 2b; subserotypes 1a and 1b, and subserotypes 3a and 3b, were genetically more closely related than indicated by serotyping alone.

Ye et al. [9] previously analyzed 37 serotype X and 69 serotype 1a, 2a, 2b, 3a, 4a, 5b, and Y isolates from China; all belonged to ST91 (later renamed ST 100), and concluded that S. flexneri epidemics in China have been caused by a single epidemic clone, ST 100. In this study, 60 isolates of S. flexneri from Guizhou Province separated into 20 STs based on a 15 loci MLST scheme; 18 of the STs were novel. The most common STs from 1972 to 1982 were ST 18, ST 126 and ST 129, however, ST 100 and ST 126 appeared between 2008 to 2010. Our results suggested that the predominant ST was ST 129 from 1972 to 1982, while ST 100 was the the predominant ST during 2008 to 2010, and the predominant ST is consistent with ST of isolates from other provinces of China [9]. Isolates belonging to the same serotype clustered in accordance with the year of isolation using all three genotyping approaches. MST indicated that the 20 STs were divided into 2 CCs, CC 18 and CC 126, and 4 singletons. CC 18 contained isolates expressing serotypes 1a, 2a, 2b, 4b, X and Y, while all the isolates in CC 126 belonged to serotype 3a. In addition, both the cluster tree and the MST, based on MLST data, showed that the isolates of serotype 4av (1973GZ03) was distant from the isolates belonging to other serotypes.

PFGE is a broadly applicable typing method with a high degree of intra- and inter laboratory reproducibility when standardized protocols are followed [23]. It has been shown to be a powerful tool in the laboratory for discriminating Shigella isolates during an outbreak [27]. In this study, PFGE discriminated the 60 isolates of S. flexneri into 39 unique PFGE patterns. Isolates belonging to the same serotypes mainly clustered together; the most closely related isolates were temporally associated with one another was well, suggesting that some drift was associated within each serotype over time. For instance, isolates of serotype 2a in cluster B2 isolated during the period of 1972 to 1982 and 2008 to 2010 were closely clustered, respectively, and similar clustering characteristics were observed for the isolates belonging to serotype 1a (cluster B1) and 3a (cluster A1), but some isolates, such as 2008GZ03, 1982GZ03 and 2009GZ23, were assigned irrespective their isolation time. Additionally, isolates with similar geographic origin were also often grouped by PFGE together as they tended to express the same serotypes.

MLVA is a prominent typing tool which has been used for characterizing S. flexneri; it has also been a useful tool for phylogenetic analysis [10]. In the present study, the S. flexneri isolates were discriminated into 44 different MTs and showed a low (approximately 20%) coefficient of similarity, indicating the high discriminatory power of the MLVA method. This finding is consistent with a previous study showing MLVA exhibited a discriminatory power greater than PFGE [10]. For most of isolates belonging to serotypes 1a, 2a and 3a, MLVA results correlated with serotyping. However, isolates of serotype X and Y were associated with serotype 1a and 2a isolates, respectively, and the serotype 4av isolate was closely related to serotype 3a isolates; this is similar to results observed previously [10]. Similar to PFGE, the majority of isolates belonging to the same serotype were temporally and geographically related.

In this study, isolates 1973GZ03 and 1978GZ01were serologically identified as serotype 4av and Yv respectively. These serotypes differ from 4a and Y because they react with monoclonal antibody MASF IV-1 [1315]. In theory, serotype 4av and Yv should have originated from serotype 4a and Y, primarily differing only in the acquisition of a 6.8k plasmid carrying a phospheantransferase gene (opt) responding for the MASF IV-1 antigenic determinant in these isolates [1315]. In this study, isolates of serotype of 4av (1973GZ03) and Yv (1978GZ01) were typed as ST 124 and ST 128 with MLST, respectively. MLST showed that 4av isolate were genetically distinct from isolates belonging to other serotypes, while Yv isolates were relatively close to isolates of serotype 2a. Molecular analysis indicates that isolate 1978GZ01 carries a dysfunctional gtrII gene within genome, and hence it is genetically similar to isolates from serotype 2a (unpublished data). This would explain the observed similarity between isolate 1978GZ01 and the majority of serotype 2a isolates.

Conclusions

In conclusion, phenotypic and molecular profiles of 60 S. flexneri isolates recovered in Guizhou between 1972 to 1982 and 2008 to 2010 were analysed. Nine serotypes (1a, 2a, 3a, 1b, 2b, X, Y, 4av and Yv) were identified, and the predominant serotype has changed from 1a to 2a in Guizhou Province. MLST differentiated the isolates into 20 sequence types (STs); 18 were novel. Four STs, ST 129, ST 100, ST 126 and ST 18, were most abundant, accounting for 65% of the isolates. The predominant ST was ST 129 in 1972 to 1982, while ST 100 was the predominant ST during 2008 to 2010. Thirty-nine NotI-PFGE (pulsotypes, PTs) were observed; eight PTs were represented by more than one isolate with six isolates sharing the PT 13 profile. MLVA analysis recognized 44 different types (MTs); seven MTs were represented by more than one isolate and MT 1 was most commonly encountered. Correlation between genetic relationships and serotypes was observed among the isolates studied; the majority of isolates belonging to the same serotype from different years clustered together based on the molecular data. These clustered isolates were also from similar geographical origins.

ReferencesNiyogiSK (2005) Shigellosis. J Microbiol43: 133143 1588008815880088IzumiyaH, TadaY, ItoK, Morita-IshiharaT, OhnishiM, et al (2009) Characterization of Shigella sonnei isolates from travel-associated cases in Japan. Journal of medical microbiology58: 14861491 doi: 10.1099/jmm.0.011809-01958990019589900von SeidleinL, KimDR, AliM, LeeH, WangX, et al (2006) A multicentre study of Shigella diarrhoea in six Asian countries: disease burden, clinical manifestations, and microbiology. PLoS medicine3: e353 1696812416968124QiuS, WangY, XuX, LiP, HaoR, et al (2013) Multidrug-resistant atypical variants of Shigella flexneri in China. Emerging infectious diseases19: 11471150 doi: 10.3201/eid1907.1112212376375423763754XiaS, XuB, HuangL, ZhaoJY, RanL, et al (2011) Prevalence and characterization of human Shigella infections in Henan Province, China, in 2006. Journal of clinical microbiology49: 232242 doi: 10.1128/JCM.01508-102106829121068291WangXY, TaoF, XiaoD, LeeH, DeenJ, et al (2006) Trend and disease burden of bacillary dysentery in China (1991–2000). Bulletin of the World Health Organization84: 561568 1687823016878230WeiX, TianK, YouL, MaQ, LiuY, et al (2012) Epidemic Characteristics and Etiological Analysis of Bacillary Ib’sentery in Guizhou Province During the Period of 2007–2010. Practical Preventive Medicine19: 1185.SunQ, LanR, WangY, ZhaoA, ZhangS, et al (2011) Development of a multiplex PCR assay targeting O-antigen modification genes for molecular serotyping of Shigella flexneri. Journal of clinical microbiology49: 37663770 doi: 10.1128/JCM.01259-112188097421880974YeC, LanR, XiaS, ZhangJ, SunQ, et al (2010) Emergence of a new multidrug-resistant serotype X variant in an epidemic clone of Shigella flexneri. Journal of clinical microbiology48: 419426 doi: 10.1128/JCM.00614-091995527319955273WangYW, WatanabeH, PhungDC, TungSK, LeeYS, et al (2009) Multilocus variable-number tandem repeat analysis for molecular typing and phylogenetic analysis of Shigella flexneri. BMC microbiology9: 278 doi: 10.1186/1471-2180-9-2782004211920042119ChiouCS, WatanabeH, WangYW, WangWL, TerajimaJ, et al (2009) Utility of multilocus variable-number tandem-repeat analysis as a molecular tool for phylogenetic analysis of Shigella sonnei. Journal of clinical microbiology47: 11491154 doi: 10.1128/JCM.01607-081922509519225095Hyytia-TreesE, SmoleSC, FieldsPA, SwaminathanB, RibotEM (2006) Second generation subtyping: a proposed PulseNet protocol for multiple-locus variable-number tandem repeat analysis of Shiga toxin-producing Escherichia coli O157 (STEC O157). Foodborne pathogens and disease3: 118131 1660298716602987SunQ, LanR, WangJ, XiaS, WangY, et al (2013) Identification and characterization of a novel Shigella flexneri serotype Yv in China. PloS one8: e70238 doi: 10.1371/journal.pone.00702382393617223936172SunQ, KnirelYA, LanR, WangJ, SenchenkovaSN, et al (2012) A novel plasmid-encoded serotype conversion mechanism through addition of phosphoethanolamine to the O-antigen of Shigella flexneri. PloS one7: e46095 doi: 10.1371/journal.pone.00460952304994723049947KnirelYA, LanR, SenchenkovaSN, WangJ, ShashkovAS, et al (2013) O-antigen structure of Shigella flexneri serotype Yv and effect of the lpt-O gene variation on phosphoethanolamine modification of S. flexneri O-antigens. Glycobiology23: 475485 doi: 10.1093/glycob/cws2222328300023283000Zhu LiTQ, ZouW, ZhangM (2013) Analysis on serotype and antimicrobial resistance of Shigella during 2010–2011 in Suzhou. Jiangsu Journal of Preventive Medicine24: 3.Ru WHL, ZhaoJ (2011) Serotype distribution and drug resistance analysis of shigella in Henan province between 2001 and 2010. Chinese Journal of Practical Medicine38: 4.Zang Z LJ, Xu Y, Wu D (2009) Analysis of Community-acquired Shigella Infection and Drug Resistance in Xian Area 2004–2008Journal of Modern Laboratory Medicine24: 3.Liu JMY, WangL, CuiJ, WangY, JiY (2011) Analysis on Change of Serum Types and Drug Resistance of Shigella in Xicheng District of Beijing from 2008 to 2010. Occupation and Health27: 3.Chen WWJ, YuW, YuY (2006) Study on the distribution and the level of resistance to drugs of Shigella in Jiading District in 2004. Modern preventive medicine33: 2.CaoY, WeiD, KamaraIL, ChenW (2012) Multi-Locus Sequence Typing (MLST) and Repetitive Extragenic Palindromic Polymerase Chain Reaction (REP-PCR), characterization of shigella spp. over two decades in Tianjin China. International journal of molecular epidemiology and genetics3: 321332 2320518423205184ChoiSY, JeonYS, LeeJH, ChoiB, MoonSH, et al (2007) Multilocus sequence typing analysis of Shigella flexneri isolates collected in Asian countries. Journal of medical microbiology56: 14601466 1796534517965345RibotEM, FairMA, GautomR, CameronDN, HunterSB, et al (2006) Standardization of pulsed-field gel electrophoresis protocols for the subtyping of Escherichia coli O157:H7, Salmonella, and Shigella for PulseNet. Foodborne pathogens and disease3: 5967 1660298016602980GorgeO, LopezS, HilaireV, LisantiO, RamisseV, et al (2008) Selection and validation of a multilocus variable-number tandem-repeat analysis panel for typing Shigella spp. Journal of clinical microbiology46: 10261036 doi: 10.1128/JCM.02027-071821621418216214IslamMA, HuqM, NabiA, TalukdarPK, AhmedD, et al (2013) Occurrence and characterization of multidrug-resistant New Delhi metallo-beta-lactamase-1-producing bacteria isolated between 2003 and 2010 in Bangladesh. Journal of medical microbiology62: 6268 doi: 10.1099/jmm.0.048066-02301919123019191AhmedAM, ShimamotoT (2014) Isolation and molecular characterization of Salmonella enterica, Escherichia coli O157:H7 and Shigella spp. from meat and dairy products in Egypt. International journal of food microbiology168–169: 5762 doi: 10.1016/j.ijfoodmicro.2013.10.01424239976