YajC, a predicted membrane protein, promotes Enterococcus faecium biofilm formation in vitro and in a rat endocarditis model.

Biofilm formation is a critical step in the pathogenesis of difficult-to-treat Gram-positive bacterial infections. We identified that YajC, a conserved membrane protein in bacteria, plays a role in biofilm formation of the clinically relevant Enterococcus faecium strain E1162. Deletion of yajC conferred significantly impaired biofilm formation in vitro and was attenuated in a rat endocarditis model. Mass spectrometry analysis of supernatants of washed ΔyajC cells revealed increased amounts in cytoplasmic and cell-surface-located proteins, including biofilm-associated proteins, suggesting that proteins on the surface of the yajC mutant are only loosely attached. In Streptococcus mutans YajC has been identified in complex with proteins of two cotranslational membrane protein-insertion pathways; the signal recognition particle (SRP)-SecYEG-YajC-YidC1 and the SRP-YajC-YidC2 pathway, but its function is unknown. In S. mutans mutation of yidC1 and yidC2 resulted in impaired protein insertion in the cell membrane and secretion in the supernatant. The E. faecium genome contains all homologous genes encoding for the cotranslational membrane protein-insertion pathways. By combining the studies in S. mutans and E. faecium, we propose that YajC is involved in the stabilization of the SRP-SecYEG-YajC-YidC1 and SRP-YajC-Yid2 pathway or plays a role in retaining proteins for proper docking to the YidC insertases for translocation in and over the membrane.

Colonization of vancomycin-resistant Enterococcus faecium in human-derived colonic epithelium: unraveling the transcriptional dynamics of host-enterococcal interactions.

Enterococcus faecium is an opportunistic pathogen able to colonize the intestines of hospitalized patients. This initial colonization is an important step in the downstream pathogenesis, which includes outgrowth of the intestinal microbiota and potential infection of the host. The impact of intestinal overgrowth on host-enterococcal interactions is not well understood. We therefore applied a RNAseq approach in order to unravel the transcriptional dynamics of E. faecium upon co-culturing with human derived colonic epithelium. Co-cultures of colonic epithelium with a hospital-associated vancomycin resistant (vanA-type) E. faecium (VRE) showed that VRE resided on top of the colonic epithelium when analyzed by microscopy. RNAseq revealed that exposure to the colonic epithelium resulted in upregulation of 238 VRE genes compared to the control condition, including genes implicated in pili expression, conjugation (plasmid_2), genes related to sugar uptake, and biofilm formation (chromosome). In total, 260 were downregulated, including the vanA operon located on plasmid_3. Pathway analysis revealed an overall switch in metabolism to amino acid scavenging and reduction. In summary, our study demonstrates that co-culturing of VRE with human colonic epithelium promotes an elaborate gene response in VRE, enhancing our insight in host-E. faecium interactions, which might facilitate the design of novel anti-infectivity strategies.

PlasmidEC and gplas2: an optimized short-read approach to predict and reconstruct antibiotic resistance plasmids in Escherichia coli.

Accurate reconstruction of Escherichia coli antibiotic resistance gene (ARG) plasmids from Illumina sequencing data has proven to be a challenge with current bioinformatic tools. In this work, we present an improved method to reconstruct E. coli plasmids using short reads. We developed plasmidEC, an ensemble classifier that identifies plasmid-derived contigs by combining the output of three different binary classification tools. We showed that plasmidEC is especially suited to classify contigs derived from ARG plasmids with a high recall of 0.941. Additionally, we optimized gplas, a graph-based tool that bins plasmid-predicted contigs into distinct plasmid predictions. Gplas2 is more effective at recovering plasmids with large sequencing coverage variations and can be combined with the output of any binary classifier. The combination of plasmidEC with gplas2 showed a high completeness (median=0.818) and F1-Score (median=0.812) when reconstructing ARG plasmids and exceeded the binning capacity of the reference-based method MOB-suite. In the absence of long-read data, our method offers an excellent alternative to reconstruct ARG plasmids in E. coli.

Metagenomic assembly is the main bottleneck in the identification of mobile genetic elements.

