Growth and adaptation of microorganisms on the cheese surface

Microbiome Assembly in Fermented Foods

For thousands of years, humans have enjoyed the novel flavors, increased shelf-life, and nutritional benefits that microbes provide in fermented foods and beverages. Recent sequencing surveys of ferments have mapped patterns of microbial diversity across space, time, and production practices. But a mechanistic understanding of how fermented food microbiomes assemble has only recently begun to emerge. Using three foods as case studies (surface-ripened cheese, sourdough starters, and fermented vegetables), we use an ecological and evolutionary framework to identify how microbial communities assemble in ferments. By combining in situ sequencing surveys with in vitro models, we are beginning to understand how dispersal, selection, diversification, and drift generate the diversity of fermented food communities. Most food producers are unaware of the ecological processes occurring in their production environments, but the theory and models of ecology and evolution can provide new approaches for managing fermented food microbiomes, from farm to ferment.

Longitudinal, Multi-Platform Metagenomics Yields a High-Quality Genomic Catalog and Guides an In Vitro Model for Cheese Communities

Microbiomes are intricately intertwined with human health, geochemical cycles, and food production. While many microbiomes of interest are highly complex and experimentally intractable, cheese rind microbiomes have proven to be powerful model systems for the study of microbial interactions. To provide a more comprehensive view of the genomic potential and temporal dynamics of cheese rind communities, we combined longitudinal, multi-platform metagenomics of three ripening washed-rind cheeses with whole-genome sequencing of community isolates. Sequencing-based approaches revealed a highly reproducible microbial succession in each cheese and the coexistence of closely related Psychrobacter species and enabled the prediction of plasmid and phage diversity and their host associations. In combination with culture-based approaches, we established a genomic catalog and a paired 16-member in vitro washed-rind cheese system. The combination of multi-platform metagenomic time-series data and an in vitro model provides a rich resource for further investigation of cheese rind microbiomes both computationally and experimentally. IMPORTANCE Metagenome sequencing can provide great insights into microbiome composition and function and help researchers develop testable hypotheses. Model microbiomes, such as those composed of cheese rind bacteria and fungi, allow the testing of these hypotheses in a controlled manner. Here, we first generated an extensive longitudinal metagenomic data set. This data set reveals successional dynamics, yields a phyla-spanning bacterial genomic catalog, associates mobile genetic elements with their hosts, and provides insights into functional enrichment of Psychrobacter in the cheese environment. Next, we show that members of the washed-rind cheese microbiome lend themselves to in vitro community reconstruction. This paired metagenomic data and in vitro system can thus be used as a platform for generating and testing hypotheses related to the dynamics within, and the functions associated with, cheese rind microbiomes.

Causality Verification for the Correlation between the Presence of Nonstarter Bacteria and Flavor Characteristics in Soft-Type Ripened Cheeses

Flavor characteristics of ripened cheese are established by various bacteria, such as lactic acid bacteria, Actinobacteria, and Proteobacteria, which spontaneously develop during the cheese-manufacturing process. We previously revealed the relationship between bacterial microbiota and flavor components in soft-type ripened cheeses by using a multiomics approach that combined metagenomics and metabolomics; however, we could not establish a causal relationship. This study aimed to substantiate the causal nature of the correlations revealed by the multiomics approach by using cheese-ripening tests with single isolate inoculation. The bacterial diversity and composition in surface mold-ripened cheeses from Japan and France varied, depending on the differences between the milks (pasteurized or raw), cheese positions (core or rind), and manufacturers. Although the volatile compounds did not clearly reflect the distinctive characteristics of the cheese samples, nonstarter lactic acid bacteria, Actinobacteria, and Proteobacteria positively correlated with ketones and sulfur compounds, as evidenced by a Spearman's correlation analysis. Cheese-ripening tests conducted after inoculation with single bacterial strains belonging to the above-mentioned taxa confirmed that these bacteria formed volatile compounds, in agreement with the correlations observed. In particular, various flavor compounds, such as acids, esters, ketones, and sulfur compounds, were detected in cheese inoculated with Pseudoalteromonas sp. TS-4-4 strain. These findings provide important insights into the role of nonstarter bacteria in the development of cheese flavor and into the effectiveness of the multiomics approach in screening for bacteria that can improve the quality of cheese products. IMPORTANCE Our previous study revealed that the existence of various bacteria, such as lactic acid bacteria, Actinobacteria, and Proteobacteria, clearly correlated with the abundance of flavor components, such as volatile compounds, in soft-type ripened cheeses via a multiomics approach that used 16S rRNA gene amplicon sequencing and headspace gas chromatography-mass spectrometry. However, this approach only showed correlations derived from statistical analyses rather than causal relationships. Therefore, in the present study, we performed cheese-ripening tests using nonstarter bacteria to substantiate the correlations revealed by the multiomics approach in soft-type ripened cheese. Our results suggest the capability of nonstarter bacteria, such as Proteobacteria, to impart flavor to cheese and the effectiveness of the multiomics approach in screening for microbial isolates that can improve the quality of cheese. Overall, our research provides new insights into the importance of bacteria in cheese production.

Analysis of Microbiota Persistence in Quebec's Terroir Cheese Using a Metabarcoding Approach

Environmental short amplicon sequencing, or metabarcoding, is commonly used to characterize the bacterial and fungal microbiota of cheese. Comparisons between different metabarcoding studies are complicated by the use of different gene markers. Here, we systematically compare different metabarcoding molecular targets using V3–V4 and V6–V8 regions of the bacterial 16S rDNA and fungal ITS1 and ITS2 regions. Taxonomic profiles varied depending on the molecular markers used. Based on data quality and detection capacity of the markers toward microorganisms usually associated with the dairy environment, the ribosomal regions V3–V4 and ITS2 were selected and further used to evaluate variability in the microbial ecosystem of terroir cheeses from the province of Quebec in Canada. Both fungal and bacterial ecosystem profiles were described for 32 different ready-to-eat bloomy-, washed- and natural-rind specialty cheese varieties. Among them, 15 were studied over two different production years. Using the Bray–Curtis dissimilarity index as an indicator of microbial shifts, we found that most variations could be explained by either a voluntary change in starter or ripening culture composition, or by changes in the cheesemaking technology. Overall, our results suggest the persistence of the microbiota between the two years studied—these data aid understanding of cheese microbiota composition and persistence during cheese ripening.

Metatranscriptomic Analyses Unravel Dynamic Changes in the Microbial and Metabolic Transcriptional Profiles in Artisanal Austrian Hard-Cheeses During Ripening

