Mutational signature in colorectal cancer caused by genotoxic pks+ E. coli.

Various species of the intestinal microbiota have been associated with the development of colorectal cancer1,2, but it has not been demonstrated that bacteria have a direct role in the occurrence of oncogenic mutations. Escherichia coli can carry the pathogenicity island pks, which encodes a set of enzymes that synthesize colibactin3. This compound is believed to alkylate DNA on adenine residues4,5 and induces double-strand breaks in cultured cells3. Here we expose human intestinal organoids to genotoxic pks+ E. coli by repeated luminal injection over five months. Whole-genome sequencing of clonal organoids before and after this exposure revealed a distinct mutational signature that was absent from organoids injected with isogenic pks-mutant bacteria. The same mutational signature was detected in a subset of 5,876 human cancer genomes from two independent cohorts, predominantly in colorectal cancer. Our study describes a distinct mutational signature in colorectal cancer and implies that the underlying mutational process results directly from past exposure to bacteria carrying the colibactin-producing pks pathogenicity island.

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.

gplas: a comprehensive tool for plasmid analysis using short-read graphs.

SUMMARY: Plasmids can horizontally transmit genetic traits, enabling rapid bacterial adaptation to new environments and hosts. Short-read whole-genome sequencing data are often applied to large-scale bacterial comparative genomics projects but the reconstruction of plasmids from these data is facing severe limitations, such as the inability to distinguish plasmids from each other in a bacterial genome. We developed gplas, a new approach to reliably separate plasmid contigs into discrete components using sequence composition, coverage, assembly graph information and network partitioning based on a pruned network of plasmid unitigs. Gplas facilitates the analysis of large numbers of bacterial isolates and allows a detailed analysis of plasmid epidemiology based solely on short-read sequence data. AVAILABILITY AND IMPLEMENTATION: Gplas is written in R, Bash and uses a Snakemake pipeline as a workflow management system. Gplas is available under the GNU General Public License v3.0 at https://gitlab.com/sirarredondo/gplas.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Conditionally essential genes for survival during starvation in Enterococcus faecium E745.

BACKGROUND: The nosocomial pathogen Enterococcus faecium can survive for prolonged periods of time on surfaces in the absence of nutrients. This trait is thought to contribute to the ability of E. faecium to spread among patients in hospitals. There is currently a lack of data on the mechanisms that are responsible for the ability of E. faecium to survive in the absence of nutrients. RESULTS: We performed a high-throughput transposon mutant library screening (Tn-seq) to identify genes that have a role in long-term survival during incubation in phosphate-buffered saline (PBS) at 20 °C. A total of 24 genes were identified by Tn-seq to contribute to survival in PBS, with functions associated with the general stress response, DNA repair, metabolism, and membrane homeostasis. The gene which was quantitatively most important for survival in PBS was usp (locus tag: EfmE745_02439), which is predicted to encode a 17.4 kDa universal stress protein. After generating a targeted deletion mutant in usp, we were able to confirm that usp significantly contributes to survival in PBS and this defect was restored by in trans complementation. The usp gene is present in 99% of a set of 1644 E. faecium genomes that collectively span the diversity of the species. CONCLUSIONS: We postulate that usp is a key determinant for the remarkable environmental robustness of E. faecium. Further mechanistic studies into usp and other genes identified in this study may shed further light on the mechanisms by which E. faecium can survive in the absence of nutrients for prolonged periods of time.

Evolution of Colistin Resistance in the Klebsiella pneumoniae Complex Follows Multiple Evolutionary Trajectories with Variable Effects on Fitness and Virulence Characteristics.