Antimicrobial resistance genes (ARG) are commonly found on acquired mobile genetic elements (MGEs) such as plasmids or transposons. Understanding the spread of resistance genes associated with mobile elements (mARGs) across different hosts and environments requires linking ARGs to the existing mobile reservoir within bacterial communities. However, reconstructing mARGs in metagenomic data from diverse ecosystems poses computational challenges, including genome fragment reconstruction (assembly), high-throughput annotation of MGEs, and identification of their association with ARGs. Recently, several bioinformatics tools have been developed to identify assembled fragments of plasmids, phages, and insertion sequence (IS) elements in metagenomic data. These methods can help in understanding the dissemination of mARGs. To streamline the process of identifying mARGs in multiple samples, we combined these tools in an automated high-throughput open-source pipeline, MetaMobilePicker, that identifies ARGs associated with plasmids, IS elements and phages, starting from short metagenomic sequencing reads. This pipeline was used to identify these three elements on a simplified simulated metagenome dataset, comprising whole genome sequences from seven clinically relevant bacterial species containing 55 ARGs, nine plasmids and five phages. The results demonstrated moderate precision for the identification of plasmids (0.57) and phages (0.71), and moderate sensitivity of identification of IS elements (0.58) and ARGs (0.70). In this study, we aim to assess the main causes of this moderate performance of the MGE prediction tools in a comprehensive manner. We conducted a systematic benchmark, considering metagenomic read coverage, contig length cutoffs and investigating the performance of the classification algorithms. Our analysis revealed that the metagenomic assembly process is the primary bottleneck when linking ARGs to identified MGEs in short-read metagenomics sequencing experiments rather than ARGs and MGEs identification by the different tools.

Gut microbiome dynamics in index patients colonized with extended-spectrum beta-lactamase (ESBL)-producing Enterobacterales after hospital discharge and their household contacts.

IMPORTANCE: Colonization with extended-spectrum beta-lactamase-producing Enterobacterales (ESBL-PE) often precedes infections and is therefore considered as a great threat for public health. Here, we studied the gut microbiome dynamics in eight index patients colonized with ESBL-PE after hospital discharge and the impact of exposure to this index patient on the gut microbiome dynamics of their household contacts. We showed that the microbiome composition from index patients is different from their household contacts upon hospital discharge and that, in some of the index patients, their microbiome composition over time shifted toward the composition of their household contacts. In contrast, household contacts showed a stable microbiome composition over time irrespective of low-level extended-spectrum beta-lactamase-producing Escherichia coli (ESBL-Ec) or extended-spectrum beta-lactamase-producing Klebsiella pneumoniae (ESBL-Kp) gut colonization, suggesting that, in healthy microbiomes, colonization resistance is able to prevent ESBL-PE expansion.

The population-level impact of Enterococcus faecalis genetics on intestinal colonization and extraintestinal infection.

IMPORTANCE: Enterococcus faecalis causes life-threatening invasive hospital- and community-associated infections that are usually associated with multidrug resistance globally. Although E. faecalis infections cause opportunistic infections typically associated with antibiotic use, immunocompromised immune status, and other factors, they also possess an arsenal of virulence factors crucial for their pathogenicity. Despite this, the relative contribution of these virulence factors and other genetic changes to the pathogenicity of E. faecalis strains remain poorly understood. Here, we investigated whether specific genomic changes in the genome of E. faecalis isolates influence its pathogenicity-infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection and intestinal colonization. Our findings indicate that E. faecalis genetics partially influence the infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection, possibly due to gut-to-bloodstream translocation, highlighting the potential substantial role of host and environmental factors, including gut microbiota, on the opportunistic pathogenic lifestyle of this bacterium.

Unraveling Enterococcus susceptibility to quaternary ammonium compounds: genes, phenotypes, and the impact of environmental conditions.

