Comparative genomics of Lactobacillus sakei with emphasis on strains from meat.

Lactobacillus sakei is a lactic acid bacterium important in food microbiology mainly due to its ability to ferment and preserve meat. The genome sequence of L. sakei strain 23K has revealed specialized metabolic capacities that reflect the bacterium's adaption to meat products, and that differentiate it from other LAB. An extensive genomic diversity analysis was conducted to elucidate the core features of the species, and to provide a better comprehension of niche adaptation of the organism. Here, we describe the genomic comparison of 18 strains of L. sakei originating mainly from processed meat against the 23K strain by comparative genome hybridization. Pulsed field gel electrophoresis was used to estimate the genome sizes of the strains, which varied from 1.880 to 2.175 Mb, and the 23K genome was among the smallest. Consequently, a large part of the genome of this strain belongs to a common gene pool invariant in this species. The majority of genes important in adaption to meat products, the ability to flexibly use meat components, and robustness during meat processing and storage were conserved, such as genes involved in nucleoside scavenging, catabolism of arginine, and the ability to cope with changing redox and oxygen levels, which is indicative of the role these genes play in niche specialization within the L. sakei species. Moreover, an additional set of sequenced L. sakei genes beyond the 23K genome was present on the microarray used, and it was demonstrated that all the strains carry remnants of or complete bacteriocin operons. The genomic divergence corresponded mainly to five regions in the 23K genome, which showed features consistent with horizontal gene transfer. Carbohydrate-fermentation profiles of the strains were evaluated in light of the CGH data, and for most substrates, the genotypes were consistent with the phenotypes. We have demonstrated a highly conserved organization of the L. sakei genomes investigated, and the 23K strain is a suitable model organism to study core features of the L. sakei species.

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.

Detection of divergent genes in microbial aCGH experiments.

BACKGROUND: Array-based comparative genome hybridization (aCGH) is a tool for rapid comparison of genomes from different bacterial strains. The purpose of such analysis is to detect highly divergent or absent genes in a sample strain compared to an index strain. Development of methods for analyzing aCGH data has primarily focused on copy number abberations in cancer research. In microbial aCGH analyses, genes are typically ranked by log-ratios, and classification into divergent or present is done by choosing a cutoff log-ratio, either manually or by statistics calculated from the log-ratio distribution. As experimental settings vary considerably, it is not possible to develop a classical discriminant or statistical learning approach. METHODS: We introduce a more efficient method for analyzing microbial aCGH data using a finite mixture model and a data rotation scheme. Using the average posterior probabilities from the model fitted to log-ratios before and after rotation, we get a score for each gene, and demonstrate its advantages for ranking and detecting divergent genes with enlarged specificity and sensitivity. RESULTS: The procedure is tested and compared to other approaches on simulated data sets, as well as on four experimental validation data sets for aCGH analysis on fully sequenced strains of Staphylococcus aureus and Streptococcus pneumoniae. CONCLUSION: When tested on simulated data as well as on four different experimental validation data sets from experiments with only fully sequenced strains, our procedure out-competes the standard procedures of using a simple log-ratio cutoff for classification into present and divergent genes.

The Response of Enterococcus faecalis V583 to Chloramphenicol Treatment.

Many Enterococcus faecalis strains display tolerance or resistance to many antibiotics, but genes that contribute to the resistance cannot be specified. The multiresistant E. faecalis V583, for which the complete genome sequence is available, survives and grows in media containing relatively high levels of chloramphenicol. No specific genes coding for chloramphenicol resistance has been recognized in V583. We used microarrays to identify genes and mechanisms behind the tolerance to chloramphenicol in V583, by comparison of cells treated with subinhibitory concentrations of chloramphenicol and untreated V583 cells. During a time course experiment, more than 600 genes were significantly differentially transcribed. Since chloramphenicol affects protein synthesis in bacteria, many genes involved in protein synthesis, for example, genes for ribosomal proteins, were induced. Genes involved in amino acid biosynthesis, for example, genes for tRNA synthetases and energy metabolism were downregulated, mainly. Among the upregulated genes were EF1732 and EF1733, which code for potential chloramphenicol transporters. Efflux of drug out of the cells may be one mechanism used by V583 to overcome the effect of chloramphenicol.

Mixture models as a method to find present and divergent genes in comparative genomic hybridization studies on bacteria.

Comparative genomic hybridization (CGH) using microarrays is performed on bacteria in order to test for genomic diversity within various bacterial species. The microarrays used for CGH are based on the genome of a fully sequenced bacterium strain, denoted reference strain. Labelled DNA fragments from a sample strain of interest and from the reference strain are hybridized to the array. Based on the obtained ratio intensities and the total intensities of the signals, each gene is classified as either present (one copy or multiple copies) or divergent (zero copies). In this paper mixture models with different number of components are tted on different combinations of variables and compared with each other. The study shows that mixture models fitted on both the ratio intensities and the total intensities including the replicates for each gene improve, compared to previously published methods, the results for several of the data sets tested. Some summaries of the data sets are proposed as a guide for the choice of model and the choice of number of components. The models are applied on data from CGH experiments with the bacteria Staphylococcus aureus and

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.

