The Genome of the Plant-Associated Lactic Acid Bacterium Lactococcus lactis KF147 Harbors a Hybrid NRPS-PKS System Conserved in Strains of the Dental Cariogenic Streptococcus mutans.

Lactococcus lactis subsp. lactis KF147 as a non-dairy strain from lactic acid bacteria (LAB) can inhabit plant tissues. It can grow on complex carbohydrates derived from plant cell walls. Its genome size is one of the largest among the sequenced lactococcal strains, possessing many genes that do not have homologues in the published genome sequences of dairy-associated L. lactis strains. In silico analysis has identified a gene cluster encoding a hybrid NRPS-PKS system (composed of non-ribosomal peptide synthetases and polyketide synthases) in the L. lactis KF147 genome, as first example of a LAB possessing such hybrid mega-enzymes. Hybrid systems produce hybrid NRP-PK secondary metabolites (natural products) in a wide variety of bacteria, fungi, and plants. In the hybrid NRPS-PKS system of L. lactis KF147, a total of 21 NRPS and 8 PKS domains were identified that are arranged into 6 NRPS modules, 3 PKS modules, and two single functional domains (trans-acyl-transferase "transAT" and thioesterase). We found homologous hybrid systems having similar gene, module, and domain organization in six other L. lactis strains and 25 strains of the dental cariogenic Streptococcus mutans. This study mainly aimed to predict the structure and function of the hybrid NRP-PK product of L. lactis KF147 using comparative genomics techniques, and included a detailed analysis of the regulatory system. Various bioinformatical approaches were used to predict the substrate specificity of the six A domains and the iterative transAT domain. Functional conservation of the A domains within different-niche-associated strains supported the prediction of the primary core structure of the putative hybrid natural product to be Leu-DLeu-Asp-DAsn-Gly-MC-MC-MC-DAsp (MC = Malonyl-CoA). Oxidative stress resistance and biofilm formation are the most probable functions of this hybrid system. The need for such a system in two different niches is argued, as an adaptation of L. lactis and S. mutans to adhere to plant tissues and human teeth, respectively, in an oxidative environment.

Lactobacillus plantarum Strains Can Enhance Human Mucosal and Systemic Immunity and Prevent Non-steroidal Anti-inflammatory Drug Induced Reduction in T Regulatory Cells.

Orally ingested bacteria interact with intestinal mucosa and may impact immunity. However, insights in mechanisms involved are limited. In this randomized placebo-controlled cross-over trial, healthy human subjects were given Lactobacillus plantarum supplementation (strain TIFN101, CIP104448, or WCFS1) or placebo for 7 days. To determine whether L. plantarum can enhance immune response, we compared the effects of three stains on systemic and gut mucosal immunity, by among others assessing memory responses against tetanus toxoid (TT)-antigen, and mucosal gene transcription, in human volunteers during induction of mild immune stressor in the intestine, by giving a commonly used enteropathic drug, indomethacin [non-steroidal anti-inflammatory drug (NSAID)]. Systemic effects of the interventions were studies in peripheral blood samples. NSAID was found to induce a reduction in serum CD4+/Foxp3 regulatory cells, which was prevented by L. plantarum TIFN101. T-cell polarization experiments showed L. plantarum TIFN101 to enhance responses against TT-antigen, which indicates stimulation of memory responses by this strain. Cell extracts of the specific L. plantarum strains provoked responses after WCFS1 and TIFN101 consumption, indicating stimulation of immune responses against the specific bacteria. Mucosal immunomodulatory effects were studied in duodenal biopsies. In small intestinal mucosa, TIFN101 upregulated genes associated with maintenance of T- and B-cell function and antigen presentation. Furthermore, L. plantarum TIFN101 and WCFS1 downregulated immunological pathways involved in antigen presentation and shared downregulation of snoRNAs, which may suggest cellular destabilization, but may also be an indicator of tissue repair. Full sequencing of the L. plantarum strains revealed possible gene clusters that might be responsible for the differential biological effects of the bacteria on host immunity. In conclusion, the impact of oral consumption L. plantarum on host immunity is strain dependent and involves responses against bacterial cell components. Some strains may enhance specific responses against pathogens by enhancing antigen presentation and leukocyte maintenance in mucosa. In future studies and clinical settings, caution should be taken in selecting beneficial bacteria as closely related strains can have different effects. Our data show that specific bacterial strains can prevent immune stress induced by commonly consumed painkillers such as NSAID and can have enhancing beneficial effects on immunity of consumers by stimulating antigen presentation and memory responses.