Vorarlberger Bergkäse (VB) is an artisanal Austrian washed-rind hard cheese produced from alpine cows' raw milk without the addition of ripening cultures. Ripening time is a key factor in VB, as it strongly influences the microbial communities present in the cheeses and the organoleptic properties of the product. In this study, the microbial and metabolic transcriptional profiles in VB rinds at different ripening times were investigated. VB products before (30 days of ripening) and after (90 days of ripening) selling were selected, RNA was extracted and subjected to shotgun metatranscriptomic sequencing. The analysis revealed some of the previously described abundant bacterial taxa of Brevibacterium, Corynebacterium, Halomonas, Psychrobacter, and Staphylococcus to be highly active in VB rinds. Additionally, the investigation of most important metabolic pathways in cheese ripening clearly showed differences in the gene transcription profiles and the active microbiota between the two ripening points investigated. At 30 days of ripening, metabolic events related with the degradation of residual lactose, lactate, citrate, proteolysis, and lipolysis were significantly more transcribed and mainly associated with Staphylococcus. On the other hand, genes involved in the degradation of smaller compounds derived from previous metabolism (i.e., metabolism of free amino acids and fatty acids) were significantly more expressed in VB rinds with 90 of ripening, and mainly associated with Brevibacterium and Corynebacterium. These latter metabolic activities are responsible of the generation of compounds, such as methanethiol and 2,3-butanediol, that are very important for the flavor and aroma characteristics of cheeses. This study shows the dynamic changes in the gene transcriptional profiles associated with energy substrates metabolism and the generation of organoleptic compounds during VB ripening and uncovers bacterial taxa as key drivers of the ripening process. These taxa might be the target for future studies toward an accelerated cheese ripening and the enhancement of its organoleptic properties.

A Review on the General Cheese Processing Technology, Flavor Biochemical Pathways and the Influence of Yeasts in Cheese

Cheese has a long history and this naturally fermented dairy product contains a range of distinctive flavors. Microorganisms in variety cheeses are an essential component and play important roles during both cheese production and ripening. However, cheeses from different countries are still handmade, the processing technology is diverse, the microbial community structure is complex and the cheese flavor fluctuates greatly. Therefore, studying the general processing technology and relationship between microbial structure and flavor formation in cheese is the key to solving the unstable quality and standardized production of cheese flavor on basis of maintaining the flavor of cheese. This paper reviews the research progress on the general processing technology and key control points of natural cheese, the biochemical pathways for production of flavor compounds in cheeses, the diversity and the role of yeasts in cheese. Combined with the development of modern detection technology, the evolution of microbial structure, population evolution and flavor correlation in cheese from different countries was analyzed, which is of great significance for the search for core functional yeast microorganisms and the industrialization prospect of traditional fermented cheese.

Unraveling the world of halophilic and halotolerant bacteria in cheese by combining cultural, genomic and metagenomic approaches

Halophilic/halotolerant bacteria are generally assumed to live in natural environments, although they may also be found in foods such as cheese and seafood. These salt-loving bacteria have been occasionally characterized in cheese, and studies on their ecological and technological functions are still scarce. We therefore selected 13 traditional cheeses to systematically characterize these microorganisms in their rinds via cultural, genomic and metagenomic methods. Using different salt-based media, we identified 35 strains with unique 16S rRNA and rpoB gene sequences, whose whole genome was sequenced. Twenty are Gram-positive species including notably Brevibacterium aurantiacum (6) and Staphylococcus equorum (3), which are also frequently added as starters. ANI and pan-genomic analyses confirm the high genetic diversity of B. aurantiacum and reveal the presence of two subspecies in S. equorum, as well as the genetic proximity of several cheese strains to bovine isolates. Additionally, we isolated 15 Gram-negative strains, potentially defining ten new species of halophilic/halotolerant cheese bacteria, in particular for the genera Halomonas and Psychrobacter. The use of all the genomes sequenced in this study as a reference to complement those existing in the databases allowed us to study the representativeness of 66 species of halophilic/halotolerant bacteria in 74 cheese rind metagenomes. While Gram-positive strains may flourish in the different types of technologies, Gram-negative species are particularly abundant in cheeses with high moisture, such as washed-rind cheeses. Finally, analyses of co-occurrences reveal assemblies, including the frequent coexistence of several species of the same genus, forming moderately complex ecosystems with functional redundancies that probably ensure stable cheese development.

Evaluation of the Relationships Between Microbiota and Metabolites in Soft-Type Ripened Cheese Using an Integrated Omics Approach

Cheese ripening is effected by various microorganisms and results in the characteristic flavors of cheese. Owing to the complexity of the microbiota involved, the relationship between microorganisms and components during ripening remains unclear. In this study, metagenomics and metabolomics were integrated to reveal these relationships in three kinds of surface mold-ripened cheeses and two kinds of bacterial smear-ripened cheeses. The microbiota is broadly divided into two groups to correspond with different cheese types. Furthermore, surface mold-ripened cheese showed similar microbiota regardless of the cheese variety, whereas bacterial smear-ripened cheese showed specific microbiota characterized by marine bacteria (MB) and halophilic and alkaliphilic lactic acid bacteria for each cheese variety. In the metabolite analysis, volatile compounds suggested differences in cheese types, although organic acids and free amino acids could not determine the cheese characteristics. On the other hand, Spearman correlation analysis revealed that the abundance of specific bacteria was related to the formation of specific organic acids, free amino acids, and volatile compounds. In particular, MB was positively correlated with esters and pyrazines, indicating their contribution to cheese quality. These methodologies and results further our understanding of microorganisms and allow us to select useful strains for cheese ripening.

Molecular Characterization and Biofilm Formation Study of Contaminant Bacteria Isolated from Domiaty and Hungarian Cheeses in Jeddah City

