Longitudinal genomic surveillance of multidrug-resistant Escherichia coli carriage in critical care patients

Mira El Chaar, Yaralynn Khoury, Gavin M. Douglas, Samir El Kazzi, Tamima Jisr, Shatha Soussi, Georgi Merhi, Rima A. Moghnieh, B. Jesse Shapiro
doi: https://doi.org/10.1101/2023.08.12.23293895

Abstract

Background Colonization with multidrug resistant E. coli strains causes a substantial health burden in hospitalized patients. We performed a longitudinal genomics study to investigate the colonization of resistant E. coli strains in critically ill patients, and to identify evolutionary changes and strain replacement events within patients.

Methods Patients were admitted to the intensive care unit and haematology wards at a major hospital in Lebanon. Perianal swabs were collected from participants on admission and during hospitalization, which were screened for extended-spectrum beta-lactamases and carbapenem-resistant Enterobacterales. We performed whole-genome sequencing and analysis on E. coli strains isolated from patients at multiple time points.

Results The E. coli isolates were genetically diverse, with 13 sequence types (STs) identified among 21 isolates sequenced. Five patients were colonized by ST131 encoding CTX-M-27, a type of beta-lactamase gene not previously been observed in Lebanon. Among the eight patients whose resident E. coli strains were tracked over time, five harbored the same E. coli strain with relatively few mutations over the 5 to 10 days of hospitalization. The other three patients were colonized by different E. coli strains over time.

Conclusion Our study provides evidence of strain diversity within patients during their hospitalization. While strains varied in their antimicrobial resistance profiles, the number of resistance genes did not increase over time. We also show that ST131 encoding CTX-M-27, which appears to be emerging as a globally important multidrug resistant E. coli strain, is also prevalent among critical care patients and deserves further monitoring.

1. Introduction

Multidrug-resistant Gram-negative bacteria represent an important threat in hospital settings (1–3). During hospitalization, patients may be colonized with bacteria that develop resistance due to changes in gene expression, or that evolve resistance via point mutation or horizontal gene transfer (4–6). In hospitalized patients, resistant strains’ survival and replication may depend on the selection pressure exerted by antibiotics (6). When resistant bacteria colonize the gastrointestinal tract, this provides an opportunity for transfer of resistance genes among pathogens and gut microbiome commensal bacteria (7–9). Multidrug-resistant Escherichia coli are a burden on healthcare systems, and are often responsible for treatment failures in patients (10–12). However, the success of the globally prevalent Escherichia coli sequence type 131 (ST131) is not easily explained simply by its antibiotic resistance profile and likely involves diverse colonization and virulence factors (13, 14). Resistant bacteria disseminate among critically-ill patients in hospitals, causing life-threatening infections (7, 15). The pathogenicity of E. coli is multifaceted, including genetic and environment factors allowing E. coli to expand its range of infection beyond the intestine, which increases disease severity (16–19).

Mixed infections of distinct E. coli strains within a single patient have been observed previously, and these strains can also evolve within patients (20). Studies tracking E. coli diversity within patients over time are still relatively rare and have focused mainly on extraintestinal infections. For example, the transition of E. coli from urine to blood usually involves very few genetic changes, but occasionally involves colonization by genetically distant strains (21). In acute infections, extraintestinal E. coli can rapidly evolve hypermutator phenotypes, generating dozens to hundreds of mutations that could provide adaptation to new tissue types (22). Over time scales of a few years, up to 32% of samples U.S. military personnel were colonized by multiple distinct E. coli phylogroups, indicating co-existence of distinct strains or replacement of one strain by another in the gut (23). In a UK hospital study, over 25% of patients were colonized by multiple distinct strains (19). Whether such strain dynamics occur during intestinal colonization over shorter time scales, or within hospitalized patients, remains unclear. We hypothesize that the hospital environment provides a diverse pool of strains for co-colonization or strain replacements, and may also select for antibiotic resistance genes as critically ill patients are treated with broad-spectrum antibiotics.

In Lebanon, the prevalence of extended-spectrum beta-lactamase (ESBL) producing E. coli has increased over time, from 12% in 2005 to 25% in 2012 (24). E. coli was the most frequently isolated bacterium from bloodstream infections (45.6%), with 79.6% of the isolates producing ESBLs (25). Most studies to date have described the prevalence of resistant E. coli strains isolated from hospitals or community settings in Lebanon (26–33). Others described their prevalence in animals and the environment (22, 34–42). Here, we analysed whole genome sequence data using a longitudinal approach to describe the colonization and evolution of multidrug resistant E. coli in critically ill patients in Lebanon.

2. Methods

2.1. Study design

A prospective observational cohort study was conducted over a six-month period from June 2021 to December 2021 at Makassed General Hospital in Lebanon. This is a 200 bed hospital located in a heavily populated area of Beirut with a medically underserved population. It serves nearly 15,000 inpatients per year that are mostly from middle to low socioeconomic status groups.

2.2. Sampling, bacterial culture and data collection

A standard protocol for infection and prevention control to screen for any patient admitted to the intensive care unit and to the hematology oncology unit has previously been established at Makassed General Hospital. A perianal swab was collected from each patient and screened for extended-spectrum beta-lactamases and carbapenem-resistant Enterobacterales (ESBL and CRE, respectively). In total, 144 patients were admitted to the intensive care and oncology units; they were recruited for this study and followed during their hospitalization. Only patients that were hospitalized for more than three days were included. In total, 97 were dropped because of their early discharge, death, or transfer to another hospital ward, leaving 47 remaining patients. Perianal swabs were collected by healthcare workers at admission and during hospitalization. A minimum of five days was required between the first, second or third collection time. A standardized questionnaire for each patient was filled by the resident consultee. The questionnaire included information about age, gender, cause of hospitalization, duration of hospitalization, antibiotics exposure (past and during hospitalization), types of antibiotics used, and bacterial infection during hospitalization.

2.3. Ethical consideration

This study was approved by the institutional review board of Makassed General Hospital. The ethical committee reviewed and approved the study protocol. An informed consent was provided by the participants that have a decision-making capacity. Participants lacking capacity were enrolled following discussion with their relatives that are involved in their care.

2.4. Bacterial Culture