Draft Genome Sequences of 24 Lactococcus lactis Strains.

The lactic acid bacterium Lactococcus lactis is widely used for the production of fermented dairy products. Here, we present the draft genome sequences of 24 L. lactis strains isolated from different environments and geographic locations.

Draft Genome Sequences of 11 Lactococcus lactis subsp. cremoris Strains.

The lactic acid bacterium Lactococcus lactis is widely used for the fermentation of dairy products. Here, we present the draft genome sequences of 11 L. lactis subsp. cremoris strains isolated from different environments.

Draft Genome Sequence of Lactobacillus delbrueckii subsp. bulgaricus LBB.B5.

Lactobacillus delbrueckii subsp. bulgaricus LBB.B5 originates from homemade Bulgarian yogurt and was selected for its ability to form a strong association with Streptococcus thermophilus The genome sequence will facilitate elucidating the genetic background behind the contribution of LBB.B5 to the taste and aroma of yogurt and its exceptional protocooperation with S. thermophilus.

Linking Bacillus cereus Genotypes and Carbohydrate Utilization Capacity.

We characterised carbohydrate utilisation of 20 newly sequenced Bacillus cereus strains isolated from food products and food processing environments and two laboratory strains, B. cereus ATCC 10987 and B. cereus ATCC 14579. Subsequently, genome sequences of these strains were analysed together with 11 additional B. cereus reference genomes to provide an overview of the different types of carbohydrate transporters and utilization systems found in B. cereus strains. The combined application of API tests, defined growth media experiments and comparative genomics enabled us to link the carbohydrate utilisation capacity of 22 B. cereus strains with their genome content and in some cases to the panC phylogenetic grouping. A core set of carbohydrates including glucose, fructose, maltose, trehalose, N-acetyl-glucosamine, and ribose could be used by all strains, whereas utilisation of other carbohydrates like xylose, galactose, and lactose, and typical host-derived carbohydrates such as fucose, mannose, N-acetyl-galactosamine and inositol is limited to a subset of strains. Finally, the roles of selected carbohydrate transporters and utilisation systems in specific niches such as soil, foods and the human host are discussed.

Draft Genome Sequence of Streptococcus thermophilus C106, a Dairy Isolate from an Artisanal Cheese Produced in the Countryside of Ireland.

The lactic acid bacterium Streptococcus thermophilus is widely used for the fermentation of dairy products. Here, we present the draft genome sequence of S. thermophilus C106 isolated from an artisanal cheese produced in the countryside of Ireland.

A Novel Quality Measure and Correction Procedure for the Annotation of Microbial Translation Initiation Sites.

The identification of translation initiation sites (TISs) constitutes an important aspect of sequence-based genome analysis. An erroneous TIS annotation can impair the identification of regulatory elements and N-terminal signal peptides, and also may flaw the determination of descent, for any particular gene. We have formulated a reference-free method to score the TIS annotation quality. The method is based on a comparison of the observed and expected distribution of all TISs in a particular genome given prior gene-calling. We have assessed the TIS annotations for all available NCBI RefSeq microbial genomes and found that approximately 87% is of appropriate quality, whereas 13% needs substantial improvement. We have analyzed a number of factors that could affect TIS annotation quality such as GC-content, taxonomy, the fraction of genes with a Shine-Dalgarno sequence and the year of publication. The analysis showed that only the first factor has a clear effect. We have then formulated a straightforward Principle Component Analysis-based TIS identification strategy to self-organize and score potential TISs. The strategy is independent of reference data and a priori calculations. A representative set of 277 genomes was subjected to the analysis and we found a clear increase in TIS annotation quality for the genomes with a low quality score. The PCA-based annotation was also compared with annotation with the current tool of reference, Prodigal. The comparison for the model genome of Escherichia coli K12 showed that both methods supplement each other and that prediction agreement can be used as an indicator of a correct TIS annotation. Importantly, the data suggest that the addition of a PCA-based strategy to a Prodigal prediction can be used to 'flag' TIS annotations for re-evaluation and in addition can be used to evaluate a given annotation in case a Prodigal annotation is lacking.

