- 1Laboratório de Enterobactérias, Departamento de Microbiologia, Imunologia e Parasitologia, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
- 2Laboratório Experimental de Patogenicidade de Enterobactérias, Departamento de Microbiologia, Imunologia e Parasitologia, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
- 3Laboratório de Biologia Molecular de Micobactérias, Departamento de Microbiologia, Imunologia e Parasitologia, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
Despite being in the era of advanced technology, the world still suffers from old infectious diseases, both intestinal and extraintestinal, where Escherichia coli plays a major role as the etiological agent. Atypical enteropathogenic E. coli (aEPEC) is one of six intestinal pathogenic E. coli pathotypes and one of the major agents causing diarrhea in low- and middle-income countries like Brazil. In this work, we have investigated to what extent a collection of aEPEC isolated from the intestinal tract of children has incorporated virulence traits involved in the development of extraintestinal infections. The phylogenetic origin and the presence of extraintestinal pathogenic E. coli (ExPEC) -related Pathogenicity Islands (PAIs) were evaluated by PCR for a collection of 111 aEPEC isolated from stool. Additionally, they were screened by PCR for the presence of specific ExPEC virulence factors. Phenotypically evaluated for bacteriocin and hemolysin production and assessed for serum resistance. Finally, four strains were sequenced and had their genome characterized. Most of the strains originated from phylogroup B1 (48.6%) and A (36.3%), followed by groups B2 (13.5%) and E (8.1%). About half of the aEPEC strains presented markers for pathogenicity islands originally described in uropathogenic E. coli (UPEC), PAI IV536 being the most prevalent. Many aEPEC strains presented the virulence genetic markers that are the hallmark of ExPEC. Besides, many strains produced bacteriocins and hemolysins and survived in human serum. Five strains fulfilled the molecular criteria to be classified as ExPEC and one as UPEC, highlighting the existence of hybrid genotypes among aEPEC strains. Three non-phylogenetic-related hybrid strains were chosen for further experiments. These strains were lethal in the Galleria mellonella model for ExPEC virulence, and the comparative analysis of their genomes revealed they belong to different EPEC/ EHEC global clonal groups. Overall, this study reports the presence of many attributes of virulence of ExPEC in a comprehensive collection of aEPEC strains. The data presented here indicate the existence of genotypic hybrid aEPEC/ExPEC and aEPEC/UPEC pathogens, suggesting that they can express both intestinal and extraintestinal virulence determinants in humans. Therefore, the consequences of their colonization and infection are more concerning and potentially life-threatening.
1 Introduction
Diarrhea is still the second leading cause of death in children under 5 years old (WHO, 2017). According to the World Health Organization (WHO), in the past, diarrhea-related deaths were due to severe dehydration and fluid loss; however, currently, septic bacterial infections are pointed out to account for an increasing proportion of all diarrhea-associated deaths (WHO, 2017).
The mortality and morbidity rates of Gram-negative septic infections are high, and E. coli is a leading cause of bloodstream infections worldwide. Some virulence factors (VFs) present in mobile genetic elements can be incorporated into diverse genomic backgrounds, enabling pathogens to cause different diseases. In this sense, the emergence of hetero-pathogenic or hybrid pathogenic E. coli emphasizes the possibility of some strains causing diarrhea and systemic infections in the same host (Mariani-Kurkdjian et al., 2014; Kessler et al., 2015; Mandomando et al., 2020; Santos et al., 2020a). Therefore, a portion of the septic infections related to diarrhea caused by E. coli can be associated with strains that have caused diarrhea and then, taking advantage of the intestinal barrier impairments, translocate and reach extraintestinal sites, causing bloodstream infection in the host. This pathogen can be either an Intestinal Pathogenic E. coli (IPEC) that received VFs from Extraintestinal Pathogenic E. coli (ExPEC) or vice versa.
Enteropathogenic Escherichia coli (EPEC) is one of the six IPEC pathotypes and is a leading cause of diarrhea in low- to middle-income countries like Brazil (Afset et al., 2004; Aranda et al., 2007; Gomes et al., 2016).
The EPEC pathotype causes diarrhea due to the presence of a pathogenicity island named Locus of enterocyte effacement (LEE) (Denamur et al., 2021). The LEE contains several genes that encode virulence factors (VFs), which work in coordination with additional genes in the genome to induce bacterial intimate adherence to enterocytes, the cell’s actin remodeling, and the rise of pedestal-like structures below the adhered bacteria, culminating in the microvilli elimination in the region surrounding the bacterial attachment (Trabulsi et al., 2002; Hernandes et al., 2009; Gomes et al., 2016). The EPEC pathotype can be subdivided into two subgroups based on the expression of an adhesin responsible for facilitating adherence onto the enterocyte, named bundle-forming pilus (BFP), which is present in strains classified as typical EPEC (tEPEC) and absent in those classified as atypical EPEC (aEPEC) (Trabulsi et al., 2002; Kaper et al., 2004; Gomes et al., 2016). In Brazil, aEPEC strains represent about 50% of the diarrhea caused by E. coli (Ori et al., 2019). Moreover, aEPEC pathotype genome is recognized to be permissive to the incorporation of DNA by horizontal gene transfer (HGT) events (Afset et al., 2008), which leads to a very heterogeneous array of virulence genetic markers displayed by the aEPEC strains. Recent reports have highlighted aEPEC isolates from extraintestinal infections, including bloodstream infections that may or may not be associated with diarrhea (Bratoeva et al., 1994; Kessler et al., 2015; Valiatti et al., 2020; Nascimento et al., 2021, 2022).
In this work, we evaluated diverse aspects related to the ExPEC pathotype in an aEPEC collection isolated from stools in diverse epidemiological studies to search for the presence of ExPEC-recognized virulence traits that could reveal the existence of hybrid intestinal pathogenic strains with the potential to cause both intestinal and extraintestinal infections.
2 Materials and methods
2.1 Bacterial strains
In this study, 111 previously published aEPEC strains isolated during epidemiological studies conducted in Brazil were revisited (Vieira et al., 2001, 2010; Gomes et al., 2004). Seventy-eight strains (70.3%) were isolated from diarrhea cases and 33 (29.7%) from asymptomatic individuals. The EPEC pathotype was determined by the presence of the eae (intimin) gene and the absence of all the other genes related to any other IPEC pathotype, while the atypical subclassification was determined by the lack of the bfpA (BFP) gene (Vieira et al., 2001, 2010; Gomes et al., 2004).
2.2 Bacteriocin production
The production of bacteriocins was assessed as described earlier (dos Santos et al., 2017). Briefly: strains were grown on three different solid culture media: Nutrient agar (NA), Tryptone yeast extract (TYE) (Difco Laboratories, Sparks-MD), and TYE plus 1% trypsin (Gibco®) (TYE-T). The indicator strain E. coli C600, grown overnight (ON) in Tryptic Soy broth (TSB) (Difco), was diluted 1:100 in sterile saline solution (0.85% NaCl) and spread on the surface of the TYE, NA, and TYE-T media plates to obtain a bacterial lawn. Three microliters of ON TSB cultures of aEPEC strains were seeded over the indicator lawn in equidistant spots. E. coli strains GC138 and HB101 were used as positive and negative controls, respectively (dos Santos et al., 2017). After ON incubation at 36°C, the strains producing bacteriocin appeared surrounded by an inhibition zone of the indicator strain in TYE and NA. As a counterproof, the presence of trypsin in the TYE-T would abolish the bacteriocin effect, resulting in non-inhibited growth of the indicator strain by the bacteriocin-producer strains.
The strains that were positive in the assay were submitted to PCR to evaluate whether the bacteriocin that they produced belongs to bacteriocin types B, E1, E2, Ia, Ib, M, and microcin V using primers (Supplementary Table S1) and conditions described by Šmajs et al. (2010).
2.3 Hemolytic activity
The hemolytic activity was assayed by inoculating 3 μL of an ON TSB culture in Columbia agar medium (Kasvi, Pinhais, PR) supplemented with 5% (v:v) of washed and defibrinated sheep blood and 10 mM CaCl2 (Beutin et al., 1989). Hemolysis was detected after 3 h and 24 h of incubation at 37°C by a clear zone surrounding the inoculated strain. E. coli strains U4-41 and EDL 933 were the positive controls for alpha-hemolysin and enterohemolysin production, respectively. Klebsiella spp. was the negative control. The presence of the hlyA and ehxA genes was searched by PCR (Supplementary Table S2) for all hemolytic strains as previously described (Schmidt et al., 1995; Yamamoto et al., 1995; Johnson and Stell, 2000).
2.4 Serum resistance assay
The serum resistance assay was done using a commercially acquired pool of human sera (Sigma-Aldrich, MO, United States). A logarithmic phase culture grown in TSB was washed twice using phosphate-buffered saline (PBS). The bacterial concentration was adjusted to 1 × 107 CFU/mL in PBS. The 100 μL of bacterial suspension was mixed with normal serum (1:1 v) and incubated at 37°C for 24 h. Another bacteria-serum mixture was prepared using Complement-inactivated serum obtained by heating at 56°C for 30 min. At 0 and 24 h, 5 μL of the mixtures were plated onto MacConkey agar to detect strains that survived serum bactericidal activity during the evaluated period. The E. coli K-12 strain MG1655 was used as a serum-sensitive control, and the strain E. coli J96 (ExPEC prototype) was used as a serum-resistant control.
The genes iss and traT, considered involved in serum resistance, were searched by PCR in all aEPEC strains as previously described (Johnson and Stell, 2000; Dezfulian et al., 2003) using specific primers (Supplementary Table S2).
2.5 Phylogenetic origin, ExPEC and UPEC molecular classification, presence of pathogenicity islands, and clonal relationship
The phylogenetic origin of aEPEC strains was determined using the Clermont quadruplex PCR scheme (Clermont et al., 2013) followed by confirmation (for phylogroups A, C, D, and E), which enables the classification of E. coli into one of seven phylogroups (A, B1, B2, C, D, E, and F).
The aEPEC strains were molecularly classified as ExPEC by the presence of at least two genes among papC, afaBCIII, sfaDE, kpsMTII, and iutA (Johnson et al., 2003), and as UPEC when presenting at least three among vat, chuA, fyuA, and yfcV (Spurbeck et al., 2012). The PCR reactions were performed using primers and conditions described in Supplementary Table S3 (Le Bouguenec et al., 1992; Johnson et al., 1997; Johnson and Stell, 2000; Spurbeck et al., 2012).
Eight pathogenicity islands (PAIs) present in ExPEC prototype strains CFT073 (ICFT073 and IICFT073), J96 (IJ96 and IIJ96), and 536 (I536, II536, III536, and IV536) were searched using primers and conditions (Supplementary Table S4) described by Sabaté et al. (2006) with minor modifications (da Silva et al., 2017).
The clonal relationship of aEPEC strains was assessed by RAPD using primers 1247 (AAGAGCCCGT) and 1283 (GCGATCCCCA) as previously described for E. coli (Nielsen et al., 2014). To determine a clonal relationship, the amplification pattern obtained for both primers was evaluated by the BioNumerics program version 7.6.3 (Applied Maths, Sint-Martens-Latem, Belgium). The similarity of 90% among the amplification profiles for both primers was established as a cutoff to determine the clonal relationship among the strains evaluated.
All PCR reactions were performed using GoTaq® Green Master Mix (Promega, WI, United States) and freshly boiled bacteria as templates. The PCR products were evaluated after electrophoresis in agarose gel.
2.6 Characterization of strains presenting the aEPEC/ExPEC and aEPEC/UPEC hybrid genotypes
2.6.1 Whole genome sequencing of the potentially hybrid aEPEC strains
The aEPEC strains 0811–4, 1551–3, 2071–1, and 3712–3, presenting molecular characteristics of the ExPEC pathotype (including UPEC), here designated as hybrid aEPEC, were chosen to have their genome sequenced by the Microbes NG sequencing service1 as previously published (Santos et al., 2021). Briefly, the genomes were sequenced using the Illumina sequencing platform to obtain 2 × 250 bp paired reads. Each genome was assembled using the SPAdes software (version 3.7), and contigs were annotated using Prokka (version 1.11). Genome analyses were conducted at the Center for Genomic Epidemiology (CGE)2 using their services for the identification of virulence genes (VirulenceFinder version 2.0) (Joensen et al., 2014; Malberg Tetzschner et al., 2020), serotype (SeroTypeFinder version 2.0) (Camacho et al., 2009; Joensen et al., 2015), antibiotic resistance genes (ResFinder version 4.1) (Zankari et al., 2017; Bortolaia et al., 2020), plasmids (PlasmidFinder version 2.0) (Carattoli et al., 2014), fumC and fimH typing (CHTyper version 1.0) (Roer et al., 2018), and sequence type determination (MLST version 2.0) (Wirth et al., 2006), following the Warwick scheme.
The presence of SNPs and number of SNP clusters in clonal-related strains were evaluated using the pathogen detection tool at the National Database of Antibiotic Resistant Organisms (NDARO)3 (Cherry, 2017).
2.6.2 Clonal relationship between hybrid aEPEC strains and EPEC/EHEC clonal groups
Two phylogenetic trees were built to evaluate the clonal relationship of the hybrid aEPEC/ExPEC strains. One tree, using E. coli genomes previously published (Hazen et al., 2016; Ingle et al., 2016; Hernandes et al., 2020), was used to identify to which EPEC/EHEC clonal group the hybrid aEPEC strains belong. The other tree was built using the genomes of the aEPEC strains 0811–4, 1551–3, and 3712–3 (this study), and 93 E. coli published genomes as previously described (Santos et al., 2021). The 93 E. coli public genomes were selected using the Similar Genome Finder Service and the Phylogenetic Tree Building Service at Bacterial and Viral Bioinformatics Resource Center–BV-BRC4 (formerly PATRIC) (Wattam et al., 2017) and details of the strains used are available in Supplementary Table S5. The genomes were downloaded from NCBI and analyzed using VirulenceFinder to determine the virulence genes considered hallmarks of EPEC, STEC, ExPEC, and UPEC. The trees’ final layout was finalized in iTOL v6 (Letunic and Bork, 2021).
2.7 Virulence in the Galleria mellonella model
The virulence of aEPEC strains carrying ExPEC characteristics was investigated “in vivo” using the Galleria mellonella model following procedures previously described (Santos et al., 2020b). An inoculum of 10 μL of the bacterial suspension containing 1 × 107 CFU/mL was injected into the larvae prolegs using a hypodermic insulin syringe with a 31-gauge needle (Uniqmed, Chungcheongnam-do, KOR). The CFU injected inoculum was confirmed by plate count. The tests were performed in three independent assays, using five larvae for each strain tested. The E. coli strains J96 and MG1655 were used as positive and negative controls, respectively, and PBS injections were used as the procedure control. The E. coli JPN15, an EAF plasmid-cured EPEC strain derivative from O127: H6 strain E2348/69 (Levine et al., 1985) was used as an aEPEC control.
2.8 Statistical analyses
The Kaplan–Meier survival curve was used for survival analysis, and the differences between the groups were determined by the log-rank (Mantel-Cox) and Gehan-Breslow-Wilcoxon tests. The threshold for statistical significance was a p-value < 0.05. All analyses were performed using GraphPad Prism 9.2 (GraphPad Software, LLC).
3 Results
3.1 Frequency of bacteriocinogenic and hemolytic strains among aEPEC
In this work, we have searched for the production of bacteriocins, hemolysins, and serum resistance, which are recognized as contributors to the ExPEC virulence pathotype (Johnson and Russo, 2018) in aEPEC strains.
Production of bacteriocins in TYE and NA was detected in 24 (21.6%) aEPEC strains, with results of both media in accordance. All were confirmed producers by the test using the same media supplemented with trypsin. The bacteriocinogenic strains were submitted to molecular typing, and the bacteriocin of 16 (67%) were successfully genetically typed. The bacteriocins Ia and Ib were the most frequently identified and occurred together in six strains, while microcin V was identified in only one strain (Table 1). The bacteriocin types E1 and E2 were identified in three strains each, and the bacteriocin M was identified in another three strains, co-occurring with bacteriocin B in two of them (Table 1).
Twenty-one (18.9%) aEPEC strains were hemolytic when cultivated in Blood agar. For three aEPEC strains, hemolysis was detected after 3 h of incubation and in the others after 6 or 18 h of incubation. The hemolytic strains were evaluated by PCR for the presence of the genes hlyA (α-hemolysin) and ehxA (enterohemolysin), confirming the occurrence of one of these genes in 19 strains. Seventeen of the hemolytic aEPEC strains harbored hlyA and two ehxA. The remaining two strains, although hemolytic, were not detectable by the molecular method used.
3.2 aEPEC strains can survive in human sera
The ability to survive the action of the Complement system active in normal human serum is a main virulence property that enables extraintestinal pathogens to succeed in passing through the bloodstream and gain access to extraintestinal sites.
The majority of aEPEC strains (n = 95, 85.6%) studied were viable after the 24-h challenge in the presence of human sera at a concentration of 50% in PBS. Among the resistant strains, only 13 survived during the challenge, reaching a low concentration in human sera, while the others kept high concentrations at the end of the assay. Sixteen strains were susceptible to human sera effects and did not survive the challenge. Thirty-eight (34.2%) strains were positive for traT, and only four of them did not display the serum-resistant phenotype. No strain was positive for the genes iss and kpsMTII with the primers used.
3.3 The phylogenetic origin, the putative presence of ExPEC PAIs, and genetic virulence markers in aEPEC
The phylogenetic origin of the 111 aEPEC strains as determined by PCR revealed that phylogroup B1 was the most frequent (54 strains—48.7%), followed by groups A (32 strains—28.8%), B2 (15 strains—13.5%), and E (9 strains—8.1%). Phylogroup D presented only one strain (0.9%). No strain from phylogroups C and F was identified. No difference was detected in phylogenetic distribution comparing strains isolated from cases and controls (Supplementary Table S6).
Searching the PAIs by PCR, it was found that about 55% of the aEPEC strains presented PAIs originally described in UPEC. Table 2 shows the distribution of these PAIs among the aEPEC strains according to their phylogroups. Five out of eight PAI markers searched were found in aEPEC, with only one PAI occurring in each strain, except for PAIs IV536, and ICFT073, which appeared combined in two strains, and PAIs II536, and IJ96 in another two (Table 2). The PAIs most frequently detected were IV536 (38 strains—34.2%), IJ96 (10 strains—9%), and ICFT073 (9 strains—8.1%), comprising 55 out of 61 PAI-positive strains. No strain presented PAI markers for islands IIJ96, I536, and III536.
Virulence factors related to ExPEC intrinsic virulence (afa/dra, sfa/foc, pap, kpsMTII, iuc/iut) and uropathogenic potential (vat, fyuA, chuA, yfcV) were searched. Many aEPEC strains presented at least one of these virulence genetic markers that are considered the hallmark of the capacity of the strains to cause systemic or urinary tract infections in murine models. Although many of these genes were spread in a variety of aEPEC strains (Table 3), only five could be classified as ExPEC (strains 0811–4, 2071–1, 1482–11, 2332–7, and 1551–3), and another one as UPEC (strain 3712–3) (Table 3), highlighting the existence of hybrid genotypes among aEPEC strains.