The aim was to study the microbiological quality of Domiaty and Hungarian cheeses, molecular identification and biofilm formation of some selected contaminant bacteria. Samples were collected from two M and P big markets in Jeddah City through the period from February to October 2018, nine visits for two types of natural cheese. Results showed that the total bacterial counts (CFU/ml) from Domiaty cheese from two markets (M and P) were 0.1 x 105, 8 x 105 and 1 x 10 5 CFU/ml respectively (3 visits of M market) and 4 x 106, 0.4 x 106, 6.5 x 103, 1 x 103, 0.1 x 103 and 0.1 x 103 CFU/ml respectively (six samples from 6 visits from P market). Results showed that the total bacterial counts (CFU/ml) from Hungarian cheese were 1.5 x 10 5, 1x 10 4, 11 x 10 4 and 4 x10 6 CFU/ml respectively from (4 visits of M market) and 0.18 x 104, 3 x 106, 22 x 106, 6 x 106 and 5 x 104 CFU/ml respectively (5 visits from P market).Different bacterial isolates from cheese were identified by morphology and biochemical test. Bacterial isolates from cheeses were identified by VITEK MS as follow: Serratia liquefaciens (D6-1, D6-2, D14-1, D13-1 and D13-2), and Pseudomonas fluorescens (D14-2) were isolated from Domiaty cheese while Enterococcus faecium (H11-2), Serratia liquefaciens (H15-1) and Streptococcus thermophilus (H14-1) were isolated from Hungarian cheese. Some selected bacterial isolates were identified by 16S rRNA. Isolates were belong to MK757978 (Raoultilla terrigena (D15-1)), MK757979 (Bacillus cereus (D16-1)), MK757980 (Enterococcus faecalis (H10-2)), MK757982 (Enterococcus fiscalism (H11-1)), MK757981 (Serratia liquefactions (H13-1)), MK757984 (Anoxybacillus flavithermus (H17-1). All bacterial isolates have been tested for the formation of biofilm using a Tissue Culture Plate (TCP). Results revealed 12.5% and 46.15% of high biofilm formation respectively for bacterial isolates of Domiaty and Hungarian cheeses.

Austrian Raw-Milk Hard-Cheese Ripening Involves Successional Dynamics of Non-Inoculated Bacteria and Fungi

Cheese ripening involves successional changes of the rind microbial composition that harbors a key role on the quality and safety of the final products. In this study, we analyzed the evolution of the rind microbiota (bacteria and fungi) throughout the ripening of Austrian Vorarlberger Bergkäse (VB), an artisanal surface-ripened cheese, by using quantitative and qualitative approaches. The real-time quantitative PCR results revealed that bacteria were more abundant than fungi in VB rinds throughout ripening, although both kingdoms were abundant along the process. The qualitative investigation was performed by high-throughput gene-targeted (amplicon) sequencing. The results showed dynamic changes of the rind microbiota throughout ripening. In the fresh products, VB rinds were dominated by Staphylococcus equorum and Candida. At early ripening times (14–30 days) Psychrobacter and Debaryomyces flourished, although their high abundance was limited to these time points. At the latest ripening times (90–160 days), VB rinds were dominated by S. equorum, Brevibacterium, Corynebacterium, and Scopulariopsis. Strong correlations were shown for specific bacteria and fungi linked to specific ripening periods. This study deepens our understanding of VB ripening and highlights different bacteria and fungi associated to specific ripening periods which may influence the organoleptic properties of the final products.

Strain-Level Diversity Impacts Cheese Rind Microbiome Assembly and Function

Diversification can generate genomic and phenotypic strain-level diversity within microbial species. This microdiversity is widely recognized in populations, but the community-level consequences of microbial strain-level diversity are poorly characterized. Using the cheese rind model system, we tested whether strain diversity across microbiomes from distinct geographic regions impacts assembly dynamics and functional outputs. We first isolated the same three bacterial species (Staphylococcus equorum, Brevibacterium auranticum, and Brachybacterium alimentarium) from nine cheeses produced in different regions of the United States and Europe to construct nine synthetic microbial communities consisting of distinct strains of the same three bacterial species. Comparative genomics identified distinct phylogenetic clusters and significant variation in genome content across the nine synthetic communities. When we assembled each synthetic community with initially identical compositions, community structure diverged over time, resulting in communities with different dominant taxa. The taxonomically identical communities showed differing responses to abiotic (high salt) and biotic (the fungus Penicillium) perturbations, with some communities showing no response and others substantially shifting in composition. Functional differences were also observed across the nine communities, with significant variation in pigment production (light yellow to orange) and in composition of volatile organic compound profiles emitted from the rinds (nutty to sulfury).IMPORTANCE Our work demonstrated that the specific microbial strains used to construct a microbiome could impact the species composition, perturbation responses, and functional outputs of that system. These findings suggest that 16S rRNA gene taxonomic profiles alone may have limited potential to predict the dynamics of microbial communities because they usually do not capture strain-level diversity. Observations from our synthetic communities also suggest that strain-level diversity has the potential to drive variability in the aesthetics and quality of surface-ripened cheeses.

Attenuated Lactococcus lactis and Surface Bacteria as Tools for Conditioning the Microbiota and Driving the Ripening of Semisoft Caciotta Cheese

This study aimed at establishing the effects of attenuated starters and surface bacteria on various features of caciotta cheese. The cheese undergoes a ripening period during which the house microbiota contaminates the surface. Conventional cheese (the control cheese [CC]) is made using only primary starters. Primary starters and attenuated (i.e., unable to grow and synthesize lactic acid) Lactococcus lactis (Lc. lactis) subsp. lactis were used to produce caciotta cheese without (ATT cheese) or with an inoculum of surface bacteria: (i) Leuconostoc lactis (Le. lactis) (LL cheese), (ii) Vibrio casei (VC cheese), (iii) Staphylococcus equorum (SE cheese), (iv) Brochothrix thermosphacta (BX cheese), and (v) a mixture of all four (MIX cheese). Attenuated Lc. lactis increased microbial diversity during cheese ripening. At the core, attenuated starter mainly increased indigenous lactococci and Lactobacillus delbrueckii group bacteria. At the surface, the main effect was on Macrococcus caseolyticus Autochthonous Le. lactis strains took advantage of the attenuated starter, becoming dominant. Adjunct Le. lactis positively affected Lactobacillus sakei group bacteria on the LL cheese surface. Adjunct V. casei, S. equorum, and B. thermosphacta did not become dominant. Surfaces of VC, SE, and BX cheeses mainly harbored Staphylococcus succinus Peptidase activities were higher in cheeses made with attenuated starter than in CC, which had the lowest concentration of free amino acids. Based on the enzymatic activities of adjunct Le. lactis, LL and MIX cheeses exhibited the highest glutamate dehydrogenase, cystathionine-γ-lyase, and esterase activities. As shown by multivariate statistical analyses, LL and MIX cheeses showed the highest similarity for microbiological and biochemical features. LL and MIX cheeses received the highest scores for overall sensory acceptability.IMPORTANCE This study provides in-depth knowledge of the effects of attenuated starters and surface bacterial strains on the microbiota and related metabolic activities during cheese ripening. The use of attenuated Lc. lactis strongly impacted the microbiota assembly of caciotta cheese. This led to improved biochemical and sensory features compared to conventional cheese. Among surface bacterial strains, Le. lactis played a key role in the metabolic activities involved in cheese ripening. This resulted in an improvement of the sensory quality of caciotta cheese. The use of attenuated lactic acid bacteria and the surface adjunct Le. lactis could be a useful biotechnology to improve the flavor formation of caciotta cheese.

Experimental Microbiomes: Models Not to Scale

Low-cost, high-throughput nucleic acid sequencing ushered the field of microbial ecology into a new era in which the microbial composition of nearly every conceivable environment on the planet is under examination. However, static "screenshots" derived from sequence-only approaches belie the underlying complexity of the microbe-microbe and microbe-host interactions occurring within these systems. Reductionist experimental models are essential to identify the microbes involved in interactions and to characterize the molecular mechanisms that manifest as complex host and environmental phenomena. Herein, we focus on three models (Bacillus-Streptomyces, Aliivibrio fischeri-Hawaiian bobtail squid, and gnotobiotic mice) at various levels of taxonomic complexity and experimental control used to gain molecular insight into microbe-mediated interactions. We argue that when studying microbial communities, it is crucial to consider the scope of questions that experimental systems are suited to address, especially for researchers beginning new projects. Therefore, we highlight practical applications, limitations, and tradeoffs inherent to each model.

Mobilome of Brevibacterium aurantiacum Sheds Light on Its Genetic Diversity and Its Adaptation to Smear-Ripened Cheeses

Brevibacterium aurantiacum is an actinobacterium that confers key organoleptic properties to washed-rind cheeses during the ripening process. Although this industrially relevant species has been gaining an increasing attention in the past years, its genome plasticity is still understudied due to the unavailability of complete genomic sequences. To add insights on the mobilome of this group, we sequenced the complete genomes of five dairy Brevibacterium strains and one non-dairy strain using PacBio RSII. We performed phylogenetic and pan-genome analyses, including comparisons with other publicly available Brevibacterium genomic sequences. Our phylogenetic analysis revealed that these five dairy strains, previously identified as Brevibacterium linens, belong instead to the B. aurantiacum species. A high number of transposases and integrases were observed in the Brevibacterium spp. strains. In addition, we identified 14 and 12 new insertion sequences (IS) in B. aurantiacum and B. linens genomes, respectively. Several stretches of homologous DNA sequences were also found between B. aurantiacum and other cheese rind actinobacteria, suggesting horizontal gene transfer (HGT). A HGT region from an iRon Uptake/Siderophore Transport Island (RUSTI) and an iron uptake composite transposon were found in five B. aurantiacum genomes. These findings suggest that low iron availability in milk is a driving force in the adaptation of this bacterial species to this niche. Moreover, the exchange of iron uptake systems suggests cooperative evolution between cheese rind actinobacteria. We also demonstrated that the integrative and conjugative element BreLI (Brevibacterium Lanthipeptide Island) can excise from B. aurantiacum SMQ-1417 chromosome. Our comparative genomic analysis suggests that mobile genetic elements played an important role into the adaptation of B. aurantiacum to cheese ecosystems.

Alcohol dehydrogenases from Proteus mirabilis contribute to alcoholic flavor

BACKGROUND

Cheese ripening involves a complex series of metabolic reactions and numerous concomitant secondary transformations. Alcohol dehydrogenase (ADH) converts aldehydes into their corresponding alcohols, which enrich cheese aroma.

RESULTS

In this study, we identified five ADH genes in Proteus mirabilis JN458, and these genes were overexpressed and characterized in Escherichia coli BL21 (DE3). The optimum pH was 7.0 for the purified recombinant ADH-1, ADH-2, and ADH-3 and 8.0 for ADH-4 and ADH-5. The optimum temperature was 40 °C for ADH-1, ADH-3, and ADH-5 and 45 °C for ADH-2 and ADH-4. The K_m value of ADH-1, ADH-2, and ADH-3 was 34.45, 16.90, and 10.01 µmol L^-1 for phenylacetaldehyde, respectively. The K_m value of ADH-4 and ADH-5 was 14.81 and 24.62 µmol L^-1 for 2-methylbutanal, respectively.

CONCLUSION

Proteus species play important roles during cheese ripening. The results of our study are important for further research on cheese flavor and for quality control during cheese production. © 2019 Society of Chemical Industry.

Brevibacterium from Austrian hard cheese harbor a putative histamine catabolism pathway and a plasmid for adaptation to the cheese environment

The genus Brevibacterium harbors many members important for cheese ripening. We performed real-time quantitative PCR (qPCR) to determine the abundance of Brevibacterium on rinds of Vorarlberger Bergkäse, an Austrian artisanal washed-rind hard cheese, over 160 days of ripening. Our results show that Brevibacterium are abundant on Vorarlberger Bergkäse rinds throughout the ripening time. To elucidate the impact of Brevibacterium on cheese production, we analysed the genomes of three cheese rind isolates, L261, S111, and S22. L261 belongs to Brevibacterium aurantiacum, whereas S111 and S22 represent novel species within the genus Brevibacterium based on 16S rRNA gene similarity and average nucleotide identity. Our comparative genomic analysis showed that important cheese ripening enzymes are conserved among the genus Brevibacterium. Strain S22 harbors a 22 kb circular plasmid which encodes putative iron and hydroxymethylpyrimidine/thiamine transporters. Histamine formation in fermented foods can cause histamine intoxication. We revealed the presence of a putative metabolic pathway for histamine degradation. Growth experiments showed that the three Brevibacterium strains can utilize histamine as the sole carbon source. The capability to utilize histamine, possibly encoded by the putative histamine degradation pathway, highlights the importance of Brevibacterium as key cheese ripening cultures beyond their contribution to cheese flavor production.

Design of microbial consortia for the fermentation of pea-protein-enriched emulsions

In order to encourage Western populations to increase their consumption of vegetables, we suggest turning legumes into novel, healthy foods by applying an old, previously widespread method of food preservation: fermentation. In the present study, a total of 55 strains from different microbial species (isolated from cheese or plants) were investigated for their ability to: (i) grow on a emulsion containing 100% pea proteins and no carbohydrates or on a 50:50 pea:milk protein emulsion containing lactose, (ii) increase aroma quality and reduce sensory off-flavors; and (iii) compete against endogenous micro organisms. The presence of carbohydrates in the mixed pea:milk emulsion markedly influenced the fermentation by strongly reducing the pH through lactic fermentation, whereas the absence of carbohydrates in the pea emulsion promoted alkaline or neutral fermentation. Lactic acid bacteria assigned to Lactobacillus plantarum, Lactobacillus rhamnosus, Lactococcus lactis and Lactobacillus casei species grew well in both the pea and pea:milk emulsions. Most of the fungal strains tested (particularly those belonging to the Mucor and Geotrichum genera) were also able to grow on both emulsions. Although most Actinobacteria and Proteobacteria did not compete with endogenous microbiota (Bacillus), some species such as Hafnia alvei, Acinetobacter johnsonii and Glutamicibacter arilaitensis grew strongly and appeared to restrict the development of the endogenous microbiota when the pea emulsion was inoculated with a combination of three to nine strains. In the mixed emulsions, lactic fermentation inhibited Actinobacteria and Proteobacteria (e.g. Brevibacterium casei, Corynebacterium casei, Staphylococcus lentus) to the greatest extent but also inhibited Bacillus (e.g. Bacillus subtilis and Bacillus licheniformis). Overall, this procedure enabled us to select two microbial consortia able to colonize pea-based products and positively influence the release of volatile compounds by generating a roasted/grilled aroma for the 100% pea emulsion, and a fruity, lactic aroma for the 50:50 pea:milk emulsion. Moreover, the fermentation in the pea-based emulsions reduced the level of hexanal, which otherwise leads to an undesired green pea aroma. Our present results show how the assembly of multiple microbial cultures can help to develop an innovative food product.

Geotrichum candidum gene expression and metabolite accumulation inside the cells reflect the strain oxidative stress sensitivity and ability to produce flavour compounds

Geotrichum candidum is a fungus-like yeast widely used as a starter culture for cheese ripening for its proteolytic and lipolytic activities and its contribution to the cheese flavours. The sequenced strain G. candidum CLIB 918 was isolated from cheese Pont-L’Evêque. This strain's ability to produce volatile compounds was compared to the ability of a known strong sulphur compound producer G. candidum strain (Gc203). The aminotransferase-coding genes BAT2 and ARO8 were identified to be involved in methionine catabolism. The production of volatile compounds indicated that the sequenced strain was a moderate producer compared to the strong producer strain. The major volatile compounds were produced from sulphur amino acid, branched-chain amino acid and fatty acid metabolisms. Metabolite content of the cells showed that the ability of the strain to produce volatile compounds was inversely proportional to its ability to store amino acids inside the cells. Reduced glutathione, hypotaurine and taurine intracellular concentrations and volatile fatty aldehyde production indicated the role of oxidative stress sensitivity in flavour production. The increase in expression of several genes in a Reblochon-type cheese at the end of ripening confirmed that oxygen and iron were key factors regulating cheese flavour production.

Characterization of Bacterial Communities in Mexican Artisanal Raw Milk "Bola de Ocosingo" Cheese by High-Throughput Sequencing

The dynamics of bacteria community of "Bola de Ocosingo" cheese, a Mexican artisanal raw milk cheese was investigated by high-throughput sequencing (454 pyrosequencing). Dairy samples (raw milk, curd, cheese at 50 and 110 days of ripening) were collected at dry (March-June) and rainy season (August-November) from three producers located in Chiapas, Mexico. In general, raw milk contained high bacterial diversity which was reduced throughout cheese manufacture. However, in two productions an important increase during cheese ripening was observed probably due to cross-contamination. Species such as Streptococcus thermophilus, Lactococcus lactis, Lactobacillus helveticus, L. delbrueckii and L. plantarum from which potential probiotic strains may be obtained, predominated during processing, varying its prevalence from one producer to another. Furthermore, low proportions of Escherichia coli/Shigella flexnerii were detected in almost all processes, however, could not be recovered by traditional methodology, indicating presence of non-cultivable cells. This work provides insights into bacteria communities of Bola de Ocosingo cheese for starter culture development, many of which are reported to provide health related benefits, and the usefulness of high-throughput sequencing to evidence cross-contamination during processing.

Deciphering intra-species bacterial diversity of meat and seafood spoilage microbiota using gyrB amplicon sequencing: A comparative analysis with 16S rDNA V3-V4 amplicon sequencing

Meat and seafood spoilage ecosystems harbor extensive bacterial genomic diversity that is mainly found within a small number of species but within a large number of strains with different spoilage metabolic potential. To decipher the intraspecies diversity of such microbiota, traditional metagenetic analysis using the 16S rRNA gene is inadequate. We therefore assessed the potential benefit of an alternative genetic marker, gyrB, which encodes the subunit B of DNA gyrase, a type II DNA topoisomerase. A comparison between 16S rDNA-based (V3-V4) amplicon sequencing and gyrB-based amplicon sequencing was carried out in five types of meat and seafood products, with five mock communities serving as quality controls. Our results revealed that bacterial richness in these mock communities and food samples was estimated with higher accuracy using gyrB than using16S rDNA. However, for Firmicutes species, 35% of putative gyrB reads were actually identified as sequences of a gyrB paralog, parE, which encodes subunit B of topoisomerase IV; we therefore constructed a reference database of published sequences of both gyrB and pare for use in all subsequent analyses. Despite this co-amplification, the deviation between relative sequencing quantification and absolute qPCR quantification was comparable to that observed for 16S rDNA for all the tested species. This confirms that gyrB can be used successfully alongside 16S rDNA to determine the species composition (richness and evenness) of food microbiota. The major benefit of gyrB sequencing is its potential for improving taxonomic assignment and for further investigating OTU richness at the subspecies level, thus allowing more accurate discrimination of samples. Indeed, 80% of the reads of the 16S rDNA dataset were represented by thirteen 16S rDNA-based OTUs that could not be assigned at the species-level. Instead, these same clades corresponded to 44 gyrB-based OTUs, which differentiated various lineages down to the subspecies level. The increased ability of gyrB-based analyses to track and trace phylogenetically different groups of strains will generate improved resolution and more reliable results for studies of the strains implicated in food processes.

Cheese brines from Danish dairies reveal a complex microbiota comprising several halotolerant bacteria and yeasts

The Danish Danbo cheese is a surface ripened semi-hard cheese, which before ripening is submerged in brine for up to 24 h. The brining is required in order to obtain the structural and organoleptic properties of the cheeses. Likewise, the content of NaCl in the cheese will influence especially the surface microbiota being of significant importance for flavour development and prevention of microbial spoilage. Even though the microbiota on cheese surfaces have been studied extensively, limited knowledge is available on the occurrence of microorganisms in cheese brine. The aim of the present study was to investigate by both culture-dependent and -independent techniques the brine microbiota in four Danish dairies producing Danbo cheese. The pH of the brines varied from 5.1 to 5.6 with a dry matter content from 20 to 27% (w/w). The content of lactate varied from 4.1 to 10.8 g/L and free amino acids from 65 to 224 mg/L. Bacteria were isolated on five different media with NaCl contents of 0.85-23.0% (w/v) NaCl. The highest count of 6.3 log CFU/mL was obtained on TSA added 4% (w/v) NaCl. For yeasts, the highest count was 3.7 log CFU/mL on MYGP added 8% (w/v) NaCl. A total of 31 bacterial and eight eukaryotic species were isolated including several halotolerant and/or halophilic species. Among bacteria, counts of ≥6.0 log CFU/mL were obtained for Tetragenococcus muriaticus and Psychrobacter celer, while counts between ≥4.5 and < 6.0 log CFU/mL were obtained for Lactococcus lactis, Staphylococcus equorum, Staphylococcus hominis, Chromohalobacter beijerinckii, Chromohalobacter japonicus and Microbacterium maritypicum. Among yeasts, counts of ≥3.5 log CFU/mL were only obtained for Debaryomyces hansenii. By amplicon-based high-throughput sequencing of 16S rRNA gene and ITS2 regions for bacteria and eukaryotes respectively, brines from the same dairy clustered together indicating the uniqueness of the dairy brine microbiota. To a great extent the results obtained by amplicon sequencing fitted with the culture-dependent technique though each of the two methodologies identified unique genera/species. Dairy brine handling procedures as e.g. microfiltration were found to influence the brine microbiota. The current study proves the occurrence of a specific dairy brine microbiota including several halotolerant and/or halophilic species most likely of sea salt origin. The importance of these species during especially the initial stages of cheese ripening and their influence on cheese quality and safety need to be investigated. Likewise, optimised brine handling procedures and microbial cultures are required to ensure an optimal brine microbiota.

Coproporphyrin III Produced by the Bacterium Glutamicibacter arilaitensis Binds Zinc and Is Upregulated by Fungi in Cheese Rinds

Bacterium-fungus interactions play key roles in the assembly of cheese rind microbial communities, but the molecular mechanisms underlying these interactions are poorly characterized. Moreover, millions of people around the world enjoy eating cheeses and cheese rinds, but our understanding of the diversity of microbial metabolites ingested during cheese consumption is limited. The discovery of zinc coproporphyrin III as the cause of pink pigment production by Glutamicibacter arilaitensis suggests that secretion of this molecule is important for microbial acquisition of trace metals.

Microbial communities of fermented food microbiomes typically exhibit predictable patterns of microbial succession. However, the biochemical mechanisms that control the diversity and dynamics of these communities are not well described. Interactions between bacteria and fungi may be one mechanism controlling the development of cheese rind microbiomes. This study characterizes a specific bacterium-fungus interaction previously discovered on cheese rinds between the bacterium Glutamicibacter arilaitensis (formerly Arthrobacter arilaitensis) and fungi of the genus Penicillium and identifies the specialized metabolites produced during cocultures. G. arilaitensis was previously shown to produce an unknown pink pigment in response to the presence of Penicillium. Using a combination of mass spectrometry, nuclear magnetic resonance (NMR), and transcriptome sequencing (RNA-seq), we determined that this pigment production is associated with production of coproporphyrin III. The discovery that coproporphyrin III preferentially bound zinc over other trace metals found in cheese curds highlights the value of using analytical chemistry to confirm identity of predicted chemical species.

IMPORTANCE Bacterium-fungus interactions play key roles in the assembly of cheese rind microbial communities, but the molecular mechanisms underlying these interactions are poorly characterized. Moreover, millions of people around the world enjoy eating cheeses and cheese rinds, but our understanding of the diversity of microbial metabolites ingested during cheese consumption is limited. The discovery of zinc coproporphyrin III as the cause of pink pigment production by Glutamicibacter arilaitensis suggests that secretion of this molecule is important for microbial acquisition of trace metals.

Author Video: An author video summary of this article is available.

Omics-Based Insights into Flavor Development and Microbial Succession within Surface-Ripened Cheese

In this study, a young Cheddar curd was used to produce two types of surface-ripened cheese, using two commercial smear-culture mixes of yeasts and bacteria. Whole-metagenome shotgun sequencing was used to screen the microbial population within the smear-culture mixes and on the cheese surface, with comparisons of microorganisms at both the species and the strain level. The use of two smear mixes resulted in the development of distinct microbiotas on the surfaces of the two test cheeses. In one case, most of the species inoculated on the cheese established themselves successfully on the surface during ripening, while in the other, some of the species inoculated were not detected during ripening and the most dominant bacterial species, Glutamicibacter arilaitensis, was not a constituent of the culture mix. Generally, yeast species, such as Debaryomyces hansenii and Geotrichum candidum, were dominant during the first stage of ripening but were overtaken by bacterial species, such as Brevibacterium linens and G. arilaitensis, in the later stages. Using correlation analysis, it was possible to associate individual microorganisms with volatile compounds detected by gas chromatography-mass spectrometry in the cheese surface. Specifically, D. hansenii correlated with the production of alcohols and carboxylic acids, G. arilaitensis with alcohols, carboxylic acids and ketones, and B. linens and G. candidum with sulfur compounds. In addition, metagenomic sequencing was used to analyze the metabolic potential of the microbial populations on the surfaces of the test cheeses, revealing a high relative abundance of metagenomic clusters associated with the modification of color, variation of pH, and flavor development. IMPORTANCE Fermented foods, in particular, surface-ripened cheese, represent a model to explain the metabolic interactions which regulate microbial succession in complex environments. This study explains the role of individual species in a heterogeneous microbial environment, i.e., the exterior of surface-ripened cheese. Through whole-metagenome shotgun sequencing, it was possible to investigate the metabolic potential of the resident microorganisms and show how variations in the microbial populations influence important aspects of cheese ripening, especially flavor development. Overall, in addition to providing fundamental insights, this research has considerable industrial relevance relating to the production of fermented food with specific qualities.

Autochthonous facility-specific microbiota dominates washed-rind Austrian hard cheese surfaces and its production environment

Cheese ripening involves the succession of complex microbial communities that are responsible for the organoleptic properties of the final products. The food processing environment can act as a source of natural microbial inoculation, especially in traditionally manufactured products. Austrian Vorarlberger Bergkäse (VB) is an artisanal washed-rind hard cheese produced in the western part of Austria without the addition of external ripening cultures. Here, the composition of the bacterial communities present on VB rinds and on different processing surfaces from two ripening cellars was assessed by near full length 16S rRNA gene amplification, cloning and sequencing. Non-inoculated aerobic bacteria dominated all surfaces in this study. VB production conditions (long ripening time, high salt concentration and low temperatures) favor the growth of psychro- and halotolerant bacteria. Several bacterial groups, such as coryneforms, Staphylococcus equorum and Halomonas dominated VB and were also found on most environmental surfaces. Analysis of OTUs shared between different surfaces suggests that VB rind bacteria are inoculated naturally during the ripening from the processing environment and that cheese surfaces exert selective pressure on these communities, as only those bacteria better adapted flourished on VB rinds. This study analyzed VB processing environment microbiota and its relationship with VB rinds for the first time, elucidating that the processing environment and the cheese microbiota should be considered as microbiologically linked ecosystems with the goal of better defining the events that take place during cheese maturation.

Comparative genomic analysis of Brevibacterium strains: insights into key genetic determinants involved in adaptation to the cheese habitat

Background

Brevibacterium strains are widely used for the manufacturing of surface-ripened cheeses, contributing to the breakdown of lipids and proteins and producing volatile sulfur compounds and red-orange pigments. The objective of the present study was to perform comparative genomic analyses in order to better understand the mechanisms involved in their ability to grow on the cheese surface and the differences between the strains.

Results

The genomes of 23 Brevibacterium strains, including twelve strains isolated from cheeses, were compared for their gene repertoire involved in salt tolerance, iron acquisition, bacteriocin production and the ability to use the energy compounds present in cheeses. All or almost all the genomes encode the enzymes involved in ethanol, acetate, lactate, 4-aminobutyrate and glycerol catabolism, and in the synthesis of the osmoprotectants ectoine, glycine-betaine and trehalose. Most of the genomes contain two contiguous genes encoding extracellular proteases, one of which was previously characterized for its activity on caseins. Genes encoding a secreted triacylglycerol lipase or involved in the catabolism of galactose and D-galactonate or in the synthesis of a hydroxamate-type siderophore are present in part of the genomes. Numerous Fe³⁺/siderophore ABC transport components are present, part of them resulting from horizontal gene transfers. Two cheese-associated strains have also acquired catecholate-type siderophore biosynthesis gene clusters by horizontal gene transfer. Predicted bacteriocin biosynthesis genes are present in most of the strains, and one of the corresponding gene clusters is located in a probable conjugative transposon that was only found in cheese-associated strains.

Conclusions

Brevibacterium strains show differences in their gene repertoire potentially involved in the ability to grow on the cheese surface. Part of these differences can be explained by different phylogenetic positions or by horizontal gene transfer events. Some of the distinguishing features concern biotic interactions with other strains such as the secretion of proteases and triacylglycerol lipases, and competition for iron or bacteriocin production. In the future, it would be interesting to take the properties deduced from genomic analyses into account in order to improve the screening and selection of Brevibacterium strains, and their association with other ripening culture components.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4322-1) contains supplementary material, which is available to authorized users.

