Population Dynamics of Lactobacillus helveticus in Swiss Gruyère-Type Cheese Manufactured With Natural Whey Cultures

Intraspecific microdiversity and ecological drivers of lactic acid bacteria in naturally fermented milk ecosystem

Traditional fermented milks are produced by inoculating technique, which selects well-adapted microorganisms that have been passed on through generations. Few reports have used naturally fermented milks as model ecosystems to investigate the mechanism of formation of intra-species microbial diversity. Here, we isolated and whole-genome-sequenced a total of 717 lactic acid bacterial isolates obtained from 12 independent naturally fermented milks collect from 12 regions across five countries. We further analyzed the within-sample intra-species phylogenies of 214 Lactobacillus helveticus isolates, 97 Lactococcus lactis subsp. lactis isolates, and 325 Lactobacillus delbrueckii subsp. bulgaricus isolates. We observed a high degree of intra-species genomic and functional gene diversity within-/between-sample(s). Single nucleotide polymorphism-based phylogenetic reconstruction revealed great within-sample intra-species heterogeneity, evolving from multiple lineages. Further phylogenetic reconstruction (presence-absence gene profile) revealed within-sample inter-clade functional diversity (based on carbohydrate-active enzyme- and peptidase-encoding genes) in all three investigated species/subspecies. By identifying and mapping clade-specific genes of intra-sample clades of the three species/subspecies to the respective fermented milk metagenome, we found extensive potential inter-/intra-species horizontal gene transfer events. Finally, the microbial composition of the samples is closely linked to the nucleotide diversity of the respective species/subspecies. Overall, our results contribute to the conservation of lactic acid bacteria resources, providing ecological insights into the microbial ecosystem of naturally fermented dairy products.

Selective enrichment of the raw milk microbiota in cheese production: Concept of a natural adjunct milk culture

In cheese production, microorganisms are usually added at the beginning of the process as primary starters to drive curd acidification, while secondary microorganisms, with other pro-technological features important for cheese ripening, are added as selected cultures. This research aimed to investigate the possibilities of influencing and selecting the raw milk microbiota using artisanal traditional methods, providing a simple method to produce a natural supplementary culture. We investigated the production of an enriched raw milk whey culture (eRWC), a natural adjunct microbial culture produced from mixing an enriched raw milk (eRM) with a natural whey culture (NWC). The raw milk was enriched by spontaneous fermentation for 21 d at 10°C. Three milk enrichment protocols were tested: heat treatment before incubation, heat treatment plus salt addition, and no treatment. The eRMs were then co-fermented with NWC (ratio of 1:10) at 38°C for 6 h (young eRWC) and 22 h (old eRWC). Microbial diversity during cultures' preparation was evaluated through the determination of colony forming units on selective growth media, and next-generation sequencing (16S rRNA gene amplicon sequencing). The enrichment step increased the streptococci and lactobacilli but reduced microbial richness and diversity of the eRMs. Although the lactic acid bacteria viable count was not significantly different between the eRWCs, they harbored higher microbial richness and diversity than NWC. Natural adjunct cultures were then tested in cheese making trials, following the microbial development, and assessing the chemical quality of the 120 d ripened cheeses. The use of eRWCs slowed the curd acidification in the first hours of cheese making but the pH 24 h after production settled to equal values for all the cheeses. Although the use of diverse eRWCs contributed to having a richer and more diverse microbiota in the early stages of cheese making, their effect decreased over time during ripening, showing an inferior effect to the raw milk microbiota. Even if more research is needed, the optimization of such a tool could be an alternative to the practice of isolating, geno-pheno-typing, and formulating mixed-defined-strain adjunct cultures that require knowledge and facilities not always available for artisanal cheese makers.

Multifactorial Microvariability of the Italian Raw Milk Cheese Microbiota and Implication for Current Regulatory Scheme

