Submitted Successfully!
To reward your contribution, here is a gift for you: A free trial for our video production service.
Thank you for your contribution! You can also upload a video entry or images related to this topic.
Version Summary Created by Modification Content Size Created at Operation
1 + 5526 word(s) 5526 2022-01-17 09:46:33 |
2 format corrected. Meta information modification 5526 2022-01-18 02:57:01 |

Video Upload Options

We provide professional Video Production Services to translate complex research into visually appealing presentations. Would you like to try it?

Confirm

Are you sure to Delete?
Cite
If you have any further questions, please contact Encyclopedia Editorial Office.
Anastasiou, R. Omics Approaches to Assess Flavor Development in Cheese. Encyclopedia. Available online: https://encyclopedia.pub/entry/18342 (accessed on 19 November 2024).
Anastasiou R. Omics Approaches to Assess Flavor Development in Cheese. Encyclopedia. Available at: https://encyclopedia.pub/entry/18342. Accessed November 19, 2024.
Anastasiou, Rania. "Omics Approaches to Assess Flavor Development in Cheese" Encyclopedia, https://encyclopedia.pub/entry/18342 (accessed November 19, 2024).
Anastasiou, R. (2022, January 17). Omics Approaches to Assess Flavor Development in Cheese. In Encyclopedia. https://encyclopedia.pub/entry/18342
Anastasiou, Rania. "Omics Approaches to Assess Flavor Development in Cheese." Encyclopedia. Web. 17 January, 2022.
Omics Approaches to Assess Flavor Development in Cheese
Edit

Cheese is characterized by a rich and complex microbiota that plays a vital role during both production and ripening, contributing significantly to the safety, quality, and sensory characteristics of the final product. In this context, it is vital to explore the microbiota composition and understand its dynamics and evolution during cheese manufacturing and ripening. Application of high-throughput DNA sequencing technologies have facilitated the more accurate identification of the cheese microbiome, detailed study of its potential functionality, and its contribution to the development of specific organoleptic properties. These technologies include amplicon sequencing, whole-metagenome shotgun sequencing, metatranscriptomics, and, most recently, metabolomics. In recent years, however, the application of multiple meta-omics approaches along with data integration analysis, which was enabled by advanced computational and bioinformatics tools, paved the way to better comprehension of the cheese ripening process, revealing significant associations between the cheese microbiota and metabolites, as well as their impact on cheese flavor and quality. 

cheese flavor omics lactic acid bacteria yeasts cheese microbiome metagenomics metatranscriptomics metabolomics

1. Omics Insights into Flavor Formation in Cheese

1.1. Genomics

The first omics technique to appear, genomics, focused on sequencing and annotation of the complete genome of an organism in order to reveal the genetic structure and to predict gene functions that are encoded within the genome. Since the publication of the long-awaited whole genome sequence of the first LAB strain in 2001, namely L. lactis subsp. lactis IL1403 [1], the number of sequenced LAB genomes has grown exponentially, due to the advent of HTS, with the genome sequences of thousands of LAB species and strains currently being deposited in public databases [2].
The study of Bolotin et al. [1] was the first genomic milestone in LAB research and laid the groundwork for a better understanding of the many aspects of bacterial physiology, metabolic pathways, and regulatory mechanisms, as well as how phenotypes are affected by genetic variations. In the following years, numerous genomics studies were conducted regarding strains of the “laboratory workhorse” species L. lactis [3][4][5][6][7]. As mentioned before, S. thermophilus, L. bulgaricus, and L. helveticus are also the main SLAB used in cheese production along with L. lactis. Therefore, due to their industrial and economic importance, efforts are constantly being made to fully sequence and annotate genomes of these taxa [8][9][10][11]. In recent years though, genomics studies included NSLAB as well, due to their considerable impact on flavor development in cheese [12][13][14][15][16][17][18][19][20][21][22][23].

Proteins Associated with Flavor Formation in Cheese

Proteolysis, lipolysis, and AAs/FAs catabolism have all been studied in detail before whole genome sequencing (WGS) appeared. However, the available genome sequences allowed an in-depth analysis of the respective genes, so as to understand the genetic instability of several traits and unravel strain-specific differences. From this perspective, a comparative genomics analysis of all the proteins involved in the proteolytic system of 22 completely sequenced LAB strains was performed, including the cell-wall-bound proteinase, peptide transporters, and peptidases [24]. Based on the results, the genomes of the Lactobacillus acidophilus, Lactobacillus johnsonii, Lactobacillus gasseri, L. bulgaricus, and L. helveticus strains studied encoded a relatively higher number of proteolysis-related genes. Furthermore, the cell-wall-bound proteinase PrtP was solely identified in the chromosomes of L. acidophilus, L. johnsonii, L. bulgaricus, L. casei, L. rhamnosus, and S. thermophilus strains, as well as in the plasmid of L. lactis subsp. cremoris SK11. On the other hand, endopeptidases PepE/PepG and proline peptidases PepI, PepR, and PepL were absent in Lactococcus and Streptococcus strains, while aminopeptidases PepC, PepN, and PepM and proline peptidases PepX and PepQ were present in all genomes analyzed [24]. More recently, comparative genomics analysis among 213 assemblies, of which 175 belonged to Lactobacillus species and 38 to associated genera, was performed [25]. Regarding the metabolic potential of the 213 strains analyzed, the authors found genes for 60 cell envelope proteinases, which are important for cleaving casein during growth in milk and thus, contribute to cheese flavor, ranging in length from 1097 to 2270 amino acids. In addition, a broad repertoire of glycoside hydrolases and glycosyltransferases was identified, which are both important in carbohydrate metabolism [25]. Therefore, genome-scale metabolic models have been constructed and applied for the in silico prediction of the metabolic patterns of LAB strains under various conditions [26][27][28][29][30][31]. It should be noted, however, that the accuracy of these models depends on the quality of the genome sequencing and the correct annotation.
Furthermore, genomes of the non-LAB genera Propionibacterium and Brevibacterium have also been sequenced and annotated, as they also contribute to the flavor of certain types of cheeses. Genomics analysis of Propionibacterium assemblies mainly focuses on the identification of genes that are involved in the two key metabolic pathways for the propionate production, i.e., the Wood–Werkman and the tricarboxylic acid cycles, the amino acid catabolic pathways, which result in the formation of volatile compounds, and the detection of esterases involved in the formation of FFAs and esters [32][33][34]. Similarly, WGS studies of Brevibacterium assemblies reported, among others, the gene repertoire responsible for the catabolism of lactose, galactose, citrate, lipids, proteins, and amino acids, which contribute to the flavor, texture, and appearance of cheeses [35][36][37][38].
Apart from bacteria, several fungal species are also used as starters in internally and surface-ripened cheeses, such as Roquefort and Camembert, respectively, with a vital role in the flavor formation of the final product. Albeit their importance, there is still a limited number of sequenced genomes in the NCBI database, e.g., only five and six partially sequenced assemblies for P. roqueforti and P. camemberti, respectively, which are two of the main fungal species used in cheese production. However, most of the WGS studies performed focused on the phylogeny and not the metabolic pathways regarding flavor formation [39][40][41][42][43].
Linking genotype to phenotype is of paramount importance to an in-depth understanding of the technological potential of a strain and therefore, to a better selection of candidates with flavor-forming metabolic potential to be used as starter or adjunct cultures in cheese production. However, as cheese ripening depends on a complex microbial community, the need for metagenomics analysis quickly arose.

1.2. Metagenomics

Metagenomics encompasses two different HTS approaches, namely amplicon sequencing and shotgun metagenome sequencing. In amplicon sequencing, a highly conserved marker gene or genome’s region of DNA directly extracted from a food microbial community, e.g., that of cheese, is amplified and sequenced. These markers are of taxonomic relevance, with the 16S rRNA gene being used for the identification of bacteria taxa and the 18S rRNA gene together with the ITS DNA region for the yeasts/fungal taxa. On the other hand, in shotgun metagenomics, the extracted DNA is fully sequenced in a non-targeted manner, thus providing not only taxonomical identification results, but also information on the metabolic potential of the microbial community by reconstructing metabolic pathways [44].
From this perspective, a shotgun metagenomics study was recently performed in 25 Cotija cheese samples [45]. This artisanal Mexican cheese is produced by raw milk without the addition of starter cultures and is ripened in an open environment. Therefore, the organoleptic characteristics of Cotija cheese depend on a range of biotic and abiotic factors. Taxonomic identification revealed that the bacterial microbiota is mainly composed of Firmicutes, followed by Actinobacteria and Proteobacteria. The authors were able to reconstruct the genome assemblies of the three dominant species detected, i.e., L. plantarum, Leuc. mesenteroides, and W. paramesenteroides, which accounted for more than 80% of the total bacterial sequences. Gene functional annotation related to Cotija cheese flavor resulted in the identification of genes involved in the catabolism of phenylalanine, branched-chain amino acids, and fatty acids [45].

1.3. Transcriptomics

