Changes in bovine milk bacterial microbiome from healthy and subclinical mastitis affected animals of the Girolando, Gyr, Guzera, and Holstein breeds.

Raw milk samples were collected from 200 dairy cows belonging to Girolando 1/2, Gyr, Guzera, and Holstein breeds, and the bacterial diversity was explored using 16S rRNA amplicon sequencing. SCC analysis showed that 69 animals were classified as affected with subclinical mastitis. The milk bacterial microbiome was dominated by Firmicutes, Proteobacteria, and Actinobacteria, with an increase of Firmicutes in animals with subclinical mastitis and Proteobacteria in healthy animals. At the family and genus level, the milk bacterial microbiome was dominated by Staphylococcus, Acinetobacter, Pseudomonas, members of the family Enterobacteriaceae, Lactococcus, Aerococcus, members of the family Rhizobiaceae, Anaerobacillus, Streptococcus, members of the family Intrasporangiaceae, members of the family Planococcaceae, Corynebacterium, Nocardioides, and Chryseobacterium. Significant differences in alpha and beta diversity analysis suggest an effect of udder health status and breed on the composition of raw bovine milk microbiota. LEfSe analysis showed 45 and 51 discriminative taxonomic biomarkers associated with udder health status and with one of the four breeds respectively, suggesting an effect of subclinical mastitis and breed on the microbiota of milk in cattle.