The increasing prevalence of multidrug-resistant Klebsiella pneumoniae has led to a resurgence in the use of colistin as a last-resort drug. Colistin is a cationic antibiotic that selectively acts on Gram-negative bacteria through electrostatic interactions with anionic phosphate groups of the lipid A moiety of lipopolysaccharides (LPSs). Colistin resistance in K. pneumoniae is mediated through loss of these phosphate groups, their modification by cationic groups, and by the hydroxylation of acyl groups of lipid A. Here, we study the in vitro evolutionary trajectories toward colistin resistance in four clinical K. pneumoniae complex strains and their impact on fitness and virulence characteristics. Through population sequencing during in vitro evolution, we found that colistin resistance develops through a combination of single nucleotide polymorphisms, insertions and deletions, and the integration of insertion sequence elements, affecting genes associated with LPS biosynthesis and modification and capsule structures. Colistin resistance decreased the maximum growth rate of one K. pneumoniaesensu stricto strain, but not those of the other three K. pneumoniae complex strains. Colistin-resistant strains had lipid A modified through hydroxylation, palmitoylation, and l-Ara4N addition. K. pneumoniaesensu stricto strains exhibited cross-resistance to LL-37, in contrast to the Klebsiella variicola subsp. variicola strain. Virulence, as determined in a Caenorhabditis elegans survival assay, was increased in two colistin-resistant strains. Our study suggests that nosocomial K. pneumoniae complex strains can rapidly develop colistin resistance through diverse evolutionary trajectories upon exposure to colistin. This effectively shortens the life span of this last-resort antibiotic for the treatment of infections with multidrug-resistant Klebsiella.

Microevolution of acquired colistin resistance in Enterobacteriaceae from ICU patients receiving selective decontamination of the digestive tract.

BACKGROUND: Colistin is an antibiotic that targets the LPS molecules present in the membranes of Gram-negative bacteria. It is used as a last-resort drug to treat infections with MDR strains. Colistin is also used in selective decontamination of the digestive tract (SDD), a prophylactic therapy used in patients hospitalized in ICUs to selectively eradicate opportunistic pathogens in the oropharyngeal and gut microbiota. OBJECTIVES: To unravel the mechanisms of acquired colistin resistance in Gram-negative opportunistic pathogens obtained from SDD-treated patients. RESULTS: Routine surveillance of 428 SDD-treated patients resulted in 13 strains with acquired colistin resistance (Escherichia coli, n = 9; Klebsiella aerogenes, n = 3; Enterobacter asburiae, n = 1) from 5 patients. Genome sequence analysis showed that these isolates represented multiple distinct colistin-resistant clones but that colistin-resistant strains within the same patient were clonally related. We identified previously described mechanisms that lead to colistin resistance, i.e. a G53 substitution in the response regulator PmrA/BasR and the acquisition of the mobile colistin resistance gene mcr-1.1, but we also observed novel variants of basR with an 18 bp deletion and a G19E substitution in the sensor histidine kinase BasS. We experimentally confirmed that these variants contribute to reduced colistin susceptibility. In a single patient, we observed that colistin resistance in a single E. coli clone evolved through two unique variants in basRS. CONCLUSIONS: We show that prophylactic use of colistin during SDD can select for colistin resistance in species that are not intrinsically colistin resistant. This highlights the importance of continued surveillance for strains with acquired colistin resistance in patients treated with SDD.

Phylogenomics of Enterococcus faecalis from wild birds: new insights into host-associated differences in core and accessory genomes of the species.

Wild birds have been suggested to be reservoirs of antimicrobial resistant and/or pathogenic Enterococcus faecalis (Efs) strains, but the scarcity of studies and available sequences limit our understanding of the population structure of the species in these hosts. Here, we analysed the clonal and plasmid diversity of 97 Efs isolates from wild migratory birds. We found a high diversity, with most sequence types (STs) being firstly described here, while others were found in other hosts including some predominant in poultry. We found that pheromone-responsive plasmids predominate in wild bird Efs while 35% of the isolates entirely lack plasmids. Then, to better understand the ecology of the species, the whole genome of fivestrains with known STs (ST82, ST170, ST16 and ST55) were sequenced and compared with all the Efs genomes available in public databases. Using several methods to analyse core and accessory genomes (AccNET, PLACNET, hierBAPS and PANINI), we detected differences in the accessory genome of some lineages (e.g. ST82) demonstrating specific associations with birds. Conversely, the genomes of other Efs lineages exhibited divergence in core and accessory genomes, reflecting different adaptive trajectories in various hosts. This pangenome divergence, horizontal gene transfer events and occasional epidemic peaks could explain the population structure of the species.

