Abstract
Objectives:
The aim of this study was to investigate the occurrence and resistance gene content of class 1 and 2 integrons among Shigella spp. and to study the genetic diversity of isolates using the pulsed-field gel electrophoresis (PFGE) method.
Methods and Results:
A total of 32 Shigella spp. were identified from 700 stool samples of patients with diarrhea from two provinces in Iran. S. sonnei (70.8%) and S. flexneri (62.5%) were the most frequent serogroups in Tehran and Razavi Khorasan provinces, respectively. Class 2 integrons were more frequent among Shigella spp. in comparison with class 1 integrons. Three different resistance gene arrays were identified among class 1 integrons. Dihydrofolate reductase (dfrA) gene cassette was detected in 78.9% of total integrons (class 1 and 2). PFGE analysis revealed clonal dissemination (62.5%) of a single clone with identical class 2 resistance gene content in Tehran province. Comparison of our Shigella pulsotypes with those published from other countries showed similar pulsotypes in India and Korea, with identical resistance profiles, which suggests dissemination of this (these) clone(s) in Asian countries.
Conclusions:
Class 2 integrons were found to be predominant among our Shigella spp. This reflects the need to monitor the acquisition and dissemination of different resistant gene cassettes among integrons. Comparison of PFGE pattern through standard procedures promoted the molecular epidemiological surveys and identification of clonal isolates in Iran and other Asian countries.
Introduction
B
Four serogroups of Shigella including serogroup A (S. dysenteriae), serogroup B (S. flexneri), serogroup C (S. boydii), and serogroup D (S. sonnei) are responsible for shigellosis, all of which except S. sonnei are further divided into subserogroups. Prevalence of different serogroups is variable worldwide.
The main virulence factors of Shigella are carried by the inv plasmid (a 200-kb virulence plasmid that contains the essential genetic information for invasion). The ipaH gene that is involved in Shigella invasion is encoded on both the inv plasmid and bacterial chromosome. This would make it possible to detect Shigella spp. through polymerase chain reaction (PCR) amplification of the ipaH genetic marker (Antoine et al., 2010).
Spread of antimicrobial resistance has increased since the 1960s and is a major problem in global public health. Shigella was the first organism that has been shown to harbor transferable antimicrobial resistance factors (Tsutomu, 1963). Mobile genetic elements such as plasmids and transposons are reported in Shigella spp. (Terajima et al., 2004). Integrons that are also reported in Shigella spp. are not mobile but can exchange exogenous DNA, known as gene cassettes, by site-specific recombination mechanisms. The most gene cassettes that have been identified inside integrons are those conferring resistance to antimicrobials (Ahmed et al., 2006).
Typing of Shigella has been achieved by both phenotypic and genotypic methods. The methodology in the first step is based on conventional phenotypic methods such as serotyping and antimicrobial susceptibility, but their low discriminatory power has promoted the development of various genotyping methods (Navia et al., 2005). Several genotyping methods have been developed for Shigella spp., among which pulsed-field gel electrophoresis (PFGE) has been recognized to be a powerful tool for discriminating the clonal relationship among Shigella strains in different laboratories (Ribot et al., 2006).
The aim of this study was to determine the serogroup distribution and resistance gene content of two classes of integrons among Shigella isolates and to investigate the genetic diversity of integron harboring isolates using the PFGE genotyping method.
Materials and Methods
Bacterial strains
Thirty-two strains of Shigella spp. from Tehran (24) (Iran Central) and Razavi Khorasan (8) (Northeast of Iran) were obtained from 700 patients with acute diarrhea from two provinces in Iran (Tehran=500, Razavi Khorasan=200) during an 18-month period from April 2009 to November 2010 with no evidence of outbreak. The prevalence of shigellosis was 4.8% and 4% during the study period in two provinces, respectively. Biochemical screening tests were done to characterize the Shigella spp. isolates from selective agar plates including MacConkey agar, Salmonella-Shigella agar and XLD agar (Merck, Hamburg, Germany). The identity of isolates was confirmed using a genus-specific primer set that amplifies within the ipaH regions. Serogrouping of isolates was performed using monospecific antisera (Baharafshan Institute of Research & Development, Tehran, Iran). S. flexneri ATCC 9290 and S. sonnei ATCC 1202 were used as controls in each assay. S. boydii and S. dysenteriae control strains were kindly provided by the World Health Organization.
Antimicrobial susceptibility
Antimicrobial susceptibility testing was performed using the disc diffusion method according to CLSI guidelines (Wikler, 2006) for ampicillin (10 μg), co-trimoxasol (25 μg), tetracycline (30 μg), chloramphenicol (30 μg), ciprofloxacin (5 μg), trimethoprim (5 μg), and streptomycin (10 μg) (Wikler, 2006). Escherichia coli ATCC 25922 was used as a control in antimicrobial susceptibility testing.
Integron characterization and sequencing of resistance-encoding gene cassettes
The total DNA was extracted by boiling and used as a template in PCR assays for ipaH gene, class1 and 2 integrons. The primer sequences used in this study are shown in Table 1.
PCR was performed in a reaction mixture with total volume of 25 μL, containing 2.5 μL 10x Taq polymerase buffer, 0.3 μL dNTPs (10 mmol L−1), 1 U Taq DNA polymerase, 0.6 μL MgCl2 (50 mmol L−1) and 0.3 mol L−1 from each primer. PCR was done as follows: initial denaturation step at 94°C for 5 min followed by 30 cycles consisting of denaturation (94°C for 1 min), annealing (56°C for 1 min), and extension (72°C for 1 min), followed by a final extension step at 72°C for 3 min. PCR products were purified using QIAquick Gel Extraction Kit (Qiagen), and direct sequencing of internal variable regions (gene cassettes) of class1 and 2 integrons was done using ABI 3730X capillary sequencer (Genfanavaran, Macrogen, Seoul, Korea).
GenBank accession numbers
One representative of each resistance gene cassette array inside class 1 integrons was deposited in the GenBank database. The assigned GenBank accession numbers are JX035724, JX035725, JX035726, and JX035727.
Pulsed-field gel electrophoresis (PFGE)
PulseNet standardized protocol was used for subtyping of Shigella spp. isolates (Ribot et al., 2006). In brief, bacteria from culture plates were adjusted to absorbance values of 0.8–1.0 at a wavelength of 610 nm in a cell suspension buffer (100 mmol L−1 Tris: 100 mmol l−1 EDTA, pH 8.0) after which agarose plugs were prepared using SeaKem Gold agarose (Lonza, Rockland, ME) and proteinase K. Cells in the agarose plugs were lysed by treatment with a lysis solution 50 mmol L−1 Tris, 50 mmol L−1 EDTA (pH 8.0), 1% sarcosine, and 0.5 mg of proteinase K for 1 h at 54°C. Six stages of washing steps were performed, twice with sterile ultrapure water and four times with Tris-EDTA (TE) buffer (10 mmol L−1 Tris, 1 mmol L−1 EDTA, pH 8.0). Forty units of XbaI restriction enzyme (Roche Diagnostic GmbH, Mannheim, Germany) was applied for plugs in a freshly prepared buffer and incubated for 4 h. XbaI digested Salmonella enterica serotype Braenderup H9812 plugs were used as DNA molecular weight size marker. The electrophoresis was performed with CHEF Mapper XA System (Bio-Rad) and consisted of 200 V at 14°C for 18 h with the increasing pulsed time from 2.16 s to 54.17 s. PFGE patterns were analyzed by Gel Compare II version 4.0 software (Applied Maths, Sint-Matenslatem, Belgium) and the patterns were compared by using the Dice coefficient and UPGMA (unweighted pair group method with arithmetic averages) clustering. A dendrogram was constructed using an optimization value of 0.50% and a position tolerance of 1.0%. Isolates with 80% similarity were considered as a single cluster. Assigning the types and interpretation of PFGE generated patterns was performed according to the guidelines set by Tenover and colleagues (Tenover et al., 1995).
Results
Serogrouping
Serogrouping testing was done for all 32 strains, and the results showed the predominance of S. sonnei (70.8%) and S. flexneri (20%) in Tehran, while in Razavi Khorasan S. flexneri was 62.5%. S. boydii constituted 8.3% and 12.5% of total Shigella strains in two provinces, respectively, while no S. dysenteriae was detected in any of the two provinces (Table 2).
WHO, World Health Organization; ATCC, American Type Culture Collection.
Antimicrobial susceptibility analysis
Forty-five percent of isolates were multidrug resistant (resistant to three classes of antimicrobial agents or more). The foremost resistance profile was SXT/TE/TMP (34.5%) (Table 3). The most prevalent resistance was seen to trimethoprim and co-trimoxazol (93.7%) and tetracycline (87.5%), while only one isolate (3.1%) was resistant to ciprofloxacin. Resistance to streptomycin, ampicillin, and chloramphenicol was 43.7%, 37.5%, and 21.8%, respectively.
Antimicrobials were classified as aminoglycosides including Streptomycin (STR) and Ampicillin (AMP), tetracyclines including Tetracycline (TE), amphenicols including Chloramphenicol (CHL), fluoroquinolones including Ciprofloxacin (CIP), trimethoprim alone (TMP), or in combination with sulfametoxazol (SXT).
MDR, multidrug-resistant to ≥3 classes of antimicrobial agents.
PCR assay and sequencing analysis of integrons
All of the isolates produced positive results in ipaH amplification. Among the total 32 strains, 13 (40.6%) carried class 1 integron (int +) using primers int1-F, int1-R (5'-conserved region) with DNA bands of 900 bp. PCR amplification of internal variable region using in-F, in-R primers produced three different products of approximately 750, 1500, and 1600 sizes in two (15.4%), one (7.7%), and three (23.1%) of int + isolates, respectively. No product was obtained for seven (54%) of the int + isolates, indicative of an empty integron with no resistance gene cassette. Sequence analysis of variable region indicated the presence of dihydrofolate reductase type I (dfrA7), aminoglycoside adenyl transferase-chloramphenicol acetyltransferase (aadB-catB3), and dihydrofolate reductase-aminoglycoside adenyltransferase (dfrA17-aadA5) resistance gene cassettes among the isolates, which correspond to 750, 1500, and 1600 bp PCR products, respectively.
Among the total 32 isolates, 25 (78.12%) harbored class 2 integron using primers hep-F, hep-R, which produced two PCR products of 1500 and 2300 bp in 9 (36%) and 16 (64%) of isolates, respectively. Sequence analysis showed that these two amplification bands corresponded to dihydrofolate reductase-streptothricin acetyltransferase (dfrA1-sat1) and dihydrofolate reductase-streptothricin acetyltransferase-aminoglycoside adenyltransferase (dfrA1-sat1-aadA1) genes, respectively.
PFGE analysis
The PFGE patterns were compared using Dice coefficient and UPGMA analysis and 11 major clusters were obtained that differed by more than three bands and further divided to 21 subclusters according to criteria described earlier (Tenover et al., 1995). Cluster analysis of the isolates revealed a single major cluster (subtypes B1–B4) constituting 46.8% of the total isolates. Although all of these strains belonged to serogroup D and were isolated from Tehran province, they showed different antimicrobial resistance profiles and resistance gene cassettes of class 1 and 2 integrons. One isolate from Razavi Khorasan province fell to this cluster with strains from Tehran (Fig. 1).

Unweighted pair group method with arithmetic averages dendrogram showing banding patterns of XbaI digested 36 Shigella isolates. The isolates are presented in comparison with their phenotypic resistance profile and integron resistance gene content. Main clusters are labeled to the right of the dendrogram.
Discussion
Previous studies have shown the prevalence of shigellosis to be 11.4%, 8.8%, and 14% during 2001–2006 in different provinces of Shiraz (central), Golestan, and Mazandaran (north), respectively (Farshad et al., 2006; Ghaemi et al., 2007; Savadkoohi et al., 2007; Pourakbari et al., 2010). In this study, the prevalence of shigellosis was shown to be 4% and 4.8% in Tehran and Razavi Khorasan, respectively. Furthermore, S. sonnei and S. flexneri were identified to be predominant species in these two provinces. These studies were conducted from 2006 to 2008, and they identified S. sonnei and S. flexneri as the dominant isolates at different time intervals (Farshad et al., 2006; Rahbar et al., 2007). The reasons for the fluctuation between S. sonnei or S. flexneri could be due to (1) travelers from neighboring countries and/or (2) overall antimicrobial use inside the country. Similar serogroups have been reported from neighboring countries such as Saudi Arabia (44%), Pakistan (58%), and Jordan (65%) (Kagalwalla et al., 1992; Rawashdeh et al., 1994; Zafar et al., 2009). On the other hand, S. sonnei was continuously reported to be the predominant strain in the United States, Canada, and the developed countries (Farshad et al., 2006).
Multidrug-resistant phenotype was seen in 45% of total isolates, with the SXTr/TEr/TMPr being the foremost resistance pattern. Many other reports also indicated the dominance of resistance to streptomycin, sulfamethoxazole, and trimethoprim among Shigella spp. (Seol et al., 2006; Vrints et al., 2009).
From 32 strains in this study, 13 (S. sonnei and S. flexneri) were positive for class 1 integron (int1+), while the dfrA gene cassette was detected in 78.9% of total integrons. The results suggest the stability of this gene cassette in two classes of integrins, which may reflect the selective pressure, exerted over a long period, by the use of trimethoprim in clinics. The gene cassette arrays of aadA-dfrA detected in 31.6% of integrons may reflect co-transfer of resistance genes due to the genetic linkage of dfrA and aadA cassettes. Our previous published data indicated the predominance of this gene cassette arrays among both EPEC and non-EPEC isolates, which provide evidence of interspecies transfer of integrons through horizontal gene transfer (Najibi et al., 2012). Our data are also consistent with the previous reports worldwide on the predominance of dfrA and aadA gene cassettes among Enterobacteriaceae (Chang et al., 2007; Ajiboye et al., 2009; Scaletsky et al., 2010).
Among 32 isolates of Shigella spp., 25 isolates (S. sonnei, S. flexneri, S. boydii) carried class 2 integrons. Two different class 2 integrons (1.5 and 2.3 kb in sizes) were identified among our isolates. These two arrays (dfrA1-sat1 and dfrA1-sat1-aadA1) are identical to those reported elsewhere, although the occurrence of a 1.5-kb cassette array among our collection was less than the 2.3-kb gene cassette. This is in contrast to the previous reports from Iran, Japan, Belgium, and Korea (Ahmed et al., 2006; Seol et al., 2006; Ranjbar et al., 2007; Vrints et al., 2009).
All streptomycin-resistant isolates contained one member of aadA/B gene cassette in class 1 or 2 integron, which emphasizes the significance of integrons in spread of aadA/B resistance gene family. The occurrence of class 1 integron was less than class 2 and showed three different gene cassette arrays. Since past years, several other reports have indicated the elevated occurrence of class 2 integrons in Shigella in comparison to class 1 (McIver et al., 2002; Oh et al., 2003; Mammina et al., 2005).
Some of our Shigella isolates (6.2%) harbored both class 1 and 2 integrons simultaneously, which reveals the capability of these species in acquiring additional resistance gene cassettes.
Several typing methods have been developed to assess the clonal relatedness of bacterial species; however, PFGE has been selected as a method of choice for Shigella spp. typing (Woodward et al., 2005). The results of PFGE analysis showed a single major cluster (cluster B) with 15 members constituting 46.8% of total isolates, all of which (except three isolates) uniformly showed identical class 2 integron resistance gene content. All of the isolates in this cluster were isolated from Tehran province (except 1 isolate from Razavi Khorasan), which suggests the clonal dissemination of a single clone throughout this province. Isolates from Razavi Khorasan showed more dispersed pulsotypes, which may be due to pilgrimage ceremonies and passenger trafficking to Afghanistan in this province. One interesting finding, however, was the absence of class 1 integron in the isolate.
Comparing the S. sonnei PFGE patterns of this study with the patterns in Korea and India has shown high pulsotype similarity and identical antimicrobial patterns, which suggests (1) dissemination or importation of the same clones in Asian countries, probably through human travel and/or migration; and (2) appearance of similar genotypes in different countries through genomic rearrangement events (Oh et al., 2003; Pazhani et al., 2008). Our previously published data have shown identical clonal isolates of other enteric pathogens (i.e., Vibrio cholerae in Iran and other Asian countries including India, Nepal, and Pakistan) (Bakhshi and Pourshafie, 2009).
Shigella standard strains fell into separate pulsotypes with no integrons within their genome.
Conclusion
In conclusion, our study focused on the analysis of class 1 and class 2 integron gene cassettes, while class 2 integrons were found to be predominant among our Shigella spp. No clear relationship was defined between specific pulsotypes and integron content of isolates. The isolates with identical pulsotypes showed different resistance profile and integron content. Global comparison of PFGE patterns through PulseNet standard procedures promoted the epidemiological investigations that have led to identification of clonal isolates in Iran and other Asian countries.
Footnotes
Acknowledgments
We thank the Research Council of Tarbiat Modares University for supporting the project. Many thanks to the World Health Organization for providing Shigella standard strains.
Disclosure Statement
No competing financial interests exist.
