Metagenomic Analysis of Dairy Bacteriophages: Extraction Method and Pilot Study on Whey Samples Derived from Using Undefined and Defined Mesophilic Starter Cultures

Lactococcal phage-host profiling through binding studies between cell wall polysaccharide types and Skunavirus receptor-binding proteins

Dairy fermentations using mesophilic starter cultures rely on the activity of specific lactic acid bacteria (LAB) such as Lactococcus lactis and Lactococcus cremoris for the acidification of milk. This biotechnological process can be affected by bacteriophage infection of LAB starter strains, which may result in delayed or even failed fermentations. Most studied lactococcal phages commence infection with the binding of a tail-associated receptor-binding protein (RBP) to a host cell surface-exposed cell wall polysaccharide (CWPS). In the present study, phage prevalence and diversity in whey samples originating from fermentations performed in various European countries employing undefined mesophilic starter cultures were investigated using phageome analysis. The range of Skunavirus RBP genotypes present in the phageomes and associated RBP-CWPS binding abilities were evaluated, resulting in the refinement and expansion of the Skunavirus RBP grouping system and the identification of several heretofore unknown Skunavirus RBP (sub)groups. These findings substantially expand our knowledge on lactococcal Skunavirus RBP diversity and their binding specificity towards CWPS receptor structures, thereby improving the predictability of fermentation outcomes and robustness of starter culture rotations and blends.

Dynamics of the viral community on the surface of a French smear-ripened cheese during maturation and persistence across production years

The surface of smear-ripened cheeses constitutes a dynamic microbial ecosystem resulting from the successive development of different microbial groups such as lactic acid bacteria, fungi, and ripening bacteria. Recent studies indicate that a viral community, mainly composed of bacteriophages, also represents a common and substantial part of the cheese microbiome. However, the composition of this community, its temporal variations, and associations between bacteriophages and their hosts remain poorly characterized. Here, we studied a French smear-ripened cheese by both viral metagenomics and 16S metabarcoding approaches to assess both the succession of phages and bacterial communities on the cheese surface during cheese ripening and their temporal variations in ready-to-eat cheeses over the years of production. We observed a clear transition of the phage community structure during ripening with a decreased relative abundance of viral species (vOTUs) associated with Lactococcus phages, which were replaced by vOTUs associated with phages infecting ripening bacteria such as Brevibacterium, Glutamicibacter, Pseudoalteromonas, and Vibrio. The dynamics of the phage community was strongly associated with bacterial successions observed on the cheese surface. Finally, while some variations in the distribution of phages were observed in ready-to-eat cheeses produced at different dates spanning more than 4 years of production, the most abundant phages were detected throughout. This result revealed the long-term persistence of the dominant phages in the cheese production environment. Together, these findings offer novel perspectives on the ecology of bacteriophages in smear-ripened cheese and emphasize the significance of incorporating bacteriophages in the microbial ecology studies of fermented foods.IMPORTANCEThe succession of diverse microbial populations is critical for ensuring the production of high-quality cheese. We observed a temporal succession of phages on the surface of a smear-ripened cheese, with new phage communities showing up when ripening bacteria start covering this surface. Interestingly, the final phage community of this cheese is also consistent over large periods of time, as the same bacteriophages were found in cheese products from the same manufacturer made over 4 years. This research highlights the importance of considering these bacteriophages when studying the microbial life of fermented foods like cheese.

The bacteriome-coupled phage communities continuously contract and shift to orchestrate the traditional rice vinegar fermentation

Amounts of microbiome studies have uncovered the microbial communities of traditional food fermentations, while in which the phageome development with time is poorly understood. Here, we conducted a study to decipher both phageome and bacteriome of the traditional rice vinegar fermentation. The vinegar phageomes showed significant differences in the alpha diversity, network density and clustering coefficient over time. Peduoviridae had the highest relative abundance. Moreover, the phageome negatively correlated to the cognate bacteriome in alpha diversity, and undergone constantly contracting and shifting across the temporal scale. Nevertheless, 257 core virial clusters (VCs) persistently occurred with time whatever the significant impacts imposed by the varied physiochemical properties. Glycoside hydrolase (GH) and glycosyltransferase (GT) families genes displayed the higher abundances across all samples. Intriguingly, diversely structuring of toxin-antitoxin systems (TAs) and CRISPR-Cas arrays were frequently harbored by phage genomes. Their divergent organization and encoding attributes underlie the multiple biological roles in modulation of network and/or contest of phage community as well as bacterial host community. This phageome-wide mapping will fuel the current insights of phage community ecology in other traditional fermented ecosystems that are challenging to decipher.