Quaternary ammonium compounds (QACs) have been extensively used in the community, healthcare facilities, and food chain, in concentrations between 20 and 30,000 mg/L. Enterococcus faecalis and Enterococcus faecium are ubiquitous in these settings and are recognized as nosocomial pathogens worldwide, but QACs' activity against strains from diverse epidemiological and genomic backgrounds remained largely unexplored. We evaluated the role of Enterococcus isolates from different sources, years, and clonal lineages as hosts of QACs tolerance genes and their susceptibility to QACs in optimal, single-stress and cross-stress growth conditions. Only 1% of the Enterococcus isolates included in this study and 0.5% of publicly available Enterococcus genomes carried qacA/B, qacC, qacG, qacJ, qacZ, qrg, bcrABC or oqxAB genes, shared with >60 species of Bacillota, Pseudomonadota, Actinomycetota, or Spirochaetota. These genes were generally found within close proximity of antibiotics and/or metals resistance genes. The minimum inhibitory concentrations (MIC) and minimum bactericidal concentrations (MBC) of benzalkonium chloride (BC) and didecyldimethylammonium chloride ranged between 0.5 and 4 mg/L (microdilution: 37°C/20 h/pH = 7/aerobiosis) for 210 E. faecalis and E. faecium isolates (two isolates carrying qacZ). Modified growth conditions (e.g., 22°C/pH = 5) increased MICBC/MBCBC (maximum of eightfold and MBCBC = 16 mg/L) and changed bacterial growth kinetics under BC toward later stationary phases in both species, including in isolates without QACs tolerance genes. In conclusion, Enterococcus are susceptible to in-use QACs concentrations and rarely carry QACs tolerance genes. However, their potential gene exchange with different microbiota, the decreased susceptibility to QACs under specific environmental conditions, and the presence of subinhibitory QACs concentrations in various settings may contribute to the selection of particular strains and, thus, require a One Health strategy to maintain QACs effectiveness. IMPORTANCE Despite the increasing use of quaternary ammonium compounds (QACs), the susceptibility of pathogens to these antimicrobials remains largely unknown. Enterococcus faecium and Enterococcus faecalis are susceptible to in-use QACs concentrations and are not main hosts of QACs tolerance genes but participate in gene transfer pathways with diverse bacterial taxa exposed to these biocides. Moreover, QACs tolerance genes often share the same genetic contexts with antibiotics and/or metals resistance genes, raising concerns about potential co-selection events. E. faecium and E. faecalis showed increased tolerance to benzalkonium chloride under specific environmental conditions (22°C, pH = 5), suggesting that strains might be selected in settings where they occur along with subinhibitory QACs concentrations. Transcriptomic studies investigating the cellular mechanisms of Enterococcus adaptation to QACs tolerance, along with longitudinal metadata analysis of tolerant populations dynamics under the influence of diverse environmental factors, are essential and should be prioritized within a One Health strategy.

Gut Colonization by ESBL-Producing Escherichia coli in Dogs Is Associated with a Distinct Microbiome and Resistome Composition.

The gut microbiome of humans and animals acts as a reservoir of extended-spectrum beta-lactamase-producing Escherichia coli (ESBL-EC). Dogs are known for having a high prevalence of ESBL-EC in their gut microbiota, although their ESBL-EC carrier status often shifts over time. We hypothesized that the gut microbiome composition of dogs is implicated in ESBL-EC colonization status. Therefore, we assessed whether ESBL-EC carriage in dogs is associated with changes in the gut microbiome and resistome. Fecal samples were collected longitudinally from 57 companion dogs in the Netherlands every 2 weeks for a total of 6 weeks (n = 4 samples/dog). Carriage of ESBL-EC was determined through selective culturing and PCR and in line with previous studies, we observed a high prevalence of ESBL-EC carriage in dogs. Using 16s rRNA gene profiling we found significant associations between detected ESBL-EC carriage and an increased abundance of Clostridium sensu stricto 1, Enterococcus, Lactococcus, and the shared genera of Escherichia-Shigella in the dog microbiome. A resistome capture sequencing approach (ResCap) furthermore, revealed associations between detected ESBL-EC carriage and the increased abundance of the antimicrobial resistance genes: cmlA, dfrA, dhfR, floR, and sul3. In summary, our study showed that ESBL-EC carriage is associated with a distinct microbiome and resistome composition. IMPORTANCE The gut microbiome of humans and animals is an important source of multidrug resistant pathogens, including beta-lactamase-producing Escherichia coli (ESBL-EC). In this study, we assessed if the carriage of ESBL-EC in dogs was associated with changes in gut composition of bacteria and antimicrobial resistant genes (ARGs). Therefore, stool samples from 57 dogs were collected every 2 weeks for a total of 6 weeks. Sixty eight percent of the dogs carried ESBL-EC during at least one of the time points analyzed. By investigating the gut microbiome and resistome composition, we observed specific changes at time points when dogs were colonized with ESBL-EC compared to time points whenESBL-EC were not detected. In conclusion, our study highlights the importance to study the microbial diversity in companion animals, as gut colonization of particular antimicrobial resistant bacteria might be an indication of a changed microbial composition that is associated with the selection of particular ARGs.

