Comprehensive molecular, genomic and phenotypic analysis of a major clone of Enterococcus faecalis MLST ST40.

BACKGROUND: Enterococcus faecalis is a multifaceted microorganism known to act as a beneficial intestinal commensal bacterium. It is also a dreaded nosocomial pathogen causing life-threatening infections in hospitalised patients. Isolates of a distinct MLST type ST40 represent the most frequent strain type of this species, distributed worldwide and originating from various sources (animal, human, environmental) and different conditions (colonisation/infection). Since enterococci are known to be highly recombinogenic we determined to analyse the microevolution and niche adaptation of this highly distributed clonal type. RESULTS: We compared a set of 42 ST40 isolates by assessing key molecular determinants, performing whole genome sequencing (WGS) and a number of phenotypic assays including resistance profiling, formation of biofilm and utilisation of carbon sources. We generated the first circular closed reference genome of an E. faecalis isolate D32 of animal origin and compared it with the genomes of other reference strains. D32 was used as a template for detailed WGS comparisons of high-quality draft genomes of 14 ST40 isolates. Genomic and phylogenetic analyses suggest a high level of similarity regarding the core genome, also demonstrated by similar carbon utilisation patterns. Distribution of known and putative virulence-associated genes did not differentiate between ST40 strains from a commensal and clinical background or an animal or human source. Further analyses of mobile genetic elements (MGE) revealed genomic diversity owed to: (1) a modularly structured pathogenicity island; (2) a site-specifically integrated and previously unknown genomic island of 138 kb in two strains putatively involved in exopolysaccharide synthesis; and (3) isolate-specific plasmid and phage patterns. Moreover, we used different cell-biological and animal experiments to compare the isolate D32 with a closely related ST40 endocarditis isolate whose draft genome sequence was also generated. D32 generally showed a greater capacity of adherence to human cell lines and an increased pathogenic potential in various animal models in combination with an even faster growth in vivo (not in vitro). CONCLUSION: Molecular, genomic and phenotypic analysis of representative isolates of a major clone of E. faecalis MLST ST40 revealed new insights into the microbiology of a commensal bacterium which can turn into a conditional pathogen.

Using a genome-scale metabolic model of Enterococcus faecalis V583 to assess amino acid uptake and its impact on central metabolism.

Increasing antibiotic resistance in pathogenic bacteria necessitates the development of new medication strategies. Interfering with the metabolic network of the pathogen can provide novel drug targets but simultaneously requires a deeper and more detailed organism-specific understanding of the metabolism, which is often surprisingly sparse. In light of this, we reconstructed a genome-scale metabolic model of the pathogen Enterococcus faecalis V583. The manually curated metabolic network comprises 642 metabolites and 706 reactions. We experimentally determined metabolic profiles of E. faecalis grown in chemically defined medium in an anaerobic chemostat setup at different dilution rates and calculated the net uptake and product fluxes to constrain the model. We computed growth-associated energy and maintenance parameters and studied flux distributions through the metabolic network. Amino acid auxotrophies were identified experimentally for model validation and revealed seven essential amino acids. In addition, the important metabolic hub of glutamine/glutamate was altered by constructing a glutamine synthetase knockout mutant. The metabolic profile showed a slight shift in the fermentation pattern toward ethanol production and increased uptake rates of multiple amino acids, especially l-glutamine and l-glutamate. The model was used to understand the altered flux distributions in the mutant and provided an explanation for the experimentally observed redirection of the metabolic flux. We further highlighted the importance of gene-regulatory effects on the redirection of the metabolic fluxes upon perturbation. The genome-scale metabolic model presented here includes gene-protein-reaction associations, allowing a further use for biotechnological applications, for studying essential genes, proteins, or reactions, and the search for novel drug targets.

Rotation gene set testing for longitudinal expression data.

