INTRODUCTION
Urinary tract infections (UTIs) are among the most prevalent bacterial diseases, posing a considerable burden to public health, particularly in the face of rising antimicrobial resistance (
1). Risk factors for UTIs include pregnancy, anatomical, and functional abnormalities of the urinary tract and the use of indwelling urinary catheters (
2). Catheters are the major risk factor for healthcare-acquired infections, and catheter-associated UTIs (CAUTIs) remain a prominent health challenge for patients in intensive care (
3).
Gram-negative pathogens have traditionally been considered the primary causative agents of UTIs (
4); however, recent studies reveal an increasing prevalence of
Enterococcus species in urinary infections (
5,
6). Historically, uropathogenic
Escherichia coli (UPEC) have been responsible for 75%–85% of all UTIs. Although Gram-positive bacteria constitute less than 15% of community-acquired UTIs, they are more frequently associated with nosocomial infections (
7,
8).
Enterococcus, naturally resistant to many first-line antibiotics and capable of forming biofilms, contributes to immune evasion and treatment failure (
9). Certain
Enterococcus strains can acquire high-level resistance to ?-lactams, aminoglycosides, glycopeptides, and even combined antibiotic therapies. In Scandinavia, vancomycin (VAN)-resistant
Enterococcus are rarely found but are causing outbreaks at an increasing and alarming speed (
10–12). Ampicillin (AMP)-resistant
E. faecium is gradually more common in healthcare-associated infections (
13). In the United States,
Enterococcus spp. were responsible for 12% of all CAUTIs between 2006 and 2007 (
14). This is likely due to their ability to elicit proinflammatory responses in the bladder and form biofilm (
15), which enables them to persist and cause chronic infections.
Bacterial populations are inherently heterogeneous, providing selective advantages during environmental changes and profoundly influencing clinical outcomes (
16). Recent studies have demonstrated that
E. faecalis exhibits heterogeneity in adhesion and biofilm formation. These key virulence properties may contribute to prolonged hospital stays and treatment failure (
17), but it remains unknown to what extent heterogeneity or polymicrobial colonization affects patient outcomes (
18).
In this study, we focus on bacterial urine colonizers among catheterized critically ill patients with severe acute respiratory syndrome coronavirus 2 infection. Early studies from similar cohorts have revealed that severely ill patients have a significantly higher risk of acquiring bacterial co-infections, especially with resistant strains, due to factors such as increased antibiotic use, utilization of urinary catheters, and the administration of immunosuppressive drugs (
19,
20). A 2022 study in Madrid investigated 87 COVID-19 patients, with 89.6% having acquired UTIs and 67.9% being related to CAUTIs.
Enterococcus was identified as the dominant genus, representing 47.4% of the isolates in the study cohort (
5). However, the role of heterogeneity in these conditions remains unknown.
Distinguishing between species and clinically relevant strains of
Enterococcus can be challenging, and their rapid emergence and importance demand a fast and accurate screening method (
21,
22). In this study, we utilized the PhenePlate-RF (PhP-RF) system, a method for strain screening that provides a biochemical fingerprint for multiple clinical
Enterococcus species, first introduced by Kühn et al. (
23). This method is fast, highly reproducible, and possesses strong discriminatory power. The PhP-RF system has been successfully used and validated for screening
Enterococcus in multiple environments, including the food chain and sewage water (
24). However, its application to clinical
Enterococcus isolates has been rare (
22). Our study employs the PhP-RF system to investigate the heterogeneity and dynamics of clinical
Enterococcus during urine colonization.
DISCUSSION
Differentiating species and clinically relevant strains of Enterococcus, often unobserved by conventional screening methods, is challenging and requires a fast and accurate screening scheme. In this study, we used the PhP-RF system, which provides a biochemical fingerprint for various clinical Enterococcus species, to investigate the heterogeneity of clinical Enterococcus during catheterized urine colonization.
Our study revealed considerable heterogeneity among
E. faecium isolates. Although underexplored, particularly in the context of urinary infections, it has previously been demonstrated that
E. faecium infections may involve multiple clones simultaneously (
26,
27). Most isolates in our study belonged to the well-characterized clonal complex 17 (CC17, formerly known as clade A1), encompassing strains from ST127 (HWP004) and ST80 (HWP143, HWP199, HWP030, HWP072, HWP083, and HWP095). CC17 is widely recognized as the most significant cluster of
E. faecium, often linked to enhanced biofilm formation, increased virulence, and AMR in healthcare settings (
28,
29). Interestingly, HWP143:4:C was identified as ST22, which is not typically associated with pathogenicity (
30), but has been proposed as an ancestral lineage of the CC17 cluster (
31). Additionally, HWP006 isolates were classified as ST192, a sequence type recently identified as an emerging clone primarily in hospital settings in Germany (
32).
Most isolates were classified into distinct PhP types. Patient HWP143 exhibited clear subpopulations due to the presence of multiple ST types. In other strains, particularly those from longitudinal samples, metabolic variations may suggest metabolic adaptations to the bladder/catheter environment.
Enterococci are generally thought to lack the necessary enzymes to metabolize urine-based carbon and nitrogen sources, such as creatine and urea (
33,
34). However, patients in ICUs with conditions such as diabetes or acute kidney injury (like those in our cohort) may excrete unaltered carbohydrates, including glucose, lactose, and L-arabinose (
35), and this could speculatively vary further in the microenvironments of a catheter. Mannitol is sometimes administered as a diuretic in clinical settings (
36); however, this was not the case in our cohort. Isolates from both HWP004 and HWP030 exhibited increased metabolic activity of L-arabinose over time. For patient HWP143, the majority of the
E. faecium population (81%) belonged to the non-pathogenic ST22, suggesting that the presence of CC17 most likely would have been overlooked by conventional screening methods (
6).
In patient HWP030, a metabolic shift was observed between collection days 5 and 6, despite genetic similarity across the isolates, with further alterations leading to a loss of CIP resistance in 93% of the population by day 8. Across
E. faecium, there were variations in AMR within patient isolates, specifically against TOB (HWP004 and HWP006), PIP/TZP (HWP004), and CIP (HWP030), as well as differences in AMR genes (HWP004 and HWP143). Although none of these antibiotics are first-line treatments for
E. faecium, this finding adds to the growing body of evidence of AMR variability among heterogeneous populations within an infection (
37). Indeed, even when the dominant populations were initially resistant, the resistant phenotype rapidly diminished in subsequent collection days, likely due to the absence of selective pressure, at least to a prevalence below our detection threshold.
Plasmid carriage varied among isolates from the same patient, independent of changes in AMR. On collection day 11, the dominant
E. faecium population of HWP004 (represented by HWP004:11E) carried rep1 but lacked repUS15, while a minor subpopulation (represented by HWP004:11L) carried both plasmids. By day 12, rep1 was absent in two of our strains (HWP004:12A/F) but reappeared in all strains by day 13, potentially maintained by isolate HWP004:12B. Upon closer investigation of HWP004:12B, it was observed that the two replicon-plasmids were located on what seemed to be a hybrid plasmid (
Fig. S4), but this was not investigated further. Previous studies have identified rep1 plasmids as conjugative plasmids specific to
E. faecalis, associated with AMR and virulence, and typically absent in
E. faecium (
Table S2) (
38,
39). Interestingly, all five HWP004:11
E. durans isolates carried this plasmid. This was not the only occurrence of the rep1 plasmid replicon, as it was additionally detected in HWP143:4C. Another notable finding was the presence of rep6 in HWP072:4A. Although this plasmid is small and cryptic, it has been proposed as an
E. faecalis-specific plasmid (
38). For isolates in patients HWP030, HWP083, and HWP199, significant nucleotide variations were observed in rep11a (suggested to be a toxin-related plasmid,
Table S2) and rep14a (small and cryptic), as these plasmid replicons showed deviations in identity level between collection days, indicating potential environmental adaptation through plasmid rearrangement.
E. faecalis isolates exhibited limited heterogeneity when assessing ST types within a patient but demonstrated notable diversity in both ST type and metabolic profile across different patients.
E. faecalis ST types are generally less homogeneous, likely due to their widespread ecological distribution, which contributes to the limited understanding of the significance of specific ST types (
40). The predominant ST type identified was ST179 (HWP003, HWP012, HWP099, and HWP085), which, along with ST16 (HWP028 and one isolate from HWP051), belongs to clonal complex 58 (
41). CC58 has been frequently associated with clinical infections, particularly in ICUs, and with AMR (
41,
42). Other identified ST types included ST538 (HWP004 and CC241), ST287 (one isolate from HWP051), ST26 (HWP078), ST81 (HWP116), and ST30 (HWP167), all ST types of which have previously been linked to clinical infections (
43–47). The stand-out example of heterogeneity among
E. faecalis was observed in patient HWP051, who was colonized by both ST16 (94% subpopulation, represented by HWP051:1A) and ST287 (6% subpopulation, represented by HWP051:1L) simultaneously. While prior research on
E. faecalis heterogeneity has often focused on strain properties, different patient samples, or outbreak scenarios (
17,
48,
49), the occurrence of mixed enterococcal site infections has not been extensively documented. The rarity of reporting mixed infections might be due to the challenges in distinguishing species, as well as the fact that screenings rarely encompass as many colonies as our current study for the investigation. Indeed, in our previous screen where 10 colonies were selected, and the batch was later assessed by MALDI-TOF, we did not detect
E. faecalis HWP004, even though
E. faecalis had been isolated from the blood of the same patient during our urine collection period (
6) (patient A in cited paper). This underscores that more pathogenic strains might be missed due to the lack of thorough investigations.
In contrast to
E. faecium,
E. faecalis isolates did not metabolize L-arabinose, melibiose, or raffinose, but efficiently utilized other carbon sources such as lactose, melezitose, inositol, sorbitol, mannitol, sucrose, amygdalin, and gluconate. These findings align with established knowledge of
E. faecalis metabolism (
50), reinforcing the species-specific metabolic pathways that have been previously documented. Notably, HWP003:5A lacked the ability to degrade melezitose, while HWP167:1A showed deficiencies in metabolizing inositol and gluconate. The isolates from patient HWP051 exhibited distinct metabolic profiles, with key differences in their ability to degrade melezitose, inositol, mannitol, amygdalin, and gluconate. Strains from patient HWP028 demonstrated increased metabolic activity for mannitol, sucrose, amygdalin, and gluconate over time, suggesting the potential importance of these carbon sources during
E. faecalis urine or catheter colonization (
6,
51,
52).
Most
E. faecalis isolates harbored the repUS43 plasmid, a conjugative plasmid previously identified in only one of our
E. faecium isolates. The isolate
E. faecalis HWP006:9L seemed to carry the repUS15 plasmid, suggesting possible HGT from
E. faecium. Interestingly, this plasmid was also found in one
E. durans isolate from the same isolation day. The repUS15 plasmid has recently been described as a carrier for high-level aminoglycoside resistance and as a potential carrier for
vanA-mediated VAN resistance (
53,
54). HWP078:1A possessed the highest number of plasmid replicons (six), including one (rep7a) not typically associated with enterococcal plasmids (
Table S2) (
38). Additionally,
E. faecalis HWP116:6A was the only isolate exhibiting phenotypic resistance to TOB, AMP, CIP, PIP, and TZP, despite carrying few AMR genes. HWP116:6A only carried the repUS43 plasmid replicon, which on further investigation turned out to be chromosomally integrated (
Fig. S4). This case of multidrug-resistant
E. faecalis is particularly concerning given the general characterization of
E. faecalis as more virulent but less resistant than
E. faecium (
55). Unlike
E. faecium,
E. faecalis plasmids did not exhibit nucleotide changes over time. We could not identify any longitudinal changes in AMR or AMR genes; the location of these genes remains unexplored in the current study. Two of the
E. faecalis strains were not found to carry any plasmid at all.
E. durans isolates, though less common, shared several metabolic features with
E. faecium, such as poor activity for melezitose and inositol metabolism and an inability to metabolize raffinose. Unlike
E. faecalis and many
E. faecium strains, none of the
E. durans isolates could metabolize sorbitol. Two
E. faecium strains were initially misidentified as
E. durans by MALDI-TOF but were confirmed as
E. faecium through sequencing. This misidentification could only partly be resolved by PhP-RF and is consistent with previous reports highlighting the similarity between
E. durans and
E. faecium (
56). Despite being prevalent for only one day, these findings are significant as
E. durans has been associated with higher rates of progression from urinary tract infections to bacteremia (
57,
58). Interestingly, speculated events of HGT might have occurred in patients co-colonized with
E. durans (HWP004 and HWP006). However, these patients also harbored
E. faecium and
E. faecalis with undetermined ST types, making it impossible to distinguish contributing species. While further molecular and genomic investigations are required to prove the phenomenon, our
E. durans HWP006:9A might be a case of interspecies plasmid conjugation.
The
E. durans strains from patient HWP004, isolated on the same day as
E. faecium HWP004:11E/L and
E. faecalis HWP004:11D, resembled the
E. faecium HWP004:11E/L phenotype, except for
E. durans HWP004:11A, which could not metabolize mannitol or sucrose. Isolates of
E. durans from patient HWP004 also exhibited a similar AMR phenotype to
E. faecium. In contrast,
E. durans HWP006 isolates displayed unique metabolic capabilities and carried more resistance genes than isolates from patient HWP004. However, no phenotypic AMR was observed, similar to co-isolated
E. faecalis isolates. Phylogenetically,
E. durans strains from patient HWP004 were metabolically similar, while the HWP006 isolates stood out both from each other and from the other
E. durans isolates, underscoring significant heterogeneity among
E. durans. Notably, all strains carried the rep1 plasmid, which has previously been suggested to be associated with
E. faecalis (
38), indicating a potential similarity between these two species.
We acknowledge several limitations in this study. Colony-forming unit (CFU) counts from frozen urine are not as optimal as using freshly collected urine. Despite our initial control experiments, the freezing procedure may have impacted different species of Enterococcus to varying degrees. This challenge is encountered by researchers working with biobanks and in microbiome research. Enterococcus spp. from complex clinical catheter samples are often polymicrobial. Although colonies were restreaked for purity, species were confirmed by MALDI-TOF, and basic local alignment search tool (BLAST) was performed on the same confirmed colony, we acknowledge that low-level contamination cannot be entirely excluded. The close phenotypic similarity in Enterococcus, combined with horizontal gene transfer, makes it inherently difficult to distinguish low-level contamination from genuine interspecies exchange. MLST classification was also inconclusive for two isolates due to incomplete or ambiguous allele calls. We therefore interpret such findings as hypothesis-generating rather than conclusive. We hope that the techniques employed in this study will help guide future research in developing effective methods for assessing bacteriuria.
Conclusion
This research provides significant insights into the heterogeneity and metabolic diversity of Enterococcus species during CAUTIs, particularly E. faecium and E. durans, across various dimensions. Our findings demonstrate considerable genetic and metabolic diversity among E. faecium isolates, with most belonging to CC17, a group known for its enhanced virulence and antibiotic resistance. Unlike earlier studies that focused primarily on single lineages, we identified the presence of multiple lineages and ST types within individual patients, including ST22, an ancestral lineage of CC17, and ST192, an emerging clone in hospital settings. Notably, metabolic adaptations were observed, such as increased metabolism of L-arabinose, and significant shifts in antibiotic resistance patterns over time. Our work reveals a previously underappreciated metabolic diversity.
E. faecium isolates exhibited variability in plasmid content and types. For instance, the rep1 plasmid, typically associated with E. faecalis, was found in both E. faecium and E. durans, hinting at HGT. Variations in plasmid nucleotide sequences over time were also observed, which could reflect environmental adaptation. E. faecalis isolates showed less genetic heterogeneity but displayed notable diversity in metabolic profiles across different patients. Evidence of mixed colonization with different ST types and variations in antibiotic resistance was found. Plasmids such as repUS43 and repUS15, plasmids not typically associated with E. faecalis, were identified, further highlighting the complexity of these infections. E. durans isolates, though less common, exhibited metabolic similarities to E. faecium. The misidentification of E. durans as E. faecium by MALDI-TOF, later confirmed through sequencing, underscores the need for accurate identification methods.
This research reveals that PhP-RF can be used to distinguish primarily clinical E. faecium and E. faecalis from each other, and that E. durans might metabolically group among E. faecium isolates. We have given insight into the dynamic nature of Enterococcus from urinary catheters, emphasizing significant variability in genetic and metabolic profiles both within individual patients and over time. The study offers a novel understanding of plasmid diversity among these species, which may have important implications for our understanding of antibiotic resistance and virulence. These findings underscore the critical need to consider metabolic adaptations and genetic heterogeneity, particularly in clinical settings where conventional screening methods may overlook the presence of multiple clones and lineages when studying Enterococcus infections, as they will come to impact disease severity and treatment outcome.