Strain-specific genomic diversity in the Mycobacterium tuberculosis complex (MTBC) is an important factor in pathogenesis that may affect virulence, transmissibility, host response and emergence of drug resistance. Several systems have been proposed to classify MTBC strains into distinct lineages and families. Here, we investigate single-nucleotide polymorphisms (SNPs) as robust (stable) markers of genetic variation for phylogenetic analysis. We identify ~92k SNP across a global collection of 1,601 genomes. The SNP-based phylogeny is consistent with the gold-standard regions of difference (RD) classification system. Of the ~7k strain-specific SNPs identified, 62 markers are proposed to discriminate known circulating strains. This SNP-based barcode is the first to cover all main lineages, and classifies a greater number of sublineages than current alternatives. It may be used to classify clinical isolates to evaluate tools to control the disease, including therapeutics and vaccines whose effectiveness may vary by strain type.

Mycobacterium tuberculosis drug resistance (DR) challenges effective tuberculosis disease control. Current molecular tests examine limited numbers of mutations, and although whole genome sequencing approaches could fully characterise DR, data complexity has restricted their clinical application. A library (1,325 mutations) predictive of DR for 15 anti-tuberculosis drugs was compiled and validated for 11 of them using genomic-phenotypic data from 792 strains. A rapid online ‘TB-Profiler’ tool was developed to report DR and strain-type profiles directly from raw sequences. Using our DR mutation library, in silico diagnostic accuracy was superior to some commercial diagnostics and alternative databases. The library will facilitate sequence-based drug-susceptibility testing.

To characterize the genetic determinants of resistance to antituberculosis drugs, we performed a genome-wide association study (GWAS) of 6,465 Mycobacterium tuberculosis clinical isolates from more than 30 countries. A GWAS approach within a mixed-regression framework was followed by a phylogenetics-based test for independent mutations. In addition to mutations in established and recently described resistance-associated genes, novel mutations were discovered for resistance to cycloserine, ethionamide and para-aminosalicylic acid. The capacity to detect mutations associated with resistance to ethionamide, pyrazinamide, capreomycin, cycloserine and para-aminosalicylic acid was enhanced by inclusion of insertions and deletions. Odds ratios for mutations within candidate genes were found to reflect levels of resistance. New epistatic relationships between candidate drug-resistance-associated genes were identified. Findings also suggest the involvement of efflux pumps (drrA and Rv2688c) in the emergence of resistance. This study will inform the design of new diagnostic tests and expedite the investigation of resistance and compensatory epistatic mechanisms.

Multidrug resistant (MDR) tuberculosis is caused by Mycobacterium tuberculosis resistant to isoniazid and rifampicin, the two most effective drugs used in tuberculosis therapy. Here, we investigated the mechanism by which resistance towards isoniazid develops and how overexpression of efflux pumps favors accumulation of mutations in isoniazid targets, thus establishing a MDR phenotype. The study was based on the in vitro induction of an isoniazid resistant phenotype by prolonged serial exposure of M. tuberculosis strains to the critical concentration of isoniazid employed for determination of drug susceptibility testing in clinical isolates. Results show that susceptible and rifampicin monoresistant strains exposed to this concentration become resistant to isoniazid after three weeks; and that resistance observed for the majority of these strains could be reduced by means of efflux pumps inhibitors. RT-qPCR assessment of efflux pump genes expression showed overexpression of all tested genes. Enhanced real-time efflux of ethidium bromide, a common efflux pump substrate, was also observed, showing a clear relation between overexpression of the genes and increased efflux pump function. Further exposure to isoniazid resulted in the selection and stabilization of spontaneous mutations and deletions in the katG gene along with sustained increased efflux activity. Together, results demonstrate the relevance of efflux pumps as one of the factors of isoniazid resistance in M. tuberculosis. These results support the hypothesis that activity of efflux pumps allows the maintenance of an isoniazid resistant population in a sub-optimally treated patient from which isoniazid genetically resistant mutants emerge. Therefore, the use of inhibitors of efflux should be considered in the development of new therapeutic strategies for preventing the emergence of MDR-TB during treatment.

On the basis of population genomic and phylogeographic analyses of 1669 Mycobacterium tuberculosis lineage 4 (L4) genomes, we find that dispersal of L4 has been completely dominated by historical migrations out of Europe. We demonstrate an intimate temporal relationship between European colonial expansion into Africa and the Americas and the spread of L4 tuberculosis (TB). Markedly, in the age of antibiotics, mutations conferring antimicrobial resistance overwhelmingly emerged locally (at the level of nations), with minimal cross-border transmission of resistance. The latter finding was found to reflect the relatively recent emergence of these mutations, as a similar degree of local restriction was observed for susceptible variants emerging on comparable time scales. The restricted international transmission of drug-resistant TB suggests that containment efforts at the level of individual countries could be successful.

The purpose of this study was to determine the levels of isoniazid and ethionamide resistance and to identify associated mutations in endemic multidrug-resistant (MDR) strains of Mycobacterium tuberculosis from the Lisbon metropolitan area, Portugal. Seventeen clinical MDR tuberculosis (TB) strains were characterized by standard and semi-quantitative drug susceptibility testing to assess the level of isoniazid and ethionamide resistance. The genes katG, inhA, ethA and ndh were screened for mutations. All strains were genotyped by 24 loci mycobacterial interspersed repetitive unit–variable number of tandem repeats (MIRU-VNTR) analysis. All strains showed high-level resistance to both isoniazid (>1 mg/L) and ethionamide (>25 mg/L). MIRU-VNTR typing revealed the presence of two main clusters, Lisboa3 and Q1, in 16/17 strains, all of which showed the C−15T mutation in the promoter region of the inhA gene. The 16 strains belong to the Latino-American-Mediterranean (LAM) genotype and the other strain belongs to the Beijing genotype. Sequencing of the inhA open reading frame revealed that the 16 strains also had mutations in the structural region of the gene, leading to the S94A substitution in 9 strains and the I194T substitution in 7 strains. The results reveal that the presence of a mutation in the inhA regulatory region together with a mutation in the inhA coding region can lead to the development of high-level isoniazid resistance and cross-resistance to ethionamide among the MDR-TB strains circulating in Lisbon. This mutational pattern also hints to a possible involvement of strain-specific factors that could be a feature of the Portuguese MDR-TB strains where the LAM family is the major circulating genotype.