The insights into the phage communities of fermented foods in the age of viral metagenomics

Phages play a critical role in the assembly and regulation of fermented food microbiome through lysis and lysogenic lifestyle, which in turn affects the yield and quality of fermented foods. Therefore, it is important to investigate and characterize the diversity and function of phages under complex microbial communities and nutrient substrate conditions to provide novel insights into the regulation of traditional spontaneous fermentation. Viral metagenomics has gradually garnered increasing attention in fermented food research to elucidate phage functions and characterize the interactions between phages and the microbial community. Advances in this technology have uncovered a wide range of phages associated with the production of traditional fermented foods and beverages. This paper reviews the common methods of viral metagenomics applied in fermented food research, and summarizes the ecological functions of phages in traditional fermented foods. In the future, combining viral metagenomics with culturable methods and metagenomics will broaden the scope of research on fermented food systems, revealing the complex role of phages and intricate phage-bacterium interactions.

Identifying ARG-carrying bacteriophages in a lake replenished by reclaimed water using deep learning techniques

As important mobile genetic elements, phages support the spread of antibiotic resistance genes (ARGs). Previous analyses of metaviromes or metagenome-assembled genomes (MAGs) failed to assess the extent of ARGs transferred by phages, particularly in the generation of antibiotic pathogens. Therefore, we have developed a bioinformatic pipeline that utilizes deep learning techniques to identify ARG-carrying phages and predict their hosts, with a special focus on pathogens. Using this method, we discovered that the predominant types of ARGs carried by temperate phages in a typical landscape lake, which is fully replenished by reclaimed water, were related to multidrug resistance and β-lactam antibiotics. MAGs containing virulent factors (VFs) were predicted to serve as hosts for these ARG-carrying phages, which suggests that the phages may have the potential to transfer ARGs. In silico analysis showed a significant positive correlation between temperate phages and host pathogens (R = 0.503, p < 0.001), which was later confirmed by qPCR. Interestingly, these MAGs were found to be more abundant than those containing both ARGs and VFs, especially in December and March. Seasonal variations were observed in the abundance of phages harboring ARGs (from 5.62 % to 21.02 %) and chromosomes harboring ARGs (from 18.01 % to 30.94 %). In contrast, the abundance of plasmids harboring ARGs remained unchanged. In summary, this study leverages deep learning to analyze phage-transferred ARGs and demonstrates an alternative method to track the production of potential antibiotic-resistant pathogens by metagenomics that can be extended to microbiological risk assessment.

A Metagenomics Approach to Enumerate Bacteriophages in a Food Niche

Dairy fermentation relies on the activity of starter cultures composed primarily of lactic acid bacteria such as Lactococcus and Streptococcus thermophilus strains to produce consistent, high-quality products. Bacteriophages are a constant threat to the industry, often causing slowed or failed fermentation resulting in significant economic losses. To ensure the continuation of reliable fermentation practices, it is important to detect and monitor the phage populations impacting different starter cultures. This has traditionally been done primarily through culture-dependent methods but has since expanded into viral metagenomics. Here we outline a protocol for a targeted virome extraction from a dairy whey sample, followed by subsequent sequencing and phageome analysis of the sample.

Longitudinal Study of Lactococcus Phages in a Canadian Cheese Factory

The presence of virulent phages is closely monitored during cheese manufacturing, as these bacterial viruses can significantly slow down the milk fermentation process and lead to low-quality cheeses. From 2001 to 2020, whey samples from cheddar cheese production in a Canadian factory were monitored for the presence of virulent phages capable of infecting proprietary strains of Lactococcus cremoris and Lactococcus lactis used in starter cultures. Phages were successfully isolated from 932 whey samples using standard plaque assays and several industrial Lactococcus strains as hosts. A multiplex PCR assay assigned 97% of these phage isolates to the Skunavirus genus, 2% to the P335 group, and 1% to the Ceduovirus genus. DNA restriction profiles and a multilocus sequence typing (MLST) scheme distinguished at least 241 unique lactococcal phages from these isolates. While most phages were isolated only once, 93 of them (out of 241, 39%) were isolated multiple times. Phage GL7 was isolated 132 times from 2006 to 2020, demonstrating that phages can persist in a cheese factory for long periods of time. Phylogenetic analysis of MLST sequences showed that phages could be clustered based on their bacterial hosts rather than their year of isolation. Host range analysis showed that Skunavirus phages exhibited a very narrow host range, whereas some Ceduovirus and P335 phages had a broader host range. Overall, the host range information was useful in improving the starter culture rotation by identifying phage-unrelated strains and helped mitigating the risk of fermentation failure due to virulent phages. IMPORTANCE Although lactococcal phages have been observed in cheese production settings for almost a century, few longitudinal studies have been performed. This 20-year study describes the close monitoring of dairy lactococcal phages in a cheddar cheese factory. Routine monitoring was conducted by factory staff, and when whey samples were found to inhibit industrial starter cultures under laboratory conditions, they were sent to an academic research laboratory for phage isolation and characterization. This led to a collection of at least 241 unique lactococcal phages, which were characterized through PCR typing and MLST profiling. Phages of the Skunavirus genus were by far the most dominant. Most phages lysed a small subset of the Lactococcus strains. These findings guided the industrial partner in adapting the starter culture schedule by using phage-unrelated strains in starter cultures and removing some strains from the starter rotation. This phage control strategy could be adapted for other large-scale bacterial fermentation processes.

