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.

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.

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.

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.

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.

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.

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.

Complete genome sequence of Geobacillus thermoglucosidans TNO-09.020, a thermophilic sporeformer associated with a dairy-processing environment.

Thermophilic spore-forming bacteria are a common cause of contamination in dairy products. We isolated the thermophilic strain Geobacillus thermoglucosidans TNO-09.020 from a milk processing plant and report the complete genome of a dairy plant isolate consisting of a single chromosome of 3.75 Mb.

Computational analysis of cysteine and methionine metabolism and its regulation in dairy starter and related bacteria.

Sulfuric volatile compounds derived from cysteine and methionine provide many dairy products with a characteristic odor and taste. To better understand and control the environmental dependencies of sulfuric volatile compound formation by the dairy starter bacteria, we have used the available genome sequence and experimental information to systematically evaluate the presence of the key enzymes and to reconstruct the general modes of transcription regulation for the corresponding genes. The genomic organization of the key genes is suggestive of a subdivision of the reaction network into five modules, where we observed distinct differences in the modular composition between the families Lactobacillaceae, Enterococcaceae, and Leuconostocaceae, on the one hand, and the family Streptococcaceae, on the other. These differences are mirrored by the way in which transcription regulation of the genes is structured in these families. In the Lactobacillaceae, Enterococcaceae, and Leuconostocaceae, the main shared mode of transcription regulation is methionine (Met) T-box-mediated regulation. In addition, the gene metK, encoding S-adenosylmethionine (SAM) synthetase, is controlled via the S(MK) box (SAM). The S(MK) box is also found upstream of metK in species of the family Streptococcaceae. However, the transcription control of the other modules is mediated via three different LysR-family regulators, MetR/MtaR (methionine), CmbR (O-acetyl[homo]serine), and HomR (O-acetylhomoserine). Redefinition of the associated DNA-binding motifs helped to identify/disentangle the related regulons, which appeared to perfectly match the proposed subdivision of the reaction network.

Complete resequencing and reannotation of the Lactobacillus plantarum WCFS1 genome.

There is growing interest in the beneficial effects of Lactobacillus plantarum on human health. The genome of L. plantarum WCFS1, first sequenced in 2001, was resequenced using Solexa technology. We identified 116 nucleotide corrections and improved function prediction for nearly 1,200 proteins, with a focus on metabolic functions and cell surface-associated proteins.

Genome sequence of the naturally plasmid-free Lactobacillus plantarum strain NC8 (CCUG 61730).

Lactobacillus plantarum is a highly versatile lactic acid bacterium found in various ecological niches, such as fermented vegetable, meat, and dairy products and the gastrointestinal tract. We sequenced the genome of L. plantarum NC8, a naturally plasmid-free strain, which has been used as a model strain in many laboratories worldwide.

PhenoLink--a web-tool for linking phenotype to ~omics data for bacteria: application to gene-trait matching for Lactobacillus plantarum strains.

BACKGROUND: Linking phenotypes to high-throughput molecular biology information generated by ~omics technologies allows revealing cellular mechanisms underlying an organism's phenotype. ~Omics datasets are often very large and noisy with many features (e.g., genes, metabolite abundances). Thus, associating phenotypes to ~omics data requires an approach that is robust to noise and can handle large and diverse data sets. RESULTS: We developed a web-tool PhenoLink (http://bamics2.cmbi.ru.nl/websoftware/phenolink/) that links phenotype to ~omics data sets using well-established as well new techniques. PhenoLink imputes missing values and preprocesses input data (i) to decrease inherent noise in the data and (ii) to counterbalance pitfalls of the Random Forest algorithm, on which feature (e.g., gene) selection is based. Preprocessed data is used in feature (e.g., gene) selection to identify relations to phenotypes. We applied PhenoLink to identify gene-phenotype relations based on the presence/absence of 2847 genes in 42 Lactobacillus plantarum strains and phenotypic measurements of these strains in several experimental conditions, including growth on sugars and nitrogen-dioxide production. Genes were ranked based on their importance (predictive value) to correctly predict the phenotype of a given strain. In addition to known gene to phenotype relations we also found novel relations. CONCLUSIONS: PhenoLink is an easily accessible web-tool to facilitate identifying relations from large and often noisy phenotype and ~omics datasets. Visualization of links to phenotypes offered in PhenoLink allows prioritizing links, finding relations between features, finding relations between phenotypes, and identifying outliers in phenotype data. PhenoLink can be used to uncover phenotype links to a multitude of ~omics data, e.g., gene presence/absence (determined by e.g.: CGH or next-generation sequencing), gene expression (determined by e.g.: microarrays or RNA-seq), or metabolite abundance (determined by e.g.: GC-MS).

