Lactococcus lactis Diversity in Undefined Mixed Dairy Starter Cultures as Revealed by Comparative Genome Analyses and Targeted Amplicon Sequencing of epsD

Database selection for shotgun metaproteomic of low-diversity dairy microbiomes

The metaproteomics field has recently gained more and more interest as a valuable tool for studying both the taxonomy and function of microbiomes, including those used in food fermentations. One crucial step in the metaproteomics pipeline is selecting a database to obtain high-quality taxonomical and functional information from microbial communities. One of the best strategies described for building protein databases is using sample-specific or study-specific protein databases obtained from metagenomic sequencing. While this is true for high-diversity microbiomes (such as gut and soil), there is still a lack of validation for different database construction strategies in low-diversity microbiomes, such as those found in fermented dairy products where starter cultures containing few species are used. In this study, we assessed the performance of various database construction strategies applied to metaproteomics on two low-diversity microbiomes obtained from cheese production using commercial starter cultures and analyzed by LC-MS/MS. Substantial differences were detected between the strategies, and the best performance in terms of the number of peptides and proteins identified from the spectra was achieved by metagenomic-derived databases. However, extensive databases constructed from a high number of available online genomes obtained a similar taxonomical and functional annotation of the metaproteome compared to the metagenomic-derived databases. Our results indicate that, in the case of low-diversity dairy microbiomes, the use of publically available genomes to construct protein databases can be considered as an alternative to metagenome-derived databases.

A multifaceted investigation of lactococcal strain diversity in undefined mesophilic starter cultures

The dairy fermentation industry relies on the activity of lactic acid bacteria in robust starter cultures to accomplish milk acidification. Maintenance of the composition of these starter cultures, whether defined or undefined, is essential to ensure consistent and high-quality fermentation end products. To date, limited information exists regarding the microbial composition of undefined starter culture systems. Here, we describe a culture-based analysis combined with a metagenomics approach to evaluate the composition of two undefined mesophilic starter cultures. In addition, we describe a qPCR-based genotype detection assay, which is capable of discerning nine distinct lactococcal genotypes to characterize these undefined starter cultures, and which can be applied to monitor compositional changes in an undefined starter culture during a fermentation.

IMPORTANCE

This study reports on the development of a combined culture-based analysis and metagenomics approach to evaluate the composition of two undefined mesophilic starter cultures. In addition, a novel qPCR-based genotype detection assay, capable of discerning nine distinct lactococcal genotypes (based on lactococcal cell wall polysaccharide biosynthesis gene clusters), was used to monitor compositional changes in an undefined starter culture following phage attack. These analytical approaches facilitate a multifaceted assessment of starter culture compositional stability during milk fermentation, which has become an important QC aspect due to the increasing demand for consistent and high-quality dairy products.

Towards the diversification of lactococcal starter and non-starter species in mesophilic dairy culture systems

Lactococcus is one of the earliest identified fermentative bacterial genera and among its member species, the dairy-associated Lactococcus lactis and Lactococcus cremoris are undoubtedly the best studied. These two species are believed to have evolved from plant-associated lactococci and through genome decay and acquisition of plasmids, have adapted to the dairy niche. The past decade has witnessed a surge of activity in novel lactococcal species identification from insect, plant and animal sources. Currently, 22 Lactococcus species are described and in this review, we summarise the genome characteristics of and phylogenetic relationships among these species. Furthermore, we explore the role of mobile elements including plasmids and bacteriophages in the diversification of lactococcal species. The pace of identification of novel lactococcal species suggests that the number of lactococcal species is likely to continue to grow. With additional sequence data for the emerging species, it will be possible to perform pathogenicity/virulence risk evaluations and generate extensive insights into the niche adaptation strategies through which they have evolved.

Antibiotic Resistance/Susceptibility Profiles of Staphylococcus equorum Strains from Cheese, and Genome Analysis for Antibiotic Resistance Genes