Genomics is widely applied to study the technological potential of a strain; however, it can only in silico predict its functional potential. Therefore, transcriptomics (RNA-seq), although a relatively new field of research, has gained much attention as it provides an accurate profile of a strain’s functional activity at a given time point.
Regarding cheese microbiota, the majority of the transcriptomic studies have been performed on L. lactis, due to its importance in cheese production [2][46][47][48][49][50][51][52][53][54]. In particular, the transcriptional responses of L. lactis have been investigated during different stresses, such as cold, heat, acid, osmotic, oxidative, and starvation [55][56][57][58][59][60][61], or during growth in media with different carbon sources [62][63][64]. Apart from L. lactis, transcriptomics studies have also been performed in other important cheese bacterial [48][65][66][67][68][69][70][71], yeasts [72][73][74], and fungal species [75][76][77][78][79][80][81]. Most interestingly, Dalmasso et al. [82] reported that P. freudenreichii CIRM-BIA1T adapted from warm (28 °C) to cold storage (4 °C) was able to express genes involved in the formation of important cheese flavor compounds through the catabolism of branched-chain amino acids. Similarly, the transcriptome profile of P. freudenreichii CIRM-BIA138 was analyzed during the adaptation of the strain to starvation [83], while that of P. freudenreichii ITG P14 during different stresses, i.e., cold, heat, and starvation [84]. Moreover, the transcriptome profile of L. rhamnosus PR1019 was recently evaluated in a cheese-like medium during carbon source starvation [85]. The analysis revealed that the strain was able to adapt under these conditions using alternative metabolic pathways, such as pyruvate degradation and ribose catabolism.

1.4. Metatranscriptomics

In contrast to the plethora of transcriptomics studies in single microbes, a limited number of metatranscriptomics projects have been performed to assess the entire gene expression of a food microbial community, such as that of cheese. Among the key problems are the short half-life of mRNA, the large amount of mRNA needed for analysis, the experimental design, and, most importantly, the selection of the proper sampling points and some technology-specific limitations, including the available bioinformatics tools and workflows for data processing and analysis [86]. The first comprehensive metatranscriptomics study was performed by Lessard et al. [87]. The authors monitored the metatranscriptome profiles of G. candidum and P. camemberti in the rind of an industrial Canadian Camembert-type cheese. Based on the functional annotation performed, transcripts related to energy metabolism, such as glycolysis/gluconeogenesis, the pentose phosphate pathway, the tricarboxylic acid cycle, and oxidative phosphorylation, were identified. Furthermore, lyases involved in the production of volatile sulfur compounds during methionine catabolism as well as transcripts related to cabbage sulfur aroma development and ammonia production were also detected [87].
More recently, Monnet et al. [88] evaluated the functional activity of the microbiome in the rind of a French Reblochon-type cheese, which was produced by S. thermophilus, L. bulgaricus, G. candidum, and D. hansenii. Metatranscriptomics analysis revealed that during ripening, only minor changes occurred in the LAB, such as an upregulation of genes involved in lipid and carbohydrate metabolism, while in yeasts, a significant upregulation of genes related to amino acid catabolism occurred from day 14 to day 35 (end of ripening), suggesting their contribution to flavor formation [88].
However, discrepancies observed between mRNA levels and protein abundance, which have been attributed to the variable half-lives of mRNA, variable rates of protein synthesis, and possible post-translational modifications of proteins [89], make the analysis of expressed proteins a valuable tool for investigating the functional activity of microbiota.

1.5. Metaproteomics

Metaproteomics is the large-scale characterization of the entire protein complement of microbial communities in a biological system of increased complexity at a given time point [90]. The progress of proteomics has been driven by the development of new technologies for peptide/protein separation, isotope labeling for quantification, and bioinformatics data analysis, while the analysis and quantification of proteins has been revolutionized by MS-based methods [91]. There are several challenges that must be overcome to address a metaproteomic study in animal food products. Among them are the presence of raw material proteins that interfere with the detection of microbial proteins, the complexity of the microbiota-expressed proteins in a fermented food present at concentration ranges that may vary dramatically, and finally, the presence of certain highly abundant proteins that are often not interesting for metaproteomic analysis, while the microbiota proteins may be several orders of magnitude less abundant. For these reasons, large-scale metaproteomic studies are limited. However, they are gradually gaining attention in the field of fermented food research, including cheese [92][93], to assess the functional diversity of microbial communities. To our knowledge, there are no metaproteomics studies on cheese flavor formation and development. However, the proteomic approach applied in a survey of bacterial proteins released in Emmental cheese revealed functional groups of proteins involved in proteolysis and glycolysis, pathways that influence the organoleptic characteristics of cheese [94].

1.6. Metabolomics

1.6.1. Analytical Techniques

Metabolomics is a relatively new field of omics research dealing with the simultaneous and high-throughput identification and quantification of low molecular mass (<1500 Da) metabolites that are not genetically encoded and are produced and modified by the metabolism of living organisms, such as organic acids, carbohydrates, amino acids, peptides, nucleic acids, vitamins, polyphenols, alkaloids, and minerals [95]. Metabolomic analyses have been classified as targeted (specific) or untargeted (non-selective) analyses.
The vast chemical diversity of the compounds in the metabolome requires efficient metabolite extraction, chromatographic separation, mass spectral detection, identification, quantification, and multivariate data analysis. Various analytical techniques are used for separation, such as liquid chromatography (LC), including high-performance LC (HPLC) or ultra-performance LC (UPLC), gas chromatography (GC), and capillary electrophoresis (CE), coupled to mass spectrometry (MS). MS-based platforms predominate because of both their ability to identify a wide range of compounds and their high throughput capacity. Apart from MS, other commonly used detection techniques are nuclear magnetic resonance (NMR) and near infrared spectrometry (NIR). However, none of the individual analytical methods are capable of effectively analyzing all the metabolites of a sample.

1.6.2. Metabolomics to Assess Flavor Formation in Cheese

Metabolomic studies in cheese provide insights on specific metabolites produced by cheese microbiota, with special interest in flavor compounds, and have significantly enabled characterization of the metabolite diversity in cheese, as well as the factors affecting it. According to Ochi et al. [96], GC/TOF-MS fingerprinting of hydrophilic low-molecular-mass compounds can be the basis of prediction models for sensory attributes, such as “rich flavor” and “sour flavor” in ripened Cheddar and Gouda, while specific amino acids [97] or volatile compounds, such as hexanoic acid, heptanoic acid, octanoic acid, 2-decenal, and acetoin [98], may be used as ripening markers in Cheddar cheese.
The cheese metabolome is affected by the degree of LAB autolysis. Grana-Padano cheeses produced in two different dairies and with different volatile profiles, as shown by solid phase micro-extraction (SPME) GC-MS, were analyzed with respect to their total culturable lactic microbiota and starter lysis. Remarkably, it was shown that increased complexity of microbial origin volatiles, such as ketones, alcohols, hydrocarbons, acetic acid, and propionic acid, was associated with the complex microbiota composition, with NSLAB (mainly L. rhamnosus/L. casei) being dominant. On the other hand, intracellular enzymes released due to SLAB cell lysis were involved in a higher content of FFAs, benzaldehyde, and organic acids, such as pyroglutamic and citric acid [99].
LC high-resolution MS (LC–HRMS) metabolic fingerprinting in model cheese revealed that the cheese metabolome is influenced by the spatial distribution of the starter L. lactis colonies throughout ripening [100]. By varying the time of renneting (at 0 h, e.g., simultaneously, or after 8 h from starter inoculation), they generated two different spatial distributions of immobilized bacterial colonies (few big colonies spread away from each other, or numerous small colonies close to each other, respectively), and identified 26 metabolites, including amino acids, organic acids, vitamins, nucleotides, and proteolysis products, being more abundant in small-colony cheeses. Moreover, according to Le Boucher et al., bacterial cells forming small colonies use the same metabolic pathways but display higher metabolic activity than big colonies, resulting in higher concentrations of metabolites being accumulated due to proteolysis or carbohydrate catabolism [101].
Investigation of the flavor compounds produced by the whole cheese microbiota or the individual contribution of SLAB or NSLAB strains is of great interest for cheese producers. Studies of this type are often performed in a cheese-based medium, mimicking the cheese environment, thus providing the nutrients and precursors present in cheese during ripening.
Sgarbi et al. explored the volatile compounds produced by different L. casei and L. rhamnosus strains grown on either a cheese-based medium (CBM) or on a LAB cell lysate-based medium (LCM). Volatile analysis by SPME GC-MS showed differences between NSLAB strains grown in CBM (where pyruvate could be a common precursor for all compounds produced) and LCM (where in the absence of lactate and citrate, NSLAB strains used mainly FFAs and FAAs), thus providing a better understanding of how NSLAB-produced volatile flavor compounds contribute to the development of cheese flavor during ripening, a first step toward the selection of wild NSLAB possessing a specific aromatic profile, for use as adjunct culture [102].
Similarly, Pogačić et al. evaluated the potential of Lactobacillus spp. and Leuconostoc spp. in the production of aroma compounds by incubating single strains in a curd-based slurry medium [103]. Analysis of volatiles by GC–MS revealed strain to strain variation. Acetoin, diacetyl, acids, and esters were mainly produced by L. rhamnosus and L. paracasei, while Leuconostoc spp. were major producers of alcohols and esters. The same curd-based slurry medium has been used by Pogačić et al. to screen LAB, Actinobacteria, P. freudenreichii, and Hafnia alvei strains for their ability to produce aroma compounds. Forty-nine out of 52 aroma compounds identified differed in their abundance among the bacteria [104].
A cheese-based medium was also used by Guarrasi et al. to grow SLAB and NSLAB strains isolated from Caciocavallo Palermitano cheese to assess their contribution to the volatile organic compounds (VOCs) production in ripened cheese [105]. By applying head space (HS)-SPME GC-MS to the fermented substrates, they found that strains of L. delbrueckii, L. casei, L. paracasei, L. rhamnosus, and Enterococcus gallinarum mainly influenced the development of the characteristic ripened-cheese aromatic compounds.
Furthermore, according to Yee et al., strains of dairy propionibacteria exhibited large inter-species and intra-species diversity in their ability to produce different aroma compounds when grown in a cheese curd-based medium [106]. GC-MS data differentiated P. freudenreichii strains from each other, as well as from Propionibacterium acidipropionici strains. For the same compound, differences between strains of the same species were as high as ~500-fold, with P. freudenreichii strains harboring the widest potential to produce cheese aroma compounds, making it a potential candidate to modulate cheese flavor.
Except cheese-based medium, cheese model systems are important tools that facilitate the study of the impact of microbial metabolism on cheese sensorial characteristics formation since they can be prepared under controlled microbiological conditions, are more economical, reproducible, and easier to obtain.
Miniature cheeses were prepared by Ruggirello et al. [107] to evaluate L. lactis commercial starter cultures’ viability and ability to produce aroma compounds during ripening. Culture-independent and -dependent approaches, as well as HS-SPME GC-MS, revealed the persistence of L. lactis in cheese throughout cheese making and ripening, although in a metabolically active hypothetical viable but not culturable (VBNC) state, since the expression of L. lactis cystathionine β-lyase (metC) and α-acetolactate synthase (als) genes, involved in the biosynthesis of sulfur aroma compounds and diacetyl/acetoin, respectively, was partially associated with acetoin, diacetyl, 2,3-butanediol, and dimethyl disulfide production in ripened cheeses.
Irlinger et al. [108] applied GC-MS for metabolic profiling of the Gram-negative bacteria Psychrobacter celer and H. alvei when grown on a mini smear cheese. Both bacteria were able to colonize the cheese surface and compete in the microbial community, modifying at the same time the aromatic content of cheeses, highlighting the fact that less abundant microorganisms can have a significant impact on cheese flavor.
In a recent study, Suzuki-Iwashima et al. [109] investigated the combined effects of LAB starters and the white fungus P. camemberti on the production of volatile compounds during the ripening of a model white mold surface-ripened cheese. Metabolomics analysis by GC-MS showed that the early ripening period was characterized by metabolites, such as lactose, galactose, lactic acid, ethanol, diacetyl, acetoin, ethyl acetate, and sulfur compounds (dimethyl sulfide and dimethyl disulfide), that derived from carbohydrate metabolism of LAB, while fungal metabolism of proteins (with amino acid-derived compounds including branched aldehydes, such as 3-methyl butanal), and of fatty acids and proteins (with methyl ketones, fatty acids and amino acids) characterized the intermediate- and the late-ripening stages, respectively.