Raw milk cheese manufactory is strictly regulated in Europe by the Protected Designation of Origin (PDO) quality scheme, which protects indigenous food products based on geographical and biotechnological features. This study encompassed the collection of 128 raw milk cheese samples across Italy to investigate the resident microbiome correlated to current PDO specifications. Shotgun metagenomic approaches highlighted how the microbial communities are primarily linked to each cheesemaking site and consequently to the use of site-specific Natural Whey Cultures (NWCs), defined by a multifactorial set of local environmental factors rather than solely by cheese type or geographical origin that guide the current PDO specification. Moreover, in-depth functional characterization of Cheese Community State Types (CCSTs) and comparative genomics efforts, including metagenomically assembled genomes (MAGs) of the dominant microbial taxa, revealed NWCs-related unique enzymatic profiles impacting the organoleptic features of the produced cheeses and availability of bioactive compounds to consumers, with putative health implications. Thus, these results highlighted the need for a profound rethinking of the current PDO designation with a focus on the production site-specific microbial metabolism to understand and guarantee the organoleptic features of the final product recognized as PDO. IMPORTANCE The Protected Designation of Origin (PDO) guarantees the traceability of food production processes, and that the production takes place in a well-defined restricted geographical area. Nevertheless, the organoleptic qualities of the same dairy products, i.e., cheeses under the same PDO denomination, differ between manufacturers. The final product's flavor and qualitative aspects can be related to the resident microbial population, not considered by the PDO denomination. Here, we analyzed a complete set of different Italian cheeses produced from raw milk through shotgun sequencing in order to study the variability of the different microbial profiles resident in Italian PDO cheeses. Furthermore, an in-depth functional analysis, along with a comparative genomic analysis, was performed in order to correlate the taxonomic information with the organoleptic properties of the final product. This analysis made it possible to highlight how the PDO denomination should be revisited to understand the effect that Natural Whey Cultures (NWCs), used in the traditional production of raw milk cheese and unique to each manufacturer, impacts on the organoleptic features of the final product.

Insights on the bacterial composition of Parmigiano Reggiano Natural Whey Starter by a culture-dependent and 16S rRNA metabarcoding portrait

Natural whey starters (NWS) are undefined bacterial communities produced daily from whey of the previous cheese-making round, by application of high temperature. As a result, in any dairy plant, NWS are continuously evolving, undefined mixtures of several strains and/or species of lactic acid bacteria, whose composition and performance strongly depend on the selective pressure acting during incubation. While NWS is critical to assure consistency to cheese-making process, little is known about the composition, functional features, and plant-to-plant fluctuations. Here, we integrated 16S rRNA metabarcoding and culture-dependent methods to profile bacterial communities of 10 NWS sampled in the production area of Parmigiano Reggiano cheese. 16S rRNA metabarcoding analysis revealed two main NWS community types, namely NWS type-H and NWS type-D. Lactobacillus helveticus was more abundant in NWS type-H, whilst Lactobacillus delbrueckii/St. thermophilus in NWS type-D, respectively. Based on the prediction of metagenome functions, NWS type-H samples were enriched in functional pathways related to galactose catabolism and purine metabolism, while NWS type-D in pathways related to aromatic and branched chain amino acid biosynthesis, which are flavor compound precursors. Culture-dependent approaches revealed low cultivability of individual colonies as axenic cultures and high genetic diversity in the pool of cultivable survivors. Co-culturing experiments showed that fermentative performance decreases by reducing the bacterial complexity of inoculum, suggesting that biotic interactions and cross-feeding relationships could take place in NWS communities, assuring phenotypic robustness. Even though our data cannot directly predict these ecological interactions, this study provides the basis for experiments targeted at understanding how selective regime affects composition, bacterial interaction, and fermentative performance in NWS.

Effect of using mycotoxin-detoxifying agents in dairy cattle feed on natural whey starter biodiversity

Natural whey cultures (NWC) are undefined multiple-strain bacterial starter communities that can be affected by even small changes along the entire dairy chain. We applied a multidisciplinary approach to investigate how the addition of 2 mycotoxin-detoxifying agents [sodium smectite and lignocellulose-based material (B1); leonardite and betaine (B2)] to cow diets modified the microbiota of the NWC in manufacture of a Grana-like cheese. Microbiological and flow cytometry analyses showed that the content and viability of lactic acid bacteria (LAB) and the total whey microbiota were not affected by the detoxifying agents, and Streptococcus thermophilus, Lactobacillus helveticus, and Limosilactobacillus fermentum were the dominant taxa. Random amplified polymorphic DNA-PCR fingerprinting and metagenomic analysis highlighted differences in the bacterial community of the NWC and in the relative abundance of Bacteroidetes that increased when B1 and B2 were included in the diet. Two of 6 St. thermophilus biotypes were detected only in control samples; conversely, none of the Lb. helveticus biotypes found in control samples were isolated from B1 and B2. In vitro tests showed that the 2 binders did not significantly affect the development of St. thermophilus, but they stimulated the growth of Lb. helveticus strains recovered only from B1 and B2 NWC. The addition of binders in cow feed can affect the LAB biotypes present in NWC.

