Efflux pump gene single-nucleotide variants associated with resistance in Mycobacterium tuberculosis isolates with discrepant drug genotypes

Introduction Single-nucleotide variants (SNVs) in Mycobacterium tuberculosis (M. tuberculosis) genomes can predict multidrug resistance (MDR) but not all phenotype-genotype correlations can be explained. We investigated SNVs in efflux pumps (EPs) in the context of M. tuberculosis drug resistance. Methods We analysed 2221 M. tuberculosis genomes from 1432 susceptible and 200 MDR, 172 pre-extensively drug resistant (XDR) and 417 XDR isolates. Analysis of 47 EP genes was conducted using MTB-VCF, an in-house bioinformatics pipeline. SNVs were categorized according to their SIFT/Polyphen scores. Resistance genotypes were also called using the TB-Profiler tool. Results Genome comparisons between susceptible and drug resistant (DR) isolates identified 418 unique SNVs in EP of which; 53.5% were in MDR, 68.9% in pre-XDR and 61.3% in XDR isolates. Twenty EPs had unique SNVs with a high SIFT/PolyPhen score, comprising 38 unique SNVs. Sixteen SNVs across 12 EP genes were significantly associated with drug resistance and enriched in pre-XDR and XDR strains. These comprised 12 previously reported SNVs (in Rv0191, Rv0507, Rv0676, Rv1217, Rv1218, Rv1273, Rv1458, Rv1819, and Rv2688) and 4 novel SNVs (in Rv1877 and Rv2333). We investigated their presence in genomes of 52 MDR isolates with phenotype-genotype discrepancies to rifampicin (RIF), isoniazid (INH), or fluoroquinolones. SNVs associated with RIF and INH (Rv1217_1218, Rv1819, Rv0450, Rv1458, Rv3827, Rv0507, Rv0676, Rv1273, and Rv2333), and with fluoroquinolone (Rv2688) resistance were present in these discrepant strains. Conclusions Considering SNVs in EPs as part of M. tuberculosis genome-based resistance interpretation may add value, especially in evaluation of XDR resistance in strains with phenotype-genotype discrepancies.