1.6.3. Metabolomics Shed Light on the Ability of Adjunct Cultures to Produce Flavor Compounds

Starter bacteria combined with adjunct LAB strains selected for desirable metabolic potential is a common tool used to control and accelerate cheese ripening [110], enhance and improve flavor intensity [111], especially in cheeses made with pasteurized milk [112], and accelerate flavor development [113]. Metabolomics approaches have been a valuable tool to assess the contribution of the potential adjuncts on cheese flavor development [89][114][115].
So far, several studies have been performed evaluating not only the contribution of potential adjunct cultures to cheese flavor development, but also their persistence, technological performance, together with interactions with the starter cultures in cheese making experiments. In a recent study, Stefanovic et al. [21], selected three L. paracasei strains based on their proteolytic enzyme activities and ability to produce flavor compounds in cheese model systems [116] to be used as adjunct cultures in Cheddar cheese manufacture. The authors found that the L. paracasei strains contributed to the development and diversification of flavor-related compounds in short-aged cheeses. Different adjunct cultures did not influence the gross cheese composition, nor primary or secondary proteolysis or lipolysis. However, cheese volatile analysis by GC-MS showed variation in long-chain aldehydes, acids, and esters that originated from the metabolism of FFAs, suggesting that starter lipolytic activity produced the primary metabolites, which were further metabolized by the adjuncts into flavor-contributing compounds.
Bancalari et al. [117] evaluated the use of a wild L. paracasei strain, selected for its ability to produce in vitro acetoin and diacetyl, as adjunct culture to enhance the flavor formation in Caciotta-type cheese. Indeed, the adjunct strain was able to develop in curd and cheese, producing higher amounts of volatile compounds and organic acids as monitored through SPME GC-MS, thus differentiating the experimental Caciotta with respect to the control cheese. Moreover, Belkheir et al. [118] prepared a Tetilla-type cheese using as adjunct cultures a high-diacetyl producer L. plantarum strain together with a peptidolytic L. brevis strain producing volatile sulfur compounds, to study the volatile profile and sensory characteristics of the model cheese and to compare it with the PDO Tetilla cheese. Volatile analysis with SPME GC-MS showed an increased abundance of acetic acid, hexanoic acid, ethyl butanoate, and ethyl hexanoate, as well as higher scores for flavor preference in the cheese made with the two adjuncts than in the control cheese, highlighting the fact that the use of selected adjunct strains would differentiate the cheese sensory properties. Interestingly, when Kocuria varians (Micrococcaceae) and Y. lipolytica, selected for their proteolytic and lipolytic activities, were used as adjunct cultures, for the manufacture of experimental Tetilla cheese from pasteurized cow’s milk, Centeno et al. [119] were able to recover the traditional flavor and sensory characteristics of raw-milk PDO Tetilla cheese. The volatile profiles of cheese manufactured with both adjuncts, detected by GC-MS, showed enhanced formation of fatty acids, esters, and sulfur compounds, thus the modifying the flavor profile of the experimental Tetilla cheese, which was considered very similar to good-quality artisanal raw-milk cheese.
In a potentially probiotic Caciotta cheese, industrially produced with autochthonous putative probiotic Lactobacillus and Kluyveromyces strains as adjunct cultures, Pisano et al. [120] investigated the adjunct’s influence on the cheese chemical and microbiological composition and sensory properties. Cheese metabolome characterization by means of 1H NMR spectroscopy (for amino acids, organic acids, and carbohydrates) together with HPLC-diode array detector/evaporative light scattering detector (DAD/ELSD) for the cholesterol, α-tocopherol, and fatty acid composition, highlighted significant variations in the cheese metabolome both in terms of the ripening time and strain combination, with Kluyveromyces and Lactobacillus strains surviving the manufacturing process and retaining their viability till the end of ripening, suggesting that Caciotta cheese can be used as a carrier for probiotic bacteria delivery.

1.6.4. Metabolomics as a Means of Cheese Authentication

Effective and reliable analytical methods are of paramount importance to securing the authenticity of PDO cheeses, aiming to protect both the product value and consumers. To this end, metabolomics-based approaches represent a powerful method to discriminate fraudulent varieties of a given food product [121][122][123].
Pisano et al. [124] analyzed the polar metabolite profiles by GC-MS, together with the predominant cultivable microbiota from buffalo and cow Mozzarella, in order to discriminate them. PDO buffalo Mozzarella exhibited a higher microbial diversity together with less psychrotrophic bacteria, while cow Mozzarella showed the highest counts of S. thermophilus, originating from the commercial starter culture. Furthermore, the polar metabolites reflected differences in the production protocols and microbiota complexity of theses cheeses, suggesting that the polar metabolite profile can be a promising tool to characterize and verify the authenticity of Italian buffalo Mozzarella. Moreover, Rocchetti et al. [125] characterized low-molecular-mass metabolites based on ultra-high-pressure liquid chromatography coupled with quadrupole time-of-flight MS (UHPLC/QTOF-MS), aiming to reveal differences between genuine PDO and non-PDO Grana Padano cheeses. Amino acids, oligopeptides, and fatty acids were the biomarkers with the highest discriminatory power.
Another authenticity problem involves adulteration related to non-declared processing methods, such as in the case of Fiore Sardo (FS) cheese, where the use of raw ovine milk is mandatory. Caboni et al. studied the polar low-molecular-mass metabolites, by GC-MS, aiming to discriminate FS cheese produced from raw or thermized milk. FAAs and saccharides were the metabolites that mostly changed, suggesting the polar low-molecular-mass metabolites as a potential biomarker for detecting milk thermization in ovine PDO cheeses [126].

1.7. Integration

Studies combining multiple-omics approaches, integrating genomics, transcriptomics, together with metabolomics, that is, a systems biology approach, provide a detailed picture of the cheese microbiota dynamics, as well as information on potential microbial interactions and their contribution with respect to cheese sensorial characteristics development. An overview of the most comprehensive studies combining multiple omics approaches to study flavor development in cheese can be found in Supplementary Table S1. In this section, integrating omics studies, which have been performed not only in mature cheese, but also during cheese ripening, are discussed, to gain deeper insights on microbial succession and metabolite production.

1.7.1. Amplicon Metagenomics-Shotgun Metagenomics

The first pioneering study incorporating different -omics techniques was performed by Wolfe et al. [127]. The authors analyzed the rind microbial communities of 137 cheese samples, including 61 natural cheeses that were left undisturbed during the aging process, 52 washed with 20% w/w NaCl, and 24 bloomy cheese rinds, from different geographic regions, milk types, and milk treatments. Using amplicon-based metagenomics analysis, only 14 bacterial and 10 fungal genera were identified at higher than 1% average abundance. Interestingly, the microbiota was found to be associated with the rind type (natural, washed, and bloomy) and moisture instead of the geographic origin. Additionally, shotgun metagenomics was performed to assess the functional potential of the microbiota based on the different rind types. In particular, shotgun data identified genes involved in several metabolic pathways associated with flavor formation, such as the catabolism of cysteine and methionine, which are known to contribute to the production of volatile sulfur compounds, and that of valine, leucine, and isoleucine, which provide sweaty and putrid aromas. It should be noted that the majority of these genes were detected in the washed rind cheeses. Furthermore, the halotolerant γ-proteobacteria genus Pseudoalteromonas, originally associated with marine environments, was found for the first time in cheese microbiota and, in particular, in the natural and bloomy cheese rinds. The shotgun metagenomics results identified a few cold-adapted enzymes produced by Pseudoalteromonas spp. that participate in lipolysis and proteolysis. Therefore, the presence of Pseudoalteromonas, as part of the cheese microbial community, could be considered beneficial, as it can contribute to the development of flavor compounds in cheeses during ripening and storage at low temperatures [127].

