Large plasmidome of dairy Lactococcus lactis subsp. lactis biovar diacetylactis FM03P encodes technological functions and appears highly unstable

Discovering genetic determinants for cell-to-cell adhesion in two prevalent conjugative lactococcal plasmids

Plasmids pNP40 and pUC11B encode two prevalent yet divergent conjugation systems, which have been characterized in detail recently. Here, we report the elucidation of the putative adhesins of the pNP40 and pUC11B conjugation systems, encoded by traAd and trsAd, respectively. Despite their significant sequence divergence, TraAd and TrsAd represent the most conserved component between the pNP40- and the pUC11B-encoded conjugation systems and share similar peptidoglycan-hydrolase domains. Protein structure prediction using AlphaFold2 highlighted the structural similarities between their predicted domains, as well as the potential homo-dimeric state of both proteins. Expression of the putative surface adhesins resulted in a cell clumping phenotype not only among cells expressing these surface adhesins but also between adhesin-expressing and non-producing cells. Furthermore, mutant derivatives of plasmids pNP40 or pUC11B carrying a mutation in traAd or trsAd, respectively, were shown to act as efficient donors provided the corresponding recipient expresses either traAd or trsAd, thus demonstrating in trans reciprocal complementarity of these proteins in conjugation systems.

The stressostat: A novel approach in adaptive laboratory evolution to improve end-product resistance

End-product inhibition in pH-controlled batch cultures, is the major limiting factor for bacterial biomass formation in the starter culture industry as well as in many other biotechnological processes. Adaptive laboratory evolution (ALE) has emerged over the past decades as a powerful tool for phenotype optimization, but none of the existing ALE methods could select for improved end-product resistance. Therefore, we developed the stressostat (STress Resistance Evolution in Substrate Surplus) as a novel continuous ALE method. Stressostat cultivation applies end-product concentrations as constant evolutionary pressure on microorganisms in the presence of substrate surplus. In this study, we improved the lactate resistance of Lactococcus lactis FM03P in 35 days of stressostat cultivations. The lactate concentrations increased over time from 530 to 675 mM, indicating the successful selection for variants with improved lactate resistance. Thirty-four variants were isolated and grouped into four clusters based on their growth rates at high lactate concentrations. In the high-throughput screening without pH control, most isolated variants could grow at high lactate concentrations (870-928 mM), while the wild type was completely inhibited. The variants grew slower than wild type at low lactate media indicating possible evolutionary trade-off. However, in pH-controlled batch cultivations, most variants produced more biomass than the wild type. In conclusion, stressostat cultivation is a valuable method to obtain L. lactis variants with improved end-product resistance and further characterization is needed to elucidate underlying resistance mechanisms and potential industrial applications.

Slc11 Synapomorphy: A Conserved 3D Framework Articulating Carrier Conformation Switch

Transmembrane carriers of the Slc11 family catalyze proton (H+)-dependent uptake of divalent metal ions (Me2+) such as manganese and iron-vital elements coveted during infection. The Slc11 mechanism of high-affinity Me2+ cell import is selective and conserved between prokaryotic (MntH) and eukaryotic (Nramp) homologs, though processes coupling the use of the proton motive force to Me2+ uptake evolved repeatedly. Adding bacterial piracy of Nramp genes spread in distinct environmental niches suggests selective gain of function that may benefit opportunistic pathogens. To better understand Slc11 evolution, Alphafold (AF2)/Colabfold (CF) 3D predictions for bacterial sequences from sister clades of eukaryotic descent (MCb and MCg) were compared using both native and mutant templates. AF2/CF model an array of native MCb intermediates spanning the transition from outwardly open (OO) to inwardly open (IO) carriers. In silico mutagenesis targeting (i) a set of (evolutionarily coupled) sites that may define Slc11 function (putative synapomorphy) and (ii) residues from networked communities evolving during MCb transition indicates that Slc11 synapomorphy primarily instructs a Me2+-selective conformation switch which unlocks carrier inner gate and contributes to Me2+ binding site occlusion and outer gate locking. Inner gate opening apparently proceeds from interaction between transmembrane helix (h) h5, h8 and h1a. MCg1 xenologs revealed marked differences in carrier shape and plasticity, owing partly to an altered intramolecular H+ network. Yet, targeting Slc11 synapomorphy also converted MCg1 IO models to an OO state, apparently mobilizing the same residues to control gates. But MCg1 response to mutagenesis differed, with extensive divergence within this clade correlating with MCb-like modeling properties. Notably, MCg1 divergent epistasis marks the emergence of the genus Bordetella-Achromobacter. Slc11 synapomorphy localizes to the 3D areas that deviate least among MCb and MCg1 models (either IO or OO) implying that it constitutes a 3D network of residues articulating a Me2+-selective carrier conformation switch which is maintained in fast-evolving clades at the cost of divergent epistatic interactions impacting carrier shape and dynamics.

