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

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.

Relationship between the Dynamics of Gross Composition, Free Fatty Acids and Biogenic Amines, and Microbial Shifts during the Ripening of Raw Ewe Milk-Derived Idiazabal Cheese

This study reports for the first time the relationship between bacterial succession, characterized by high-throughput sequencing (sequencing of V3-V4 16S rRNA regions), and the evolution of gross composition, free fatty acids (FFAs) and biogenic amines (BAs) during cheese ripening. Specifically, Idiazabal PDO cheese, a raw ewe milk-derived semi-hard o hard cheese, was analysed. Altogether, 8 gross parameters were monitored (pH, dry matter, protein, fat, Ca, Mg, P and NaCl) and 21 FFAs and 8 BAs were detected. The ripening time influenced the concentration of most physico-chemical parameters, whereas the producer mainly affected the gross composition and FFAs. Through an O2PLS approach, the non-starter lactic acid bacteria Lactobacillus, Enterococcus and Streptococcus were reported as positively related to the evolution of gross composition and FFAs release, while only Lactobacillus was positively related to BAs production. Several environmental or non-desirable bacteria showed negative correlations, which could indicate the negative impact of gross composition on their growth, the antimicrobial effect of FFAs and/or the metabolic use of FFAs by these genera, and their ability to degrade BAs. Nonetheless, Obesumbacterium and Chromohalobacter were positively associated with the synthesis of FFAs and BAs, respectively. This research work provides novel information that may contribute to the understanding of possible functional relationships between bacterial communities and the evolution of several cheese quality and safety parameters.

Bioinformatic Prediction and Characterization of Proteins in Porphyra dentata by Shotgun Proteomics

Porphyra dentata is an edible red seaweed with high nutritional value. It is widely cultivated and consumed in East Asia and has vast economic benefits. Studies have found that P. dentata is rich in bioactive substances and is a potential natural resource. In this study, label-free shotgun proteomics was first applied to identify and characterize different harvest proteins in P. dentata. A total of 13,046 different peptides were identified and 419 co-expression target proteins were characterized. Bioinformatics was used to study protein characteristics, functional expression, and interaction of two important functional annotations, amino acid, and carbohydrate metabolism. Potential bioactive peptides, protein structure, and potential ligand conformations were predicted, and the results suggest that bioactive peptides may be utilized as high-quality active fermentation substances and potential targets for drug production. Our research integrated the global protein database, the first time bioinformatic analysis of the P. dentata proteome during different harvest periods, improves the information database construction and provides a framework for future research based on a comprehensive understanding.

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.

Evaluation of bacterial diversity of traditional cheese in Tarbagatay Prefecture, China, and its correlation with cheese quality

In Xinjiang, China, traditional handmade cheese is made from fresh milk under natural environmental conditions and is a common fermented dairy product in the region. Due to differences in production methods between regions, the research conducted on the bacterial diversity of traditional handmade cheese is not comprehensive. Hence, little is known about the relationship between bacteria and cheese quality. Therefore, in this study, cheese samples from Tarbagatay Prefecture, Xinjiang, were chosen for investigation. The bacteria in 17 cheese samples were analyzed by sequencing 16S rRNA using Illumina MiSeq technology. The results showed that there were two dominant bacterial phyla and six dominant bacterial genera in the cheeses. Of these, Lactobacillus and Lactococcus displayed the most significant positive correlation with cheese quality. This study provides data to support the improvement of traditional cheese quality via microbial diversity and lays a foundation for the industrialization of traditional cheese production.

Evaluating the Authenticity of the Raw-Milk Cheese Fontina (PDO) with Respect to Similar Cheeses

The implementation of quality assurance schemes for the assessment of PDO food authenticity is an issue involving manufacturers, traders, retailers and consumers. In this respect, reliable analytical methods are needed to integrate paper-trailing information. The feasibility of distinguishing the Italian Fontina PDO cheese from the generic Fontal cheese was preliminarily evaluated on a set of commercial samples by measuring selected parameters (pH, alkaline phosphatase activity, content of copper, volatiles, extent of proteolysis) related to the different manufacturing processes. The relative profile of free amino acids proved to be a promising tool. A new set of 41 samples of Fontina PDO cheese was collected at representative dairies within the recognized production area and analyzed for free amino acids. A chemometric model of Fontina PDO cheese was built based on the mean content and standard deviation of 15 free amino acids. On this basis, all of the PDO samples were correctly identified, whereas all of the Fontal cheeses were recognized as different 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.