Comparative genome analysis of central nitrogen metabolism and its control by GlnR in the class Bacilli.

BACKGROUND: The assimilation of nitrogen in bacteria is achieved through only a few metabolic conversions between alpha-ketoglutarate, glutamate and glutamine. The enzymes that catalyze these conversions are glutamine synthetase, glutaminase, glutamate dehydrogenase and glutamine alpha-ketoglutarate aminotransferase. In low-GC Gram-positive bacteria the transcriptional control over the levels of the related enzymes is mediated by four regulators: GlnR, TnrA, GltC and CodY. We have analyzed the genomes of all species belonging to the taxonomic families Bacillaceae, Listeriaceae, Staphylococcaceae, Lactobacillaceae, Leuconostocaceae and Streptococcaceae to determine the diversity in central nitrogen metabolism and reconstructed the regulation by GlnR. RESULTS: Although we observed a substantial difference in the extent of central nitrogen metabolism in the various species, the basic GlnR regulon was remarkably constant and appeared not affected by the presence or absence of the other three main regulators. We found a conserved regulatory association of GlnR with glutamine synthetase (glnRA operon), and the transport of ammonium (amtB-glnK) and glutamine/glutamate (i.e. via glnQHMP, glnPHQ, gltT, alsT). In addition less-conserved associations were found with, for instance, glutamate dehydrogenase in Streptococcaceae, purine catabolism and the reduction of nitrite in Bacillaceae, and aspartate/asparagine deamination in Lactobacillaceae. CONCLUSIONS: Our analyses imply GlnR-mediated regulation in constraining the import of ammonia/amino-containing compounds and the production of intracellular ammonia under conditions of high nitrogen availability. Such a role fits with the intrinsic need for tight control of ammonia levels to limit futile cycling.

Comparative analyses imply that the enigmatic Sigma factor 54 is a central controller of the bacterial exterior.

BACKGROUND: Sigma-54 is a central regulator in many pathogenic bacteria and has been linked to a multitude of cellular processes like nitrogen assimilation and important functional traits such as motility, virulence, and biofilm formation. Until now it has remained obscure whether these phenomena and the control by Sigma-54 share an underlying theme. RESULTS: We have uncovered the commonality by performing a range of comparative genome analyses. A) The presence of Sigma-54 and its associated activators was determined for all sequenced prokaryotes. We observed a phylum-dependent distribution that is suggestive of an evolutionary relationship between Sigma-54 and lipopolysaccharide and flagellar biosynthesis. B) All Sigma-54 activators were identified and annotated. The relation with phosphotransfer-mediated signaling (TCS and PTS) and the transport and assimilation of carboxylates and nitrogen containing metabolites was substantiated. C) The function annotations, that were represented within the genomic context of all genes encoding Sigma-54, its activators and its promoters, were analyzed for intra-phylum representation and inter-phylum conservation. Promoters were localized using a straightforward scoring strategy that was formulated to identify similar motifs. We found clear highly-represented and conserved genetic associations with genes that concern the transport and biosynthesis of the metabolic intermediates of exopolysaccharides, flagella, lipids, lipopolysaccharides, lipoproteins and peptidoglycan. CONCLUSION: Our analyses directly implicate Sigma-54 as a central player in the control over the processes that involve the physical interaction of an organism with its environment like in the colonization of a host (virulence) or the formation of biofilm.

Novel sigmaB regulation modules of Gram-positive bacteria involve the use of complex hybrid histidine kinases.

A common bacterial strategy to cope with stressful conditions is the activation of alternative sigma factors that control specific regulons enabling targeted responses. In the human pathogen Bacillus cereus, activation of the major stress-responsive sigma factor σ(B) is controlled by a signalling route that involves the multi-sensor hybrid histidine kinase RsbK. RsbK-type kinases are not restricted to the B. cereus group, but occur in a wide variety of other bacterial species, including members of the the low-GC Gram-positive genera Geobacillus and Paenibacillus as well as the high-GC actinobacteria. Genome context and protein sequence analyses of 118 RsbK homologues revealed extreme variability in N-terminal sensory as well as C-terminal regulatory domains and suggested that RsbK-type kinases are subject to complex fine-tuning systems, including sensitization and desensitization via methylation and demethylation within the helical domain preceding the H-box. The RsbK-mediated stress-responsive sigma factor activation mechanism that has evolved in B. cereus and the other species differs markedly from the extensively studied and highly conserved RsbRST-mediated σ(B) activation route found in Bacillus subtilis and other low-GC Gram-positive bacteria. Implications for future research on sigma factor control mechanisms are presented and current knowledge gaps are briefly discussed.