1.7.2. Amplicon Metagenomics-Metabolomics-Metatranscriptomics

De Pasquale et al. [128] used Fiore Sardo, Pecorino Siciliano, and Pecorino Toscano cheeses as hard cheese model systems to study the spatial distribution of metabolically active microbiota and its effect on secondary proteolysis and VOC production, since these cheeses present a decreasing NaCl gradient from the surface to the center and an opposite moisture trend, properties that affect the microbiota distribution. By combining 16S rRNA gene pyrosequencing (targeting RNA) and Purge and Trap coupled with GC-MS (PT GC-MS), they found that in all cheese varieties, the poorest VOC profile was detected in the core region, due to the low oxygen availability, with high levels of alcohols originating from aldehydes and methyl-ketones reduction. Mesophilic lactobacilli (predominantly L. plantarum) positively correlated to alcohols, aldehydes, methyl and branched esters, and sulfur compounds. Thermophilic LAB in Pecorino Siciliano (including L. delbrueckii and S. thermophilus) positively correlated with the total concentration of FAAs, in different cheese regions, as well as with alcohols and related esters, while Brevibacterium sp. present on the surface of Pecorino Toscano correlated with alcohols, aldehydes, ketones, esters, and sulfur compounds.
Recently, Turri et al. [129] characterized the mature Historic Rebel (HR) cheese, an Italian heritage cheese produced from raw cow milk in the Alps. Microbiota diversity, assessed by 16S rRNA gene amplicon sequencing, revealed that the core microbiota comprising Streptococcus, Lactobacillus, Lactococcus, Leuconostoc, and Pediococcus genera correlated positively with the VOCs hexanal, 2-heptanal, 3-hydroxybutan-2-one (or acetoin), and ethanol respectively, as determined by SPME GC-MS. Moreover, a lipidomics approach was included by applying Dynamic Headspace (DHS) GC-MS to analyze the terpene fraction and the polyunsaturated fatty acids composition of HR cheese, parameters that were closely related with pasture vegetation and feeding, respectively, and contribute to the richness of cheese flavor.
In a sophisticated study [73], metabolite analysis by Ultra HPLC-MS (UHPLC-MS) and HPLC-UV in combination with dual RNA-seq analysis were applied in ripened lab-scale cheeses to provide insight into the metabolic interactions in a simple synthetic community composed of three species commonly used for the production of smear-ripened cheese. A strong mutualistic interaction between Brevibacterium aurantiacum and H. alvei was proposed, according to which, proteases and lipases secreted by B. aurantiacum liberate energy compounds from caseins and triglycerides that stimulate H. alvei, which in turn produces siderophore that increases iron availability for B. aurantiacum. Furthermore, the proteolytic activity of B. aurantiacum led to increased methionine catabolism in H. alvei producing methanethiol, a precursor for a wide variety of volatile sulfur compounds that contribute to cheese flavor.

1.7.3. Shotgun Metagenomics-Metatranscriptomics

Recently, Duru et al. [130] and DeFilippis et al. [131] studied the effect of modifying the cheese ripening temperature on microbial community structure and function. More specifically, in a thorough study, the metagenome and metatranscriptome profiles of the semi-hard Swiss-type Maasdam cheese were studied during warm (20 °C) and cold (5 °C) room ripening, using L. lactis subsp. lactis and L. lactis subsp. cremoris as starter strains and as adjunct cultures strains of P. freudenreichii subsp. shermanii, L. rhamnosus, and L. helveticus [130]. The authors constructed four genomes (one genome per species) from the shotgun data to near completeness (higher than 97%), and based on the mean DNA read coverage, L. lactis was found to be the dominant species. Annotation of genome assemblies and pathways reconstruction identified genes required for FFAs biosynthesis in all genomes, proteolytic enzymes only in the LAB, and lipolytic enzymes in the genomes of L. lactis, L. rhamnosus, and P. freudenreichii. Furthermore, genes for valine catabolism were found in L. lactis and P. freudenreichii, while those for methionine and cysteine catabolism in L. lactis, L. rhamnosus, and L. helveticus genomes. Therefore, all species are important for the flavor formation of Maasdam cheese to a different extent. RNA-seq analysis confirmed the dominance of L. lactis, as more than 85% of the transcript reads mapped uniquely to this species. Moreover, metatranscriptomic data showed that L. lactis was metabolically active despite the ripening temperature. However, this was not the case for the other species, as genes related to the central metabolism were downregulated during cold room ripening, suggesting that fewer flavor compounds were produced. On the other hand, according to De Filippis et al. [131], elevation of the ripening temperature of Caciocavallo Silano cheese from standard (16 °C) to experimental (20 °C) temperatures directly affected microbiota diversity and metabolism. More specifically, 16S rRNA amplicon and shotgun metatranscriptome sequencing revealed an increased relative abundance of NSLAB in cheese ripened at higher temperatures, together with differential expression of 651 genes. Furthermore, overexpression of proteolysis, lipolysis, amino, and fatty acid catabolism-related genes at 20 °C correlated with increased production of cheese VOCs, as found by SPME GC-MS, and significantly increased the cheese maturation rate.
The functional potential of the microbiome in an experimental surface-ripened cheese was recently assessed by Dugat-Bony et al. [132]. Using reference genomes, the assembly of the shotgun metagenomics data resulted in the genome construction of the nine microbes used for the production of the cheese, i.e., six bacteria species (L. lactis, Staphylococcus equorum, Corynebacterium casei, H. alvei, B. aurantiacum, and Arthrobacter arilaitensis), and three yeast species (D. hansenii, G. candidum, and K. lactis). In addition to the shotgun metagenomics, the authors also employed metatranscriptomics to link metabolically related transcript reads to the microbiome at specific time points during cheese ripening. RNA-seq analysis revealed that enzymes involved in lactose fermentation were found to be expressed by L. lactis and K. lactis. Moreover, lactate degradation was attributed to D. hansenii and G. candidum, due to the high levels of lactate dehydrogenase transcripts detected in these species. Furthermore, although several transcript reads associated with protein and lipid metabolism were mapped to L. lactis, it was found that G. candidum was the main contributor to proteolysis and lipolysis. In addition, G. candidum was also found to be a key microbe regarding the catabolism of amino acids and thus, the production of flavor compounds throughout ripening. However, it should be noted that transcripts related to amino acid catabolism also mapped to the L. lactis genome, mainly at the early stage of ripening, as well as to the C. casei and H. alvei genomes at the end of ripening [132].