Table 3. Frequency of virulence genetic markersa, and molecular classification of ExPEC and UPEC in 111 aEPEC strains.
Regarding the presence of VFs associated with PAIs, except for fyuA and PAI IV, the other VFs did not hold this association in most strains. This was pointed out by the absence or rare co-occurrence of pap (related to PAI IICFT073, PAI II536, and PAI IJ96), sfa (related to PAI III536), kps (related to PAI ICFT073), and hlyA (related to PAI IJ96) among the aEPEC strains evaluated.
Table 4 presents a summary of the virulence properties that were detected in the aEPEC strains studied, comparing the frequency of each property by phylogroup.
Not surprisingly, the occurrence of the virulence properties followed the phylogenetic distribution, and most of the traits occurred in phylogroups A and B1, which are the most frequent among aEPEC strains. Serum-resistant strains were identified in all phylogroups, with phylogroup E concentrating the higher frequency of resistant strains than phylogroups A, B1, and B2 (Table 4). The five aEPEC/ExPEC hybrid strains (0811–4, 2071–1, 1482–11, 2332–7, and 1551–3) belonged to phylogroup B1, and the hybrid aEPEC/UPEC (3712–3) belonged to phylogroup B2.
3.4 Clonal relationship and genomic sequencing of hybrid aEPEC strains
The RAPD PCR methodology was used to evaluate the clonal relationship among the aEPEC strains, especially the six aEPEC strains identified as possessing hybrid genotypes, according to the amplification obtained using the primers 1247 and 1283 and considering a cutoff of 90% of similarity for both PCRs. RAPD highlighted the high heterogenicity of the aEPEC strains evaluated, with most not being grouped in a cluster (Supplementary Figure S1).
The hybrid aEPEC strains 0811–4, 1482–11, and 2332–7 are closely related, indicating a clonal relationship. The strain 2071–1 presents a similarity of >95% using primer 1247 but a similarity between 87 and 90% with primer 1283, which raised doubts about its relationship with the other three aEPEC strains. On the other hand, 1551–3 and 3712–3 were not closely related either to each other, or to the first group (Figure 1).

