Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 02 October 2025

Sec. Antimicrobials, Resistance and Chemotherapy

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1656938

This article is part of the Research TopicAddressing Antimicrobial Resistance: Multimodal Strategies for Low-Resource EnvironmentsView all 5 articles

Comprehensive molecular and phenotypic profiling of uropathogenic Escherichia coli in a Honduran healthcare setting: virulence, resistance and phylogeny

Bryan Ortiz&#x;Bryan Ortiz1Manuel G. Ballesteros-Monrreal&#x;Manuel G. Ballesteros-Monrreal2Celeste Galindo,Celeste Galindo1,3Daniel RiveraDaniel Rivera1Luis RiveraLuis Rivera1Fernando PrezFernando Pérez1Victoria MaldonadoVictoria Maldonado1Elixia ValleElixia Valle1Ana EstradaAna Estrada1Pablo Mendez-Pfeiffer,Pablo Mendez-Pfeiffer2,4Lourdes EnríquezLourdes Enríquez1Dora ValenciaDora Valencia2Gustavo Fontecha
Gustavo Fontecha1*
  • 1Instituto de Investigaciones en Microbiología, Facultad de Ciencias, Universidad Nacional Autónoma de Honduras, Tegucigalpa, Honduras
  • 2Departamento de Ciencias Químico-Biológicas y Agropecuarias, Universidad de Sonora, Hermosillo, Mexico
  • 3Instituto Hondureño de Seguridad Social, Tegucigalpa, Honduras
  • 4Departamento de Bionanotecnología, Centro de Nanociencias y Nanotecnología, Universidad Nacional Autónoma de México, Km 107 Carretera Tijuana-Ensenada, Ensenada, Mexico

Introduction: Urinary tract infections (UTIs) are a major public health concern, further complicated by the rise of multidrug-resistant bacterial strains. Uropathogenic Escherichia coli (UPEC) is the primary causative agent of UTIs, notable for its genetic diversity and its ability to acquire both virulence and antimicrobial resistance determinants.

Methods: This study conducted a comprehensive phenotypic and molecular characterization of 126 UPEC strains isolated from a tertiary care hospital in Tegucigalpa, Honduras. Seventeen virulence genes were screened, antimicrobial susceptibility to 17 antibiotics was assessed, phylogenetic grouping was performed, and potential clonal relationships were analyzed using ERIC-PCR.

Results: Strains isolated from male patients exhibited significantly higher virulence gene counts (mean: 10.48 vs. 8.06; p = 0.0029), resistance indices (RI = 0.46 vs. 0.27; p < 0.0001), and multidrug resistance rates (88% vs. 63%; p = 0.009) compared to those from female patients. Extended-spectrum β-lactamase (ESBL) production was observed in 42% of isolates, with a higher prevalence in males (59%; p = 0.049). Phylogroup B2 was the most frequent (29%) and was significantly associated with virulence genes papG-II, hlyA, cnf-1, fyuA, and iucD. Despite high genetic heterogeneity observed through ERIC-PCR, clonal clusters sharing similar phylogroups, virulence profiles, and resistance phenotypes were identified. A weak but significant correlation was found between virulence and resistance indices (r = 0.1796; p = 0.0442).

Discussion: This study provides the first in-depth molecular and phenotypic characterization of UPEC in Honduras. The detection of highly virulent and multidrug-resistant strains underscores the need to reinforce local molecular surveillance and to revise empirical treatment guidelines based on local epidemiological data.

1 Introduction

Urinary tract infections (UTIs) are a major global public health concern, affecting over 150 million individuals each year (Öztürk and Murt, 2020). These infections are associated with considerable morbidity and encompass a broad spectrum of clinical manifestations, ranging from asymptomatic bacteriuria and cystitis to severe pyelonephritis, which may progress to septic shock and potentially fatal multiorgan failure (Agarwal et al., 2012; McLellan and Hunstad, 2016). In 2019 alone, an estimated 400 million UTI cases were reported worldwide, resulting in over 230,000 deaths and 520,200 disability-adjusted life years (DALYs)—a nearly 70% increase compared to 1990 (Yang et al., 2022).

Escherichia coli is a Gram-negative bacterium that typically resides as a commensal in the gastrointestinal tracts of humans and animals (Lagerstrom and Hadly, 2023; Lagerstrom and Hadly, 1948). However, its high genomic plasticity enables it to acquire exogenous genetic elements, such as plasmids, transposons, insertion sequences, and integrons (Yousuf et al., 2020; Tenaillon et al., 2010). This capacity promotes the emergence of strains with enhanced virulence and antimicrobial resistance, posing a significant threat to public health (Lagerstrom and Hadly, 1948; Whelan et al., 2023). Based on infection site, clinical presentation, and genotypic traits, E. coli strains are broadly classified into two categories: intestinal pathotypes, which cause diarrheal disease, and extraintestinal pathotypes (Regua-Mangia et al., 2010). The latter include uropathogenic E. coli (UPEC), neonatal meningitis-associated E. coli (NMEC), sepsis-associated E. coli (SEPEC), avian pathogenic E. coli (APEC), and mammary pathogenic E. coli (MPEC) (Sora et al., 2021).

E. coli strains are currently categorized into eight major phylogenetic groups (A, B1, B2, C, D, E, F, and G), which reflect their evolutionary relationships and pathogenic potential. In addition, five cryptic clades (I–V) have been described, corresponding to phylogenetically related but distinct Escherichia species (Clermont et al., 2000; Clermont et al., 2013; Clermont et al., 2019). Most extraintestinal pathogenic strains belong to phylogroups B2 and D (Halaji et al., 2022), whereas commensal strains are more commonly associated with phylogroups A and B1 (da Silva et al., 2017).

UPEC is the predominant cause of UTIs in both community and hospital settings, accounting for approximately 80% of community-acquired and over 50% of healthcare-associated infections (Kang et al., 2018; Xia et al., 2017; Kucheria et al., 2005). It is responsible for more than 75% of uncomplicated UTIs and about 65% of complicated cases (Whelan et al., 2023; Medina and Castillo-Pino, 2019). This high prevalence is closely linked to the virulence potential of phylogroups B2 and D, which harbor strains capable of expressing a wide range of virulence factors (Halaji et al., 2022; Royer et al., 2023). These include adhesins, toxins, iron acquisition systems, and immune evasion mechanisms, all of which facilitate adherence, invasion, and persistence within the urinary tract epithelium, contributing to recurrent and difficult-to-treat infections (Alam Parvez and Rahman, 2018; Behzadi and Carevic, 2019; Ballesteros-Monrreal et al., 2021).

In addition to their pathogenicity, UPEC strains often exhibit multidrug resistance (MDR), phenotypes, representing a growing public health challenge (Akhand and Karicheri, 2022). MDR infections are projected to cause an increasing number of deaths globally by 2050 (de Kraker et al., 2016). Furthermore, as with virulence traits, antimicrobial resistance profiles in UPEC vary considerably by geographic region, complicating empirical treatment decisions (Yang et al., 2022; Ballesteros-Monrreal et al., 2021; Ramakrishnan et al., 2022). This highlights the urgent need for local and regional epidemiological studies to better understand UPEC’s resistance and virulence patterns, thereby informing and updating clinical management guidelines.

To date, no epidemiological surveillance study in Honduras has addressed the molecular characterization of UPEC in clinical settings. Therefore, this study aimed to characterize the antimicrobial resistance, virulence gene distribution, and phylogenetic profiles of UPEC isolates recovered from a major tertiary hospital in Tegucigalpa, Honduras. These findings are intended to support clinicians and public health stakeholders by providing data for more effective management of UTIs in the region.

2 Materials and methods

2.1 Clinical isolates

Urine cultures phenotypically identified as Escherichia coli with colony counts ≥100,000 CFU/mL, processed between July and August 2022 at a tertiary-level hospital of the Honduran Social Security Institute (IHSS) in Tegucigalpa, were donated to this study. The isolates were provided along with anonymized sociodemographic data (sex, age, and hospitalization status). No additional clinical information, such as treatment history, underlying diseases, or urinary cathether use, was available. Only molecular characterization of the isolates was performed. No direct contact with patients occurred at any stage of the study.

2.2 Phenotypic identification and antimicrobial susceptibility testing

Bacterial identification was performed using the BD Phoenix™ Identification System for Gram-negative bacilli, following the manufacturer’s instructions. Briefly, one to two colonies from a 24-h culture were suspended in Phoenix ID broth to a 0.5 McFarland turbidity standard. From this standardized suspension, 25 μL were transferred to Phoenix AST broth, which was supplemented with a drop of Phoenix AST indicator (containing methylene blue and resazurin) to detect bacterial growth. The mixture was then loaded into BD Phoenix™ NMIC-500 panels (catalog number BD 449023), which were sealed, registered, and processed using the BD Phoenix™ M50 system. Results were analyzed using Epicenter data management software (version 6.61A, BD Diagnostic Systems) after incubation for 6 to 16 h.

Antibiotics were grouped into 12 categories, according to the classification proposed by Magiorakos et al. (2012) (1) Aminoglycosides: Amikacin (AMK), Gentamicin (GM) (2) Sulfonamides: Trimethoprim/Sulfamethoxazole (TSX); (3) Nitrofurans: Nitrofurantoin (MAC); (4) Fluoroquinolones: Ciprofloxacin (CIP), Levofloxacin (LVX); (5) Penicillins with β-lactamase inhibitors: Amoxicillin/Clavulanic acid (AMC), Piperacillin/Tazobactam (TZP); (6) Penicillins: Ampicillin (AMP); (7) First-generation cephalosporins: Cefuroxime (CFX); (8) Second-generation cephalosporins: Cefazolin (CZ); (9) Third-generation cephalosporins: Ceftriaxone (CRO); (10) Fourth-generation cephalosporins: Cefepime (FEP); (11) Cephamycin: Cefoxitin (FOX); and (12) Carbapenems: Ertapenem (ETP), Imipenem (IMP), Meropenem (MEM). Based on the resistance profile of each isolate, strains were classified as follows: Non-multidrug-resistant (NMDR): Non-susceptible to two or fewer categories; Multidrug-resistant (MDR): Non-susceptible to three or more categories; extensively drug-resistant (XDR): Susceptible to only one of the tested categories; Pandrug-resistant (PDR): Non-susceptible to all tested categories (Magiorakos et al., 2012).

The prevalence of resistance for each antibiotic was graphically represented as a percentage. Additionally, a resistance index (RI) was calculated for each isolate using the formula previously reported by Karim et al. (2023),

RI=ab

Where “a” represents the number of antibiotics to which a given isolate was resistant, and “b” corresponds to the total number of antibiotics tested in this study. The resulting RI values were subsequently used for comparative analyses based on patient sex, age group, phylogroup classification, and other variables.

2.3 DNA extraction and molecular identification

Phenotypically confirmed E. coli isolates were inoculated in Luria-Bertani (LB) broth (Liofilchem®) and incubated at 37 °C for 24 h. Genomic DNA was extracted using the Wizard® Genomic DNA Purification Kit (Promega, Madison, WI, USA), following the manufacturer’s protocol. DNA extracts were eluted in 100 μL of elution buffer and stored at −20 °C until further use.

Molecular identification of E. coli was performed by PCR amplification of the ybbW gene, which encodes an allantoin receptor reported as highly specific to this species (Walker et al., 2017). PCR reactions were conducted in a final volume of 25 μL containing 12.5 μL of 2 × PCR Master Mix (Promega Corp. Madison, WI, USA), 1 μL of each primer (10 μM), 8 μL of nuclease-free water, and 1.5 μL of DNA (40 ng/μL). Thermocycling conditions were as follows: initial denaturation at 95 °C for 4 min; 37 cycles of denaturation at 95 °C for 70 s, annealing at 60 °C for 70 s, and extension at 72 °C for 70 s; followed by a final extension step at 72 °C for 10 min. All PCR products were resolved on 2% agarose gels stained with ethidium bromide and visualized using the Gel Doc™ EZ Imaging System. Primer sequences used for amplification are listed in Supplementary Table 1.