Abundance and potential contribution of Gram-negative cheese rind bacteria from Austrian artisanal hard cheeses

Many different Gram-negative bacteria have been shown to be present on cheese rinds. Their contribution to cheese ripening is however, only partially understood until now. Here, cheese rind samples were taken from Vorarlberger Bergkäse (VB), an artisanal hard washed-rind cheese from Austria. Ripening cellars of two cheese production facilities in Austria were sampled at the day of production and after 14, 30, 90 and 160days of ripening. To obtain insights into the possible contribution of Advenella, Psychrobacter, and Psychroflexus to cheese ripening, we sequenced and analyzed the genomes of one strain of each genus isolated from VB cheese rinds. Additionally, quantitative PCRs (qPCRs) were performed to follow the abundance of Advenella, Psychrobacter, and Psychroflexus on VB rinds during ripening in both facilities. qPCR results showed that Psychrobacter was most abundant on cheese rinds and the abundance of Advenella decreased throughout the first month of ripening and increased significantly after 30days of ripening (p<0.01). Psychrobacter and Psychroflexus increased significantly during the first 30 ripening days (p<0.01), and decreased to their initial abundance during the rest of the ripening time (p<0.05). Genome sequencing resulted in 17 to 27 contigs with assembly sizes of 2.7 Mbp for Psychroflexus, 3 Mbp for Psychrobacter, and 4.3 Mbp for Advenella. Our results reveal that each genome harbors enzymes shown to be important for cheese ripening in other bacteria such as: Cystathionine/Methionine beta or gamma-lyases, many proteases and peptidases (including proline iminopeptidases), aminotransferases, and lipases. Thus, all three isolates have the potential to contribute positively to cheese ripening. In conclusion, the three species quantified were stable community members throughout the ripening process and their abundance on cheese rinds together with the results from genome sequencing suggest an important contribution of these bacteria to cheese ripening.