Integrated molecular approaches for fermented food microbiome research

Molecular technologies, including high-throughput sequencing, have expanded our perception of the microbial world. Unprecedented insights into the composition and function of microbial communities have generated large interest, with numerous landmark studies published in recent years relating the important roles of microbiomes and the environment-especially diet and nutrition-in human, animal, and global health. As such, food microbiomes represent an important cross-over between the environment and host. This is especially true of fermented food microbiomes, which actively introduce microbial metabolites and, to a lesser extent, live microbes into the human gut. Here, we discuss the history of fermented foods, and examine how molecular approaches have advanced research of these fermented foods over the past decade. We highlight how various molecular approaches have helped us to understand the ways in which microbes shape the qualities of these products, and we summarize the impacts of consuming fermented foods on the gut. Finally, we explore how advances in bioinformatics could be leveraged to enhance our understanding of fermented foods. This review highlights how integrated molecular approaches are changing our understanding of the microbial communities associated with food fermentation, the creation of unique food products, and their influences on the human microbiome and health.

High Level of Interaction between Phages and Bacteria in an Artisanal Raw Milk Cheese Microbial Community

Microbial starter cultures are used in the production of many cheeses around the world, such as Parmigiano-Reggiano, in Italy, Époisses, in France, and Canastra, in Brazil, providing many of the unique features of these cheeses. Bacteriophages (phages) are ubiquitous and well known to modulate the structure of bacterial communities, and recent data indicate that cheeses contain a high abundance of naturally occurring phages. Here, we analyze the viral and bacterial metagenomes of Canastra cheese: a traditional artisanal Brazilian cheese produced using an endogenous starter culture and raw milk. Over 1,200 viral operational taxonomic units were recovered using both isolated viral-like particles and complete metagenomic DNA. Common viral families identified included Siphoviridae and Myoviridae, with 40% of putative phage genomes unidentified at the family level of classification. We observed very high phage diversity, which varied greatly across different cheese producers, with 28% of phage genomes detected in only one producer. Several metagenome-assembled genomes were recovered for lactic acid-producing bacteria, as well as nonstarter bacterial species, and we identified several phage-bacterium interactions, at the strain level of resolution, varying across distinct cheese producers. We postulate that at least one bacterial strain detected could be endogenous and unique to the Canastra cheese-producing region in Brazil and that its growth seems to be modulated by autochthonous phages present in this artisanal production system. This phage-host relationship is likely to influence the fermentation dynamics and ultimately the sensorial profile of these cheeses, with implications for other similar cheese production systems around the world. IMPORTANCE Our work demonstrated a dynamic yet stable microbial ecosystem during cheese production using an endogenous starter culture. This was observed across several distinct producers and was marked by genomic evidence of continued phage-bacterium interactions, such as the presence of bacterial defense mechanisms. Furthermore, we provide evidence of unique microbial signatures for each individual cheese producer studied in the region, a fact that may have profound consequences on product traceability. This was the first effort to describe and understand the bacteriophage composition and ecological dynamics within the Brazilian Canastra cheese production system. The study of this prototypical backslopping production system provides a solid background for further mechanistic studies of the production of many cheeses around the world.

The Application of Metagenomics to Study Microbial Communities and Develop Desirable Traits in Fermented Foods

The microbial communities present within fermented foods are diverse and dynamic, producing a variety of metabolites responsible for the fermentation processes, imparting characteristic organoleptic qualities and health-promoting traits, and maintaining microbiological safety of fermented foods. In this context, it is crucial to study these microbial communities to characterise fermented foods and the production processes involved. High Throughput Sequencing (HTS)-based methods such as metagenomics enable microbial community studies through amplicon and shotgun sequencing approaches. As the field constantly develops, sequencing technologies are becoming more accessible, affordable and accurate with a further shift from short read to long read sequencing being observed. Metagenomics is enjoying wide-spread application in fermented food studies and in recent years is also being employed in concert with synthetic biology techniques to help tackle problems with the large amounts of waste generated in the food sector. This review presents an introduction to current sequencing technologies and the benefits of their application in fermented foods.