2.4 Virulence gene profiling

Seventeen virulence genes associated with UPEC pathogenesis in the urinary tract were analyzed. These genes were functionally grouped into five categories: (1) Adhesion: fimH, papG-II, papC, sfaS/focC; (2) Motility: fliCD; (3) Toxin production: hlyA, cnf-1, vat, sat; (4) Immune evasion: agn43, traT, kpsM; and (5) Iron acquisition: feoB, fyuA, iucD, iha, iutA. Gene detection was performed via multiplex PCR following the protocol described by Ballesteros-Monrreal et al. (2021). PCR reactions were carried out in a final volume of 25 μL, containing 12.5 μL of 2 × PCR Master Mix (Promega Corp. Madison, WI, USA), 1 μL of each primer (10 μM), and 2 μL of genomic DNA. Positive controls included DNA from E. coli strains CFT073, ATCC 25922, and UPEC 39, as described in prior publications (Ballesteros-Monrreal et al., 2021). Primer sequences, amplicon sizes, and annealing temperatures (Tm, °C) are provided in Supplementary Table 1.

2.5 Phylogenetic group determination

Phylogenetic group assignment was performed using Clermont et al. (2000, 2013, 2019) methods, based on PCR detection of the arpA, chuA, yjaA, trpA, ybgD, and cfaB genes, as well as the TspE4. C2 DNA fragment. PCR reactions were conducted following previously described protocols (Mendoza et al., 2023).

2.6 Eric PCR and clonality estimation

Clonality among all isolates was assessed through the analysis of consensus sequences of repetitive intergenic regions (ERIC), a technique widely validated and recommended for molecular typing of E. coli (Alsultan and Elhadi, 2022; López-Ramírez et al., 2018; Dishan et al., 2023; Ardakani and Ranjbar, 2016; Mahmoud et al., 2020). The methodology followed that described by Movahedi et al. (2021). PCR reactions were carried out in a final volume of 25 μL containing 12.5 μL of 2 × PCR Master Mix (Promega Corp., Madison, WI, USA), 8.5 μL of nuclease-free water, 2 μL of genomic DNA, and 1 μL of each ERIC primer. Thermocycling conditions were as follows: initial denaturation at 94 °C for 5 min; 40 cycles of denaturation at 94 °C for 1 min, annealing at 53 °C for 40 s, and extension at 72 °C for 1 min; followed by a final extension at 72 °C for 5 min. Band sizes, in relation to the molecular weight marker, as well as the similarity matrix for dendrogram construction, were generated using the GelJ v2.0 software.1 Cluster analysis was performed using the unweighted pair group method with arithmetic mean (UPGMA) and the DICE similarity coefficient. To avoid overestimation of clonality, isolates were classified as “suspected clones” only when they shared ≥90% similarity in ERIC-PCR patterns and exhibited identical or highly similar resistance, virulence, and phylogroup profiles. This integrative approach ensured a more stringent and reliable interpretation of genetic relatedness. Clonal distribution was further analyzed in relation to hospital wards and patient origin (hospitalized vs. non-hospitalized).

2.7 Statistical analyses and graphical representations

A comprehensive statistical analysis approach was employed to ensure the reliability and interpretability of results. Descriptive statistics were initially calculated, including means, medians, and interquartile ranges, according to the distribution of the data. Assumptions of normality and homogeneity of variances were evaluated using the Shapiro–Wilk and Bartlett tests, respectively. Due to violations of these assumptions in some groups, the non-parametric Kruskal–Wallis test was used for overall group comparisons. When significant differences were observed, post hoc analyses were conducted using Dunn’s test with Bonferroni correction to determine which specific groups differed. For categorical variables such as sex, phylogenetic group, or gene presence/absence, Fisher’s exact test was used. Additionally, odds ratios (ORs) and relative risks (RRs) with 95% confidence intervals were calculated to assess the strength of associations between categorical variables. To explore potential associations between virulence factors and antimicrobial resistance, Pearson correlation analysis was performed. A p < 0.05 was considered statistically significant in all analyses.

All statistical analyses were performed using GraphPad Prism version 7.0 (GraphPad Software, San Diego, CA, USA)2 and JASP v0.19.3.3 Graphical representations were generated using GraphPad Prism and RStudio with R version 4.5.0 (2025-04-11). Data visualization in R was conducted using the ggplot2, pheatmap, and magick packages (R Core Team, 2025). ERIC-PCR dendrograms were visualized and edited using the interactive iTOL environment, version 7.2 (Letunic and Bork, 2021).

3 Results

3.1 Patients and clinical isolates

Between July and August 2022, a total of 126 clinical isolates were identified as Escherichia coli through both biochemical and molecular methods. Of these, 99 isolates (78.6%) were obtained from female patients, whose ages ranged from 0 to 88 years, while the remaining 27 isolates (21.4%) were from male patients aged between 0 and 85 years.

Analysis of age group distribution revealed that the 13–50 years category was the most prevalent among both sexes, accounting for 41% of female and 59% of male patients. This was followed by patients older than 50 years, comprising 38% of females and 22% of males. However, no statistically significant association was found between age groups and patient sex (p > 0.05).

3.2 Virulence genes

3.2.1 Virulence genes by sex and age group

Following molecular confirmation of E. coli isolates, the presence of 17 virulence-associated genes was assessed. Isolates from female patients carried an average of 8.06 ± 3.26 virulence genes (range: 3–15), while those from male patients showed a significantly higher average of 10.48 ± 2.69 genes (range: 7–15). This difference was statistically significant (p = 0.0029), suggesting a higher virulence gene burden in isolates from male patients.

When stratified by age group, isolates from females under the age of five harbored the highest average number of virulence genes (9.43 ± 3.85), although this trend did not reach statistical significance (p > 0.05). Among male patients, the 6–12 years age group showed the highest virulence gene load, with an average of 13.3 ± 3.05 genes, a statistically significant difference when compared to the 13–50 and >50 years age groups (p < 0.05).

Regarding the prevalence of individual genes, the most frequently detected across both sexes were the type 1 fimbrial adhesin gene fimH (96%) and the iron acquisition gene feoB (100%) (Figure 1A). Notably, a higher overall prevalence of virulence genes was observed in isolates from male patients. However, statistically significant differences (p < 0.05) were observed only for the following genes: papG-II, papC, hlyA, cnf-1, kpsM, fyuA, iucD, and iha (Figure 1B), which may indicate a potential sex-specific association of certain virulence factors with infections in male patients.

Figure 1
Two bar charts display the prevalence of various genes. Chart A shows overall prevalence percentages for genes in a sample of 126, with fimH and feoB at 98 percent and 100 percent, respectively. Chart B compares the gene prevalence between women (pink, n=99) and men (blue, n=27). The fimH gene shows 100 percent prevalence in men and 97 percent in women. Significant differences are marked with asterisks; for example, papG_II is significantly more prevalent in men. Charts highlight gender-based differences in gene prevalence.

Figure 1. Prevalence of virulence-associated genes (VAGs) in clinical UPEC isolates. (A) Overall prevalence of VAGs. (B) Prevalence of VAGs stratified by patient sex. fimH, Fimbrial adhesin of type 1 pilus; sfaD/focC, S fimbriae minor subunit / F1C fimbriae chaperone; papG-II, Type P pilus adhesin allele 2; papC, Type P pilus chaperone; fliCD, Flagellin subunit / flagellar cap; sat, Autotransporter secreted toxin; hlyA, α-hemolysin; vat, Vacuolating autotransporter toxin; cnf-1, Cytotoxic necrotizing factor 1; kpsM, Capsular antigen variant; agn43, Antigen 43; traT, Serum resistance protein; fyuA, Yersiniabactin receptor; iucD, Aerobactin biosynthesis gene; iha, Bifunctional enterobactin receptor/adhesin protein; iutA, Ferric aerobactin receptor; feoB, Ferrous iron transport protein B.

3.2.2 Virulence genes and clinical setting

Potential differences in the distribution of virulence genes according to the clinical origin of the isolates—hospital-acquired vs. community-acquired—were investigated. Although no statistically significant differences were found in the overall prevalence of virulence genes between these patient groups (p > 0.05) (Supplementary Figure 1), noteworthy pathogenic trends emerged. Isolates from outpatients more frequently carried genes associated with iron acquisition, serum resistance, and renal adherence, whereas hospital-acquired isolates exhibited a higher frequency of genes linked to motility and toxigenicity.

The distribution of virulence genes was further analyzed by hospital ward of origin. Among isolates from female patients, the General Medicine, Gynecology, Outpatient Clinics, Internal Medicine, and Emergency wards showed a higher prevalence of each virulence gene (Figure 2A). However, these differences were not statistically significant (p > 0.05). Similarly, in male patients, a higher prevalence of virulence genes was noted in isolates from General Medicine, Urology, Pediatrics, Emergency, and Outpatient wards (Figure 2B), yet no statistically significant differences were detected (p > 0.05).

Figure 2
Heatmaps and box plots of virulence indices are shown. Image A displays a heatmap of virulence genes across different clinical areas, categorized by adherence, immune evasion, iron uptake, motility, and toxicity. Image B provides a similar heatmap for a smaller set of clinical areas. Image C features box plots of the index virulence for various hospital zones, indicating the spread and median values. Image D shows smaller box plots for selected hospital zones, highlighting lower indices compared to Image C.

Figure 2. Distribution of virulence-associated genes (VAGs) by hospital ward. (A) Female patients; (B) male patients. Virulence indices by hospital ward in (C) female patients; (D) male patients. GYN, Gynecology; INT, Internal medicine; EXT CON, External consultation; EMERG, Emergency; PED, Pediatrics; MED, Medicine; URO, Urology. In the box plots, each box represents the interquartile range (IQR), i.e., the interval between the 25th percentile (Q1) and the 75th percentile (Q3), encompassing the central 50% of the data. The line within the box indicates the median, while the whiskers extend to the minimum and maximum values within 1.5 times the IQR.

To provide a more integrated assessment of the virulence gene load per strain, individual virulence indices were calculated by dividing the number of detected genes by the total number evaluated (n = 17). These indices were compared across hospital wards using the non-parametric Kruskal–Wallis test, followed by Dunn’s post hoc test with Bonferroni correction. Boxplots were generated to visualize the distribution of virulence indices by ward for female (Figure 2C) and male (Figure 2D) patients. Each plot displays the median, interquartile range (IQR), and outliers, allowing for the evaluation of dispersion patterns and relative virulence burden across clinical units. To ensure robust statistical analysis, wards with fewer than five isolates were grouped under the category “Others.”

Among female patients, although no statistically significant differences were observed between wards (p > 0.05), notable trends were identified. Isolates from the Medicine, Pediatrics, and Outpatient wards showed proportionally higher virulence indices compared to those from Gynecology and Internal Medicine. However, no pairwise comparison reached statistical significance (p > 0.99) (Figure 2C). Similarly, in male patients, no statistically significant differences were observed (p > 0.05). Nonetheless, the highest virulence indices were found in isolates from wards classified under “Others” and from the Emergency ward, while the lowest indices were recorded in Urology and Internal Medicine (Figure 2D). Although these differences were not statistically significant, the observed trends may reflect a differential distribution of more virulent strains across specific clinical areas.

3.2.3 Virulence profiles

