archcid-65855
archcid-65855

article information

Archives of Clinical Infectious Diseases: Vol.14, issue 3; e65855
published online: June 15, 2019
article type: Research Article
received: January 3, 2018
revised: March 29, 2019
accepted: May 5, 2019
DOI: https://doi.org/10.5812/archcid.65855
Crossmark

Multi Locus VNTR (MLVA) Typing and Detection of the OI-122 Pathogenicity Island in Typical and Atypical Enteropathogenic Escherichia coli Isolated from Children with Acute Diarrhea

authors:

avatar Mojtaba Memariani ORCID 1 , avatar Shahin Najar-Peerayeh 1 , * , avatar Taghi Zahraei Salehi ORCID 2

Department of Bacteriology, Faculty of Medical Sciences, Tarbiat Modares University, Tehran, Iran
Department of Microbiology, Faculty of Veterinary Medicine, University of Tehran, Tehran, Iran

how to cite: Memariani M, Najar-Peerayeh S, Zahraei Salehi T. Multi Locus VNTR (MLVA) Typing and Detection of the OI-122 Pathogenicity Island in Typical and Atypical Enteropathogenic Escherichia coli Isolated from Children with Acute Diarrhea. Arch Clin Infect Dis. 2019;14(3):e65855. https://doi.org/10.5812/archcid.65855.

Abstract

Background:

Pathogenicity O island 122 (OI-122) associated with the severity of Enteropathogenic Escherichia coli (EPEC) disease.

Objectives:

The aim of this study was the characterization of OI-122 in Iranian EPEC isolates.

Methods:

The distribution patterns of OI-122 were investigated in 42 EPEC strains by detection of OI-122 genes and serogrouping, multi locus VNTR (MLVA) typing, and phylogenetic grouping.

Results:

The complete OI-122 was identified in 18.75% of tEPEC and in 3.8% of aEPEC strains, an incomplete OI-122 with various combinations of genes was found in 75% of tEPEC and 30.76% of strains aEPEC, and OI-122 genes were absent in 6.25% of tEPEC and 65.3% of aEPEC strains. nleB (52.4%) was the most frequently detected, and the prevalence of pagC, efa1/lifA, and sen genes were 19%, 16.7%, and 14.3%, respectively. The most common phylogenetic group among tEPEC and aEPEC strains were B1 (68.75% and 80.76%, respectively). The common serogroups among EPEC isolates were O128, O111, O55, O127, and O44, and 45.2% of the isolates were untypeable. Forty-one MLVA types among 42 EPEC strains were classified into 7 clonal complexes.

Conclusions:

The diverse distribution of OI-122 genes among MLVA clonal groups provides evidence for dynamic evolution regarding the OI-122 pathogenicity island in tEPEC and aEPEC strains. The results also indicate that the acquisition of OI-122 gene contents occurs in a modular manner.

1. Background

Escherichia coli strains are among the major causes of diarrhea particularly in infants in developing countries (1). The virulence factors associated with E. coli mediated diarrhea have been used for classification of strains as distinct diarrhoeagenic pathovar. Enteropathogenic Escherichia coli (EPEC) is one of the diarrhoeagenic pathovars that belongs to an attaching and effacing (A/E) lesion forming family of pathogens; this family includes enterohaemorrhagic E. coli (EHEC), rabbit diarrhoeagenic E. coli (RDEC), Escherichia albertii, and the murine pathogen Citrobacter rodentium (2). The A/E lesions are comprised of bacteria intimately bound to the host cell surface leading to damage of microvilli and pedestal formation. The pathogenicity island locus of enterocyte effacement (LEE) required for the production of these lesions. The core region of LEE is highly conserved, which encodes the type III secretion system, whereas a higher diversity observes among genes encoding effector proteins. EPEC strains can be classified into typical EPEC (tEPEC) and atypical EPEC (aEPEC) subgroups based on the presence or absence of the type IV-like bundle-forming pili (BFP) that encoded in a ca. 90-kb E. coli adherence factor plasmid (pEAF) (3). Besides LEE pathogenicity island (PAI), various PAIs have been found on chromosomes of EPEC strains, such as pathogenicity island OI-122.