Anal swabs were cultured on three different culture media plates: MacConkey agar without antibiotics, MacConkey agar with ertapenem (0.5 mg/L) to isolate carbapenamase-producing bacteria, and MacConkey agar with ceftriaxone (4 μg/ml) for isolation of ESBL-producing bacteria. All plates were incubated at 37 °C for 24 h. Following bacterial growth, two colonies of the same color, morphology and shape were picked and pooled, then taxonomically classified using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF). Antimicrobial susceptibility using the disk diffusion method was determined for all E. coli isolates for the following antibiotics: amoxicillin (20 μg), amoxicillin/clavulanic acid (20/10 μg), cefepime (30 μg), ceftriaxone (30 μg), piperacillin-tazobactam (100/10 µg), ceftazidime-avibactam (20/10 μg), cephalothin (30 μg), ertapenem (10 μg), fosfomycin (200 μg), trimethoprim/sulfamethoxazole (1.25/23.75 μg), ciprofloxacin (5 μg), colistin (10μg), imipenem (10 μg), nitrofurantoin (300 μg), amikacin (30 μg), tetracyclin (30 μg), and gentamicin (10 µg).

2.5. Whole genome sequencing

Genomic DNA was extracted from the E. coli isolates using QIAamp DNA Mini kit (Qiagen) following the manufacturer’s guidelines. DNA concentrations were measured using Qubit (Thermo Fisher Scientific) with the Qubit dsDNA HS Assay Kit (Thermo Fisher Scientific). Sequencing libraries were made with Nextera XT DNA Library Preparation Kit (Illumina) according to manufacturer’s instructions. Illumina sequencing was performed on a Novaseq- 6000 producing paired-end 2 × 150 bp reads, with an average base quality (Phred score) of 36 and G + C content of 54.1%.

2.6. Genomic Analysis

Reads were processed with Trimmomatic v0.39 (43) to remove adaptor sequences and quality filter the reads; we specified the TruSeq3 adapters and the quality filtering parameters of LEADING:3, TRAILING:3, SLIDINGWINDOW:4:15, and MINLEN:30. The assembly was performed using SPAdes Genome Assembler v3.15.4 with default parameters (44). Assemblies were annotated with Prokka v1.14.5 (45) with the kingdom specified as bacteria. The resulting annotated genes were processed by Panaroo v1.3.0 to create the core genome alignment and pangenome (46). This was run with the ‘strict’ clean mode parameter and with the ‘clustal’ option specified for the core genome alignment. Plots were generated using Roary v3.12.0 including the pangenome frequency plot, a presence and absence matrix against a tree and a pie chart of the pangenome, breaking down the core, soft core, shell and cloud (47). The phylogenetic tree along with the pangenome of the isolates were visualized with Phandango v1.1.0 (48). A maximum-likelihood phylogeny was generated based on the core genome alignment from Panaroo input into RaxML (46, 49). The resulting tree was visualised with the Interactive Tree of Life (iTOL v6) (50). SRST2 (v0.2.0) was used to identify MLST types and serotype (51). Phylogroups were classified based on the ClermonTyping method (52). FastANI v1.33 was used to calculate the Average Nucleotide Identity (ANI) of orthologous genomic regions between pairs of genome assemblies (53). The presence of antimicrobial resistance genes, putative virulence factors, and plasmid replicons were studied using ABRicate with the ResFinder database, Virulence Factor database, and PlasmidFinder database (54). We used Snippy v4.6.0 to compare E. coli genomes from the same patient and identify single nucleotide variants (SNVs). The reference genome used was the E. coli strain identified at T0 (at hospital admission) or at T1 (the first collection time during hospitalization) (55). The major allele frequency distribution was also based on the identified variants, including SNVs and insertions/deletions.

2.7. Data Availability

All sequence data has been deposited in DDBJ/ENA/GenBank under BioProject PRJNA962847 and accession numbers SAMN35447827, SAMN35449847, SAMN35450267, SAMN35690035, SAMN35690036, SAMN35690056, SAMN35690090, SAMN35690097, SAMN35690110, SAMN35709899, SAMN35709920, SAMN35711041, SAMN35711058, SAMN35711154, SAMN35713710, SAMN35714296, SAMN35716244, SAMN35728356, SAMN35728381, SAMN35728629, SAMN35731198 and SAMN35731315.

3. Results

3.1. Study participants and microbial carriage

In total, 20 (43%) patients were eligible for this study, all of whom were colonized by antibiotic resistant isolates. Throughout their hospitalization, 14 (70%) patients were colonized with E. coli, one with Raoultella ornithinolytica and one with Klebsiella pneumoniae. Four patients were colonized by multiple bacterial species during their hospitalization. Here, we focus on the E. coli-colonized patients. The median age of patients infected with E. coli was 67 years (range = 17 - 88 years), eight were female (57%) and six (43%) were male. These 14 patients were admitted to different hospital wards: the intensive care unit (N=7), the haematology and oncology ward (N=3) and the isolation ward (N=4). The median duration of hospitalization for these patients was 29 days (range = 9 - 73 days). The cause of admission was either pneumonia, chemotherapy, bleeding, or sepsis (Table S1). In total, 11 out of 14 patients had an indwelling catheter during hospitalization and 13 patients were under corticosteroid therapy. The patients were exposed mostly to broad-spectrum antibiotics during their hospitalization, including: meropenem which was used in 50% of patients; colistin and piperacillin/tazobactam in 50%; amikacin, ceftriaxone and ceftazidime/avibactam in 50%; and levofloxacin in 29% of patients. Vancomycin was used in 43% of patients; trimethoprim-sulfamethoxazole, teicoplanin, gentamycin and clarithromycin were used in less than 20% of patients. The median time between first sample collection (T0) and second collection (T1) was 9 days (range = 5 - 28 days). During hospitalization, 7/14 (50%) of patients were documented to have infectious episodes (Table S1), including respiratory tract infections (N=5), followed by bacteremia (N=5), urinary tract (N=3) and skin and soft tissue (N=3).

3.2. Genomic typing and phylogeny of E. coli isolates

In total, 12 out of 14 patients infected with E. coli were selected for further analysis, of which eight patients had isolates from multiple time points, yielding 22 whole-genome sequences in total. Two patients (MGH 3A/B and MGH 12A) and two isolates (MGH 2A and MGH 6B) were dropped since their assembly sizes exceeded 8 Mbp, much larger than expected for an E. coli genome. Two of these assemblies also had median SNV major allele frequencies in the 0.93-0.97 range, whereas all others were very close to 1 (Figure S1). This suggests that some of the unexpectedly large assemblies could have been due to contamination, or due to pooling colonies of two distinct E. coli strains. The remaining 22 genomes considered for further analysis can be considered single strains.