Extensive horizontal gene transfer in cheese-associated bacteria

Acquisition of genes through horizontal gene transfer (HGT) allows microbes to rapidly gain new capabilities and adapt to new or changing environments. Identifying widespread HGT regions within multispecies microbiomes can pinpoint the molecular mechanisms that play key roles in microbiome assembly. We sought to identify horizontally transferred genes within a model microbiome, the cheese rind. Comparing 31 newly sequenced and 134 previously sequenced bacterial isolates from cheese rinds, we identified over 200 putative horizontally transferred genomic regions containing 4733 protein coding genes. The largest of these regions are enriched for genes involved in siderophore acquisition, and are widely distributed in cheese rinds in both Europe and the US. These results suggest that HGT is prevalent in cheese rind microbiomes, and that identification of genes that are frequently transferred in a particular environment may provide insight into the selective forces shaping microbial communities.

DOI:http://dx.doi.org/10.7554/eLife.22144.001

From the depths of the ocean to the lining of the human gut, almost every environment on Earth is home to a unique community of microorganisms referred to as a microbiome. Within these communities, unrelated microorganisms can exchange genetic information through a process known as horizontal gene transfer. For example, genes linked to antibiotic resistance are often transferred between different microorganisms, which can create increasingly drug resistant microbes and has important implications for human health.