Control of Bacterial Phenotype and Chromosomal Gene Expression by Single Plasmids of Lactococcus lactis IL594

Plasmid-free Lactococcus lactis IL1403 is one of the best-characterized representatives of lactic acid bacteria (LAB), intensively used in broad microbiology worldwide. Its parent strain, L. lactis IL594, contains seven plasmids (pIL1-pIL7) with resolved DNA sequences and an indicated role for overall plasmid load in enhancing host-adaptive potential. To determine how individual plasmids manipulate the expression of phenotypes and chromosomal genes, we conducted global comparative phenotypic analyses combined with transcriptomic studies in plasmid-free L. lactis IL1403, multiplasmid L. lactis IL594, and its single-plasmid derivatives. The presence of pIL2, pIL4, and pIL5 led to the most pronounced phenotypic differences in the metabolism of several carbon sources, including some β-glycosides and organic acids. The pIL5 plasmid also contributed to increased tolerance to some antimicrobial compounds and heavy metal ions, especially those in the toxic cation group. Comparative transcriptomics showed significant variation in the expression levels of up to 189 chromosomal genes due to the presence of single plasmids and 435 unique chromosomal genes that were resultant of the activity of all plasmids, which may suggest that the observed phenotypic changes are not only the result of a direct action of their own genes but also originate from indirect actions through crosstalk between plasmids and the chromosome. The data obtained here indicate that plasmid maintenance leads to the development of important mechanisms of global gene regulation that provide changes in the central metabolic pathways and adaptive properties of L. lactis and suggest the possibility of a similar phenomenon among other groups of bacteria.

Incorporating plasmid biology and metagenomics into a holistic model of the human gut microbiome

The human gut microbiome is often described as the collection of bacteria, archaea, fungi, protists, and viruses associated with an individual, with no acknowledgement of the plasmid constituents. However, like viruses, plasmids are autonomous intracellular replicating entities that can influence the genotype and phenotype of their host and mediate trans-kingdom interactions. Plasmids are frequently noted as vehicles for horizontal gene transfer and for spreading antibiotic resistance, yet their multifaceted contribution to mutualistic and antagonistic interactions within the human microbiome and impact on human health is overlooked. In this review, we highlight the importance of plasmids and their biological properties as overlooked components of microbiomes. Subsequent human microbiome studies should include dedicated analyses of plasmids, particularly as a holistic understanding of human-microbial interactions is required before effective and safe interventions can be implemented to improve human well-being.

The Presence of Plasmids in Lactococcus lactis IL594 Determines Changes in the Host Phenotype and Expression of Chromosomal Genes

The L. lactis IL594 strain contains seven plasmids (pIL1 to pIL7) and is the parental strain of the plasmid-free L. lactis IL1403, one of the most studied lactic acid bacteria (LAB) strain. The genetic sequences of pIL1 to pIL7 plasmids have been recently described, however the knowledge of global changes in host phenotype and transcriptome remains poor. In the present study, global phenotypic analyses were combined with transcriptomic studies to evaluate a potential influence of plasmidic genes on overall gene expression in industrially important L. lactis strains. High-throughput screening of phenotypes differences revealed pronounced phenotypic differences in favor of IL594 during the metabolism of some C-sources, including lactose and β-glucosides. A plasmids-bearing strain presented increased resistance to unfavorable growth conditions, including the presence of heavy metal ions and antimicrobial compounds. Global comparative transcriptomic study of L. lactis strains revealed variation in the expression of over 370 of chromosomal genes caused by plasmids presence. The general trend presented upregulated energy metabolism and biosynthetic genes, differentially expressed regulators, prophages and cell resistance proteins. Our findings suggest that plasmids maintenance leads to significant perturbation in global gene regulation that provides change in central metabolic pathways and adaptive properties of the IL594 cells.

