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.

Comparative Genome Analysis of Lactococcus lactis Indicates Niche Adaptation and Resolves Genotype/Phenotype Disparity.

Lactococcus lactis is one of the most important micro-organisms in the dairy industry for the fermentation of cheese and buttermilk. Besides the conversion of lactose to lactate it is responsible for product properties such as flavor and texture, which are determined by volatile metabolites, proteolytic activity and exopolysaccharide production. While the species Lactococcus lactis consists of the two subspecies lactis and cremoris their taxonomic position is confused by a group of strains that, despite of a cremoris genotype, display a lactis phenotype. Here we compared and analyzed the (draft) genomes of 43 L. lactis strains, of which 19 are of dairy and 24 are of non-dairy origin. Machine-learning algorithms facilitated the identification of orthologous groups of protein sequences (OGs) that are predictors for either the taxonomic position or the source of isolation. This allowed the unambiguous categorization of the genotype/phenotype disparity of ssp. lactis and ssp. cremoris strains. A detailed analysis of phenotypic properties including plasmid-encoded genes indicates evolutionary changes during niche adaptations. The results are consistent with the hypothesis that dairy isolates evolved from plant isolates. The analysis further suggests that genomes of cremoris phenotype strains are so eroded that they are restricted to a dairy environment. Overall the genome comparison of a diverse set of strains allowed the identification of niche and subspecies specific genes. This explains evolutionary relationships and will aid the identification and selection of industrial starter cultures.

Putative antibiotic resistance genes present in extant Bacillus licheniformis and Bacillus paralicheniformis strains are probably intrinsic and part of the ancient resistome.

Whole-genome sequencing and phenotypic testing of 104 strains of Bacillus licheniformis and Bacillus paralicheniformis from a variety of sources and time periods was used to characterize the genetic background and evolution of (putative) antimicrobial resistance mechanisms. Core proteins were identified in draft genomes and a phylogenetic analysis based on single amino acid polymorphisms allowed the species to be separated into two phylogenetically distinct clades with one outlier. Putative antimicrobial resistance genes were identified and mapped. A chromosomal ermD gene was found at the same location in all B. paralichenformis and in 27% of B. licheniformis genomes. Erythromycin resistance correlated very well with the presence of ermD. The putative streptomycin resistance genes, aph and aadK, were found in the chromosome of all strains as adjacent loci. Variations in amino acid sequence did not correlate with streptomycin susceptibility although the species were less susceptible than other Bacillus species. A putative chloramphenicol resistance gene (cat), encoding a novel chloramphenicol acetyltransferase protein was also found in the chromosome of all strains. Strains encoding a truncated CAT protein were sensitive to chloramphenicol. For all four resistance genes, the diversity and genetic context followed the overall phylogenetic relationship. No potentially mobile genetic elements were detected in their vicinity. Moreover, the genes were only distantly related to previously-described cat, aph, aad and erm genes present on mobile genetic elements or in other species. Thus, these genes are suggested to be intrinsic to B. licheniformis and B. paralicheniformis and part of their ancient resistomes. Since there is no evidence supporting horizontal transmission, these genes are not expected to add to the pool of antibiotic resistance elements considered to pose a risk to human or animal health. Whole-genome based phylogenetic and sequence analysis, combined with phenotypic testing, is proposed to be suitable for determining intrinsic resistance and evolutionary relationships.

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.

Plasmid Complement of Lactococcus lactis NCDO712 Reveals a Novel Pilus Gene Cluster.

Lactococcus lactis MG1363 is an important gram-positive model organism. It is a plasmid-free and phage-cured derivative of strain NCDO712. Plasmid-cured strains facilitate studies on molecular biological aspects, but many properties which make L. lactis an important organism in the dairy industry are plasmid encoded. We sequenced the total DNA of strain NCDO712 and, contrary to earlier reports, revealed that the strain carries 6 rather than 5 plasmids. A new 50-kb plasmid, designated pNZ712, encodes functional nisin immunity (nisCIP) and copper resistance (lcoRSABC). The copper resistance could be used as a marker for the conjugation of pNZ712 to L. lactis MG1614. A genome comparison with the plasmid cured daughter strain MG1363 showed that the number of single nucleotide polymorphisms that accumulated in the laboratory since the strains diverted more than 30 years ago is limited to 11 of which only 5 lead to amino acid changes. The 16-kb plasmid pSH74 was found to contain a novel 8-kb pilus gene cluster spaCB-spaA-srtC1-srtC2, which is predicted to encode a pilin tip protein SpaC, a pilus basal subunit SpaB, and a pilus backbone protein SpaA. The sortases SrtC1/SrtC2 are most likely involved in pilus polymerization while the chromosomally encoded SrtA could act to anchor the pilus to peptidoglycan in the cell wall. Overexpression of the pilus gene cluster from a multi-copy plasmid in L. lactis MG1363 resulted in cell chaining, aggregation, rapid sedimentation and increased conjugation efficiency of the cells. Electron microscopy showed that the over-expression of the pilus gene cluster leads to appendices on the cell surfaces. A deletion of the gene encoding the putative basal protein spaB, by truncating spaCB, led to more pilus-like structures on the cell surface, but cell aggregation and cell chaining were no longer observed. This is consistent with the prediction that spaB is involved in the anchoring of the pili to the cell.

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.