These 22 E. coli genomes were genetically diverse, with 13 STs (ST10, ST44, ST69, ST131, ST224, ST227, ST 540, ST648, ST1286, ST1431, ST1491) identified (Figure 1). The most common was ST131 which was carried by 5/12 patients (41%).

Figure 1:
Figure 1: Distribution of antimicrobial resistance genes across the E. coli genomes recovered from patients.

Genetic determinants of resistance are grouped according to their corresponding antimicrobial classes, which are colour coded. Sequence type profiles, serotypes, phylogroup, and Fim type are indicated for each isolate. The time differences between T0/T1 or T1/T2 are indicated under Collection Timeline. *MGH 2B, MGH 5B, MGH 6A and MGH 14A strains were sequenced only at T0 or T1. The plasmid replicons identified in each isolate were included in the table. The patients are differentiated by numbers and each letter next to the number represent the different time points of collection: A (T0), B (T1) and C (T2)

We inferred the maximum-likelihood phylogeny of these isolates using a core-gene alignment. We included E. coli genomes previously sequenced from clinical samples collected in Lebanon (31, 32). Based on this phylogeny, we identified 5 phylogroups. In total, five isolates were phylogroup C and eight were phylogroup A. The ST131 strains (within phylogroup B2) were the most frequent (N=7). Phylogroups F and D were each represented by one isolate. In total, 11 serotypes were identified with the most common being O25:H4 and O9:H30 (Figure 2).

Figure 2:
Figure 2: Maximum likelihood phylogenetic tree of 54 E. coli isolates from clinical samples in Lebanon.

The tree includes genomes sequenced in the current study (indicated by asterisks). The phylogroup of the samples are shown in different colors.

Of the six patients sampled at two time points 5-11 days apart, five patients (MGH 1, MGH 8, MGH 9, MGH 10, MGH 11) had pairs of isolates that clustered very closely on the phylogeny (Figure 2) and had identical profiles of AMR genes (Figure 1). This is consistent with E. coli lineages persisting in patients between sampling points, although we cannot exclude transmission events of near-identical genomes among patients; for example, the two isolates from patient MGH 8 were very closely related to isolates from patients MGH 4, 6, and 14, suggesting possibly recent transmission of ST131 strains. In contrast, the two isolates from patient MGH 13 sampled seven days apart came from entirely different phylogroups (C and D; Figure 2) and had distinct AMR gene and plasmid profiles (Figure 1). This suggests either a mixed infection in patient MGH 13 (with phylogroup C sampled first and D sampled second), or a strain replacement that occurred between time points.

3.3. Prevalence and pangenome variation of ST 131 strains

Among the 12 patients with sequenced isolates, five harbored the ST131 type. Of these, four strains belonged to the O25:H4 serotype and fimH 30 variant and one strain to the O16:H5 serotype and fimH 41 variant. ST131 formed a distinct clade on the phylogenetic tree, classified as phylogroup B2 (Figure 2). All ST131 genomes encode blaCTX-M-27, blaEC-5, and the IncF plasmid, and 80% of them also encode blaOXA-181. They included type 1 Fimbriae, P fimbriae, and enterotoxin. All ST131 genomes were highly genetically similar, with average nucleotide identity values ranging from 99.5% to 99.9% (Figure S2). A pangenome analysis revealed 3308 core genes present in all isolates, 85 ‘soft core’ genes, 2520 ‘shell genes,’ and 2653 ‘cloud’ genes. The gene presence/absence patterns are shown alongside the phylogeny and strain typing information (Figure 3). Notably, ST131 strains with the O25:H4 serotype have a unique pattern of gene presence/absence, which is distinct from other ST131 strains and from other phylogroups.

Figure 3:
Figure 3: Gene presence/absence matrix from pangenome analysis of 22 E. coli isolates.

The Pangenome and maximum likelihood tree and distribution of accessory genes were visualized using Phandango. Blue and white represent the presence and absence of genes, respectively.

3.4. Antimicrobial susceptibility profiles and genetic determinants of AMR

Among the 14 patients carrying E. coli, 13 carried multidrug-resistant (MDR) strains. Only one patient (MGH 11) carried a strain that was resistant to cephalosporins, yet it was susceptible to all other antibiotics (Table S2).

Susceptibility testing was performed on 29 E. coli isolates collected from patients at different time points. Identical strains (defined as ANI of 99.9% or more) isolated from the same patient at different time points were omitted from the analysis (N=7). Among the remaining 22 E. coli strains, 55% (N=12/22) were ESBL, 5% (N=1/22) CRE, and 41% (N=9/22) ESBL/CRE. All patients enrolled in this study carried resistant strains upon admission. A total of 95% (N=21/22) were resistant or intermediate to cefotaxime, 86% to ceftriaxone (N=19/22), 54% (N=12/22) to ceftazidime, 63% (N=14/22) to ciprofloxacin, 50% (N=11/22) to ertapenem, 32% (N=7/22) to imipenem and meropenem, and 10% (N=2/22) to ceftazidime/avibactam (Table S2).

The genomes of these isolates were screened for known genetic determinants of AMR using the ABRicate pipeline. More than 46 different acquired AMR genes were detected; β-lactamases (N=19 genes) and aminoglycosidases (N=9 genes) were among the most common, and patients contained 2-16 resistance genes, with a median of 11 per patient (Figure 1). Among the 12 patients with non-identical E. coli strains, Class A, C and D beta-lactamases were mostly identified, including blaCTX-M-15 (N=5/12, 41%), blaCTX-M-27 (N=5/12, 41%), blaEC-5 (N=5/12, 41%), blaEC-15 (N=4/12, 33%), blaOXA-181(N=4/12, 33%). The beta-lactamases blaNDM-5 and blaOXA-244 were seen in one and two patients respectively. Several other AMR genes were also detected (Figure 2), with qnrS1 (encoding quinolone resistance) being quite prevalent (N=9/12, 75%).

3.5. Plasmid and virulence factor identification

The ABRicate pipeline was also used to identify the incompatibility groups for the recovered plasmid replicons (Figure S3). IncFI was the most common plasmid (73%; N= 16/22), followed by IncFII (63%; M= 14/22), IncX (41%, N= 9/22), IncI (23%; N= 5/22), IncB (27%; N=6/22) and IncH (14%; N= 3/22).