Dairy streptococcal cell wall and exopolysaccharide genome diversity

The large-scale and high-intensity application of Streptococcus thermophilus species in milk fermentation processes is associated with a persistent threat of (bacterio)phage infection. Phage infection of starter cultures may cause inconsistent, slow or even failed fermentations with consequent diminished product quality and/or output. The phage life cycle commences with the recognition of, and binding to, a specific host-encoded and surface-exposed receptor, which in the case of S. thermophilus can be the rhamnose-glucose polysaccharide (RGP; specified by the rgp gene cluster) or exopolysaccharide (EPS; specified by the eps gene cluster). The genomic diversity of 23 S . thermophilus strains isolated from unpasteurized dairy products was evaluated, including a detailed analysis of the rgp and eps loci. In the present study, five novel eps genotypes were identified while variations of currently recognized rgp gene cluster types were also observed. Furthermore, the diversity of rgp genotypes amongst retrieved isolates positively correlated with phage diversity based on phageome analysis of eight representative dairy products. Our findings therefore substantially expand our knowledge on S. thermophilus’ strain and phage diversity in (artisanal) dairy products and highlight the merit of phageome analysis of artisanal and traditional fermented foods as a sensitive marker of dominant microbiota involved in the fermentation.

An Adenosine Triphosphate- Dependent 5'-3' DNA Helicase From sk1-Like Lactococcus lactis F13 Phage

Here, we describe functional characterization of an early gene (gp46) product of a virulent Lactococcus lactis sk1-like phage, vB_Llc_bIBBF13 (abbr. F13). The GP46_F13 protein carries a catalytically active RecA-like domain belonging to the P-loop NTPase superfamily. It also retains features characteristic for ATPases forming oligomers. In order to elucidate its detailed molecular function, we cloned and overexpressed the gp46 gene in Escherichia coli. Purified GP46_F13 protein binds to DNA and exhibits DNA unwinding activity on branched substrates in the presence of adenosine triphosphate (ATP). Size exclusion chromatography with multi-angle light scattering (SEC-MALS) experiments demonstrate that GP46_F13 forms oligomers, and further pull-down assays show that GP46_F13 interacts with host proteins involved in replication (i.e., DnaK, DnaJ, topoisomerase I, and single-strand binding protein). Taking together the localization of the gene and the obtained results, GP46_F13 is the first protein encoded in the early-expressed gene region with helicase activity that has been identified among lytic L. lactis phages up to date.

Virome studies of food production systems: time for 'farm to fork' analyses

The food industry is under increasing pressure to produce high quality, traceable and minimally processed foods that are produced using sustainable approaches and ingredients. In line with the latter, there is an increased pressure for plant-based products to replace animal-derived products. Until recently, research efforts have mainly focused on dairy and meat products owing to their economic importance. The shift towards plant-based diets and food production requires a corresponding shift in research efforts to define the microbial requirements for and composition of (novel) plant-based foods, the (micro)organisms that are beneficial to such production systems, and the abundance and role of (bacterio)phages in shaping the microbial landscape of these foods. In this review, we explore current efforts in the area of virome analysis of foods and food production environments and highlight the need for more unified approaches to understand the contribution of phages in food safety and quality, and to develop novel tools to enhance the traceability of foods.

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.

Bacteriophage-mediated manipulation of the gut microbiome – promises and presents limitations

Gut microbiome (GM) composition and function are linked to human health and disease, and routes for manipulating the GM have become an area of intense research. Due to its high treatment efficacy, the use of fecal microbiota transplantation (FMT) is generally accepted as a promising experimental treatment for patients suffering from GM imbalances (dysbiosis), e.g. caused by recurrent Clostridioides difficile infections (rCDI). Mounting evidence suggests that bacteriophages (phages) play a key role in successful FMT treatment by restoring the dysbiotic bacterial GM. As a refinement to FMT, removing the bacterial component of donor feces by sterile filtration, also referred to as fecal virome transplantation (FVT), decreases the risk of invasive infections caused by bacteria. However, eukaryotic viruses and prophage-encoded virulence factors remain a safety issue. Recent in vivo studies show how cascading effects are initiated when phage communities are transferred to the gut by e.g. FVT, which leads to changes in the GM composition, host metabolome, and improve host health such as alleviating symptoms of obesity and type-2-diabetes (T2D). In this review, we discuss the promises and limitations of FVT along with the perspectives of using FVT to treat various diseases associated with GM dysbiosis.