Figure 1. Dendrogram of six aEPEC presenting hybrid genotypes. Amplification profiles obtained in RAPD-PCR using primer 1247 (A) and primer 1283 (B) were evaluated by the BioNumerics and used to build the dendrogram. Note the proximity among strains 0811–4, 2071–1, 1482–11, and 2332–7. Strains 1551–3 and 3712–3 were less correlated. The red dashed line represents the cutoff of 90% similarity in the amplification pattern used to determine the clonality of the strains.
Considering the clonal relationship among 0811–4, 1482–11, and 2332–7, the aEPEC 0811–4 and the more distant strain 2071–1 were chosen to represent this group, being submitted together with 1551–3 and 3712–3 to genome sequencing.
The genome sequences revealed that the strains had predicted genome sizes ranging from 4,784,160 to 5,328,743 bp, with strains 0811–4 and 2071–1 displaying larger predicted genome sizes (Table 5). Three STs were identified among the sequenced strains: 0811–4 and 2071–1, which were clonally related, as suggested by RAPD, belong to ST381; strain 1551–5 belongs to ST29, and strain 3712–3 to ST1092 (Table 5).

Table 5. Genome predicted size, metrics of the WGS, and MLST of aEPEC hybrid strains evaluated in this study.
A phylogenetic analysis conducted with 0811–4 and 2071–1 confirmed that they are clonally related (Supplementary Figure S2). Additionally, the pathogen detection tool showed that they belong to the same SNP cluster PDS000052783.4, differing by only 18 SNPs. Therefore, only 0811–4 was taken as a reference for the four aEPEC/ExPEC strains from ST381 for other evaluations.
A second phylogenetic analysis, including all hybrid strains, showed that these three strains belong to different EPEC/EHEC clonal groups, with 1551–3 belonging to the EHEC 2 clonal group, 0811–4 to EPEC 14, and 3712–3 to EPEC 4 (Figure 2).

Figure 2. EPEC and EHEC clonal lineages. The whole genome sequence of the hybrid strains was evaluated to identify to which EPEC/EHEC lineage they belong. Each one of the hybrids belongs to a different clonal lineage: EPEC 4 (3712–3), EPEC 14 (0811–4), and EHEC 2 (1551–3). The phylogenetic tree was built with the codon tree method, evaluating 1,000 single-copy CDS (amino acid and nucleic acid sequences) in an RA × ML matrix, using the genomes of previously published strains (Hazen et al., 2016; Hernandes et al., 2020) and the hybrids aEPEC/ExPEC (0811–4 and 1551–3) and aEPEC/UPEC (3712–3). Each strip segment represents one of the EPEC or EHEC clonal lineages. The hybrid strains are highlighted in bold. Escherichia fergusonii was used as an out-group to root the tree. Bootstrap above 80 was displayed in the tree.
The assembled genomes were submitted to VFDB, and the VF content was compared (Figure 3). It is noticeable that the VF content of the strains is slightly different, although 0811–4 and 1551–3 shared many VFs. In common, besides the LEE, all strains bear the genes related to the expression of CFA/I, ECP, curli, EhaB, UpaG, and Yersiniabactin (Figure 3). In addition, a diverse combination of auto transporters, including SPATES (Vat, Tsh, and Pic) and iron acquisition systems (Aerobactin, Salmonella iron transport system, and Heme receptor) was found (Figure 3). Additionally, the profile of VF genes used for molecular classification was identified as complete, including all the genes of operons, and, except for 0811–4, the strains did not share the intimin type or serotype. Strain 0811–3 harbors intimin type epsilon-1, while 1551–3 harbors epsilon-2, and 3712–3 intimin type zeta (Table 6).

Figure 3. VFs of hybrid aEPEC identified. The VFs identified using VFDB were used to build this figure. The filled squares represent the presence of the genes, while the empty squares represent the absence. Only genes identified in at least one strain were displayed in this figure.
The finding of hybrid aEPEC strains presenting UPEC or ExPEC molecular signatures drove us to search for similar genomic characteristics in a larger group of aEPEC strains that have their genomes accessible in the NCBI. The similar genomes, identified using MASH/MinHash, helped the construction of an additional phylogenetic tree that was added with the VFs used to molecular define the EPEC, EHEC, ExPEC, and UPEC pathotypes.
The tree’s analyses showed that all EPEC genomes closely related to 3712–3 were also hybrid aEPEC or EHEC bearing at least three of the required VFs to be classified as UPEC (Figure 4), suggesting that at least some of the strains that compose the EPEC4 clonal group (see Figure 2) are potentially extraintestinal pathogens. In contrast, the other two hybrid aEPEC in this study were unique in presenting ExPEC VFs within their clusters. Yet, it is worth mentioning that 0811–4 is clonal related with other three hybrid aEPEC strains (1482–11, 2071–1, and 2332–7) evaluated in the present study, pointing out that the spread of aEPEC/ExPEC hybrid genotypes is a reality among diarrheagenic strains isolated in Brazil.