We screened the isolates for virulence genes which are often encoded by E. coli, and detected a wide variety (Figure S4). Some virulence genes were commonly detected, including ompA (100%) which is an outer membrane protein required for conjugation. Iron regulatory proteins were detected in 86% (N= 19/22) of isolates.

We recovered some of the known markers of intestinal and extraintestinal virulence in some isolates, such as the invasion and evasion factors (kpsM and kpsD [41%, N= 9/22]) and adherence factors (fim [86%, N= 19/22] and pap genes [41%, N= 9/22]) that are associated with intestinal and extraintestinal infections in humans. One patient (MGH 9) harboured an enteroaggregative E. coli strain with an anti-aggregation gene, and another patient (MGH 5) an enterotoxigenic E. coli with heat-labile (LT) and heat-stable (ST) enterotoxins genes. Other virulence genes detected in the study included chu, which codes for an outer membrane hemin receptor (41%, N= 9/22), and the enterotoxin gene, senB (32%, N= 7/22).

3.6. Strain dynamics within patients

A comparative genomic analysis was performed to compare the strains isolated from patients at different time points. Among the eight patients who were followed over time, five (MGH 1, MGH 8, MGH 9, MGH 10 and MGH 11) had very closely-related E. coli isolates at T0 and T1 (Figure 1, Figure 2). The duration of hospitalization for these patients ranged between 5 and 11 days, and genomes differed by an average of 16 SNVs (range 10-22 SNVs). In contrast, three patients (MGH 4, MGH 7 and MGH 13) were colonized by a highly divergent strain (different by > 1000 SNVs or <97% ANI; Figure S2) between time points, leading to marked changes in AMR gene profiles (Figure 1). In two patients, (MGH 4 and MGH 7), there was a putative strain replacement between T0 and T1, with the same strain persisting between T1 and T2, differing by a few SNVs (Figure 4). Comparing only closely-related pairs of genomes from the same patient yielded no evidence for a molecular clock (Spearman’s correlation between number of SNVs and number of days separating isolates, rho = 0.24, p = 0.60). The lack of clock signal could be due to a small sample size, short duration of sampling, or due to sampling only one isolate per patient per time point.

Figure 4:
Figure 4: Summary of E. coli genetic diversity and putative strain replacements within patients over time.

Each patient is shown on a separate row, with time indicated along the horizontal axis. E. coli genomes are indicated with rectangles or triangles, connected by arrows showing the number of days between samples and their average nucleotide identity (ANI). Genomes sharing a high ANI (99.9% or more) are shown with the same shape, and those with lower ANI (different strains) are shown in different shapes. The total number of genetic variations within each patient are indicated on the right: INS (insertion sequences), and complex variants (the combination of SNVs and multiple nucleotide polymorphism).

Of the 120 genetic variants (SNVs or complex variants of multiple mismatches nucleotides or insertions/deletions) observed between closely-related isolates from the same patient, half (N=60/120) were located in non-coding regions and 29% (N= 35/120) were synonymous changes. Most occurred in hypothetical proteins (N=18/35, 51%) (Table S3). Missense variants, or insertions were observed in 21% of the total variants detected (N=25/120). These were detected in hypothetical proteins (N= 12/25, 48%), transposases (N=5/25, 20%), a putative protein YjdJ (N=1/25, 4%), antigen 43 (N=2/25, 8%), cytoskeleton bundling-enhancing antitoxin (N=1/25, 4%), ompF outer membrane porin F (N=1/25, 4%) and D-alanine--D-alanine ligase A (N=1/25, 4%) (Table S3). We next tested the hypothesis that the frequency of known AMR genes or virulence factors (VFs) might increase within patients over time, due to exposure to antibiotics and other selective pressures in the critical care hospital environment. Although certain patients (e.g. MGH 13) acquired several AMR genes and VFs over time, there was no consistent trend across patients (Table S3, Figure 5). Changes were most pronounced in patients MGH 4, MGH 7 and MGH 13, who experienced a putative strain replacement between time points (Figure 4, Figure 5).

Figure 5:
Figure 5: Summary of genes acquired or lost in E. coli strains isolated from eight patients at two different points.

Gene gains (values above 0) or losses (below 0) over time (between T0 and T1) are indicated for AMR genes, virulence factors (VF), plasmids and insertion sequences (IS).

Discussion

The carriage rate of resistant E. coli isolates has been increasing in both healthcare and community settings. A recent systematic review has shown that at least one in five inpatients worldwide were carriers of ESBL and that the Eastern Mediterranean, which includes Lebanon, had the highest carriage rate (45.6%) (56). Such high prevalence was confirmed in our study while screening patients for resistant Gram-negative bacteria: 70% of patients were colonized with a multi-drug resistant E. coli isolate upon admission and during their hospitalization. All patients were critically ill and were given broad-spectrum antibiotics during their hospitalization. Despite the limited number of patients enrolled, we were able to identify that rectal colonization by ST131 subgroup fimH30-O25b, clade C1-M27, harboring blaCTX-M-27 was prevalent, sampled in 33% (N=4/12) of patients. One additional patient harbored a closely related isolate: ST131 subgroup fimH41-O25b, clade A-M27, harboring blaCTX-M-27. Comparative analyses of the five ST131 isolates showed that these strains share an average nucleotide identity of more than 99% and harbor identical virulence factors including type 1 fimbriae (FimH30) and the secreted autotransporter toxin, both of which are involved in uropathogenesis by mediating human bladder epithelial cell adhesion, invasion and biofilm formation (57–60). This strain is also resistant to cephalosporins, carbapenems, quinolones and tetracyclines. ST131 E. coli isolates differ from most other MDR E. coli by being more pathogenic causing often urinary tract infections (61).

