Niels O. Verhulst

Background The skin is inhabited by a variety of micro‐organisms, with bacteria representing the predominant taxon of the skin microbiome. In sheep, the skin bacterial community of healthy animals has been addressed in few studies, only with culture‐based methods or sequencing of cloned amplicons. Objectives The objectives of this study were to determine the sheep skin bacterial community composition by using metabarcoding for a detailed characterisation and to determine the effect of body part, breed and environment. Materials and Methods Overall, 267 samples were taken from 89 adult female sheep, belonging to three different breeds and kept on nine different farms in Switzerland. From every individual, one sample each was taken from belly, left ear and left leg and metabarcoding of the 16S rRNA V3–V4 hypervariable region was performed. Results The main phyla identified were Actinobacteriota, Firmicutes, Proteobacteria and Bacteriodota. The alpha diversity as determined by Shannon's diversity index was significantly different between sheep from different farms. Beta diversity analysis by principal coordinate analysis (PCoA) showed clustering of the samples by farm and body site, while breed had only a marginal influence. A sparse partial least squares discriminant analysis (sPLS‐DA) revealed seven main groups of operational taxonomic units (OTUs) of which groups of OTUs were specific for some farms. Conclusions and Clinical Relevance These findings indicate that environment has a larger influence on skin microbial variability than breed, although the sampled breeds, the most abundant ones in Switzerland, are phenotypically similar. Future studies on the sheep skin microbiome may lead to novel insights in skin diseases and prevention.