Identification of Antimicrobial Resistance-Associated Genes through Whole Genome Sequencing of Mycoplasma bovis Isolates with Different Antimicrobial Resistances

Antimicrobial resistance (AMR) in Mycoplasma bovis has been previously associated with topoisomerase and ribosomal gene mutations rather than specific resistance-conferring genes. Using whole genome sequencing (WGS) to identify potential new AMR mechanisms for M. bovis, it was found that a 2019 clinical isolate with high MIC (2019-043682) for fluoroquinolones, macrolides, lincosamides, pleuromutilins and tetracyclines had a new core genome multilocus sequencing (cgMLST) type (ST10-like) and 91% sequence similarity to the published genome of M. bovis PG45. Closely related to PG45, a 1982 isolate (1982-M6152) shared the same cgMLST type (ST17), 97.2% sequence similarity and low MIC results. Known and potential AMR- associated genetic events were identified through multiple sequence alignment of the three genomes. Isolate 2019-043682 had 507 genes with non-synonymous mutations (NSMs) and 67 genes disrupted. Isolate 1982-M6152 had 81 NSMs and 20 disruptions. Using functional roles and known mechanisms of antimicrobials, a 55 gene subset was assessed for AMR potential. Seventeen were previously identified from other bacteria as sites of AMR mutation, 38 shared similar functions to them, and 11 contained gene-disrupting mutations. This study indicated that M. bovis may obtain high AMR characteristics by mutating or disrupting other functional genes, in addition to topoisomerases and ribosomal genes.