Gene set analysis methods are popular tools for identifying differentially expressed gene sets in microarray data. Most existing methods use a permutation test to assess significance for each gene set. The permutation test's assumption of exchangeable samples is often not satisfied for time-series data and complex experimental designs, and in addition it requires a certain number of samples to compute p-values accurately. The method presented here uses a rotation test rather than a permutation test to assess significance. The rotation test can compute accurate p-values also for very small sample sizes. The method can handle complex designs and is particularly suited for longitudinal microarray data where the samples may have complex correlation structures. Dependencies between genes, modeled with the use of gene networks, are incorporated in the estimation of correlations between samples. In addition, the method can test for both gene sets that are differentially expressed and gene sets that show strong time trends. We show on simulated longitudinal data that the ability to identify important gene sets may be improved by taking the correlation structure between samples into account. Applied to real data, the method identifies both gene sets with constant expression and gene sets with strong time trends.

Enterococcus faecalis grows on ascorbic acid.

We show that Enterococcus faecalis can utilize ascorbate for fermentative growth. In chemically defined media, growth yield was limited by the supply of amino acids, and the cells showed a much higher demand for amino acids than when they were grown on glucose.

Smoothing gene expression data with network information improves consistency of regulated genes.

Gene set analysis methods have become a widely used tool for including prior biological knowledge in the statistical analysis of gene expression data. Advantages of these methods include increased sensitivity, easier interpretation and more conformity in the results. However, gene set methods do not employ all the available information about gene relations. Genes are arranged in complex networks where the network distances contain detailed information about inter-gene dependencies. We propose a method that uses gene networks to smooth gene expression data with the aim of reducing the number of false positives and identify important subnetworks. Gene dependencies are extracted from the network topology and are used to smooth genewise test statistics. To find the optimal degree of smoothing, we propose using a criterion that considers the correlation between the network and the data. The network smoothing is shown to improve the ability to identify important genes in simulated data. Applied to a real data set, the smoothing accentuates parts of the network with a high density of differentially expressed genes.

Adaptation of Brucella melitensis Antimicrobial Susceptibility Testing to the ISO 20776 Standard and Validation of the Method.

Brucellosis, mainly caused by Brucella (B.) melitensis, is associated with a risk of chronification and relapses. Antimicrobial susceptibility testing (AST) standards for B. melitensis are not available, and the agent is not yet listed in the EUCAST breakpoint tables. CLSI recommendations for B. melitensis exist, but they do not fulfill the requirements of the ISO 20776 standard regarding the culture medium and the incubation conditions. Under the third EU Health Programme, laboratories specializing in the diagnostics of highly pathogenic bacteria in their respective countries formed a working group within a Joint Action aiming to develop a suitable method for the AST of B. melitensis. Under the supervision of EUCAST representatives, this working group adapted the CLSI M45 document to the ISO 20776 standard after testing and validation. These adaptations included the comparison of various culture media, culture conditions and AST methods. A Standard Operation Procedure was derived and an interlaboratory validation was performed in order to evaluate the method. The results showed pros and cons for both of the two methods but also indicate that it is not necessary to abandon Mueller-Hinton without additives for the AST of B. melitensis.

Comparative genomic analysis reveals significant enrichment of mobile genetic elements and genes encoding surface structure-proteins in hospital-associated clonal complex 2 Enterococcus faecalis.