Complete Genome Sequence of Lactobacillus helveticus JCM 1004, an Aminopeptidase-Producing Lactic Acid Bacterium

We report the complete genome sequence of Lactobacillus helveticus JCM 1004, an aminopeptidase-producing lactic acid bacterium. The genome consists of a circular chromosome which comprises 2,261,280 bp, with a G+C content of 37.56%. The genome was predicted to harbor 13 rRNA genes, 64 tRNA genes, and 2,462 protein-coding sequences.

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.

Microbial Interactions within the Cheese Ecosystem and Their Application to Improve Quality and Safety

The cheese microbiota comprises a consortium of prokaryotic, eukaryotic and viral populations, among which lactic acid bacteria (LAB) are majority components with a prominent role during manufacturing and ripening. The assortment, numbers and proportions of LAB and other microbial biotypes making up the microbiota of cheese are affected by a range of biotic and abiotic factors. Cooperative and competitive interactions between distinct members of the microbiota may occur, with rheological, organoleptic and safety implications for ripened cheese. However, the mechanistic details of these interactions, and their functional consequences, are largely unknown. Acquiring such knowledge is important if we are to predict when fermentations will be successful and understand the causes of technological failures. The experimental use of "synthetic" microbial communities might help throw light on the dynamics of different cheese microbiota components and the interplay between them. Although synthetic communities cannot reproduce entirely the natural microbial diversity in cheese, they could help reveal basic principles governing the interactions between microbial types and perhaps allow multi-species microbial communities to be developed as functional starters. By occupying the whole ecosystem taxonomically and functionally, microbiota-based cultures might be expected to be more resilient and efficient than conventional starters in the development of unique sensorial properties.

Probiotic potential of γ-aminobutyric acid (GABA)-producing yeast and its influence on the quality of cheese

Kazakh cheese is a traditional dairy product in Xinjiang, China. To study the function and potential probiotic characteristics of yeast in Kazakh cheese and its contribution to cheese fermentation, we screened the γ-aminobutyric acid (GABA)-producing yeasts Pichia kudriavzevii 1-21, Kluyveromyces marxianus B13-5, Saccharomyces cerevisiae DL6-20, and Kluyveromyces lactis DY1-10. We investigated the potential probiotic properties of these strains and their use in cheese fermentation (cheeses designated CSP, CSM, CSS, and CSI, respectively); a control with no added yeast was designated CS. The results showed that the 4 yeast strains all showed high self-polymerization (2- and 24-h autoaggregation capacity of >80 and 90%, respectively), hydrophobicity (40-92% variation, low hydrophobicity in xylene, but within the range of probiotics), and the ability to survive the gastrointestinal tract (survival rate >75% after simulation), indicating the probiotic ability of the strains in vitro. The GABA production capacity of the CSM cheese increased (to 95.6 mg/100 g), but its protein content did not change significantly, and amino acid degradation was obvious. The GABA production capacity of the CSS cheese decreased (to 450 mg/kg); its protein content declined, and its amino acid content increased. Except for water and protein, we found no obvious differences in most physical and chemical indicators. Kluyveromyces marxianus B13-5 helped to form the desired texture. Multivariate statistical analysis showed that fermentation of the cheese with the 4 yeasts improved the production of esters and alcohols. The CSS cheese had good aroma production performance, because S. cerevisiae DL6-20 produced high concentrations of isoamyl alcohol, hexanoic acid ethyl ester, benzyl alcohol, octanoic acid ethyl ester, 3-hydroxy-2-butanone, and hexanoic acid; the content of 2-methyl-propanoic acid was low. Compared with the CSP cheese, the CSI and CSM cheeses had a fruitier aroma and a milder odor, but the CSI and CSM cheeses had high concentrations of ethyl acetate, butanoic acid, ethyl ester, 3-methyl-1-butanol-acetate, ethyl hexanoate, ethyl octanoate, acetic acid 2-phenylethyl ester, and ethyl lactate; concentrations of 3-methyl-butanoic acid, propanoic acid, acetic acid, and butanoic acid were low. The CSP cheese had stronger acid-producing ability. The order of fragrance production performance was CSS > CSI, CSM > CSP > CS. Research into the fermentation mechanisms of GABA-producing yeast in cheese will provide a theoretical basis for the quality control and industrial production of Kazakh cheese.

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.

