ORIGINAL RESEARCH article
High Prevalence of Virulence Genes in Specific Genotypes of Atypical Enteropathogenic Escherichia coli
- 1State Key Laboratory of Infectious Disease Prevention and Control, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
- 2Longgang Center for Disease Control and Prevention, Shenzhen, China
- 3Shandong Center for Disease Control and Prevention, Jinan, China
- 4Zigong Center for Disease Control and Prevention, Zigong, China
Atypical enteropathogenic Escherichia coli (aEPEC) strains are emerging enteropathogens that have been detected worldwide. A collection of 228 aEPEC strains (121 from diarrheal patients, 27 from healthy carriers, 47 from animals and 33 from raw meats) were investigated for serotypes, virulence gene profiles and phylogenetic relationships. Sixty-six O serogroups were identified. Serogroup O51 was the most prevalent, followed by O119, O26 and O76. For the 20 virulence genes detected, statistically significant differences were observed in the overall prevalence of efa1 (lifA), nleB, nleE, set/ent, paa, and ehxA genes among strains from diarrheal patients, healthy carriers, animals and raw meats, respectively. Strains from diarrheal patients had significantly higher levels of efa1 (lifA) (29.8 vs. 0%, P = 0.0002), nleB (41.3 vs. 7.4%, P = 0.0004), nleE (43.8 vs. 7.4%, P = 0.0002) and set/ent (41.3 vs. 7.4%, P = 0.0004) genes than strains obtained from healthy carriers. The paa gene was identified more often in isolates from raw meats (63.6 vs. 14.8%, P < 0.0001), animals (42.6 vs. 14.8%, P < 0.0122), and diarrheal patients (36.4 vs. 14.8%, P < 0.0225) than in strains obtained from healthy carriers. The ehxA gene was detected more frequently in strains from raw meats than in strains from diarrheal patients (27.3 vs. 2.5%, P = 0.0000) and healthy carriers (27.3 vs. 7.4%, P = 0.0474). The phylogenetic marker, yjaA, was more frequently observed in strains among healthy carriers than in diarrheal patient strains. Among the 228 aEPEC strains, 79 sequence types (STs) were identified. The prominent STs, which comprised strains carrying the four OI-122 genes and lpfA, were ST40, ST328, and ST29. Overall, the results indicate that aEPEC strains isolated in China are highly heterogeneous. aEPEC strains that are potentially more pathogenic appear to be related to specific STs or clonal complexes and serotypes. The high prevalence of diarrhea-associated genes in animal or raw meat strains suggests a zoonotic transmission pathway for potentially human pathogenic aEPEC.
Globally, one in ten childhood deaths are due to diarrheal disease among children under 5 years old, leading to about 800,000 fatalities worldwide annually, most of which occur in areas of sub-Saharan Africa and south Asia (Kotloff et al., 2013). Enteropathogenic Escherichia coli (EPEC), the first pathogenic type of E. coli to be associated with human disease, is one of the most prevalent pathogen infecting children worldwide (Ochoa and Contreras, 2011). EPEC strains are able to form attaching and effacing (A/E) lesions through intimate adherence and effacement of the intestinal microvilli of the host. A/E lesion production is dependent on the pathogenicity island LEE (locus of enterocyte effacement), which is integral to EPEC pathogenicity. Based on the presence or absence of the EPEC adherence factor plasmid (pEAF), EPEC may be further classified as typical EPEC (tEPEC) or atypical EPEC (aEPEC), respectively (Croxen et al., 2013). In the last 10 years, epidemiological studies of tEPEC and aEPEC have shown a much higher prevalence of aEPEC in both developed and developing countries, and aEPEC has emerged as an important pathogen (Hu and Torres, 2015). We, and other investigators, have clearly shown its high prevalence in China.
The majority of tEPEC strains produce virulence traits on LEE and pEAF, while aEPEC strains mostly have both additional and heterogeneous virulence properties (Hu and Torres, 2015). Bando et al. found that aEPEC strains showed genetic similarity with other diarrheagenic E. coli (DEC) pathotypes and contained virulence factors derived from other DEC pathotypes. So they inferred that aEPEC strains might have a particular genetic background that allows the acquisition and expression of virulent genes derived from other pathotypes (Bando et al., 2009). As a consequence, aEPEC strains are phylogenetically heterogeneous and carry virulence factors of other pathogenic E. coli more often than tEPEC strains (Hernandes et al., 2009). However, the virulence determinants linked with aEPEC infection have yet to be determined.
Besides the LEE region, another pathogenicity island (PAI), OI-122 has been described in aEPEC. This PAI harbors genes efa1 (lifA) that encodes lymphocyte inhibitory factor (LifA) (Abu-Median et al., 2006), set/ent that encodes a homologous Shigella flexneri enterotoxin (Nataro et al., 1995), and nleB and nleE genes that encode proteins, which inhibit pro-inflammatory signaling (Newton et al., 2010). The simultaneous presence of all the OI-122 genes i.e., efa1 (lifA), set/ent, nleB, and nleE, is statistically associated with diarrhea compared to healthy controls (Afset et al., 2006; Mercado et al., 2016).
In addition to the LEE-encoded factors, other virulence determinants may contribute to the pathogenesis of aEPEC. These include astA, which encodes the enteroaggregative heat-stable toxin 1 (EAST1), and has been shown to be significantly associated with diarrhea in a case-control study in Brazil (Dulguer et al., 2003) and a waterborne outbreak of diarrhea in Japan (Yatsuyanagi et al., 2003). Two toxins can cause apoptosis: cytolethal distending toxin (CDT) which causes characteristic and irreversible cell cycle arrest (Jinadasa et al., 2011), and subtilase cytotoxin (SubAB) which triggers endoplasmic reticulum stress signaling pathways leading to apoptosis (Paton and Paton, 2010). Autotransporters are often associated with virulence functions such as adherence, aggregation, invasion, biofilm formation, and toxicity (Abreu et al., 2013). Pathogenesis is mediated by the plasmid encoded toxin (Pet) (Ruiz et al., 2014), the Pic serine protease (Abreu et al., 2016), and an extracellular serine protease (EspP) (Brunder et al., 1999). Antiaggregation protein (dispersin), encoded by the aap gene, acts by dispersing the pathogen through the mucus layer produced by the intestinal epithelial cells, and is translocated by the antiaggregation transporter protein encoded by aat (Sheikh et al., 2002). Other putative virulence factors such as an enterohemolysin (EhxA) and a catalase/peroxidase (KatP), encoded by the genes located in a megaplasmid of Shiga toxin-producing E. coli (STEC) have also been reported (Kobayashi et al., 2013). The phylogenetic marker gene yjaA, was reported to be negatively associated with diarrhea (Wang et al., 2013).
The first step in gastrointestinal infection is the initial attachment to the surface of the host intestinal epithelium (Kalita et al., 2014). Several adherence-related factors have been described as being present in EPEC or STEC isolates, such as, the enterohemorrhagic E. coli (EHEC) autotransporter C, encoded by ehaC (Abreu et al., 2013); the Iha virulence factor, encoded by iha (Tarr et al., 2000; Bardiau et al., 2009); ToxB, encoded by toxB on the pO157 plasmid (Tatsuno et al., 2001); an autoagglutinating adhesin, Saa (Paton et al., 2001); porcine attaching and effacing-associated adhesin (Paa), encoded by paa (Batisson et al., 2003; Maluta et al., 2014); long polar fimbria (LPF) (Bardiau et al., 2010).
Since tEPEC strains are rarely isolated from animals, humans are generally considered to be the main reservoir of tEPEC. Whereas aEPEC are present in both healthy and diseased animals and humans (Hernandes et al., 2009; Hu and Torres, 2015). In the current investigation, we characterized 228 aEPEC strains recovered from multiple sources (diarrheal patients, healthy carriers, animals and raw meats), and attempted to investigate their molecular characteristics, phylogenetic relationship, and understand the potential factors involved in human infection.
Materials and Methods
In a previous study, a total of 3401 specimens from different sources were collected from 2006 to 2015 in seven geographical regions in China. In which, stool samples of diarrheal patients were collected during clinical treatment at sentinel hospitals, fecal samples of healthy humans (without diarrhea in recent 2 weeks) were collected during routine physical examination, and stool samples of animals (cattle, pig, live chicken, bird, marmota, ochtona) and raw meat samples (chicken meat, beef, pork, mutton) were collected in routine surveys (Xu et al., 2016). Briefly, 1.5 ml of each enrichment [cultured overnight in modified Tryptone Soya Broth (mTSB) (Oxoid, UK)] was centrifuged. Then, 150 μl of the lysis buffer (100 mM NaCl, 10 mM Tris–HCl [pH 8.3], 1 mM EDTA [pH 9.0], 1% Triton X-100) was added, resuspended, boiled and centrifuged. The released DNA were subjected to PCR assay for eae gene (Hu and Torres, 2015). The eae-positive enrichment culture were inoculated onto CHROMagar™ ECC plate (CHROMagar, France) and further processed for isolation of aEPEC as previously described (Xu et al., 2016). The eae-positive colonies were subcultured on Luria-Bertani (Oxoid, UK) plates and incubated for further confirmation by standard biochemical tests. Isolates with eae+, stx1/stx2−, and bfpA− were confirmed as aEPEC (Hu and Torres, 2015). Only one isolate from each sample was kept for further analysis.
Totally, 228 aEPEC isolates were identified during the epidemiological studies. 121, 27, and 47 strains were isolated from the fecal samples of diarrheal patients, healthy carriers and animals, respectively. The remaining 33 isolates were recovered from raw meats. Bacteria were stored at −80°C before being subcultured aerobically at 37°C on Luria-Bertani (LB) agar at 37°C.
Extraction of DNA from aEPEC Isolates
To prepare template DNA for PCR, the distilled water-boiling method was employed (Ooka et al., 2009). A single colony from the aEPEC isolates was suspended in 100 μl of deionized water and boiled for 10 min. After being centrifuged at 10,000 × g for 5 min, the supernatant containing DNA was collected, and 1 μl of the supernatant was used as template DNA.
Serotyping of aEPEC Isolates
The O serogroups were screened by the PCR-based system arranged by Iguchi et al. (2015). The PCR results were confirmed using complete E. coli O antisera (Statens Serum Institute, Denmark). The isolates were defined as O-untypable if they did not react with any O antisera. H typing was performed by PCR amplification and sequencing of the fliC gene using primers fliC-F (5′-ATGGCACAAGTCATTAATACCCAAC-3′) and fliC-R (5′-CTAACCCTGCAGCAGAGACA-3′) (Fields et al., 1997). These were then compared to the known variants of fliC gene on the SerotypeFinder database (https://cge.cbs.dtu.dk/services/SerotypeFinder/) (Joensen et al., 2015). The isolate was considered as H-untypable if fliC was negative by PCR.
Detection of Virulence Genes by PCR
The 228 aEPEC strains were investigated by PCR for the presence of putative virulence genes or adhesion genes (ehxA, subAB, nleB, nleE, set/ent, espP, katP, astA, pet, aat, ehaC, pic, aap and cdt, iha, efa1, lfpA, saa, toxB, and paa) and a phylogenetic marker gene yja. Three variants (lfpAO113, lfpAO157-154, lfpAO157-141) of the LPF encoding gene (lpfA) and two variants (cdt1, cdt2) of the cdt gene were included. PCR amplifications were performed in a thermal cycler (SensoQuest Labcycler, Germany). PCR cycling was performed according to the following conditions: initial denaturation at 95°C for 5 min, 30 cycles, each of denaturation at 95°C for 1 min, then annealing at the corresponding annealing temperature of 1 min and nucleotide extension at 72°C for 1 min. All primers, sizes of amplified fragments and PCR conditions used for this screening and corresponding references are shown in Table 1. Amplified products were analyzed by 1.5% agarose gel electrophoresis.
Multilocus Sequence Typing (MLST)
MLST was carried out on seven conserved housekeeping genes (adk, fumC, gyrB, icd, mdh, purA, and recA) according to the scheme of the E. coli MLST database (http://mlst.warwick.ac.uk/mlst/mlst/mlst/dbs/Ecoli/). The PCR products were purified with the QIAquick PCR purification kit (Qiagen, Germany), and double-strand sequenced using the ABI 3730 Automated DNA Analyzer (Applied Biosystems, USA). Each of the seven gene loci was assigned an allele number upon submission of the sequences to the E. coli MLST database. Allelic sequences with new variations gained new allele numbers. The allelic profile was used to generate a specific sequence type (ST) for each isolate.
For each strain, seven gene sequences were concatenated in the following order: adk, fumC, gyrB, icd, mdh, purA, and recA, to generate a 3423-bp DNA sequence for phylogenetic analysis. The concatamers were then aligned using the ClustalW program from the MEGA 6 software (http://www.megasoftware.net/). A Neighbor-Joining tree was constructed based on the maximum composite likelihood model with 1,000 bootstrap re-samplings using MEGA 6. To identify closely related genotypes, a minimum spanning tree (MST) based on different STs was constructed using Bionumerics software, version 4.6 (Applied Maths, Belgium).
The presence of each virulence gene in aEPEC isolates from different sources were compared with each other by a two-tailed chi-square test or a Fisher's exact test using Epi Info software, version 3.5.3 (Maldonado et al., 2014). P-values < 0.05 were considered statistically significant.
Serogroups and Serotypes
The 228 aEPEC isolates belonged to 66 O serogroups with the exception of 17 (7.5%) strains being O-untypable (ONT) (Table 2, Table S1). Serogroup O51 was the most prevalent and was identified in 20 (8.8%) isolates, followed by O119 and O26, which were identified in 12 (5.3%) and 10 (4.4%) isolates, respectively. The remaining serogroups consisted of between 1 and 9 aEPEC strains. In diarrheal patients, animals and raw meats, the most common serogroups detected were O51, O119, and O76, respectively. No predominant serogroup was observed in strains obtained from healthy carriers.
Twenty-five different H-types were detected among the aEPEC strains (Table 2, Table S1). The most common were the flagellar types H11 (22 strains, 9.6%), H21 (20 strains, 8.8%), and H7 (20 strains, 8.8%). The H-type was not identified in 11 strains (4.8%). The most frequently detected H-types in diarrheal patients, animals and raw meats were H19, H21, and H7, respectively. No prominent H-type was observed in strains from healthy carriers.
Virulence Gene Profiles
With the exception of saa, subAB, espP, aat, aap, and pic, all of the genes investigated were detected (Table 3, Table S1). The most frequently detected virulence gene among the 228 aEPEC isolates was ehaC (n = 183, 80.3%). The isolates from diarrheal patients, healthy carriers, animals, and raw meats exhibited a high prevalence of ehaC, with a detection rate of 82.6, 59.3, 83.0, and 84.8%, respectively. Significant difference in the overall presence of ehaC among the four different sources (P = 0.0350) was observed. ehaC was identified more often in isolates from diarrheal patients (82.6 vs. 59.3%, P = 0.0066), animals (83.0 vs. 59.3%, P = 0.0158) and raw meats (84.8 vs. 59.3%, P = 0.0159) than in isolates from healthy carriers.
Among the adhesion genes investigated, lpfA exhibited the highest prevalence in animal strains (66.0%), followed by diarrheal patient strains (54.5%), raw meat strains (48.5%), and healthy carrier strains (40.7%). However, we could not detect a significant difference in the overall presence of lpfA. The paa gene was identified with the highest prevalence in raw meat strains (63.6%), followed by animal strains (42.6%), diarrheal patient strains (36.4%), and healthy carrier strains (14.8%). There was a statistically significant difference in the overall presence of paa between the four different sources (P = 0.0013). paa was identified more often in isolates from raw meats (63.6 vs. 14.8%, P < 0.0001), animals (42.6 vs. 14.8%, P < 0.0122) and diarrheal patients (36.4 vs. 14.8%, P < 0.0225) than in strains from healthy carriers. Compared to lpfA or paa, the prevalence of iha and toxB was lower in aEPEC strains (15.8 and 8.8%, respectively).
We also investigated four genes located on OI-122 (efa1 (lifA), set/ent, nleB, and nleE), and the phylogenetic marker yjaA gene. Of the 121 strains from diarrheal patients, 36 were efa1 (lifA)-positive (29.8%), 50 were nleB-positive (41.3%), 53 were nleE-positive (43.8%), 50 were set/ent-positive (41.3%), and 44 were yjaA-positive (36.4%). Of the 27 strains from healthy individuals, none were efa1 (lifA)-positive, two strains were nleB-, nleE-, and set/ent-positive (7.4%), respectively, and 17 were yjaA-positive (63.0%). Among the 47 strains from animals, 14 were efa1 (lifA)-positive (29.8%), 22 were nleB-positive (46.8%), 24 were nleE-positive (51.1%), 21 were set/ent-positive (44.7%), and 18 were yjaA-positive (38.3%). Of the 33 strains from raw meats, 9 were efa1 (lifA)-positive (27.3%), 22 were nleB-positive (66.7%), 20 were nleE-positive (60.1%), 22 were set/ent-positive (66.7%), and 14 were yjaA-positive (42.4%). The complete OI-122 was detected in 35 (28.9%), 0, 13 (27.7%), and 6 (18.2%) strains from diarrheal patients, healthy carriers, animals and raw meats, respectively (Table S1). In contrast, incomplete OI-122 was observed in 19 (15.7%), 2 (7.4%), 12 (25.5%), and 18 (54.5%) strains from diarrheal patients, healthy carriers, animals and raw meats, respectively. aEPEC strains carrying OI-O122 genes with different combinations but lacking efa1 (lifA) were also found. There were statistically significant differences in the overall presence of efa1 (lifA), nleB, nleE, and set/ent genes among the four different sources (P = 0.0130, P = 0.0001, P = 0.0002, and P = 0.0001, respectively). Strains from diarrheal patients had significantly more efa1 (lifA) (29.8% vs. 0, P = 0.0002), nleB (41.3 vs. 7.4%, P = 0.0004), nleE (43.8 vs. 7.4%, P = 0.0002) and set/ent (41.3 vs. 7.4%, P = 0.0004) genes than strains obtained from healthy carriers. Meanwhile, efa1 (lifA), nleB, nleE, and set/ent genes were each found to be more prevalent in aEPEC strains isolated from animals or raw meats than strains from healthy carriers. yjaA was more frequently observed among healthy carrier strains than in diarrheal patient isolates (63.0 vs. 36.4%, P = 0.0105) and in animal isolates (63.0 vs. 38.3%, P = 0.0354).
Among the virulence genes linked to the EHEC pO157 plasmid, ehxA and katP were detected at high prevalence in raw meat strains (27.3 and 18.2%, respectively), while these genes were detected at low prevalence in diarrheal patient strains (2.5 and 7.4%, respectively). There was a statistically significant difference in the overall presence of ehxA among the four different sources (P = 0.0001). It was detected more often in isolates from raw meats than in strains from diarrheal patients (27.3 vs. 2.5%, P = 0.0000) and healthy carriers (27.3 vs. 7.4%, P = 0.0474). However, no significant difference was observed between raw meat strains and animal strains. In contrast, cdt was present at higher prevalence in animal strains (12.8%) and healthy carrier strains (11.1%) than in raw meat strains (6.1%) and diarrheal patient strains (0.8%). Furthermore, astA and pet were observed at higher prevalences in healthy carrier strains (18.5 and 29.6%, respectively) than in diarrheal patient strains (both were 10.7%).
Multilocus Sequence Typing
Seventy-nine different sequence types (STs) were detected among 228 aEPEC isolates (Figure 1, Table S1). In diarrheal patients, healthy carriers, animals and raw meats, 53, 17, 26, and 15 STs were detected, respectively. The most frequently detected ST was ST40 (20/228, 8.8%), with 12 strains from diarrheal patients, seven from animals and one from a healthy carrier. The second predominant ST, ST10, was identified from eight diarrheal patient isolates, four raw meat isolates, three healthy carrier isolates, and one animal isolate. ST29 and ST328 comprised 14 and 13 isolates, respectively. Another 35 STs comprised 2-10 isolates, while 40 STs had only one single isolate. Among 25 STs that included three or more isolates, four STs (ST328, ST378, ST2178, and ST2346) were detected only in diarrheal patient isolates, and four STs (ST342, ST302, ST21, and ST4059) were present only in animal isolates.
Figure 1. Phylogenetic relationships and virulence profiles of the 79 STs among the 228 aEPEC isolates. An unrooted phylogenetic tree was constructed by the neighbor-joining algorithm based on the Maximum Composite Likelihood model of nucleotide substitution. Bootstrap values greater than 90% based on 1000 replications are given at the internal nodes. STs highlighted in gray were the most prominent STs containing aEPEC isolates and harbored genes located on OI-122 (nleB, nleE, set/ent, or efa1 (lifA)) and lpfA.
We constructed an MST containing 79 STs obtained in our study (Figure 2, Table S1). They were clustered into 34 singletons and 14 non-overlapping groups or clonal complexes (CC). Most STs differed from each other by 2 or more alleles. Only two STs (i.e., ST10 and ST29) contained isolates from all sources (diarrheal patients, health carriers, animals and raw meats).
Figure 2. Minimal spanning tree of the 228 aEPEC isolates based on STs. Each circle represents a given ST with size proportional to the number of isolates. The colors for the slices of the pie represent different sources of isolates: diarrheal patient in green, animal in red, healthy carrier in yellow and raw meat in purple. The numbers on connecting lines show the number of allelic difference between two STs. The number in the circle is the ST number.
Prevalence of Virulence Genes among Different STs
There were quite differences in the number of virulence genes among different isolates, ranging from 0 to 12 genes per isolate. Two isolates from diarrheal patients belonging to ST11 (O157:H7) and ST29 (O7:H11) contained 12 and 11 virulence genes, respectively (Figure 1, Table S1).
Genes located on the pathogenicity island OI-122 (nleB, nleE, set/ent, or efa1 (lifA)) were identified in isolates belonging to 33 different STs. Complete OI-122 was detected in 54 isolates belonging to 16 different STs. The most prominent STs were ST40 (15 isolates), ST328 (13 isolates), and ST29 (8 isolates). The most common clonal complexes were CC278 (16 strains), CC40 (15 strains), and CC29 (11 isolates). The presence of the OI-122 genes was associated with different combinations of other virulence genes. Among the 54 isolates harboring the complete OI-122, 51 strains were ehaC positive with 42 strains being lpfA positive, while paa, ehxA, iha, toxB, katP, astA, and cdt were present in 26, 9, 9, 6, 6, 4, and 1 isolates, respectively (Figure 1, Table S1).
Three variants of the lpfA gene (lpfAO113, lpfAO157/OI–154, and lpfAO157/OI–141) were detected in 40 different STs, which contained 124 isolates. The most prominent STs were ST40 (all carrying lpfA), ST328 (all carrying lpfA) and ST29 (all except one carrying the lpfA gene). In contrast, for the 32 isolates contained in CC10, the lpfA gene was only detected in four isolates without any of the OI-122 genes.
The present study is the first detailed molecular profiling study of aEPEC strains from China. Serotyping results showed that 66 serogroups were identified, most of which were non-classical serogroups (Hernandes et al., 2009). The WHO has recognized 12 O serogroups as the classical EPEC serogroups, i.e., O26, O55, O86, O111, O114, O119, O125, O126, O127, O128, O142, and O158. Nine classical serotypes were found in our study, which contained 38 (16.7%) aEPEC isolates. They were not collectively more prevalent than the remaining serogroups. In fact, the most frequently detected serogroup of aEPEC reported so far is O51, followed by O145, O26, O55 and O111, and O119 (Hernandes et al., 2009). Globally, there are various studies showing that aEPEC has different serogrous, such as O76 (Møller-Stray et al., 2012), O127 (Hao et al., 2012), O119, O26, and O88 (Jenkins et al., 2003). In this study, 7.4% of aEPEC isolates were O untypable, which is much lower than that reported by others (Blanco et al., 2006; Monaghan et al., 2013; Zhang et al., 2016). Blanco et al. found that 110 EPEC strains were typed to 43 O serogroups and 69 serotypes with 44 new serotypes and only 13% could be assigned to classical EPEC serotypes (Blanco et al., 2006). So, the vast diversity of serotypes suggested that serotyping alone is insufficient for the detection of aEPEC. To better characterize aEPEC strains, it is necessary to combine the phenotypic and genotypic information.
Virulence gene profiling is essential in determining the pathogenic potential of these strains. In our investigation, statistically meaningful differences were observed among the distributions of OI-122 genes efa1 (lifA), nleB, nleE, set/ent and the phylogenetic marker gene yjaA. Afset et al. observed that efa1 (lifA) was the gene with the strongest statistical association with diarrhea (Afset et al., 2006). The other three genes were also found to be linked to diarrhea, while gene yjaA was negatively associated with diarrhea. Our results showed that OI-122 genes efa1 (lifA), nleB, nleE, set/ent, and the phylogenetic marker gene yjaA were associated with diarrhea among the aEPEC isolates studied, which is consistent with the results detected by Afset and colleagues (Afset et al., 2006). Furthermore, efa1 (lifA), nleB, nleE, and set/ent were found to be more prevalent in aEPEC strains isolated from animals or raw meats than in strains from healthy carriers. These findings indicated that aEPEC strains from different sources might carry different virulence gene profiles. The presence of aEPEC strains carrying OI-122 genes in animals and raw meats might pose a zoonotic risk to humans. This finding also supports the concept that animals are the sources of aEPEC infection in humans (Wang et al., 2013).
The paa gene is frequently detected in porcine EPEC. Paa is involved in intimate attachment of the bacteria to enterocytes and causes post-weaning diarrhea in pigs (Batisson et al., 2003). In accordance with previous studies, the positive association between paa and diarrhea was also observed in our study. Balière et al. (2016) observed a significantly higher (60%) occurrence of paa in EPEC isolates from environment samples, which is similar to the frequency we observed in raw meat isolates (63.6%) in this study. paa was much more frequently detected in strains isolated from raw meats than from other sources, and was especially higher than strains from healthy carriers (14.8%). Similarly, the gene ehxA that encodes the EHEC hemolysin, was also found at a higher frequency among aEPEC from raw meats, although we could not confirm the positive association of ehxA with diarrhea (Afset et al., 2006; Scaletsky et al., 2009) in this study. In a previous study, an outbreak caused by aEPEC following the consumption of contaminated food or water was reported (Wedley et al., 2013). Thus, the high presence of aEPEC possessing paa or ehxA in raw meats indicates that these isolates could cause disease if ingested by humans, and therefore, remains a serious public health threat to consumers.
In a previous study, the virulence marker astA was significantly associated with diarrhea (Dulguer et al., 2003). In Japan, a waterborne outbreak of diarrhea was found to be associated with aEPEC strains carrying the astA gene (Yatsuyanagi et al., 2003). However, a negative statistical association with diarrhea has also been reported (Afset et al., 2006). In the current investigation assessing 228 isolates, no correlation was confirmed, in agreement with another study (Scaletsky et al., 2009).
Autotransporter (AT) proteins are often associated with virulence attributes such as adherence, biofilm formation, and toxicity and these proteins have been identified in both tEPEC and aEPEC (Abreu et al., 2013). Several studies have investigated the presence of AT genes in EPEC (Abreu et al., 2013; Balière et al., 2016). In the present study, four AT genes (ehaC, espP, pet, and pic) were detected. Among these genes, the most frequent gene observed in our isolates was ehaC, which was detected at frequencies similar to those observed by Abreu and colleagues (Abreu et al., 2013). However, it was found more often in isolates from diarrheal patients, animals and raw meats than in isolates from healthy carriers. Our hypothesis is that gene ehaC might play a role in aEPEC infection.
Afset et al. assessed the phylogenetic group and association with diarrhea of 56 aEPEC strains and found that phylogenetic groups B1 and D were weakly associated with diarrhea (Afset et al., 2008). Wang et al. observed that phylogroup B1 was detected more frequently among patient strains than healthy carrier strains (Wang et al., 2013). Furthermore, they also demonstrated that the virulence group Ia comprising the efa1 (lifA) gene and/or lpfA genes but without the yjaA gene, was also detected more frequently in clinical isolates compared with healthy carriers. A gene cluster encoding enzymes involved in O-antigen synthesis was found to be significantly associated with lethal infections (Donnenberg et al., 2015). The pathogenesis of aEPEC also seems to be related to the serotypes (Hu and Torres, 2015). In our investigation, a different method of MLST was employed for phylogenetic group determination. There was a trend for the retention of complete OI-122 and the relationship with specific STs or clonal complexes and serotypes. Genes of OI-122 or lpfA were more frequently detected in ST40 strains with O119:H21 or O109, ST328 strains with O51:H7 or O88:H25, ST795 strains with O76:H7, and ST29 or ST21 strains with O26. These aEPEC strains that were isolated from diarrheal patients, animals and raw meats harbored genes reported to be associated with human diarrhea and might be considered more potentially pathogenic. It has been reported that many virulence genes are present on mobile elements, such as, plasmids, PAIs, or phages and are highly interchangeable among different bacterial strains through horizontal transfer (Hacker and Kaper, 2000). The results that some virulence genes are present more often in specific STs and serotypes than other members suggest that some STs and serotypes of aEPEC might have different genomic backgrounds which permit the acquisition and expression of virulence factors. However, the reason for this is unknown. The detection of such aEPEC strains by subgroup demands further attention. Molecular assessment combining phylogenetic analysis and virulence gene profiling might be employed to establish the pathogenic potential of aEPEC strains.
However, in the present study, there were fewer isolates from raw meats and healthy carriers than from diarrheal patients. Unequal sample sizes and bias in data collection might exist. Nevertheless, this was the most comprehensive data analyzing aEPEC isolates in China. Also, the results did show the vast diversity of the aEPEC isolates from different sources in China.
In conclusion, aEPEC strains in China are heterogeneous in serotypes, virulence gene profiles and phylogenies. It is possible that different aEPEC strains might have different pathogenic potential. The aEPEC strains which were more potentially pathogenic seemed to be related to specific STs or clonal complexes and serotypes.
YXu and YX designed the project; YJ, BH, HW, RF, and SF carried out the sampling work; YXu, XB, SF, and RF carried out the experiments and generated data; YXu and YX analyzed data and drafted the manuscript. All authors have read and approved the final version of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the National Natural Science Foundation of China , the State Key Laboratory for Infectious Disease Prevention and Control [2015SKLID504], and the National Basic Research Priorities Program of China [2015CB554201].
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fcimb.2017.00109/full#supplementary-material
Abreu, A. G., Abe, C. M., Nune, K. O., Moraes, C. T., Chavez-Dueñas, L., Navarro-Garcia, F., et al. (2016). The serine protease Pic as a virulence factor of atypical enteropathogenic Escherichia coli. Gut. Microb. 7, 115–125. doi: 10.1080/19490976.2015.1136775
Abreu, A. G., Bueris, V., Porangaba, T. M., Sircili, M. P., Navarro-Garcia, F., and Elias, W. P. (2013). Autotransporter protein-encoding genes of diarrheagenic Escherichia coli are found in both typical and atypical enteropathogenic E. coli strains. Appl. Environ. Microbiol. 79, 411–414. doi: 10.1128/AEM.02635-12
Abu-Median, A. B., van Diemen, P. M., Dziva, F., Vlisidou, I., Wallis, T. S., and Stevens, M. P. (2006). Functional analysis of lymphostatin homologues in enterohaemorrhagic Escherichia coli. FEMS Microbiol. Lett. 258, 43–49. doi: 10.1111/j.1574-6968.2006.00198.x
Afset, J. E., Anderssen, E., Bruant, G., Harel, J., Wieler, L., and Bergh, K. (2008). Phylogenetic backgrounds and virulence profiles of atypical enteropathogenic Escherichia coli strains from a case-control study using multilocus sequence typing and DNA microarray analysis. J. Clin. Microbiol. 46, 2280–2290. doi: 10.1128/JCM.01752-07
Afset, J. E., Bruant, G., Brousseau, R., Harel, J., Anderssen, E., Bevanger, L., et al. (2006). Identification of virulence genes linked with diarrhea due to atypical enteropathogenic Escherichia coli by DNA microarray analysis and PCR. J. Clin. Microbiol. 44, 3703–3711. doi: 10.1128/JCM.00429-06
Bai, X., Zhao, A., Lan, R., Xin, Y., Xie, H., Meng, Q., et al. (2013). Shiga toxin-producing Escherichia coli in yaks (Bos grunniens) from the Qinghai-Tibetan Plateau, China. PLoS ONE 8:e65537. doi: 10.1371/journal.pone.0065537
Balière, C., Rincé, A., Delannoy, S., Fach, P., and Gourmelon, M. (2016). Molecular profiling of Shiga toxin-producing Escherichia coli and Enteropathogenic E. coli strains isolated from French coastal environments. Appl. Environ. Microbiol. 82, 3913–3927. doi: 10.1128/AEM.00271-16
Bando, S. Y., Andrade, F. B., Guth, B. E., Elias, W. P., Moreira-Filho, C. A., and Pestana De Castro, A. F. (2009). Atypical enteropathogenic Escherichia coli genomic background allows the acquisition of non-EPEC virulence factors. FEMS Microbiol. Lett. 299, 22–30. doi: 10.1111/j.1574-6968.2009.01735.x
Bardiau, M., Labrozzo, S., and Mainil, J. G. (2009). Putative adhesins of enteropathogenic and enterohemorrhagic Escherichia coli of serogroup O26 isolated from humans and cattle. J. Clin. Microbiol. 47, 2090–2096. doi: 10.1128/JCM.02048-08
Batisson, I., Guimond, M. P., Girard, F., An, H., Zhu, C., Oswald, E., et al. (2003). Characterization of the novel factor paa involved in the early steps of the adhesion mechanism of attaching and effacing Escherichia coli. Infect. Immun. 71, 4516–4525. doi: 10.1128/IAI.71.8.4516-4525.2003
Blanco, M., Blanco, J. E., Dahbi, G., Alonso, M. P., Mora, A., Coira, M. A., et al. (2006). Identification of two new intimin types in atypical enteropathogenic Escherichia coli. Int. Microbiol. 9, 103–110.
Brunder, W., Schmidt, H., and Karch, H. (1996). KatP, a novel catalase-peroxidase encoded by the large plasmid of enterohaemorrhagic Escherichia coli O157:H7. Microbiology 142(Pt 11), 3305–3315. doi: 10.1099/13500872-142-11-3305
Brunder, W., Schmidt, H., Frosch, M., and Karch, H. (1999). The large plasmids of Shiga-toxin-producing Escherichia coli (STEC) are highly variable genetic elements. Microbiology 145, 10. doi: 10.1099/13500872-145-5-1005
Clermont, O., Bonacorsi, S., and Bingen, E. (2000). Rapid and simple determination of the Escherichia coli phylogenetic group. Appl. Environ. Microbiol. 66, 4555–4558. doi: 10.1128/AEM.66.10.4555-4558.2000
Croxen, M. A., Law, R. J., Scholz, R., Keeney, K. M., Wlodarska, M., and Finlay, B. B. (2013). Recent advances in understanding enteric pathogenic Escherichia coli. Clin. Microbiol. Rev. 26, 822–880. doi: 10.1128/CMR.00022-13
Donnenberg, M. S., Hazen, T. H., Farag, T. H., Panchalingam, S., Antonio, M., Hossain, A., et al. (2015). Bacterial factors associated with lethal outcome of enteropathogenic Escherichia coli infection: genomic case-control studies. PLoS Negl. Trop. Dis. 9:e0003791. doi: 10.1371/journal.pntd.0003791
Dulguer, M. V., Fabbricotti, S. H., Bando, S. Y., Moreira-Filho, C. A., Fagundes-Neto, U., and Scaletsky, I. C. (2003). Atypical enteropathogenic Escherichia coli strains: phenotypic and genetic profiling reveals a strong association between enteroaggregative E. coli heat-stable enterotoxin and diarrhea. J. Infect. Dis. 188, 1685–1694. doi: 10.1086/379666
Fields, P. I., Blom, K., Hughes, H. J., Helsel, L. O., Feng, P., and Swaminathan, B. (1997). Molecular characterization of the gene encoding H antigen in Escherichia coli and development of a PCR-restriction fragment length polymorphism test for identification of E. coli O157:H7 and O157:NM. J. Clin. Microbiol. 35, 1066–1070.
Hao, R., Qiu, S., Wang, Y., Yang, G., Su, W., Song, L., et al. (2012). Quinolone-resistant Escherichia coli O127a:K63 serotype with an extended-spectrum-beta-lactamase phenotype from a food poisoning outbreak in China. J. Clin. Microbiol. 50, 2450–2451. doi: 10.1128/JCM.00276-12
Hernandes, R. T., Elias, W. P., Vieira, M. A., and Gomes, T. A. (2009). An overview of atypical enteropathogenic Escherichia coli. FEMS Microbiol. Lett. 297, 137–149. doi: 10.1111/j.1574-6968.2009.01664.x
Iguchi, A., Iyoda, S., Seto, K., Morita-Ishihara, T., Scheutz, F., Ohnishi, M., et al. (2015). Escherichia coli O-Genotyping PCR: a comprehensive and practical platform for molecular O serogrouping. J. Clin. Microbiol. 53, 2427–2432. doi: 10.1128/JCM.00321-15
Jenkins, C., Lawson, A. J., Cheasty, T., Willshaw, G. A., Wright, P., Dougan, G., et al. (2003). Subtyping intimin genes from enteropathogenic Escherichia coli associated with outbreaks and sporadic cases in the United Kingdom and Eire. Mol. Cell. Probes 17, 149–156. doi: 10.1016/S0890-8508(03)00046-X
Jinadasa, R. N., Bloom, S. E., Weiss, R. S., and Duhamel, G. E. (2011). Cytolethal distending toxin: a conserved bacterial genotoxin that blocks cell cycle progression, leading to apoptosis of a broad range of mammalian cell lineages. Microbiology 157, 1851–1875. doi: 10.1099/mic.0.049536-0
Joensen, K. G., Tetzschner, A. M., Iguchi, A., Aarestrup, F. M., and Scheutz, F. (2015). Rapid and easy in silico serotyping of Escherichia coli isolates by use of whole-genome sequencing data. J. Clin. Microbiol. 53, 2410–2426. doi: 10.1128/JCM.00008-15
Kobayashi, N., Lee, K., Yamazaki, A., Saito, S., Furukawa, I., Kono, T., et al. (2013). Virulence gene profiles and population genetic analysis for exploration of pathogenic serogroups of Shiga toxin-producing Escherichia coli. J. Clin. Microbiol. 51, 4022–4028. doi: 10.1128/JCM.01598-13
Kotloff, K. L., Nataro, J. P., Blackwelder, W. C., Nasrin, D., Farag, T. H., Panchalingam, S., et al. (2013). Burden and aetiology of diarrhoeal disease in infants and young children in developing countries (the Global Enteric Multicenter Study, GEMS): a prospective, case-control study. Lancet 382, 209–222. doi: 10.1016/S0140-6736(13)60844-2
Maldonado, N. A., Múnera, M. I., López, J. A., Sierra, P., Robledo, C., Robledo, J., et al. (2014). Trends in antibiotic resistance in Medellin and municipalities of the Metropolitan area between 2007 and 2012: results of six years of surveillance. Biomedica 34, 433–446. doi: 10.7705/biomedica.v34i3.1658
Maluta, R. P., Fairbrother, J. M., Stella, A. E., Rigobelo, E. C., Martinez, R., and de Ávila, F. A. (2014). Potentially pathogenic Escherichia coli in healthy, pasture-raised sheep on farms and at the abattoir in Brazil. Vet. Microbiol. 169, 89–95. doi: 10.1016/j.vetmic.2013.12.013
Mercado, E. H., Piscoche, C., Contreras, C., Durand, D., Riveros, M., Ruiz, J., et al. (2016). Pathogenicity island O-122 in enteropathogenic Escherichia coli strains is associated with diarrhea severity in children from Lima Peru. Int. J. Med. Microbiol. 306, 231–236. doi: 10.1016/j.ijmm.2016.05.005
Møller-Stray, J., Eriksen, H. M., Bruheim, T., Kapperud, G., Lindstedt, B. A., Skeie, Å., et al. (2012). Two outbreaks of diarrhoea in nurseries in Norway after farm visits, April to May 2009. Euro Surveill. 17, 20321–20327.
Monaghan, Á., Byrne, B., Fanning, S., Sweeney, T., McDowell, D., and Bolton, D. J. (2013). Serotypes and virulence profiles of atypical enteropathogenic Escherichia coli (EPEC) isolated from bovine farms and abattoirs. J. Appl. Microbiol. 114, 595–603. doi: 10.1111/jam.12064
Nataro, J. P., Seriwatana, J., Fasano, A., Maneval, D. R., Guers, L. D., Noriega, F., et al. (1995). Identification and cloning of a novel plasmid-encoded enterotoxin of enteroinvasive Escherichia coli and Shigella strains. Infect. Immun. 63, 4721–4728.
Newton, H. J., Pearson, J. S., Badea, L., Kelly, M., Lucas, M., Holloway, G., et al. (2010). The type III effectors NleE and NleB from enteropathogenic E. coli and OspZ from Shigella block nuclear translocation of NF-κB p65. PLoS Pathog 6:e1000898. doi: 10.1371/journal.ppat.1000898
Nicholls, L., Grant, T. H., and Robins-Browne, R. M. (2000). Identification of a novel genetic locus that is required for in vitro adhesion of a clinical isolate of enterohaemorrhagic Escherichia coli to epithelial cells. Mol. Microbiol. 35, 275–288. doi: 10.1046/j.1365-2958.2000.01690.x
Ooka, T., Terajima, J., Kusumoto, M., Iguchi, A., Kurokawa, K., Ogura, Y., et al. (2009). Development of a multiplex PCR-based rapid typing method for enterohemorrhagic Escherichia coli O157 strains. J. Clin. Microbiol. 47, 2888–2894. doi: 10.1128/JCM.00792-09
Paton, A. W., and Paton, J. C. (2002). Direct detection and characterization of Shiga toxigenic Escherichia coli by multiplex PCR for stx1, stx2, eae, ehxA, and saa. J. Clin. Microbiol. 40, 271–274. doi: 10.1128/JCM.40.1.271-274.2002
Paton, A. W., Srimanote, P., Woodrow, M. C., and Paton, J. C. (2001). Characterization of Saa, a novel autoagglutinating adhesin produced by locus of enterocyte effacement-negative Shiga-toxigenic Escherichia coli strains that are virulent for humans. Infect. Immun. 69, 6999–7009. doi: 10.1128/IAI.69.11.6999-7009.2001
Patzi-Vargas, S., Zaidi, M. B., Perez-Martinez, I., León-Cen, M., Michel-Ayala, A., Chaussabel, D., et al. (2015). Diarrheagenic Escherichia coli carrying supplementary virulence genes are an important cause of moderate to severe diarrhoeal disease in Mexico. PLoS Negl. Trop. Dis. 9:e0003510. doi: 10.1371/journal.pntd.0003510
Ruiz, R. C., Melo, K. C., Rossato, S. S., Barbosa, C. M., Corrêa, L. M., Elias, W. P., et al. (2014). Atypical enteropathogenic Escherichia coli secretes plasmid encoded toxin. Biomed Res. Int. 2014:896235. doi: 10.1155/2014/896235
Scaletsky, I. C., Aranda, K. R., Souza, T. B., Silva, N. P., and Morais, M. B. (2009). Evidence of pathogenic subgroups among atypical enteropathogenic Escherichia coli strains. J. Clin. Microbiol. 47, 3756–3759. doi: 10.1128/JCM.01599-09
Schmidt, H., Zhang, W. L., Hemmrich, U., Jelacic, S., Brunder, W., Tarr, P. I., et al. (2001). Identification and characterization of a novel genomic island integrated at selC in locus of enterocyte effacement-negative, Shiga toxin-producing Escherichia coli. Infect. Immun. 69, 6863–6873. doi: 10.1128/IAI.69.11.6863-6873.2001
Sheikh, J., Czeczulin, J. R., Harrington, S., Hicks, S., Henderson, I. R., Le Bouguénec, C., et al. (2002). A novel dispersin in enteroaggregative Escherichia coli. J. Clin. Invest. 110, 1329–1337. doi: 10.1172/JCI16172
Szalo, I. M., Goffaux, F., Pirson, V., Piérard, D., Ball, H., and Mainil, J. (2002). Presence in bovine enteropathogenic (EPEC) and enterohaemorrhagic (EHEC) Escherichia coli of genes encoding for putative adhesins of human EHEC strains. Res. Microbiol. 153, 653–658. doi: 10.1016/S0923-2508(02)01379-7
Tarr, C. L., Large, T. M., Moeller, C. L., Lacher, D. W., Tarr, P. I., Acheson, D. W., et al. (2002). Molecular characterization of a serotype O121:H19 clone, a distinct Shiga toxin-producing clone of pathogenic Escherichia coli. Infect. Immun. 70, 6853–6859. doi: 10.1128/IAI.70.12.6853-6859.2002
Tarr, P. I., Bilge, S. S., Vary, J. C. Jr., Jelacic, S., Habeeb, R. L., Ward, T. R., et al. (2000). Iha: a novel Escherichia coli O157:H7 adherence-conferring molecule encoded on a recently acquired chromosomal island of conserved structure. Infect. Immun. 68, 1400–1407. doi: 10.1128/IAI.68.3.1400-1407.2000
Tatsuno, I., Horie, M., Abe, H., Miki, T., Makino, K., Shinagawa, H., et al. (2001). toxB gene on pO157 of enterohemorrhagic Escherichia coli O157:H7 is required for full epithelial cell adherence phenotype. Infect. Immun. 69, 6660–6669. doi: 10.1128/IAI.69.11.6660-6669.2001
Toma, C., Martínez Espinosa, E., Song, T., Miliwebsky, E., Chinen, I., Iyoda, S., et al. (2004). Distribution of putative adhesins in different seropathotypes of Shiga toxin-producing Escherichia coli. J. Clin. Microbiol. 42, 4937–4946. doi: 10.1128/JCM.42.11.4937-4946.2004
Wang, L., Wakushima, M., Aota, T., Yoshida, Y., Kita, T., Maehara, T., et al. (2013). Specific properties of enteropathogenic Escherichia coli isolates from diarrheal patients and comparison to strains from foods and fecal specimens from cattle, swine, and healthy carriers in Osaka City, Japan. Appl. Environ. Microbiol. 79, 1232–1240. doi: 10.1128/AEM.03380-12
Wedley, A. L., Elajnef, H. M., and Fletcher, J. N. (2013). Characterization of a novel EAST-negative enteropathogenic E. coli strain implicated in a food-borne outbreak of diarrhoea in adults. APMIS 121, 494–502. doi: 10.1111/apm.12024
Xu, Y., Bai, X., Zhao, A., Zhang, W., Ba, P., Liu, K., et al. (2016). Genetic diversity of intimin gene of atypical enteropathogenic Escherichia coli isolated from human, animals and raw meats in China. PLoS ONE 11:e0152571. doi: 10.1371/journal.pone.0152571
Yatsuyanagi, J., Saito, S., Miyajima, Y., Amano, K., and Enomoto, K. (2003). Characterization of atypical enteropathogenic Escherichia coli strains harboring the astA gene that were associated with a waterborne outbreak of diarrhea in Japan. J. Clin. Microbiol. 41, 2033–2039. doi: 10.1128/JCM.41.5.2033-2039.2003
Keywords: E. coli, EPEC, serotyping, MLST, virulence gene
Citation: Xu Y, Bai X, Jin Y, Hu B, Wang H, Sun H, Fan R, Fu S and Xiong Y (2017) High Prevalence of Virulence Genes in Specific Genotypes of Atypical Enteropathogenic Escherichia coli. Front. Cell. Infect. Microbiol. 7:109. doi: 10.3389/fcimb.2017.00109
Received: 21 January 2017; Accepted: 17 March 2017;
Published: 04 April 2017.
Edited by:Philip R. Hardwidge, Kansas State University, USA
Reviewed by:Shannon D. Manning, Michigan State University, USA
Jorge Blanco, Universidade de Santiago de Compostela, Spain
Copyright © 2017 Xu, Bai, Jin, Hu, Wang, Sun, Fan, Fu and Xiong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yanwen Xiong, firstname.lastname@example.org