Identification of food-grade subtilisins as gluten-degrading enzymes to treat celiac disease.

Gluten are proline- and glutamine-rich proteins present in wheat, barley, and rye and contain the immunogenic sequences that drive celiac disease (CD). Rothia mucilaginosa, an oral microbial colonizer, can cleave these gluten epitopes. The aim was to isolate and identify the enzymes and evaluate their potential as novel enzyme therapeutics for CD. The membrane-associated R. mucilaginosa proteins were extracted and separated by DEAE chromatography. Enzyme activities were monitored with paranitroanilide-derivatized and fluorescence resonance energy transfer (FRET) peptide substrates, and by gliadin zymography. Epitope elimination was determined in R5 and G12 ELISAs. The gliadin-degrading Rothia enzymes were identified by LC-ESI-MS/MS as hypothetical proteins ROTMU0001_0241 (C6R5V9_9MICC), ROTMU0001_0243 (C6R5W1_9MICC), and ROTMU0001_240 (C6R5V8_9MICC). A search with the Basic Local Alignment Search Tool revealed that these are subtilisin-like serine proteases belonging to the peptidase S8 family. Alignment of the major Rothia subtilisins indicated that all contain the catalytic triad with Asp (D), His (H), and Ser (S) in the D-H-S order. They cleaved succinyl-Ala-Ala-Pro-Phe-paranitroanilide, a substrate for subtilisin with Pro in the P2 position, as in Tyr-Pro-Gln and Leu-Pro-Tyr in gluten, which are also cleaved. Consistently, FRET substrates of gliadin immunogenic epitopes comprising Xaa-Pro-Xaa motives were rapidly hydrolyzed. The Rothia subtilisins and two subtilisins from Bacillus licheniformis, subtilisin A and the food-grade Nattokinase, efficiently degraded the immunogenic gliadin-derived 33-mer peptide and the immunodominant epitopes recognized by the R5 and G12 antibodies. This study identified Rothia and food-grade Bacillus subtilisins as promising new candidates for enzyme therapeutics in CD.

Comparative Genomics of Iron-Transporting Systems in Bacillus cereus Strains and Impact of Iron Sources on Growth and Biofilm Formation.

Iron is an important element for bacterial viability, however it is not readily available in most environments. We studied the ability of 20 undomesticated food isolates of Bacillus cereus and two reference strains for capacity to use different (complex) iron sources for growth and biofilm formation. Studies were performed in media containing the iron scavenger 2,2-Bipyridine. Transcriptome analysis using B. cereus ATCC 10987 indeed showed upregulation of predicted iron transporters in the presence of 2,2-Bipyridine, confirming that iron was depleted upon its addition. Next, the impact of iron sources on growth performance of the 22 strains was assessed and correlations between growth stimulation and presence of putative iron transporter systems in the genome sequences were analyzed. All 22 strains effectively used Fe citrate and FeCl3 for growth, and possessed genes for biosynthesis of the siderophore bacillibactin, whereas seven strains lacked genes for synthesis of petrobactin. Hemoglobin could be used by all strains with the exception of one strain that lacked functional petrobactin and IlsA systems. Hemin could be used by the majority of the tested strains (19 of 22). Notably, transferrin, ferritin, and lactoferrin were not commonly used by B. cereus for growth, as these iron sources could be used by 6, 3, and 2 strains, respectively. Furthermore, biofilm formation was found to be affected by the type of iron source used, including stimulation of biofilms at liquid-air interphase (FeCl3 and Fe citrate) and formation of submerged type biofilms (hemin and lactoferrin). Our results show strain variability in the genome-encoded repertoire of iron-transporting systems and differences in efficacy to use complex iron sources for growth and biofilm formation. These features may affect B. cereus survival and persistence in specific niches.

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.

Nearly Complete Genome Sequence of Lactobacillus plantarum Strain NIZO2877.

Lactobacillus plantarum is a versatile bacterial species that is isolated mostly from foods. Here, we present the first genome sequence of L. plantarum strain NIZO2877 isolated from a hot dog in Vietnam. Its two contigs represent a nearly complete genome sequence.

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.

Resequencing of the Lactobacillus plantarum Strain WJL Genome.

Lactobacillus plantarum strain WJL is a symbiont isolated from the Drosophila melanogaster gut. The genome of L. plantarum WJL, first sequenced in 2013, was resequenced and rescaffolded in this study. A combination of Sanger and Illumina sequencing allowed us to reduce the number of contigs from 102 to 13. This work contributes to a better understanding of the genome and function of this organism.

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.