Fig. 5

Covariance between microbiome composition and wMel abundance. A PCoA analysis of microbiome variation identified a primary axis of variation (PC1: 58.73%) that corresponded with the eigenvector for wMel which was the largest observed vector. Other large eigenvectors representing highly abundant taxa were largely orthogonal to Axis 1. B Axis 1 loadings were significantly correlated with wMel abundance across individuals but not other microbiome species. C PCoA analysis of microbiome variation identified a secondary axis of variation (PC2: 10.19%) which corresponded with the eigenvectors for B. subtilis, M. haemolytica, and to a lesser extent wMel. D Axis 2 loadings were significantly correlated with B. subtilis and M. haemolytica abundance and anticorrelated with wMel abundance across individuals