Psychroflexus aurantiacus sp. nov., isolated from soil in the Yellow River Delta wetlands

A Gram-stain-negative, strictly aerobic, non-motile, orange-coloured bacterium, designated YR1-1^T, was isolated from a soil sample collected from the Yellow River Delta wetlands (PR China). Growth was observed at a salinity of 1.0-15.0 % NaCl, 4-45 °C and pH 6.0-9.0. The results of phylogenetic analysis based on the 16S rRNA gene sequences indicated that YR1-1^T represented a member of the genus Psychroflexus, with the highest sequence similarity to Psychroflexus sediminis YIM-C238^T (97.9 %), followed by Psychroflexus aestuariivivens (97.1 %) and Psychroflexus torquis (96.4 %). The average nucleotide identity and digital DNA-DNA hybridization values between YR1-1^T and other closely related type strains of species of the genus Psychroflexus were 68.7-86.3% and 17.8-30.9 %. The genome of the strain was 2 899 374 bp in length with 39.8 % DNA G+C content. The predominant fatty acids (>10 %) were iso-C_15 : 0 and anteiso-C_15 : 0. The major respiratory quinone was menaquinone-6 (MK-6) and the major polar lipids were phosphatidylethanolamine, phospholipid, diphosphatidylglycerol, two unidentified aminolipids and four unidentified lipids. The combined genotypic and phenotypic data indicate that YR1-1^T represents a novel species within the genus Psychroflexus, for which the name Psychroflexus aurantiacus sp. nov., is proposed. The type strain is YR1-1^T (=KCTC 72794^T=CGMCC 1.17458^T).

The transcriptome of Listeria monocytogenes during co-cultivation with cheese rind bacteria suggests adaptation by induction of ethanolamine and 1,2-propanediol catabolism pathway genes

The survival of Listeria (L.) monocytogenes in foods and food production environments (FPE) is dependent on several genes that increase tolerance to stressors; this includes competing with intrinsic bacteria. We aimed to uncover genes that are differentially expressed (DE) in L. monocytogenes sequence type (ST) 121 strain 6179 when co-cultured with cheese rind bacteria. L. monocytogenes was cultivated in broth or on plates with either a Psychrobacter or Brevibacterium isolate from cheese rinds. RNA was extracted from co-cultures in broth after two or 12 hours and from plates after 24 and 72 hours. Broth co-cultivations with Brevibacterium or Psychrobacter yielded up to 392 and 601 DE genes, while plate co-cultivations significantly affected the expression of up to 190 and 485 L. monocytogenes genes, respectively. Notably, the transcription of virulence genes encoding the Listeria adhesion protein and Listeriolysin O were induced during plate and broth co-cultivations. The expression of several systems under the control of the global stress gene regulator, σB, increased during co-cultivation. A cobalamin-dependent gene cluster, responsible for the catabolism of ethanolamine and 1,2-propanediol, was upregulated in both broth and plate co-cultures conditions. Finally, a small non-coding (nc)RNA, Rli47, was induced after 72 hours of co-cultivation on plates and accounted for 50-90% of the total reads mapped to L. monocytogenes. A recent study has shown that Rli47 may contribute to L. monocytogenes stress survival by slowing growth during stress conditions through the suppression of branch-chained amino acid biosynthesis. We hypothesize that Rli47 may have an impactful role in the response of L. monocytogenes to co-cultivation by regulating a complex network of metabolic and virulence mechanisms.

Psychroflexus maritimus sp. nov., isolated from coastal sediment