Figure 4. Phylogenetic tree comparing the VFs of the aEPEC hybrid strains and 93 E. coli strains with similar genomes available at NCBI. All E. coli strains displaying a similar genome with aEPEC/UPEC strain 3712–3 also harbor VFs related to the UPEC molecular classification. The VFs used for strain comparison were those associated with the molecular classification of EPEC (in blue: bfpA, eae, and tir), STEC (in purple: stx), ExPEC (in red: pap, sfa, afa, kpsMT, and iuc), and UPEC (in green: fyuA, chuA, vat, and yfcV) and identified by VirulenceFinder. Bootstrap ≥ 80 were displayed in the tree as red dots. The three hybrid strains sequenced in the present study were highlighted in bold.
3.5 The hybrid aEPEC tend to be virulent in the G. mellonella model
In this study, the G. mellonella model was used to evaluate the pathogenic potential of three of those above-mentioned aEPEC strains that presented genotypic and phenotypic characteristics of ExPEC, meaning that they represented a hybrid class of pathogenic strains.
The results presented in Figure 5 show that all hybrid strains evaluated (aEPEC/ExPEC 0811–4 and 1551–3, and aEPEC/UPEC 3712–3) were capable of killing more larvae than the negative-control strains MG1655 and JPN15, a variant of prototype EPEC strain E2348/69 devoid of the virulence plasmid. However, the survival rate of larvae for the three hybrid strains was higher when compared to the prototype J96, being higher for 1551–3 (p = 0.0012) than for 0811–4 and 3712–3 (p = 0.01 and 0.02, respectively). There was no statistical difference in virulence among the three hybrid strains in this model (p < 0.05). These data suggest that the inoculum of 1 × 105 CFU per larvae can differentiate diarrheagenic EPEC strains from those harboring a hybrid genotype. However, the inoculum used cannot distinguish the genomic diversity of the aEPEC hybrid strains.