The ST131 strains carrying CTX-M-15 have been reported to cause many infections globally. Recently, the the CTX-M-27-producing clade C1 of E. coli ST131 has emerged, and is thought to have epidemic potential both in community and healthcare setting (62–64). To the best of our knowledge, our study reports the first description of E. coli ST131 clade C1-M27 circulating in Lebanon. ST131 in Lebanon was first described in 2016 in animals; however, it was associated with blaCTX-M-15 (39). In our study, we have also shown the dominance of blaCTX-M-15 in different strains. This gene generally colonizes the gastrointestinal tract of farmed cattle and birds, as well as raw meat intended for human consumption which is a part of the Lebanese weekly diet, and was also identified in surface water (27, 39, 42). Further surveillance should be implemented in Lebanon to understand the transmission of blaCTX-M-15 and blaCTX-M-27 from animals to humans. We also identified blaOXA-181 in the ST131 clone. The blaOXA-181 gene is mainly found in E. coli and K. pneumoniae. It was first reported in 2007 (65) and was subsequently identified in several countries (66–69). However, it is not often isolated in ST131 but mostly linked with ST410 and ST1284 (70–72). The presence of both blaOXA-181 and blaCTX-M27 genes in ST131 was first observed in a young man with a war-related wound E. coli infection, which escalated to a series of recurrent infections over three months; this strain was shown to harbor a blaCTX-M 27 gene that was transferred from Morganella morganii (73). This exemplifies the ability of E. coli to rapidly acquire resistance genes from other species and highlights the need for continuous surveillance of gene transfer and resistance evolution in different E. coli sequence types.

To date, few studies have investigated within-patient diversity of E. coli in the same patient over time (74). In this pilot study, we sequenced E. coli genomes upon admission and during hospitalization. A limitation of our study is that we only sequenced two pooled isolates per patient per time point, making it difficult to distinguish strain replacement events from persistent mixed infections. Future efforts should ideally sequence multiple genomes from each sample to better characterize within-patient diversity and improve inference of transmission events (19, 60–62). In addition, using short read sequencing, it was difficult to confidently link resistance genes to plasmids; long-read sequencing could be used to fully assemble plasmids and make these linkages with confidence. Finally, our study only sampled perianal swabs. Sampling additional body sites and tissues could help establish E. coli transmission routes within patients.

Despite the limitations, we identified three patients with a possible strain replacement over a few days of hospitalization. In two of these patients sampled at a third time point, the new strain was retained. While we cannot exclude stable co-colonization of these different strains over time, which could have been missed by sequencing two pooled genomes per sample, our study opens the possibility that strain replacement events could plausibly occur during hospitalization. We also hypothesized that the frequency of AMR genes would increase over time in our patients, who were all treated with antibiotics. Although patients and strains both varied widely in their AMR and virulence gene content, there was no evidence for increasing frequency of AMR genes over time. Our work has therefore shown no major changes in SNV occurrence that accounts for disease susceptibility and resistance; this could be due to a lack of power to detect an effect in our small cohort, or could also suggest that resistant strains are already circulating and that resistance is transmitted upon infection rather than evolving within patients.

In conclusion, our results provide evidence for the recent emergence of ST131 subgroup fimH30-O25b, clade C1-M27, in Lebanon and reinforces the need for continuous genomic surveillance of this clone among patients in critical care units. Our study enables longitudinal stain tracking of E. coli strains, their colonization dynamics and their diversity within hospitals over time.

Data Availability

All data produced in the present work are contained in the manuscript

Supplementary Figures legends

Figure S1: Boxplots displaying the major allele frequency distributions (i.e., frequency of higher allele, whether alternative or reference) based on single-nucleotide variants and insertions/deletions. These variants were called by Snippy, based on mapped reads to a reference Escherichia coli genome. Samples are separated by whether they passed an independent quality control step of whether the whole-genome assembly was >= 8 megabase pairs (Mbp), indicating a potential mixture of strains or other contamination. The number of variants per sample ranges from 8,971-82,636 (mean: 45,920). The box within each boxplots represents the median (centre line) and the 25th and 75th percentiles as the box edges. The whiskers represent 1.5 * inter-quartile range (i.e., the range between the 25th and 75th percentiles), or to the most extreme value. Points indicate outliers outside this range.

Figure S2: Heatmap summary of Average Nucleotide Identity between each of the 22 strains sequenced in this study. Values range from 96% ANI to 100% ANI.

Figure S3: Heatmap showing the plasmid types identified in each of 22 E. coli genomes. The x axis represents the plasmid types.

Figure S4: Heatmap showing the virulence factors identified in each of 22 E. coli genomes. The x axis represents the virulent genes identified using ABRicate.

Supplementary Tables legends

Table S1: Demographic and Clinical characteristics of 14 patients enrolled in the study. The table includes demographic data in addition to the hospitalization status, risk factors, and antibiotic use during hospitalization. N/A refers to Not Available.

Table S2: Antimicrobial susceptibility testing results of 32 isolates identified in 15 patients using the disk diffusion method.

Table S3: General characteristics of genetic variants identified in seven patients at two different time points.

