Tuberculosis (TB) remains one of the most common infectious diseases caused by Mycobacterium tuberculosis complex (MTBC). To panoramically analyze MTBC's genomic methylation, we completed the genomes of 12 MTBC strains (Mycobacterium bovis; M. bovis BCG; M. microti; M. africanum; M. tuberculosis H37Rv; H37Ra; and 6 M. tuberculosis clinical isolates) belonging to different lineages and characterized their methylomes using single-molecule real-time (SMRT) technology. We identified three ^m6A sequence motifs and their corresponding methyltransferase (MTase) genes, including the reported mamA, hsdM and a newly discovered mamB. We also experimentally verified the methylated motifs and functions of HsdM and MamB. Our analysis indicated the MTase activities varied between 12 strains due to mutations/deletions. Furthermore, through measuring ‘the methylated-motif-site ratio’ and ‘the methylated-read ratio’, we explored the methylation status of each modified site and sequence-read to obtain the ‘precision methylome’ of the MTBC strains, which enabled intricate analysis of MTase activity at whole-genome scale. Most unmodified sites overlapped with transcription-factor binding-regions, which might protect these sites from methylation. Overall, our findings show enormous potential for the SMRT platform to investigate the precise character of methylome, and significantly enhance our understanding of the function of DNA MTase.

Tuberculosis (TB) has surpassed HIV as the leading infectious disease killer worldwide since 2014. The main pathogen, Mycobacterium tuberculosis (Mtb), contains ~4,000 genes that account for ~90% of the genome. However, it is still unclear which of these genes are primary/secondary, which are responsible for generality/individuality, and which interconvert during evolution. Here we utilized a pan-genomic analysis of 36 Mtb genomes to address these questions. We identified 3,679 Mtb core (i.e., primary) genes, determining their phenotypic generality (e.g., virulence, slow growth, dormancy). We also observed 1,122 dispensable and 964 strain-specific secondary genes, reflecting partially shared and lineage-/strain-specific individualities. Among which, five L2 lineage-specific genes might be related to the increased virulence of the L2 lineage. Notably, we discovered 28 Mtb “Super Core Genes” (SCGs: more than a copy in at least 90% strains), which might be of increased importance, and reflected the “super phenotype generality.” Most SCGs encode PE/PPE, virulence factors, antigens, and transposases, and have been verified as playing crucial roles in Mtb pathogenicity. Further investigation of the 28 SCGs demonstrated the interconversion among SCGs, single-copy core, dispensable, and strain-specific genes through copy number variations (CNVs) during evolution; different mutations on different copies highlight the delicate adaptive-evolution regulation amongst Mtb lineages. This reflects that the importance of genes varied through CNVs, which might be driven by selective pressure from environment/host-adaptation. In addition, compared with Mycobacterium bovis (Mbo), Mtb possesses 48 specific single core genes that partially reflect the differences between Mtb and Mbo individuality.

Bone and Joint tuberculosis (BJTB) constitutes about 10% of total extra-pulmonary TB cases. Since the BJTB is a paucibacillary condition, there has been no systematic study on the bacterial characterization, especially the epidemiological feature. Here we collected the mycobacterial clinical isolates, analyzed the clinical features and the bacteriological characteristics from 113 BJTB cases reported in China. The mean age of the cases was 40.33 years while most of the patients fell into the 20–29 year age group; local pain was the most common onset symptom of BJTB cases; mean time from symptom onset to BJTB diagnosis was 13.16 months. 31 isolates were defined as drug resistant, including 15 multidrug resistant (MDR) and 2 extensively drug resistant (XDR) isolates according to the drug susceptibility test outcomes; after spoligotyping, 87.6% (99/113) isolates were categorized as Beijing family. In contrast to the isolates from pulmonary tuberculosis patients, here the MIRU-VNTR assay did not find anything significant. A prolonged time span for BJTB diagnosis highlights the requirement of paying further attention to BJTB infection in China. This study provides essential insights into the demographic and microbial characteristics of BJTB cases in China.

p -Aminosalicylic acid (PAS) is an important compound for treating multidrug-resistant tuberculosis (TB). Previous studies showed that thyA mutations are often related to PAS resistance in clinical isolates. We performed a systematic analysis of isolate genotypes and detected mutations in three folate pathway genes ( folC , thyA , and ribD ) in 61.1% (127/208) of PAS-resistant isolates, including 11 double mutants. This result expands our knowledge about the distribution and frequency of mutations related to PAS resistance in mycobacterial clinical isolates.