Characterization of blaKPC-2 and blaNDM-1 Plasmids of a K. pneumoniae ST11 Outbreak Clone.

The most common resistance mechanism to carbapenems is the production of carbapenemases. In 2021, the Pan American Health Organization warned of the emergence and increase in new carbapenemase combinations in Enterobacterales in Latin America. In this study, we characterized four Klebsiella pneumoniae isolates harboring blaKPC and blaNDM from an outbreak during the COVID-19 pandemic in a Brazilian hospital. We assessed their plasmids' transference ability, fitness effects, and relative copy number in different hosts. The K. pneumoniae BHKPC93 and BHKPC104 strains were selected for whole genome sequencing (WGS) based on their pulsed-field gel electrophoresis profile. The WGS revealed that both isolates belong to ST11, and 20 resistance genes were identified in each isolate, including blaKPC-2 and blaNDM-1. The blaKPC gene was present on a ~56 Kbp IncN plasmid and the blaNDM-1 gene on a ~102 Kbp IncC plasmid, along with five other resistance genes. Although the blaNDM plasmid contained genes for conjugational transfer, only the blaKPC plasmid conjugated to E. coli J53, without apparent fitness effects. The minimum inhibitory concentrations (MICs) of meropenem/imipenem against BHKPC93 and BHKPC104 were 128/64 and 256/128 mg/L, respectively. Although the meropenem and imipenem MICs against E. coli J53 transconjugants carrying the blaKPC gene were 2 mg/L, this was a substantial increment in the MIC relative to the original J53 strain. The blaKPC plasmid copy number was higher in K. pneumoniae BHKPC93 and BHKPC104 than in E. coli and higher than that of the blaNDM plasmids. In conclusion, two ST11 K. pneumoniae isolates that were part of a hospital outbreak co-harbored blaKPC-2 and blaNDM-1. The blaKPC-harboring IncN plasmid has been circulating in this hospital since at least 2015, and its high copy number might have contributed to the conjugative transfer of this particular plasmid to an E. coli host. The observation that the blaKPC-containing plasmid had a lower copy number in this E. coli strain may explain why this plasmid did not confer phenotypic resistance against meropenem and imipenem.

Enterococcus spp. from chicken meat collected 20 years apart overcome multiple stresses occurring in the poultry production chain: Antibiotics, copper and acids.

Poultry meat has been a vehicle of antibiotic resistant bacteria and genes. Yet, the diversity of selective pressures associated with their maintenance in the poultry-production chain remains poorly explored. We evaluated the susceptibility of Enterococcus spp. from chicken meat collected 20 years apart to antibiotics, metals, acidic pH and peracetic acid-PAA. Contemporary chicken-meat samples (n = 53 batches, each including a pool of neck skin from 10 single carcasses) were collected in a slaughterhouse facility using PAA as disinfectant (March-August 2018, North of Portugal). Broilers were raised in intensive farms (n = 29) using CuSO4 and organic acids as feed additives. Data were compared with that of 67 samples recovered in the same region during 1999-2001. All 2018 samples had multidrug resistant-MDR isolates, with >45 % carrying Enterococcus faecalis, Enterococcus faecium or Enterococcus gallinarum resistant to tetracycline, erythromycin, ampicillin, quinupristin-dalfopristin, ciprofloxacin, chloramphenicol or aminoglycosides. Resistance rates were similar (P > 0.05) to those of 1999-2001 samples for all but five antibiotics. The decrease of samples carrying vancomycin-resistant isolates from 46 % to 0 % between 1999-2001 and 2018 was the most striking difference. Isolates from both periods were similarly susceptible to acid pH [minimum-growth pH (4.5-5.0), minimum-survival pH (3.0-4.0)] and to PAA (MIC90 = 100-120 mg/L/MBC90 = 140-160 mg/L; below concentrations used in slaughterhouse). Copper tolerance genes (tcrB and/or cueO) were respectively detected in 21 % and 4 % of 2018 and 1999-2001 samples. The tcrB gene was only detected in E. faecalis (MICCuSO4 > 12 mM), and their genomes were compared with other international ones of chicken origin (PATRIC database), revealing a polyclonal population and a plasmid or chromosomal location for tcrB. The tcrB plasmids shared diverse genetic modules, including multiple antimicrobial resistance genes (e.g. to tetracyclines, chloramphenicol, macrolide-lincosamide-streptogramin B-MLSB, aminoglycosides, bacitracin, coccidiostats). When in chromosome, the tcrB gene was co-located closely to merA (mercury) genes. Chicken meat remains an important vehicle of MDR Enterococcus spp. able to survive under diverse stresses (e.g. copper, acid) potentially contributing to these bacteria maintenance and flux among animal-environment-humans.

