Comprehensive Virulence Gene Profiling of Bovine Non-aureus Staphylococci Based on Whole-Genome Sequencing Data

Non-aureus staphylococci (NAS) are the most frequently isolated pathogens from intramammary infection (IMI) in dairy cattle. Virulence factors (VFs) and mechanisms by which NAS cause IMI are not fully known. Herein, we analyzed the distribution of 191 VFs in 441 genomes of 25 NAS species, after classifying VFs into functional categories: adherence (n = 28), exoenzymes (n = 21), immune evasion (n = 20), iron metabolism (n = 29), and toxins (n = 93). In addition to establishing VF gene profiles, associations of VF genes between and among functional categories were computed, revealing distinctive patterns of association among VFs for various NAS species. Associations were also computed for low, medium, and high somatic cell count (SCC) and clinical mastitis (CM) isolates, demonstrating distinctive patterns of associations for low SCC and CM isolates, but no differences between high SCC and CM isolates. To determine whether VF distributions had any association with SCC or CM, various clustering approaches, including complete linkages, Ward clustering, and t-distributed stochastic neighbor embedding, were applied. However, no clustering of isolates representing low SCC, medium SCC, or high SCC or CM was identified. Regression analysis to test for associations with individual VF functional categories demonstrated that each additional toxin and host immune evasion gene increased the odds of having high SCC or CM, although an overall increase in the number of VFs was not associated with increased SCC or occurrence of CM. In conclusion, we established comprehensive VF gene profiling, determined VF gene distributions and associations, calculated pathogenic potentials of all NAS species, and detected no clear link between VF genes and mastitis. IMPORTANCE Non-aureus staphylococci (NAS) are the most frequently isolated pathogens from milk in dairy cattle worldwide. The virulence factors (VFs) and mechanisms by which these bacteria cause udder infection are not fully known. We determined the distribution and associations of 191 VFs in 25 NAS species and investigated the relationship between VFs and disease. Although the overall number of VFs was not associated with disease severity, increasing numbers of toxin and host immune evasion genes specifically were associated with more severe disease outcomes. These findings suggest that the development of disease and the interactions of VFs with the host are complex and determined by the interplay of genes rather than just the presence of virulence genes. Together, our results provide foundational genetic knowledge to other researchers to design and conduct further experiments, focusing on understanding the synergy between VFs and roles of individual NAS species in IMI and characterizing species-specific effects on udder health.