Genome-Wide Association Study of Nucleotide Variants Associated with Resistance to Nine Antimicrobials in Mycoplasma bovis

Antimicrobial resistance (AMR) studies of Mycoplasma bovis have generally focused on specific loci versus using a genome-wide association study (GWAS) approach. A GWAS approach, using two different models, was applied to 194 Mycoplasma bovis genomes. Both a fixed effects linear model (FEM) and a linear mixed model (LMM) identified associations between nucleotide variants (NVs) and antimicrobial susceptibility testing (AST) phenotypes. The AMR phenotypes represented fluoroquinolones, tetracyclines, phenicols, and macrolides. Both models identified known and novel NVs associated (Bonferroni adjusted p < 0.05) with AMR. Fluoroquinolone resistance was associated with multiple NVs, including previously identified mutations in gyrA and parC. NVs in the 30S ribosomal protein 16S were associated with tetracycline resistance, whereas NVs in 5S rRNA, 23S rRNA, and 50S ribosomal proteins were associated with phenicol and macrolide resistance. For all antimicrobial classes, resistance was associated with NVs in genes coding for ABC transporters and other membrane proteins, tRNA-ligases, peptidases, and transposases, suggesting a NV-based multifactorial model of AMR in M. bovis. This study was the largest collection of North American M. bovis isolates used with a GWAS for the sole purpose of identifying novel and non-antimicrobial-target NVs associated with AMR.