Cultivable microbial diversity, peptide profiles, and bio-functional properties in Parmigiano Reggiano cheese

Introduction Lactic acid bacteria (LAB) communities shape the sensorial and functional properties of artisanal hard-cooked and long-ripened cheeses made with raw bovine milk like Parmigiano Reggiano (PR) cheese. While patterns of microbial evolution have been well studied in PR cheese, there is a lack of information about how this microbial diversity affects the metabolic and functional properties of PR cheese. Methods To fill this information gap, we characterized the cultivable fraction of natural whey starter (NWS) and PR cheeses at different ripening times, both at the species and strain level, and investigated the possible correlation between microbial composition and the evolution of peptide profiles over cheese ripening. Results and discussion The results showed that NWS was a complex community of several biotypes belonging to a few species, namely, Streptococcus thermophilus , Lactobacillus helveticus , and Lactobacillus delbrueckii subsp. lactis . A new species-specific PCR assay was successful in discriminating the cheese-associated species Lacticaseibacillus casei , Lacticaseibacillus paracasei , Lacticaseibacillus rhamnosus , and Lacticaseibacillus zeae . Based on the resolved patterns of species and biotype distribution, Lcb. paracasei and Lcb. zeae were most frequently isolated after 24 and 30 months of ripening, while the number of biotypes was inversely related to the ripening time. Peptidomics analysis revealed more than 520 peptides in cheese samples. To the best of our knowledge, this is the most comprehensive survey of peptides in PR cheese. Most of them were from β-caseins, which represent the best substrate for LAB cell-envelope proteases. The abundance of peptides from β-casein 38–88 region continuously increased during ripening. Remarkably, this region contains precursors for the anti-hypertensive lactotripeptides VPP and IPP, as well as for β-casomorphins. We found that the ripening time strongly affects bioactive peptide profiles and that the occurrence of Lcb. zeae species is positively linked to the incidence of eight anti-hypertensive peptides. This result highlighted how the presence of specific LAB species is likely a pivotal factor in determining PR functional properties.