Horizontal gene transfer has been studied for almost 100 years, but examining it directly is challenging because, almost by definition, it requires studying a community of microbes rather than one microbe in isolation. As such, researchers are looking for simple models of microbial communities that can be easily manipulated in experiments.

Bonham et al. have now turned to the outer surface of cheese, also known as cheese rind, to better understand horizontal gene transfer. As a model system, the cheese rind microbiome is relatively simple to work with because cheese rind is easy to replicate in the laboratory, and the microbes growing on cheese can be grown on their own or in combinations with other microbes. By comparing the genetic material of 165 cheese-associated bacteria to one another, Bonham et al. identified over 4,000 genes that were shared between the bacteria, including several large clusters of genes that were shared by many species.

Many of the identified genes (about 23% to be precise) appear to help the microorganisms acquire nutrients that are known to be in short supply on the surface of cheese surface, including iron. Bacteria typically use specialized molecules called siderophores to scavenge for iron and uptake systems to carry the iron-bound siderophore back into the cell. Notably, only the genes associated with the uptake systems were found in some of the shared gene clusters. This finding suggests that horizontal gene transfer has allowed some microbes to “cheat” and take up iron-bound siderophores without expending energy to produce the siderophores themselves.

Using the cheese rind microbiome as a model system, it becomes possible to explore how horizontal gene transfer works in more detail than before. A better understanding of this process can then be applied to other important microbiomes, including those where genes conferring antibiotic resistance are commonly exchanged.