Preferential deletion events in the direct repeat locus of Mycobacterium tuberculosis.

The "Harlingen" IS6110 restriction fragment length polymorphism (RFLP) cluster has linked over 100 tuberculosis cases in The Netherlands since 1993. Four Mycobacterium tuberculosis isolates that were epidemiologically linked to this cluster had different spoligotype patterns, as well as slightly divergent IS6110 profiles, compared to the majority of the isolates. Sequencing of the direct repeat (DR) locus revealed sequence polymorphisms at the putative deletion sites. These deletion footprints provided evidence for independent deletions of the central region of the DR locus in three isolates, while the different genotype of the fourth isolate was explained by transmission. Our finding suggests that convergent deletions in the DR locus occur frequently. However, deletion footprints are not suitable to detect convergent deletions in the DR because they seem to be exceptional. Deletion footprints in the DR were not described previously, and we did not observe them in any public M. tuberculosis complex sequences. We conclude that preferential deletions in the DR loci of closely related strains are usually an unnoted event that interferes with clustering of closely related strains.

PanCGHweb: a web tool for genotype calling in pangenome CGH data.

UNLABELLED: A pangenome is the total of genes present in strains of the same species. Pangenome microarrays allow determining the genomic content of bacterial strains more accurately than conventional comparative genome hybridization microarrays. PanCGHweb is the first tool that effectively calls genotype based on pangenome microarray data. AVAILABILITY: PanCGHweb, the web tool is accessible from: http://bamics2.cmbi.ru.nl/websoftware/pancgh/.

Molecular description and industrial potential of Tn6098 conjugative transfer conferring alpha-galactoside metabolism in Lactococcus lactis.

A novel 51-kb conjugative transposon of Lactococcus lactis, designated Tn6098, encoding the capacity to utilize α-galactosides such as raffinose and stachyose, was identified and characterized. Alpha-galactosides are a dominant carbon source in many plant-derived foods. Most dairy lactococcus strains are unable to use α-galactosides as a growth substrate, yet many of these strains are known to have beneficial industrial traits. Conjugal transfer of Tn6098 was demonstrated from the plant-derived donor strain L. lactis KF147 to the recipient L. lactis NZ4501, a derivative of the dairy model strain L. lactis MG1363. The integration of Tn6098 into the genome of the recipient strain was confirmed by Illumina sequencing of the transconjugant L. lactis NIZO3921. The molecular structure of the integration site was confirmed by a PCR product spanning the insertion site. A 15-bp direct repeat sequence (TTATACCATAATTAC) is present on either side of Tn6098 in the chromosome of L. lactis KF147. One copy of this sequence is also present in the L. lactis MG1363 chromosome and represents the sole integration site. Phenotypic characterization of all strains showed that the transconjugant has not only acquired the ability to grow well in soy milk, a substrate rich in α-galactosides, but also has retained the flavor-forming capabilities of the recipient strain L. lactis MG1363. This study demonstrates how (induced) conjugation can be used to exploit the beneficial industrial traits of industrial dairy lactic acid bacteria in fermentation of plant-derived substrates.

Genomic diversity and versatility of Lactobacillus plantarum, a natural metabolic engineer.

In the past decade it has become clear that the lactic acid bacterium Lactobacillus plantarum occupies a diverse range of environmental niches and has an enormous diversity in phenotypic properties, metabolic capacity and industrial applications. In this review, we describe how genome sequencing, comparative genome hybridization and comparative genomics has provided insight into the underlying genomic diversity and versatility of L. plantarum. One of the main features appears to be genomic life-style islands consisting of numerous functional gene cassettes, in particular for carbohydrates utilization, which can be acquired, shuffled, substituted or deleted in response to niche requirements. In this sense, L. plantarum can be considered a "natural metabolic engineer".

Complete genome sequence of Lactococcus lactis subsp. lactis KF147, a plant-associated lactic acid bacterium.

Lactococcus lactis is a lactic acid bacterium used in the production of many fermented dairy products. We report the complete genome sequence of L. lactis subsp. lactis KF147, a nondairy strain isolated from mung bean sprouts. The circular chromosome of 2,598,144 bp, the largest among the sequenced lactococcal strains, encodes many properties related to adaptation to the plant environment.