In food, bacteria carrying antibiotic resistance genes could play a prominent role in the spread of resistance. Staphylococcus equorum populations can become large in a number of fermented foods, yet the antibiotic resistance properties of this species have been little studied. In this work, the resistance/susceptibility (R/S) profile of S. equorum strains (n = 30) from cheese to 16 antibiotics was determined by broth microdilution. The minimum inhibitory concentration (MIC) for all antibiotics was low in most strains, although higher MICs compatible with acquired genes were also noted. Genome analysis of 13 strains showed the S. equorum resistome to be composed of intrinsic mechanisms, acquired mutations, and acquired genes. As such, a plasmidic cat gene providing resistance to chloramphenicol was found in one strain; this was able to provide resistance to Staphylococcus aureus after electroporation. An msr(A) polymorphic gene was identified in five strains. The Mrs(A) variants were associated with variable resistance to erythromycin. However, the genetic data did not always correlate with the phenotype. As such, all strains harbored a polymorphic fosB/fosD gene, although only one acquired copy was associated with strong resistance to fosfomycin. Similarly, a plasmid-associated blaR1-blaZI operon encoding a penicillinase system was identified in five ampicillin- and penicillin G-susceptible strains. Identified genes not associated with phenotypic resistance further included mph(C) in two strains and norA in all strains. The antibiotic R/S status and gene content of S. equorum strains intended to be employed in food systems should be carefully determined.

The rotation of primary starter culture mixtures results in batch-to-batch variations during Gouda cheese production

Industrial production of Gouda cheeses mostly relies on a rotated use of different mixed-strain lactic acid bacteria starter cultures to avoid phage infections. However, it is unknown how the application of these different starter culture mixtures affect the organoleptic properties of the final cheeses. Therefore, the present study assessed the impact of three different starter culture mixtures on the batch-to-batch variations among Gouda cheeses from 23 different batch productions in the same dairy company. Both the cores and rinds of all these cheeses were investigated after 36, 45, 75, and 100 weeks of ripening by metagenetics based on high-throughput full-length 16S rRNA gene sequencing accompanied with an amplicon sequence variant (ASV) approach as well as metabolite target analysis of non-volatile and volatile organic compounds. Up to 75 weeks of ripening, the acidifying Lactococcus cremoris and Lactococcus lactis were the most abundant bacterial species in the cheese cores. The relative abundance of Leuconostoc pseudomesenteroides was significantly different for each starter culture mixture. This impacted the concentrations of some key metabolites, such as acetoin produced from citrate, and the relative abundance of non-starter lactic acid bacteria (NSLAB). Cheeses with the least Leuc. pseudomesenteroides contained more NSLAB, such as Lacticaseibacillus paracasei that was taken over by Tetragenococcus halophilus and Loigolactobacillus rennini upon ripening time. Taken together, the results indicated a minor role of leuconostocs in aroma formation but a major impact on the growth of NSLAB. The relative abundance of T. halophilus (high) and Loil. rennini (low) increased with ripening time from rind to core. Two main ASV clusters of T. halophilus could be distinguished, which were differently correlated with some metabolites, both beneficial (regarding aroma formation) and undesirable ones (biogenic amines). A well-chosen T. halophilus strain could be a candidate adjunct culture for Gouda cheese production.

Dynamics of Starter and Non-Starter Lactic Acid Bacteria Populations in Long-Ripened Cheddar Cheese Using Propidium Monoazide (PMA) Treatment

The microbial community of industrially produced Canadian Cheddar cheese was examined from curd to ripened cheese at 30-32 months using a combination of viable plate counts of SLAB (GM17) and NSLAB (MRSv), qPCR and 16S rRNA gene amplicon sequencing. Cell treatment with propidium monoazide excluded DNA of permeable cells from amplification. The proportion of permeable cells of both Lactococcus spp. and Lacticaseibacillus spp. was highest at 3-6 months. While most remaining Lacticaseibacillus spp. cells were intact during later ripening stages, a consistent population of permeable Lactococcus spp. cells was maintained over the 32-month period. While Lactococcus sequence variants were significant biomarkers for viable cheese curd communities at 0-1 m, Lacticaseibacillus was identified as a distinctive biomarker for cheeses from 7 to 20 months. From 24 to 32 months, Lacticaseibacillus was replaced in significance by four genera (Pediococcus and Latilactobacillus at 24 m and at 30-32 m, Secundilactobacillus and Paucilactobacillus). These results underscore the importance of monitoring potential defects in cheeses aged over 24 months, which could be diagnosed early through microbial DNA profiling to minimize potential waste of product. Future perspectives include correlating volatile flavor compounds with microbial community composition as well as the investigation of intra-species diversity.

