Genome-wide analysis of Mycobacterium tuberculosis polymorphisms reveals lineage-specific associations with drug resistance

Background

Continuing evolution of the Mycobacterium tuberculosis (Mtb) complex genomes associated with resistance to anti-tuberculosis drugs is threatening tuberculosis disease control efforts. Both multi- and extensively drug resistant Mtb (MDR and XDR, respectively) are increasing in prevalence, but the full set of Mtb genes involved are not known. There is a need for increased sensitivity of genome-wide approaches in order to elucidate the genetic basis of anti-microbial drug resistance and gain a more detailed understanding of Mtb genome evolution in a context of widespread antimicrobial therapy. Population structure within the Mtb complex, due to clonal expansion, lack of lateral gene transfer and low levels of recombination between lineages, may be reducing statistical power to detect drug resistance associated variants.

Results

To investigate the effect of lineage-specific effects on the identification of drug resistance associations, we applied genome-wide association study (GWAS) and convergence-based (PhyC) methods to multiple drug resistance phenotypes of a global dataset of Mtb lineages 2 and 4, using both lineage-wise and combined approaches. We identify both well-established drug resistance variants and novel associations; uniquely identifying associations for both lineage-specific and -combined GWAS analyses. We report 17 potential novel associations between antimicrobial resistance phenotypes and Mtb genomic variants.

Conclusions

For GWAS, both lineage-specific and -combined analyses are useful, whereas PhyC may perform better in contexts of greater diversity. Unique associations with XDR in lineage-specific analyses provide evidence of diverging evolutionary trajectories between lineages 2 and 4 in response to antimicrobial drug therapy.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5615-3) contains supplementary material, which is available to authorized users.

Longitudinal whole genome analysis of pre and post drug treatment Mycobacterium tuberculosis isolates reveals progressive steps to drug resistance

Tuberculosis (TB) is one of the leading causes of death due to an infectious disease in the world. Understanding the mechanisms of drug resistance has become pivotal in the detection and treatment of newly emerging resistant TB cases. We have analyzed three pairs of Mycobacterium tuberculosis strains pre- and post-drug treatment to identify mutations involved in the progression of resistance to the drugs rifampicin and isoniazid. In the rifampicin resistant strain, we confirmed a mutation in rpoB (S450L) that is known to confer resistance to rifampicin. We discovered a novel L101R mutation in the katG gene of an isoniazid resistant strain, which may directly contribute to isoniazid resistance due to the proximity of the mutation to the katG isoniazid-activating site. Another isoniazid resistant strain had a rare mutation in the start codon of katG. We also identified a number of mutations in each longitudinal pair, such as toxin-antitoxin mutations that may influence the progression towards resistance or may play a role in compensatory fitness. These findings improve our knowledge of drug resistance progression during therapy and provide a methodology to monitor longitudinal strains using whole genome sequencing, polymorphism comparison, and functional annotation.

Whole genome sequencing reveals genomic heterogeneity and antibiotic purification in Mycobacterium tuberculosis isolates

BACKGROUND

Whole genome sequencing has revolutionised the interrogation of mycobacterial genomes. Recent studies have reported conflicting findings on the genomic stability of Mycobacterium tuberculosis during the evolution of drug resistance. In an age where whole genome sequencing is increasingly relied upon for defining the structure of bacterial genomes, it is important to investigate the reliability of next generation sequencing to identify clonal variants present in a minor percentage of the population. This study aimed to define a reliable cut-off for identification of low frequency sequence variants and to subsequently investigate genetic heterogeneity and the evolution of drug resistance in M. tuberculosis.

METHODS

Genomic DNA was isolated from single colonies from 14 rifampicin mono-resistant M. tuberculosis isolates, as well as the primary cultures and follow up MDR cultures from two of these patients. The whole genomes of the M. tuberculosis isolates were sequenced using either the Illumina MiSeq or Illumina HiSeq platforms. Sequences were analysed with an in-house pipeline.

RESULTS