BACKGROUND: Enterococci rank among the leading causes of nosocomial infections. The failure to identify pathogen-specific genes in Enterococcus faecalis has led to a hypothesis where the virulence of different strains may be linked to strain-specific genes, and where the combined endeavor of the different gene-sets result in the ability to cause infection. Population structure studies by multilocus sequence typing have defined distinct clonal complexes (CC) of E. faecalis enriched in hospitalized patients (CC2, CC9, CC28 and CC40). RESULTS: In the present study, we have used a comparative genomic approach to investigate gene content in 63 E. faecalis strains, with a special focus on CC2. Statistical analysis using Fisher's exact test revealed 252 significantly enriched genes among CC2-strains. The majority of these genes were located within the previously defined mobile elements phage03 (n = 51), efaB5 (n = 34) and a vanB associated genomic island (n = 55). Moreover, a CC2-enriched genomic islet (EF3217 to -27), encoding a putative phage related element within the V583 genome, was identified. From the draft genomes of CC2-strains HH22 and TX0104, we also identified a CC2-enriched non-V583 locus associated with the E. faecalis pathogenicity island (PAI). Interestingly, surface related structures (including MSCRAMMs, internalin-like and WxL protein-coding genes) implicated in virulence were significantly overrepresented (9.1%; p = 0.036, Fisher's exact test) among the CC2-enriched genes. CONCLUSION: In conclusion, we have identified a set of genes with potential roles in adaptation or persistence in the hospital environment, and that might contribute to the ability of CC2 E. faecalis isolates to cause disease.

Rotation testing in gene set enrichment analysis for small direct comparison experiments.

Gene Set Enrichment Analysis (GSEA) is a method for analysing gene expression data with a focus on a priori defined gene sets. The permutation test generally used in GSEA for testing the significance of gene set enrichment involves permutation of a phenotype vector and is developed for data from an indirect comparison design, i.e. unpaired data. In some studies the samples representing two phenotypes are paired, e.g. samples taken from a patient before and after treatment, or if samples representing two phenotypes are hybridised to the same two-channel array (direct comparison design). In this paper we will focus on data from direct comparison experiments, but the methods can be applied to paired data in general. For these types of data, a standard permutation test for paired data that randomly re-signs samples can be used. However, if the sample size is very small, which is often the case for a direct comparison design, a permutation test will give very imprecise estimates of the p-values. Here we propose using a rotation test rather than a permutation test for estimation of significance in GSEA of direct comparison data with a limited number of samples. Our proposed rotation test makes GSEA applicable to direct comparison data with few samples, by depending on rotations of the data instead of permutations. The rotation test is a generalisation of the permutation test, and can in addition be used on indirect comparison data and for testing significance of other types of test statistics outside the GSEA framework.

Improved analysis of bacterial CGH data beyond the log-ratio paradigm.

BACKGROUND: Existing methods for analyzing bacterial CGH data from two-color arrays are based on log-ratios only, a paradigm inherited from expression studies. We propose an alternative approach, where microarray signals are used in a different way and sequence identity is predicted using a supervised learning approach. RESULTS: A data set containing 32 hybridizations of sequenced versus sequenced genomes have been used to test and compare methods. A ROC-analysis has been performed to illustrate the ability to rank probes with respect to Present/Absent calls. Classification into Present and Absent is compared with that of a gaussian mixture model. CONCLUSION: The results indicate our proposed method is an improvement of existing methods with respect to ranking and classification of probes, especially for multi-genome arrays.

Comparative genomics of Enterococcus faecalis from healthy Norwegian infants.

BACKGROUND: Enterococcus faecalis, traditionally considered a harmless commensal of the intestinal tract, is now ranked among the leading causes of nosocomial infections. In an attempt to gain insight into the genetic make-up of commensal E. faecalis, we have studied genomic variation in a collection of community-derived E. faecalis isolated from the feces of Norwegian infants. RESULTS: The E. faecalis isolates were first sequence typed by multilocus sequence typing (MLST) and characterized with respect to antibiotic resistance and properties associated with virulence. A subset of the isolates was compared to the vancomycin resistant strain E. faecalis V583 (V583) by whole genome microarray comparison (comparative genomic hybridization (CGH)). Several of the putative enterococcal virulence factors were found to be highly prevalent among the commensal baby isolates. The genomic variation as observed by CGH was less between isolates displaying the same MLST sequence type than between isolates belonging to different evolutionary lineages. CONCLUSION: The variations in gene content observed among the investigated commensal E. faecalis is comparable to the genetic variation previously reported among strains of various origins thought to be representative of the major E. faecalis lineages. Previous MLST analysis of E. faecalis have identified so-called high-risk enterococcal clonal complexes (HiRECC), defined as genetically distinct subpopulations, epidemiologically associated with enterococcal infections. The observed correlation between CGH and MLST presented here, may offer a method for the identification of lineage-specific genes, and may therefore add clues on how to distinguish pathogenic from commensal E. faecalis. In this work, information on the core genome of E. faecalis is also substantially extended.