Needle in a Whey-Stack: PhRACS as a Discovery Tool for Unknown Phage-Host Combinations

The field of metagenomics has rapidly expanded to become the go-to method for complex microbial community analyses. However, there is currently no straightforward route from metagenomics to traditional culture-based methods of strain isolation, particularly in (bacterio)phage biology, leading to an investigative bottleneck. Here, we describe a method that exploits specific phage receptor binding protein (RBP)-host cell surface receptor interaction enabling isolation of phage-host combinations from an environmental sample. The method was successfully applied to two complex sample types-a dairy-derived whey sample and an infant fecal sample, enabling retrieval of specific and culturable phage hosts. IMPORTANCE PhRACS aims to bridge the current divide between in silico genetic analyses (i.e., phageomic studies) and traditional culture-based methodology. Through the labeling of specific bacterial hosts with fluorescently tagged recombinant phage receptor binding proteins and the isolation of tagged cells using flow cytometry, PhRACS allows the full potential of phageomic data to be realized in the wet laboratory.

Genetic diversity of Lactobacillus delbrueckii isolated from raw milk in Hokkaido, Japan

Lactic acid bacteria (LAB) play important roles in acid production and flavor formation in fermented dairy products. Lactic acid bacteria strains with distinct characteristics confer unique features to products. Diverse LAB have been identified in raw milk and traditional fermented milk prepared from raw milk. However, little is known about LAB in raw milk in Japan. To preserve diverse LAB as potential starters or probiotics for future use, we have isolated and identified various kinds of LAB from raw milk produced in Japan. In this study, we focused on Lactobacillus delbrueckii, one of the most important species in the dairy industry. We identified L. delbrueckii subspecies isolated from raw milk in Hokkaido, Japan, by analyzing intraspecific diversity using 4 distinct methods, hsp60 cluster analysis, multilocus sequence analysis, core-genome analysis, and whole-genome analysis based on average nucleotide identity. The subspecies distribution and a new dominant subset of L. delbrueckii from raw milk in Japan were revealed. The discovery of new strains with different genotypes is important for understanding the geographic distribution and characteristics of the bacteria and further their use as a microbial resource with the potential to express unconventional flavors and functionalities. The strains identified in this study may have practical applications in the development of fermented dairy products.

Massive Survey on Bacterial-Bacteriophages Biodiversity and Quality of Natural Whey Starter Cultures in Trentingrana Cheese Production

This study focused on the microbial and bacteriophages identification and characterization in cheese-production facilities that use natural whey starter (NWS) cultures for Trentingrana production. Bacterial and phage screening was carried out on cooked not acidified whey and NWS samples isolated from six dairy factories, for 4 consecutive days in four different months. By means of a combined approach, using plate counts, bacterial isolation, and metataxonomic analysis Lactobacillus helveticus was found occurring as the dominant species in NWS cultures and Levilactobacillus brevis as codominant in the cheese factories where the temperature of NWS production was mainly lower than 40°C, suggesting that the variability in the parameters of the NWS culture preparation could differently modulate the bacterial species in NWS cultures. Using turbidity test approach on 303 bacterial isolates from the NWS cultures, 120 distinct phages were identified. L. helveticus phage contamination of NWS cultures was revealed in most of the analyzed samples, but despite the great recovery of bacteriophage contamination cases, the microbial quality of NWS cultures was high. Our results support the presence of natural bacteriophage resistance mechanisms in L. helveticus. The use of NWS cultures probably creates an ideal environment for the proliferation of different L. helveticus strains balanced with their phages without a clear dominance. It is evident, from this study, that the presence of a high biodiversity of NWS bacterial strains is relevant to avoid phages dominance in NWS cultures and consequently to keep a good acidification ability.

Metataxonomic and metagenomic approaches for the study of undefined strain starters for cheese manufacture