Group IIA-Secreted Phospholipase A2 in Human Serum Kills Commensal but Not Clinical Enterococcus faecium Isolates.

Human innate immunity employs cellular and humoral mechanisms to facilitate rapid killing of invading bacteria. The direct killing of bacteria by human serum is attributed mainly to the activity of the complement system, which forms pores in Gram-negative bacteria. Although Gram-positive bacteria are considered resistant to killing by serum, we uncover here that normal human serum effectively kills Enterococcus faecium Comparison of a well-characterized collection of commensal and clinical E. faecium isolates revealed that human serum specifically kills commensal E. faecium strains isolated from normal gut microbiota but not clinical isolates. Inhibitor studies show that the human group IIA secreted phospholipase A2 (hGIIA), but not complement, is responsible for killing of commensal E. faecium strains in human normal serum. This is remarkable since the hGIIA concentration in "noninflamed" serum was considered too low to be bactericidal against Gram-positive bacteria. Mechanistic studies showed that serum hGIIA specifically causes permeabilization of commensal E. faecium membranes. Altogether, we find that a normal concentration of hGIIA in serum effectively kills commensal E. faecium and that resistance of clinical E. faecium to hGIIA could have contributed to the ability of these strains to become opportunistic pathogens in hospitalized patients.

Identification of a Novel Genomic Island Associated with vanD-Type Vancomycin Resistance in Six Dutch Vancomycin-Resistant Enterococcus faecium Isolates.

Genomic comparison of the first six Dutch vanD-type vancomycin-resistant Enterococcus faecium (VRE) isolates with four vanD gene clusters from other enterococcal species and anaerobic gut commensals revealed that the vanD gene cluster was located on a genomic island of variable size. Phylogenetic inferences revealed that the Dutch VRE isolates were genetically not closely related and that genetic variation of the vanD-containing genomic island was not species specific, suggesting that this island is transferred horizontally between enterococci and anaerobic gut commensals.

RNA-seq and Tn-seq reveal fitness determinants of vancomycin-resistant Enterococcus faecium during growth in human serum.

BACKGROUND: The Gram-positive bacterium Enterococcus faecium is a commensal of the human gastrointestinal tract and a frequent cause of bloodstream infections in hospitalized patients. The mechanisms by which E. faecium can survive and grow in blood during an infection have not yet been characterized. Here, we identify genes that contribute to growth of E. faecium in human serum through transcriptome profiling (RNA-seq) and a high-throughput transposon mutant library sequencing approach (Tn-seq). RESULTS: We first sequenced the genome of E. faecium E745, a vancomycin-resistant clinical isolate, using a combination of short- and long read sequencing, revealing a 2,765,010 nt chromosome and 6 plasmids, with sizes ranging between 9.3 kbp and 223.7 kbp. We then compared the transcriptome of E. faecium E745 during exponential growth in rich medium and in human serum by RNA-seq. This analysis revealed that 27.8% of genes on the E. faecium E745 genome were differentially expressed in these two conditions. A gene cluster with a role in purine biosynthesis was among the most upregulated genes in E. faecium E745 upon growth in serum. The E. faecium E745 transposon mutant library was then used to identify genes that were specifically required for growth of E. faecium in serum. Genes involved in de novo nucleotide biosynthesis (including pyrK_2, pyrF, purD, purH) and a gene encoding a phosphotransferase system subunit (manY_2) were thus identified to be contributing to E. faecium growth in human serum. Transposon mutants in pyrK_2, pyrF, purD, purH and manY_2 were isolated from the library and their impaired growth in human serum was confirmed. In addition, the pyrK_2 and manY_2 mutants were tested for their virulence in an intravenous zebrafish infection model and exhibited significantly attenuated virulence compared to E. faecium E745. CONCLUSIONS: Genes involved in carbohydrate metabolism and nucleotide biosynthesis of E. faecium are essential for growth in human serum and contribute to the pathogenesis of this organism. These genes may serve as targets for the development of novel anti-infectives for the treatment of E. faecium bloodstream infections.