The analysis of virulence profiles revealed high diversity, with a total of 92 distinct profiles distributed between male and female patients (Supplementary material 1A). Among female patients, 70 unique profiles were identified; of these, 14 were shared among 44 isolates (Supplementary material 1B). In male patients, 22 distinct profiles were detected, two of which were shared among seven isolates (Supplementary material 1C).

3.2.4 Co-ocurrence of virulence genes

Several virulence genes were found to co-occur within the same clinical isolates. Statistical analysis was conducted to assess the significance of this co-occurrence (Supplementary material 1D). Significant co-occurrence was observed for multiple genes, particularly agn43, cnf-1, fyuA, hlyA, iucD, kpsM, papC, papG-II, sat, sfaD, traT, and vat, suggesting potential genetic linkage, likely mediated by mobile genetic elements.

3.3 Antimicrobial resistance

3.3.1 Overall resistance by sex

Resistance to 17 antibiotics, including those recommended for UTI treatment in Honduras, was evaluated. Overall, the resistance index (RI) in isolates from female patients (n = 99) was 0.27 ± 0.20 (median: 0.2353), while isolates from male patients (n = 27) showed a significantly higher RI, with a mean of 0.46 ± 0.22 and a median of 0.5294. Since data from the female group did not follow a normal distribution (p < 0.0001; Shapiro–Wilk test), the nonparametric Mann–Whitney test was applied, revealing a statistically significant difference between groups (p < 0.0001). This result suggests a higher burden of antimicrobial resistance in isolates from male patients compared to those from female patients.

3.3.2 Resistance by age group and multidrug resistance

Comparison of resistance indices by age group within each sex showed no statistically significant differences. However, males over 50 years of age exhibited significantly higher RIs than females under 5 years (p = 0.0108), females aged 13 to 50 years (p = 0.0077), and females over 50 years (p = 0.0410). These findings suggest a trend toward increased antimicrobial resistance in older male patients, although caution is warranted due to limited sample sizes in some subgroups, which may reduce statistical power and underestimate actual differences, particularly among males under 50.

Based on antibiotic resistance profiles, 67% of the 126 isolates were classified as MDR (multidrug-resistant), 2% as XDR (extensively drug-resistant), and 31% as NMDR (non-multidrug-resistant). The highest relative proportion of MDR isolates was observed in male patients (88%) compared to 63% in female patients. This difference was statistically significant (p = 0.009).

3.3.3 Prevalence of antibiotic resistance

The frequency of resistance to key antimicrobial categories included in the study was evaluated. Among the 126 clinical isolates analyzed, the highest prevalence of resistance was observed against penicillins (AMP), sulfonamides (TSX), fluoroquinolones (LVX, CIP), and second- to fourth-generation cephalosporins (CFX, CRO, FEP) (Figure 3A). Resistance levels were consistently higher in isolates from male patients compared to those from female patients across all tested antimicrobials. However, statistically significant differences were observed only for aminoglycosides (GM, p = 0.0246), penicillins (AMP, p = 0.0151), β-lactam/β-lactamase inhibitor combinations (AMC, p = 0.0044), first- to fourth-generation cephalosporins (CZ, p = 0.0044; CRO, p = 0.0078; FEP, p = 0.0018; CFX, p = 0.0177), carbapenems (IMP, p = 0.0446; MEM, p = 0.0307), nitrofurans (MAC, p = 0.0327), and fluoroquinolones (CIP, p = 0.0085; LVX, p = 0.0084) (Figure 3B).

Figure 3
Two bar charts depict antibiotic resistance prevalence. Chart A shows overall prevalence across various antibiotics, with AMP having the highest at eighty percent. Chart B compares prevalence between women and men for each antibiotic. Notably, men have higher resistance to GM and CIP, while women have a higher prevalence for GM and AMP. Both charts indicate variation in resistance percentages across different antibiotics.

Figure 3. Prevalence of resistance to the 17 antibiotics evaluated. (A) Overall prevalence; (B) Prevalence according to the sex of the patients. AMK, Amikacin; GM, Gentamicin; TSX, Trimethoprim-Sulfamethoxazole; AMC, Amoxicillin-Clavulanic Acid; AMP, Ampicillin; CZ, Cefazolin; FEP, Cefepime; FOX, Cefoxitin; CRO, Ceftriaxone; CFX, Cefuroxime; ETP, Ertapenem; IMP, Imipenem; MEM, Meropenem; TZP, Piperacillin-Tazobactam; MAC, Nitrofurantoin; CIP, Ciprofloxacin; LVX, Levofloxacin.

The observed differences were supported by odds ratios (OR) and relative risk (RR) values indicating a higher probability of resistance in male patients, with confidence intervals that in all cases excluded the null value (Supplementary material 2A).

When comparing resistance by age and sex, adult males over 50 years of age exhibited the highest resistance burden, with rates exceeding 50% for several antibiotics, including gentamicin (GM), ampicillin (AMP), amoxicillin-clavulanic acid (AMC), cephalosporins, and fluoroquinolones. Elevated resistance levels were also observed in males aged 13–50 years, although slightly lower (Supplementary material 2A). In contrast, isolates from pediatric males showed low prevalence of resistance, with no cases reported for carbapenems and only one for fluoroquinolones. Among adult and pediatric females, resistance levels were moderate overall; however, in girls under 5 years and those aged 6–12, very high prevalence was observed, particularly for trimethoprim/sulfamethoxazole (up to 100%), amoxicillin-clavulanic acid (40%), and ampicillin (80%) (Supplementary material 2B).

3.3.4 ESBL, carbapenemase production, and multidrug resistance

According to the Phoenix™ system results, 53 (42.06%) of the 126 isolates were extended-spectrum β-lactamase (ESBL) producers, while 10 (7.93%) produced carbapenemases. Sex-based analysis revealed that 37% of the isolates from female patients and 59% of those from male patients were ESBL producers, a statistically significant difference (p = 0.049). In contrast, carbapenemase production was similar in both groups, with 8% in females and 7% in males, showing no significant difference (p > 0.05). Notably, a considerable proportion of ESBL-producing strains were observed in female pediatric patients, but not in male infants (Supplementary material 2B).

3.3.5 Co-ocurrence of antibiotic resistance

As multiple isolates showed simultaneous resistance to various antimicrobial classes, Pearson correlation analysis was conducted to identify significant statistical associations among resistance profiles, including phenotypic profiles of ESBL and carbapenemase (CAR) production.

This analysis revealed positive associations between multiple pairs of antibiotics, particularly a strong correlation between ceftriaxone (CRO) and ciprofloxacin (CIP) (r = 0.433, p < 0.001), suggesting concurrent resistance to β-lactams and fluoroquinolones. High correlations were also found between CRO and other cephalosporins such as cefepime (FEP) (r = 0.87, p < 0.001) and cefuroxime (CFX) (r = 0.873, p < 0.001), as well as with ertapenem (ETP) (r = 0.28, p < 0.001) and meropenem (MEM) (r = 0.209, p = 0.009). CIP also showed significant correlations with GM, AMC, AMP, FEP, and nitrofurantoin (MAC), among others (r values between 0.20 and 0.43; p < 0.01) (Supplementary material 2C).

The Pearson correlation analysis also demonstrated a highly significant association between the ESBL phenotype and resistance to cephalosporins such as CRO (r = 0.886), FEP (r = 0.821), and CFX (r = 0.794), all with p < 0.001. Positive correlations were also found between ESBL and resistance to fluoroquinolones such as CIP (r = 0.355, p < 0.001), as well as piperacillin-tazobactam (TZP) (r = 0.282) and MAC (r = 0.187). Regarding the CAR phenotype, strong correlations were identified with ETP (r = 0.891), imipenem (IMP) (r = 0.198), and MEM (r = 0.449), all statistically significant (p < 0.05).

3.3.6 Virulence vs. resistance

Additionally, when comparing the virulence and resistance indices of the 126 isolates, a weak but statistically significant positive correlation was observed between the two variables (r = 0.1796, p = 0.0442). This finding suggests that some isolates simultaneously carry a high virulence gene load and an extensive antimicrobial resistance profile.

3.4 Phylogenetic groups

3.4.1 Phylogroups in the general population

The phylogenetic groups were identified, revealing a diverse distribution among the 126 strains evaluated. The most prevalent phylogroup was B2, present in 29% of the isolates, followed by groups A and D, both with a prevalence of 21%. Lower proportions were observed for phylogroups B1 (10%), Clade I (6%), F (4%), G (2%), and C (1%). Additionally, 5% of the isolates could not be assigned to any known phylogroup and were classified as “Unknown Phylogroup” (Figure 4A).

Figure 4
Bar charts compare prevalence percentages of phylogenetic groups. Chart A shows overall prevalence, with group B2 having the highest at 29 percent. Chart B distinguishes gender, showing men have a higher prevalence in group B2 at 59 percent, compared to women's 21 percent.

Figure 4. Prevalence of E. coli phylogenetic groups. (A) Prevalence of each phylogenetic group among the total number of clinical isolates analyzed; (B) Prevalence according to the sex of the patient. “?”: Unknown phylogroup.

3.4.2 Phylogroups by sex

When analyzing distribution by sex, a greater phylogenetic diversity was observed in isolates from female patients compared to those from males. In female samples, the most represented phylogroups were A and D, each with a prevalence of 24%, followed by B2 at 21%. In contrast, male isolates showed a marked predominance of phylogroup B2, accounting for 59%, a statistically significant difference compared to females. Additionally, phylogroup C was exclusively identified in isolates from male patients (Figure 4B). The higher prevalence of phylogroup B2 among isolates from male patients aligns with the increased virulence gene burden observed in this group.

3.4.3 Phylogroups by hospital ward

After characterizing the overall prevalence of phylogroups among the clinical isolates, their distribution across different hospital wards was explored in order to identify potential associations with the clinical environment of origin. To this end, isolates were stratified according to the patient’s sex and the hospital ward of origin, and the relative proportions of each phylogroup within each group were calculated (Figures 5A,B). Although apparent differences in distribution were observed, no statistically significant associations were found between phylogenetic groups and hospital areas (p > 0.05).

Figure 5
Two heatmaps display the prevalence of phylogenetic groups across various hospital zones. Image A shows shades of purple, indicating prevalence levels, with different hospital zones on the vertical axis and phylogenetic groups on the horizontal axis. Image B, similar in layout, uses shades of blue and gray to depict prevalence variations. Each section includes numerical values corresponding to the prevalence level.

Figure 5. Distribution of phylogenetic groups by hospital ward. (A) Distribution of phylogroups by hospital ward in female patients; (B) Distribution of phylogroups by hospital ward in male patients. “?”: Unknown phylogroup.

3.4.4 Virulence and phylogenetic groups

Differences in the total number of virulence genes among the various phylogenetic groups were assessed using the virulence index (VI). Due to the small sample sizes of phylogroups C (n = 1) and G (n = 3), these were grouped under the category “Others” to avoid statistical bias.

The analysis revealed statistically significant differences among phylogroups (p < 0.05). In particular, phylogroup B2 showed a significantly higher VI compared to phylogroups A (p = 0.0495) and B1 (p = 0.0406), reinforcing its association with a higher genetic virulence load. While no significant differences were observed with other groups, it is noteworthy that phylogroup F displayed an even higher mean than B2. However, this difference was not statistically significant, likely due to the small sample size of group F (n = 5), which limits the statistical power of the analysis. Similarly, other groups with low representation, such as Clade I or those classified as “Others,” showed variable medians without reaching statistical significance (Supplementary Figure 2). These findings suggest that, although phylogroup B2 harbors the highest proportion of strains with multiple virulence genes, the pathogenic potential of other phylogroups, such as F, may be underestimated and should be further assessed in studies with larger sample sizes. When examining the distribution of virulence genes by phylogenetic group, similar patterns were observed among several phylogroups, particularly B2, D, and those classified as unknown (Figure 6B). This trend was consistent in isolates from female patients but was not observed in those from male patients (Supplementary Figures 3, 4, respectively).