Undefined strain starters are used for the production of many traditional and artisanal cheeses. Composition of undefined starters depends on several factors, and the diversity in strains and species significantly affects cheese quality and features. Culture-dependent approaches have long been used for the microbial profiling and functionalities of undefined cultures but underestimate their diversity due to culturability biases. Recently, culture-independent methods, based on high-throughput sequencing (HTS), have been preferred, with a significant boost in resolution power and sensitivity level. Amplicon targeted (AT) metagenomics, based on 16S rRNA sequencing, returned a larger microbiota diversity at genus and, sometimes, at species levels for artisanal starters of several PDO cheeses, but was inappropriate for populations with high strain diversity, and other gene targets were tested in AT approaches. Shotgun metagenomics (total DNA) and metatranscriptomics (total RNA), although are more powerful in depicting diversity and functionality of undefined cultures, have been rarely applied because of some limitations (e.g., high costs and laboriousness, need for bioinformatics skills). The advantages of HTS technologies are undoubted, but some hurdles need to be still overcame (e.g., resolution power, discrepancy between active and inactive cells, robust analytic pipelines, cost and time reduction for integrated approaches) so that HTS become routinary and convenient for defining complexity, microbial interactions (including host-phage relationships) and evolution in cheeses of undefined starters.

Complete Genome Sequence of Lactococcus lactis subsp. lactis bv. diacetylactis SD96

The genome of Lactococcus lactis subsp. lactis bv. diacetylactis SD96, a strain used for cheese production, is presented. SD96 is refractory to phage attack, which is a desired property for starter bacteria. Its 10 plasmids provide industrially important traits, such as lactose and citrate metabolism, proteolytic activity, and phage resistance.

Unprecedented Diversity of Lactococcal Group 936 Bacteriophages Revealed by Amplicon Sequencing of the Portal Protein Gene

Lactococcus lactis is one of the most important bacteria in dairy fermentations, being used in the production of cheese and buttermilk. The processes are vulnerable to phage attacks, and undefined mixtures of lactococcal strains are often used to reduce the risk of bacteriophage caused fermentation failure. Other preventive measures include culture rotation to prevent phage build-up and phage monitoring. Phage diversity, rather than quantity, is the largest threat to fermentations using undefined mixed starter cultures. We have developed a method for culture independent diversity analysis of lytic bacteriophages of the 936 group, the phages most commonly found in dairies. Using, as a target, a highly variable region of the portal protein gene, we demonstrate an unprecedented diversity and the presence of new 936 phages in samples taken from cheese production. The method should be useful to the dairy industry and starter culture manufacturers in their efforts to reduce phage problems.

A Diverse Range of Human Gut Bacteria Have the Potential To Metabolize the Dietary Component Gallic Acid

The human gut microbiota contains a broad variety of bacteria that possess functional genes, with resultant metabolites that affect human physiology and therefore health. Dietary gallates are phenolic components that are present in many foods and beverages and are regarded as having health-promoting attributes. However, the potential for metabolism of these phenolic compounds by the human microbiota remains largely unknown. The emergence of high-throughput sequencing (HTS) technologies allows this issue to be addressed. In this study, HTS was used to assess the incidence of gallate-decarboxylating bacteria within the gut microbiota of healthy individuals for whom bacterial diversity was previously determined to be high. This process was facilitated by the design and application of degenerate PCR primers to amplify a region encoding the catalytic C subunit of gallate decarboxylase (LpdC) from total metagenomic DNA extracted from human fecal samples. HTS resulted in the generation of a total of 3,261,967 sequence reads and revealed that the primary gallate-decarboxylating microbial phyla in the intestinal microbiota were Firmicutes (74.6%), Proteobacteria (17.6%), and Actinobacteria (7.8%). These reads corresponded to 53 genera, i.e., 47% of the bacterial genera detected previously in these samples. Among these genera, Anaerostipes and Klebsiella accounted for the majority of reads (40%). The usefulness of the HTS-lpdC method was demonstrated by the production of pyrogallol from gallic acid, as expected for functional gallate decarboxylases, among representative strains belonging to species identified in the human gut microbiota by this method.IMPORTANCE Despite the increasing wealth of sequencing data, the health contributions of many bacteria found in the human gut microbiota have yet to be elucidated. This study applies a novel experimental approach to predict the ability of gut microbes to carry out a specific metabolic activity, i.e., gallate metabolism. The study showed that, while gallate-decarboxylating bacteria represented 47% of the bacterial genera detected previously in the same human fecal samples, no gallate decarboxylase homologs were identified from representatives of Bacteroidetes The presence of functional gallate decarboxylases was demonstrated in representative Proteobacteria and Firmicutes strains from the human microbiota, an observation that could be of considerable relevance to the in vivo production of pyrogallol, a physiologically important bioactive compound.