The Two-Component System ChtRS Contributes to Chlorhexidine Tolerance in Enterococcus faecium.

Enterococcus faecium is one of the primary causes of nosocomial infections. Disinfectants are commonly used to prevent infections with multidrug-resistant E. faecium in hospitals. Worryingly, E. faecium strains that exhibit tolerance to disinfectants have already been described. We aimed to identify and characterize E. faecium genes that contribute to tolerance to the disinfectant chlorhexidine (CHX). We used a transposon mutant library, constructed in a multidrug-resistant E. faecium bloodstream isolate, to perform a genome-wide screen to identify genetic determinants involved in tolerance to CHX. We identified a putative two-component system (2CS), composed of a putative sensor histidine kinase (ChtS) and a cognate DNA-binding response regulator (ChtR), which contributed to CHX tolerance in E. faecium Targeted chtR and chtS deletion mutants exhibited compromised growth in the presence of CHX. Growth of the chtR and chtS mutants was also affected in the presence of the antibiotic bacitracin. The CHX- and bacitracin-tolerant phenotype of E. faecium E1162 was linked to a unique, nonsynonymous single nucleotide polymorphism in chtR Transmission electron microscopy showed that upon challenge with CHX, the ΔchtR and ΔchtS mutants failed to divide properly and formed long chains. Normal growth and cell morphology were restored when the mutations were complemented in trans Morphological abnormalities were also observed upon exposure of the ΔchtR and ΔchtS mutants to bacitracin. The tolerance to both chlorhexidine and bacitracin provided by ChtRS in E. faecium highlights the overlap between responses to disinfectants and antibiotics and the potential for the development of cross-tolerance for these classes of antimicrobials.

Dissemination of cephalosporin resistance genes between Escherichia coli strains from farm animals and humans by specific plasmid lineages.

Third-generation cephalosporins are a class of ß-lactam antibiotics that are often used for the treatment of human infections caused by Gram-negative bacteria, especially Escherichia coli. Worryingly, the incidence of human infections caused by third-generation cephalosporin-resistant E. coli is increasing worldwide. Recent studies have suggested that these E. coli strains, and their antibiotic resistance genes, can spread from food-producing animals, via the food-chain, to humans. However, these studies used traditional typing methods, which may not have provided sufficient resolution to reliably assess the relatedness of these strains. We therefore used whole-genome sequencing (WGS) to study the relatedness of cephalosporin-resistant E. coli from humans, chicken meat, poultry and pigs. One strain collection included pairs of human and poultry-associated strains that had previously been considered to be identical based on Multi-Locus Sequence Typing, plasmid typing and antibiotic resistance gene sequencing. The second collection included isolates from farmers and their pigs. WGS analysis revealed considerable heterogeneity between human and poultry-associated isolates. The most closely related pairs of strains from both sources carried 1263 Single-Nucleotide Polymorphisms (SNPs) per Mbp core genome. In contrast, epidemiologically linked strains from humans and pigs differed by only 1.8 SNPs per Mbp core genome. WGS-based plasmid reconstructions revealed three distinct plasmid lineages (IncI1- and IncK-type) that carried cephalosporin resistance genes of the Extended-Spectrum Beta-Lactamase (ESBL)- and AmpC-types. The plasmid backbones within each lineage were virtually identical and were shared by genetically unrelated human and animal isolates. Plasmid reconstructions from short-read sequencing data were validated by long-read DNA sequencing for two strains. Our findings failed to demonstrate evidence for recent clonal transmission of cephalosporin-resistant E. coli strains from poultry to humans, as has been suggested based on traditional, low-resolution typing methods. Instead, our data suggest that cephalosporin resistance genes are mainly disseminated in animals and humans via distinct plasmids.

Genome-wide Screening Identifies Phosphotransferase System Permease BepA to Be Involved in Enterococcus faecium Endocarditis and Biofilm Formation.