CiVi: circular genome visualization with unique features to analyze sequence elements.

UNLABELLED: We have developed CiVi, a user-friendly web-based tool to create custom circular maps to aid the analysis of microbial genomes and sequence elements. Sequence related data such as gene-name, COG class, PFAM domain, GC%, and subcellular location can be comprehensively viewed. Quantitative gene-related data (e.g. expression ratios or read counts) as well as predicted sequence elements (e.g. regulatory sequences) can be uploaded and visualized. CiVi accommodates the analysis of genomic elements by allowing a visual interpretation in the context of: (i) their genome-wide distribution, (ii) provided experimental data and (iii) the local orientation and location with respect to neighboring genes. CiVi thus enables both experts and non-experts to conveniently integrate public genome data with the results of genome analyses in circular genome maps suitable for publication. CONTACT: L.Overmars@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. AVAILABILITY AND IMPLEMENTATION: CiVi is freely available at http://civi.cmbi.ru.nl.

Draft Genome Sequence of Lactobacillus plantarum CMPG5300, a Human Vaginal Isolate.

The draft genome of a highly auto-aggregating Lactobacillus plantarum strain isolated from a human vagina is reported. The peculiar phenotype also provides an adhesive and co-aggregative potential with various pathogens, which could be of significance in the vaginal niche. Detailed genome analysis could aid in identifying the adhesins of the strain.

Safety of the surrogate microorganism Enterococcus faecium NRRL B-2354 for use in thermal process validation.

Enterococcus faecium NRRL B-2354 is a surrogate microorganism used in place of pathogens for validation of thermal processing technologies and systems. We evaluated the safety of strain NRRL B-2354 based on its genomic and functional characteristics. The genome of E. faecium NRRL B-2354 was sequenced and found to comprise a 2,635,572-bp chromosome and a 214,319-bp megaplasmid. A total of 2,639 coding sequences were identified, including 45 genes unique to this strain. Hierarchical clustering of the NRRL B-2354 genome with 126 other E. faecium genomes as well as pbp5 locus comparisons and multilocus sequence typing (MLST) showed that the genotype of this strain is most similar to commensal, or community-associated, strains of this species. E. faecium NRRL B-2354 lacks antibiotic resistance genes, and both NRRL B-2354 and its clonal relative ATCC 8459 are sensitive to clinically relevant antibiotics. This organism also lacks, or contains nonfunctional copies of, enterococcal virulence genes including acm, cyl, the ebp operon, esp, gelE, hyl, IS16, and associated phenotypes. It does contain scm, sagA, efaA, and pilA, although either these genes were not expressed or their roles in enterococcal virulence are not well understood. Compared with the clinical strains TX0082 and 1,231,502, E. faecium NRRL B-2354 was more resistant to acidic conditions (pH 2.4) and high temperatures (60°C) and was able to grow in 8% ethanol. These findings support the continued use of E. faecium NRRL B-2354 in thermal process validation of food products.

Combining chemoinformatics with bioinformatics: in silico prediction of bacterial flavor-forming pathways by a chemical systems biology approach "reverse pathway engineering".