DOI:http://dx.doi.org/10.7554/eLife.22144.002

Use of a metagenetic approach to monitor the bacterial microbiota of "Tomme d'Orchies" cheese during the ripening process

The study of microbial ecosystems in artisanal foodstuffs is important to complete in order to unveil its diversity. The number of studies performed on dairy products has increased during the last decade, particularly those performed on milk and cheese derivative products. In this work, we investigated the bacterial content of "Tomme d'Orchies" cheese, an artisanal pressed and uncooked French cheese. To this end, a metagenetic analysis, using Illumina technology, was utilized on samples taken from the surface and core of the cheese at 0, 1, 3, 14 and 21days of ripening process. In addition to the classical microbiota found in cheese, various strains likely from environmental origin were identified. A large difference between the surface and the core content was observed within samples withdrawn during the ripening process. The main species encountered in the core of the cheese were Lactococcus spp. and Streptococcus spp., with an inversion of this ratio during the ripening process. Less than 2.5% of the whole population was composed of strains issued from environmental origin, as Lactobacillales, Corynebacterium and Brevibacterium. In the core, about 85% of the microbiota was attributed to the starters used for the cheese making. In turn, the microbiota of the surface contained less than 30% of these starters and interestingly displayed more diversity. The predominant genus was Corynebacterium sp., likely originating from the environment. The less abundant microbiota of the surface was composed of Bifidobacteria, Brevibacterium and Micrococcales. To summarize, the "Tomme d'Orchies" cheese displayed a high diversity of bacterial species, especially on the surface, and this diversity is assumed to arise from the production environment and subsequent ripening process.

