A comprehensive, large-scale analysis of “terroir” cheese and milk microbiota reveals profiles strongly shaped by both geographical and human factors

An exhaustive analysis was performed on more than 2000 microbiotas from French Protected Designation of Origin (PDO) cheeses covering most cheese families produced throughout the world. Thanks to a complete and accurate set of associated metadata, we have carried out a deep analysis of the ecological drivers of microbial communities in milk and “terroir” cheeses. We show that bacterial and fungal microbiota from milk differed significantly across dairy species while sharing a core microbiome consisting of four microbial species. By contrast, no microbial species was detected in all ripened cheese samples. Our network analysis suggested that the cheese microbiota was organized into independent network modules. These network modules comprised mainly species with an overall relative abundance lower than 1%, showing that the most abundant species were not those with the most interactions. Species assemblages differed depending on human drivers, dairy species, and geographical area, thus demonstrating the contribution of regional know-how to shaping the cheese microbiota. Finally, an extensive analysis at the milk-to-cheese batch level showed that a high proportion of cheese taxa were derived from milk under the influence of the dairy species and PDO.