A Gram-stain-negative, strictly aerobic, non-flagellated, gliding, rod-shaped bacterium, designated C1^T was isolated from the coastal sediment of Xiaoshi Island, Weihai, China. Phylogenetic analysis of the 16S rRNA gene and the genome sequence of the newly isolated strain revealed that it belonged to the genus Psychroflexus within the family Flavobacteriaceae. The 16S rRNA gene sequence similarities between strain C1^T and the type strains of Psychroflexus members ranged from 91.6 to 94.9%. The AAI, POCP, ANI and dDDH values between strain C1^T and Psychroflexus torquis ATCC 700755^T were 67.0%, 53.1%, 70.5% and 20.2%, respectively. The DNA G+C content was 34.0 mol%. Strain C1^T grows optimally at 28-30 ℃, at pH 7.5-8.0 and with 3.0-5.0% (w/v) NaCl, and its colonies were orange-colored, convex and circular on the MA plate. Positive for hydrolysis of Tween 20 and catalase activities. The dominant respiratory quinone was menaquinone-6, and the major fatty acids were iso-C_17:0 3-OH and iso-C_15:0. The polar lipids of strain C1^T consisted of phosphatidylethanolamine and three unidentified lipids. On the basis of phylogenetic, phenotypic and chemotaxonomic data, it is considered that strain C1^T represents a novel species within the genus Psychroflexus, for which, the name Psychroflexus maritimus sp. nov., is proposed. The type strain is C1^T (= MCCC 1H00415^T = KCTC 72796^T).

Microbial communities of a variety of cheeses and comparison between core and rind region of cheeses

Understanding the microbial community of cheese is important in the dairy industry, as the microbiota contributes to the safety, quality, and physicochemical and sensory properties of cheese. In this study, the microbial compositions of different cheeses (Cheddar, provolone, and Swiss cheese) and cheese locations (core, rind, and mixed) collected from the Arbuthnot Dairy Center at Oregon State University were analyzed using 16S rRNA gene amplicon sequencing with the Illumina MiSeq platform (Illumina, San Diego, CA). A total of 225 operational taxonomic units were identified from the 4,675,187 sequencing reads generated. Streptococcus was observed to be the most abundant organism in provolone (72 to 85%) and Swiss (60 to 67%), whereas Lactococcus spp. were found to dominate Cheddar cheese (27 to 76%). Species richness varied significantly by cheese. According to alpha diversity analysis, porter-soaked Cheddar cheese exhibited the highest microbial richness, whereas smoked provolone cheese showed the lowest. Rind regions of each cheese changed color through smoking and soaking for the beverage process. In addition, the microbial diversity of the rind region was higher than the core region because smoking and soaking processes directly contacted the rind region of each cheese. The microbial communities of the samples clustered by cheese, indicated that, within a given type of cheese, microbial compositions were very similar. Moreover, 34 operational taxonomic units were identified as biomarkers for different types of cheese through the linear discriminant analysis effect size method. Last, both carbohydrate and AA metabolites comprised more than 40% of the total functional annotated genes from 9 varieties of cheese samples. This study provides insight into the microbial composition of different types of cheese, as well as various locations within a cheese, which is applicable to its safety and sensory quality.

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.

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.

The role of the surface smear microbiome in the development of defective smear on surface-ripened red-smear cheese

The complex smear microbiota colonizing the surface of red-smear cheese fundamentally impacts the ripening process, appearance and shelf life of cheese. To decipher the prokaryotic composition of the cheese smear microbiome, the surface of a semi-hard surface ripened cheese was studied post-ripening by culture-based and culture-independent molecular approaches. The aim was to detect potential bacterial alterations in the composition of the cheese smear microbiota resulting from cheese storage in vacuum film-prepackaging, which is often accompanied by the development of a surface smear defect. Next-generation sequencing of amplified 16S rRNA gene fragments revealed an unexpected high diversity of a total of 132 different genera from the domains Bacteria and Archaea on the cheese surface. Beside typical smear organisms, our study revealed the presence of several microorganisms so far not associated with cheese, but related to milk, farm and cheese dairy environments. A 16S ribosomal RNA based analysis from total RNA identified the major metabolically active populations in the cheese surface smear as Actinobacteria of the genera Corynebacterium, Brevibacterium, Brachybacterium and Agrococcus. Comparison of data on a higher phylogenetic level revealed distinct differences in the composition of the cheese smear microbiome from the different samples. While the proportions of Proteobacteria and Bacteroidetes were increased in the smear of prepacked samples and in particular in defective smear, staphylococci showed an opposite trend and turned out to be strongly decreased in defective smear. In conclusion, next-generation sequencing of amplified 16S rRNA genes and 16S rRNA from total RNA extracts provided a much deeper insight into the bacterial composition of the cheese smear microbiota. The observed shifts in the microbial composition of samples from defect surface smear suggest that certain members of the Proteobacteria contribute to the observed negative organoleptic properties of the surface smear of cheese after prepacking in plastic foil.