Development of a High-Throughput Microfluidic qPCR System for the Quantitative Determination of Quality-Relevant Bacteria in Cheese

The composition of the cheese microbiome has an important impact on the sensorial quality and safety of cheese. Therefore, much effort has been made to investigate the microbial community composition of cheese. Quantitative real-time polymerase chain reaction (qPCR) is a well-established method for detecting and quantifying bacteria. High-throughput qPCR (HT-qPCR) using microfluidics brings further advantages by providing fast results and by decreasing the cost per sample. We have developed a HT-qPCR approach for the rapid and cost-efficient quantification of microbial species in cheese by designing qPCR assays targeting 24 species/subspecies commonly found in cheese. Primer pairs were evaluated on the Biomark (Fluidigm) microfluidic HT-qPCR system using DNA from single strains and from artificial mock communities. The qPCR assays worked efficiently under identical PCR conditions, and the validation showed satisfying inclusivity, exclusivity, and amplification efficiencies. Preliminary results obtained from the HT-qPCR analysis of DNA samples of model cheeses made with the addition of adjunct cultures confirmed the potential of the microfluidic HT-qPCR system to screen for selected bacterial species in the cheese microbiome. HT-qPCR data of DNA samples of two downgraded commercial cheeses showed that this approach provides valuable information that can help to identify the microbial origin of quality defects. This newly developed HT-qPCR system is a promising approach that will allow simultaneous monitoring of quality-relevant species in fermented foods with high bacterial diversity, thereby opening up new perspectives for the control and assurance of high product quality.

Microbiota and Metabolite Profiling Combined With Integrative Analysis for Differentiating Cheeses of Varying Ripening Ages

Cheese maturation and flavor development results from complex interactions between milk substrates, cheese microbiota and their metabolites. In this study, bacterial 16S rRNA-gene sequencing, untargeted metabolomics (gas chromatography-mass spectrometry) and data integration analyses were used to characterize and differentiate commercial Cheddar cheeses of varying maturity made by the same and different manufacturers. Microbiota and metabolite compositions varied between cheeses of different ages and brands, and could be used to distinguish the cheeses. Individual amino acids and carboxylic acids were positively correlated with the ripening age for some brands. Integration and Random Forest analyses revealed numerous associations between specific bacteria and metabolites including a previously undescribed positive correlation between Thermus and phenylalanine and a negative correlation between Streptococcus and cholesterol. Together these results suggest that multi-omics analyses has the potential to be used for better understanding the relationships between cheese microbiota and metabolites during ripening and for discovering biomarkers for validating cheese age and brand authenticity.

Evaluation of bacterial communities of Grana Padano cheese by DNA metabarcoding and DNA fingerprinting analysis

The composition of the bacterial community of Grana Padano (GP) cheese was evaluated by an amplicon-based metagenomic approach (DNA metabarcoding) and RAPD-PCR fingerprinting. One hundred eighteen cheeses, which included 118 dairies located in the production area of GP, were collected. Two hundred fifty-four OTUs were detected, of which 82 were further discriminated between dominant (32 OTUs; > 1% total reads) and subdominant (50 OTUs; between 0.1% and 1% total reads) taxa. Lactobacillus (L.) delbrueckii, Lacticaseibacillus (Lact.) rhamnosus, Lact. casei, Limosilactobacillus fermentum, Lactococcus (Lc.) raffinolactis, L. helveticus, Streptococcus thermophilus, and Lc. lactis were the major dominant taxa ('core microbiota'). The origin of samples significantly impacted on both richness, evenness, and the relative abundance of bacterial species, with peculiar pattern distribution among the five GP production regions. A differential analysis allowed to find bacterial species significantly associated with specific region pairings. The analysis of pattern similarity among RAPD-PCR profiles highlighted the presence of a 'core' community banding pattern present in all the GP samples, which was strictly associated with the core microbiota highlighted by DNA metabarcoding. A trend to group samples according to the five production regions was also observed. This study widened our knowledge on the bacterial composition and ecology of Grana Padano cheese.

Ecology of Lactobacilli Present in Italian Cheeses Produced from Raw Milk