Phenotypic and Safety Assessment of the Cheese Strain Lactiplantibacillus plantarum LL441, and Sequence Analysis of its Complete Genome and Plasmidome

This work describes the phenotypic typing and complete genome analysis of LL441, a dairy Lactiplantibacillus plantarum strain. LL441 utilized a large range of carbohydrates and showed strong activity of some carbohydrate-degrading enzymes. The strain grew slowly in milk and produced acids and ketones along with other volatile compounds. The genome of LL441 included eight circular molecules, the bacterial chromosome, and seven plasmids (pLL441-1 through pLL441-7), ranging in size from 8.7 to 53.3 kbp. Genome analysis revealed vast arrays of genes involved in carbohydrate utilization and flavor formation in milk, as well as genes providing acid and bile resistance. No genes coding for virulence traits or pathogenicity factors were detected. Chromosome and plasmids were packed with insertion sequence (IS) elements. Plasmids were also abundant in genes encoding heavy metal resistance traits and plasmid maintenance functions. Technologically relevant phenotypes linked to plasmids, such as the production of plantaricin C (pLL441-1), lactose utilization (pLL441-2), and bacteriophage resistance (pLL441-4), were also identified. The absence of acquired antibiotic resistance and of phenotypes and genes of concern suggests L. plantarum LL441 be safe. The strain might therefore have a use as a starter or starter component in dairy and other food fermentations or as a probiotic.

Food grade microbial synthesis of the butter aroma compound butanedione using engineered and non-engineered Lactococcus lactis

The design-build-test-learn (DBTL) cycle has been implemented in metabolic engineering processes for optimizing the production of valuable compounds, including food ingredients. However, the use of recombinant microorganisms for producing food ingredients is associated with different challenges, e.g., in the EU, a content of more than 0.9% of such ingredients requires to be labeled. Therefore, we propose to expand the DBTL cycle and use the "learn" module to guide the development of non-engineered strains for clean label production. Here, we demonstrate how this approach can be used to generate engineered and natural cell factories able to produce the valuable food flavor compound - butanedione (diacetyl). Through comprehensive rerouting of the metabolism of Lactococcus lactis MG1363 and re-installment of the capacity to metabolize lactose and dairy protein, we managed to achieve a high titer of diacetyl (6.7 g/L) in pure dairy waste. Based on learnings from the engineering efforts, we successfully achieved the production of diacetyl without using recombinant DNA technology. We accomplish the latter by process optimization and by relying on high-throughput screening using a microfluidic system. Our results demonstrate the great potential that lies in combining metabolic engineering and natural approaches for achieving efficient production of food ingredients.

Directed Recovery and Molecular Characterization of Antibiotic Resistance Plasmids from Cheese Bacteria

Resistance to antimicrobials is a growing problem of worldwide concern. Plasmids are thought to be major drivers of antibiotic resistance spread. The present work reports a simple way to recover replicative plasmids conferring antibiotic resistance from the bacteria in cheese. Purified plasmid DNA from colonies grown in the presence of tetracycline and erythromycin was introduced into plasmid-free strains of Lactococcus lactis, Lactiplantibacillus plantarum and Lacticaseibacillus casei. Following antibiotic selection, the plasmids from resistant transformants were isolated, analyzed by restriction enzyme digestion, and sequenced. Seven patterns were obtained for the tetracycline-resistant colonies, five from L. lactis, and one each from the lactobacilli strains, as well as a single digestion profile for the erythromycin-resistant transformants obtained in L. lactis. Sequence analysis respectively identified tet(S) and ermB in the tetracycline- and erythromycin-resistance plasmids from L. lactis. No dedicated resistance genes were detected in plasmids conferring tetracycline resistance to L. casei and L. plantarum. The present results highlight the usefulness of the proposed methodology for isolating functional plasmids that confer antibiotic resistance to LAB species, widen our knowledge of antibiotic resistance in the bacteria that inhabit cheese, and emphasize the leading role of plasmids in the spread of resistance genes via the food chain.