Species-Specific Patterns of Gut Metabolic Modules in Dutch Individuals with Different Dietary Habits.

Diet is an important determinant of the human gut microbiome. Here, we analyzed fecal metagenomes of Dutch adults following omnivorous, pescatarian, vegan, and vegetarian diets. We compared the taxonomic composition of individuals from our study with publicly available gut metagenomes from westernized and non-westernized societies. We observed that, despite long-term transition to diets rich in plant fibers (vegan or vegetarian), the microbiomes of these were typical of westernized populations, and similar in composition to omnivores. Although there were no major differences in metabolic modules, we identified differences in the species that contributed to particular functions, such as carbohydrate degradation and short-chain fatty acid metabolism. Overall, this study shows functional redundancy of the microbiomes among westernized populations, which is independent of long-term individual dietary habits. IMPORTANCE Diet is an important modulator of the human gut microbiome, which is susceptible to increased consumption of plant fibers in vegan or vegetarian lifestyles. To investigate this, we compared the gut microbiome of Dutch adults following omnivorous, pescatarian, vegan and vegetarian diets. We did not observe major differences in the gut microbiome composition and function between individuals with different dietary habits. However, we observed differences in the species that contribute to the core functions of the gut microbiome. Our study thus emphasizes the need to better understand the species-specific functional changes associated with dietary habits in the human gut microbiome.

Evolution of Chlorhexidine Susceptibility and of the EfrEF Operon among Enterococcus faecalis from Diverse Environments, Clones, and Time Spans.

Chlorhexidine (CHX) is widely used to control the spread of pathogens (e.g., human/animal clinical settings, ambulatory care, food industry). Enterococcus faecalis, a major nosocomial pathogen, is broadly distributed in diverse hosts and environments facilitating its exposure to CHX over the years. Nevertheless, CHX activity against E. faecalis is understudied. Our goal was to assess CHX activity and the variability of ChlR-EfrEF proteins (associated with CHX tolerance) among 673 field isolates and 1,784 E. faecalis genomes from the PATRIC database from different sources, time spans, clonal lineages, and antibiotic-resistance profiles. The CHX MIC (MICCHX) and minimum bactericidal concentration (MBCCHX) against E. faecalis presented normal distributions (0.5 to 64 mg/L). However, more CHX-tolerant isolates were detected in the food chain and recent human infections, suggesting an adaptability of E. faecalis populations in settings where CHX is heavily used. Heterogeneity in ChlR-EfrEF sequences was identified, with isolates harboring incomplete ChlR-EfrEF proteins, particularly the EfrE identified in the ST40 clonal lineage, showing low MICCHX (≤1mg/L). Distinct ST40-E. faecalis subpopulations carrying truncated and nontruncated EfrE were detected, with the former being predominant in human isolates. This study provides a new insight about CHX susceptibility and ChlR-EfrEF variability within diverse E. faecalis populations. The MICCHX/MBCCHX of more tolerant E. faecalis (MICCHX = 8 mg/L; MBCCHX = 64 mg/L) remain lower than in-use concentrations of CHX (≥500 mg/L). However, increased CHX use, combined with concentration gradients occurring in diverse environments, potentially selecting multidrug-resistant strains with different CHX susceptibilities, signals the importance of monitoring the trends of E. faecalis CHX tolerance within a One Health approach. IMPORTANCE Chlorhexidine (CHX) is a disinfectant and antiseptic used since the 1950s and included in the World Health Organization's list of essential medicines. It has been widely applied in hospitals, the community, the food industry, animal husbandry and pets. CHX tolerance in Enterococcus faecalis, a ubiquitous bacterium and one of the leading causes of human hospital-acquired infections, remains underexplored. Our study provides novel and comprehensive insights about CHX susceptibility within the E. faecalis population structure context, revealing more CHX-tolerant subpopulations from the food chain and recent human infections. We further show a detailed analysis of the genetic diversity of the efrEF operon (previously associated with E. faecalis CHX tolerance) and its correlation with CHX phenotypes. The recent strains with a higher tolerance to CHX and the multiple sources where bacteria are exposed to this biocide alert us to the need for the continuous monitoring of E. faecalis adaptation toward CHX tolerance within a One Health approach.