OI-122 is a 23-kb pathogenicity island with three distinct modules in O157:H7 strain EDL 933. Module 1 carries pagC gene, module 2 carries sen, nleB and nleE genes and the third module contains efa1/lifA and efa2 genes, these modules are separated from each other by various IS elements (4). However, it was reported in various forms, complete with three distinct modules and partially with one or two modules (5). It has been shown that the presence of the OI-122 genes (i.e. efa1/lifA, sen, pagC, nleB, and nleE) in aEPEC is statistically linked with severe diseases and outbreaks (6). Other studies have found that OI-122 was more prevalent in tEPEC than in aEPEC and significantly associated with diarrhea in aEPEC strains (7, 8). In our previous study, we investigated the prevalence of the EPEC strain in clinical E. coli isolates from 363 stool samples of children from 2011 to 2013. Our results showed the presence of the EPEC strains in 42 of the isolates (9). The present study was designed to characterize the pathogenicity O island 122 (OI-122) among EPEC strains and to investigate the epidemiological relationship using MLVA.

2. Objectives

The aim of this study was to investigate the distribution patterns of this genomic island in Iranian tEPEC and aEPEC strains by detecting OI-122 genes and serogrouping, multi locus VNTR (MLVA) typing, and phylogenetic grouping.

3. Methods

3.1. Bacterial Strains

Forty-two EPEC strains, defined as eae positive E. coli without stx1 and stx2 genes, were collected from stool samples from children with diarrhea (≤ 10 years old), in the study conducted during the period from 2011 to 2013 in the four pediatric hospitals in Tehran, Iran (9). The 16 typical (bfpA+) and 26 atypical (bfpA-) EPEC isolates were stored at -70°C until further analysis.

3.2. Detection of OI-122 Genes

All of the EPEC isolates were screened by PCR for the presence of OI-122 genes (i.e. efa/lifA, sen, pagC, and nleB) using primers described elsewhere (7-11). Genomic DNA of E. coli strain EDL933 was used as positive control (Table 1).

Table 1.

Primers Used in This Study for Detection of Phylogenetic Groups, the Virulence Genes and MLVA Typing

GenePrimer Sequence (5’ to 3’)Size of Amplicon, bpTemperature, °C
chuAF: GAC GAA CCA ACG GTC AGG AT27955°C
R: TGC CGC CAG TAC CAA AGA CA
yjaAF: TGA AGT GTC AGG AGA CGC TG21155°C
R: ATG GAG AAT GCG TTC CTC AAC
TspE4C2F: GAG TAA TGT CGG GGC ATT CA15255°C
R: CGC GCC AAC AAA GTA TTA CG
efa/lifAF: AGA ATG GAA GAT CAC ACC30955°C
R: ATA ATG CCT TTC ATC CAC AC
nleBF: GCT TTC ACC GAT AAG GAC AAC27255°C
R: TCG CCA TCA ACA AAA ATA CC
pagCF: ATG AGT GGT TCA ACA CTG52054°C
R: CCA ACT CCA ACA GTA AAT CC
senF: GGA TGG AAC CAT ACC TGG55058°C
R: CGC AAT CAA TTG CTA ATG C
ms06F: AAA CGG GAG AGC CGG TTA TT3955°C
R: TGT TGG TAC AAC GGC TCC TG
ms07F: GTC AGT TCG CCC AGA CAC AG3955°C
R: CGG TGT CAG CAA ATC CAG AG
ms09F: GTG CCA TCG GGC AAA ATT AG17955°C
R: CCG ATA AGG GAG CAG GCT AGT
ms11F: GAA ACA GGC CCA GGC TAC AC9655°C
R: CTG GCG CTG GTT ATG GGT AT
ms21F: GCT GAT GGC GAA GGA GAA GA14155°C
R: GGG AGT ATG CGG TCA AAA GC
ms23F: GCT CCG CTG ATT GAC TCC TT37555°C
R: AAC TGG CGG CGT TTA TCA AG
ms32F: GAG ATT GCC GAA GTG TTG C10155°C
R: AAC TGG CGG CGT TTA TCA AG