Survey of genomic diversity among Enterococcus faecalis strains by microarray-based comparative genomic hybridization.

We have compared nine Enterococcus faecalis strains with E. faecalis V583 by comparative genomic hybridization using microarrays (CGH). The strains used in this study (the "test" strains) originated from various environments. CGH is a powerful and promising tool for obtaining novel information on genome diversity in bacteria. By CGH, one obtains clues about which genes are present or divergent in the strains, compared to a reference strain (here, V583). The information obtained by CGH is important from both ecological and systematic points of view. CGH of E. faecalis showed considerable diversity in gene content: Compared to V583, the percentage of divergent genes in the test strains varied from 15% to 23%, and 154 genes were divergent in all strains. The main variation was found in regions corresponding to exogenously acquired or mobile DNA in V583. Antibiotic resistance genes, virulence factors, and integrated plasmid genes dominated among the divergent genes. The strains examined showed various contents of genes corresponding to the pTEF1, pTEF2, and pTEF3 genes in V583. The extensive transport and metabolic capabilities of V583 appeared similar in the test strains; CGH indicated that the ability to transport and metabolize various carbohydrates was similar in the test strains (verified by API 50 CH assays). The contents of genes related to stress tolerance appeared similar in V583 and the nine test strains, supporting the view of E. faecalis as an organism able to resist harsh conditions.

Regression as a method to predict copy numbers in comparative genomic hybridization studies on bacteria.

Comparative genomic hybridizations (CGH) using microarrays are performed with bacteria in order to determine the level of genomic similarity between various strains. The microarrays applied in CGH experiments are constructed on the basis of the genome sequence of one strain, which is used as a control, or reference, in each experiment. A strain being compared with the known strain is called the unknown strain. The ratios of fluorescent intensities obtained from the spots on the microarrays can be used to determine which genes are divergent in the unknown strain, as well as to predict the copy number of actual genes in the unknown strain. In this paper, we focus on the prediction of gene copy number based on data from CGH experiments. We assumed a linear connection between the log2 of the copy number and the observed log2-ratios, then predictors based on the factor analysis model and the linear random model were proposed in an attempt to identify the copy numbers. These predictors were compared to using the ratio of the intensities directly. Simulations indicated that the proposed predictors improved the prediction of the copy number in most situations. The predictors were applied on CGH data obtained from experiments with Enterococcus faecalis strains in order to determine copy number of relevant genes in five different strains.

Complete sequence of the enterocin Q-encoding plasmid pCIZ2 from the multiple bacteriocin producer Enterococcus faecium L50 and genetic characterization of enterocin Q production and immunity.

The locations of the genetic determinants for enterocin L50 (EntL50A and EntL50B), enterocin Q (EntQ), and enterocin P (EntP) in the multiple bacteriocin producer Enterococcus faecium strain L50 were determined. These bacteriocin genes occur at different locations; entL50AB (encoding EntL50A and EntL50B) are on the 50-kb plasmid pCIZ1, entqA (encoding EntQ) is on the 7.4-kb plasmid pCIZ2, and entP (encoding EntP) is on the chromosome. The complete nucleotide sequence of pCIZ2 was determined to be 7,383 bp long and contains 10 putative open reading frames (ORFs) organized in three distinct regions. The first region contains three ORFs: entqA preceded by two divergently oriented genes, entqB and entqC. EntqB shows high levels of similarity to bacterial ATP-binding cassette (ABC) transporters, while EntqC displays no significant similarity to any known protein. The second region encompasses four ORFs (orf4 to orf7), and ORF4 and ORF5 display high levels of similarity to mobilization proteins from E. faecium and Enterococcus faecalis. In addition, features resembling a transfer origin region (oriT) were found in the promoter area of orf4. The third region contains three ORFs (orf8 to orf10), and ORF8 and ORF9 exhibit similarity to the replication initiator protein RepE from E. faecalis and to RepB proteins, respectively. To clarify the minimum requirement for EntQ synthesis, we subcloned and heterologously expressed a 2,371-bp fragment from pCIZ2 that encompasses only the entqA, entqB, and entqC genes in Lactobacillus sakei, and we demonstrated that this fragment is sufficient for EntQ production. Moreover, we also obtained experimental results indicating that EntqB is involved in ABC transporter-mediated EntQ secretion, while EntqC confers immunity to this bacteriocin.

Transcriptional response of Enterococcus faecalis V583 to erythromycin.

