Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.

Single-nucleotide polymorphism in two representative multidrug-resistant Mycobacterium bovis isolates collected from patients in a Spanish hospital harboring a human infection outbreak

Mycobacterium bovis is the etiological agent of tuberculosis in domestic and wild animals. Its involvement as a human pathogen has been highlighted again with the recent descriptions of transmission through dairy products (18), reactivation or primary infection in human immunodeficiency virus-infected patients (5), and association with meat industry workers, animal keepers, or hunters (3). Strains resistant to antituberculous drugs (M. bovis is naturally resistant to pyrazinamide) pose an additional risk (2). Several studies have demonstrated that mutations in target genes are associated with resistance to antituberculous drugs (4, 7, 10, 11, 16). However, most of them have been developed in Mycobacterium tuberculosis strains and limited data are available regarding M. bovis isolates. The aim of this study was to characterize by sequencing the main genes involved in antibiotic resistance in two multidrug-resistant (MDR) M. bovis isolates in a human outbreak detected in a hospital in Madrid that subsequently spread to several countries (5, 6, 15). The isolates were resistant to 11 drugs, but only their rpoB and katG genes have been analyzed so far (1, 14). We studied the first (93/R1) and last (95/R4) M. bovis isolates of this nosocomial outbreak, characterized by spoligotyping as SB0426 (hexacode 63-5F-5E-7F-FF-60 in the database at www.mbovis.org) (1, 13). Several genes involved in resistance to isoniazid (katG, ahpC, inhA, and the oxyR-ahpC intergenic region), rifampin (rpoB), streptomycin (rrs, rpsL), ethambutol (embB), and quinolones (gyrA) were studied. These genes, or fragments of genes, were amplified and sequenced as previously described (12). The sequence analysis revealed polymorphisms in five (ahpC, rpoB, rpsL, embB, and gyrA) out of nine analyzed genes (Table 1). Nucleotide substitutions in four genes cause a change in the encoded amino acid. Two additional synonymous mutations in ahpC and rpsL differentiated the first and last isolates from the outbreak.