Model-based integration of genomics and metabolomics reveals SNP functionality in Mycobacterium tuberculosis

Human tuberculosis is caused by members of the Mycobacterium tuberculosis complex (MTBC) and presents variable disease outcomes. The variation has primarily been attributed to host and environmental factors, but recent evidence indicates an additional role of genetic diversity among MTBC clinical strains. Here, we used metabolomics to unravel the potential role of genetic variations in conferring strain-specific adaptive capacity and vulnerability. To systematically identify functionality of single nucleotide polymorphisms (SNPs), we developed a constraint-based approach that integrates metabolomic and genomic data. Model-based predictions were systematically tested against independent metabolome data; they correctly classified SNP effects in pyruvate kinase and suggested a genetic basis for strain-specific sensitivity to the antibiotic para-aminosalicylic acid. Our method is broadly applicable to mutations in enzyme-encoding genes across microbial life, opening new possibilities for identifying strain-specific metabolic vulnerabilities that could lead to more selective treatment strategies.