Enterococcus faecium is a common cause of nosocomial infections, of which infective endocarditis is associated with substantial mortality. In this study, we used a microarray-based transposon mapping (M-TraM) approach to evaluate a rat endocarditis model and identified a gene, originally annotated as "fruA" and renamed "bepA," putatively encoding a carbohydrate phosphotransferase system (PTS) permease (biofilm and endocarditis-associated permease A [BepA]), as important in infective endocarditis. This gene is highly enriched in E. faecium clinical isolates and absent in commensal isolates that are not associated with infection. Confirmation of the phenotype was established in a competition experiment of wild-type and a markerless bepA mutant in a rat endocarditis model. In addition, deletion of bepA impaired biofilm formation in vitro in the presence of 100% human serum and metabolism of ß-methyl-D-glucoside. ß-glucoside metabolism has been linked to the metabolism of glycosaminoglycans that are exposed on injured heart valves, where bacteria attach and form vegetations. Therefore, we propose that the PTS permease BepA is directly implicated in E. faecium pathogenesis.

Whole-Genome Multilocus Sequence Typing of Extended-Spectrum-Beta-Lactamase-Producing Enterobacteriaceae.

Molecular typing has become indispensable in the detection of nosocomial transmission of bacterial pathogens and the identification of sources and routes of transmission in outbreak settings, but current methods are labor-intensive, are difficult to standardize, or have limited resolution. Whole-genome multilocus sequence typing (wgMLST) has emerged as a whole-genome sequencing (WGS)-based gene-by-gene typing method that may overcome these limitations and has been applied successfully for several species in outbreak settings. In this study, genus-, genetic-complex-, and species-specific wgMLST schemes were developed for Citrobacter spp., the Enterobacter cloacae complex, Escherichia coli, Klebsiella oxytoca, and Klebsiella pneumoniae and used to type a national collection of 1,798 extended-spectrum-beta-lactamase-producing Enterobacteriaceae (ESBL-E) isolates obtained from patients in Dutch hospitals. Genus-, genetic-complex-, and species-specific thresholds for genetic distance that accurately distinguish between epidemiologically related and unrelated isolates were defined for Citrobacter spp., the E. cloacae complex, E. coli, and K. pneumoniae wgMLST was shown to have higher discriminatory power and typeability than in silico MLST. In conclusion, the wgMLST schemes developed in this study facilitate high-resolution WGS-based typing of the most prevalent ESBL-producing species in clinical practice and may contribute to further elucidation of the complex epidemiology of antimicrobial-resistant Enterobacteriaceae wgMLST opens up possibilities for the creation of a Web-accessible database for the global surveillance of ESBL-producing bacterial clones.

Multilevel population genetic analysis of vanA and vanB Enterococcus faecium causing nosocomial outbreaks in 27 countries (1986-2012).

OBJECTIVES: Vancomycin-resistant Enterococcus faecium (VREfm) have been increasingly reported since the 1980s. Despite the high number of published studies about VRE epidemiology, the dynamics and evolvability of these microorganisms are still not fully understood. A multilevel population genetic analysis of VREfm outbreak strains since 1986, representing the first comprehensive characterization of plasmid content in E. faecium, was performed to provide a detailed view of potential transmissible units. METHODS: From a comprehensive MeSH search, we identified VREfm strains causing hospital outbreaks (1986-2012). In total, 53 VanA and 18 VanB isolates (27 countries, 5 continents) were analysed and 82 vancomycin-susceptible E. faecium (VSEfm) were included for comparison. Clonal relatedness was established by PFGE and MLST (goeBURST/Bayesian Analysis of Population Structure, BAPS). Characterization of van transposons (PCR mapping, RFLP, sequencing), plasmids (transfer, ClaI-RFLP, PCR typing of relaxases, replication-initiation proteins and toxin-antitoxin systems, hybridization, sequencing), bacteriocins and virulence determinants (PCR, hybridization, sequencing) was performed. RESULTS: VREfm were mainly associated with major human lineages ST17, ST18 and ST78. VREfm and VSEfm harboured plasmids of different families [RCR, small theta plasmids, RepA_N (pRUM/pLG1) and Inc18] able to yield mosaic elements. Tn1546-vanA was mainly located on pRUM/Axe-Txe (USA) and Inc18-pIP186 (Europe) plasmids. The VanB2 type (Tn5382/Tn1549) was predominant among VanB strains (chromosome and plasmids). CONCLUSIONS: Both strains and plasmids contributed to the spread and persistence of vancomycin resistance among E. faecium. Horizontal gene transfer events among genetic elements from different clonal lineages (same or different species) result in chimeras with different stability and host range, complicating the surveillance of epidemic plasmids.