Figure 6
Panel A shows seven radar charts, each representing different genetic profiles labeled as PG with varying sample sizes. Each chart features different genes plotted in a circular format. Panel B displays four radar charts overlaid in different colors, showing comparative genetic expressions across PG groups. Each gene is placed at distinct radial axes with markings for reference.

Figure 6. Virulence and phylogenetic groups. (A) Distribution of virulence genes by phylogenetic group; (B) Phylogroups showing similarities in the distribution of virulence genes. The presence of asterisks indicates a statistically significant association (p < 0.05).

Additionally, an analysis was conducted to determine whether specific virulence genes were more prevalent in a given phylogenetic group. Overall, phylogroup B2 displayed the highest number of statistically significant associations, with a high prevalence of virulence genes. Notably, positive associations were found for the genes papG-II, papC, sat, hlyA, vat, cnf-1, kpsM, fyuA, iucD, and iha, all of which were detected more frequently in isolates belonging to this phylogroup (Figure 6A). Odds ratio (OR) and relative risk (RR) values (Supplementary material 3) confirmed these associations, indicating that B2 strains not only exhibited a higher prevalence of these genes but also had a significantly greater likelihood of harboring them compared to other phylogroups. Some genes, such as sat (OR = 4.52; RR = 1.68), cnf-1 (OR = 4.51; RR = 3.43), vat (OR = 4.3; RR = 3.75), hlyA (OR = 3.96; RR = 1.57), and kpsM (OR = 3.25; RR = 1.68), showed the highest values, indicating a strong association with key properties such as cytotoxicity and immune evasion.

In contrast, phylogroups A and B1 showed a lower prevalence of most of the aforementioned virulence genes, with statistically significant differences. These findings support the association of phylogroup B2 with virulence factors related to cytotoxicity and adhesion at both bladder and renal levels, and are consistent with previous evidence describing phylogroups A and B1 as lower-virulence or commensal groups.

3.4.5 Resistance and phylogenetic group

Differences in the resistance index (RI) among the various phylogenetic groups were evaluated using the non-parametric Kruskal–Wallis test, as the assumption of normality was not met. Although visual differences were observed in the range of RI values for certain phylogroups—particularly for group A (in males), B2 (in all isolates), and the unknown group (in females)—the analyses conducted on the overall population, as well as those stratified by sex, did not show statistically significant differences (p = 0.176) (Figures 7AC). These findings suggest that, based on the available data, no systematic difference in resistance index can be confirmed among the phylogenetic groups evaluated.

Figure 7
Bar graphs labeled A, B, and C comparing the resistance index (RI) across different phylogenetic groups labeled with symbols: ?, B2, CI, A, D, G, B1, and F. Each graph shows varied heights of bars and distribution of data points, representing the RI for each group, with error bars indicating variability. Graph A has all groups, B lacks group G, and C shows groups A, D, B1, B2, and F.

Figure 7. Distribution of the total number of antibiotics to which Escherichia coli isolates are resistant, according to their phylogenetic group. (A) All isolates (n = 126); (B) Isolates from female patients (n = 99); (C) Isolates from male patients (n = 27). In each graph, bars indicate the median and interquartile ranges (IQR).

When evaluating the association between resistance to specific antibiotics and phylogenetic groups, statistically significant differences (p < 0.05) were observed in certain cases. Among the isolates obtained from male patients, a higher prevalence of resistance to cefoxitin (FOX), ertapenem (ETP), and imipenem (IMP) was recorded, each with a 67% prevalence in phylogenetic group A, compared to other groups. In female patients, a greater number of carbapenemase-producing isolates (62.5%) was identified in group A, as well as higher resistance to gentamicin (GM) in group B2 (57%), to ampicillin (AMP) in group D (96%), and to levofloxacin (LVX) and ciprofloxacin (CIP) in the undetermined phylogenetic group (100%). All these associations were statistically significant (p < 0.05).

3.5 Eric PCR

To evaluate the potential clonal relationship among clinical uropathogenic E. coli isolates, a clustering analysis was performed based on profiles obtained by ERIC-PCR. This analysis was organized at three levels: by hospital ward, among non-hospitalized patients, and among hospitalized patients. At the hospital ward level, possible clonal relationships were identified among isolates from eight clinical departments. The wards with the highest number of clusters were Emergency (5 clusters; 13 strains), Internal Medicine (7 clusters; 14 strains), Gynecology (2 clusters; 6 strains), General Medicine (2 clusters; 4 strains), Pathological Puerperium (1 cluster; 2 strains), Endocrinology (1 cluster; 2 strains), Pediatrics (1 cluster; 2 strains), Urology (2 clusters; 4 strains), and Outpatients (1 cluster; 2 strains) (Supplementary Figures 5A–P). Among hospitalized patients, 58 isolates were grouped into 20 clusters (Supplementary Figures O,P), while in non-hospitalized patients, 39 isolates were grouped into 16 clusters (Supplementary Figure 5P).

Despite the high similarity in banding patterns obtained through ERIC-PCR, we observed variations in resistance and virulence profiles, and even in phylogenetic grouping. In this context, it is important to consider that ERIC sequences do not differentiate specific virulence or antibiotic resistance genes, which can be acquired through horizontal gene transfer events. Similarly, phylogenetic groups are determined by specific genes that cannot be distinguished via ERIC-PCR. However, we identified some isolates with features suggestive of possible clonality: they showed identical ERIC profiles according to our cut-off criteria, along with high similarity in both virulence and antibiotic resistance profiles. These isolates were detected by hospital ward (Figure 8A), in the group of hospitalized patients (Figure 8B), and the outpatient group (Figure 8C). Figure 8 visually summarizes some of the clusters with the highest similarity, integrating genomic, phenotypic, and epidemiological profiles.

Figure 8
Three phylogenetic trees labeled A, B, and C, display genetic relationships among UPEC strains. Each tree includes a matrix of red and blue circles indicating the presence or absence of genetic features, respectively. The accompanying tables list strain details, conditions, gender, and age groups. Star icons denote specific attributes at the far right. Tree scales vary, with A at scale 1, while B and C are at scale 0.1.

Figure 8. Clinical isolates with the highest similarity, according to clustering criteria, virulence profiles, and antibiotic resistance phenotypes. (A) Analysis by hospital ward; (B) Analysis in hospitalized patients; (C) Analysis in outpatients. Isolates shown in blue exhibit the highest similarity in virulence, resistance, and phylogenetic groups.

4 Discussion

Urinary tract infections (UTIs) pose a significant public health challenge due to their high incidence and the growing prevalence of multidrug-resistant bacterial strains (Magiorakos et al., 2012; Flores-Mireles et al., 2015). Uropathogenic Escherichia coli (UPEC) is the main causative agent, noted for its high capacity to acquire both virulence and antibiotic resistance genes (Ballesteros-Monrreal et al., 2021; Akhand and Karicheri, 2022). In Honduras, information on UPEC is scarce; therefore, characterizing local strains is essential to improve microbiological diagnosis, empirical treatment, and epidemiological surveillance. In this study, 126 clinical isolates confirmed as E. coli were analyzed, of which 99 (78.6%) were obtained from women and 27 (21.4%) from men. This higher proportion in women aligns with the literature, which points to anatomical and hormonal factors as predisposing elements for UTIs (Flores-Mireles et al., 2015; Raphael et al., 2024). These factors include the shorter female urethra compared to that of males, as well as the proximity of the anus and vagina, which facilitates ascending colonization (He et al., 2025; Deltourbe et al., 2022).

One of the findings in this study was that most isolates came from patients aged 13 to 50, followed by patients over 50, with no significant differences between sexes. In the 13–50 age group, events such as sexual activity, use of contraceptives, pregnancy, and in some cases prostatic hyperplasia occur more frequently—factors that are all associated with a higher risk of UTI (He et al., 2025; Deltourbe et al., 2022). Although no detailed clinical information was collected from patients in this study, these factors have been widely described as risk elements for disease development. Therefore, their influence cannot be ruled out and should be considered in future research that integrates clinical, microbiological, and demographic variables more comprehensively.

4.1 Virulence determinants

On the other hand, the analysis of virulence genes revealed significant differences between clinical isolates from male and female patients. Overall, strains from male patients exhibited a significantly higher virulence gene burden, with an average of 10.48 genes per strain, compared to 8.06 genes in female patients (p = 0.0029). This difference suggests that infectious episodes in male patients may be associated with strains possessing a greater pathogenic potential, which aligns with previous studies that have described a higher requirement for adaptability to colonize the male urinary tract, particularly in the presence of obstructive factors or urinary catheterization (Ruiz et al., 2002; AP et al., 2017).

Regarding age group distribution, results in females showed a higher mean number of genes in patients under 5 years old; however, this difference did not reach statistical significance. In males, the 6-to-12-year-old group stood out with a significantly higher genetic load than that observed in adolescents and adults (p < 0.05). This finding suggests that, in specific pediatric subgroups, pathogenesis may involve a higher number of virulence factors, which could translate into more aggressive clinical presentations or increased risk of recurrence. The elevated virulence gene load in isolates from pediatric patients has been previously reported by other authors (Hernández-Chiñas et al., 2019; Luna-Pineda et al., 2018), and although the sample size in our study limits definitive conclusions, these data suggest that the pediatric population, particularly school-aged boys, represents a risk group that warrants closer clinical and microbiological monitoring. Implementing follow-up strategies in this cohort could be key to detecting emerging virulence patterns, preventing complications, and timely adjusting therapeutic decisions.

Among the most prevalent genes, fimH (96%) and feoB (100%) were identified, both essential for initial adherence to the uroepithelium and iron acquisition, respectively, as well as for urinary tract colonization, particularly in nutrient-limited environments (Alam Parvez and Rahman, 2018). The high frequency of these genes is consistent with their central role in the pathophysiology of UPEC and has been widely reported (Ballesteros-Monrreal et al., 2021; Zhao et al., 2009; Yun et al., 2014; Tabasi et al., 2015; Sarshar et al., 2020). However, when stratified by sex, significant differences were observed in the prevalence of genes with more specialized functions. In males, there was a higher presence of genes related to renal-level adherence (papG-II, papC), toxins (hlyA, cnf-1), capsule formation (kpsM), and iron acquisition systems (fyuA, iucD, iha), all of which are associated with increased persistence and tissue damage in complicated infections (Ballesteros-Monrreal et al., 2021).

These results are consistent with previous reports: Kudinha et al. (2013) described a virulence gradient with higher prevalence of virulence factors in pyelonephritis compared to cystitis and healthy controls; Firoozeh et al. (2014) showed that pyelonephritis isolates harbor a greater diversity of virulence genes, particularly pap, traT, aer, hly, and PAIs; Hyun et al. (2022) identified sex-based differences with male isolates enriched in toxin genes such as hlyA and cnf1; and Hsiao et al. (2024) demonstrated that phylogroup B2 predominates in UTI and bacteremic cases, with genes such as papGII enriched in invasive infections (Hyun et al., 2022; Hsiao et al., 2024; Firoozeh et al., 2014; Kudinha et al., 2013). Collectively, these findings support our observation that isolates from male patients harbor more virulence determinants and are also associated with complicated clinical outcomes.