Lifestyle, metabolism and environmental adaptation in Lactococcus lactis

Lactococcus lactis serves as a paradigm organism for the lactic acid bacteria (LAB). Extensive research into the molecular biology, metabolism and physiology of several model strains of this species has been fundamental for our understanding of the LAB. Genomic studies have provided new insights into the species L. lactis, including the resolution of the genetic basis of its subspecies division, as well as the control mechanisms involved in the fine-tuning of growth rate and energy metabolism. In addition, it has enabled novel approaches to study lactococcal lifestyle adaptations to the dairy application environment, including its adjustment to near-zero growth rates that are particularly relevant in the context of cheese ripening. This review highlights various insights in these areas and exemplifies the strength of combining experimental evolution with functional genomics and bacterial physiology research to expand our fundamental understanding of the L. lactis lifestyle under different environmental conditions.

The large plasmidome of Lactococcus lactis subsp. lactis bv. diacetylactis S50 confers its biotechnological properties

Plasmids are autonomous episomally replicating genetic elements, which carry backbone genes important for the replication and maintenance within their host, and accessory genes that might confer an advantage to their host under specific selective pressure in its ecological niche. The genome of dairy isolate L. lactis subsp. lactis bv. diacetylactis S50 was sequenced using the PacBio SMRT Cell Seq-RSII platform and revealed to possess one of the largest plasmidomes among L. lactis strains studied so far, harboring six plasmids: pS6 (5553 bp), pS7a (7308 bp), pS7b (7266 bp), pS19 (19,027 bp), pS74 (74,256 bp) and pS127 (127,002 bp) in total representing 8.9% of genome size (240,412 bp). Based on predicted plasmid replication proteins and origins it appears that all six plasmids replicate via the theta-type mechanism. The two the largest plasmids (pS74 and pS127), carry a number of genes known to be important for growth and survival in the dairy environment. These genes encode technological functions such as bacteriocin production, protein degradation, magnesium and cobalt/nickel transporters, selenium binding, exopolysaccharides (EPS) production, bacteriophage and stress resistance. Beside genes for replication, the small plasmids (pS6, pS7a, pS7a, and pS19) also carry genes important for mobilization and host survival such as type I restriction-modification (R-M) system, metal transporters, enzymes and transcriptional regulators. All plasmids in S50 strain are mobilizable, containing an oriT sequences, while pS127 is self-conjugative and allows for mobilization of the other plasmids. Small plasmids are prone to structural and segregational instability, while pS127 appeared to be segregationally stable thanks to the possession of two partition systems. The main characteristic of plasmid pS74 is EPS production, while plasmid pS127 is characterized by proteinase and multiple bacteriocins, tra locus, phage abortive systems and metal transporters. In addition to LcnA and LcnB, plasmid pS127 encodes several bacteriocin-pheromone molecules and a new bacteriocin named LcnS50, with narrow spectrum of action limited to lactococci, that has been successfully cloned and heterologously expressed.

Comparative Genomics of Streptococcus thermophilus Support Important Traits Concerning the Evolution, Biology and Technological Properties of the Species

Streptococcus thermophilus is a major starter for the dairy industry with great economic importance. In this study we analyzed 23 fully sequenced genomes of S. thermophilus to highlight novel aspects of the evolution, biology and technological properties of this species. Pan/core genome analysis revealed that the species has an important number of conserved genes and that the pan genome is probably going to be closed soon. According to whole genome phylogeny and average nucleotide identity (ANI) analysis, most S. thermophilus strains were grouped in two major clusters (i.e., clusters A and B). More specifically, cluster A includes strains with chromosomes above 1.83 Mbp, while cluster B includes chromosomes below this threshold. This observation suggests that strains belonging to the two clusters may be differentiated by gene gain or gene loss events. Furthermore, certain strains of cluster A could be further subdivided in subgroups, i.e., subgroup I (ASCC 1275, DGCC 7710, KLDS SM, MN-BM-A02, and ND07), II (MN-BM-A01 and MN-ZLW-002), III (LMD-9 and SMQ-301), and IV (APC151 and ND03). In cluster B certain strains formed one distinct subgroup, i.e., subgroup I (CNRZ1066, CS8, EPS, and S9). Clusters and subgroups observed for S. thermophilus indicate the existence of lineages within the species, an observation which was further supported to a variable degree by the distribution and/or the architecture of several genomic traits. These would include exopolysaccharide (EPS) gene clusters, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs)-CRISPR associated (Cas) systems, as well as restriction-modification (R-M) systems and genomic islands (GIs). Of note, the histidine biosynthetic cluster was found present in all cluster A strains (plus strain NCTC12958^T) but was absent from all strains in cluster B. Other loci related to lactose/galactose catabolism and urea metabolism, aminopeptidases, the majority of amino acid and peptide transporters, as well as amino acid biosynthetic pathways were found to be conserved in all strains suggesting their central role for the species. Our study highlights the necessity of sequencing and analyzing more S. thermophilus complete genomes to further elucidate important aspects of strain diversity within this starter culture that may be related to its application in the dairy industry.