Biotic Interactions Shape the Ecological Distributions of Staphylococcus Species

Many metagenomic sequencing studies have observed the presence of closely related bacterial species or genotypes in the same microbiome. Previous attempts to explain these patterns of microdiversity have focused on the abiotic environment, but few have considered how biotic interactions could drive patterns of microbiome diversity. We dissected the patterns, processes, and mechanisms shaping the ecological distributions of three closely related Staphylococcus species in cheese rind biofilms. Paradoxically, the most abundant species (S. equorum) is the slowest colonizer and weakest competitor based on growth and competition assays in the laboratory. Through in vitro community reconstructions, we determined that biotic interactions with neighboring fungi help resolve this paradox. Species-specific stimulation of the poor competitor by fungi of the genus Scopulariopsis allows S. equorum to dominate communities in vitro as it does in situ. Results of comparative genomic and transcriptomic experiments indicate that iron utilization pathways, including a homolog of the S. aureus staphyloferrin B siderophore operon pathway, are potential molecular mechanisms underlying Staphylococcus-Scopulariopsis interactions. Our integrated approach demonstrates that fungi can structure the ecological distributions of closely related bacterial species, and the data highlight the importance of bacterium-fungus interactions in attempts to design and manipulate microbiomes.

Decades of culture-based studies and more recent metagenomic studies have demonstrated that bacterial species in agriculture, medicine, industry, and nature are unevenly distributed across time and space. The ecological processes and molecular mechanisms that shape these distributions are not well understood because it is challenging to connect in situ patterns of diversity with mechanistic in vitro studies in the laboratory. Using tractable cheese rind biofilms and a focus on coagulase-negative Staphylococcus (CNS) species, we demonstrate that fungi can mediate the ecological distributions of closely related bacterial species. One of the Staphylococcus species studied, S. saprophyticus, is a common cause of urinary tract infections. By identifying processes that control the abundance of undesirable CNS species, cheese producers will have more precise control on the safety and quality of their products. More generally, Staphylococcus species frequently co-occur with fungi in mammalian microbiomes, and similar bacterium-fungus interactions may structure bacterial diversity in these systems.

Microbiological quality of smear-ripened cheeses stored in different temperature regimes

The purpose of this work was to study smear-ripened cheese, especially its microbiological quality. Samples were stored in different temperature conditions. The first group (A) of samples was stored in a refrigerator at 6 °C. The second group (B) of samples was stored at 6 °C for 21 days, next at -18 °C for 7 days and 7 days at 6 °C. The third group (C) of samples was stored  at  6 °C before the date of minimum durability, next 7 weeks at -18 °C and after that at 6 °C for 7 days. I have observed  lactic  acid  bacteria, Brevibacterium  linens, coliforms,  psychrotrophic  organisms, Escherichia  coli, moulds  and yeast. The number of Brevibacterium linens was higher (p <0.05) at cool ing/freezing for 1 week (log CFU.g-1) than cooling (log CFU.g-1) and cooling/freezing for 7 weeks after 35 days. A higher (p <0.05) number of psychrotrophic microorganisms was recorded at the end of the monitoring in samples stored in a refrigerator (A/49) in comparison with cheeses stored using cooling/freezing  regime  for  one  week  (B/49). Among  the  samples  stored  at  these  temperature  regimes,  there  were  no statistically significant differences (p >0.05) in the numbers of LAB, coliforms, E. coli, moulds and yeast, neither at the end of the DMD nor at storage for 49 or 91 days, respectively.

Adaptive Horizontal Gene Transfers between Multiple Cheese-Associated Fungi

Domestication is an excellent model for studies of adaptation because it involves recent and strong selection on a few, identified traits [1–5]. Few studies have focused on the domestication of fungi, with notable exceptions [6–11], despite their importance to bioindustry [12] and to a general understanding of adaptation in eukaryotes [5]. Penicillium fungi are ubiquitous molds among which two distantly related species have been independently selected for cheese making—P. roqueforti for blue cheeses like Roquefort and P. camemberti for soft cheeses like Camembert. The selected traits include morphology, aromatic profile, lipolytic and proteolytic activities, and ability to grow at low temperatures, in a matrix containing bacterial and fungal competitors [13–15]. By comparing the genomes of ten Penicillium species, we show that adaptation to cheese was associated with multiple recent horizontal transfers of large genomic regions carrying crucial metabolic genes. We identified seven horizontally transferred regions (HTRs) spanning more than 10 kb each, flanked by specific transposable elements, and displaying nearly 100% identity between distant Penicillium species. Two HTRs carried genes with functions involved in the utilization of cheese nutrients or competition and were found nearly identical in multiple strains and species of cheese-associated Penicillium fungi, indicating recent selective sweeps; they were experimentally associated with faster growth and greater competitiveness on cheese and contained genes highly expressed in the early stage of cheese maturation. These findings have industrial and food safety implications and improve our understanding of the processes of adaptation to rapid environmental changes.

•

New HTRs are found in cheese fungi

•

HTRs are flanked by specific transposable elements

•

HTRs have spread in cheese-associated fungi through recent selective sweeps

•

Experiments link two HTRs to growth and competitive advantages on cheese