Regarding the virulence profiles of the isolates, high variability was observed, with a total of 92 distinct profiles: 14 were shared among multiple isolates from female patients, and only 2 among male isolates. Notably, isolates with shared profiles were not concentrated in a specific hospital ward or age group, suggesting a random and unfocused distribution of genetically similar strains within the hospital environment, with no clear evidence of clonal clustering within this dataset. Additionally, the low frequency of shared profiles is consistent with the high heterogeneity observed in the ERIC-PCR patterns, suggesting that most cases correspond to independent infectious events and a high diversity in infection sources. Nevertheless, ERIC-PCR results revealed two isolates that showed high similarity in their virulence characteristics, resistance, phylogenetic group, and hospital origin (Figure 8C). This finding underscores the importance of implementing continuous molecular surveillance to identify potential outbreaks, common transmission routes, or the emergence of successful clones with predominant virulence profiles.

Likewise, the high heterogeneity observed in virulence profiles may be related to the acquisition and maintenance of genes through horizontal transfer mechanisms, such as plasmids and pathogenicity islands (PAIs), reflecting the remarkable adaptability of uropathogenic isolates to diverse ecological niches and selective pressures. In our analysis, the significant co-occurrence of genes such as papG-II, papC, sfaD/focC, vat, fyuA, kpsM, and iucD is consistent with the findings of Ballesteros-Monrreal et al. (2021), who identified similar clusters in uropathogenic E. coli strains. More recently, Ballesteros-Monrreal et al. (2024) corroborated the presence of these PAIs through molecular characterization, suggesting that the observed gene clustering in our isolates likely reflects the functional and structural organization of these genetic elements rather than random distribution. This finding supports the hypothesis that the simultaneous presence of multiple virulence factors may be mediated by shared genomic structures. Moreover, considering that genes such as iutA and iucD have previously been associated with ColV-type plasmids (Johnson et al., 2006; Johnson and Nolan, 2009), it is reasonable to assume that some of our isolates may carry plasmids with horizontal dissemination capacity. This aspect is particularly relevant in hospital settings, where the circulation of strains harboring PAIs and plasmids contributes not only to pathogenic potential but also to the spread of antimicrobial resistance. Future studies should aim to confirm the presence and characterize the structure of such plasmids to better understand their role in the dissemination of resistance and virulence determinants. Ideally through whole-genome sequencing (WGS) to resolve their genetic context and evolutionary dynamics.

4.2 Antimicrobial resistance

Regarding the characterization of the isolates based on antibiotic resistance, the analysis of the resistance index (RI) revealed significant differences between clinical isolates according to the patient’s sex. Isolates from male patients exhibited a significantly higher antimicrobial resistance burden, with an average RI of 0.46, compared to 0.27 in female-derived isolates (p < 0.0001). Additionally, 69% of all isolates were classified as MDR (including XDR), with the highest proportion of MDR found in male patients (88%; p = 0.009). This difference reinforces the previously observed trend in virulence gene burden, suggesting that strains isolated from male patients not only possess greater pathogenic potential but also exhibit increased antimicrobial resistance. However, one of the limitations of our study is the smaller male sample size (n = 27), which restricts the possibility of drawing definitive conclusions in this subgroup, although the observed trends are consistent with international evidence (Ruiz et al., 2002; AP et al., 2017; Khanal et al., 2024).

These results are consistent with findings by Khanal et al. (2024), who analyzed over 85,000 uropathogenic E. coli isolates in Australia and found a high prevalence of MDR strains in both hospital and community settings, with a significantly higher proportion in males compared to females (16.5% vs. 12.8% in hospitals; 6.4% vs. 5.2% in the community; p < 0.001). Although the Australian study also lacked individual clinical metadata, the authors attributed these differences to the fact that urinary tract infections in males are more frequently associated with conditions that complicate the clinical course, such as benign prostatic hyperplasia, urolithiasis, use of urinary catheters, history of urological instrumentation, or prior surgical interventions (Khanal et al., 2024). These conditions lead to urinary flow obstruction, increased risk of persistent colonization, and the need for prolonged or broad-spectrum antibiotic treatments, favoring the selection of resistant strains. Similar to our findings, Khanal et al. (2024) noted that, despite UTIs being more frequent in females, the strains isolated from males showed higher multidrug resistance rates (Khanal et al., 2024). This disparity suggests that UTIs in males may be caused by strains with more complex virulence and resistance profiles. This finding underscores the need to adapt empirical treatments not only to the type of infection but also to the epidemiological features associated with patient sex.

Among the 17 antibiotics tested, the highest prevalence of resistance in the 126 strains was observed for AMP (80%), TSX (67%), fluoroquinolones (CIP and LVX), and several second- to fourth-generation cephalosporins. When stratified by sex, isolates from male patients showed higher resistance rates for all evaluated antibiotics, with statistically significant differences observed for aminoglycosides, simple and combination penicillins, cephalosporins, carbapenems, nitrofurans, and fluoroquinolones. Similarly, the ESBL phenotype was more prevalent in this group (59%; p = 0.049). These results are particularly relevant when contrasted with currently recommended empirical treatment regimens in Honduras (Honduras, 2019; Ucros, 2009; Giovanni Erazo Trimarchi et al., 2015; Honduras, 2016). For outpatient pediatric patients, first-line treatment includes TSX and AMC, whereas in hospital settings, AMP, GM, AMK, CRO, or ampicillin-sulbactam (AMS) are commonly used (Johnson et al., 2006; Johnson and Nolan, 2009; Ucros, 2009). However, our findings reveal a high prevalence of resistance to TSX (60–100%), AMP (83–85%), and AMC (27–29%) in pediatric isolates of both sexes (with resistance reaching 66.7% in males and 40% in females), as well as to GM (33.3%) in males and CRO (33.3%) in females (Supplementary material 2B), suggesting potentially limited efficacy of these therapeutic options in the pediatric population. This situation is particularly concerning in children under 5 years of age, not only due to the frequent use of these antibiotics as first-line treatments but also due to the detection of a high burden of virulence genes. The coexistence of these factors may compromise therapeutic success and worsen clinical prognosis. Additionally, among pediatric female patients, ESBL- and CAR-producing isolates were observed, whereas none were detected in male children (Supplementary material 2B). This difference may be explained by the higher frequency of UTIs in females, which increases exposure to antibiotic treatments and consequently the selective pressure favoring the emergence of MDR strains.

In adult outpatients, commonly used antibiotics include AMC, TSX, CIP, LVX, and CRO; in hospitalized patients, fluoroquinolones, CRO, AMS, TZP, and carbapenems such as ertapenem (ETP) are used (Giovanni Erazo Trimarchi et al., 2015; Loscalzo et al., 2022). In our study, high resistance prevalence was observed among adult patients against these antimicrobials, particularly TSX (60–83.3%), CIP (39–93.75%), and CRO (31–83%) (Supplementary material 2B). This resistance pattern is not unique to our setting: studies conducted in the USA and Australia have reported elevated resistance rates to these drugs, especially TSX, more commonly observed in patients with recurrent UTIs, anatomical abnormalities of the urinary tract, or recent antibiotic treatment history, which has raised concerns regarding its use as a first-line empirical therapy (Khanal et al., 2024; Wesolek et al., 2022). These findings also align with a systematic review and meta-analysis conducted in Mexico, which reported high resistance to ampicillin (79%), third-generation cephalosporins (40–50%), quinolones and fluoroquinolones (45–61%), as well as TSX (61.9%), along with a considerable prevalence of ESBL-producing isolates (39.46%) (Ballesteros-Monrreal et al., 2023). This body of evidence underscores the urgent need to reassess empirical treatment guidelines for UTIs, considering both selective pressure and local resistance contexts, which may vary significantly depending on geographic region.

The high prevalence of resistance to these antibiotics may be linked to well-established genetic mechanisms. For TSX, resistance has been associated with dfrA and sul genes, commonly found in mobile integrons (Poey et al., 2024). In the case of CIP, resistance can result from chromosomal mutations in gyrA and parC, as well as plasmid-mediated mechanisms involving efflux pumps (qepA), target protection (qnrB, qnrA, qnrS), or enzymatic modification [aac(6′)-Ib-cr] (Malekzadegan et al., 2019). Regarding CRO, one of the most frequent mechanisms is the production of extended-spectrum β-lactamases (ESBLs), which confer resistance to multiple cephalosporins. Notably, many of these determinants can coexist within conjugative plasmids that harbor resistance genes for various antimicrobial classes (Salah et al., 2019; Juraschek et al., 2022), facilitating co-selection phenomena.

The co-occurrence of antimicrobial resistance phenotypes in our isolates is unlikely to be random, but rather reflects their organization within mobile genetic elements such as integrons and conjugative plasmids, which frequently harbor multiple resistance determinants within the same replicon. This genetic architecture facilitates horizontal gene transfer and promotes co-selection, whereby exposure to a single antibiotic may indirectly maintain resistance to several unrelated drugs. Similar patterns have been described in Enterobacteriaceae worldwide, underscoring the role of these mobile elements in the persistence and dissemination of multidrug resistance. Future studies using whole-genome sequencing (WGS) will be essential to characterize the specific plasmids and integrons circulating in Honduras and to resolve the genetic context underlying these associations.

Although the presence of these genes was not directly assessed, a strong correlation was observed between resistance to ceftriaxone and ciprofloxacin (r = 0.433; p < 0.001), suggesting possible co-transmission of these phenotypes. Additionally, significant associations were found between the ESBL phenotype and resistance to both cephalosporins and fluoroquinolones (r = 0.355; p < 0.001), potentially indicating the presence of mobile genetic elements that integrate and disseminate gene clusters associated with multidrug resistance in the studied population, posing a considerable therapeutic threat, particularly in settings where these antimicrobials are frequently used empirically.

This pattern of co-dissemination of resistance mechanisms has been previously documented in recent studies conducted in Mexico. In clinical isolates from women with UTIs, Ballesteros-Monrreal et al. (2024) reported a high prevalence of blaCTX-M-2 and aac(6′)-Ib genes in isolates with simultaneous resistance phenotypes to β-lactams, fluoroquinolones, aminoglycosides, and ESBL production. Similarly, Marquez-Salazar et al. (2025) found that more than 95% of community-acquired E. coli isolates resistant to both β-lactams and fluoroquinolones carried ESBL and plasmid-mediated quinolone resistance (PMQR) genes—mainly aac(6′)-Ib-cr, qnrB, and qepA—pointing to the presence of conjugative plasmids carrying multiple resistance genes. The coexistence of these phenotypes may offer clues regarding potential resistance genotypes.

A limitation of our study is that we did not specifically detect these genetic determinants, highlighting the need for future investigations aimed at identifying genes associated with mobile genetic elements, as well as the circulating plasmids and integrons in E. coli strains from Honduras. Such research is essential to anticipate therapeutic efficacy crises and guide rational antibiotic use policies. Collectively, the findings of this study reveal a concerning landscape of antimicrobial resistance in uropathogenic E. coli strains in the Honduran context, necessitating a reevaluation of current therapeutic protocols and the strengthening of local surveillance. Integrating phenotypic resistance data with epidemiological and molecular information should be a priority for clinical decision-making and the formulation of antibiotic stewardship policies. In this regard, the use of personalized antibiograms, continuous training of medical personnel, and the incorporation of local data into clinical guidelines are key pillars for containing the spread of multidrug-resistant strains and preserving the effectiveness of available antimicrobials.

4.3 Phylogeny and clonality

The phylogenetic characterization of the 126 clinical E. coli isolates revealed a heterogeneous distribution, with a predominance of phylogroup B2 (29%), followed by A and D (21% each)—lineages commonly associated with extraintestinal infections. We also identified phylogroups B1, Clade I, F, G, and C, along with 5% of non-typeable strains. Notably, phylogroup G, recently incorporated into the scheme proposed by Clermont et al. (2019), was detected. Although rare, its clinical relevance cannot be ruled out due to its shared characteristics with B2. The application of this updated scheme may explain the lower proportion of non-typeable strains compared to previous studies. Nevertheless, the persistence of unclassified isolates suggests the possible presence of phylogenetic lineages not yet described in Honduras.

