Enterococcus carry broad-host QACs tolerance genes widely spread across different bacterial taxa
The presence of QACs tolerance genes was investigated in Enterococcus from diverse epidemiological backgrounds, including 210 Enterococcus isolates from human and non-human sources (105 E. faecium and 105 E. faecalis) and 22,428 genomes publicly available at the NCBI database.
Only 1% of the isolates, which included one ST17
E. faecium and one ST6
E. faecalis carrying the
qacZ gene, were found to carry QACs tolerance genes. Additionally, among the public genomes analyzed, the presence of QACs tolerance genes (
qacA/B,
qacC,
qacG,
qacJ,
qacZ,
qrg,
bcrABC, and
oqxAB) was observed in only 0.5% of the genomes (
n = 117), which encompassed 93
E. faecalis, 23
E. faecium, and 1
Enterococcus lactis genomes (
Fig. 1; Table S1). The isolates ST17
E. faecium-E241 and ST6
E. faecalis-V583 were collected from a hospital sewage (2002) and a hospitalized patient (1987), respectively (
21,
30). The 117
Enterococcus isolates from public genomes were collected between 1987 and 2020 and belonged to diverse clonal lineages (30 STs among 93
E. faecalis; 13 STs among 23
E. faecium) and sources (Table S1). Among these, the
qac and
qrg genes were mostly detected in humans, including clinical and surveillance isolates (97%,
n = 69/71, among isolates with identifiable sources) (
P ≤ 0.01), whereas
bcrABC and
oqxAB were more prevalent in
Enterococcus strains from the food chain environment (68%,
n = 28/41) (
P ≤ 0.01) (
Fig. 1; Table S1).
QACs tolerance genes identical to the ones found in
Enterococcus strains were shared with bacterial species from Bacillota, Pseudomonadota, Actinomycetota, and Spirochaetota (former Firmicutes, Proteobacteria, Actinobacteria, and Spirochaetes, respectively) phyla (
Fig. 1) (
32). The
qacA/B,
qacC, and
qacJ were detected in 20 genera (20, 38, and 4 species, respectively), mainly in
Staphylococcus spp. (97%,
n = 2,691/2,777, among bacteria carrying
qac genes) (
P ≤ 0.01) and particularly in
S. aureus (58%,
n = 1,565/2,691 of the
Staphylococcus genomes carrying
qac genes) (
Fig. 1). The predominant genera sharing an identical
qrg with
Enterococcus strains was
Streptococcus sp. (98%,
n = 183/187) (
P ≤ 0.01), among the five genera (14 species) in which the gene was found. The
qacG and
qacZ were only found in
Enterococcus, according to our sequence selection criteria. As observed for
Enterococcus isolates, most strains carrying
qac and
qrg genes had a human origin (88%,
n = 2,387/2,699, among isolates with identifiable sources) (
P ≤ 0.01) (
Fig. 1). The
bcrABC gene cluster was identified in six genera (10 species), and
oqxAB gene cluster was identified in four genera (five species), mainly in
Listeria monocytogenes (97%,
n = 2,975/3,060, among those carrying
bcrABC) (
P ≤ 0.01) and
Escherichia coli (78%,
n = 80/102, among those carrying
oqxAB) (
P ≤ 0.01), respectively (
Fig. 1). Of note,
bcrABC and
oqxAB were, as mentioned above for enterococci, more prevalent in food chain isolates (84%,
n = 2,645/2,942) (
P ≤ 0.01) (
Fig. 1).
Genetic contexts of QACs tolerance genes are diverse and enriched in antimicrobial resistance genes
The co-occurrence of QACs tolerance genes with other antimicrobial resistance genes (e.g., antibiotics, metals) in the same genetic contexts and the similarity of these QACs tolerance genetic contexts between Enterococcus and other bacterial taxa were evaluated.
Considering the few complete genomes included in the analysis for the different taxa, namely of Enterococcus spp., most QACs tolerance genes were found on plasmids of different sizes (n = 18/19; 2–129 kb) rather than in the chromosome (n = 1/19) (Fig. S1). For the most part, the genetic contexts of QACs tolerance genes studied presented high variability among the different bacterial taxa (Fig. S1-A-G), with exception of some isolates sharing few genes. This was the case for E. faecium C132 and Staphylococcus aureus CC1-1 or Staphylococcus warneri Ani-LG-025, sharing the qacA/B and beta-lactam resistance genes, or E. faecalis C32 and Staphylococcus capitis 15–101 or Staphylococcus massiliensis P3, sharing the qacA/B and copper tolerance genes (Fig. S1-A), among others. Enterococcus spp. from the clinical and environmental settings shared qacZ genetic contexts carrying aac(6′)-Ie-aph(2′′)-Ia, coding for aminoglycoside resistance, insertion sequences, and recombinases (Fig. S1-D). Also, qrg genetic contexts were shared between E. faecalis and Streptococcus spp., mostly presenting insertion sequences, recombinases, and hypothetical proteins (Fig. S1-E). Finally, E. faecalis from different epidemiological backgrounds shared bcrABC genes, insertion sequences, genes coding for recombinases, hypothetical proteins, or bacteriocin associated (Fig. S1-F).
Regardless of the strains’ source, geographical region, or date of isolation, we observed that different genetic determinants for antibiotic resistance and metal tolerance were located adjacent to QACs tolerance genes. Of note, several of the genetic contexts analyzed (n = 32) in diverse bacteria, including Enterococcus spp., harbored genes conferring resistance to aminoglycosides, beta-lactams, macrolides, lincosamides, streptogramin B, mupirocin, bleomycin, tetracycline, chloramphenicol/florfenicol, trimethoprim, fosfomycin, and/or sulfonamide (Fig. S1-A,B,C,D,F,G). Metal tolerance genes, namely to copper, cadmium, zinc, or arsenic, were also detected within the vicinity of qacA/B, qacC, qacJ, qrg, and bcrABC genes (Fig. S1-A,B,C,E,F). Additionally, the genetic contexts analyzed were highly enriched in insertion sequences, recombinases, and replication associated proteins.
QACs susceptibility assays of E. faecalis and E. faecium
The susceptibility of 105 E. faecalis and 105 E. faecium isolates (including two isolates carrying qacZ), with diverse epidemiological and clonal backgrounds (Table S2), to BC and DDAC was determined.
The MIC
BC in standard conditions were similar for the
E. faecalis and
E. faecium studied (MIC
50 = 2 mg/L and MIC
90 = 2 mg/L for both) (
Fig. 2). The highest MIC
BC of 4 mg/L was observed in 1
E. faecalis and 9
E. faecium recovered from different sources, years, and clonal lineages. Of these, most were MDR (
n = 8/10; resistant to three or more antibiotics from different families), and two of them harbored the gene
qacZ (ST17
E. faecium-E241 and ST6
E. faecalis-V583). MBC
BC distributions for both species also showed similar MBC
50 = 2 mg/L and MBC
90 = 4 mg/L (
Fig. 2). Likewise, the 43% (
n = 45/105) of
E. faecalis and 19% (
n = 20/105) of
E. faecium isolates showing the highest MBC
BC of 4 mg/L comprised, in both cases, isolates from different epidemiological and genetic backgrounds. On the contrary,
E. faecium (
n = 8) with the lowest MBC
BC of 1 mg/L were mostly from the food chain (
n = 7/8), belonging to different STs and years. The value of MIC
BC and MBC
BC of the control strain
E. faecalis ATCC 29212 varied between 1-2 mg/L and 2-4 mg/L, respectively.
Both species were susceptible to lower concentrations of DDAC compared to BC (
P ≤ 0.01) (
Fig. 2). For DDAC, the MIC
50 was 1 mg/L for both species, and the MIC
90 was 1 mg/L or 2 mg/L for
E. faecalis and
E. faecium, respectively (
Fig. 2). Strains with an MIC
DDAC of 2 mg/L (
n = 8
E.
faecalis and
n = 29
E. faecium) were diverse and comprised most of the
Enterococcus with the highest MIC
BC (
n = 7/10, including the two isolates carrying
qacZ). Finally, the MBC
50 was 2 mg/L for both species, and the MBC
90 was 4 mg/L and 2 mg/L for
E. faecalis and
E. faecium, respectively. Again, isolates showing the highest MBC
DDAC of 4 mg/L or the lowest MBC
DDAC of 1 mg/L (18%,
n = 19/105, and 4%,
n = 4/105, of the
E. faecalis; 4%,
n = 4/105, and 23%,
n = 24/105, of the
E. faecium populations tested, respectively) were associated with different sources, years, and STs. The MIC
DDAC and MBC
DDAC of the control strain
E. faecalis ATCC 29212 varied between 0.5-1 mg/L and 1-2 mg/L, respectively.
Considering all sources,
E. faecalis seem to be more tolerant to bactericidal concentrations of the QACs tested (higher MBC) than
E. faecium (
P ≤ 0.01), although the latter presented a significantly higher MIC
DDAC (
P ≤ 0.01) (
Fig. 2). The MIC and MBC distributions of the isolates tested were analyzed separately by source and time span (5-year intervals) (Fig. S2), with the following significant differences among them: the MIC
DDAC and MIC
BC were higher for
E. faecium recovered from human infections compared to those from the food chain (
P ≤ 0.01), but a significant increasing trend in the MIC
BC over the years was detected in
E. faecium isolates from the food chain (
P ≤ 0.01) (Fig. S2). QACs susceptibility among MDR and non-MDR
E. faecalis and
E. faecium was similar (
P > 0.01) (Fig. S3). BC and DDAC MIC and MBC values for vancomycin- or linezolid-resistant isolates varied within the ranges described for the whole population (Fig. S3).
Different growth conditions occurring in the environment affect Enterococcus spp. susceptibility to QACs
To assess the activity of QACs on
Enterococcus growth or survival under diverse environmental conditions that mimic real contexts, MIC
BC and MBC
BC were determined in anaerobic conditions (e.g., occurring in sewage) (
37,
38), room temperature (22°C; representing abiotic surfaces in healthcare or agri-food sectors) (
39), and mildly acidic pH (pH = 5; resembling human skin and fresh fruits/vegetables contaminated with QACs residues) (
40 – 42). BC was the QAC chosen for the susceptibility assays with modified conditions, as both species demonstrated higher tolerance and phenotype variability (wider MIC and MBC ranges) to this compound compared to DDAC under standard conditions. MIC
BC and MBC
BC distributions for the modified conditions studied are shown in
Fig. 3.
The analysis revealed similar results between aerobic (standard) and anaerobic conditions (
P > 0.01), but increased tolerance to BC at room temperature (22°C) and/or mildly acidic pH (pH = 5) conditions for
E. faecalis and
E. faecium strains (
P ≤ 0.01) (
Fig. 3). In anaerobic conditions, MIC
BC or MBC
BC increased by no more than twofold in a few strains (six
E. faecalis and six
E. faecium) belonging to diverse sources, dates, and STs. The MBC
BC of
qacZ + ST17
E.
faecium E241 increased from 4 to 8 mg/L, outside the range obtained in standard conditions. Also, the MIC
BC and MBC
BC increased from 4 to 8 mg/L in a human
E. faecium strain (ST412), recovered from a clinical infection in 2011, that did not contain any of the known QACs tolerance genes.
Contrastingly, four to eightfold MIC
BC and MBC
BC increases were observed at 22°C and/or pH = 5 conditions in diverse strains of both species. When isolates were tested at pH = 5, the MIC
BC increased significantly for
E. faecalis (
P ≤ 0.01) but not for
E. faecium (
P > 0.01), although MBC
BC were higher for both (
P ≤ 0.01) (
Fig. 3). The MIC
BC of one
E. faecalis increased fourfold (from 1 to 4 mg/L; urban wastewater treatment plant from Tunisia, ST23, 2014) and the MIC
BC of the
qacZ + E. faecalis V583 increased from 4 to 8 mg/L. The highest MBC
BC obtained among the investigated
E. faecalis strains was also 8 mg/L, observed in nine strains exhibiting a twofold (
n = 8; from several sources and STs) or fourfold rise (
n = 1; poultry carcass, ST843, 2018). In the case of
E. faecium isolates, the highest MBC
BC at pH = 5 was 16 mg/L, corresponding to
qacZ + E. faecium E241 and one strain isolated from the faeces of a long-term care facility patient in 2016 (ST262). It is noteworthy that these strains exhibited an MBC
BC of 4 mg/L under standard conditions.
An increased tolerance to BC was more pronounced at 22°C than at pH = 5 for E. faecium, reflected both in MICBC (n = 3 isolates with fourfold and n = 2 isolates with eightfold increases at 22°C, compared to standard conditions) and MBCBC (n = 6 isolates with fourfold and n = 2 isolates with eightfold increases at 22°C) values. The highest MBCBC of 16 mg/L was detected in two food chain isolates (ready-to-eat salad and poultry carcass, ST12 and ST352, 2010–2019), while the remaining strains had an MBCBC of 8 mg/L (n = 18). Among E. faecalis, in all but one of the isolates, the MBCBC increased to 8 mg/L [following twofold (n = 17) or fourfold (n = 4) MBCBC increases] at 22°C, whereas the MICBC did not change (P > 0.01) for most strains (n = 19).
The effect of the combination of the test conditions (growth at 22°C and pH = 5) was also studied, excluding anaerobiosis since it did not significantly influence QACs susceptibility. Indeed, higher MICBC or MBCBC were observed for both species under 22°C and pH = 5 stress than when each condition was tested separately (P ≤ 0.01). MBCBC increased from 1 to 4 mg/L (standard conditions) to 8 mg/L at 22°C (93%; n = 21/22 E. faecalis and n = 18/20 E. faecium) and 16 mg/L (100%; n = 7/7 E. faecalis and n = 7/7 E. faecium) for the combination of 22°C and pH = 5 for most isolates, regardless of the species, epidemiological or genetic background. MICBC also reached the highest values of 16 mg/L in four E. faecium strains from human infection and colonization, hospital sewage (E. faecium ST17 carrying qacZ), and food chain (2002–2020), and of 8 mg/L in four E. faecalis from human infection (n = 2, including an ST6 strain carrying qacZ), food chain, and environmental origins (1987–2019).
Two
E. faecalis (V583 with
qacZ; S37-25 without
qacZ) and two
E. faecium (E241 with
qacZ; F1651 without
qacZ) were included in kinetic assays, showing changes in their growth dynamics under the modified conditions studied compared with the standard ones (Fig. S4 and S5). Such modifications in the growth curves were strain specific, with
E. faecalis isolates showing more similarities among each other compared to the
E. faecium isolates. This was observed in the kinetic assays conducted both with subinhibitory concentrations of BC (0.5 mg/L for
E. faecium F1651 and 1 mg/L for
E. faecium E241,
E. faecalis V583, and S37-25) as well as without BC (
Fig. 4; Fig. S4 and S5). Comparing the growth kinetics of each strain under BC plus modified conditions with the standard ones (
Fig. 4), different bacterial adaptations were observed, including extended lag phases and, in several cases, slower exponential growth, resulting in a delayed entry into the stationary phase. Also, for the
E. faecium F1651 and the two
E. faecalis, their growth with BC plus 22°C or 22°C + pH = 5 surpassed that occurring in standard conditions during the time of the assay, which may suggest that these modified conditions are better for bacterial multiplication under BC stress.
Specifically, at pH = 5 or 22°C, we described earlier that both
E. faecium and
E. faecalis increased their MBC
BC (two to eightfold) compared to standard conditions, which may be explained by the longer exponential growth phase in such modified conditions (
Fig. 4). Additionally, the highest MBC
BC increase observed for
E. faecium (eightfold) may be related to their exponential growth for even longer periods compared to
E. faecalis, which enter much earlier into the stationary phase (<20 hours) (
Fig. 4). Of note, for the combination of 22°C + pH = 5 conditions, growth kinetics for all
E. faecium and
E. faecalis strains included in the analysis have shown that they were still in the exponential multiplication phase after 24 hours. This may explain the higher MBC
BC observed for both species under the combination of 22°C and pH = 5 stress than for any other growth condition tested, with all isolates reaching the highest MBC
BC of 16 mg/L. No clear differences were detected in isolates carrying the
qacZ gene (
E. faecium E241;
E. faecalis V583) or not (
E. faecium F1651;
E. faecalis S37-25) in any condition tested.