External validation of WGS-based antimicrobial susceptibility prediction tools, KOVER-AMR and ResFinder 4.1, for Escherichia coli clinical isolates.

OBJECTIVE: To externally validate whole genome sequence-antimicrobial susceptibility testing phenotype prediction tools KOVER-AMR and ResFinder 4.1 for Escherichia coli clinical isolates from Dutch routine care. METHODS: A random sample of 234 E. coli and 283 third generation cephalosporin-resistant E. coli isolates from urine and blood were collected (2014-2017). Culture-antimicrobial susceptibility testing was performed using VITEK 2 and BD Phoenix. Sequences were used as input for KOVER-AMR-SCM, KOVER-AMR-CART, and ResFinder 4.1. The concordance, major error rate (MER), and very major error rate (VMER) were calculated, with subsequent comparison to U.S. Food and Drug Administration (FDA) criteria (MER ≤3% and VMER with a 95% confidence interval ≤1.5-≤7.5%). RESULTS: ResFinder 4.1 performed better than KOVER-AMR-models; however, neither tool achieved overall (V)MERs below FDA criteria. KOVER-AMR-SCM, KOVER-AMR-CART, and ResFinder 4.1, MER (cumulative all antimicrobials) were: 5.1% (4.4-5.9), 4.3% (3.6-5.0), and 5.1% (4.5-5.8), respectively. MERs ≤3% were achieved for 6 (SCM) and 5 (CART) of the 11 tested antimicrobials for KOVER-AMR-models and for 9/13 antimicrobials tested with ResFinder 4.1. KOVER-AMR-SCM, KOVER-AMR-CART, and ResFinder 4.1 cumulative VMERs were: 26% (24-28), 29% (27-31), and 11% (9.2-12). VMERs with a 95% CI ≤ 1.5-≤7.5 were only achieved for 4/13 tested antimicrobials with ResFinder 4.1. DISCUSSION: In this study, whole genome sequence-antimicrobial susceptibility testing phenotype prediction tools KOVER-AMR and ResFinder 4.1 did not meet the FDA criteria needed for clinical diagnostic use in 517 E. coli clinical isolates from Dutch routine care. The tested tools should be further improved before they can be used for clinical decision making.

The interplay between community and hospital Enterococcus faecium clones within health-care settings: a genomic analysis.

BACKGROUND: The genomic relationships among Enterococcus faecium isolates are the subject of ongoing research that seeks to clarify the origins of observed lineages and the extent of horizontal gene transfer between them, and to robustly identify links between genotypes and phenotypes. E faecium is considered to form distinct groups-A and B-corresponding to isolates derived from patients who were hospitalised (A) and isolates from humans in the community (B). The additional separation of A into the so-called clades A1 and A2 remains an area of uncertainty. We aimed to investigate the relationships between A1 and non-A1 groups and explore the potential role of non-A1 isolates in shaping the population structure of hospital E faecium. METHODS: We collected short-read sequence data from invited groups that had previously published E faecium genome data. This hospital-based isolate collection could be separated into three groups (or clades, A1, A2, and B) by augmenting the study genomes with published sequences derived from human samples representing the previously defined genomic clusters. We performed phylogenetic analyses, by constructing maximum-likelihood phylogenetic trees, and identified historical recombination events. We assessed the pan-genome, did resistome analysis, and examined the genomic data to identify mobile genetic elements. Each genome underwent chromosome painting by use of ChromoPainter within FineSTRUCTURE software to assess ancestry and identify hybrid groups. We further assessed highly admixed regions to infer recombination directionality. FINDINGS: We assembled a collection of 1095 hospital E faecium sequences from 34 countries, further augmented by 33 published sequences. 997 (88%) of 1128 genomes clustered as A1, 92 (8%) as A2, and 39 (4%) as B. We showed that A1 probably emerged as a clone from within A2 and that, because of ongoing gene flow, hospital isolates currently identified as A2 represent a genetic continuum between A1 and community E faecium. This interchange of genetic material between isolates from different groups results in the emergence of hybrid genomes between clusters. Of the 1128 genomes, 49 (4%) hybrid genomes were identified: 33 previously labelled as A2 and 16 previously labelled as A1. These interactions were fuelled by a directional pattern of recombination mediated by mobile genetic elements. By contrast, the contribution of B group genetic material to A1 was limited to a few small regions of the genome and appeared to be driven by genomic sweep events. INTERPRETATION: A2 and B isolates coming into the hospital form an important reservoir for ongoing A1 adaptation, suggesting that effective long-term control of the effect of E faecium could benefit from strategies to reduce these genomic interactions, such as a focus on reducing the acquisition of hospital A1 strains by patients entering the hospital. FUNDING: Wellcome Trust.