References

  1. 1.↵
    Jean SS, Harnod D, Hsueh PR. Global Threat of Carbapenem-Resistant Gram-Negative Bacteria. Front Cell Infect Microbiol. 2022;12:823684.
    OpenUrl
  2. 2.
    Haque M, Sartelli M, McKimm J, Abu Bakar M. Health care-associated infections - an overview. Infect Drug Resist. 2018;11:2321–33.
    OpenUrlCrossRefPubMed
  3. 3.↵
    European Antimicrobial Resistance C. The burden of bacterial antimicrobial resistance in the WHO European region in 2019: a cross-country systematic analysis. Lancet Public Health. 2022;7(11):e897–e913.
    OpenUrl
  4. 4.↵
    Woodford N, Ellington MJ. The emergence of antibiotic resistance by mutation. Clin Microbiol Infect. 2007;13(1):5–18.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.
    van Hoek AH, Mevius D, Guerra B, Mullany P, Roberts AP, Aarts HJ. Acquired antibiotic resistance genes: an overview. Front Microbiol. 2011;2:203.
    OpenUrlCrossRefPubMed
  6. 6.↵
    Munita JM, Arias CA. Mechanisms of Antibiotic Resistance. Microbiol Spectr. 2016;4(2).
  7. 7.↵
    Kantele A, Kuenzli E, Dunn SJ, Dance DAB, Newton PN, Davong V, et al. Dynamics of intestinal multidrug-resistant bacteria colonisation contracted by visitors to a high-endemic setting: a prospective, daily, real-time sampling study. Lancet Microbe. 2021;2(4):e151–e8.
    OpenUrl
  8. 8.
    Brito IL. Examining horizontal gene transfer in microbial communities. Nat Rev Microbiol. 2021;19(7):442–53.
    OpenUrlCrossRef
  9. 9.↵
    Lerner A, Matthias T, Aminov R. Potential Effects of Horizontal Gene Exchange in the Human Gut. Front Immunol. 2017;8:1630.
    OpenUrlCrossRef
  10. 10.↵
    Manges AR, Geum HM, Guo A, Edens TJ, Fibke CD, Pitout JDD. Global Extraintestinal Pathogenic Escherichia coli (ExPEC) Lineages. Clin Microbiol Rev. 2019;32(3).
  11. 11.
    Mills EG, Martin MJ, Luo TL, Ong AC, Maybank R, Corey BW, et al. A one-year genomic investigation of Escherichia coli epidemiology and nosocomial spread at a large US healthcare network. Genome Med. 2022;14(1):147.
    OpenUrl
  12. 12.↵
    Matsumoto T, Mikamo H, Ohge H, Yanagihara K, Weerdenburg E, Go O, et al. Distribution of extraintestinal pathogenic Escherichia coli O-serotypes and antibiotic resistance in blood isolates collected from patients in a surveillance study in Japan. J Infect Chemother. 2022;28(11):1445–51.
    OpenUrl
  13. 13.↵
    McNally A, Kallonen T, Connor C, Abudahab K, Aanensen DM, Horner C, et al. Diversification of Colonization Factors in a Multidrug-Resistant Escherichia coli Lineage Evolving under Negative Frequency-Dependent Selection. mBio. 2019;10(2).
  14. 14.↵
    Pitout JDD, Finn TJ. The evolutionary puzzle of Escherichia coli ST131. Infect Genet Evol. 2020;81:104265.
    OpenUrlCrossRef
  15. 15.↵
    Prado V, Hernandez-Tejero M, Mucke MM, Marco F, Gu W, Amoros A, et al. Rectal colonization by resistant bacteria increases the risk of infection by the colonizing strain in critically ill patients with cirrhosis. J Hepatol. 2022;76(5):1079–89.
    OpenUrl
  16. 16.↵
    Rozwadowski M, Gawel D. Molecular Factors and Mechanisms Driving Multidrug Resistance in Uropathogenic Escherichia coli-An Update. Genes (Basel). 2022;13(8).
  17. 17.
    Daga AP, Koga VL, Soncini JGM, de Matos CM, Perugini MRE, Pelisson M, et al. Escherichia coli Bloodstream Infections in Patients at a University Hospital: Virulence Factors and Clinical Characteristics. Front Cell Infect Microbiol. 2019;9:191.
    OpenUrl
  18. 18.
    Nagarjuna D, Mittal G, Dhanda RS, Verma PK, Gaind R, Yadav M. Faecal Escherichia coli isolates show potential to cause endogenous infection in patients admitted to the ICU in a tertiary care hospital. New Microbes New Infect. 2015;7:57–66.
    OpenUrl
  19. 19.↵
    Ludden C, Coll F, Gouliouris T, Restif O, Blane B, Blackwell GA, et al. Defining nosocomial transmission of Escherichia coli and antimicrobial resistance genes: a genomic surveillance study. Lancet Microbe. 2021;2(9):e472–e80.
    OpenUrl
  20. 20.↵
    Levert M, Zamfir O, Clermont O, Bouvet O, Lespinats S, Hipeaux MC, et al. Molecular and evolutionary bases of within-patient genotypic and phenotypic diversity in Escherichia coli extraintestinal infections. PLoS Pathog. 2010;6(9):e1001125.
    OpenUrlCrossRefPubMed
  21. 21.↵
    McNally A, Alhashash F, Collins M, Alqasim A, Paszckiewicz K, Weston V, Diggle M. Genomic analysis of extra-intestinal pathogenic Escherichia coli urosepsis. Clin Microbiol Infect. 2013;19(8):E328–34.
    OpenUrlCrossRefPubMed
  22. 22.↵
    Bridier-Nahmias A, Launay A, Bleibtreu A, Magnan M, Walewski V, Chatel J, et al. Escherichia coli Genomic Diversity within Extraintestinal Acute Infections Argues for Adaptive Evolution at Play. mSphere. 2021;6(1).
  23. 23.↵
    Mende K, Beckius ML, Zera WC, Yu X, Cheatle KA, Aggarwal D, et al. Phenotypic and genotypic changes over time and across facilities of serial colonizing and infecting Escherichia coli isolates recovered from injured service members. J Clin Microbiol. 2014;52(11):3869–77.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    Daoud Z, Salem Sokhn E, Masri K, Matar GM, Doron S. Escherichia coli Isolated from Urinary Tract Infections of Lebanese Patients between 2005 and 2012: Epidemiology and Profiles of Resistance. Front Med (Lausanne). 2015;2:26.
    OpenUrl
  25. 25.↵
    Haddad S, Jabbour JF, Hindy JR, Makki M, Sabbagh A, Nayfeh M, et al. Bacterial bloodstream infections and patterns of resistance in patients with haematological malignancies at a tertiary centre in Lebanon over 10 years. J Glob Antimicrob Resist. 2021;27:228–35.
    OpenUrl
  26. 26.↵
    Obeid A, Maliha P, Abdallah S, Akl E, Deeb M, El Moussawi H, et al. ESBL-producing Escherichia coli and Klebsiella pneumoniae in two major Lebanese hospitals: molecular epidemiology and correlation with consumption. J Infect Dev Ctries. 2018;12(2.1):16S.
    OpenUrl
  27. 27.↵
    Al-Mir H, Osman M, Drapeau A, Hamze M, Madec JY, Haenni M. Spread of ESC-, carbapenem- and colistin-resistant Escherichia coli clones and plasmids within and between food workers in Lebanon. J Antimicrob Chemother. 2021;76(12):3135–43.
    OpenUrl
  28. 28.
    Moussally M, Zahreddine N, Kazma J, Ahmadieh R, Kan SS, Kanafan ZA. Prevalence of antibiotic-resistant organisms among hospitalized patients at a tertiary care center in Lebanon, 2010-2018. J Infect Public Health. 2021;14(1):12–6.
    OpenUrl
  29. 29.
    Rosenthal VD, Belkebir S, Zand F, Afeef M, Tanzi VL, Al-Abdely HM, et al. Six-year multicenter study on short-term peripheral venous catheters-related bloodstream infection rates in 246 intensive units of 83 hospitals in 52 cities of 14 countries of Middle East: Bahrain, Egypt, Iran, Jordan, Kingdom of Saudi Arabia, Kuwait, Lebanon, Morocco, Pakistan, Palestine, Sudan, Tunisia, Turkey, and United Arab Emirates-International Nosocomial Infection Control Consortium (INICC) findings. J Infect Public Health. 2020;13(8):1134–41.
    OpenUrlCrossRef
  30. 30.
    Moghnieh RA, Kanafani ZA, Tabaja HZ, Sharara SL, Awad LS, Kanj SS. Epidemiology of common resistant bacterial pathogens in the countries of the Arab League. Lancet Infect Dis. 2018;18(12):e379–e94.
    OpenUrl
  31. 31.↵
    Tokajian S, Salloum T, Eisen JA, Jospin G, Farra A, Mokhbat JE, Coil DA. Genomic attributes of extended-spectrum beta-lactamase-producing Escherichia coli isolated from patients in Lebanon. Future Microbiol. 2017;12:213–26.
    OpenUrl
  32. 32.↵
    Dagher C, Salloum T, Alousi S, Arabaghian H, Araj GF, Tokajian S. Molecular characterization of Carbapenem resistant Escherichia coli recovered from a tertiary hospital in Lebanon. PLoS One. 2018;13(9):e0203323.
    OpenUrl
  33. 33.↵
    Khoury ME, Salloum T, Jisr T, Chaar ME, Tokajian S. Whole-Genome Sequence Analysis of Carbapenem-Resistant Enterobacteriaceae Recovered from Hospitalized Patients. J Glob Antimicrob Resist. 2023.
  34. 34.↵
    Tokajian S, Moghnieh R, Salloum T, Arabaghian H, Alousi S, Moussa J, et al. Extended-spectrum beta-lactamase-producing Escherichia coli in wastewaters and refugee camp in Lebanon. Future Microbiol. 2018;13:81–95.
    OpenUrl
  35. 35.
    Kassem, II, Mann D, Li S, Deng X. Draft genome sequences and resistome analysis of multidrug-resistant mcr-1-harbouring Escherichia coli isolated from pre-harvest poultry in Lebanon. J Glob Antimicrob Resist. 2021;25:114–6.
    OpenUrl
  36. 36.
    Dandachi I, Fayad E, Sleiman A, Daoud Z, Rolain JM. Dissemination of Multidrug-Resistant and mcr-1 Gram-Negative Bacilli in Broilers, Farm Workers, and the Surrounding Environment in Lebanon. Microb Drug Resist. 2020;26(4):368–77.
    OpenUrl
  37. 37.
    Dandachi I, Sokhn ES, Dahdouh EA, Azar E, El-Bazzal B, Rolain JM, Daoud Z. Prevalence and Characterization of Multi-Drug-Resistant Gram-Negative Bacilli Isolated From Lebanese Poultry: A Nationwide Study. Front Microbiol. 2018;9:550.
  38. 38.
    Diab M, Hamze M, Bonnet R, Saras E, Madec JY, Haenni M. Extended-spectrum beta-lactamase (ESBL)- and carbapenemase-producing Enterobacteriaceae in water sources in Lebanon. Vet Microbiol. 2018;217:97–103.
    OpenUrl
  39. 39.↵
    Diab M, Hamze M, Madec JY, Haenni M. High Prevalence of Non-ST131 CTX-M-15-Producing Escherichia coli in Healthy Cattle in Lebanon. Microb Drug Resist. 2017;23(2):261–6.
    OpenUrlCrossRef
  40. 40.
    Kassem, II, Osman M, Hassan J, Sulaiman A, Mann D, Esseili MA, et al. First Report of the Mobile Colistin Resistance Gene, mcr-1.26, in Multidrug-Resistant Escherichia coli Isolated from Retail Chicken Meat. J Glob Antimicrob Resist. 2023.
  41. 41.
    Mikhayel M, Leclercq SO, Sarkis DK, Doublet B. Occurrence of the Colistin Resistance Gene mcr-1 and Additional Antibiotic Resistance Genes in ESBL/AmpC-Producing Escherichia coli from Poultry in Lebanon: A Nationwide Survey. Microbiol Spectr. 2021;9(2):e0002521.
    OpenUrl
  42. 42.↵
    Moussa J, Abboud E, Tokajian S. The dissemination of antimicrobial resistance determinants in surface water sources in Lebanon. FEMS Microbiol Ecol. 2021;97(9).
  43. 43.↵
    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
    OpenUrlCrossRefPubMed
  45. 45.↵
    Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    Tonkin-Hill G, MacAlasdair N, Ruis C, Weimann A, Horesh G, Lees JA, et al. Producing polished prokaryotic pangenomes with the Panaroo pipeline. Genome Biol. 2020;21(1):180.
    OpenUrlCrossRefPubMed
  47. 47.↵
    Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MT, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31(22):3691–3.
    OpenUrlCrossRefPubMed
  48. 48.↵
    Hadfield J, Croucher NJ, Goater RJ, Abudahab K, Aanensen DM, Harris SR. Phandango: an interactive viewer for bacterial population genomics. Bioinformatics. 2018;34(2):292–3.
    OpenUrl
  49. 49.↵
    Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–W6.
    OpenUrlCrossRefPubMed
  51. 51.↵
    Inouye M, Dashnow H, Raven LA, Schultz MB, Pope BJ, Tomita T, et al. SRST2: Rapid genomic surveillance for public health and hospital microbiology labs. Genome Med. 2014;6(11):90.
    OpenUrlCrossRefPubMed
  52. 52.↵
    Beghain J, Bridier-Nahmias A, Le Nagard H, Denamur E, Clermont O. ClermonTyping: an easy-to-use and accurate in silico method for Escherichia genus strain phylotyping. Microb Genom. 2018;4(7).
  53. 53.↵
    Jain C, Rodriguez RL, Phillippy AM, Konstantinidis KT, Aluru S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun. 2018;9(1):5114.
    OpenUrlCrossRefPubMed
  54. 54.↵
    T S. ABRicate GitHub2020 [Available from: https://github.com/tseemann/abricate
  55. 55.↵
    T. S. Snippy: fast bacterial variant calling from NGS reads. 2020.
  56. 56.↵
    Bezabih YM, Bezabih A, Dion M, Batard E, Teka S, Obole A, et al. Comparison of the global prevalence and trend of human intestinal carriage of ESBL-producing Escherichia coli between healthcare and community settings: a systematic review and meta-analysis. JAC Antimicrob Resist. 2022;4(3):dlac048.
    OpenUrl
  57. 57.↵
    Dautin N. Serine protease autotransporters of enterobacteriaceae (SPATEs): biogenesis and function. Toxins (Basel). 2010;2(6):1179–206.
    OpenUrl
  58. 58.
    Qin J, Wilson KA, Sarkar S, Heras B, O’Mara ML, Totsika M. Conserved FimH mutations in the global Escherichia coli ST131 multi-drug resistant lineage weaken interdomain interactions and alter adhesin function. Comput Struct Biotechnol J. 2022;20:4532–41.
    OpenUrl
  59. 59.
    Spaulding CN, Klein RD, Ruer S, Kau AL, Schreiber HL, Cusumano ZT, et al. Selective depletion of uropathogenic E. coli from the gut by a FimH antagonist. Nature. 2017;546(7659):528–32.
    OpenUrlCrossRefPubMed
  60. 60.↵
    Poole NM, Green SI, Rajan A, Vela LE, Zeng XL, Estes MK, Maresso AW. Role for FimH in Extraintestinal Pathogenic Escherichia coli Invasion and Translocation through the Intestinal Epithelium. Infect Immun. 2017;85(11).
  61. 61.↵
    Nicolas-Chanoine MH, Bertrand X, Madec JY. Escherichia coli ST131, an intriguing clonal group. Clin Microbiol Rev. 2014;27(3):543–74.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    Peirano G, Pitout JDD. Extended-Spectrum beta-Lactamase-Producing Enterobacteriaceae: Update on Molecular Epidemiology and Treatment Options. Drugs. 2019;79(14):1529–41.
    OpenUrlCrossRef
  63. 63.
    Matsumura Y, Johnson JR, Yamamoto M, Nagao M, Tanaka M, Takakura S, et al. CTX-M-27- and CTX-M-14-producing, ciprofloxacin-resistant Escherichia coli of the H30 subclonal group within ST131 drive a Japanese regional ESBL epidemic. J Antimicrob Chemother. 2015;70(6):1639–49.
    OpenUrlCrossRefPubMed
  64. 64.↵
    Ghosh H, Doijad S, Falgenhauer L, Fritzenwanker M, Imirzalioglu C, Chakraborty T. bla(CTX-M-27)-Encoding Escherichia coli Sequence Type 131 Lineage C1-M27 Clone in Clinical Isolates, Germany. Emerg Infect Dis. 2017;23(10):1754–6.
    OpenUrlCrossRef
  65. 65.↵
    Castanheira M, Deshpande LM, Mathai D, Bell JM, Jones RN, Mendes RE. Early dissemination of NDM-1- and OXA-181-producing Enterobacteriaceae in Indian hospitals: report from the SENTRY Antimicrobial Surveillance Program, 2006-2007. Antimicrob Agents Chemother. 2011;55(3):1274–8.
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    Kim JS, Hong CK, Park SH, Jin YH, Han S, Kim HS, et al. Emergence of NDM-4 and OXA-181 carbapenemase-producing Klebsiella pneumoniae. J Glob Antimicrob Resist. 2020;20:332–3.
    OpenUrl
  67. 67.
    Nigg A, Brilhante M, Dazio V, Clement M, Collaud A, Gobeli Brawand S, et al. Shedding of OXA-181 carbapenemase-producing Escherichia coli from companion animals after hospitalisation in Switzerland: an outbreak in 2018. Euro Surveill. 2019;24(39).
  68. 68.
    Abd El Ghany M, Sharaf H, Al-Agamy MH, Shibl A, Hill-Cawthorne GA, Hong PY. Genomic characterization of NDM-1 and 5, and OXA-181 carbapenemases in uropathogenic Escherichia coli isolates from Riyadh, Saudi Arabia. PLoS One. 2018;13(8):e0201613.
    OpenUrl
  69. 69.↵
    Ahn K, Hwang GY, Kim YK, Kim HY, Jeong HR, Hong JS, Uh Y. Nosocomial Outbreak Caused by NDM-5 and OXA-181 Carbapenemase Co-producing Escherichia coli. Infect Chemother. 2019;51(2):177–82.
    OpenUrl
  70. 70.↵
    Peirano G, Chen L, Nobrega D, Finn TJ, Kreiswirth BN, DeVinney R, Pitout JDD. Genomic Epidemiology of Global Carbapenemase-Producing Escherichia coli, 2015-2017. Emerg Infect Dis. 2022;28(5):924–31.
    OpenUrl
  71. 71.
    Mataseje LF, Boyd DA, Fuller J, Haldane D, Hoang L, Lefebvre B, et al. Characterization of OXA-48- like carbapenemase producers in Canada, 2011-14. J Antimicrob Chemother. 2018;73(3):626–33.
    OpenUrl
  72. 72.↵
    Pitout JDD, Peirano G, Kock MM, Strydom KA, Matsumura Y. The Global Ascendency of OXA-48- Type Carbapenemases. Clin Microbiol Rev. 2019;33(1).
  73. 73.↵
    McGann P, Snesrud E, Ong AC, Appalla L, Koren M, Kwak YI, et al. War wound treatment complications due to transfer of an IncN plasmid harboring bla(OXA-181) from Morganella morganii to CTX-M-27-producing sequence type 131 Escherichia coli. Antimicrob Agents Chemother. 2015;59(6):3556–62.
    OpenUrlAbstract/FREE Full Text
  74. 74.↵
    Richter TKS, Hazen TH, Lam D, Coles CL, Seidman JC, You Y, et al. Temporal Variability of Escherichia coli Diversity in the Gastrointestinal Tracts of Tanzanian Children with and without Exposure to Antibiotics. mSphere. 2018;3(6).
Back to top
Posted August 16, 2023.
Download PDF

Supplementary Material

Data/Code
Email
Longitudinal genomic surveillance of multidrug-resistant Escherichia coli carriage in critical care patients
Mira El Chaar, Yaralynn Khoury, Gavin M. Douglas, Samir El Kazzi, Tamima Jisr, Shatha Soussi, Georgi Merhi, Rima A. Moghnieh, B. Jesse Shapiro
medRxiv 2023.08.12.23293895; doi: https://doi.org/10.1101/2023.08.12.23293895
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area