Regarding sex-specific distribution, greater phylogenetic diversity was observed among female patients, with predominance of groups A and D, whereas phylogroup B2 was significantly more prevalent in male patients (59%), potentially associated with a higher burden of virulence and resistance genes, as well as anatomical and clinical factors specific to males. Although no significant differences were found among hospital wards, we identified a trend toward the concentration of B2 strains in departments such as Internal Medicine and Urology, possibly reflecting their association with more complex clinical presentations. Virulence index analysis showed significantly higher values in phylogroup B2 compared to A and B1, with associations to key genes related to adhesion, cytotoxicity, immune evasion, and iron acquisition, such as papG-II, papC, sat, hlyA, vat, cnf-1, kpsM, fyuA, iucD, and iha. Although phylogroup F exhibited a high median index, the low frequency of isolates precludes definitive conclusions but suggests an underestimated pathogenic potential that warrants further investigation.

Regarding antimicrobial resistance, no significant differences were observed in the overall index across phylogroups, although sex-stratified patterns emerged: in males, phylogroup A showed higher resistance to carbapenems; in females, B2 exhibited elevated resistance to gentamicin, D to ampicillin, and non-typeable strains to fluoroquinolones. These findings point to the existence of phylogroup-specific profiles possibly modulated by distinct selective pressures. Our findings align with those of Ballesteros-Monrreal et al. (2021), who also reported high phylogenetic diversity in clinical isolates from women, with 25% of strains being non-typeable. In that study, phylogroup B2 harbored the highest virulence gene load (iroN 91%, papC 64%, fyuA 100%), whereas B1 harbored the greatest antimicrobial resistance. Similarly, Hemati et al. (2024) identified B2 as the most prevalent group (32%), followed by E, C, and B1, and reported that over half of the isolates could not be phylogenetically assigned. Kara et al. (2024) also highlighted B2 as the dominant phylogroup (45.5%), followed by E (22.4%) and 12% non-typeable strains. Complementarily, Wang et al. (2023) characterized 844 clinical isolates and confirmed the predominance of phylogroup B2 (64.8%), with the highest average number of virulence genes (10.1 per isolate). Moreover, although phylogroups C, E, and F were less frequent, they displayed high multidrug resistance rates, supporting the importance of monitoring minority lineages that may be acquiring clinically relevant traits.

4.4 Implications and future directions

Collectively, our data identify phylogroup B2 as the predominant uropathogenic lineage in this cohort, supported by its association with multiple virulence genes and prominent role in resistance profiles. These findings highlight the need to strengthen molecular surveillance of uropathogenic E. coli in response to increasing genetic diversity and the potential emergence of clinically significant new lineages. This study presents one of the first comprehensive molecular and phenotypic characterizations of UPEC in Honduras. Our findings demonstrate a high prevalence of virulence factors, MDR, and ESBL production among isolates from a tertiary-care hospital, with a particular impact on pediatric and male patient populations. The observed resistance to commonly prescribed antibiotics, including ampicillin, trimethoprim-sulfamethoxazole, fluoroquinolones, and cephalosporins, raises significant concerns regarding the continued effectiveness of empiric treatment protocols recommended in national guidelines.

Phylogenetic analysis revealed substantial genetic diversity, with phylogroup B2 emerging as the predominant lineage, associated with increased virulence and resistance. The ERIC-PCR profiles, in combination with overlapping resistance and virulence patterns, suggest potential clonal relationships and highlight the urgent need for enhanced molecular surveillance strategies. In a context where data on urinary tract infections remain limited, our results provide critical evidence to guide clinical decision-making, support antimicrobial stewardship initiatives, and reinforce microbiological surveillance at the local level.

Another limitation of our study is the short sampling period (2 months), which restricts epidemiological generalizability. Future investigations should extend the collection period to capture seasonal variation and temporal dynamics in UPEC populations. Recognizing this limitation underscores the importance of complementing cross-sectional molecular studies with longitudinal surveillance to fully appreciate the epidemiological behavior of uropathogenic E. coli in our setting.

Finally, one of the future goals of our research is the complete sequencing of the genomes of isolates exhibiting co-occurrence of virulence and resistance genes. Whole-genome sequencing (WGS) represents a more robust approach compared to PCR-based detection, as it allows for the comprehensive reconstruction of the genomic context in which these determinants are embedded. Through WGS, it will be possible to confirm the presence of mobilizable genetic elements such as conjugative plasmids, transposons, and pathogenicity islands, and to elucidate their structural organization and evolutionary dynamics. This level of resolution is crucial to assess their potential for horizontal transfer within the hospital setting, where the circulation of highly adaptable clones represents a major epidemiological threat. In addition, characterizing these isolates in such depth provides a unique opportunity to move beyond descriptive surveillance and into translational applications. WGS-based analysis will enable the identification of the most relevant virulence factors associated with urinary tract infections in Honduras, including those that directly contribute to host colonization, immune evasion, or tissue damage. Understanding these genetic traits in a population-specific context may guide the design of novel preventive or therapeutic strategies, either by prioritizing the development of vaccines, antivirulence drugs, or adjuvant therapies that target key determinants of pathogenicity as has been observed in previously published studies focused on other relevant models (Tiwari et al., 2014; Khan et al., 2023; Sohail et al., 2025; Bunduki et al., 2021). Consequently, the genomic dissection of local UPEC populations does not only represent a step toward a more precise understanding of their epidemiology but also opens avenues for the identification of clinically actionable targets that could improve the management of UTIs in vulnerable groups.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

This study involved the analysis of phenotypic and sociodemographic data previously obtained by the hospital, as well as the molecular characterization of Escherichia coli strains donated by a clinical laboratory. No patient contact occurred, and no identifiable personal data were accessed. Therefore, ethical approval and informed consent were not required.

Author contributions

BO: Formal analysis, Writing – original draft, Conceptualization, Data curation, Investigation, Writing – review & editing. MB-M: Formal analysis, Writing – review & editing, Writing – original draft, Investigation, Data curation, Conceptualization. CG: Data curation, Investigation, Writing – review & editing. DR: Writing – review & editing, Data curation, Investigation. LR: Investigation, Writing – review & editing. FP: Writing – review & editing, Investigation. VM: Investigation, Writing – review & editing. EV: Writing – review & editing, Investigation. AE: Writing – review & editing, Investigation. PM-P: Writing – review & editing, Formal analysis, Visualization. LE: Data curation, Investigation, Writing – review & editing. DV: Investigation, Writing – review & editing. GF: Conceptualization, Supervision, Funding acquisition, 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 partially funded by DICIHT-UNAH.

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 Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.1656938/full#supplementary-material

Footnotes

References

Agarwal, J., Srivastava, S., and Singh, M. (2012). Pathogenomics of uropathogenic Escherichia coli. Indian J. Med. Microbiol. 30, 141–149. doi: 10.4103/0255-0857.96657

PubMed Abstract | Crossref Full Text | Google Scholar

Akhand, A., and Karicheri, R. (2022). Antibiotic resistance among Uropathogenic Escherichia Coli. Eur. J. Molecular Clin. Medicine. 9, 4914–4920. Available at: https://link.gale.com/apps/doc/A765922960/HRCA?u=anon~193888f7&sid=googleScholar&xid=431b774c

Google Scholar

Alam Parvez, M. S., and Rahman, D. (2018). “Virulence factors of Uropathogenic E. coli” in Microbiology of urinary tract infections-microbial agents and predisposing factors. ed. P. Behzadi (Rijeka: IntechOpen).

Google Scholar

Alsultan, A., and Elhadi, N. (2022). Evaluation of ERIC-PCR method for determining genetic diversity among Escherichia coli isolated from human and retail imported frozen shrimp and beef. Int. J. Food Contam. 9:12. doi: 10.1186/s40550-022-00098-1

Crossref Full Text | Google Scholar

AP, d. S. d.-S., de Sousa, V. S., Martins, N., da Silva Dias, R. C., Bonelli, R. R., Riley, L. W., et al. (2017). Escherichia coli sequence type 73 as a cause of community acquired urinary tract infection in men and women in Rio de Janeiro, Brazil. Diagn. Microbiol. Infect. Dis. 88, 69–74. doi: 10.1016/j.diagmicrobio.2017.01.024

Crossref Full Text | Google Scholar

Ardakani, M. A., and Ranjbar, R. (2016). Molecular typing of uropathogenic E. coli strains by the ERIC-PCR method. Electron. Physician 8, 2291–2295. doi: 10.19082/2291

PubMed Abstract | Crossref Full Text | Google Scholar

Ballesteros-Monrreal, M. G., Arenas-Hernández, M. M., Barrios-Villa, E., Juarez, J., Álvarez-Ainza, M. L., Taboada, P., et al. (2021). Bacterial morphotypes as important trait for uropathogenic E. coli diagnostic; a virulence-phenotype-phylogeny study. Microorganisms. 9:2381. doi: 10.3390/microorganisms9112381

PubMed Abstract | Crossref Full Text | Google Scholar

Ballesteros-Monrreal, M. G., Mendez-Pfeiffer, P., Barrios-Villa, E., Arenas-Hernández, M. M., Enciso-Martínez, Y., Sepúlveda-Moreno, C. O., et al. (2023). Uropathogenic Escherichia coli in Mexico, an overview of virulence and resistance determinants: systematic review and meta-analysis. Arch. Med. Res. 54, 247–260. doi: 10.1016/j.arcmed.2023.01.001

PubMed Abstract | Crossref Full Text | Google Scholar

Ballesteros-Monrreal, M. G., Mendez-Pfeiffer, P., Ortíz, B., Bolado-Martínez, E., Álvarez-Ainza, M. L., Enciso-Martínez, Y., et al. (2024). Uropathogenic E. Coli and hybrid Pathotypes in Mexican women with urinary tract infections: a comprehensive molecular and phenotypic overview. Curr. Issues Mol. Biol. 46, 5909–5928. doi: 10.3390/cimb46060353

PubMed Abstract | Crossref Full Text | Google Scholar

Behzadi, P., and Carevic, B. (2019). Microbiology of urinary tract infections: Microbial agents and predisposing factors. London: BoD–Books on Demand.

Google Scholar

Bunduki, G. K., Heinz, E., Phiri, V. S., Noah, P., Feasey, N., and Musaya, J. (2021). Virulence factors and antimicrobial resistance of uropathogenic Escherichia coli (UPEC) isolated from urinary tract infections: a systematic review and meta-analysis. BMC Infect. Dis. 21, 1–13. doi: 10.1186/s12879-021-06435-7

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Clermont, O., Dixit, O. V. A., Vangchhia, B., Condamine, B., Dion, S., Bridier-Nahmias, A., et al. (2019). Characterization and rapid identification of phylogroup G in Escherichia coli, a lineage with high virulence and antibiotic resistance potential. Environ. Microbiol. 21, 3107–3117. doi: 10.1111/1462-2920.14713

PubMed Abstract | Crossref Full Text | Google Scholar

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, 1–8. doi: 10.1186/s12866-017-0979-5

Crossref Full Text | Google Scholar

de Kraker, M. E., Stewardson, A. J., and Harbarth, S. (2016). Will 10 million people die a year due to antimicrobial resistance by 2050? PLoS Med. 13:e1002184. doi: 10.1371/journal.pmed.1002184

PubMed Abstract | Crossref Full Text | Google Scholar

Deltourbe, L., Mariano, L. L., Hreha, T. N., Hunstad, D. A., and Ingersoll, M. A. (2022). The impact of biological sex on diseases of the urinary tract. Mucosal Immunol. 15, 857–866. doi: 10.1038/s41385-022-00549-0

PubMed Abstract | Crossref Full Text | Google Scholar