The impact of host metapopulation structure on the population genetics of colonizing bacteria.

Many key bacterial pathogens are frequently carried asymptomatically, and the emergence and spread of these opportunistic pathogens can be driven, or mitigated, via demographic changes within the host population. These inter-host transmission dynamics combine with basic evolutionary parameters such as rates of mutation and recombination, population size and selection, to shape the genetic diversity within bacterial populations. Whilst many studies have focused on how molecular processes underpin bacterial population structure, the impact of host migration and the connectivity of the local populations has received far less attention. A stochastic neutral model incorporating heightened local transmission has been previously shown to fit closely with genetic data for several bacterial species. However, this model did not incorporate transmission limiting population stratification, nor the possibility of migration of strains between subpopulations, which we address here by presenting an extended model. We study the consequences of migration in terms of shared genetic variation and show by simulation that the previously used summary statistic, the allelic mismatch distribution, can be insensitive to even large changes in microepidemic and migration rates. Using likelihood-free inference with genotype network topological summaries we fit a simpler model to commensal and hospital samples from the common nosocomial pathogens Staphylococcus aureus, Staphylococcus epidermidis, Enterococcus faecalis and Enterococcus faecium. Only the hospital data for E. faecium display clearly marked deviations from the model predictions which may be attributable to its adaptation to the hospital environment.

Growth condition-dependent cell surface proteome analysis of Enterococcus faecium.