3.3. Phylogenetic Grouping

Phylogroup determination for EPEC isolate was accomplished using the previously described PCR methodology (12). Briefly, the chuA and yjaA genes and DNA fragment TSPE4.C2 amplified by PCR. The phylogenetic group A (negative for all three genes), B1 (TSPE4.C2 positive), B2 (both chuA and yjaA genes positive) and D (chuA positive) was detected based on presence or absence of the chuA and yjaA genes and DNA fragment TSPE4.C2.

3.4. Serogrouping

The slide agglutination test was performed using EPEC O-specific polyvalent antisera (Mast, UK) according to manufacturer’s instructions. The polyvalent antisera consist of three separated pools, able to react with the following serogroups: poly group 2 (O26, O55, O111, O119, and O126), poly group 3 (O86, O114, O125, O127 and O128) and poly group 4 (O44, O112, O124 and O142). The agglutinated isolates further tested by monovalent O antisera.

3.5. MLVA Genotyping

MLVA genotyping were performed with the primer for seven loci (ms06, ms07, ms09, ms11, ms21, ms23, and ms32) (13). Repeat units were imported into Microsoft Excel 2007. The minimum spanning tree (MST) was constructed with a categorical coefficient based on allelic profiles of the strains (www.MLVAplus.net). Clonal complexes (CCs) are defined as a group of allelic profiles in which every profile shares at least 5 loci in common with at least one other member of the group (14). Genetic diversity of each VNTR locus was determined by calculating the Simpson’s index of diversity (SID) using the V-DICE software (Health Protection Agency, London, UK; http://www.hpa-bioinfotools.org.uk/cgi-bin/DICI/DICI.pl). To evaluate congruence between MLVA CCs and the other typing methods, the adjusted Wallace coefficient (AW) was calculated (www.comparingpartitions.info).

3.6. Statistical Analysis

Fisher exact test was used for data analysis. A P value of < 0.01 was considered statistically significant.

4. Results

4.1. Phylogenetic Group Determination

The most common phylogenetic group among tEPEC and aEPEC strains were B1 (68.75% and 80.76%, respectively). 4.6% of the aEPEC strains belonged to phylogenetic group A, 7.6% strains belonged to phylogenetic group D and one isolate belonged to phylogenetic group B2. 25% of the tEPEC strains belonged to phylogenetic group D and one isolate belonged to phylogenetic group B2 (Table 2).

Table 2.

Comparison of Phylogenetic Classification, OI-122 Genes, Serotype, MLVA Typing for Each of the 42 EPEC Strains Identified in a Study Among Children with Acute Diarrhea

Strainefa/lifAsennleBpagCPhylogroupSerogroupMLVA TypeHospital
tEPEC Strains
TMU29nnpnB1O128M1B
TMU53nnpnB1ONTM2M
TMU75nnnnB1ONTM3B
TMU112pnppB1O127M4A
TMU135nnpnB1O127M5M
TMU145nnpnDONTM6A
TMU155pnpnB1O128M7A
TMU157nppnDONTM8M
TMU164npppDONTM9M
TMU168nnppB1ONTM10A
TMU172npppB2O111M11B
TMU192ppppB1ONTM12A
TMU226ppppB1ONTM13B
TMU337pnnnDO44M14H
TMU346pnpnB1O111M15H
TMU391nnpnB1O111M16H
aEPEC Strains
TMU16nnnnAONTM17M
TMU19nnpnDONTM18M
TMU45nnnnB1ONTM19A
TMU49nnnnB1ONTM20M
TMU81nnpnB1O111M21M
TMU93nnpnB1O55M22B
TMU106nnpnB1O128M23M
TMU211nnnnB1O128M23M
TMU120nnppB1O55M24A
TMU130nnnnB1O111M25A
TMU132pnppDONTM26B
TMU180nnnnB1O128M27A
TMU188nnnnB1O128M28A
TMU194nnnnB1O111M29M
TMU202nnnnB1O44M30A
TMU243nnpnB1O128M31B
TMU252nnnnAONTM32H
TMU261nnnnB1ONTM33B
TMU264nnnnB1ONTM34B
TMU268nnnnB1ONTM35B
TMU278npnnB2ONTM36B
TMU329nnnnB1O55M37A
TMU353nnpnB1O127M38H
TMU362nnnnB1O55M39A
TMU369nnnnB1ONTM40A
TMU380nnnnB1O55M41A

Abbreviation: ONT; O-serogroup non-typeable.

4.2. Serogrouping

The most common serogroups among EPEC isolates were found to be the members of O128 (n = 7, 16.7%) and O111 (n = 6, 14.3%), followed by O55 (n = 5, 11.9%), O127 (n = 3, 7%), and O44 (n = 2, 4.8%). However, the majority of the isolates (n = 19, 45.2%) were untypeable (ONT) with the antisera used (Table 2).

4.3. MLVA Types

The MLVA identified 41 different allelic profiles (MLVA types) among 42 EPEC isolates. The Simpson's diversity of index was 0.975 that indicated the high discriminatory power for typing of EPEC isolates. An existence of identical size alleles for at least 5 VNTR markers, corresponding to the 71% genetic similarity cut-off was used to define MLVA clonal complexes (CCs) in this study. Therefore, 23 (54.8%) isolates were classified into 7 CCs, while the remaining isolates (45.2%) were assigned singletons. MLVA-CC1 was the largest lineage (clone) consisting of 7 MLVA types (Figure 1).

Minimum spanning tree (MST) based on MLVA data of EPEC isolates. The circles indicate the MLVA types and the size of the circles indicate the number of isolates. The number of loci that differ between two MLVA types is indicated on the lines connecting the MLVA types. The colours of the circles indicate the serogroups of the isolates and clonal complexes (CCs) are indicated by grey halos. The red dotted circles represent tEPEC isolates.
Minimum spanning tree (MST) based on MLVA data of EPEC isolates. The circles indicate the MLVA types and the size of the circles indicate the number of isolates. The number of loci that differ between two MLVA types is indicated on the lines connecting the MLVA types. The colours of the circles indicate the serogroups of the isolates and clonal complexes (CCs) are indicated by grey halos. The red dotted circles represent tEPEC isolates.
Figure 1.

Minimum spanning tree (MST) based on MLVA data of EPEC isolates. The circles indicate the MLVA types and the size of the circles indicate the number of isolates. The number of loci that differ between two MLVA types is indicated on the lines connecting the MLVA types. The colours of the circles indicate the serogroups of the isolates and clonal complexes (CCs) are indicated by grey halos. The red dotted circles represent tEPEC isolates.

MLVA analysis of O-serogroups showed that the O128 strains were found in CC1 and O55 strains were grouped in CC3, the remaining CCs were contained of other O serogroups and untypeable strains (Figure 1).

4.4. OI-122 Genes

The individual OI-122 genes were detected in 93.7% of tEPEC and 42.3% of aEPEC isolates (Table 2). The distribution of individual OI-122 modules is shown in Table 3.The complete OI-122 was identified in 3 (18.75%) of tEPEC and in 1 (3.84%) of aEPEC strains. Regarding the prevalence of the modules, module 2 was the most prevalent, followed by 1 and 3. OI-122 modules were not found in 18 (42.8%) isolates.

Table 3.

The Diversity of OI-122 Modules Among in EPEC Strains in This Study

OI-122 ModulesStrain
tEPEC, n = 16aEPEC, n = 26Total, n = 42
Complete OI-122314
Module 1---
Module 26713
Module 31-1
Module 1+2314
Module 1+3---
Module 2+32-2
Absent OI-12211718

nleB (n = 22, 52.4%) was the most frequently detected, and the prevalence of pagC, efa1/lifA, and sen genes were 19% (n = 8), 16.7% (n = 7), and 14.3%, respectively. There was a strong association between nleB and tEPEC strains (P < 0.0001). Furthermore, it was shown that efa/lifA (P = 0.008), sen (P = 0.023), and pagC (P = 0.025) had statistically significant associations with tEPEC strains (Table 2).

The EPEC strains were classified into three different groups by component analysis of the OI-122 genes in the study (Table 3). The first group included strains containing a complete OI-122 module (3 strains in tEPEC and 1 in aEPEC), an incomplete OI-122 module carrying strains classified in second group (12 strains in tEPEC and 8 in aEPEC) and the third group lack the OI-122 module (1 strains in tEPEC and 17 in aEPEC). The OI-122 genes were distributed among strains with different serogroups or even untypeable EPEC strains and there were not any clonal relationships between OI-122 genes and distinct O serogroups.

5. Discussion

The main goal of this study was to determine the prevalence and distribution of pathogenicity island OI-122 in tEPEC and aEPEC. Four principle results emerged from our data: (1) the OI-122 modules were found in all tEPEC (except one isolate), but the levels of distribution varied among the OI-122 modules observed within isolates, (2) more aEPEC (57.7%) lake of OI-122 genes and the OI-122 modules was prevalent in tEPEC than in aEPEC isolates, (3) there was statistical evidence for the existence of nleB, efa/lifA, sen and pagC genes in tEPEC genotype, and (4) there was a high degree of genetic heterogeneity among OI-122 carrying EPEC strains. These results were consistent with those previously reported; suggesting that the acquisition of complete OI-122 island, OI-122 module, and OI-122 genes are evidence of horizontal gene transfer and the evolutionary dynamics among EPEC strains (5, 15). While obtained of the OI-122 genes are associated with severe diseases and outbreaks (6), our results could not indicate the association between OI-122 module and severity of disease because all of the strains analyzed in this study were isolated from children with acute diarrhea.

A complete OI-122 gene (carrying efa1/lifA, pagC, sen, and nleB) was more prevalent in tEPEC (18.75%) than in aEPEC (3.84%). Similar results were reported by Vieira et al. and Salvador et al. (5, 15).

The efa1/lifA gene (coding for EHEC factor for adherence protein) was absent in our aEPEC strains with one exception, this gene encodes a lymphostatin, which has been shown to have adhesive potentials and inhibitory effects in lymphocyte functions (4-6). An efa1/lifA gene was found in 37.5% of tEPEC and in 3.8% of aEPEC strains in this study, the significant association with diarrhea was reported for the efa1/lifA gene (6). The efa1/lifA frequency in aEPEC strains in Japan, Brazil, Norway, Australia and New Zealand were 32.5%, 30.4%, and 28.8% respectively, which it is higher than the results in the current study for aEPEC strains (5-7).

Module 2 was the prevalent module among tEPEC and aEPEC strains in this study, which it is consistent with previous reports (5-7). This module encodes nleB (virulence factor) and sen (a putative enterotoxin) genes (16, 17). The nleB detected in the majority (87.5%) of the tEPEC strains examined, suggesting that it represents a stable acquisition of the positive clonal lineages.

The pagC gene (immunogen and bacterial survival factor within macrophages) from module 1 was present in strains carrying a complete OI-122 with all three modules and strains carrying the combination of this module and module 2. This module plays an important role (immunogen and bacterial survival factor within macrophages) in O157: H7 infections, (18, 19) but its role in EPEC infection is unknown. The heterogeneous distribution of OI-122 modules was observed in EPEC strains in this work and that of other researchers, (5-7) probably due to the acquisition of OI-122 gene contents in a modular manner.

Most of the strains (76.19%) in this study belonged to the B1 phylogenetic group and only two isolates were the B2 phylogenetic group. This typing method classifies strains according to the presence or absence of the three genes chuA, yjaA, and tspE4C2, and the only B2 group is positive for a yjaA gene (12). The significant negative association with diarrhea for the yjaA gene reported in the study of Afset et al. (6) and suggested that the presence of the yjaA gene may be a marker for low diarrheagenic potential aEPEC strains. The yjaA gene as the negative marker for diarrhea may be one possible reason for the association of B2 phylogenetic group with extraintestinal E. coli infections (20-22).

The common serogroups were O128, O111, O55 and O44 in EPEC in this study, and 50% of tEPEC and 42.3% of aEPEC strains were untypeable with the traditional O antisera. The presence of the high number of untypeable strains in this study indicates that O serogrouping is not the valuable method in the diagnosis of EPEC infections.

The extensive heterogeneity among the tEPEC and aEPEC strains in this study was found by the MLVA analysis; differences in the MLVA type were identified between strains belonging to the same phylogenetic group or O serogroups. However, two strains shared a same (M23) MLVA type (Table 2), these strains were isolated from a sister and brother, which are belong to the B1 phylogenetic group and O128 serogroup of aEPEC. Interestingly, while the OI-122 genes were absent in one of these strains, the other strain harbored the nleB gene. MLVA typing identified seven clonal complexes (7CC) in EPEC strains, but clonal relationships were not found with OI-122 genes. The MLST results of aEPEC isolates from China also showed very heterogeneity (23). Staples et al. (24) analyzed EPEC isolates by MLST and MLVA. Their isolates were highly heterogeneous. They concluded that the human enteric EPEC population might be a complex of commensal or pathogenic strains. MLVA typing of aEPEC strains from poultry suggested that poultry strains were closer to bovine strains and were less similar to patient strains (25).

In conclusion, the results of this study provide evidence for dynamic evolution regarding the OI-122 pathogenicity island in tEPEC and aEPEC strains. The OI-122 gene contents distribution in a modular manner in this study indicating that clinical, animal and cell line studies need for an understanding of the OI-122 genes roles in pathogenicity of EPEC strains. Since no specific results on the distribution of the OI-122 pathogenicity island in EPEC strains are available in Iran, the results of this study shed light on the roles of OI-122 genes in the pathogenicity and epidemiology of EPEC strains, particularly of the aEPEC strains.

References

  • 1.

    Hu J, Torres AG. Enteropathogenic Escherichia coli: foe or innocent bystander? Clin Microbiol Infect. 2015;21(8):729-34. [PubMed ID: 25726041]. [PubMed Central ID: PMC4497942]. https://doi.org/10.1016/j.cmi.2015.01.015.

  • 2.

    Croxen MA, Finlay BB. Molecular mechanisms of Escherichia coli pathogenicity. Nat Rev Microbiol. 2010;8(1):26-38. [PubMed ID: 19966814]. https://doi.org/10.1038/nrmicro2265.

  • 3.

    Tobe T, Hayashi T, Han CG, Schoolnik GK, Ohtsubo E, Sasakawa C. Complete DNA sequence and structural analysis of the enteropathogenic Escherichia coli adherence factor plasmid. Infect Immun. 1999;67(10):5455-62. [PubMed ID: 10496929]. [PubMed Central ID: PMC96904].

  • 4.

    Bugarel M, Martin A, Fach P, Beutin L. Virulence gene profiling of enterohemorrhagic (EHEC) and enteropathogenic (EPEC) Escherichia coli strains: a basis for molecular risk assessment of typical and atypical EPEC strains. BMC Microbiol. 2011;11:142. [PubMed ID: 21689465]. [PubMed Central ID: PMC3133550]. https://doi.org/10.1186/1471-2180-11-142.

  • 5.

    Vieira MA, Salvador FA, Silva RM, Irino K, Vaz TM, Rockstroh AC, et al. Prevalence and characteristics of the O122 pathogenicity island in typical and atypical enteropathogenic Escherichia coli strains. J Clin Microbiol. 2010;48(4):1452-5. [PubMed ID: 20181917]. [PubMed Central ID: PMC2849565]. https://doi.org/10.1128/JCM.01944-09.

  • 6.

    Afset JE, Bruant G, Brousseau R, Harel J, Anderssen E, Bevanger L, et al. Identification of virulence genes linked with diarrhea due to atypical enteropathogenic Escherichia coli by DNA microarray analysis and PCR. J Clin Microbiol. 2006;44(10):3703-11. [PubMed ID: 17021100]. [PubMed Central ID: PMC1594803]. https://doi.org/10.1128/JCM.00429-06.

  • 7.

    Narimatsu H, Ogata K, Makino Y, Ito K. Distribution of non-locus of enterocyte effacement pathogenic island-related genes in Escherichia coli carrying eae from patients with diarrhea and healthy individuals in Japan. J Clin Microbiol. 2010;48(11):4107-14. [PubMed ID: 20844211]. [PubMed Central ID: PMC3020830]. https://doi.org/10.1128/JCM.00677-10.

  • 8.

    Konczy P, Ziebell K, Mascarenhas M, Choi A, Michaud C, Kropinski AM, et al. Genomic O island 122, locus for enterocyte effacement, and the evolution of virulent verocytotoxin-producing Escherichia coli. J Bacteriol. 2008;190(17):5832-40. [PubMed ID: 18586943]. [PubMed Central ID: PMC2519529]. https://doi.org/10.1128/JB.00480-08.

  • 9.

    Memariani M, Najar Peerayeh S, Shokouhi Mostafavi SK, Zahraei Salehi T. Detection of class 1 and 2 integrons among enteropathogenic Escherichia coli isolates. Arch Pediatr Infect Dis. 2014;3(2). https://doi.org/10.5812/pedinfect.16372.

  • 10.

    Scaletsky IC, Aranda KR, Souza TB, Silva NP, Morais MB. Evidence of pathogenic subgroups among atypical enteropathogenic Escherichia coli strains. J Clin Microbiol. 2009;47(11):3756-9. [PubMed ID: 19759223]. [PubMed Central ID: PMC2772600]. https://doi.org/10.1128/JCM.01599-09.

  • 11.

    Makobe CK, Sang WK, Kikuvi G, Kariuki S. Molecular characterization of virulence factors in diarrhoeagenic Escherichia coli isolates from children in Nairobi, Kenya. J Infect Dev Ctries. 2012;6(8):598-604. [PubMed ID: 22910565]. https://doi.org/10.3855/jidc.2082.

  • 12.

    Clermont O, Bonacorsi S, Bingen E. Rapid and simple determination of the Escherichia coli phylogenetic group. Appl Environ Microbiol. 2000;66(10):4555-8. [PubMed ID: 11010916]. [PubMed Central ID: PMC92342]. https://doi.org/10.1128/aem.66.10.4555-4558.2000.

  • 13.

    Gorge O, Lopez S, Hilaire V, Lisanti O, Ramisse V, Vergnaud G. Selection and validation of a multilocus variable-number tandem-repeat analysis panel for typing Shigella spp. J Clin Microbiol. 2008;46(3):1026-36. [PubMed ID: 18216214]. [PubMed Central ID: PMC2268344]. https://doi.org/10.1128/JCM.02027-07.

  • 14.

    Weniger T, Krawczyk J, Supply P, Niemann S, Harmsen D. MIRU-VNTRplus: a web tool for polyphasic genotyping of Mycobacterium tuberculosis complex bacteria. Nucleic Acids Res. 2010;38(Web Server issue):W326-31. [PubMed ID: 20457747]. [PubMed Central ID: PMC2896200]. https://doi.org/10.1093/nar/gkq351.

  • 15.

    Salvador FA, Hernandes RT, Vieira MA, Rockstroh AC, Gomes TA. Distribution of non-LEE-encoded type 3 secretion system dependent effectors in enteropathogenic Escherichia coli. Braz J Microbiol. 2014;45(3):851-5. [PubMed ID: 25477918]. [PubMed Central ID: PMC4204969]. https://doi.org/10.1590/s1517-83822014000300014.

  • 16.

    Karmali MA, Mascarenhas M, Shen S, Ziebell K, Johnson S, Reid-Smith R, et al. Association of genomic O island 122 of Escherichia coli EDL 933 with verocytotoxin-producing Escherichia coli seropathotypes that are linked to epidemic and/or serious disease. J Clin Microbiol. 2003;41(11):4930-40. [PubMed ID: 14605120]. [PubMed Central ID: PMC262514]. https://doi.org/10.1128/jcm.41.11.4930-4940.2003.

  • 17.

    Kelly M, Hart E, Mundy R, Marches O, Wiles S, Badea L, et al. Essential role of the type III secretion system effector NleB in colonization of mice by Citrobacter rodentium. Infect Immun. 2006;74(4):2328-37. [PubMed ID: 16552063]. [PubMed Central ID: PMC1418941]. https://doi.org/10.1128/IAI.74.4.2328-2337.2006.

  • 18.

    Ju W, Shen J, Toro M, Zhao S, Meng J. Distribution of pathogenicity islands OI-122, OI-43/48, and OI-57 and a high-pathogenicity island in Shiga toxin-producing Escherichia coli. Appl Environ Microbiol. 2013;79(11):3406-12. [PubMed ID: 23524679]. [PubMed Central ID: PMC3648051]. https://doi.org/10.1128/AEM.03661-12.

  • 19.

    Wickham ME, Lupp C, Mascarenhas M, Vazquez A, Coombes BK, Brown NF, et al. Bacterial genetic determinants of non-O157 STEC outbreaks and hemolytic-uremic syndrome after infection. J Infect Dis. 2006;194(6):819-27. [PubMed ID: 16941350]. https://doi.org/10.1086/506620.

  • 20.

    Navidinia M, Najar Peerayeh S, Fallah F, Bakhshi B. Phylogenetic groups and pathogenicity island markers in Escherichia coli isolated from children. Jundishapur J Microbiol. 2013;6(10). https://doi.org/10.5812/jjm.8362.

  • 21.

    Yun KW, Kim DS, Kim W, Lim IS. Molecular typing of uropathogenic Escherichia coli isolated from Korean children with urinary tract infection. Korean J Pediatr. 2015;58(1):20-7. [PubMed ID: 25729395]. [PubMed Central ID: PMC4342777]. https://doi.org/10.3345/kjp.2015.58.1.20.

  • 22.

    Shokouhi Mostafavi SK, Najar Peerayeh S, Mohabbati Mobarez A, Parizi MK. Detection of virulence factors, phylogroups, serogroups and biofilm formation among CTX-M-1 positive Escherichia coli isolated from patients with pyelonephritis. Biomed Res. 2018;29(8). https://doi.org/10.4066/biomedicalresearch.29-18-270.

  • 23.

    Xu Y, Bai X, Jin Y, Hu B, Wang H, Sun H, et al. High prevalence of virulence genes in specific genotypes of atypical enteropathogenic Escherichia coli. Front Cell Infect Microbiol. 2017;7:109. [PubMed ID: 28421169]. [PubMed Central ID: PMC5378719]. https://doi.org/10.3389/fcimb.2017.00109.

  • 24.

    Staples M, Doyle CJ, Graham RM, Jennison AV. Molecular epidemiological typing of enteropathogenic Escherichia coli strains from Australian patients. Diagn Microbiol Infect Dis. 2013;75(3):320-4. [PubMed ID: 23357294]. https://doi.org/10.1016/j.diagmicrobio.2012.11.010.

  • 25.

    Parvej MS, Nakamura H, Alam MA, Wang L, Zhang S, Emura K, et al. Host range-associated clustering based on multilocus variable-number tandem-repeat analysis, phylotypes, and virulence genes of atypical enteropathogenic Escherichia coli strains. Appl Environ Microbiol. 2019;85(6). [PubMed ID: 30658974]. [PubMed Central ID: PMC6414391]. https://doi.org/10.1128/AEM.02796-18.