Impact of long-term dietary habits on the human gut resistome in the Dutch population.

The human gut microbiome plays a central role in health and disease. Environmental factors, such as lifestyle and diet, are known to shape the gut microbiome as well as the reservoir of resistance genes that these microbes harbour; the resistome. In this study we assessed whether long-term dietary habits within a single geographical region (the Netherlands) impact the human gut resistome. Faecal samples from Dutch omnivores, pescatarians, vegetarians and vegans were analysed by metagenomic shotgun sequencing (MSS) (n = 149) and resistome capture sequencing approach (ResCap) (n = 64). Among all diet groups, 119 and 145 unique antibiotic resistance genes (ARGs) were detected by MSS or ResCap, respectively. Five or fifteen ARGs were shared between all diet groups, based on MSS and ResCap, respectively. The total number of detected ARGs by MSS or ResCap was not significantly different between the groups. MSS also revealed that vegans have a distinct microbiome composition, compared to other diet groups. Vegans had a lower abundance of Streptococcus thermophilus and Lactococcus lactis compared to pescatarians and a lower abundance of S. thermophilus when compared to omnivores. In summary, our study showed that long-term dietary habits are not associated with a specific resistome signature.

Effects of early-life antibiotics on the developing infant gut microbiome and resistome: a randomized trial.

Broad-spectrum antibiotics for suspected early-onset neonatal sepsis (sEONS) may have pronounced effects on gut microbiome development and selection of antimicrobial resistance when administered in the first week of life, during the assembly phase of the neonatal microbiome. Here, 147 infants born at ≥36 weeks of gestational age, requiring broad-spectrum antibiotics for treatment of sEONS in their first week of life were randomized 1:1:1 to receive three commonly prescribed intravenous antibiotic combinations, namely penicillin + gentamicin, co-amoxiclav + gentamicin or amoxicillin + cefotaxime (ZEBRA study, Trial Register NL4882). Average antibiotic treatment duration was 48 hours. A subset of 80 non-antibiotic treated infants from a healthy birth cohort served as controls (MUIS study, Trial Register NL3821). Rectal swabs and/or faeces were collected before and immediately after treatment, and at 1, 4 and 12 months of life. Microbiota were characterized by 16S rRNA-based sequencing and a panel of 31 antimicrobial resistance genes was tested using targeted qPCR. Confirmatory shotgun metagenomic sequencing was executed on a subset of samples. The overall gut microbial community composition and antimicrobial resistance gene profile majorly shift directly following treatment (R2 = 9.5%, adjusted p-value = 0.001 and R2 = 7.5%, adjusted p-value = 0.001, respectively) and normalize over 12 months (R2 = 1.1%, adjusted p-value = 0.03 and R2 = 0.6%, adjusted p-value = 0.23, respectively). We find a decreased abundance of Bifidobacterium spp. and increased abundance of Klebsiella and Enterococcus spp. in the antibiotic treated infants compared to controls. Amoxicillin + cefotaxime shows the largest effects on both microbial community composition and antimicrobial resistance gene profile, whereas penicillin + gentamicin exhibits the least effects. These data suggest that the choice of empirical antibiotics is relevant for adverse ecological side-effects.

Functional characterization of a gene cluster responsible for inositol catabolism associated with hospital-adapted isolates of Enterococcus faecium.