The incompleteness of genome-scale metabolic models is a major bottleneck for systems biology approaches, which are based on large numbers of metabolites as identified and quantified by metabolomics. Many of the revealed secondary metabolites and/or their derivatives, such as flavor compounds, are non-essential in metabolism, and many of their synthesis pathways are unknown. In this study, we describe a novel approach, Reverse Pathway Engineering (RPE), which combines chemoinformatics and bioinformatics analyses, to predict the "missing links" between compounds of interest and their possible metabolic precursors by providing plausible chemical and/or enzymatic reactions. We demonstrate the added-value of the approach by using flavor-forming pathways in lactic acid bacteria (LAB) as an example. Established metabolic routes leading to the formation of flavor compounds from leucine were successfully replicated. Novel reactions involved in flavor formation, i.e. the conversion of alpha-hydroxy-isocaproate to 3-methylbutanoic acid and the synthesis of dimethyl sulfide, as well as the involved enzymes were successfully predicted. These new insights into the flavor-formation mechanisms in LAB can have a significant impact on improving the control of aroma formation in fermented food products. Since the input reaction databases and compounds are highly flexible, the RPE approach can be easily extended to a broad spectrum of applications, amongst others health/disease biomarker discovery as well as synthetic biology.

Lactobacillus paracasei comparative genomics: towards species pan-genome definition and exploitation of diversity.

Lactobacillus paracasei is a member of the normal human and animal gut microbiota and is used extensively in the food industry in starter cultures for dairy products or as probiotics. With the development of low-cost, high-throughput sequencing techniques it has become feasible to sequence many different strains of one species and to determine its "pan-genome". We have sequenced the genomes of 34 different L. paracasei strains, and performed a comparative genomics analysis. We analysed genome synteny and content, focussing on the pan-genome, core genome and variable genome. Each genome was shown to contain around 2800-3100 protein-coding genes, and comparative analysis identified over 4200 ortholog groups that comprise the pan-genome of this species, of which about 1800 ortholog groups make up the conserved core. Several factors previously associated with host-microbe interactions such as pili, cell-envelope proteinase, hydrolases p40 and p75 or the capacity to produce short branched-chain fatty acids (bkd operon) are part of the L. paracasei core genome present in all analysed strains. The variome consists mainly of hypothetical proteins, phages, plasmids, transposon/conjugative elements, and known functions such as sugar metabolism, cell-surface proteins, transporters, CRISPR-associated proteins, and EPS biosynthesis proteins. An enormous variety and variability of sugar utilization gene cassettes were identified, with each strain harbouring between 25-53 cassettes, reflecting the high adaptability of L. paracasei to different niches. A phylogenomic tree was constructed based on total genome contents, and together with an analysis of horizontal gene transfer events we conclude that evolution of these L. paracasei strains is complex and not always related to niche adaptation. The results of this genome content comparison was used, together with high-throughput growth experiments on various carbohydrates, to perform gene-trait matching analysis, in order to link the distribution pattern of a specific phenotype to the presence/absence of specific sets of genes.

Classification of the adenylation and acyl-transferase activity of NRPS and PKS systems using ensembles of substrate specific hidden Markov models.

There is a growing interest in the Non-ribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs) of microbes, fungi and plants because they can produce bioactive peptides such as antibiotics. The ability to identify the substrate specificity of the enzyme's adenylation (A) and acyl-transferase (AT) domains is essential to rationally deduce or engineer new products. We here report on a Hidden Markov Model (HMM)-based ensemble method to predict the substrate specificity at high quality. We collected a new reference set of experimentally validated sequences. An initial classification based on alignment and Neighbor Joining was performed in line with most of the previously published prediction methods. We then created and tested single substrate specific HMMs and found that their use improved the correct identification significantly for A as well as for AT domains. A major advantage of the use of HMMs is that it abolishes the dependency on multiple sequence alignment and residue selection that is hampering the alignment-based clustering methods. Using our models we obtained a high prediction quality for the substrate specificity of the A domains similar to two recently published tools that make use of HMMs or Support Vector Machines (NRPSsp and NRPS predictor2, respectively). Moreover, replacement of the single substrate specific HMMs by ensembles of models caused a clear increase in prediction quality. We argue that the superiority of the ensemble over the single model is caused by the way substrate specificity evolves for the studied systems. It is likely that this also holds true for other protein domains. The ensemble predictor has been implemented in a simple web-based tool that is available at http://www.cmbi.ru.nl/NRPS-PKS-substrate-predictor/.