Dishan, A., Hizlisoy, H., Barel, M., Disli, H., Gungor, C., Ertas Onmaz, N., et al. (2023). Biofilm formation, antibiotic resistance and genotyping of Shiga toxin-producing Escherichia coli isolated from retail chicken meats. Br. Poult. Sci. 64, 63–73. doi: 10.1080/00071668.2022.2116697

PubMed Abstract | Crossref Full Text | Google Scholar

Firoozeh, F., Saffari, M., Neamati, F., and Zibaei, M. (2014). Detection of virulence genes in Escherichia coli isolated from patients with cystitis and pyelonephritis. Int. J. Infect. Dis. 29, 219–222. doi: 10.1016/j.ijid.2014.03.1393

PubMed Abstract | Crossref Full Text | Google Scholar

Flores-Mireles, A. L., Walker, J. N., Caparon, M., and Hultgren, S. J. (2015). Urinary tract infections: epidemiology, mechanisms of infection and treatment options. Nat. Rev. Microbiol. 13, 269–284. doi: 10.1038/nrmicro3432

PubMed Abstract | Crossref Full Text | Google Scholar

Giovanni Erazo Trimarchi, G. J. P., Espinoza, M., Aceituno, N., Solórzano, L., Álvarez, A., Caballero, A., et al. (2015). Manual de medicina interna Diagnostico y tratamiento. Postgrado de Medicina Interna Promoción 2015 FdCM. Tegucigalpa: Universidad Nacional Autonoma de Honduras.

Google Scholar

Halaji, M., Fayyazi, A., Rajabnia, M., Zare, D., Pournajaf, A., and Ranjbar, R. (2022). Phylogenetic group distribution of uropathogenic Escherichia coli and related antimicrobial resistance pattern: a meta-analysis and systematic review. Front. Cell. Infect. Microbiol. 12:790184. doi: 10.3389/fcimb.2022.790184

PubMed Abstract | Crossref Full Text | Google Scholar

He, Y., Zhao, J., Wang, L., Han, C., Yan, R., Zhu, P., et al. (2025). Epidemiological trends and predictions of urinary tract infections in the global burden of disease study 2021. Sci. Rep. 15:4702. doi: 10.1038/s41598-025-89240-5

PubMed Abstract | Crossref Full Text | Google Scholar

Hemati, S., Halimi, S., Jabalameli, F., Emaneini, M., and Beigverdi, R. (2024). Phylogenetic group, antibiotic resistance, virulence gene, and genetic diversity of Escherichia coli causing bloodstream infections in Iran. Front. Microbiol. 15:1426510. doi: 10.3389/fmicb.2024.1426510

PubMed Abstract | Crossref Full Text | Google Scholar

Hernández-Chiñas, U., Pérez-Ramos, A., Belmont-Monroy, L., Chávez-Berrocal, M. E., González-Villalobos, E., Navarro-Ocaña, A., et al. (2019). Characterization of auto-agglutinating and non-typeable uropathogenic Escherichia coli strains. J. Inf. Dev. Countries. 13, 465–472. doi: 10.3855/jidc.11098

PubMed Abstract | Crossref Full Text | Google Scholar

Honduras, SD. (2016). Protocolos para la atención durante la preconcepción, el embarazo, el parto, el puerperio y del neonato PT04: 2016. Tegucigalpa M.D.C.: Secretaría de Salud de Honduras. Available at: https://honduras.bvsalud.org/wp-content/uploads/2024/04/PROTOCOLOS-PARA-LA-ATENCION-DURANTE-LA-PRECONCEPCION-EL-EMBARAZO-EL-PARTO-EL-PUERPERIO-Y-DEL-NEONATO-VOLUMEN-2-ATENCION-DEL-PARTO-Y-DEL-NEONATO.pdf

Google Scholar

Honduras, SD (2019). Protocolo para la atencion del menor de 5 años en el primer nivel, PT39: 2019. Tegucigalpa M.D.C.: Secretaría de Salud de Honduras. Available at: https://honduras.bvsalud.org/wp-content/uploads/2024/04/Protocolo-de-Atencion-Integral-a-los-y-las-Adolescentes-en-el-Primer-Nivel-de-Atencion.pdf

Google Scholar

Hsiao, C.-Y., Lee, Y.-C., Shyu, D. J., Su, C.-P., Lin, M.-Y., Guo, N.-Y., et al. (2024). Molecular characterization of high and low virulent Escherichia coli clinical strains isolated from patients with urinary tract infections with or without bacteremia in southern Taiwan. Inf. Drug Resist. 17, 2389–2399. doi: 10.2147/IDR.S458925

PubMed Abstract | Crossref Full Text | Google Scholar

Hyun, M., Lee, J. Y., and Kim, H. A. (2022). Comparison of the clinical and genotypic characteristics of uropathogenic Escherichia coli strains according to sex in Korea. Microb. Drug Resist. 28, 988–996. doi: 10.1089/mdr.2021.0144

PubMed Abstract | Crossref Full Text | Google Scholar

Johnson, T. J., and Nolan, L. K. (2009). Pathogenomics of the virulence plasmids of Escherichia coli. Microbiol. Mol. Biol. Rev. 73, 750–774. doi: 10.1128/MMBR.00015-09

PubMed Abstract | Crossref Full Text | Google Scholar

Johnson, T. J., Siek, K. E., Johnson, S. J., and Nolan, L. K. (2006). DNA sequence of a ColV plasmid and prevalence of selected plasmid-encoded virulence genes among avian Escherichia coli strains. J. Bacteriol. 188, 745–758. doi: 10.1128/JB.188.2.745-758.2006

PubMed Abstract | Crossref Full Text | Google Scholar

Juraschek, K., Malekzadah, J., Malorny, B., Käsbohrer, A., Schwarz, S., Meemken, D., et al. (2022). Characterization of qnrB-carrying plasmids from ESBL-and non-ESBL-producing Escherichia coli. BMC Genomics 23:365. doi: 10.1186/s12864-022-08564-y

PubMed Abstract | Crossref Full Text | Google Scholar

Kang, C.-I., Kim, J., Park, D. W., Kim, B.-N., Ha, U.-S., Lee, S.-J., et al. (2018). Clinical practice guidelines for the antibiotic treatment of community-acquired urinary tract infections. Infection Chemother. 50, 67–100. doi: 10.3947/ic.2018.50.1.67

PubMed Abstract | Crossref Full Text | Google Scholar

Kara, A., Massaro, C., Giammanco, G. M., Alduina, R., and Boussoualim, N. (2024). Phylogenetic diversity, antibiotic resistance, and virulence of Escherichia coli strains from urinary tract infections in Algeria. Antibiotics. 13:773. doi: 10.3390/antibiotics13080773

PubMed Abstract | Crossref Full Text | Google Scholar

Karim, M. R., Zakaria, Z., Hassan, L., Mohd Faiz, N., and Ahmad, N. I. (2023). Antimicrobial resistance profiles and co-existence of multiple antimicrobial resistance genes in mcr-Harbouring Colistin-resistant Enterobacteriaceae isolates recovered from poultry and poultry meats in Malaysia. Antibiotics. 12:1060. doi: 10.3390/antibiotics12061060

PubMed Abstract | Crossref Full Text | Google Scholar

Khan, A., Sohail, S., Yaseen, S., Fatima, S., Wisal, A., Ahmed, S., et al. (2023). Exploring and targeting potential druggable antimicrobial resistance targets ArgS, SecY, and MurA in Staphylococcus sciuri with TCM inhibitors through a subtractive genomics strategy. Funct. Integr. Genomics 23:254. doi: 10.1007/s10142-023-01179-w

PubMed Abstract | Crossref Full Text | Google Scholar

Khanal, N., Cortie, C. H., Story, C., Jones, S., Mansfield, K. J., Miyakis, S., et al. (2024). Multidrug resistance in urinary E. coli higher in males compared to females. BMC Urol. 24:255. doi: 10.1186/s12894-024-01654-x

PubMed Abstract | Crossref Full Text | Google Scholar

Kucheria, R., Dasgupta, P., Sacks, S., Khan, M., and Sheerin, N. (2005). Urinary tract infections: new insights into a common problem. Postgrad. Med. J. 81, 83–86. doi: 10.1136/pgmj.2004.023036

PubMed Abstract | Crossref Full Text | Google Scholar

Kudinha, T., Johnson, J., Andrew, S., Kong, F., Anderson, P., and Gilbert, G. (2013). Distribution of phylogenetic groups, sequence type ST131, and virulence-associated traits among Escherichia coli isolates from men with pyelonephritis or cystitis and healthy controls. Clin. Microbiol. Infect. 19, E173–E180. doi: 10.1111/1469-0691.12123

PubMed Abstract | Crossref Full Text | Google Scholar

Lagerstrom, K. M., and Hadly, E. A. (1948). The under-investigated wild side of Escherichia coli: genetic diversity, pathogenicity and antimicrobial resistance in wild animals. Proc. R. Soc. B 288:20210399. doi: 10.1098/rspb.2021.0399

PubMed Abstract | Crossref Full Text | Google Scholar

Lagerstrom, K. M., and Hadly, E. A. (2023). Under-appreciated phylogroup diversity of Escherichia coli within and between animals at the urban-wildland interface. Appl. Environ. Microbiol. 89, e00142–e00123. doi: 10.1128/aem.00142-23

Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

López-Ramírez, K. L., Díaz-Maldonado, K. C., Vergara Espinoza, M. A., Santamaría-Veliz, O., Serquén-López, L. M., Canelo Olinda, B., et al. (2018). Patrón de clonalidad mediante ERIC-PCR y REP-PCR de Escherichia coli y Klebsiella pneumoniae productores de betalactamasas de espectro extendido, aisladas de pacientes con infección urinaria intrahospitalaria. Hospital Regional Lambayeque, Perú. Horizonte Médico 18, 11–18. doi: 10.24265/horizmed.2018.v18n2.03

Crossref Full Text | Google Scholar

Loscalzo, J., Fauci, A. S., Kasper, D. L., Hauser, S. L., Longo, D. L., and Jameson, J. L. (2022). Harrison's principles of internal medicine. London: McGraw Hill.

Google Scholar

Luna-Pineda, V. M., Ochoa, S. A., Cruz-Cordova, A., Cazares-Dominguez, V., Reyes-Grajeda, J. P., Flores-Oropeza, M. A., et al. (2018). Features of urinary Escherichia coli isolated from children with complicated and uncomplicated urinary tract infections in Mexico. PLoS One 13:e0204934. doi: 10.1371/journal.pone.0204934

PubMed Abstract | Crossref Full Text | Google Scholar

Magiorakos, A.-P., Srinivasan, A., Carey, R. B., Carmeli, Y., Falagas, M., Giske, C., et al. (2012). Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect. 18, 268–281. doi: 10.1111/j.1469-0691.2011.03570.x

PubMed Abstract | Crossref Full Text | Google Scholar

Mahmoud, A. T., Ibrahem, R. A., Salim, M. T., Gabr, A., and Halby, H. M. (2020). Prevalence of some virulence factors and genotyping of hospital-acquired uropathogenic Escherichia coli isolates recovered from cancer patients. J. Global Antimicrobial Resistance. 23, 211–216. doi: 10.1016/j.jgar.2020.08.003

PubMed Abstract | Crossref Full Text | Google Scholar

Malekzadegan, Y., Rastegar, E., Moradi, M., Heidari, H., and Sedigh, E.-S. H. (2019). Prevalence of quinolone-resistant uropathogenic Escherichia coli in a tertiary care hospital in South Iran. Infection Drug Resist. 12, 1683–1689. doi: 10.2147/IDR.S206966

PubMed Abstract | Crossref Full Text | Google Scholar