Enterococcus faecium is a nosocomial, multidrug-resistant pathogen. Whole genome sequence studies revealed that hospital-associated E. faecium isolates are clustered in a separate clade A1. Here, we investigated the distribution, integration site and function of a putative iol gene cluster that encodes for myo-inositol (MI) catabolism. This iol gene cluster was found as part of an ~20 kbp genetic element (iol element), integrated in ICEEfm1 close to its integrase gene in E. faecium isolate E1679. Among 1644 E. faecium isolates, ICEEfm1 was found in 789/1227 (64.3 %) clade A1 and 3/417 (0.7 %) non-clade A1 isolates. The iol element was present at a similar integration site in 180/792 (22.7 %) ICEEfm1-containing isolates. Examination of the phylogenetic tree revealed genetically closely related isolates that differed in presence/absence of ICEEfm1 and/or iol element, suggesting either independent acquisition or loss of both elements. E. faecium iol gene cluster containing isolates E1679 and E1504 were able to grow in minimal medium with only myo-inositol as carbon source, while the iolD-deficient mutant in E1504 (E1504∆iolD) lost this ability and an iol gene cluster negative recipient strain gained this ability after acquisition of ICEEfm1 by conjugation from donor strain E1679. Gene expression profiling revealed that the iol gene cluster is only expressed in the absence of other carbon sources. In an intestinal colonization mouse model the colonization ability of E1504∆iolD mutant was not affected relative to the wild-type E1504 strain. In conclusion, we describe and functionally characterise a gene cluster involved in MI catabolism that is associated with the ICEEfm1 island in hospital-associated E. faecium isolates. We were unable to show that this gene cluster provides a competitive advantage during gut colonisation in a mouse model. Therefore, to what extent this gene cluster contributes to the spread and ecological specialisation of ICEEfm1-carrying hospital-associated isolates remains to be investigated.

Recovering Escherichia coli Plasmids in the Absence of Long-Read Sequencing Data.

The incidence of infections caused by multidrug-resistant E. coli strains has risen in the past years. Antibiotic resistance in E. coli is often mediated by acquisition and maintenance of plasmids. The study of E. coli plasmid epidemiology and genomics often requires long-read sequencing information, but recently a number of tools that allow plasmid prediction from short-read data have been developed. Here, we reviewed 25 available plasmid prediction tools and categorized them into binary plasmid/chromosome classification tools and plasmid reconstruction tools. We benchmarked six tools (MOB-suite, plasmidSPAdes, gplas, FishingForPlasmids, HyAsP and SCAPP) that aim to reliably reconstruct distinct plasmids, with a special focus on plasmids carrying antibiotic resistance genes (ARGs) such as extended-spectrum beta-lactamase genes. We found that two thirds (n = 425, 66.3%) of all plasmids were correctly reconstructed by at least one of the six tools, with a range of 92 (14.58%) to 317 (50.23%) correctly predicted plasmids. However, the majority of plasmids that carried antibiotic resistance genes (n = 85, 57.8%) could not be completely recovered as distinct plasmids by any of the tools. MOB-suite was the only tool that was able to correctly reconstruct the majority of plasmids (n = 317, 50.23%), and performed best at reconstructing large plasmids (n = 166, 46.37%) and ARG-plasmids (n = 41, 27.9%), but predictions frequently contained chromosome contamination (40%). In contrast, plasmidSPAdes reconstructed the highest fraction of plasmids smaller than 18 kbp (n = 168, 61.54%). Large ARG-plasmids, however, were frequently merged with sequences derived from distinct replicons. Available bioinformatic tools can provide valuable insight into E. coli plasmids, but also have important limitations. This work will serve as a guideline for selecting the most appropriate plasmid reconstruction tool for studies focusing on E. coli plasmids in the absence of long-read sequencing data.

Apparent nosocomial adaptation of Enterococcus faecalis predates the modern hospital era.

Enterococcus faecalis is a commensal and nosocomial pathogen, which is also ubiquitous in animals and insects, representing a classical generalist microorganism. Here, we study E. faecalis isolates ranging from the pre-antibiotic era in 1936 up to 2018, covering a large set of host species including wild birds, mammals, healthy humans, and hospitalised patients. We sequence the bacterial genomes using short- and long-read techniques, and identify multiple extant hospital-associated lineages, with last common ancestors dating back as far as the 19th century. We find a population cohesively connected through homologous recombination, a metabolic flexibility despite a small genome size, and a stable large core genome. Our findings indicate that the apparent hospital adaptations found in hospital-associated E. faecalis lineages likely predate the "modern hospital" era, suggesting selection in another niche, and underlining the generalist nature of this nosocomial pathogen.