A Specific Sugar Moiety in the Lactococcus lactis Cell Wall Pellicle Is Required for Infection by CHPC971, a Member of the Rare 1706 Phage Species

Lactococcus lactis is a Gram-positive bacterium widely used as a starter culture for the production of different dairy products, especially a large variety of cheeses. Infection of lactococcal starter cultures by bacteriophages is one of the major causes of fermentation failure and often leads to production halt. Lactococcal bacteriophages belonging to the c2, 936, and P335 species are the most commonly isolated in dairy plants and have been extensively investigated in the past three decades. Information regarding bacteriophages belonging to less commonly isolated species is, on the other hand, less extensive, although these phages can also contribute to starter culture infection. Here, we report the nucleotide sequence of the newly isolated L. lactis phage CHPC971, belonging to the rare 1706 species of lactococcal phages. We investigated the nature of the host receptor recognized by the phage and collected evidence that strongly suggests that it binds to a specific sugar moiety in the cell wall pellicle of its host. An in silico analysis of the genome of phage CHPC971 identified the hypothetical genes involved in receptor binding.IMPORTANCE Gathering information on how lactococcal bacteriophages recognize their host and proliferate in the dairy environment is of vital importance for the establishment of proper starter culture rotation plans and to avoid fermentation failure and consequent great economic losses for dairy industries. We provide strong evidence on the type of receptor recognized by a newly isolated 1706-type lactococcal bacteriophage, increasing knowledge of phage-host interactions relevant to dairying. This information can help to prevent phage infection events that, so far, are hard to predict and avoid.

On the Role of Bioinformatics and Data Science in Industrial Microbiome Applications

Advances in sequencing and computational biology have drastically increased our capability to explore the taxonomic and functional compositions of microbial communities that play crucial roles in industrial processes. Correspondingly, commercial interest has risen for applications where microbial communities make important contributions. These include food production, probiotics, cosmetics, and enzyme discovery. Other commercial applications include software that takes the user's gut microbiome data as one of its inputs and outputs evidence-based, automated, and personalized diet recommendations for balanced blood sugar levels. These applications pose several bioinformatic and data science challenges that range from requiring strain-level resolution in community profiles to the integration of large datasets for predictive machine learning purposes. In this perspective, we provide our insights on such challenges by touching upon several industrial areas, and briefly discuss advances and future directions of bioinformatics and data science in microbiome research.

A comparative genomics approach for identifying host-range determinants in Streptococcus thermophilus bacteriophages

Comparative genomics has proven useful in exploring the biodiversity of phages and understanding phage-host interactions. This knowledge is particularly useful for phages infecting Streptococcus thermophilus, as they constitute a constant threat during dairy fermentations. Here, we explore the genetic diversity of S. thermophilus phages to identify genetic determinants with a signature for host specificity, which could be linked to the bacterial receptor genotype. A comparative genomic analysis was performed on 142 S. thermophilus phage genomes, 55 of which were sequenced in this study. Effectively, 94 phages were assigned to the group cos (DT1), 36 to the group pac (O1205), six to the group 5093, and six to the group 987. The core genome-based phylogeny of phages from the two dominating groups and their receptor binding protein (RBP) phylogeny corresponded to the phage host-range. A role of RBP in host recognition was confirmed by constructing a fluorescent derivative of the RBP of phage CHPC951, followed by studying the binding of the protein to the host strain. Furthermore, the RBP phylogeny of the cos group was found to correlate with the host genotype of the exocellular polysaccharide-encoding operon. These findings provide novel insights towards developing strategies to combat phage infections in dairies.

Bacteriophages in the Dairy Environment: From Enemies to Allies

The history of dairy farming goes back thousands of years, evolving from a traditional small-scale production to the industrialized manufacturing of fermented dairy products. Commercialization of milk and its derived products has been very important not only as a source of nourishment but also as an economic resource. However, the dairy industry has encountered several problems that have to be overcome to ensure the quality and safety of the final products, as well as to avoid economic losses. Within this context, it is interesting to highlight the role played by bacteriophages, or phages, viruses that infect bacteria. Indeed, bacteriophages were originally regarded as a nuisance, being responsible for fermentation failure and economic losses when infecting lactic acid bacteria, but are now considered promising antimicrobials to fight milk-borne pathogens without contributing to the increase in antibiotic resistance.