DAIRYdb: a manually curated reference database for improved taxonomy annotation of 16S rRNA gene sequences from dairy products

Background

Reads assignment to taxonomic units is a key step in microbiome analysis pipelines. To date, accurate taxonomy annotation of 16S reads, particularly at species rank, is still challenging due to the short size of read sequences and differently curated classification databases. The close phylogenetic relationship between species encountered in dairy products, however, makes it crucial to annotate species accurately to achieve sufficient phylogenetic resolution for further downstream ecological studies or for food diagnostics. Curated databases dedicated to the environment of interest are expected to improve the accuracy and resolution of taxonomy annotation.

Results

We provide a manually curated database composed of 10’290 full-length 16S rRNA gene sequences from prokaryotes tailored for dairy products analysis (https://github.com/marcomeola/DAIRYdb).

The performance of the DAIRYdb was compared with the universal databases Silva, LTP, RDP and Greengenes. The DAIRYdb significantly outperformed all other databases independently of the classification algorithm by enabling higher accurate taxonomy annotation down to the species rank. The DAIRYdb accurately annotates over 90% of the sequences of either single or paired hypervariable regions automatically.

The manually curated DAIRYdb strongly improves taxonomic annotation accuracy for microbiome studies in dairy environments. The DAIRYdb is a practical solution that enables automatization of this key step, thus facilitating the routine application of NGS microbiome analyses for microbial ecology studies and diagnostics in dairy products.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5914-8) contains supplementary material, which is available to authorized users.

Long-read based de novo assembly of low-complexity metagenome samples results in finished genomes and reveals insights into strain diversity and an active phage system

BACKGROUND

Complete and contiguous genome assemblies greatly improve the quality of subsequent systems-wide functional profiling studies and the ability to gain novel biological insights. While a de novo genome assembly of an isolated bacterial strain is in most cases straightforward, more informative data about co-existing bacteria as well as synergistic and antagonistic effects can be obtained from a direct analysis of microbial communities. However, the complexity of metagenomic samples represents a major challenge. While third generation sequencing technologies have been suggested to enable finished metagenome-assembled genomes, to our knowledge, the complete genome assembly of all dominant strains in a microbiome sample has not been demonstrated. Natural whey starter cultures (NWCs) are used in cheese production and represent low-complexity microbiomes. Previous studies of Swiss Gruyère and selected Italian hard cheeses, mostly based on amplicon metagenomics, concurred that three species generally pre-dominate: Streptococcus thermophilus, Lactobacillus helveticus and Lactobacillus delbrueckii.

RESULTS

Two NWCs from Swiss Gruyère producers were subjected to whole metagenome shotgun sequencing using the Pacific Biosciences Sequel and Illumina MiSeq platforms. In addition, longer Oxford Nanopore Technologies MinION reads had to be generated for one to resolve repeat regions. Thereby, we achieved the complete assembly of all dominant bacterial genomes from these low-complexity NWCs, which was corroborated by a 16S rRNA amplicon survey. Moreover, two distinct L. helveticus strains were successfully co-assembled from the same sample. Besides bacterial chromosomes, we could also assemble several bacterial plasmids and phages and a corresponding prophage. Biologically relevant insights were uncovered by linking the plasmids and phages to their respective host genomes using DNA methylation motifs on the plasmids and by matching prokaryotic CRISPR spacers with the corresponding protospacers on the phages. These results could only be achieved by employing long-read sequencing data able to span intragenomic as well as intergenomic repeats.