Carriage of ESBL/AmpC-producing or ciprofloxacin non-susceptible Escherichia coli and Klebsiella spp. in healthy people in Norway.

BACKGROUND: Asymptomatic carriage has been recognised as an important risk factor for infection caused by antibiotic resistant bacteria. A 14% global prevalence of Extended-Spectrum Beta-lactamase (ESBL) carriage was recently reported, but large intra-and interregional variations were observed. We investigated the faecal carriage rates of ESBL-, AmpC-producing and ciprofloxacin non-susceptible Escherichia coli and Klebsiella spp. in healthy Norwegians. METHODS: Rectal samples were obtained from 284 volunteers, together with demographic data and information on recent travel history. The rectal samples were screened by selective plating and E. coli and Klebsiella spp. identified using MALDI-TOF. Phenotypic and molecular characterization of resistant isolates was also performed. RESULTS: ESBL- or AmpC-producing E. coli and Klebsiella spp. were isolated from 4.9% and 3.2% of the study population, respectively. Carriage of ciprofloxacin non-susceptible isolates was detected in 9.9% of the volunteers. Molecular typing of ESBL/plasmid-mediated AmpC (pAmpC)-producing isolates suggested an allodemic situation rather than the dissemination of a specific clone in the Norwegian community. In concurrence with previous findings, travel to South-East Asia was associated with increased risk of carrying resistant E. coli or Klebsiella spp., highlighting the contribution of factors such as increased global mobility in erasing the boundaries between healthcare and community settings when it comes to spread of resistant bacteria. CONCLUSIONS: Overall, our study recognised Norway as a low-incidence country for faecal carriage of resistant bacteria among healthy individuals. Furthermore, our work denoted the importance of healthy humans as a reservoir for transmission of antibiotic resistant E. coli and Klebsiella spp.

Bioluminescence based biosensors for quantitative detection of enterococcal peptide-pheromone activity reveal inter-strain telesensing in vivo during polymicrobial systemic infection.

Enterococcus faecalis is a significant threat in the nosocomial setting due to the emergence of isolates that are multi-antibiotic resistant, refractory to the available therapies and equipped with a variety of pathogenicity determinants. This bacterium uses quorum-sensing systems to regulate its physiological processes, including the expression of virulence traits, to adapt and proliferate within a host. Here, we describe the construction and application of two bioluminescence-based reporter systems for the direct detection of the quorum-sensing regulated expression of (i) the gelatinase biosynthesis-activating pheromone (GBAP) and (ii) the cytolysin small subunit (CylL(S)) in natural samples. The two E. faecalis reporters conditionally expressed bioluminescence in the presence of GBAP and CylL(S) both in the supernatants of liquid cultures and in an agar-overlay assay in as little as three hours, with a high level of sensitivity. Biosensors employed to investigate the interaction between the fsr and cyl systems revealed that fsr impeded CylL(S) activity by 75%. Furthermore, we identified a clinical E. faecalis isolate that acted as a biological cheater, producing cytolysin only upon sensing CylL(S)-producers in its environment. This isolate enhanced its virulence during polymicrobial systemic infection of Galleria mellonella.

Transcriptomic and functional analysis of NaCl-induced stress in Enterococcus faecalis.