The highly autoaggregative and adhesive phenotype of the vaginal Lactobacillus plantarum strain CMPG5300 is sortase dependent.

Lactobacilli are important for the maintenance of a healthy ecosystem in the human vagina. Various mechanisms are postulated but so far are poorly substantiated by molecular studies, such as mutant analysis. Bacterial autoaggregation is an interesting phenomenon that can promote adhesion to host cells and displacement of pathogens. In this study, we report on the identification of a human vaginal isolate, Lactobacillus plantarum strain CMPG5300, which shows high autoaggregative and adhesive capacity. To investigate the importance of sortase-dependent proteins (SDPs) in these phenotypes, a gene deletion mutant was constructed for srtA, the gene encoding the housekeeping sortase that covalently anchors these SDPs to the cell surface. This mutant lost the capacity to autoaggregate, showed a decrease in adhesion to vaginal epithelial cells, and lost biofilm-forming capacity under the conditions tested. These results indicate that the housekeeping sortase SrtA of CMPG5300 is a key determinant of the peculiar surface properties of this vaginal Lactobacillus strain.

MGcV: the microbial genomic context viewer for comparative genome analysis.

BACKGROUND: Conserved gene context is used in many types of comparative genome analyses. It is used to provide leads on gene function, to guide the discovery of regulatory sequences, but also to aid in the reconstruction of metabolic networks. We present the Microbial Genomic context Viewer (MGcV), an interactive, web-based application tailored to strengthen the practice of manual comparative genome context analysis for bacteria. RESULTS: MGcV is a versatile, easy-to-use tool that renders a visualization of the genomic context of any set of selected genes, genes within a phylogenetic tree, genomic segments, or regulatory elements. It is tailored to facilitate laborious tasks such as the interactive annotation of gene function, the discovery of regulatory elements, or the sequence-based reconstruction of gene regulatory networks. We illustrate that MGcV can be used in gene function annotation by visually integrating information on prokaryotic genes, like their annotation as available from NCBI with other annotation data such as Pfam domains, sub-cellular location predictions and gene-sequence characteristics such as GC content. We also illustrate the usefulness of the interactive features that allow the graphical selection of genes to facilitate data gathering (e.g. upstream regions, ID's or annotation), in the analysis and reconstruction of transcription regulation. Moreover, putative regulatory elements and their corresponding scores or data from RNA-seq and microarray experiments can be uploaded, visualized and interpreted in (ranked-) comparative context maps. The ranked maps allow the interpretation of predicted regulatory elements and experimental data in light of each other. CONCLUSION: MGcV advances the manual comparative analysis of genes and regulatory elements by providing fast and flexible integration of gene related data combined with straightforward data retrieval. MGcV is available at http://mgcv.cmbi.ru.nl.

Genotype-phenotype matching analysis of 38 Lactococcus lactis strains using random forest methods.

BACKGROUND: Lactococcus lactis is used in dairy food fermentation and for the efficient production of industrially relevant enzymes. The genome content and different phenotypes have been determined for multiple L. lactis strains in order to understand intra-species genotype and phenotype diversity and annotate gene functions. In this study, we identified relations between gene presence and a collection of 207 phenotypes across 38 L. lactis strains of dairy and plant origin. Gene occurrence and phenotype data were used in an iterative gene selection procedure, based on the Random Forest algorithm, to identify genotype-phenotype relations. RESULTS: A total of 1388 gene-phenotype relations were found, of which some confirmed known gene-phenotype relations, such as the importance of arabinose utilization genes only for strains of plant origin. We also identified a gene cluster related to growth on melibiose, a plant disaccharide; this cluster is present only in melibiose-positive strains and can be used as a genetic marker in trait improvement. Additionally, several novel gene-phenotype relations were uncovered, for instance, genes related to arsenite resistance or arginine metabolism. CONCLUSIONS: Our results indicate that genotype-phenotype matching by integrating large data sets provides the possibility to identify gene-phenotype relations, possibly improve gene function annotation and identified relations can be used for screening bacterial culture collections for desired phenotypes. In addition to all gene-phenotype relations, we also provide coherent phenotype data for 38 Lactococcus strains assessed in 207 different phenotyping experiments, which to our knowledge is the largest to date for the Lactococcus lactis species.

Genome sequence of Lactobacillus pentosus KCA1: vaginal isolate from a healthy premenopausal woman.

The vaginal microbiota, in particular Lactobacillus species, play an important role in female health through modulation of immunity, countering pathogens and maintaining a pH below 4.7. We report the isolation and genome sequence of Lactobacillus pentosus strain KCA1 (formally known as L. plantarum) from the vagina of a healthy Nigerian woman. The genome was sequenced using Illumina GA II technology. The resulting 16,920,226 paired-end reads were assembled with the Velvet tool. Contigs were annotated using the RAST server, and manually curated. A comparative analysis with the available genomes of L. pentosus IG1 and L. plantarum WCFS1 showed that over 15% of the predicted functional activities are found only in this strain. The strain has a chromosome sequence of 3,418,159 bp with a G+C content of 46.4%, and is devoid of plasmids. Novel gene clusters or variants of known genes relative to the reference genomes were found. In particular, the strain has loci encoding additional putative mannose phosphotransferase systems. Clusters of genes include those for utilization of hydantoin, isopropylmalate, malonate, rhamnosides, and genes for assimilation of polyglycans, suggesting the metabolic versatility of L. pentosus KCA1. Loci encoding putative phage defense systems were also found including clustered regularly interspaced short palindromic repeats (CRISPRs), abortive infection (Abi) systems and toxin-antitoxin systems (TA). A putative cluster of genes for biosynthesis of a cyclic bacteriocin precursor, here designated as pentocin KCA1 (penA) were identified. These findings add crucial information for understanding the genomic and geographic diversity of vaginal lactobacilli.

Complete Genome Sequence of Anoxybacillus flavithermus TNO-09.006, a Thermophilic Sporeformer Associated with a Dairy-Processing Environment.

Spores of thermophilic spore-forming bacteria are a common cause of contamination in dairy products. We isolated the thermophilic strain Anoxybacillus flavithermus TNO-09.006 from a milk-processing plant, and we report the complete genome of this isolate consisting of a single chromosome of 2.65 Mb.

Inferring patient to patient transmission of Mycobacterium tuberculosis from whole genome sequencing data.

BACKGROUND: Mycobacterium tuberculosis is characterised by limited genomic diversity, which makes the application of whole genome sequencing particularly attractive for clinical and epidemiological investigation. However, in order to confidently infer transmission events, an accurate knowledge of the rate of change in the genome over relevant timescales is required. METHODS: We attempted to estimate a molecular clock by sequencing 199 isolates from epidemiologically linked tuberculosis cases, collected in the Netherlands spanning almost 16 years. RESULTS: Multiple analyses support an average mutation rate of ~0.3 SNPs per genome per year. However, all analyses revealed a very high degree of variation around this mean, making the confirmation of links proposed by epidemiology, and inference of novel links, difficult. Despite this, in some cases, the phylogenetic context of other strains provided evidence supporting the confident exclusion of previously inferred epidemiological links. CONCLUSIONS: This in-depth analysis of the molecular clock revealed that it is slow and variable over short time scales, which limits its usefulness in transmission studies. However, the superior resolution of whole genome sequencing can provide the phylogenetic context to allow the confident exclusion of possible transmission events previously inferred via traditional DNA fingerprinting techniques and epidemiological cluster investigation. Despite the slow generation of variation even at the whole genome level we conclude that the investigation of tuberculosis transmission will benefit greatly from routine whole genome sequencing.