CONCLUSIONS

Here, we demonstrate the feasibility of complete de novo genome assembly of all dominant strains from low-complexity NWCs based on whole metagenomics shotgun sequencing data. This allowed to gain novel biological insights and is a fundamental basis for subsequent systems-wide omics analyses, functional profiling and phenotype to genotype analysis of specific microbial communities.

Omics-based comparative analysis of putative mobile genetic elements in Lactococcus lactis

Lactococcus lactis is globally used in food fermentation. Genomics is useful to investigate speciation and differential occurrence of (un)desired gene functions, often related to mobile DNA. This study investigates L. lactis for putative chromosomal mobile genetic elements through comparative genomics, and analyses how they contribute to chromosomal variation at strain level. Our work identified 95 loci that may range over 10% of the chromosome size when including prophages, and the loci display a marked differential occurrence in the analysed strains. Analysis of differential transcriptomics data revealed how mobile genetic elements may impact the host physiology in response to conditional changes. This insight in the genetic variation of mobile genetic elements in L. lactis holds potential to further identify important functions related to food and biotechnology applications within this important species.

The Lactococcus lactis Pan-Plasmidome

Plasmids are autonomous, self-replicating, extrachromosomal genetic elements that are typically not essential for growth of their host. They may encode metabolic capabilities, which promote the maintenance of these genetic elements, and may allow adaption to specific ecological niches and consequently enhance survival. Genome sequencing of 16 Lactococcus lactis strains revealed the presence of 83 plasmids, including two megaplasmids. The limitations of Pacific Biosciences SMRT sequencing in detecting the total plasmid complement of lactococcal strains is examined, while a combined Illumina/SMRT sequencing approach is proposed to combat these issues. Comparative genome analysis of these plasmid sequences combined with other publicly available plasmid sequence data allowed the definition of the lactococcal plasmidome, and facilitated an investigation into (bio) technologically important plasmid-encoded traits such as conjugation, bacteriocin production, exopolysaccharide (EPS) production, and (bacterio) phage resistance.

Large Circular Plasmids from Groundwater Plasmidomes Span Multiple Incompatibility Groups and Are Enriched in Multimetal Resistance Genes

Naturally occurring plasmids constitute a major category of mobile genetic elements responsible for harboring and transferring genes important in survival and fitness. A targeted evaluation of plasmidomes can reveal unique adaptations required by microbial communities. We developed a model system to optimize plasmid DNA isolation procedures targeted to groundwater samples which are typically characterized by low cell density (and likely variations in the plasmid size and copy numbers). The optimized method resulted in successful identification of several hundred circular plasmids, including some large plasmids (11 plasmids more than 50 kb in size, with the largest being 1.7 Mb in size). Several interesting observations were made from the analysis of plasmid DNA isolated in this study. The plasmid pool (plasmidome) was more conserved than the corresponding microbiome distribution (16S rRNA based). The circular plasmids were diverse as represented by the presence of seven plasmid incompatibility groups. The genes carried on these groundwater plasmids were highly enriched in metal resistance. Results from this study confirmed that traits such as metal, antibiotic, and phage resistance along with toxin-antitoxin systems are encoded on abundant circular plasmids, all of which could confer novel and advantageous traits to their hosts. This study confirms the ecological role of the plasmidome in maintaining the latent capacity of a microbiome, enabling rapid adaptation to environmental stresses.IMPORTANCE Plasmidomes have been typically studied in environments abundant in bacteria, and this is the first study to explore plasmids from an environment characterized by low cell density. We specifically target groundwater, a significant source of water for human/agriculture use. We used samples from a well-studied site and identified hundreds of circular plasmids, including one of the largest sizes reported in plasmidome studies. The striking similarity of the plasmid-borne ORFs in terms of taxonomical and functional classifications across several samples suggests a conserved plasmid pool, in contrast to the observed variability in the 16S rRNA-based microbiome distribution. Additionally, the stress response to environmental factors has stronger conservation via plasmid-borne genes as marked by abundance of metal resistance genes. Last, identification of novel and diverse plasmids enriches the existing plasmid database(s) and serves as a paradigm to increase the repertoire of biological parts that are available for modifying novel environmental strains.

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.