The robust physiology of Enterococcus faecalis facilitates tolerance to various stresses. We here report the transcriptional response of E. faecalis V583 to growth in the presence of 6.5% NaCl. Among the early responses observed was an immediate down-regulation of mscL, accompanied by an up-regulation of genes predicted to be involved in uptake of extracellular potassium and glycine betaine. The high NaCl concentration also induced expression of chaperons and cell envelope related traits, such as the enterococcal polysaccharide antigen (epa) locus. Functional genetic analysis revealed reduced salt stress resistance in both epaB and epaE mutants. The reduced salt resistance phenotype associated with the epaB mutant was restored by complementation, hence demonstrating a role of Epa in the physiological robustness of E. faecalis. Furthermore, we demonstrate that Epa confers increased resistance towards multiple cell envelope stress-inducing factors. Accordingly, these findings delineate a potential link between the robust nature of E. faecalis and its ability to perform as a human pathogen, and provide a new perspective on the mechanisms by which Epa contributes to virulence. Notably, the high NaCl concentration also resulted in strict repression of the gelE-sprE operon and impaired gelatinase activity. We demonstrate that NaCl antagonize the GBAP-pheromone dependent induction in a concentration dependent manner.

Comparative genomic analysis of pathogenic and probiotic Enterococcus faecalis isolates, and their transcriptional responses to growth in human urine.

Urinary tract infection (UTI) is the most common infection caused by enterococci, and Enterococcus faecalis accounts for the majority of enterococcal infections. Although a number of virulence related traits have been established, no comprehensive genomic or transcriptomic studies have been conducted to investigate how to distinguish pathogenic from non-pathogenic E. faecalis in their ability to cause UTI. In order to identify potential genetic traits or gene regulatory features that distinguish pathogenic from non-pathogenic E. faecalis with respect to UTI, we have performed comparative genomic analysis, and investigated growth capacity and transcriptome profiling in human urine in vitro. Six strains of different origins were cultivated and all grew readily in human urine. The three strains chosen for transcriptional analysis showed an overall similar response with respect to energy and nitrogen metabolism, stress mechanism, cell envelope modifications, and trace metal acquisition. Our results suggest that citrate and aspartate are significant for growth of E. faecalis in human urine, and manganese appear to be a limiting factor. The majority of virulence factors were either not differentially regulated or down-regulated. Notably, a significant up-regulation of genes involved in biofilm formation was observed. Strains from different origins have similar capacity to grow in human urine. The overall similar transcriptional responses between the two pathogenic and the probiotic strain suggest that the pathogenic potential of a certain E. faecalis strain may to a great extent be determined by presence of fitness and virulence factors, rather than the level of expression of such traits.

Transcriptional responses of Enterococcus faecalis V583 to bovine bile and sodium dodecyl sulfate.

Resistance to bile is a prerequisite property of the gastrointestinal bacterial flora. Bile acids are powerful detergents, and resistance to sodium dodecyl sulfate (SDS) has therefore often been considered relevant to studies of bile resistance. We have studied the effects of bovine bile (BB) and SDS on Enterococcus faecalis V583 by traditional growth studies and microarrays. Transcriptional responses were studied by time course experiments. In the presence of BB (V583-BB) or SDS (V583-SDS), 308 and 209 genes were identified as differentially expressed at one or more time points, respectively. In V583 treated with both BB and SDS (V583-BB-SDS), 254 genes showed differential expression. Detergents exert their toxic effects primarily on the microbial membrane. The enrichment of differentially transcribed genes that encode proteins with membrane-associated functions and/or locations indicates a major impact of all three treatments on the integrity and functionality of the cell membrane. Two gene clusters involved in fatty acid biosynthesis were repressed in V583-BB and V583-BB-SDS and partly induced in V583-SDS. Furthermore, two EmrB/QacA family drug resistance transporters and a vacuolar-type ATPase were induced in V583-BB and V583-BB-SDS. None of the putative bile salt hydrolase homologs in V583 showed differential expression during the bile treatments. The transcriptional profile of V583-BB-SDS was qualitatively more similar to the response in V583-BB than to that in V583-SDS, suggesting that the presence of bile suppresses the effects of SDS in V583-BB-SDS. The overall results presented here indicate that different mechanisms are involved in detergent resistance in E. faecalis.