A transcriptional profile of Enterococcus faecalis V583 (V583) treated with erythromycin is presented. This is the first study describing a complete transcriptional profile of Enterococcus. E. faecalis is a common and nonvirulent bacterium in many natural environments, but also an important cause of nosocomial infections. We have used a genome-wide microarray based on the genome sequence of V583 to study gene expression in cells exposed to erythromycin. V583 is resistant to relatively high concentrations of erythromycin, but growth is retarded by the treatment. The effect of erythromycin treatment on V583 was studied by a time course experiment; samples were extracted at five time points over a period of 90 min. A drastic change in gene transcription was seen with the erythromycin-treated cells compared to the untreated cells. Altogether, 260 genes were down-regulated at one or more time points, while 340 genes were up-regulated. Genes encoding hypothetical proteins and genes encoding transport and binding proteins were the two most dominating groups of differentially expressed genes. The gene encoding ermB (EFA0007) was expressed, but not differentially, which indicated that other genes are important for the survival and growth maintenance of V583 treated with erythromycin. One of these genes is a putative MsrC-like protein, which was up-regulated at all time points studied. Other specific genes that were found to be up-regulated were genes encoding ABC transporters and two-component regulatory systems, and these may be genes that are important for the specific response of V583 to erythromycin.

Phylogeny and functional expression of ribulose 1,5-bisphosphate carboxylase/oxygenase from the autotrophic ammonia-oxidizing bacterium Nitrosospira sp. isolate 40KI.

The autotrophic ammonia-oxidizing bacteria (AOB), which play an important role in the global nitrogen cycle, assimilate CO(2) by using ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO). Here we describe the first detailed study of RubisCO (cbb) genes and proteins from the AOB. The cbbLS genes from Nitrosospira sp. isolate 40KI were cloned and sequenced. Partial sequences of the RubisCO large subunit (CbbL) from 13 other AOB belonging to the beta and gamma subgroups of the class Proteobacteria are also presented. All except one of the beta-subgroup AOB possessed a red-like type I RubisCO with high sequence similarity to the Ralstonia eutropha enzyme. All of these new red-like RubisCOs had a unique six-amino-acid insert in CbbL. Two of the AOB, Nitrosococcus halophilus Nc4 and Nitrosomonas europaea Nm50, had a green-like RubisCO. With one exception, the phylogeny of the AOB CbbL was very similar to that of the 16S rRNA gene. The presence of a green-like RubisCO in N. europaea was surprising, as all of the other beta-subgroup AOB had red-like RubisCOs. The green-like enzyme of N. europaea Nm50 was probably acquired by horizontal gene transfer. Functional expression of Nitrosospira sp. isolate 40KI RubisCO in the chemoautotrophic host R. eutropha was demonstrated. Use of an expression vector harboring the R. eutropha cbb control region allowed regulated expression of Nitrosospira sp. isolate 40KI RubisCO in an R. eutropha cbb deletion strain. The Nitrosospira RubisCO supported autotrophic growth of R. eutropha with a doubling time of 4.6 h. This expression system may allow further functional analysis of AOB cbb genes.

Inducible bacteriocin production in Lactobacillus is regulated by differential expression of the pln operons and by two antagonizing response regulators, the activity of which is enhanced upon phosphorylation.

Expression of the five (pln) operons involved in the bacteriocin production of Lactobacillus plantarum C11 is regulated by a so-called pheromone-based signal-transducing network, in which the peptide pheromone (PlnA) induces bacteriocin production through the action of a histidine protein kinase (PlnB) and two antagonizing response regulators (PlnC as an activator and PlnD as a negative regulator). All pln-regulated promoters contain a conserved pair of direct repeats that serve as binding sites for PlnC and PlnD. In the present work, we show that the five PlnA-responsive operons are differentially expressed with regard to both timing and strength, and that the pheromone triggers a strong autoactivating loop of the regulatory unit (plnABCD) during an early stage of induction that gradually leads to enhanced activation of the other operons. The transport operon (plnGHSTUV), which is involved in the secretion of the pheromone and bacteriocins, is also expressed relatively early upon induction, but is quickly turned off soon after peak expression. Further investigation of the various promoters revealed that, although subtle differences within the promoter regions could account for the observed differential regulation, the presence of a downstream promoter-proximal sequence in one promoter was found to cause delayed peak activity. How phosphorylation regulates the activity of the pln response regulators was also accessed by direct mutagenesis at their phosphorylation sites. It was found that the two response regulators exert activity at two different levels: a low level when they are not phosphorylated and an elevated level when they are phosphorylated. The present data demonstrate that bacteriocin production in L. plantarum C11 is a highly regulated process, in which different regulatory mechanisms are applied to fine tune the timing and strength of expression of the five pln operons.