Among the bacterial genera that are used for cheese production, Lactobacillus is a key taxon of high industrial relevance that is commonly present in commercial starter cultures for dairy fermentations. Certain lactobacilli play a defining role in the development of the organoleptic features during the ripening stages of particular cheeses. We performed an in-depth 16S rRNA gene-based microbiota analysis coupled with internally transcribed spacer-mediated Lactobacillus compositional profiling of 21 common Italian cheeses produced from raw milk in order to evaluate the ecological distribution of lactobacilli associated with this food matrix. Statistical analysis of the collected data revealed the existence of putative Lactobacillus community state types (LCSTs), which consist of clusters of Lactobacillus (sub)species. Each LCST is dominated by one or two taxa that appear to represent keystone elements of an elaborate network of positive and negative interactions with minor components of the cheese microbiota. The results obtained in this study reveal the existence of peculiar cheese microbiota assemblies that represent intriguing targets for further functional studies aimed at dissecting the species-specific role of bacteria in cheese manufacturing.IMPORTANCE The microbiota is known to play a key role in the development of the organoleptic features of dairy products. Lactobacilli have been reported to represent one of the main components of the nonstarter bacterial population, i.e., bacteria that are not deliberately added to the milk, harbored by cheese, although the species-level composition of this microbial population has never been assessed in detail. In the present study, we applied a recently developed metagenomic approach that employs an internally transcribed spacer to profile the Lactobacillus population harbored by cheese produced from raw milk at the (sub)species level. The obtained data revealed the existence of particular Lactobacillus community state types consisting of clusters of Lactobacillus (sub)species that tend to cooccur in the screened cheeses. Moreover, analysis of covariances between members of this genus indicate that these taxa form an elaborate network of positive and negative interactions that define specific clusters of covariant lactobacilli.

Investigation of the Lactic Acid Bacteria in Kazak Cheese and Their Contributions to Cheese Fermentation

Kazak cheese is a traditional dairy product fermented by lactic acid bacteria (LAB) in Xinjiang. To investigate the LAB in Kazak cheese and their contributions to cheese fermentation, four representative LAB, Streptococcus thermophilus B8, Lactobacillus helveticus B6, Weissella confusa B14, and Lactobacillus rhamnosus B10, were isolated from Kazak cheese and subsequently used to ferment cheeses, which were named StC, LhC, WcC, and LrC, respectively. The result showed that most of the physical and chemical indicators had no significant difference, except for moisture and fat. W. confusa B14 was beneficial to the production of amino acids, whereas S. thermophilus B8 promoted the formation of organic acids and contributed to formation ideal texture property. Furthermore, the four cheeses all possessed a strong fruity aroma, with brandy, sweet, herbaceous, pungent, and fatty aromas being the most prominent in WcC. This is because L. helveticus B6 produced a high concentration of hexanal, nonanal, octanal, 3-methylbutanoic acid, ethyl acetate, ethyl butanoate, isoamyl acetate, and ethyl hexanoate in LhC. Research on the fermentation mechanism of LAB in cheese will provide a theoretical basis for the quality control and industrial production of Kazak cheese.

Safety and technological application of autochthonous Streptococcus thermophilus cultures in the buffalo Mozzarella cheese

Thermophilic and mesophilic lactic acid bacteria (LAB), such as Streptococcus thermophilus, Lactobacillus delbrueckii subsp. bulgaricus, Lactobacillus helveticus, and Lactococcus lactis, play a crucial role in the technological and sensory quality of Mozzarella cheese. In this study, the safety (genes encoding virulence factors and antibiotic resistance) and acidifying activity of autochthonous S. thermophilus cultures were evaluated in order to choose the most suitable strain for industrial application. The safe and good acidifying culture was tested in two buffalo Mozzarella cheese batches: Mozzarella cheeses produced with autochthonous culture (SJRP107) and commercial culture (STM5). The cultivable LAB was evaluated by culture-dependent method (plate counting) and the quantification of S. thermophilus cultures (commercial and autochthonous) were evaluated by culture-independent method RealT-qPCR (real-time quantitative polymerase chain reaction). The texture, physicochemical and proteolytic properties of the Mozzarella cheeses were similar for both batches. The nonstarter LAB count was higher during manufacture than in the storage, and the RealT-qPCR indicated the presence of S. thermophilus culture until the end of storage. S. thermophilus SJRP107 presented high potential for safety application in the production of Mozzarella cheese. Furthermore, considering the culture characteristics and their relationship with product quality, further studies could be helpful to determine their effect on the sensory characteristics of the cheese.

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.