GWAS for quantitative resistance phenotypes in Mycobacterium tuberculosis reveals resistance genes and regulatory regions

Drug resistance diagnostics that rely on the detection of resistance-related mutations could expedite patient care and TB eradication. We perform minimum inhibitory concentration testing for 12 anti-TB drugs together with Illumina whole-genome sequencing on 1452 clinical Mycobacterium tuberculosis (MTB) isolates. We evaluate genome-wide associations between mutations in MTB genes or non-coding regions and resistance, followed by validation in an independent data set of 792 patient isolates. We confirm associations at 13 non-canonical loci, with two involving non-coding regions. Promoter mutations are measured to have smaller average effects on resistance than gene body mutations. We estimate the heritability of the resistance phenotype to 11 anti-TB drugs and identify a lower than expected contribution from known resistance genes. This study highlights the complexity of the genomic mechanisms associated with the MTB resistance phenotype, including the relatively large number of potentially causal loci, and emphasizes the contribution of the non-coding portion of the genome.