The last 30 years Enterococcus faecium has become an important nosocomial pathogen in hospitals worldwide. The aim of this study was to obtain insight in the cell surface proteome of E. faecium when grown in laboratory and clinically relevant conditions. Enterococcus faecium E1162, a clinical blood stream isolate, was grown until mid-log phase in brain heart infusion medium (BHI) with, or without 0.02% bile salts, Tryptic Soy Broth with 1% glucose (TSBg) and urine, and its cell surface was "shaved" using immobilized trypsin. Peptides were identified using MS/MS. Mapping against the translated E1162 whole genome sequence identified 67 proteins that were differentially detected in different conditions. In urine, 14 proteins were significantly more and nine proteins less abundant relative to the other conditions. Growth in BHI-bile and TSBg, revealed four and six proteins, respectively, which were uniquely present in these conditions while two proteins were uniquely present in both conditions. Thus, proteolytic shaving of E. faecium cells identified differentially surface exposed proteins in different growth conditions. These proteins are of special interest as they provide more insight in the adaptive mechanisms and may serve as targets for the development of novel therapeutics against this multi-resistant emerging pathogen. All MS data have been deposited in the ProteomeXchange with identifier PXD002497 (http://proteomecentral.proteomexchange.org/dataset/PXD002497).

A genomic virulence reference map of Enterococcus faecalis reveals an important contribution of phage03-like elements in nosocomial genetic lineages to pathogenicity in a Caenorhabditis elegans infection model.

In the present study, the commensal and pathogenic host-microbe interaction of Enterococcus faecalis was explored using a Caenorhabditis elegans model system. The virulence of 28 E. faecalis isolates representing 24 multilocus sequence types (MLSTs), including human commensal and clinical isolates as well as isolates from animals and of insect origin, was investigated using C. elegans strain glp-4 (bn2ts); sek-1 (km4). This revealed that 6 E. faecalis isolates behaved in a commensal manner with no nematocidal effect, while the remaining strains showed a time to 50% lethality ranging from 47 to 120 h. Principal component analysis showed that the difference in nematocidal activity explained 94% of the variance in the data. Assessment of known virulence traits revealed that gelatinase and cytolysin production accounted for 40.8% and 36.5% of the observed pathogenicity, respectively. However, coproduction of gelatinase and cytolysin did not increase virulence additively, accounting for 50.6% of the pathogenicity and therefore indicating a significant (26.7%) saturation effect. We employed a comparative genomic analysis approach using the 28 isolates comprising a collection of 82,356 annotated coding sequences (CDS) to identify 2,325 patterns of presence or absence among the investigated strains. Univariate statistical analysis of variance (ANOVA) established that individual patterns positively correlated (n = 61) with virulence. The patterns were investigated to identify potential new virulence traits, among which we found five patterns consisting of the phage03-like gene clusters. Strains harboring phage03 showed, on average, 17% higher killing of C. elegans (P = 4.4e(-6)). The phage03 gene cluster was also present in gelatinase-and-cytolysin-negative strain E. faecalis JH2-2. Deletion of this phage element from the JH2-2 clinical strain rendered the mutant apathogenic in C. elegans, and a similar mutant of the nosocomial V583 isolate showed significantly attenuated virulence. Bioinformatics investigation indicated that, unlike other E. faecalis virulence traits, phage03-like elements were found at a higher frequency among nosocomial isolates. In conclusion, our report provides a valuable virulence map that explains enhancement in E. faecalis virulence and contributes to a deeper comprehension of the genetic mechanism leading to the transition from commensalism to a pathogenic lifestyle.

Investigating the mobilome in clinically important lineages of Enterococcus faecium and Enterococcus faecalis.

BACKGROUND: The success of Enterococcus faecium and E. faecalis evolving as multi-resistant nosocomial pathogens is associated with their ability to acquire and share adaptive traits, including antimicrobial resistance genes encoded by mobile genetic elements (MGEs). Here, we investigate this mobilome in successful hospital associated genetic lineages, E. faecium sequence type (ST)17 (n=10) and ST78 (n=10), E. faecalis ST6 (n=10) and ST40 (n=10) by DNA microarray analyses. RESULTS: The hybridization patterns of 272 representative targets including plasmid backbones (n=85), transposable elements (n=85), resistance determinants (n=67), prophages (n=29) and clustered regularly interspaced short palindromic repeats (CRISPR)-cas sequences (n=6) separated the strains according to species, and for E. faecalis also according to STs. RCR-, Rep_3-, RepA_N- and Inc18-family plasmids were highly prevalent and with the exception of Rep_3, evenly distributed between the species. There was a considerable difference in the replicon profile, with rep 17/pRUM , rep 2/pRE25 , rep 14/EFNP1 and rep 20/pLG1 dominating in E. faecium and rep 9/pCF10 , rep 2/pRE25 and rep 7 in E. faecalis strains. We observed an overall high correlation between the presence and absence of genes coding for resistance towards antibiotics, metals, biocides and their corresponding MGEs as well as their phenotypic antimicrobial susceptibility pattern. Although most IS families were represented in both E. faecalis and E. faecium, specific IS elements within these families were distributed in only one species. The prevalence of IS256-, IS3-, ISL3-, IS200/IS605-, IS110-, IS982- and IS4-transposases was significantly higher in E. faecium than E. faecalis, and that of IS110-, IS982- and IS1182-transposases in E. faecalis ST6 compared to ST40. Notably, the transposases of IS981, ISEfm1 and IS1678 that have only been reported in few enterococcal isolates were well represented in the E. faecium strains. E. faecalis ST40 strains harboured possible functional CRISPR-Cas systems, and still resistance and prophage sequences were generally well represented. CONCLUSIONS: The targeted MGEs were highly prevalent among the selected STs, underlining their potential importance in the evolution of hospital-adapted lineages of enterococci. Although the propensity of inter-species horizontal gene transfer (HGT) must be emphasized, the considerable species-specificity of these MGEs indicates a separate vertical evolution of MGEs within each species, and for E. faecalis within each ST.