Marquez-Salazar, D. A., Delgadillo-Valles, R., Hernandez-Acevedo, G. N., Barrios-Villa, E., Muñiz-Salazar, R., Lopez-Valencia, G., et al. (2025). Coexistence of PMQR and ESBL genes among clinical Escherichia coli isolates from community-acquired UTI in Mexicali, on the US-Mexico border. Braz. J. Infect. Dis. 29:104554. doi: 10.1016/j.bjid.2025.104554

PubMed Abstract | Crossref Full Text | Google Scholar

McLellan, L. K., and Hunstad, D. A. (2016). Urinary tract infection: pathogenesis and outlook. Trends Mol. Med. 22, 946–957. doi: 10.1016/j.molmed.2016.09.003

PubMed Abstract | Crossref Full Text | Google Scholar

Medina, M., and Castillo-Pino, E. (2019). An introduction to the epidemiology and burden of urinary tract infections. Ther. Adv. Urol. 11:1756287219832172. doi: 10.1177/1756287219832172

PubMed Abstract | Crossref Full Text | Google Scholar

Mendoza, K., Ortiz, B., Rivera, L., Peña, T., Chirinos-Escobar, M., Enríquez, L., et al. (2023). Monitoring of microbial contamination of groundwater in the upper Choluteca River basin, Honduras. Water. 15:2116. doi: 10.3390/w15112116

Crossref Full Text | Google Scholar

Movahedi, M., Zarei, O., Hazhirkamal, M., Karami, P., Shokoohizadeh, L., and Taheri, M. (2021). Molecular typing of Escherichia coli strains isolated from urinary tract infection by ERIC-PCR. Gene Reports. 23:101058. doi: 10.1016/j.genrep.2021.101058

Crossref Full Text | Google Scholar

Öztürk, R., and Murt, A. (2020). Epidemiology of urological infections: a global burden. World. J. Urol. 38, 2669–2679. doi: 10.1007/s00345-019-03071-4

PubMed Abstract | Crossref Full Text | Google Scholar

Poey, M. E., de Los, S. E., Aznarez, D., García-Laviña, C. X., and Laviña, M. (2024). Genetics of resistance to trimethoprim in cotrimoxazole resistant uropathogenic Escherichia coli: integrons, transposons, and single gene cassettes. Front. Microbiol. 15:1395953. doi: 10.3389/fmicb.2024.1395953

PubMed Abstract | Crossref Full Text | Google Scholar

R Core Team (2025). R: A language and environment for statistical Computing_. Austria: R Foundation for Statistical Computing V.

Google Scholar

Ramakrishnan, V., Marialouis, X. A., Al-Ansari, M. M., Al-Humaid, L., Santhanam, A., and Obulisamy, P. K. (2022). Multilocus sequence typing and ERIC-PCR fingerprinting of virulent clinical isolates of uropathogenic multidrug resistant Escherichia coli. J. King Saud Univ. Sci. 34:101874. doi: 10.1016/j.jksus.2022.101874

Crossref Full Text | Google Scholar

Raphael, E., Argante, L., Cinconze, E., Nannizzi, S., Belmont, C., Mastrangelo, C. F., et al. (2024). Incidence and recurrence of urinary tract infections caused by Uropathogenic Escherichia coli: a retrospective cohort study. Research and reports. Urology 22, 253–264. doi: 10.2147/RRU.S470605

Crossref Full Text | Google Scholar

Regua-Mangia, A. H., Irino, K., da Silva, P. R., Pimentel Bezerra, R. M., Santos Périssé, A. R., and Teixeira, L. M. (2010). Molecular characterization of uropathogenic and diarrheagenic Escherichia coli pathotypes. J. Basic Microbiol. 50, S107–S115. doi: 10.1002/jobm.200900364

Crossref Full Text | Google Scholar

Royer, G., Clermont, O., Marin, J., Condamine, B., Dion, S., Blanquart, F., et al. (2023). Epistatic interactions between the high pathogenicity island and other iron uptake systems shape Escherichia coli extra-intestinal virulence. Nat. Commun. 14:3667. doi: 10.1038/s41467-023-39428-y

PubMed Abstract | Crossref Full Text | Google Scholar

Ruiz, J., Simon, K., Horcajada, J. P., Velasco, M., Barranco, M., Roig, G., et al. (2002). Differences in virulence factors among clinical isolates of Escherichia coli causing cystitis and pyelonephritis in women and prostatitis in men. J. Clin. Microbiol. 40, 4445–4449. doi: 10.1128/JCM.40.12.4445-4449.2002

PubMed Abstract | Crossref Full Text | Google Scholar

Salah, F. D., Soubeiga, S. T., Ouattara, A. K., Sadji, A. Y., Metuor-Dabire, A., Obiri-Yeboah, D., et al. (2019). Distribution of quinolone resistance gene (qnr) in ESBL-producing Escherichia coli and Klebsiella spp. in Lomé, Togo. Antimicrob. Resist. Infect. Control 8, 1–8. doi: 10.1186/s13756-019-0552-0

Crossref Full Text | Google Scholar

Sarshar, M., Behzadi, P., Ambrosi, C., Zagaglia, C., Palamara, A. T., and Scribano, D. (2020). FimH and anti-adhesive therapeutics: a disarming strategy against uropathogens. Antibiotics. 9:397. doi: 10.3390/antibiotics9070397

PubMed Abstract | Crossref Full Text | Google Scholar

Sohail, S., Wisal, A., Rana, S., Abid Khan, R. M., Ullah, A., Khan, F.-A., et al. (2025). Druggable proteins of Alistipes reveal promising antimicrobial targets against chronic intestinal inflammation. Front. Pharmacol. 16:1555062. doi: 10.3389/fphar.2025.1555062

PubMed Abstract | Crossref Full Text | Google Scholar

Sora, V. M., Meroni, G., Martino, P. A., Soggiu, A., Bonizzi, L., and Zecconi, A. (2021). Extraintestinal pathogenic Escherichia coli: virulence factors and antibiotic resistance. Pathogens. 10:1355. doi: 10.3390/pathogens10111355

PubMed Abstract | Crossref Full Text | Google Scholar

Tabasi, M., Karam, M. R. A., Habibi, M., Yekaninejad, M. S., and Bouzari, S. (2015). Phenotypic assays to determine virulence factors of uropathogenic Escherichia coli (UPEC) isolates and their correlation with antibiotic resistance pattern. Osong. Public Health Res. Perspect. 6, 261–268. doi: 10.1016/j.phrp.2015.08.002

PubMed Abstract | Crossref Full Text | Google Scholar

Tenaillon, O., Skurnik, D., Picard, B., and Denamur, E. (2010). The population genetics of commensal Escherichia coli. Nat. Rev. Microbiol. 8, 207–217. doi: 10.1038/nrmicro2298

PubMed Abstract | Crossref Full Text | Google Scholar

Tiwari, S., da Costa, M. P., Almeida, S., Hassan, S. S., Jamal, S. B., Oliveira, A., et al. (2014). C. pseudotuberculosis Phop confers virulence and may be targeted by natural compounds. Integr. Biol. 6, 1088–1099. doi: 10.1039/C4IB00140K

PubMed Abstract | Crossref Full Text | Google Scholar

Ucros, R. S. (2009). Guías de pediatría práctica basadas en la evidencia. Guías de pediatría práctica basadas en la evidencia, XVIII–651-XVIII. Bogotá: Editorial Médica Panamericana.

Google Scholar

Walker, D. I., McQuillan, J., Taiwo, M., Parks, R., Stenton, C. A., Morgan, H., et al. (2017). A highly specific Escherichia coli qPCR and its comparison with existing methods for environmental waters. Water Res. 126, 101–110. doi: 10.1016/j.watres.2017.08.032

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, M.-C., Fan, Y.-H., Zhang, Y.-Z., Bregente, C. J. B., Lin, W.-H., Chen, C.-A., et al. (2023). Characterization of uropathogenic Escherichia coli phylogroups associated with antimicrobial resistance, virulence factor distribution, and virulence-related phenotypes. Infect. Genet. Evol. 114:105493. doi: 10.1016/j.meegid.2023.105493

PubMed Abstract | Crossref Full Text | Google Scholar

Wesolek, J. L., Wu, J. Y., Smalley, C. M., Wang, L., and Campbell, M. J. (2022). Risk factors for trimethoprim and sulfamethoxazole-resistant Escherichia coli in ED patients with urinary tract infections. Am. J. Emerg. Med. 56, 178–182. doi: 10.1016/j.ajem.2022.03.052

PubMed Abstract | Crossref Full Text | Google Scholar

Whelan, S., Lucey, B., and Finn, K. (2023). Uropathogenic Escherichia coli (UPEC)-associated urinary tract infections: the molecular basis for challenges to effective treatment. Microorganisms. 11:2169. doi: 10.3390/microorganisms11092169

PubMed Abstract | Crossref Full Text | Google Scholar

Xia, L., Liu, Y., Xia, S., Kudinha, T., Xiao, S.-n., Zhong, N.-s., et al. (2017). Prevalence of ST1193 clone and IncI1/ST16 plasmid in E-coli isolates carrying Bla CTX-M-55 gene from urinary tract infections patients in China. Sci. Rep. 7:44866. doi: 10.1038/srep44866

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, X., Chen, H., Zheng, Y., Qu, S., Wang, H., and Yi, F. (2022). Disease burden and long-term trends of urinary tract infections: a worldwide report. Front. Public Health 10:888205. doi: 10.3389/fpubh.2022.888205

PubMed Abstract | Crossref Full Text | Google Scholar

Yousuf, M., Iuliani, I., Veetil, R. T., Seshasayee, A. S. N., Sclavi, B., and Cosentino, L. M. (2020). Early fate of exogenous promoters in E. coli. Nucleic Acids Res. 48, 2348–2356. doi: 10.1093/nar/gkz1196

PubMed Abstract | Crossref Full Text | Google Scholar

Yun, K. W., Kim, H. Y., Park, H. K., Kim, W., and Lim, I. S. (2014). Virulence factors of uropathogenic Escherichia coli of urinary tract infections and asymptomatic bacteriuria in children. J. Microbiol. Immunol. Infect. 47, 455–461. doi: 10.1016/j.jmii.2013.07.010

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, L., Chen, X., Zhu, X., Yang, W., Dong, L., Xu, X., et al. (2009). Prevalence of virulence factors and antimicrobial resistance of uropathogenic Escherichia coli in Jiangsu province (China). Urology 74, 702–707. doi: 10.1016/j.urology.2009.01.042

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Uropathogenic Escherichia coli (UPEC), UPEC virulence, UPEC resistance, Honduras, urinary infection

Citation: Ortiz B, Ballesteros-Monrreal MG, Galindo C, Rivera D, Rivera L, Pérez F, Maldonado V, Valle E, Estrada A, Mendez-Pfeiffer P, Enríquez L, Valencia D and Fontecha G (2025) Comprehensive molecular and phenotypic profiling of uropathogenic Escherichia coli in a Honduran healthcare setting: virulence, resistance and phylogeny. Front. Microbiol. 16:1656938. doi: 10.3389/fmicb.2025.1656938

Received: 06 July 2025; Accepted: 17 September 2025;
Published: 02 October 2025.

Edited by:

Marwan Osman, Yale University, United States

Reviewed by:

Hassan S. S., University of Karachi, Pakistan
Eman Jassim Mohammed, Al Mustansiriya University, Iraq

Copyright © 2025 Ortiz, Ballesteros-Monrreal, Galindo, Rivera, Rivera, Pérez, Maldonado, Valle, Estrada, Mendez-Pfeiffer, Enríquez, Valencia and Fontecha. 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: Gustavo Fontecha, Z3VzdGF2by5mb250ZWNoYUB1bmFoLmVkdS5obg==

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.