Figure 5. Virulence in G. mellonella was assessed for one hybrid strain belonging to each RAPD type using the Kaplan–Meier survival curve. Each larva was injected with 1 × 105 CFU, and incubated at 37°C for 3 days. The ExPEC prototype strain J96 was used as a positive control, and PBS as the procedure control. ****p < 0.0001. Strains JPN15 and MG1655 were used as negative controls.
4 Discussion
Atypical EPEC is the main E. coli pathotype responsible for diarrhea in Brazil (Afset et al., 2004; Aranda et al., 2004; Liebchen et al., 2011; Ori et al., 2019; Beraldo et al., 2023) and other low- and middle-income countries (Jarquin et al., 2022; Kara et al., 2022; Yamani and Elhadi, 2022). Since this pathotype was defined, it has been essentially recognized as an intestinal pathogen. This concept started to be challenged by an increasing number of reports referring to the isolation of IPEC of various pathotypes from the bloodstream.
E. coli is known to be a bacterium permissive for the acquisition and stable maintenance of foreign DNA brought by various mechanisms that mediate horizontal gene transfer (HGT) (Desvaux et al., 2020). Such genome flexibility helped to understand the existence of strains presenting mixed virulence genes that used to be characteristic of specific pathotypes.
In 2015, Kessler et al. (2015) reported one case of an inpatient with diarrhea that evolved to bacteremia with multiorgan dysfunction due to an ExPEC bearing EPEC VFs. Even though it was not the first report regarding the EPEC pathotype causing extraintestinal infections (Bratoeva et al., 1994; Abe et al., 2008; Toval et al., 2014), the study contributed to the increase in the number of studies reporting the presence of EPEC strains in extraintestinal infections (Riveros et al., 2017; Valiatti et al., 2020; Munhoz et al., 2021; Nascimento et al., 2021, 2022; Yousefipour et al., 2023).
The knowledge of those versatile bacteria prompted the scientific community to consider the concept of hybrid pathogenic E. coli strains (Santos et al., 2020a) therefore intensifying the studies for a better understanding of their virulence potential and their role as a threat to public health. The importance of hybrid pathogenic E. coli bearing IPEC and ExPEC VFs relies on the possibility of these diarrheagenic bacteria translocating from the intestinal lumen, crossing the gut barrier, and reaching the bloodstream to cause systemic infection, which poses an increasing risk of patient death (Kessler et al., 2015; Mandomando et al., 2020; Santos et al., 2020a). Considering the WHO reports that inform that most diarrhea-related deaths are due to septic shock (WHO, 2017), it is important to identify what are the pathogens that are associated with both conditions and, therefore, require additional surveillance and better diagnosis to avoid fatal outcomes mainly among children.
In the present work, we studied 111 aEPEC strains isolated from patients and asymptomatic carriers, evaluating aspects related to phylogenetic origin, presence of PAIs and VFs associated with the ExPEC pathotype, to assess the extension of sharing of virulence characteristics between these two pathotypes, as well as the existence of putative hybrid pathogens among aEPEC strains.
As various studies (Hazen et al., 2016; Richter et al., 2018; Hernandes et al., 2020) have already demonstrated, the phylogenetic analyses showed that most of the aEPEC strains belonged to phylogroups B1 (48.6%) and A (28.8%).
As many as 61 (54.9%) strains possessed genetic PAI markers first described in UPEC. The presence of PAI markers in all phylogroups of aEPEC, except for phylogroup D, poorly represented by one strain, suggests that the accessibility of the aEPEC population to HGT is independent of its phylogenetic origin. Diverse PAI markers were identified in the sample studied, with PAI IV536 being the most prevalent (34.2%), and detected mainly in phylogroups B1, and A.
The PAI IV536, also known as High Pathogenicity Island (HPI), harbors the iron uptake system called yersiniabactin. Recent studies demonstrated that yersiniabactin contributes to the intrinsic virulence of extraintestinal pathogens and is involved in the capacity of E. coli to kill mice in a murine model of sepsis (Galardini et al., 2020). Besides, it has been used as a molecular marker to identify E. coli strains that are potentially uropathogenic (Spurbeck et al., 2012).
It has to be pointed out that, in our results, the detection of the PAI markers was not in complete agreement with the detection of the virulence genes they were supposed to carry according to the literature, i.e., hlyA, pap, kpsMTII, and iuc for ICFT073, pap for IICFT073, hlyA and pap for IJ96, and pap for II536 (Dobrindt et al., 2002; Sabaté et al., 2006; Desvaux et al., 2020). The only exception occurred with PAI IV536 which was associated with the genes fyuA and irp2 (Dobrindt et al., 2002; Desvaux et al., 2020) in all positive strains. These findings may be due to the presence of incomplete PAIs or even differences in the content of genes they carry. Studies will be addressed to better understand the relationship between the UPEC PAIs and the aEPEC pathotype.
Although more than 50 VFs have been described so far as playing a role in ExPEC virulence (Johnson and Russo, 2018), different studies have shown that it is possible to use a minimum set of VFs to recognize E. coli strains as potential agents of extraintestinal infections, regardless of the bacteria isolation source (Picard et al., 1999; Johnson et al., 2000, 2006, 2022; Spurbeck et al., 2012; Wijetunge et al., 2015; Mellata et al., 2018; Galardini et al., 2020). Further refined studies evaluating the bacterial virulence in vivo in the murine models for sepsis and urinary tract infection (Johnson et al., 2003; Spurbeck et al., 2012), have demonstrated that two sets of genes would be reliable to molecularly classify E. coli that can cause extraintestinal infections in humans (Johnson et al., 2003; Spurbeck et al., 2012). Accordingly, in this study, aEPEC was positive for at least one of all the gene markers used for this molecular classification, except for kpsMTII, papC, and sfaDE. This demonstrates once again that aEPEC frequently shares genetic features responsible for the pathogenicity of ExPEC. Moreover, six aEPEC reunited the right gene combinations to have genotypes identified as hybrid pathogens: five aEPEC/ExPEC, and one aEPEC/UPEC.
Besides the presence of a variety of ExPEC virulence genes, many aEPEC strains were shown to produce bacteriocins and hemolysins which are virulence traits also described as involved in ExPEC virulence, although not exclusively.
The distribution of the virulence properties evaluated among aEPEC in phylogroups showed that their occurrence followed the phylogenetic distribution, occurring more in phylogroups A and B1, which were the most frequent among aEPEC strains. Additionally, the presence of PAI markers and the identification of hybrid strains in phylogroup B1 reinforce the genomic plasticity of strains from these phylogroups, enabling them to receive additional VFs related to extraintestinal pathogenicity, and therefore, prompting them to cause extraintestinal infections.
The resistance to the action of complement present in the human serum is an essential property for pathogens to survive in the bloodstream (Miajlovic and Smith, 2014). Surprisingly, most aEPEC strains (85.6%) survived the serum resistance assay. This is a frequency similar to that found among strains isolated from bloodstream infections in clinical settings (Miajlovic and Smith, 2014), and this frequency is much higher than the one previously reported for fecal strains (Miajlovic and Smith, 2014).
Among the three genes (kpsMTII, traT, and iss) usually related to serum-resistant phenotype and searched in this survey, only traT was found in less than 36% of the serum-resistant strains, while kpsMTII and iss were not detected at all. In regard to the iss gene, it has to be mentioned that the primers used to search it are specific to type-1 iss sequenced (Dezfulian et al., 2003; Johnson et al., 2008), and therefore, it is not possible to rule out the role of the other variant genes in the serum-resistant phenotype observed among aEPEC strains analyzed. Additionally, some of the aEPEC strains studied harbor genes related to serine-protease autotransporters of Enterobacteriaceae (SPATE) (data not shown). Some of these autotransporter proteins were shown to mediate complement cleavage, resulting in serum resistance (Freire et al., 2022; Correa et al., 2024). Therefore, further studies are required to clarify the resistant phenotype identified in the aEPEC strains studied here.
Overall, considering all the differences between IPEC and ExPEC, including those related to the VFs they harbor (Riley, 2020), the simple passage of IPEC strains through the bloodstream would not pose a problem regarding the disease treatment, given that they are not supposed to possess virulence factors to survive or colonize out of the intestinal tract (Johnson and Russo, 2018; Riley, 2020; Denamur et al., 2021). However, the picture changes if these strains are hybrid pathogenic aEPEC/UPEC and aEPEC/ExPEC like those we have detected in this aEPEC survey. In this case, infected patients, mainly children, would be at higher risk of developing extraintestinal infections after colonization of the intestinal tract, which could be either symptomatic or asymptomatic. Riveros et al. (2017), evaluated bacteremia strains isolated from children under 2 years old in Peru and identified diverse IPEC pathotypes, including aEPEC, among the isolates. So far, we are not aware of reports in Brazil about comprehensive epidemiological studies evaluating the presence of IPEC strains carrying ExPEC virulence genes isolated from the bloodstream of children. Nevertheless, there is a report (Liberatore et al., 2011) of an aEPEC isolate, not genotypically hybrid, but displaying the serum-resistant phenotype, which could translocate the intestinal barrier in an experimental translocation model in rats (Liberatore et al., 2011). After the intestinal translocation, the strain was recovered from diverse organs and caused tissue hypoperfusion (Liberatore et al., 2011). This data indicates that among the strains that present serum resistance but lack known virulence markers, there are strains with extraintestinal pathogenic potential that have not yet been adequately evaluated and can potentially cause extraintestinal infections, especially in young children.
The pathogenic potential of three hybrid strains studied in this work (aEPEC/ExPEC 0811–4, 1551–3, and aEPEC/UPEC 3712–3) was proven in the G. mellonella model, which has been reported as reliable to assess the virulence of ExPEC (Alghoribi et al., 2014; Williamson et al., 2014; Ciesielczuk et al., 2015; Pereira et al., 2018). The three strains were lethal in that model, although not as efficiently as the positive ExPEC control E. coli J96. However, they were significantly more lethal than the non-hybrid aEPEC control, strain JPN15, and the non-pathogenic control, E. coli MG 1655.
The genome comparative analysis of the three hybrid strains 0811–4, 1551–3, and 3712–3, showed that they are not phylogenetically related, and accordingly, they belong to different EPEC/EHEC clonal groups, namely EPEC14, EHEC2, and EPEC4, respectively. Corroborating the above analysis, the three hybrid strains belong to different serotypes and sequence types, and all present a complete PAI LEE - the pathogenicity island of EPEC, but display different variants of intimin and its receptor.
Interestingly, the phylogenetic analyses with additional E. coli genomes showed that all EPEC strains closely related to 3712–3 are hybrid EPEC/UPEC and bear at least three of the required VFs to be classified as UPEC, suggesting that part of the strains that compose the EPEC4 clonal group are potentially extraintestinal pathogens. In accordance, a few strains identified in this cluster were isolated from extraintestinal infections. Curiously, two other strains isolated in Brazil, PA58, and BA1250, previously published (Gioia-Di Chiacchio et al., 2018; Munhoz et al., 2021), also belong to this cluster. Strain BA1250 was isolated from a diarrhea case in a different geographical region (Munhoz et al., 2021), and strain PA58 was a typical EPEC/Shiga-toxin producer E. coli (STEC) hetero-pathogenic strain isolated from a bird (Gioia-Di Chiacchio et al., 2018) that also carries UPEC molecular markers. The other two strains, 0811–4 and 1551–3, are the only strains that bear ExPEC VFs in their clusters (Figure 4). However, it is important to mention that 0811–4 has a clonal relationship with the other three aEPEC strains evaluated in the present study, pointing out that this hybrid genotype is present, and somehow often, among aEPEC strains isolated in São Paulo, Brazil.
Finally, the data herein report the detection of various ExPEC virulence attributes in a comprehensive collection of aEPEC strains. In six out of 111 aEPEC strains studied, it was found the genes fulfilling the molecular criteria for their classification as aEPEC/ExPEC or aEPEC/UPEC hybrid pathogens. In addition, these hybrid pathogens isolated from diarrheal cases were proven to be virulent in an in vivo G. mellonella model for ExPEC, suggesting that they can express both intestinal and extraintestinal virulence determinants in humans, turning the consequences of their infection and colonization more severe, and potentially life-threatening.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA871776/.
Ethics statement
The manuscript presents research on animals that do not require additional ethical approval for their study.
Author contributions
AS: Data curation, Formal analysis, Investigation, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. RSS: Formal analysis, Investigation, Validation, Visualization, Writing – review & editing. MV: Investigation, Validation, Writing – review & editing. CN: Investigation, Supervision, Validation, Writing – review & editing. MR: Formal analysis, Investigation, Validation, Visualization, Writing – review & editing. BG: Investigation, Resources, Supervision, Writing – review & editing. TG: Data curation, Funding acquisition, Resources, Writing – review & editing. RMS: Conceptualization, Data curation, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) thematic research grant number 2017/14821–7. Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) provided a scholarship to RSS under financial code 001 and a postdoctoral fellowship to AS (88882.306532/2018–01). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of interest
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.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1599350/full#supplementary-material
Footnotes
References
Abe, C. M., Salvador, F. A., Falsetti, I. N., Vieira, M. A. M., Blanco, J., Blanco, J. E., et al. (2008). Uropathogenic Escherichia coli (UPEC) strains may carry virulence properties of diarrhoeagenic E. coli. FEMS Immunol. Med. Microbiol. 52, 397–406. doi: 10.1111/j.1574-695X.2008.00388.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., Bevanger, L., Romundstad, P., and Bergh, K. (2004). Association of atypical enteropathogenic Escherichia coli (EPEC) with prolonged diarrhoea. J. Med. Microbiol. 53, 1137–1144. doi: 10.1099/jmm.0.45719-0
Alghoribi, M. F., Gibreel, T. M., Dodgson, A. R., Beatson, S. A., and Upton, M. (2014). Galleria mellonella infection model demonstrates high lethality of ST69 and ST127 Uropathogenic E. coli. PLoS One 9:e101547. doi: 10.1371/journal.pone.0101547
Aranda, K. R. S., Fabbricotti, S. H., Fagundes-Neto, U., and Scaletsky, I. C. A. (2007). Single multiplex assay to identify simultaneously enteropathogenic, enteroaggregative, enterotoxigenic, enteroinvasive and Shiga toxin-producing Escherichia coli strains in Brazilian children. FEMS Microbiol. Lett. 267, 145–150. doi: 10.1111/J.1574-6968.2006.00580.X
Aranda, K. R. S., Fagundes-Neto, U., and Scaletsky, I. C. A. (2004). Evaluation of multiplex PCRs for diagnosis of infection with diarrheagenic Escherichia coli and Shigella spp. J. Clin. Microbiol. 42, 5849–5853. doi: 10.1128/JCM.42.12.5849-5853.2004
Beraldo, L. G., Borges, C. A., Maluta, R. P., Cardozo, M. V., and de Ávila, F. A. (2023). Molecular analysis of enteropathogenic Escherichia coli (EPEC) isolates from healthy food-producing animals and humans with diarrhoea. Zoonoses Public Health 70, 117–124. doi: 10.1111/ZPH.13007
Beutin, L., Montenegro, M. A., Orskov, I., Orskov, F., Prada, J., Zimmermann, S., et al. (1989). Close association of verotoxin (Shiga-like toxin) production with enterohemolysin production in strains of Escherichia coli. J. Clin. Microbiol. 27, 2559–2564. doi: 10.1128/jcm.27.11.2559-2564.1989
Bortolaia, V., Kaas, R. S., Ruppe, E., Roberts, M. C., Schwarz, S., Cattoir, V., et al. (2020). ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 75, 3491–3500. doi: 10.1093/JAC/DKAA345
Bratoeva, M. P., Wolf, M. K., Marks, J. K., and Cantey, J. R. (1994). A case of diarrhea, bacteremia, and fever caused by a novel strain of Escherichia coli. J. Clin. Microbiol. 32, 1383–1386. doi: 10.1128/jcm.32.5.1383-1386.1994
Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC Bioinformatics 10:421. doi: 10.1186/1471-2105-10-421
Carattoli, A., Zankari, E., García-Fernández, A., Voldby Larsen, M., Lund, O., Villa, L., et al. (2014). In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 58, 3895–3903. doi: 10.1128/AAC.02412-14
Cherry, J. L. (2017). A practical exact maximum compatibility algorithm for reconstruction of recent evolutionary history. BMC Bioinformatics 18:127. doi: 10.1186/s12859-017-1520-4
Ciesielczuk, H., Betts, J., Phee, L., Doumith, M., Hope, R., Woodford, N., et al. (2015). Comparative virulence of urinary and bloodstream isolates of extra-intestinal pathogenic Escherichia coli in a galleria mellonella model. Virulence 6, 145–151. doi: 10.4161/21505594.2014.988095
Clermont, O., Christenson, J. K., Denamur, E., and Gordon, D. M. (2013). The Clermont Escherichia coli phylo-typing method revisited: improvement of specificity and detection of new phylo-groups. Environ. Microbiol. Rep. 5, 58–65. doi: 10.1111/1758-2229.12019
Correa, G. B., Freire, C. A., Dibo, M., Huerta-Cantillo, J., Navarro-Garcia, F., Barbosa, A. S., et al. (2024). Plasmid-encoded toxin of Escherichia coli cleaves complement system proteins and inhibits complement-mediated lysis in vitro. Front. Cell. Infect. Microbiol. 14:1327241. doi: 10.3389/FCIMB.2024.1327241
da Silva, L. C., de Mello Santos, A. C., and Silva, R. M. (2017). Uropathogenic Escherichia coli pathogenicity islands and other ExPEC virulence genes may contribute to the genome variability of enteroinvasive E. coli. BMC Microbiol. 17:68. doi: 10.1186/s12866-017-0979-5
Denamur, E., Clermont, O., Bonacorsi, S., and Gordon, D. (2021). The population genetics of pathogenic Escherichia coli. Nat. Rev. Microbiol. 19, 37–54. doi: 10.1038/s41579-020-0416-x
Desvaux, M., Dalmasso, G., Beyrouthy, R., Barnich, N., Delmas, J., and Bonnet, R. (2020). Pathogenicity factors of Genomic Islands in intestinal and Extraintestinal Escherichia coli. Front. Microbiol. 11:2065. doi: 10.3389/fmicb.2020.02065
Dezfulian, H., Batisson, I., Fairbrother, J. M., Lau, P. C. K., Nassar, A., Szatmari, G., et al. (2003). Presence and characterization of Extraintestinal pathogenic Escherichia coli virulence genes in F165-Positive E. coli strains isolated from diseased calves and pigs. J. Clin. Microbiol. 41, 1375–1385. doi: 10.1128/JCM.41.4.1375-1385.2003
Dobrindt, U., Blum-Oehler, G., Nagy, G., Schneider, G., Johann, A., Gottschalk, G., et al. (2002). Genetic structure and distribution of four pathogenicity islands (PAI I536 to PAI IV536) of uropathogenic Escherichia coli strain 536. Infect. Immun. 70, 6365–6372. doi: 10.1128/IAI.70.11.6365-6372.2002
dos Santos, L. F., Biscola, F. T., Gonçalves, E. M., and Guth, B. E. C. (2017). Biofilm formation, invasiveness and colicinogeny in locus of enterocyte and effacement negative O113:H21 Shigatoxigenic Escherichia coli. J. Appl. Microbiol. 122, 1101–1109. doi: 10.1111/JAM.13409
Freire, C. A., Silva, R. M., Ruiz, R. C., Pimenta, D. C., Bryant, J. A., Henderson, I. R., et al. (2022). Secreted autotransporter toxin (sat) mediates innate immune system evasion. Front. Immunol. 13:433. doi: 10.3389/fimmu.2022.844878
Galardini, M., Clermont, O., Baron, A., Busby, B., Dion, S., Schubert, S., et al. (2020). Major role of iron uptake systems in the intrinsic extra-intestinal virulence of the genus Escherichia revealed by a genome-wide association study. PLoS Genet. 16:e1009065. doi: 10.1371/journal.pgen.1009065
Gioia-Di Chiacchio, R. M., Cunha, M. P. V., de Sá, L. R. M., Davies, Y. M., Pereira, C. B. P., Martins, F. H., et al. (2018). Novel hybrid of typical Enteropathogenic Escherichia coli and Shiga-toxin-Producing E. coli (tEPEC/STEC) emerging from pet birds. Front. Microbiol. 9:2975. doi: 10.3389/fmicb.2018.02975
Gomes, T. A. T., Elias, W. P., Scaletsky, I. C. A., Guth, B. E. C., Rodrigues, J. F., Piazza, R. M. F., et al. (2016). Diarrheagenic Escherichia coli. Braz. J. Microbiol. 47, 3–30. doi: 10.1016/j.bjm.2016.10.015
Gomes, T. A. T., Irino, K., Girão, D. M., Girão, V. B. C., Guth, B. E. C., Vaz, T. M. I., et al. (2004). Emerging enteropathogenic Escherichia coli strains? Emerg. Infect. Dis. 10, 1851–1855. doi: 10.3201/EID1010.031093
Hazen, T. H., Donnenberg, M. S., Panchalingam, S., Antonio, M., Hossain, A., Mandomando, I., et al. (2016). Genomic diversity of EPEC associated with clinical presentations of differing severity. Nat. Microbiol. 1:15014. doi: 10.1038/NMICROBIOL.2015.14
Hernandes, R. T., Elias, W. P., Vieira, M. A. M., and Gomes, T. A. T. (2009). An overview of atypical enteropathogenic Escherichia coli. FEMS Microbiol. Lett. 297, 137–149. doi: 10.1111/J.1574-6968.2009.01664.X
Hernandes, R. T., Hazen, T. H., Dos Santos, L. F., Richter, T. K. S., Michalski, J. M., and Rasko, D. A. (2020). Comparative genomic analysis provides insight into the phylogeny and virulence of atypical enteropathogenic Escherichia coli strains from Brazil. PLoS Negl. Trop. Dis. 14, e0008373–e0008319. doi: 10.1371/JOURNAL.PNTD.0008373
Ingle, D. J., Tauschek, M., Edwards, D. J., Hocking, D. M., Pickard, D. J., Azzopardi, K. I., et al. (2016). Evolution of atypical enteropathogenic E. coli by repeated acquisition of LEE pathogenicity island variants. Nat. Microbiol. 1:15010. doi: 10.1038/nmicrobiol.2015.10
Jarquin, C., Morales, O., McCracken, J. P., Lopez, M. R., Lopez, B., Reyes, L., et al. (2022). Burden of Diarrhoeagenic Escherichia coli in Santa Rosa, Guatemala in active health-services surveillance during 2008-2009 and 2014-2015. Trop. Med. Int. Health 27, 408–417. doi: 10.1111/TMI.13735
Joensen, K. G., Scheutz, F., Lund, O., Hasman, H., Kaas, R. S., Nielsen, E. M., et al. (2014). Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli. J. Clin. Microbiol. 52, 1501–1510. doi: 10.1128/JCM.03617-13
Joensen, K. G., Tetzschner, A. M. 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
Johnson, J. R., Clermont, O., Menard, M., Kuskowski, M. A., Picard, B., and Denamur, E. (2006). Experimental mouse lethality of Escherichia coli isolates, in relation to accessory traits, phylogenetic group, and ecological source. J. Infect. Dis. 194, 1141–1150. doi: 10.1086/507305
Johnson, J. R., Kuskowski, M., Denamur, E., Elion, J., and Picard, B. (2000). Clonal origin, virulence factors, and virulence (multiple letters). Infect. Immun. 68, 424–425. doi: 10.1128/IAI.68.1.424-425.2000
Johnson, T. J., Miller, E. A., Flores-Figueroa, C., Munoz-Aguayo, J., Cardona, C., Fransen, K., et al. (2022). Refining the definition of the avian pathogenic Escherichia coli (APEC) pathotype through inclusion of high-risk clonal groups. Poult. Sci. 101:102009. doi: 10.1016/J.PSJ.2022.102009
Johnson, J. R., Murray, A. C., Gajewski, A., Sullivan, M., Snippes, P., Kuskowski, M. A., et al. (2003). Isolation and molecular characterization of nalidixic acid-resistant extraintestinal pathogenic Escherichia coli from retail chicken products. Antimicrob. Agents Chemother. 47, 2161–2168. doi: 10.1128/aac.47.7.2161-2168.2003
Johnson, J. R., and Russo, T. A. (2018). Molecular epidemiology of extraintestinal pathogenic Escherichia coli. EcoSal Plus 8:ESP0004-2017. doi: 10.1128/ecosalplus.ESP-0004-2017
Johnson, J. R., Stapleton, A. E., Russo, T. A., Scheutz, F., Brown, J. J., and Maslow, J. N. (1997). Characteristics and prevalence within serogroup O4 of a J96-like clonal group of uropathogenic Escherichia coli O4:H5 containing the class I and class III alleles of papG. Infect. Immun. 65, 2153–2159. doi: 10.1128/IAI.65.6.2153-2159.1997
Johnson, J. R., and Stell, A. L. (2000). Extended virulence genotypes of Escherichia coli strains from patients with urosepsis in relation to phylogeny and host compromise. J. Infect. Dis. 181, 261–272. doi: 10.1086/315217
Johnson, T. J., Wannemuehler, Y. M., and Nolan, L. K. (2008). Evolution of the iss gene in Escherichia coli. Appl. Environ. Microbiol. 74, 2360–2369. doi: 10.1128/AEM.02634-07
Kaper, J. B., Nataro, J. P., and Mobley, H. L. T. (2004). Pathogenic Escherichia coli. Nat. Rev. Microbiol. 2, 123–140. doi: 10.1038/nrmicro818
Kara, Y., Klzll, M. C., Klllç, Ö., Us, T., and Dinleyici, E. Ç. (2022). Investigation of the use of multiplex PCR in childhood diarrhea with clinical and epidemiological features. J. Trop. Pediatr. 68:fmac090. doi: 10.1093/TROPEJ/FMAC090
Kessler, R., Nisa, S., Hazen, T. H., Horneman, A., Amoroso, A., Rasko, D. A., et al. (2015). Diarrhea, bacteremia and multiorgan dysfunction due to an extraintestinal pathogenic Escherichia coli strain with enteropathogenic E. coli genes. Pathog. Dis. 73:ftv076. doi: 10.1093/femspd/ftv076
Le Bouguenec, C., Archambaud, M., and Labigne, A. (1992). Rapid and specific detection of the pap, afa, and sfa adhesin-encoding operons in uropathogenic Escherichia coli strains by polymerase chain reaction. J. Clin. Microbiol. 30, 1189–1193. doi: 10.1128/JCM.30.5.1189-1193.1992
Letunic, I., and Bork, P. (2021). Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, W293–W296. doi: 10.1093/nar/gkab301
Levine, M. M., Nataro, J. P., Karch, H., Baldini, M. M., Kaper, J. B., Black, R. E., et al. (1985). The diarrheal response of humans to some classic serotypes of enteropathogenic Escherichia coli is dependent on a plasmid encoding an enteroadhesiveness factor. J. Infect. Dis. 152, 550–559. doi: 10.1093/INFDIS/152.3.550
Liberatore, A. M. A., Moreira, F. C., Gomes, T. A. T., Menchaca-Diaz, J. L., and Koh, I. H. J. (2011). Typical and atypical enteropathogenic Escherichia coli bacterial translocation associated with tissue hypoperfusion in rats. Braz. J. Med. Biol. Res. 44, 1018–1024. doi: 10.1590/S0100-879X2011007500105
Liebchen, A., Benz, I., Mellmann, A., Karch, H., Gomes, T. A. T., Yamamoto, D., et al. (2011). Characterization of Escherichia coli strains isolated from patients with diarrhea in São Paulo, Brazil: identification of intermediate virulence factor profiles by multiplex PCR. J. Clin. Microbiol. 49, 2274–2278. doi: 10.1128/JCM.00386-11
Malberg Tetzschner, A. M., Johnson, J. R., Johnston, B. D., Lund, O., and Scheutz, F. (2020). In silico genotyping of Escherichia coli isolates for Extraintestinal virulence genes by use of whole-genome sequencing data. J. Clin. Microbiol. 58:e01269-20. doi: 10.1128/JCM.01269-20
Mandomando, I., Vubil, D., Boisen, N., Quintó, L., Ruiz, J., Sigaúque, B., et al. (2020). Escherichia coli ST131 clones harbouring AggR and AAF/V fimbriae causing bacteremia in Mozambican children: emergence of new variant of fimH27 subclone. PLoS Negl. Trop. Dis. 14:e0008274. doi: 10.1371/journal.pntd.0008274
Mariani-Kurkdjian, P., Lemaître, C., Bidet, P., Perez, D., Boggini, L., Kwon, T., et al. (2014). Haemolytic-uraemic syndrome with bacteraemia caused by a new hybrid Escherichia coli pathotype. N. Microbes N. Infect. 2, 127–131. doi: 10.1002/nmi2.49
Mellata, M., Johnson, J. R., and Curtiss, R. (2018). Escherichia coli isolates from commercial chicken meat and eggs cause sepsis, meningitis and urinary tract infection in rodent models of human infections. Zoonoses Public Health 65, 103–113. doi: 10.1111/ZPH.12376
Miajlovic, H., and Smith, S. G. (2014). Bacterial self-defence: how Escherichia coli evades serum killing. FEMS Microbiol. Lett. 354, 1–9. doi: 10.1111/1574-6968.12419
Munhoz, D. D., Santos, F. F., Mitsunari, T., Schüroff, P. A., Elias, W. P., Carvalho, E., et al. (2021). Hybrid atypical Enteropathogenic and Extraintestinal Escherichia coli (aEPEC/ExPEC) BA1250 strain: a draft genome. Pathogens 10:475. doi: 10.3390/pathogens10040475
Nascimento, J. A. S., Santos, F. F., Santos-Neto, J. F., Trovão, L. O., Valiatti, T. B., Pinaffi, I. C., et al. (2022). Molecular epidemiology and presence of hybrid pathogenic Escherichia coli among isolates from community-acquired urinary tract infection. Microorganisms 10:302. doi: 10.3390/microorganisms10020302
Nascimento, J. A. S., Santos, F. F., Valiatti, T. B., Santos-Neto, J. F., Santos, A. C. M., Cayô, R., et al. (2021). Frequency and diversity of hybrid Escherichia coli strains isolated from urinary tract infections. Microorganisms 9:693. doi: 10.3390/microorganisms9040693
Nielsen, K. L., Godfrey, P. A., Stegger, M., Andersen, P. S., Feldgarden, M., and Frimodt-Møller, N. (2014). Selection of unique Escherichia coli clones by random amplified polymorphic DNA (RAPD): evaluation by whole genome sequencing. J. Microbiol. Methods 103, 101–103. doi: 10.1016/j.mimet.2014.05.018
Ori, E. L., Takagi, E. H., Andrade, T. S., Miguel, B. T., Cergole-Novella, M. C., Guth, B. E. C., et al. (2019). Diarrhoeagenic Escherichia coli and Escherichia albertii in Brazil: Pathotypes and serotypes over a 6-year period of surveillance. Epidemiol. Infect. 147:e10. doi: 10.1017/S0950268818002595
Pereira, T. C., de Barros, P. P., de Oliveira Fugisaki, L. R., Rossoni, R. D., de Ribeiro, F. C., de Menezes, R. T., et al. (2018). Recent advances in the use of galleria mellonella model to study immune responses against human pathogens. J. Fungi 4:128. doi: 10.3390/JOF4040128
Picard, B., Garcia, J. S., Gouriou, S., Duriez, P., Brahimi, N., Bingen, E., et al. (1999). The link between phylogeny and virulence in Escherichia coli extraintestinal infection. Infect. Immun. 67, 546–553. doi: 10.1128/iai.67.2.546-553.1999
Richter, T. K. S., Hazen, T. H., Lam, D., Coles, C. L., Seidman, J. C., You, Y., et al. (2018). Temporal variability of Escherichia coli diversity in the gastrointestinal tracts of Tanzanian children with and without exposure to antibiotics. mSphere 3:e00558-18. doi: 10.1128/MSPHERE.00558-18
Riley, L. W. (2020). Distinguishing Pathovars from Nonpathovars: Escherichia coli. Microbiol Spectr 8:AME-0014-2020. doi: 10.1128/microbiolspec.AME-0014-2020
Riveros, M., García, W., García, C., Durand, D., Mercado, E., Ruiz, J., et al. (2017). Molecular and phenotypic characterization of diarrheagenic Escherichia coli strains isolated from bacteremic children. Am. J. Trop. Med. Hyg. 97, 1329–1336. doi: 10.4269/ajtmh.17-0066
Roer, L., Johannesen, T. B., Hansen, F., Stegger, M., Tchesnokova, V., Sokurenko, E., et al. (2018). CHTyper, a web tool for subtyping of extraintestinal pathogenic Escherichia coli based on the fumC and fimH alleles. J. Clin. Microbiol. 56, 63–81. doi: 10.1128/JCM.00063-18
Sabaté, M., Moreno, E., Pérez, T., Andreu, A., and Prats, G. (2006). Pathogenicity island markers in commensal and uropathogenic Escherichia coli isolates. Clin. Microbiol. Infect. 12, 880–886. doi: 10.1111/j.1469-0691.2006.01461.x
Santos, A. C. M., Fuga, B., Esposito, F., Cardoso, B., Santos, F. F., Valiatti, T. B., et al. (2021). Unveiling the virulent genotype and unusual biochemical behavior of Escherichia coli ST59. Appl. Environ. Microbiol. 87:e0074321. doi: 10.1128/AEM.00743-21
Santos, A. C. M., Santos, F. F., Silva, R. M., and Gomes, T. A. T. (2020a). Diversity of hybrid- and hetero-pathogenic Escherichia coli and their potential implication in more severe diseases. Front. Cell. Infect. Microbiol. 10:339. doi: 10.3389/FCIMB.2020.00339
Santos, A. C. M., Silva, R. M., Valiatti, T. B., Santos, F. F., Santos-Neto, J. F., Cayô, R., et al. (2020b). Virulence potential of a multidrug-resistant Escherichia coli strain belonging to the emerging clonal group ST101-B1 isolated from bloodstream infection. Microorganisms 8:827. doi: 10.3390/microorganisms8060827
Schmidt, H., Beutin, L., and Karch, H. (1995). Molecular analysis of the plasmid-encoded hemolysin of Escherichia coli O157:H7 strain EDL 933. Infect. Immun. 63, 1055–1061. doi: 10.1128/iai.63.3.1055-1061.1995
Šmajs, D., Micenková, L., Šmarda, J., Vrba, M., Ševčíková, A., Vališová, Z., et al. (2010). Bacteriocin synthesis in uropathogenic and commensal Escherichia coli: colicin E1 is a potential virulence factor. BMC Microbiol. 10:288. doi: 10.1186/1471-2180-10-288
Spurbeck, R. R., Dinh, P. C., Walk, S. T., Stapleton, A. E., Hooton, T. M., Nolan, L. K., et al. (2012). Escherichia coli isolates that carry vat, fyuA, chuA, and yfcV efficiently colonize the urinary tract. Infect. Immun. 80, 4115–4122. doi: 10.1128/IAI.00752-12
Toval, F., Köhler, C. D., Vogel, U., Wagenlehner, F., Mellmann, A., Fruth, A., et al. (2014). Characterization of Escherichia coli isolates from hospital inpatients or outpatients with urinary tract infection. J. Clin. Microbiol. 52, 407–418. doi: 10.1128/JCM.02069-13
Trabulsi, L. R., Keller, R., and Tardelli Gomes, T. A. (2002). Typical and atypical enteropathogenic Escherichia coli. Emerg. Infect. Dis. 8, 508–513. doi: 10.3201/eid0805.010385
Valiatti, T. B., Santos, F. F., Santos, A. C. M., Nascimento, J. A. S., Silva, R. M., Carvalho, E., et al. (2020). Genetic and virulence characteristics of a hybrid atypical Enteropathogenic and Uropathogenic Escherichia coli (aEPEC/UPEC) strain. Front. Cell. Infect. Microbiol. 10:492. doi: 10.3389/fcimb.2020.00492
Vieira, M. A. M., Andrade, J. R. C., Trabulsi, L. R., Rosa, A. C. P., Dias, A. M. G., Ramos, S. R. T. S., et al. (2001). Phenotypic and genotypic characteristics of Escherichia coli strains of non-enteropathogenic E. coli (EPEC) serogroups that carry EAE and lack the EPEC adherence factor and Shiga toxin DNA probe sequences. J. Infect. Dis. 183, 762–772. doi: 10.1086/318821
Vieira, M. A. M., Salvador, F. A., Silva, R. M., Irino, K., Vaz, T. M. I., Rockstroh, A. C., et al. (2010). Prevalence and characteristics of the O122 pathogenicity island in typical and atypical enteropathogenic Escherichia coli strains. J. Clin. Microbiol. 48, 1452–1455. doi: 10.1128/JCM.01944-09
Wattam, A. R., Davis, J. J., Assaf, R., Boisvert, S., Brettin, T., Bun, C., et al. (2017). Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center. Nucleic Acids Res. 45, D535–D542. doi: 10.1093/nar/gkw1017
WHO (2017). Diarrhoeal disease. World Health Organization. Available at: https://www.who.int/news-room/fact-sheets/detail/diarrhoeal-disease (Accessed January 9, 2022).
Wijetunge, D. S. S., Gongati, S., Debroy, C., Kim, K. S., Couraud, P. O., Romero, I. A., et al. (2015). Characterizing the pathotype of neonatal meningitis causing Escherichia coli (NMEC) applied microbiology. BMC Microbiol. 15, 1–15. doi: 10.1186/s12866-015-0547-9
Williamson, D. A., Mills, G., Johnson, J. R., Porter, S., and Wiles, S. (2014). In vivo correlates of molecularly inferred virulence among extraintestinal pathogenic Escherichia coli (ExPEC) in the wax moth galleria mellonella model system. Virulence 5, 388–393. doi: 10.4161/viru.27912
Wirth, T., Falush, D., Lan, R., Colles, F., Mensa, P., Wieler, L. H., et al. (2006). Sex and virulence in Escherichia coli: an evolutionary perspective. Mol. Microbiol. 60, 1136–1151. doi: 10.1111/j.1365-2958.2006.05172.x
Yamamoto, S., Terai, A., Yuri, K., Kurazono, H., Takeda, Y., and Yoshida, O. (1995). Detection of urovirulence factors in Escherichia coli by multiplex polymerase chain reaction. FEMS Immunol. Med. Microbiol. 12, 85–90. doi: 10.1016/0928-8244(95)00053-A
Yamani, L. Z., and Elhadi, N. (2022). Virulence characteristics, antibiotic resistance patterns and molecular typing of Enteropathogenic producing Escherichia coli (EPEC) isolates in Eastern Province of Saudi Arabia: 2013-2014. Infect. Drug. Resist. 15, 6763–6772. doi: 10.2147/IDR.S388956
Yousefipour, M., Rezatofighi, S. E., and Ardakani, M. R. (2023). Detection and characterization of hybrid uropathogenic Escherichia coli strains among E. coli isolates causing community-acquired urinary tract infection. J. Med. Microbiol. 72:001660. doi: 10.1099/JMM.0.001660
Zankari, E., Allesøe, R., Joensen, K. G., Cavaco, L. M., Lund, O., and Aarestrup, F. M. (2017). PointFinder: a novel web tool for WGS-based detection of antimicrobial resistance associated with chromosomal point mutations in bacterial pathogens. J. Antimicrob. Chemother. 72, 2764–2768. doi: 10.1093/jac/dkx217
Keywords: ExPEC, UPEC, aEPEC, children diarrhea, Escherichia coli, hybrid pathogenic, hetero-pathogenic
Citation: Santos ACM, Silva RS, Vieira MAM, Niero CV, Ribeiro MSF, Guth BEC, Gomes TAT and Silva RM (2025) The diverse virulence potential of atypical enteropathogenic Escherichia coli isolated from diarrhea: the emergence of a hybrid pathotype? Front. Microbiol. 16:1599350. doi: 10.3389/fmicb.2025.1599350
Edited by:
Moataz Abd El Ghany, The University of Sydney, AustraliaReviewed by:
Domonkos Sváb, HUN-REN Veterinary Medical Research Institute, HungaryMichael Olson, Snow College, United States
Copyright © 2025 Santos, Silva, Vieira, Niero, Ribeiro, Guth, Gomes and Silva. 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) and the copyright owner(s) 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: Rosa M. Silva, cm9zYS51bmlmZXNwQGdtYWlsLmNvbQ==; Ana C. M. Santos, Y2Fyb2xpbmEubWVsbG9AdW5pZmVzcC5icg==
†These authors have contributed equally to this work and share first authorship