Using next-generation sequencing in combination with Sanger sequencing and statistical analysis we defined a read frequency cut-off of 30% to identify low frequency M. tuberculosis variants with high confidence. Using this cut-off we demonstrated a high rate of genetic diversity between single colonies isolated from one population, showing that by using the current sequencing technology, single colonies are not a true reflection of the genetic diversity within a whole population and vice versa. We further showed that numerous heterogeneous variants emerge and then disappear during the evolution of isoniazid resistance within individual patients. Our findings allowed us to formulate a model for the selective bottleneck which occurs during the course of infection, acting as a genomic purification event.

CONCLUSIONS

Our study demonstrated true levels of genetic diversity within an M. tuberculosis population and showed that genetic diversity may be re-defined when a selective pressure, such as drug exposure, is imposed on M. tuberculosis populations during the course of infection. This suggests that the genome of M. tuberculosis is more dynamic than previously thought, suggesting preparedness to respond to a changing environment.

Genetic markers, genotyping methods & next generation sequencing in Mycobacterium tuberculosis

Molecular epidemiology (ME) is one of the main areas in tuberculosis research which is widely used to study the transmission epidemics and outbreaks of tubercle bacilli. It exploits the presence of various polymorphisms in the genome of the bacteria that can be widely used as genetic markers. Many DNA typing methods apply these genetic markers to differentiate various strains and to study the evolutionary relationships between them. The three widely used genotyping tools to differentiate Mycobacterium tuberculosis strains are IS6110 restriction fragment length polymorphism (RFLP), spacer oligotyping (Spoligotyping), and mycobacterial interspersed repeat units - variable number of tandem repeats (MIRU-VNTR). A new prospect towards ME was introduced with the development of whole genome sequencing (WGS) and the next generation sequencing (NGS) methods, where the entire genome is sequenced that not only helps in pointing out minute differences between the various sequences but also saves time and the cost. NGS is also found to be useful in identifying single nucleotide polymorphisms (SNPs), comparative genomics and also various aspects about transmission dynamics. These techniques enable the identification of mycobacterial strains and also facilitate the study of their phylogenetic and evolutionary traits.

Reengineering redox sensitive GFP to measure mycothiol redox potential of Mycobacterium tuberculosis during infection

Mycobacterium tuberculosis (Mtb) survives under oxidatively hostile environments encountered inside host phagocytes. To protect itself from oxidative stress, Mtb produces millimolar concentrations of mycothiol (MSH), which functions as a major cytoplasmic redox buffer. Here, we introduce a novel system for real-time imaging of mycothiol redox potential (EMSH ) within Mtb cells during infection. We demonstrate that coupling of Mtb MSH-dependent oxidoreductase (mycoredoxin-1; Mrx1) to redox-sensitive GFP (roGFP2; Mrx1-roGFP2) allowed measurement of dynamic changes in intramycobacterial EMSH with unprecedented sensitivity and specificity. Using Mrx1-roGFP2, we report the first quantitative measurements of EMSH in diverse mycobacterial species, genetic mutants, and drug-resistant patient isolates. These cellular studies reveal, for the first time, that the environment inside macrophages and sub-vacuolar compartments induces heterogeneity in EMSH of the Mtb population. Further application of this new biosensor demonstrates that treatment of Mtb infected macrophage with anti-tuberculosis (TB) drugs induces oxidative shift in EMSH , suggesting that the intramacrophage milieu and antibiotics cooperatively disrupt the MSH homeostasis to exert efficient Mtb killing. Lastly, we analyze the membrane integrity of Mtb cells with varied EMSH during infection and show that subpopulation with higher EMSH are susceptible to clinically relevant antibiotics, whereas lower EMSH promotes antibiotic tolerance. Together, these data suggest the importance of MSH redox signaling in modulating mycobacterial survival following treatment with anti-TB drugs. We anticipate that Mrx1-roGFP2 will be a major contributor to our understanding of redox biology of Mtb and will lead to novel strategies to target redox metabolism for controlling Mtb persistence.