This study aims to genomically elucidate six isolates of rapidly growing non-tuberculous mycobacteria (RGM) derived from Siamese fighting fish (Betta splendens). These isolates had previously undergone phenotypic and biochemical characterization, antibiotic susceptibility testing, and in vivo virulence assessment. Initial DNA barcoding using the 16S rRNA sequence assigned these six isolates to five different species, namely Mycobacterium chelonae (BN1983), M. cosmeticum (BN1984 and N041), M. farcinogenes (SNSK5), M. mucogenicum (BN1956), and M. senegalense (BN1985). However, the identification relied solely on the highest percent identity of the 16S rRNA gene, raising concerns about the taxonomic ambiguity of these species. Comprehensive whole genome sequencing (WGS) and extended genomic comparisons using multilocus sequence typing (MLST), average nucleotide identity (ANI), and digital DNA-DNA hybridization (dDDH) led to the reclassification of BN1985 and SNSK5 as M. conceptionense while confirming BN1983 as M. chelonae and BN1984 and N041 as M. cosmeticum. Notably, the analysis of the BN1956 isolate revealed a potential new species that is proposed here as M. mucogenicum subsp. phocaicum sp. nov. Common genes encoding "mycobacterial" virulence proteins, such as PE and PPE family proteins, MCE, and YrbE proteins, were detected in all six isolates. Two species, namely M. chelonae and M. cosmeticum, appear to have horizontally acquired T6SS-II (clpB), catalase (katA), GroEL (groel), and capsule (rmlb) from distantly related environmental bacteria such as Klebsiella sp., Neisseria sp., Clostridium sp., and Streptococcus sp. This study provides the first draft genome sequence of RGM isolates currently circulating in B. splendens and underscores the necessity of WGS for the identification and classification of mycobacterial species.
Keywords: Betta splendens; average nucleotide identity; genomic sequencing; non-tuberculous mycobacteria; whole genome.