This study assembles DNA adenine methylomes for 93 Mycobacterium tuberculosis complex (MTBC) isolates from seven lineages paired with fully-annotated, finished, de novo assembled genomes. Integrative analysis yielded four key results. First, methyltransferase allele-methylome mapping corrected methyltransferase variant effects previously obscured by reference-based variant calling. Second, heterogeneity analysis of partially active methyltransferase alleles revealed that intracellular stochastic methylation generates a mosaic of methylomes within isogenic cultures, which we formalize as ‘intercellular mosaic methylation’ (IMM). Mutation-driven IMM was nearly ubiquitous in the globally prominent Beijing sublineage. Third, promoter methylation is widespread and associated with differential expression in the ΔhsdM transcriptome, suggesting promoter HsdM-methylation directly influences transcription. Finally, comparative and functional analyses identified 351 sites hypervariable across isolates and numerous putative regulatory interactions. This multi-omic integration revealed features of methylomic variability in clinical isolates and provides a rational basis for hypothesizing the functions of DNA adenine methylation in MTBC physiology and adaptive evolution.

Background: Rifampicin (RIF) is a key first-line drug used to treat tuberculosis, a pulmonary disease caused by Mycobacterium tuberculosis. However antibiotic resistance to RIF is prevalent despite an apparent fitness cost. RIF resistance is primarily caused by mutations in the RIF resistance determining region in the rpoB gene, at the cost of slower growth in rich media. Compensatory mutations in the genes rpoA and rpoC have been shown to alleviate this fitness cost. These compensatory mutations may explain how RIF resistant strains have spread so rapidly. However, the effect of compensation on transmission is still unclear, partly because of uncertainty over which rpoABC mutations compensate for which RIF resistance markers. Objectives: We performed an association study on a globally representative set of 4309 whole genome sequenced clinical M. tuberculosis isolates to identify novel putative compensatory mutations, determine the prevalence of known and previously reported putative compensatory mutations, and determine which RIF resistance markers associate with these compensatory mutations. Results and Conclusions: Only 20.0% (216/1079) of RIF resistant isolates carried previously reported high-probability compensatory mutations, suggesting existence of other compensatory mutations. Using a strict phylogenetic approach, we identified 18 novel putative compensatory mutations in rpoC, rpoB, and rpoA. Novel and previously reported compensatory mutations were strongly associated with the RIFR marker rpoB:S450L, suggesting compensation may be specific to particular RIFR markers. These findings will aid identification of RIF-resistant M. tuberculosis strains with restored fitness. Such strains pose a greater risk of causing resistant outbreaks.

Point mutations in the rrs gene and eis promoter are known to confer resistance to second-line injectable drugs (SLIDs) amikacin (AMK), capreomycin (CAP), and kanamycin (KAN). While mutations in these canonical genes confer a majority of SLID-resistance, alternative mechanisms of resistance are not uncommon and threaten effective treatment decisions when using conventional molecular diagnostics. In total, 1184 clinical M. tuberculosis isolates from 7 countries were studied for genomic markers associated with phenotypic resistance. The markers rrs:A1401G and rrs:G1484T were associated with resistance to all three SLIDs, and three known markers in the eis promoter (eis:G-10A, eis:C-12T, and eis:C-14T) were similarly associated with kanamycin resistance (KAN-R). Among 325, 324, 270 AMK-R, CAP-R, and KAN-R isolates, 264 (81.2%), 250 (77.2%), and 249 (92.3%) harbored canonical mutations, respectively. Thirteen isolates harbored more than one canonical mutation. Canonical mutations did not account for 111 of the phenotypically resistant isolates. A gene-wise method identified three genes and promoters with mutations that, on aggregate, associated with unexplained resistance to at least one SLID. Our analysis associated whiB7 promoter mutations with KAN resistance, supporting clinical relevance for the previously demonstrated role of whiB7 overexpression in KAN resistance. We also provide evidence for the novel association of ppe51 (a gene previously associated with various antimicrobial compounds) with AMK resistance, and for the novel association of thrB with AMK and CAP resistance. The use of gene-wise association can provide additional insight, and therefore is recommended for identification of rare mechanisms of resistance when individual mutations carry insufficient statistical power.

No abstract

No abstract

Point mutations in the rrs gene and the eis promoter are known to confer resistance to the second-line injectable drugs (SLIDs) amikacin (AMK), capreomycin (CAP), and kanamycin (KAN). While mutations in these canonical genes confer the majority of SLID resistance, alternative mechanisms of resistance are not uncommon and threaten effective treatment decisions when using conventional molecular diagnostics.

No abstract

No abstract