References

  1. Bolotin, A.; Wincker, P.; Mauger, S.; Jaillon, O.; Malarme, K.; Weissenbach, J.; Ehrlich, S.D.; Sorokin, A. The complete genome sequence of the lactic acid bacterium Lactococcus lactis ssp. lactis IL1403. Genome Res. 2001, 11, 731–753.
  2. Kok, J.; van Gijtenbeek, L.A.; de Jong, A.; van der Meulen, S.B.; Solopova, A.; Kuipers, O.P. The evolution of gene regulation research in Lactococcus lactis. FEMS Microbiol. Rev. 2017, 41, S220–S243.
  3. Kelleher, P.; Bottacini, F.; Mahony, J.; Kilcawley, K.N.; van Sinderen, D. Comparative and functional genomics of the Lactococcus lactis taxon; insights into evolution and niche adaptation. BMC Genom. 2017, 18, 267.
  4. Kok, J.; Buist, G.; Zomer, A.L.; van Hijum, S.A.F.T.; Kuipers, O.P. Comparative and functional genomics of lactococci. FEMS Microbiol. Rev. 2005, 29, 411–433.
  5. Linares, D.M.; Kok, J.; Poolman, B. Genome sequences of Lactococcus lactis MG1363 (revised) and NZ9000 and comparative physiological studies. J. Bacteriol. 2010, 192, 5806–5812.
  6. Oliveira, A.P.; Nielsen, J.; Förster, J. Modeling Lactococcus lactis using a genome-scale flux model. BMC Microbiol. 2005, 5, 39.
  7. Wels, M.; Siezen, R.; van Hijum, S.; Kelly, W.J.; Bachmann, H. Comparative genome analysis of Lactococcus lactis indicates niche adaptation and resolves genotype/phenotype disparity. Front. Microbiol. 2019, 10, 4.
  8. Alexandraki, V.; Kazou, M.; Blom, J.; Pot, B.; Papadimitriou, K.; Tsakalidou, E. Comparative genomics of Streptococcus thermophilus support important traits concerning the evolution, biology and technological properties of the species. Front. Microbiol. 2019, 10, 2916.
  9. Broadbent, J.R.; Cai, H.; Larsen, R.L.; Hughes, J.E.; Welker, D.L.; De Carvalho, V.G.; Tompkins, T.A.; Ardo, Y.; Vogensen, F.; De Lorentiis, A.; et al. Genetic diversity in proteolytic enzymes and amino acid metabolism among Lactobacillus helveticus strains. J. Dairy Sci. 2011, 94, 4313–4328.
  10. Broadbent, J.R.; Hughes, J.E.; Welker, D.L.; Tompkins, T.A.; Steele, J.L. Complete genome sequence for Lactobacillus helveticus CNRZ 32, an industrial cheese starter and cheese flavor adjunct. Genome Announc. 2013, 1, e00590-13.
  11. Makarova, K.; Slesarev, A.; Wolf, Y.; Sorokin, A.; Mirkin, B.; Koonin, E.; Pavlov, A.; Pavlova, N.; Karamychev, V.; Polouchine, N.; et al. Comparative genomics of the lactic acid bacteria. Proc. Natl. Acad. Sci. USA 2006, 103, 15611–15616.
  12. Iskandar, C.F.; Cailliez-Grimal, C.; Borges, F.; Revol-Junelles, A.-M. Review of lactose and galactose metabolism in lactic acid bacteria dedicated to expert genomic annotation. Trends Food Sci. Technol. 2019, 88, 121–132.
  13. Buron-Moles, G.; Chailyan, A.; Dolejs, I.; Forster, J.; Mikš, M.H. Uncovering carbohydrate metabolism through a genotype-phenotype association study of 56 lactic acid bacteria genomes. Appl. Microbiol. Biotechnol. 2019, 103, 3135–3152.
  14. Cai, H.; Thompson, R.; Budinich, M.F.; Broadbent, J.R.; Steele, J.L. Genome sequence and comparative genome analysis of Lactobacillus casei: Insights into their niche-associated evolution. Genome Biol. Evol. 2009, 1, 239–257.
  15. Claesson, M.J.; van Sinderen, D.; O’Toole, P.W. The genus Lactobacillus—A genomic basis for understanding its diversity. FEMS Microbiol. Lett. 2007, 269, 22–28.
  16. Jiang, Y.; Zhang, J.; Zhao, X.; Zhao, W.; Yu, Z.; Chen, C.; Yang, Z. Complete genome sequencing of exopolysaccharide-producing Lactobacillus plantarum K25 provides genetic evidence for the probiotic functionality and cold endurance capacity of the strain. Biosci. Biotechnol. Biochem. 2018, 82, 1225–1233.
  17. Kazou, M.; Alexandraki, V.; Blom, J.; Pot, B.; Tsakalidou, E.; Papadimitriou, K. Comparative genomics of Lactobacillus acidipiscis ACA-DC 1533 isolated from traditional Greek Kopanisti cheese against species within the Lactobacillus salivarius clade. Front. Microbiol. 2018, 9, 1244.
  18. Liu, J.; Chan, S.H.J.; Chen, J.; Solem, C.; Jensen, P.R. Systems biology—A guide for understanding and developing improved strains of lactic acid bacteria. Front. Microbiol. 2019, 10, 876.
  19. Liu, M.; Nauta, A.; Francke, C.; Siezen Roland, J. Comparative genomics of enzymes in flavor-forming pathways from amino acids in lactic acid bacteria. Appl. Environ. Microbiol. 2008, 74, 4590–4600.
  20. Stefanovic, E.; Casey, A.; Cotter, P.; Cavanagh, D.; Fitzgerald, G.; McAuliffe, O. Draft genome sequence of Lactobacillus casei DPC6800, an isolate with the potential to diversify flavor in cheese. Genome Announc. 2016, 4, e00063-16.
  21. Stefanovic, E.; Kilcawley, K.N.; Roces, C.; Rea, M.C.; O’Sullivan, M.; Sheehan, J.J.; McAuliffe, O. Evaluation of the potential of Lactobacillus paracasei adjuncts for flavor compounds development and diversification in short-aged Cheddar cheese. Front. Microbiol. 2018, 9, 1506.
  22. Wu, C.; Huang, J.; Zhou, R. Genomics of lactic acid bacteria: Current status and potential applications. Crit. Rev. Microbiol. 2017, 43, 393–404.
  23. Zheng, J.; Ruan, L.; Sun, M.; Gänzle, M. A genomic view of lactobacilli and pediococci demonstrates that phylogeny matches ecology and physiology. Appl. Environ. Microbiol. 2015, 81, 7233–7243.
  24. Liu, M.; Bayjanov, J.R.; Renckens, B.; Nauta, A.; Siezen, R.J. The proteolytic system of lactic acid bacteria revisited: A genomic comparison. BMC Genom. 2010, 11, 36.
  25. Sun, Z.; Harris, H.M.B.; McCann, A.; Guo, C.; Argimón, S.; Zhang, W.; Yang, X.; Jeffery, I.B.; Cooney, J.C.; Kagawa, T.F.; et al. Expanding the biotechnology potential of lactobacilli through comparative genomics of 213 strains and associated genera. Nat. Commun. 2015, 6, 8322.
  26. Flahaut, N.A.; Wiersma, A.; van de Bunt, B.; Martens, D.E.; Schaap, P.J.; Sijtsma, L.; Dos Santos, V.A.; de Vos, W.M. Genome-scale metabolic model for Lactococcus lactis MG1363 and its application to the analysis of flavor formation. Appl. Microbiol. Biotechnol. 2013, 97, 8729–8739.
  27. Koduru, L.; Kim, Y.; Bang, J.; Lakshmanan, M.; Han, N.S.; Lee, D.-Y. Genome-scale modeling and transcriptome analysis of Leuconostoc mesenteroides unravel the redox governed metabolic states in obligate heterofermentative lactic acid bacteria. Sci. Rep. 2017, 7, 15721.
  28. Kristjansdottir, T.; Bosma, E.F.; Branco dos Santos, F.; Özdemir, E.; Herrgård, M.J.; França, L.; Ferreira, B.; Nielsen, A.T.; Gudmundsson, S. A metabolic reconstruction of Lactobacillus reuteri JCM 1112 and analysis of its potential as a cell factory. Microb. Cell Fact. 2019, 18, 186.
  29. Pastink, M.I.; Teusink, B.; Hols, P.; Visser, S.; de Vos, W.M.; Hugenholtz, J. Genome-scale model of Streptococcus thermophilus LMG18311 for metabolic comparison of lactic acid bacteria. Appl. Environ. Microbiol. 2009, 75, 3627–3633.
  30. Vinay-Lara, E.; Hamilton, J.J.; Stahl, B.; Broadbent, J.R.; Reed, J.L.; Steele, J.L. Genome-scale reconstruction of metabolic networks of Lactobacillus casei ATCC 334 and 12A. PLoS ONE 2014, 9, e110785.
  31. Xu, N.; Liu, J.; Ai, L.; Liu, L. Reconstruction and analysis of the genome-scale metabolic model of Lactobacillus casei LC2W. Gene 2015, 554, 140–147.
  32. Falentin, H.; Deutsch, S.-M.; Jan, G.; Loux, V.; Thierry, A.; Parayre, S.; Maillard, M.-B.; Dherbécourt, J.; Cousin, F.J.; Jardin, J.; et al. The complete genome of Propionibacterium freudenreichii CIRM-BIA1T, a hardy actinobacterium with food and probiotic applications. PLoS ONE 2010, 5, e11748.
  33. Koskinen, P.; Deptula, P.; Smolander, O.-P.; Tamene, F.; Kammonen, J.; Savijoki, K.; Paulin, L.; Piironen, V.; Auvinen, P.; Varmanen, P. Complete genome sequence of Propionibacterium freudenreichii DSM 20271(T). Stand. Genom. Sci. 2015, 10, 83.
  34. Piwowarek, K.; Lipińska, E.; Hać-Szymańczuk, E.; Kieliszek, M.; Kot, A.M. Sequencing and analysis of the genome of Propionibacterium Freudenreichii T82 strain: Importance for industry. Biomolecules 2020, 10, 348.
  35. Anast, J.M.; Dzieciol, M.; Schultz, D.L.; Wagner, M.; Mann, E.; Schmitz-Esser, S. Brevibacterium from Austrian hard cheese harbor a putative histamine catabolism pathway and a plasmid for adaptation to the cheese environment. Sci. Rep. 2019, 9, 6164.
  36. Levesque, S.; de Melo, A.G.; Labrie, S.J.; Moineau, S. Mobilome of Brevibacterium aurantiacum sheds light on its genetic diversity and its adaptation to smear-ripened cheeses. Front. Microbiol. 2019, 10, 1270.
  37. Olender, A.; Rutyna, P.; Niemcewicz, M.; Bogut, A.; Ciesielka, M.; Teresiński, G. Draft whole-genome sequence of Brevibacterium casei strain isolated from a bloodstream infection. Braz. J. Microbiol. 2020, 51, 685–689.
  38. Pham, N.-P.; Layec, S.; Dugat-Bony, E.; Vidal, M.; Irlinger, F.; Monnet, C. Comparative genomic analysis of Brevibacterium strains: Insights into key genetic determinants involved in adaptation to the cheese habitat. BMC Genom. 2017, 18, 955.
  39. Cheeseman, K.; Ropars, J.; Renault, P.; Dupont, J.; Gouzy, J.; Branca, A.; Abraham, A.-L.; Ceppi, M.; Conseiller, E.; Debuchy, R.; et al. Multiple recent horizontal transfers of a large genomic region in cheese making fungi. Nat. Commun. 2014, 5, 2876.
  40. Dumas, E.; Feurtey, A.; Rodríguez de la Vega, R.C.; Le Prieur, S.; Snirc, A.; Coton, M.; Thierry, A.; Coton, E.; Le Piver, M.; Roueyre, D.; et al. Independent domestication events in the blue-cheese fungus Penicillium roqueforti. Mol. Ecol. 2020, 29, 2639–2660.
  41. Gillot, G.; Jany, J.-L.; Coton, M.; Le Floch, G.; Debaets, S.; Ropars, J.; López-Villavicencio, M.; Dupont, J.; Branca, A.; Giraud, T.; et al. Insights into Penicillium roqueforti morphological and genetic diversity. PLoS ONE 2015, 10, e0129849.
  42. Perkins, V.; Vignola, S.; Lessard, M.-H.; Plante, P.-L.; Corbeil, J.; Dugat-Bony, E.; Frenette, M.; Labrie, S. Phenotypic and genetic characterization of the cheese ripening yeast Geotrichum candidum. Front. Microbiol. 2020, 11, 737.
  43. Ropars, J.; Didiot, E.; Rodríguez de la Vega, R.C.; Bennetot, B.; Coton, M.; Poirier, E.; Coton, E.; Snirc, A.; Le Prieur, S.; Giraud, T. Domestication of the emblematic white cheese-making fungus Penicillium camemberti and its diversification into two varieties. Curr. Biol. 2020, 30, 4441–4453.e4444.
  44. De Filippis, F.; Parente, E.; Ercolini, D. Recent past, present, and future of the food microbiome. Annu. Rev. Food Sci. Technol. 2018, 9, 589–608.
  45. Escobar-Zepeda, A.; Sanchez-Flores, A.; Quirasco Baruch, M. Metagenomic analysis of a Mexican ripened cheese reveals a unique complex microbiota. Food Microbiol. 2016, 57, 116–127.
  46. de Jong, A.; Hansen, M.E.; Kuipers, O.P.; Kilstrup, M.; Kok, J. The transcriptional and gene regulatory network of Lactococcus lactis MG1363 during growth in milk. PLoS ONE 2013, 8, e53085.
  47. Del Rio, B.; Redruello, B.; Martin, M.C.; Fernandez, M.; de Jong, A.; Kuipers, O.P.; Ladero, V.; Alvarez, M.A. Transcriptome profiling of Lactococcus lactis subsp. cremoris CECT 8666 in response to agmatine. Genom. Data 2016, 7, 112–114.
  48. Desfossés-Foucault, É.; LaPointe, G.; Roy, D. Transcription profiling of interactions between Lactococcus lactis subsp. cremoris SK11 and Lactobacillus paracasei ATCC 334 during Cheddar cheese simulation. Int. J. Food Microbiol. 2014, 178, 76–86.
  49. Dressaire, C.; Redon, E.; Gitton, C.; Loubiere, P.; Monnet, V.; Cocaign-Bousquet, M. Investigation of the adaptation of Lactococcus lactis to isoleucine starvation integrating dynamic transcriptome and proteome information. Microb. Cell Fact. 2011, 10 (Suppl. S1), S18.
  50. Gupta, J.A.; Thapa, S.; Verma, M.; Som, R.; Mukherjee, K.J. Genomics and transcriptomics analysis reveals the mechanism of isobutanol tolerance of a laboratory evolved Lactococcus lactis strain. Sci. Rep. 2020, 10, 10850.
  51. Kuipers, O.P.; de Jong, A.; Baerends, R.J.S.; van Hijum, S.A.F.T.; Zomer, A.L.; Karsens, H.A.; den Hengst, C.D.; Kramer, N.E.; Buist, G.; Kok, J. Transcriptome analysis and related databases of Lactococcus lactis. In Proceedings of the 7th Symposium on Lactic Acid Bacteria, Egmond aan Zee, The Netherlands, 1–5 September 2002; Siezen, R.J., Kok, J., Abee, T., Schasfsma, G., Eds.; Springer: Dordrecht, The Netherlands, 2002; pp. 113–122.
  52. Nouaille, S.; Even, S.; Charlier, C.; Le Loir, Y.; Cocaign-Bousquet, M.; Loubière, P. Transcriptomic response of Lactococcus lactis in mixed culture with Staphylococcus aureus. Appl. Environ. Microbiol. 2009, 75, 4473–4482.
  53. Puri, P.; Wetzel, C.; Saffert, P.; Gaston, K.W.; Russell, S.P.; Cordero Varela, J.A.; van der Vlies, P.; Zhang, G.; Limbach, P.A.; Ignatova, Z.; et al. Systematic identification of tRNAome and its dynamics in Lactococcus lactis. Mol. Microbiol. 2014, 93, 944–956.
  54. Wu, H.; Liu, J.; Miao, S.; Zhao, Y.; Zhu, H.; Qiao, M.; Saris, P.E.J.; Qiao, J. Contribution of YthA, a PspC family transcriptional regulator of Lactococcus lactis F44 acid tolerance and nisin yield: A transcriptomic approach. Appl. Environ. Microbiol. 2018, 84, e02483-17.
  55. Cretenet, M.; Laroute, V.; Ulvé, V.; Jeanson, S.; Nouaille, S.; Even, S.; Piot, M.; Girbal, L.; Le Loir, Y.; Loubière, P.; et al. Dynamic analysis of the Lactococcus lactis transcriptome in cheeses made from milk concentrated by ultrafiltration reveals multiple strategies of adaptation to stresses. Appl. Environ. Microbiol. 2011, 77, 247–257.
  56. Duwat, P.; Sourice, S.; Cesselin, B.; Lamberet, G.; Vido, K.; Gaudu, P.; Le Loir, Y.; Violet, F.; Loubière, P.; Gruss, A. Respiration capacity of the fermenting bacterium Lactococcus lactis and its positive effects on growth and survival. J. Bacteriol. 2001, 183, 4509–4516.
  57. Larsen, N.; Moslehi-Jenabian, S.; Werner, B.B.; Jensen, M.L.; Garrigues, C.; Vogensen, F.K.; Jespersen, L. Transcriptome analysis of Lactococcus lactis subsp. lactis during milk acidification as affected by dissolved oxygen and the redox potential. Int. J. Food Microbiol. 2016, 226, 5–12.
  58. Lu, J.; Cui, L.; Lin, S.; Hao, L.; Cao, N.; Yi, J.; Liu, X.; Lu, L.; Kang, Q. Short communication: Global transcriptome analysis of Lactococcus lactis ssp. lactis in response to gradient freezing. J. Dairy Sci. 2019, 102, 3933–3938.
  59. Taïbi, A.; Dabour, N.; Lamoureux, M.; Roy, D.; LaPointe, G. Comparative transcriptome analysis of Lactococcus lactis subsp. cremoris strains under conditions simulating Cheddar cheese manufacture. Int. J. Food Microbiol. 2011, 146, 263–275.
  60. van der Meulen, S.B.; de Jong, A.; Kok, J. Early transcriptome response of Lactococcus lactis to environmental stresses reveals differentially expressed small regulatory RNAs and tRNAs. Front. Microbiol. 2017, 8, 1704.
  61. Zhu, Z.; Yang, P.; Wu, Z.; Zhang, J.; Du, G. Systemic understanding of Lactococcus lactis response to acid stress using transcriptomics approaches. J. Ind. Microbiol. Biotechnol. 2019, 46, 1621–1629.
  62. Dressaire, C.; Redon, E.; Milhem, H.; Besse, P.; Loubière, P.; Cocaign-Bousquet, M. Growth rate regulated genes and their wide involvement in the Lactococcus lactis stress responses. BMC Genom. 2008, 9, 343.
  63. van der Meulen, S.B.; de Jong, A.; Kok, J. Transcriptome landscape of Lactococcus lactis reveals many novel RNAs including a small regulatory RNA involved in carbon uptake and metabolism. RNA Biol. 2016, 13, 353–366.
  64. Zomer, A.L.; Buist, G.; Larsen, R.; Kok, J.; Kuipers, O.P. Time-resolved determination of the CcpA regulon of Lactococcus lactis subsp. cremoris MG1363. J. Bacteriol. 2007, 189, 1366–1381.
  65. Chen, C.; Wang, L.; Lu, Y.; Yu, H.; Tian, H. Comparative transcriptional analysis of Lactobacillus plantarum and its ccpA-knockout mutant under galactooligosaccharides and glucose conditions. Front. Microbiol. 2019, 10, 1584.
  66. Chun, B.H.; Kim, K.H.; Jeon, H.H.; Lee, S.H.; Jeon, C.O. Pan-genomic and transcriptomic analyses of Leuconostoc mesenteroides provide insights into its genomic and metabolic features and roles in kimchi fermentation. Sci. Rep. 2017, 7, 11504.
  67. Giaretta, S.; Treu, L.; Vendramin, V.; da Silva Duarte, V.; Tarrah, A.; Campanaro, S.; Corich, V.; Giacomini, A. Comparative Transcriptomic analysis of Streptococcus thermophilus TH1436 and TH1477 showing different capability in the use of galactose. Front. Microbiol. 2018, 9, 1765.
  68. Jääskeläinen, E.; Vesterinen, S.; Parshintsev, J.; Johansson, P.; Riekkola, M.-L.; Björkroth, J. Production of buttery-odor compounds and transcriptome response in Leuconostoc gelidum subsp. gasicomitatum LMG18811T during growth on various carbon sources. Appl. Environ. Microbiol. 2015, 81, 1902–1908.
  69. Kim, S.Y.; Kim, S.A.; Jang, Y.J.; Seo, S.O.; Han, N.S. Screening of endogenous strong promoters of Leuconostoc citreum EFEL2700 based on transcriptome analysis and its application for food-grade production of β-galactosidase. J. Biotechnol. 2021, 325, 1–6.
  70. Li, J.-s.; Bi, Y.-t.; Dong, C.; Yang, J.-f.; Liang, W.-d. Transcriptome analysis of adaptive heat shock response of Streptococcus thermophilus. PLoS ONE 2011, 6, e25777.
  71. Padmanabhan, A.; Tong, Y.; Wu, Q.; Zhang, J.; Shah, N.P. Transcriptomic insights into the growth phase- and sugar-associated changes in the exopolysaccharide production of a high eps-producing Streptococcus thermophilus ASCC 1275. Front. Microbiol. 2018, 9, 1919.
  72. Malek, R.; Bonnarme, P.; Irlinger, F.; Frey-Klett, P.; Onésime, D.; Aubert, J.; Loux, V.; Beckerich, J.-M. Transcriptomic response of Debaryomyces hansenii during mixed culture in a liquid model cheese medium with Yarrowia lipolytica. Int. J. Food Microbiol. 2018, 264, 53–62.
  73. Pham, N.-P.; Landaud, S.; Lieben, P.; Bonnarme, P.; Monnet, C. Transcription profiling reveals cooperative metabolic interactions in a microbial cheese-ripening community composed of Debaryomyces hansenii, Brevibacterium aurantiacum, and Hafnia alvei. Front. Microbiol. 2019, 10, 1901.
  74. Prista, C.; Michán, C.; Miranda, I.M.; Ramos, J. The halotolerant Debaryomyces hansenii, the Cinderella of non-conventional yeasts. Yeast 2016, 33, 523–533.
  75. Becerra, M.; González-Siso, M.I.; Cerdán, M.E. A transcriptome analyis of Kluyveromyces lactis growing in cheese whey. Int. Dairy J. 2006, 16, 207–214.
  76. Castellote, J.; Fraud, S.; Irlinger, F.; Swennen, D.; Fer, F.; Bonnarme, P.; Monnet, C. Investigation of Geotrichum candidum gene expression during the ripening of Reblochon-type cheese by reverse transcription-quantitative PCR. Int. J. Food Microbiol. 2015, 194, 54–61.
  77. Pracharova, P.; Lieben, P.; Pollet, B.; Beckerich, J.M.; Bonnarme, P.; Landaud, S.; Swennen, D. Geotrichum candidum gene expression and metabolite accumulation inside the cells reflect the strain oxidative stress sensitivity and ability to produce flavour compounds. FEMS Yeast Res. 2019, 19, foy111.
  78. Punt, M.; van den Brule, T.; Teertstra, W.R.; Dijksterhuis, J.; den Besten, H.M.W.; Ohm, R.A.; Wösten, H.A.B. Impact of maturation and growth temperature on cell-size distribution, heat-resistance, compatible solute composition and transcription profiles of Penicillium roqueforti conidia. Food Res. Int. 2020, 136, 109287.
  79. Rojas-Aedo, J.F.; Gil-Durán, C.; Del-Cid, A.; Valdés, N.; Álamos, P.; Vaca, I.; García-Rico, R.O.; Levicán, G.; Tello, M.; Chávez, R. The biosynthetic gene cluster for andrastin A in Penicillium roqueforti. Front. Microbiol. 2017, 8, 813.
  80. Suleau, A.; Gourdon, P.; Reitz-Ausseur, J.; Casaregola, S. Transcriptomic analysis of extensive changes in metabolic regulation in Kluyveromyces lactis strains. Eukaryot. Cell 2006, 5, 1360–1370.
  81. Suleau, A.; Jacques, N.; Reitz-Ausseur, J.; Casaregola, S. Intraspecific gene expression variability in the yeast Kluyveromyces lactis revealed by micro-array analysis. FEMS Yeast Res. 2005, 5, 595–604.
  82. Dalmasso, M.; Aubert, J.; Briard-Bion, V.; Chuat, V.; Deutsch, S.-M.; Even, S.; Falentin, H.; Jan, G.; Jardin, J.; Maillard, M.-B.; et al. A temporal -omic study of Propionibacterium freudenreichii CIRM-BIA1T adaptation strategies in conditions mimicking cheese ripening in the cold. PLoS ONE 2012, 7, e29083.
  83. Aburjaile, F.F.; Rohmer, M.; Parrinello, H.; Maillard, M.B.; Beaucher, E.; Henry, G.; Nicolas, A.; Madec, M.N.; Thierry, A.; Parayre, S.; et al. Adaptation of Propionibacterium freudenreichii to long-term survival under gradual nutritional shortage. BMC Genom. 2016, 17, 1007.
  84. Falentin, H.; Postollec, F.; Parayre, S.; Henaff, N.; Le Bivic, P.; Richoux, R.; Thierry, A.; Sohier, D. Specific metabolic activity of ripening bacteria quantified by real-time reverse transcription PCR throughout Emmental cheese manufacture. Int. J. Food Microbiol. 2010, 144, 10–19.
  85. Lazzi, C.; Turroni, S.; Mancini, A.; Sgarbi, E.; Neviani, E.; Brigidi, P.; Gatti, M. Transcriptomic clues to understand the growth of Lactobacillus rhamnosus in cheese. BMC Microbiol. 2014, 14, 28.
  86. Shakya, M.; Lo, C.-C.; Chain, P.S.G. Advances and challenges in metatranscriptomic analysis. Front. Genet. 2019, 10, 904.
  87. Lessard, M.H.; Viel, C.; Boyle, B.; St-Gelais, D.; Labrie, S. Metatranscriptome analysis of fungal strains Penicillium camemberti and Geotrichum candidum reveal cheese matrix breakdown and potential development of sensory properties of ripened Camembert-type cheese. BMC Genom. 2014, 15, 235.
  88. Monnet, C.; Dugat-Bony, E.; Swennen, D.; Beckerich, J.-M.; Irlinger, F.; Fraud, S.; Bonnarme, P. Investigation of the activity of the microorganisms in a Reblochon-style cheese by metatranscriptomic analysis. Front. Microbiol. 2016, 7, 536.
  89. Afshari, R.; Pillidge, C.; Dias, D.; Osborn, A.; Gill, H. Cheesomics: The future pathway to understanding cheese flavour and quality. Crit. Rev. Food Sci. Nutr. 2018, 60, 33–47.
  90. Wilmes, P.; Bond, P.L. The application of two-dimensional polyacrylamide gel electrophoresis and downstream analyses to a mixed community of prokaryotic microorganisms. Environ. Microbiol. 2004, 6, 911–920.
  91. Zhang, Y.; Fonslow, B.R.; Shan, B.; Baek, M.-C.; Yates, J.R. Protein analysis by shotgun/bottom-up proteomics. Chem. Rev. 2013, 113, 2343–2394.
  92. Soggiu, A.; Piras, C.; Mortera, S.L.; Alloggio, I.; Urbani, A.; Bonizzi, L.; Roncada, P. Unravelling the effect of clostridia spores and lysozyme on microbiota dynamics in Grana Padano cheese: A metaproteomics approach. J. Proteom. 2016, 147, 21–27.
  93. Yang, L.; Fan, W.; Xu, Y. Metaproteomics insights into traditional fermented foods and beverages. Compr. Rev. Food Sci. Food Saf. 2020, 19, 2506–2529.
  94. Gagnaire, V.; Piot, M.; Camier, B.; Vissers, J.P.; Jan, G.; Léonil, J. Survey of bacterial proteins released in cheese: A proteomic approach. Int. J. Food Microbiol. 2004, 94, 185–201.
  95. Mozzi, F.; Ortiz, M.E.; Bleckwedel, J.; De Vuyst, L.; Pescuma, M. Metabolomics as a tool for the comprehensive understanding of fermented and functional foods with lactic acid bacteria. Food Res. Int. 2013, 54, 1152–1161.
  96. Ochi, H.; Bamba, T.; Naito, H.; Iwatsuki, K.; Fukusaki, E. Metabolic fingerprinting of hard and semi-hard natural cheeses using gas chromatography with flame ionization detector for practical sensory prediction modeling. J. Biosci. Bioeng. 2012, 114, 506–511.
  97. Ochi, H.; Sakai, Y.; Koishihara, H.; Abe, F.; Bamba, T.; Fukusaki, E. Monitoring the ripening process of Cheddar cheese based on hydrophilic component profiling using gas chromatography-mass spectrometry. J. Dairy Sci. 2013, 96, 7427–7441.
  98. Gan, H.H.; Yan, B.; Linforth, R.S.T.; Fisk, I.D. Development and validation of an APCI-MS/GC-MS approach for the classification and prediction of Cheddar cheese maturity. Food Chem. 2016, 190, 442–447.
  99. Lazzi, C.; Povolo, M.; Locci, F.; Bernini, V.; Neviani, E.; Gatti, M. Can the development and autolysis of lactic acid bacteria influence the cheese volatile fraction? The case of Grana Padano. Int. J. Food Microbiol. 2016, 233, 20–28.
  100. Le Boucher, C.; Courant, F.; Royer, A.L.; Jeanson, S.; Lortal, S.; Dervilly-Pinel, G.; Thierry, A.; Le Bizec, B. LC–HRMS fingerprinting as an efficient approach to highlight fine differences in cheese metabolome during ripening. Metabolomics 2015, 11, 1117–1130.
  101. Le Boucher, C.; Gagnaire, V.; Briard-Bion, V.; Jardin, J.; Maillard, M.B.; Dervilly-Pinel, G.; Le Bizec, B.; Lortal, S.; Jeanson, S.; Thierry, A. Spatial distribution of Lactococcus lactis colonies modulates the production of major metabolites during the ripening of a model cheese. Appl. Environ. Microbiol. 2016, 82, 202–210.
  102. Sgarbi, E.; Lazzi, C.; Tabanelli, G.; Gatti, M.; Neviani, E.; Gardini, F. Nonstarter lactic acid bacteria volatilomes produced using cheese components. J. Dairy Sci. 2013, 96, 4223–4234.
  103. Pogačić, T.; Maillard, M.B.; Leclerc, A.; Hervé, C.; Chuat, V.; Valence, F.; Thierry, A. Lactobacillus and Leuconostoc volatilomes in cheese conditions. Appl. Microbiol. Biotechnol. 2016, 100, 2335–2346.
  104. Pogačić, T.; Maillard, M.-B.; Leclerc, A.; Hervé, C.; Chuat, V.; Yee, A.L.; Valence, F.; Thierry, A. A methodological approach to screen diverse cheese-related bacteria for their ability to produce aroma compounds. Food Microbiol. 2015, 46, 145–153.
  105. Guarrasi, V.; Sannino, C.; Moschetti, M.; Bonanno, A.; Di Grigoli, A.; Settanni, L. The individual contribution of starter and non-starter lactic acid bacteria to the volatile organic compound composition of Caciocavallo Palermitano cheese. Int. J. Food Microbiol. 2017, 259, 35–42.
  106. Yee, A.L.; Maillard, M.B.; Roland, N.; Chuat, V.; Leclerc, A.; Pogačić, T.; Valence, F.; Thierry, A. Great interspecies and intraspecies diversity of dairy propionibacteria in the production of cheese aroma compounds. Int. J. Food Microbiol. 2014, 191, 60–68.
  107. Ruggirello, M.; Giordano, M.; Bertolino, M.; Ferrocino, I.; Cocolin, L.; Dolci, P. Study of Lactococcus lactis during advanced ripening stages of model cheeses characterized by GC-MS. Food Microbiol. 2018, 74, 132–142.
  108. Irlinger, F.; Yung, S.A.; Sarthou, A.S.; Delbès-Paus, C.; Montel, M.C.; Coton, E.; Coton, M.; Helinck, S. Ecological and aromatic impact of two Gram-negative bacteria (Psychrobacter celer and Hafnia alvei) inoculated as part of the whole microbial community of an experimental smear soft cheese. Int. J. Food Microbiol. 2012, 153, 332–338.
  109. Suzuki-Iwashima, A.; Matsuura, H.; Iwasawa, A.; Shiota, M. Metabolomics analyses of the combined effects of lactic acid bacteria and Penicillium camemberti on the generation of volatile compounds in model mold-surface-ripened cheeses. J. Biosci. Bioeng. 2020, 129, 333–347.
  110. Gobbetti, M.; Di Cagno, R.; Calasso, M.; Neviani, E.; Fox, P.F.; De Angelis, M. Drivers that establish and assembly the lactic acid bacteria biota in cheeses. Trends Food Sci. Technol. 2018, 78, 244–254.
  111. Settanni, L.; Moschetti, G. Non-starter lactic acid bacteria used to improve cheese quality and provide health benefits. Food Microbiol. 2010, 27, 691–697.
  112. Gobbetti, M.; De Angelis, M.; Di Cagno, R.; Mancini, L.; Fox, P.F. Pros and cons for using non-starter lactic acid bacteria (NSLAB) as secondary/adjunct starters for cheese ripening. Trends Food Sci. Technol. 2015, 45, 167–178.
  113. Smid, E.J.; Kleerebezem, M. Production of aroma compounds in lactic fermentations. Annu. Rev. Food Sci. Technol. 2014, 5, 313–326.
  114. Blaya, J.; Barzideh, Z.; LaPointe, G. Symposium review: Interaction of starter cultures and nonstarter lactic acid bacteria in the cheese environment1. J. Dairy Sci. 2018, 101, 3611–3629.
  115. Khattab, A.R.; Guirguis, H.A.; Tawfik, S.M.; Farag, M.A. Cheese ripening: A review on modern technologies towards flavor enhancement, process acceleration and improved quality assessment. Trends Food Sci. Technol. 2019, 88, 343–360.
  116. Stefanovic, E.; Thierry, A.; Maillard, M.B.; Bertuzzi, A.; Rea, M.C.; Fitzgerald, G.; McAuliffe, O.; Kilcawley, K.N. Strains of the Lactobacillus casei group show diverse abilities for the production of flavor compounds in 2 model systems. J. Dairy Sci. 2017, 100, 6918–6929.
  117. Bancalari, E.; Montanari, C.; Levante, A.; Alinovi, M.; Neviani, E.; Gardini, F.; Gatti, M. Lactobacillus paracasei 4341 as adjunct culture to enhance flavor in short ripened Caciotta-type cheese. Food Res. Int. 2020, 135, 109284.
  118. Belkheir, K.; Zadi-Karam, H.; Karam, N.-E.; Carballo, J.; Centeno, J.A. Effects of selected mesophilic Lactobacillus strains obtained from camel milk on the volatile and sensory profiles of a model short-ripened pressed cows’ milk cheese. Int. Dairy J. 2020, 109, 104738.
  119. Centeno, J.; Garabal, J.I.; Docampo, F.; Lorenzo, J.M.; Carballo, J. Recovering traditional raw-milk Tetilla cheese flavour and sensory attributes by using Kocuria varians and Yarrowia lipolytica adjunct cultures. Int. J. Food Microbiol. 2017, 251, 33–40.
  120. Pisano, M.B.; Rosa, A.; Putzu, D.; Cesare Marincola, F.; Mossa, V.; Viale, S.; Fadda, M.E.; Cosentino, S. Influence of autochthonous putative probiotic cultures on microbiota, lipid components and metabolome of Caciotta cheese. Front. Microbiol. 2020, 11, 583745.
  121. Kamilari, E.; Tomazou, M.; Antoniades, A.; Tsaltas, D. High throughput sequencing technologies as a new toolbox for deep analysis, characterization and potentially authentication of protection designation of origin cheeses? Int. J. Food Sci. 2019, 2019, 5837301.
  122. Checa, A.; Saurina, J. Metabolomics and PDO. In Comprehensive Analytical Chemistry; de la Guardia, M., Gonzálvez, A., Eds.; Elsevier: Amsterdam, The Netherlands, 2013; Volume 60, pp. 123–143.
  123. Danezis, G.P.; Pappas, A.C.; Tsiplakou, E.; Pappa, E.C.; Zacharioudaki, M.; Tsagkaris, A.S.; Papachristidis, C.A.; Sotirakoglou, K.; Zervas, G.; Georgiou, C.A. Authentication of Greek Protected Designation of Origin cheeses through elemental metabolomics. Int. Dairy J. 2020, 104, 104599.
  124. Pisano, M.B.; Scano, P.; Murgia, A.; Cosentino, S.; Caboni, P. Metabolomics and microbiological profile of Italian mozzarella cheese produced with buffalo and cow milk. Food Chem. 2016, 192, 618–624.
  125. Rocchetti, G.; Lucini, L.; Gallo, A.; Masoero, F.; Trevisan, M.; Giuberti, G. Untargeted metabolomics reveals differences in chemical fingerprints between PDO and non-PDO Grana Padano cheeses. Food Res. Int. 2018, 113, 407–413.
  126. Caboni, P.; Maxia, D.; Scano, P.; Addis, M.; Dedola, A.; Pes, M.; Murgia, A.; Casula, M.; Profumo, A.; Pirisi, A. A gas chromatography-mass spectrometry untargeted metabolomics approach to discriminate Fiore Sardo cheese produced from raw or thermized ovine milk. J. Dairy Sci. 2019, 102, 5005–5018.
  127. Wolfe, B.E.; Button, J.E.; Santarelli, M.; Dutton, R.J. Cheese rind communities provide tractable systems for in situ and in vitro studies of microbial diversity. Cell 2014, 158, 422–433.
  128. De Pasquale, I.; Di Cagno, R.; Buchin, S.; De Angelis, M.; Gobbetti, M. Spatial distribution of the metabolically active microbiota within Italian PDO ewes’ milk cheeses. PLoS ONE 2016, 11, e0153213.
  129. Turri, F.; Cremonesi, P.; Battelli, G.; Severgnini, M.; Brasca, M.; Gandini, G.; Pizzi, F. High biodiversity in a limited mountain area revealed in the traditional production of Historic Rebel cheese by an integrated microbiota–lipidomic approach. Sci. Rep. 2021, 11, 10374.
  130. Duru, I.C.; Laine, P.; Andreevskaya, M.; Paulin, L.; Kananen, S.; Tynkkynen, S.; Auvinen, P.; Smolander, O.P. Metagenomic and metatranscriptomic analysis of the microbial community in Swiss-type Maasdam cheese during ripening. Int. J. Food Microbiol. 2018, 281, 10–22.
  131. De Filippis, F.; Genovese, A.; Ferranti, P.; Gilbert, J.A.; Ercolini, D. Metatranscriptomics reveals temperature-driven functional changes in microbiome impacting cheese maturation rate. Sci. Rep. 2016, 6, 21871.
  132. Dugat-Bony, E.; Straub, C.; Teissandier, A.; Onésime, D.; Loux, V.; Monnet, C.; Irlinger, F.; Landaud, S.; Leclercq-Perlat, M.-N.; Bento, P.; et al. Overview of a surface-ripened cheese community functioning by meta-omics analyses. PLoS ONE 2015, 10, e0124360.
More
Information
Contributor MDPI registered users' name will be linked to their SciProfiles pages. To register with us, please refer to https://encyclopedia.pub/register :
View Times: 570
Entry Collection: Environmental Sciences
Revisions: 2 times (View History)
Update Date: 18 Jan 2022
1000/1000
ScholarVision Creations