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

Large-scale statistical analysis of Mycobacterium tuberculosis genome sequences identifies compensatory mutations associated with multi-drug resistance

Tuberculosis (TB), caused by Mycobacterium tuberculosis, has a significant impact on global health worldwide. The development of multi-drug resistant strains that are resistant to the first-line drugs isoniazid and rifampicin threatens public health security. Rifampicin and isoniazid resistance are largely underpinned by mutations in rpoB and katG respectively and are associated with fitness costs. Compensatory mutations are considered to alleviate these fitness costs and have been observed in rpoC/rpoA (rifampicin) and oxyR'-ahpC (isoniazid). We developed a framework (CompMut-TB) to detect compensatory mutations from whole genome sequences from a large dataset comprised of 18,396 M. tuberculosis samples. We performed association analysis (Fisher's exact tests) to identify pairs of mutations that are associated with drug-resistance, followed by mediation analysis to identify complementary or full mediators of drug-resistance. The analyses revealed several potential mutations in rpoC (N = 47), rpoA (N = 4), and oxyR'-ahpC (N = 7) that were considered either 'highly likely' or 'likely' to confer compensatory effects on drug-resistance, including mutations that have previously been reported and validated. Overall, we have developed the CompMut-TB framework which can assist with identifying compensatory mutations which is important for more precise genome-based profiling of drug-resistant TB strains and to further understanding of the evolutionary mechanisms that underpin drug-resistance.

Genetic Characterization and Population Structure of Drug-Resistant Mycobacterium tuberculosis Isolated from Brazilian Patients Using Whole-Genome Sequencing

The present study aimed to determine the genetic diversity of isolates of Mycobacterium tuberculosis (Mtb) from presumed drug-resistant tuberculosis patients from several states of Brazil. The isolates had been submitted to conventional drug susceptibility testing for first- and second-line drugs. Multidrug-resistant (MDR-TB) (54.8%) was the most frequent phenotypic resistance profile, in addition to an important high frequency of pre-extensive resistance (p-XDR-TB) (9.2%). Using whole-genome sequencing (WGS), we characterized 298 Mtb isolates from Brazil. Besides the analysis of genotype distribution and possible correlations between molecular and clinical data, we determined the performance of an in-house WGS pipeline with other online pipelines for Mtb lineages and drug resistance profile definitions. Sub-lineage 4.3 (52%) was the most frequent genotype, and the genomic approach revealed a p-XDR-TB level of 22.5%. We detected twenty novel mutations in three resistance genes, and six of these were observed in eight phenotypically resistant isolates. A cluster analysis of 170 isolates showed that 43.5% of the TB patients belonged to 24 genomic clusters, suggesting considerable ongoing transmission of DR-TB, including two interstate transmissions. The in-house WGS pipeline showed the best overall performance in drug resistance prediction, presenting the best accuracy values for five of the nine drugs tested. Significant associations were observed between suffering from fatal disease and genotypic p-XDR-TB (p = 0.03) and either phenotypic (p = 0.006) or genotypic (p = 0.0007) ethambutol resistance. The use of WGS analysis improved our understanding of the population structure of MTBC in Brazil and the genetic and clinical data correlations and demonstrated its utility for surveillance efforts regarding the spread of DR-TB, hopefully helping to avoid the emergence of even more resistant strains and to reduce TB incidence and mortality rates.

Comparative Study on Tuberculosis Drug Resistance and Molecular Detection Methods Among Different Mycobacterium Tuberculosis Lineages

Purpose

This study aims to compare drug resistance and detection efficacy across different Mycobacterium tuberculosis lineages, offering insights for precise treatment and molecular diagnosis.

Methods

161 strains of Mycobacterium tuberculosis (M.tb) were tested for drug resistance using Phenotypic Drug Susceptibility Testing (pDST), High-Resolution Melting analysis (HRM), and Whole Genome Sequencing (WGS) methods. The main focus was on evaluating the accuracy of different methods for detecting resistance to rifampicin (RIF), isoniazid (INH), and streptomycin (SM).

Results

Among the 161 strains of M.tb, 83.85% (135/161) were fully sensitive to RIF, INH, and SM according to pDST, and the rate of multidrug resistance was 4.35% (7/161). The drug resistance rates of lineage 2 M.tb to the three drugs (26/219, 11.87%) were significantly higher than those of non-lineage 2 M.tb (12/264, 4.45%) (P<0.05). Compared with pDST, WGS had a sensitivity of 100%, 94.12%, and 92.31% and a specificity of 100%, 99.31%, and 98.65% for RIF, INH, and SM, respectively, with no significant difference. The sensitivity of HRM for RIF, INH, and SM was 87.50%, 52.94%, and 76.92%, respectively, while the specificity was 96.08%, 99.31%, and 99.32%, respectively. The sensitivity of HRM for detecting INH resistance was significantly lower than that of pDST (P=0.039). Compared with HRM, WGS increased the sensitivity of RIF, INH, and SM by 12.50%, 41.18%, and 15.38%, respectively.

Conclusion

There are significant differences in drug resistance rates among different lineages of M.tb, with lineage 2 having higher rates of RIF, INH, and SM resistance than lineages 3 and 4. The sensitivity of HRM is far lower than that of pDST, and currently, the accuracy of HRM is not sufficient to replace pDST. WGS has no significant difference in detecting drug resistance compared with pDST but can identify new anti-tuberculosis drug-resistant mutations, providing effective guidance for clinical decision-making.

Horizontal Gene Transfer and Drug Resistance Involving Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) acquires drug resistance at a rate comparable to that of bacterial pathogens that replicate much faster and have a higher mutation rate. One explanation for this rapid acquisition of drug resistance in Mtb is that drug resistance may evolve in other fast-replicating mycobacteria and then be transferred to Mtb through horizontal gene transfer (HGT). This paper aims to address three questions. First, does HGT occur between Mtb and other mycobacterial species? Second, what genes after HGT tend to survive in the recipient genome? Third, does HGT contribute to antibiotic resistance in Mtb? I present a conceptual framework for detecting HGT and analyze 39 ribosomal protein genes, 23S and 16S ribosomal RNA genes, as well as several genes targeted by antibiotics against Mtb, from 43 genomes representing all major groups within Mycobacterium. I also included mgtC and the insertion sequence IS6110 that were previously reported to be involved in HGT. The insertion sequence IS6110 shows clearly that the Mtb complex participates in HGT. However, the horizontal transferability of genes depends on gene function, as was previously hypothesized. HGT is not observed in functionally important genes such as ribosomal protein genes, rRNA genes, and other genes chosen as drug targets. This pattern can be explained by differential selection against functionally important and unimportant genes after HGT. Functionally unimportant genes such as IS6110 are not strongly selected against, so HGT events involving such genes are visible. For functionally important genes, a horizontally transferred diverged homologue from a different species may not work as well as the native counterpart, so the HGT event involving such genes is strongly selected against and eliminated, rendering them invisible to us. In short, while HGT involving the Mtb complex occurs, antibiotic resistance in the Mtb complex arose from mutations in those drug-targeted genes within the Mtb complex and was not gained through HGT.

Molecular epidemiology and transmission dynamics of multi-drug resistant tuberculosis strains using whole genome sequencing in the Amhara region, Ethiopia

BACKGROUND

Drug resistant Mycobacterium tuberculosis prevention and care is a major challenge in Ethiopia. The World health organization has designated Ethiopia as one of the 30 high burden multi-drug resistant tuberculosis (MDR-TB) countries. There is limited information regarding genetic diversity and transmission dynamics of MDR-TB in Ethiopia.

OBJECTIVE

To investigate the molecular epidemiology and transmission dynamics of MDR-TB strains using whole genome sequence (WGS) in the Amhara region.

METHODS

Forty-five MDR-TB clinical isolates from Amhara region were collected between 2016 and 2018, and characterized using WGS and 24-loci Mycobacterium Interspersed Repetitive Units Variable Number of Tandem Repeats (MIRU-VNTR) typing. Clusters were defined based on the maximum distance of 12 single nucleotide polymorphisms (SNPs) or alleles as the upper threshold of genomic relatedness. Five or less SNPs or alleles distance or identical 24-loci VNTR typing is denoted as surrogate marker for recent transmission.

RESULTS

Forty-one of the 45 isolates were analyzed by WGS and 44% (18/41) of the isolates were distributed into 4 clusters. Of the 41 MDR-TB isolates, 58.5% were classified as lineage 4, 36.5% lineage 3 and 5% lineage 1. Overall, TUR genotype (54%) was the predominant in MDR-TB strains. 41% (17/41) of the isolates were clustered into four WGS groups and the remaining isolates were unique strains. The predominant cluster (Cluster 1) was composed of nine isolates belonging to lineage 4 and of these, four isolates were in the recent transmission links.

CONCLUSIONS

Majority of MDR-TB strain cluster and predominance of TUR lineage in the Amhara region give rise to concerns for possible ongoing transmission. Efforts to strengthen TB laboratory to advance diagnosis, intensified active case finding, and expanded contact tracing activities are needed in order to improve rapid diagnosis and initiate early treatment. This would lead to the interruption of the transmission chain and stop the spread of MDR-TB in the Amhara region.

Comparison of in silico predicted Mycobacterium tuberculosis spoligotypes and lineages from whole genome sequencing data

Bacterial strain-types in the Mycobacterium tuberculosis complex underlie tuberculosis disease, and have been associated with drug resistance, transmissibility, virulence, and host-pathogen interactions. Spoligotyping was developed as a molecular genotyping technique used to determine strain-types, though recent advances in whole genome sequencing (WGS) technology have led to their characterization using SNP-based sub-lineage nomenclature. Notwithstanding, spoligotyping remains an important tool and there is a need to study the congruence between spoligotyping-based and SNP-based sub-lineage assignation. To achieve this, an in silico spoligotype prediction method ("Spolpred2") was developed and integrated into TB-Profiler. Lineage and spoligotype predictions were generated for > 28 k isolates and the overlap between strain-types was characterized. Major spoligotype families detected were Beijing (25.6%), T (18.6%), LAM (13.1%), CAS (9.4%), and EAI (8.3%), and these broadly followed known geographic distributions. Most spoligotypes were perfectly correlated with the main MTBC lineages (L1-L7, plus animal). Conversely, at lower levels of the sub-lineage system, the relationship breaks down, with only 65% of spoligotypes being perfectly associated with a sub-lineage at the second or subsequent levels of the hierarchy. Our work supports the use of spoligotyping (membrane or WGS-based) for low-resolution surveillance, and WGS or SNP-based systems for higher-resolution studies.

In vitro modeling of isoniazid resistance mechanisms in Mycobacterium tuberculosis H37Rv

Introduction

Mycobacterium tuberculosis (MTB), the causative agent of tuberculosis, has been a global threat to human beings for several decades. Treating tuberculosis has become more difficult as the prevalence of drug-resistant tuberculosis has increased globally. Evidence suggests that the comprehensive landscape of resistance mechanisms in MTB is ambiguous. More importantly, little is known regarding the series of events connected to resistance mechanisms in MTB before exposure to anti-TB drugs, during exposure to the drugs, and finally, when the MTB becomes resistant after exposure, upon analyses of its genome.

Methods

We used the wild-type strain of MTB (H37Rv) in an in vitro model for generating induced resistance using a sub-inhibitory concentration of isoniazid, and the generated resistance-associated variants (RAVs) were identified using the whole genome sequencing method.

Results

The detection of an inhA promoter mutation (fabG1-15C>T), which results in increased production of InhA protein, was found to be a major mechanism for developing resistance to isoniazid in the first place. We observed adaptation of MTB resistance mechanisms in high isoniazid stress by alteration and abolishment of KatG due to the detection of katG S315N, the common region of mutation that confers isoniazid resistance, along with katG K414N, katG N138S, and katG A162E. Furthermore, we detected the ahpC-72C>T and ahpC 21C>A mutations, but further investigation is needed to determine their role in compensating for the loss of KatG activity.

Discussion

This suggests that increased InhA production is the main mechanism where there are low levels of isoniazid, whereas the alteration of KatG was found to be utilized in mycobacterium with a high concentration of isoniazid. Our work demonstrates that this in vitro approach of generating induced resistance could provide clinically relevant information after the fabG1-15C>T mutation, which is the common mutation found in clinical isolates. Moreover, other mutations detected in this work can also be found in clinical isolates. These findings may shed light on the impact of isoniazid in generating RAV and the resistance mechanism scenario that mycobacterium used under various isoniazid-pressuring conditions. More research is needed to understand better the role of RAV and mechanical resistance events within the mycobacterium genome in promoting a promising drug prediction platform that could lead to the right treatment for patients with MDR-TB and XDR-TB.

Genetic mutations underlying isoniazid-resistant Mycobacterium tuberculosis in Khyber Pakhtunkhwa, Pakistan

Tuberculosis, caused by Mycobacterium tuberculosis, is a major public health issue in Pakistan. Isoniazid is a first-line pro-drug that requires activation through an enzyme called catalase peroxidase, but is subject to widespread resistance, driven by mutations in katG and inhA genes and other loci with compensatory effects (e.g., ahpC). Here, we used whole genome sequencing data from 51 M. tuberculosis isolates collected from Khyber Pakhtunkhwa province (years 2016-2019; all isoniazid phenotypically resistant) to investigate the genetic diversity of mutations in isoniazid candidate genes. The most common mutations underlying resistance were katG S315T (37/51), fabG1 -15C>T (13/51; inhA promoter), and inhA -154G>A (7/51). Other less common mutations (n < 5) were also identified in katG (R128Q, V1A, W505*, A109T, D311G) and candidate compensatory genes ahpC (-54C>T, -51G>A) and oxyS (M249T). Using DynaMut2 software, the mutants exhibited various degrees of stability and flexibility on protein structures, with some katG mutations leading to a decrease in KatG protein flexibility. Overall, the characterisation of circulating isoniazid resistant-linked mutations will assist in drug resistant TB management and control activities in a highly endemic area of Pakistan.

Mycobacterium tuberculosis functional genetic diversity, altered drug sensitivity, and precision medicine

In the face of the unrelenting global burden of tuberculosis (TB), antibiotics remain our most effective tools to save lives and control the spread of Mycobacterium tuberculosis (Mtb). However, we confront a dual challenge in our use of antibiotics: simplifying and shortening the TB drug regimen while also limiting the emergence and propagation of antibiotic resistance. This task is now more feasible due to the increasing availability of bacterial genomic data at or near the point of care. These resources create an opportunity to envision how integration of bacterial genetic determinants of antibiotic response into treatment algorithms might transform TB care. Historically, Mtb drug resistance studies focused on mutations in genes encoding antibiotic targets and the resulting increases in the minimal inhibitory concentrations (MICs) above a breakpoint value. But recent progress in elucidating the effects of functional genetic diversity in Mtb has revealed various genetic loci that are associated with drug phenotypes such as low-level MIC increases and tolerance which predict the development of resistance and treatment failure. As a result, we are now poised to advance precision medicine approaches in TB treatment. By incorporating information regarding Mtb genetic characteristics into the development of drug regimens, clinical care which tailors antibiotic treatment to maximize the likelihood of success has come into reach.

Genome-wide association studies of global Mycobacterium tuberculosis resistance to 13 antimicrobials in 10,228 genomes identify new resistance mechanisms

The emergence of drug-resistant tuberculosis is a major global public health concern that threatens the ability to control the disease. Whole-genome sequencing as a tool to rapidly diagnose resistant infections can transform patient treatment and clinical practice. While resistance mechanisms are well understood for some drugs, there are likely many mechanisms yet to be uncovered, particularly for new and repurposed drugs. We sequenced 10,228 Mycobacterium tuberculosis (MTB) isolates worldwide and determined the minimum inhibitory concentration (MIC) on a grid of 2-fold concentration dilutions for 13 antimicrobials using quantitative microtiter plate assays. We performed oligopeptide- and oligonucleotide-based genome-wide association studies using linear mixed models to discover resistance-conferring mechanisms not currently catalogued. Use of MIC over binary resistance phenotypes increased sample heritability for the new and repurposed drugs by 26% to 37%, increasing our ability to detect novel associations. For all drugs, we discovered uncatalogued variants associated with MIC, including in the Rv1218c promoter binding site of the transcriptional repressor Rv1219c (isoniazid), upstream of the vapBC20 operon that cleaves 23S rRNA (linezolid) and in the region encoding an α-helix lining the active site of Cyp142 (clofazimine, all p < 10-7.7). We observed that artefactual signals of cross-resistance could be unravelled based on the relative effect size on MIC. Our study demonstrates the ability of very large-scale studies to substantially improve our knowledge of genetic variants associated with antimicrobial resistance in M. tuberculosis.

Variants in Bedaquiline-Candidate-Resistance Genes: Prevalence in Bedaquiline-Naive Patients, Effect on MIC, and Association with Mycobacterium tuberculosis Lineage

Studies have shown that variants in bedaquiline-resistance genes can occur in isolates from bedaquiline-naive patients. We assessed the prevalence of variants in all bedaquiline-candidate-resistance genes in bedaquiline-naive patients, investigated the association between these variants and lineage, and the effect on phenotype. We used whole-genome sequencing to identify variants in bedaquiline-resistance genes in isolates from 509 bedaquiline treatment naive South African tuberculosis patients. A phylogenetic tree was constructed to investigate the association with the isolate lineage background. Bedaquiline MIC was determined using the UKMYC6 microtiter assay. Variants were identified in 502 of 509 isolates (98.6%), with the highest (85%) prevalence of variants in the Rv0676c (mmpL5) gene. We identified 36 unique variants, including 19 variants not reported previously. Only four isolates had a bedaquiline MIC equal to or above the epidemiological cut-off value of 0.25 μg/mL. Phylogenetic analysis showed that 14 of the 15 variants observed more than once occurred monophyletically in one Mycobacterium tuberculosis (sub)lineage. The bedaquiline MIC differed between isolates belonging to lineage 2 and 4 (Fisher's exact test, P = 0.0004). The prevalence of variants in bedaquiline-resistance genes in isolates from bedaquiline-naive patients is high, but very few (<2%) isolates were phenotypically resistant. We found an association between variants in bedaquiline resistance genes and Mycobacterium tuberculosis (sub)lineage, resulting in a lineage-dependent difference in bedaquiline phenotype. Future studies should investigate the impact of the presence of variants on bedaquiline-resistance acquisition and treatment outcome.

Characterisation of drug-resistant Mycobacterium tuberculosis mutations and transmission in Pakistan

Tuberculosis, caused by Mycobacterium tuberculosis, is a high-burden disease in Pakistan, with multi-drug (MDR) and extensive-drug (XDR) resistance, complicating infection control. Whole genome sequencing (WGS) of M. tuberculosis is being used to infer lineages (strain-types), drug resistance mutations, and transmission patterns-all informing infection control and clinical decision making. Here we analyse WGS data on 535 M. tuberculosis isolates sourced across Pakistan between years 2003 and 2020, to understand the circulating strain-types and mutations related to 12 anti-TB drugs, as well as identify transmission clusters. Most isolates belonged to lineage 3 (n = 397; 74.2%) strain-types, and were MDR (n = 328; 61.3%) and (pre-)XDR (n = 113; 21.1%). By inferring close genomic relatedness between isolates (< 10-SNPs difference), there was evidence of M. tuberculosis transmission, with 55 clusters formed consisting of a total of 169 isolates. Three clusters consist of M. tuberculosis that are similar to isolates found outside of Pakistan. A genome-wide association analysis comparing 'transmitted' and 'non-transmitted' isolate groups, revealed the nusG gene as most significantly associated with a potential transmissible phenotype (P = 5.8 × 10^-10). Overall, our study provides important insights into M. tuberculosis genetic diversity and transmission in Pakistan, including providing information on circulating drug resistance mutations for monitoring activities and clinical decision making.

Distribution of Common and Rare Genetic Markers of Second-Line-Injectable-Drug Resistance in Mycobacterium tuberculosis Revealed by a Genome-Wide Association Study

Point mutations in the rrs gene and the eis promoter are known to confer resistance to the second-line injectable drugs (SLIDs) amikacin (AMK), capreomycin (CAP), and kanamycin (KAN). While mutations in these canonical genes confer the majority of SLID resistance, alternative mechanisms of resistance are not uncommon and threaten effective treatment decisions when using conventional molecular diagnostics. In total, 1,184 clinical Mycobacterium tuberculosis isolates from 7 countries were studied for genomic markers associated with phenotypic resistance. The markers rrs:A1401G and rrs:G1484T were associated with resistance to all three SLIDs, and three known markers in the eis promoter (eis:G-10A, eis:C-12T, and eis:C-14T) were similarly associated with kanamycin resistance (KAN-R). Among 325, 324, and 270 AMK-R, CAP-R, and KAN-R isolates, 274 (84.3%), 250 (77.2%), and 249 (92.3%) harbored canonical mutations, respectively. Thirteen isolates harbored more than one canonical mutation. Canonical mutations did not account for 103 of the phenotypically resistant isolates. A genome-wide association study identified three genes and promoters with mutations that, on aggregate, were associated with unexplained resistance to at least one SLID. Our analysis associated whiB7 5'-untranslated-region mutations with KAN resistance, supporting clinical relevance for this previously demonstrated mechanism of KAN resistance. We also provide evidence for the novel association of CAP resistance with the promoter of the Rv2680-Rv2681 operon, which encodes an exoribonuclease that may influence the binding of CAP to the ribosome. Aggregating mutations by gene can provide additional insight and therefore is recommended for identifying rare mechanisms of resistance when individual mutations carry insufficient statistical power.

An evolutionary functional genomics approach identifies novel candidate regions involved in isoniazid resistance in Mycobacterium tuberculosis

Efforts to eradicate tuberculosis are hampered by the rise and spread of antibiotic resistance. Several large-scale projects have aimed to specifically link clinical mutations to resistance phenotypes, but they were limited in both their explanatory and predictive powers. Here, we combine functional genomics and phylogenetic associations using clinical strain genomes to decipher the architecture of isoniazid resistance and search for new resistance determinants. This approach has allowed us to confirm the main target route of the antibiotic, determine the clinical relevance of redox metabolism as an isoniazid resistance mechanism and identify novel candidate genes harboring resistance mutations in strains with previously unexplained isoniazid resistance. This approach can be useful for characterizing how the tuberculosis bacilli acquire resistance to new antibiotics and how to forestall them.

Victoria Furió et al. apply functional genomics and evolutionary analyses to the study of antibiotic resistance in tuberculosis. Focusing on isoniazid resistance and using genomic data from clinical strains, they identify novel candidate genes with resistance mutations and further uncover the mechanisms underlying drug resistance.

Characterization of Mutations Associated with Streptomycin Resistance in Multidrug-Resistant Mycobacterium tuberculosis in Zambia

Streptomycin (STR) is recommended for the management of multidrug-resistant tuberculosis (MDR-TB). Streptomycin resistance-conferring mutation types and frequency are shown to be influenced by genotypes of circulating strains in a population. This study aimed to characterize the mutations in MDR-TB isolates and examine their relationship with the genotypes in Zambia. A total of 138 MDR-TB isolates stored at the University Teaching Hospital Tuberculosis Reference Laboratory in Zambia were analyzed using spoligotyping and sequencing of STR resistance-associated genes. Streptomycin resistance was observed in 65.9% (91/138) of MDR-TB isolates. Mutations in rpsL, rrs, and gidB accounted for 33%, 12.1%, and 49.5%, respectively. Amino acid substitution K43R in rpsL was strongly associated with the CAS1_Kili genotype (p < 0.0001). The combination of three genes could predict 91.2% of STR resistance. Clustering of isolates based on resistance-conferring mutations and spoligotyping was observed. The clustering of isolates suggests that the increase in STR-resistant MDR-TB in Zambia is largely due to the spread of resistant strains from inadequate treatment. Therefore, rapid detection of STR resistance genetically is recommended before its use in MDR-TB treatment in Zambia.

Structural and Genomic Insights Into Pyrazinamide Resistance in Mycobacterium tuberculosis Underlie Differences Between Ancient and Modern Lineages

Resistance to drugs used to treat tuberculosis disease (TB) continues to remain a public health burden, with missense point mutations in the underlying Mycobacterium tuberculosis bacteria described for nearly all anti-TB drugs. The post-genomics era along with advances in computational and structural biology provide opportunities to understand the interrelationships between the genetic basis and the structural consequences of M. tuberculosis mutations linked to drug resistance. Pyrazinamide (PZA) is a crucial first line antibiotic currently used in TB treatment regimens. The mutational promiscuity exhibited by the pncA gene (target for PZA) necessitates computational approaches to investigate the genetic and structural basis for PZA resistance development. We analysed 424 missense point mutations linked to PZA resistance derived from ∼35K M. tuberculosis clinical isolates sourced globally, which comprised the four main M. tuberculosis lineages (Lineage 1-4). Mutations were annotated to reflect their association with PZA resistance. Genomic measures (minor allele frequency and odds ratio), structural features (surface area, residue depth and hydrophobicity) and biophysical effects (change in stability and ligand affinity) of point mutations on pncA protein stability and ligand affinity were assessed. Missense point mutations within pncA were distributed throughout the gene, with the majority (>80%) of mutations with a destabilising effect on protomer stability and on ligand affinity. Active site residues involved in PZA binding were associated with multiple point mutations highlighting mutational diversity due to selection pressures at these functionally important sites. There were weak associations between genomic measures and biophysical effect of mutations. However, mutations associated with PZA resistance showed statistically significant differences between structural features (surface area and residue depth), but not hydrophobicity score for mutational sites. Most interestingly M. tuberculosis lineage 1 (ancient lineage) exhibited a distinct protein stability profile for mutations associated with PZA resistance, compared to modern lineages.

Phenotypic and genotypic features of the Mycobacterium tuberculosis lineage 1 subgroup in central Vietnam

Mycobacterium tuberculosis (Mtb) has different features depending on different geographic areas. We collected Mtb strains from patients with smear-positive pulmonary tuberculosis in Da Nang, central Vietnam. Using a whole genome sequencing platform, including genome assembly complemented by long-read-sequencing data, genomic characteristics were studied. Of 181 Mtb isolates, predominant Vietnamese EAI4_VNM and EAI4-like spoligotypes (31.5%), ZERO strains (5.0%), and part of EAI5 (11.1%) were included in a lineage-1 (L1) sublineage, i.e., L1.1.1.1. These strains were found less often in younger people, and they genetically clustered less frequently than other modern strains. Patients infected with ZERO strains demonstrated less lung infiltration. A region in RD2bcg spanning six loci, i.e., PE_PGRS35, cfp21, Rv1985c, Rv1986, Rv1987, and erm(37), was deleted in EAI4_VNM, EAI4-like, and ZERO strains, whereas another 118 bp deletion in furA was specific only to ZERO strains. L1.1.1.1-sublineage-specific deletions in PE_PGRS4 and PE_PGRS22 were also identified. RD900, seen in ancestral lineages, was present in majority of the L1 members. All strains without IS6110 (5.0%) had the ZERO spoligo-pattern. Distinctive features of the ancestral L1 strains provide a basis for investigation of the modern versus ancestral Mtb lineages and allow consideration of countermeasures against this heterogeneous pathogen.

Heterogeneous Streptomycin Resistance Level Among Mycobacterium tuberculosis Strains From the Same Transmission Cluster

Widespread and frequent resistance to the second-line tuberculosis (TB) medicine streptomycin, suggests ongoing transmission of low fitness cost streptomycin resistance mutations. To investigate this hypothesis, we studied a cohort of 681 individuals from a TB epidemic in Portugal. Whole-genome sequencing (WGS) analyses were combined with phenotypic growth studies in culture media and in mouse bone marrow derived macrophages. Streptomycin resistance was the most frequent resistance in the cohort accounting for 82.7% (n = 67) of the resistant Mycobacterium tuberculosis isolates. WGS of 149 clinical isolates identified 13 transmission clusters, including three clusters containing only streptomycin resistant isolates. The biggest cluster was formed by eight streptomycin resistant isolates with a maximum of five pairwise single nucleotide polymorphisms of difference. Interestingly, despite their genetic similarity, these isolates displayed different resistance levels to streptomycin, as measured both in culture media and in infected mouse bone marrow derived macrophages. The genetic bases underlying this phenotype are a combination of mutations in gid and other genes. This study suggests that specific streptomycin resistance mutations were transmitted in the cohort, with the resistant isolates evolving at the cluster level to allow low-to-high streptomycin resistance levels without a significative fitness cost. This is relevant not only to better understand transmission of streptomycin resistance in a clinical setting dominated by Lineage 4 M. tuberculosis infections, but mainly because it opens new prospects for the investigation of selection and spread of drug resistance in general.

Highly transmitted M. tuberculosis strains are more likely to evolve MDR/XDR and cause outbreaks, but what makes them highly transmitted?

Multi-Drug-Resistant strains of Mycobacterium tuberculosis (MDR-TB) are a serious obstacle to global TB eradication. While most MDR-TB strains are infrequently transmitted, a few cause large transmission clusters that contribute substantially to local MDR-TB burdens. Here we examine whether the known mutations in these strains can explain their success. Drug resistance mutations differ in fitness costs and strains can also acquire compensatory mutations (CM) to restore fitness, but some highly transmitted MDR strains have no CM. The acquisition of resistance mutations that maintain high transmissibility seems to occur by chance and are more likely in strains that are intrinsically highly transmitted and cause many cases. Modern Beijing lineage strains have caused several large outbreaks, but MDR outbreaks are also caused by ancient Beijing and lineage 4 strains, suggesting the lineage is less important than the characteristics of the individual strain. The development of fluoroquinolone resistance appears to represent another level of selection, in which strains must surmount unknown fitness costs of gyrA mutations. The genetic determinants of high transmission are poorly defined but may involve genes encoding proteins involved in molybdenum acquisition and the Esx systems. In addition, strains eliciting lower cytokine responses and producing more caseating granulomas may have advantages for transmission. Successful MDR/XDR strains generally evolve from highly transmitted drug sensitive parent strains due to selection pressures from deficiencies in local TB control programs. Until TB incidence is considerably reduced, there will likely be highly transmitted strains that develop resistance to any new antibiotic.

Compromised base excision repair pathway in Mycobacterium tuberculosis imparts superior adaptability in the host

Tuberculosis caused by Mycobacterium tuberculosis (Mtb) is a significant public health concern, exacerbated by the emergence of drug-resistant TB. To combat the host’s dynamic environment, Mtb encodes multiple DNA repair enzymes that play a critical role in maintaining genomic integrity. Mtb possesses a GC-rich genome, rendering it highly susceptible to cytosine deaminations, resulting in the occurrence of uracils in the DNA. UDGs encoded by ung and udgB initiate the repair; hence we investigated the biological impact of deleting UDGs in the adaptation of pathogen. We generated gene replacement mutants of uracil DNA glycosylases, individually (RvΔung, RvΔudgB) or together (RvΔdKO). The double KO mutant, RvΔdKO exhibited remarkably higher spontaneous mutation rate, in the presence of antibiotics. Interestingly, RvΔdKO showed higher survival rates in guinea pigs and accumulated large number of SNPs as revealed by whole-genome sequence analysis. Competition assays revealed the superior fitness of RvΔdKO over Rv, both in ex vivo and in vivo conditions. We propose that compromised DNA repair results in the accumulation of mutations, and a subset of these drives adaptation in the host. Importantly, this property allowed us to utilize RvΔdKO for the facile identification of drug targets.

Mutation in the genome of bacteria contributes to the acquisition of drug resistance. Mutations in bacteria can arise due to exposures to antibiotics, oxidative, reductive, and many other stresses that bacteria encounter in the host. Mtb has multiple DNA repair mechanisms, including a base excision repair pathway to restore the damaged genome. Here we set out to determine the impact of deleting the Uracil DNA base excision pathway on pathogen adaptability to both antibiotic and host induced stresses. Combinatorial mutant of Mtb UDGs showed higher spontaneous rates of mutations when subjected to antibiotic stress and showed higher survival levels in the guinea pig model of infection. Whole-genome sequence analysis showed significant accumulation of SNPs, suggesting that mutations providing survival advantage may have been positively selected. We also showed that double mutant of Mtb UDGs would be an excellent means to identify antibiotic targets in the bacteria. Competition experiments wherein we pitted wild type and double mutant against each other demonstrated that double mutant has a decisive edge over the wild type. Together, data suggest that the absence of a base excision repair pathway leads to higher mutations and provides a survival advantage under stress. They could be an invaluable tool for identifying targets of new antibiotics.

Early Drug Development and Evaluation of Putative Antitubercular Compounds in the -Omics Era

Tuberculosis (TB) is an infectious disease caused by the bacterium Mycobacterium tuberculosis. According to the WHO, the disease is one of the top 10 causes of death of people worldwide. Mycobacterium tuberculosis is an intracellular pathogen with an unusually thick, waxy cell wall and a complex life cycle. These factors, combined with M. tuberculosis ability to enter prolonged periods of latency, make the bacterium very difficult to eradicate. The standard treatment of TB requires 6-20months, depending on the drug susceptibility of the infecting strain. The need to take cocktails of antibiotics to treat tuberculosis effectively and the emergence of drug-resistant strains prompts the need to search for new antitubercular compounds. This review provides a perspective on how modern -omic technologies facilitate the drug discovery process for tuberculosis treatment. We discuss how methods of DNA and RNA sequencing, proteomics, and genetic manipulation of organisms increase our understanding of mechanisms of action of antibiotics and allow the evaluation of drugs. We explore the utility of mathematical modeling and modern computational analysis for the drug discovery process. Finally, we summarize how -omic technologies contribute to our understanding of the emergence of drug resistance.

Evaluation of drug susceptibility profile of Mycobacterium tuberculosis Lineage 1 from Brazil based on whole genome sequencing and phenotypic methods

BACKGROUND

The evaluation of procedures for drug susceptibility prediction of Mycobacterium tuberculosis based on genomic data against the conventional reference method test based on culture is realistic considering the scenario of growing number of tools proposals based on whole-genome sequences (WGS).

OBJECTIVES

This study aimed to evaluate drug susceptibility testing (DST) outcome based on WGS tools and the phenotypic methods performed on isolates of M. tuberculosis Lineage 1 from the state of Pará, Brazil, generally associated with low levels of drug resistance.

METHODOLOGY

Culture based DST was performed using the Proportion Method in Löwenstein-Jensen medium on 71 isolates that had been submitted to WGS. We analysed the seven main genome sequence-based tools for resistance and lineage prediction applied to M. tuberculosis and for comparison evaluation we have used the Kappa concordance test.

FINDINGS

When comparing the WGS-based tools against the DST, we observed the highest level of agreement using TB-profiler. Among the tools, TB-profiler, KvarQ and Mykrobe were those which identified the largest number of TB-MDR cases. Comparing the four most sensitive tools regarding resistance prediction, agreement was observed for 43 genomes.

MAIN CONCLUSIONS

Drug resistance profiling using next-generation sequencing offers rapid assessment of resistance-associated mutations, therefore facilitating rapid access to effective treatment.

Evaluation of the frequency of mutation genes in multidrug-resistant tuberculosis (MDR-TB) strains in Beijing, China

The aim of this study was to explore the frequency and distribution of gene mutations that are related to isoniazid (INH) and rifampin (RIF)-resistance in the strains of the multidrug-resistant tuberculosis (MDR-TB) Mycobacterium tuberculosis (M.tb) in Beijing, China. In this retrospective study, the genotypes of 173 MDR-TB strains were analysed by spoligotyping. The katG, inhA genes and the promoter region of inhA, in which genetic mutations confer INH resistance; and the rpoB gene, in which genetic mutations confer RIF resistance, were sequenced. The percentage of resistance-associated nucleotide alterations among the strains of different genotypes was also analysed. In total, 90.8% (157/173) of the MDR strains belonged to the Beijing genotype. Population characteristics were not significantly different among the strains of different genotypes. In total, 50.3% (87/173) strains had mutations at codon S315T of katG; 16.8% (29/173) of strains had mutations in the inhA promoter region; of them, 5.5% (15/173) had point mutations at -15 base (C→T) of the inhA promoter region. In total, 86.7% (150/173) strains had mutations at rpoB gene; of them, 40% (69/173) strains had mutations at codon S531L of rpoB. The frequency of mutations was not significantly higher in Beijing genotypic MDR strains than in non-Beijing genotypes. Beijing genotypic MDR-TB strains were spreading in Beijing and present a major challenge to TB control in this region. A high prevalence of katG Ser315Thr, inhA promoter region (-15C→T) and rpoB (S531L) mutations was observed. Molecular diagnostics based on gene mutations was a useful method for rapid detection of MDR-TB in Beijing, China.

Robust detection of point mutations involved in multidrug-resistant Mycobacterium tuberculosis in the presence of co-occurrent resistance markers

Tuberculosis disease is a major global public health concern and the growing prevalence of drug-resistant Mycobacterium tuberculosis is making disease control more difficult. However, the increasing application of whole-genome sequencing as a diagnostic tool is leading to the profiling of drug resistance to inform clinical practice and treatment decision making. Computational approaches for identifying established and novel resistance-conferring mutations in genomic data include genome-wide association study (GWAS) methodologies, tests for convergent evolution and machine learning techniques. These methods may be confounded by extensive co-occurrent resistance, where statistical models for a drug include unrelated mutations known to be causing resistance to other drugs. Here, we introduce a novel ‘cannibalistic’ elimination algorithm (“Hungry, Hungry SNPos”) that attempts to remove these co-occurrent resistant variants. Using an M. tuberculosis genomic dataset for the virulent Beijing strain-type (n = 3,574) with phenotypic resistance data across five drugs (isoniazid, rifampicin, ethambutol, pyrazinamide, and streptomycin), we demonstrate that this new approach is considerably more robust than traditional methods and detects resistance-associated variants too rare to be likely picked up by correlation-based techniques like GWAS.

Tuberculosis is one of the deadliest infectious diseases, being responsible for more than one million deaths per year. The causing bacteria are becoming increasingly drug-resistant, which is hampering disease control. At the same time, an unprecedented amount of bacterial whole-genome sequencing is increasingly informing clinical practice. In order to detect the genetic alterations responsible for developing drug resistance and predict resistance status from genomic data, bio-statistical methods and machine learning models have been employed. However, due to strongly overlapping drug resistance phenotypes and genotypes in multidrug-resistant datasets, the results of these correlation-based approaches frequently also contain mutations related to resistance against other drugs. In the past, this issue has often been ignored or partially resolved by either restricting the input data or in post-analysis screening—with both strategies relying on prior information. Here we present a heuristic algorithm for finding resistance-associated variants and demonstrate that it is considerably more robust towards co-occurrent resistance compared to traditional techniques. The software is available at https://github.com/julibeg/HHS.

Robust barcoding and identification of Mycobacterium tuberculosis lineages for epidemiological and clinical studies

BACKGROUND

Tuberculosis, caused by bacteria in the Mycobacterium tuberculosis complex (MTBC), is a major global public health burden. Strain-specific genomic diversity in the known lineages of MTBC is an important factor in pathogenesis that may affect virulence, transmissibility, host response and emergence of drug resistance. Fast and accurate tracking of MTBC strains is therefore crucial for infection control, and our previous work developed a 62-single nucleotide polymorphism (SNP) barcode to inform on the phylogenetic identity of 7 human lineages and 64 sub-lineages.

METHODS

To update this barcode, we analysed whole genome sequencing data from 35,298 MTBC isolates (~ 1 million SNPs) covering 9 main lineages and 3 similar animal-related species (M. tuberculosis var. bovis, M. tuberculosis var. caprae and M. tuberculosis var. orygis). The data was partitioned into training (N = 17,903, 50.7%) and test (N = 17,395, 49.3%) sets and were analysed using an integrated phylogenetic tree and population differentiation (FST) statistical approach.

RESULTS

By constructing a phylogenetic tree on the training MTBC isolates, we characterised 90 lineages or sub-lineages or species, of which 30 are new, and identified 421 robust barcoding mutations, of which a minimal set of 90 was selected that included 20 markers from the 62-SNP barcode. The barcoding SNPs (90 and 421) discriminated perfectly the 86 MTBC isolate (sub-)lineages in the test set and could accurately reconstruct the clades across the combined 35k samples.

CONCLUSIONS

The validated 90 SNPs can be used for the rapid diagnosis and tracking of MTBC strains to assist public health surveillance and control. To facilitate this, the SNP markers have now been incorporated into the TB-Profiler informatics platform ( https://github.com/jodyphelan/TBProfiler ).

Population structure of multidrug-resistant Mycobacterium tuberculosis clinical isolates in Colombia

Emergence of multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) isolates is a major public health problem that threatens progress made in tuberculosis (TB) care and control worldwide. In Colombia, the prevalence of MDR tuberculosis (MDR-TB) has increased slowly but steadily since 2001. However, the population structure of the MDR-TB strains circulating in Colombia is sparsely known. In this work, 203 MDR isolates isolated in 2012-2013 were collected, and characterized by spoligotyping, followed by 24-loci MIRU-VNTR (data available for 190 isolates). The most prevalent genotypes corresponded to SIT42/LAM9 (12.81%), SIT62/H1 (10.34%), and SIT190/Beijing (10.34%). A fine analysis showed that although the MDR strains came from 29 of the 33 departments of Colombia, the distribution of these main lineages was not at random and depended on the city of isolation (p-value <0.000001). Both LAM and Beijing lineage strains were significantly associated with MDR-TB (p-value <0.0001): LAM lineage was associated with 2 patterns of MDR, namely combined resistance to INH + Rifampin (HR), and to SHRE (Streptomycin + INH + Rifampin + Ethambutol), while the Beijing lineage strains were essentially associated with MDR (SHRE). Interestingly, distribution of genotypic lineages in function of drug resistance information (e.g. pansusceptible vs. MDR) was different in our setting as compared to other countries in Latin America. However, MIRU-VNTR patterns were unique for all strains, an observation that did not support active transmission of circulating MDR clones.

Combining structure and genomics to understand antimicrobial resistance

Antimicrobials against bacterial, viral and parasitic pathogens have transformed human and animal health. Nevertheless, their widespread use (and misuse) has led to the emergence of antimicrobial resistance (AMR) which poses a potentially catastrophic threat to public health and animal husbandry. There are several routes, both intrinsic and acquired, by which AMR can develop. One major route is through non-synonymous single nucleotide polymorphisms (nsSNPs) in coding regions. Large scale genomic studies using high-throughput sequencing data have provided powerful new ways to rapidly detect and respond to such genetic mutations linked to AMR. However, these studies are limited in their mechanistic insight. Computational tools can rapidly and inexpensively evaluate the effect of mutations on protein function and evolution. Subsequent insights can then inform experimental studies, and direct existing or new computational methods. Here we review a range of sequence and structure-based computational tools, focussing on tools successfully used to investigate mutational effect on drug targets in clinically important pathogens, particularly Mycobacterium tuberculosis. Combining genomic results with the biophysical effects of mutations can help reveal the molecular basis and consequences of resistance development. Furthermore, we summarise how the application of such a mechanistic understanding of drug resistance can be applied to limit the impact of AMR.

Challenges of Immune Response Diversity in the Human Population Concerning New Tuberculosis Diagnostics, Therapies, and Vaccines

Universal approaches to the prevention and treatment of human diseases fail to take into account profound immune diversity resulting from genetic variations across populations. Personalized or precision medicine takes into account individual lifestyle, environment, and biology (genetics and immune status) and is being adopted in several disease intervention strategies such as cancer and heart disease. However, its application in infectious diseases, particularly global diseases such as tuberculosis (TB), is far more complex and in a state of infancy. Here, we discuss the impact of human genetic variations on immune responses and how they relate to failures seen in current TB diagnostic, therapy, and vaccine approaches across populations. We offer our perspective on the challenges and potential for more refined approaches going forward.

The Genetic Background Modulates the Evolution of Fluoroquinolone-Resistance in Mycobacterium tuberculosis

Fluoroquinolones (FQ) form the backbone in experimental treatment regimens against drug-susceptible tuberculosis. However, little is known on whether the genetic variation present in natural populations of Mycobacterium tuberculosis (Mtb) affects the evolution of FQ-resistance (FQ-R). To investigate this question, we used nine genetically distinct drug-susceptible clinical isolates of Mtb and measured their frequency of resistance to the FQ ofloxacin (OFX) in vitro. We found that the Mtb genetic background led to differences in the frequency of OFX-resistance (OFX-R) that spanned two orders of magnitude and substantially modulated the observed mutational profiles for OFX-R. Further, in vitro assays showed that the genetic background also influenced the minimum inhibitory concentration and the fitness effect conferred by a given OFX-R mutation. To test the clinical relevance of our in vitro work, we surveyed the mutational profile for FQ-R in publicly available genomic sequences from clinical Mtb isolates, and found substantial Mtb lineage-dependent variability. Comparison of the clinical and the in vitro mutational profiles for FQ-R showed that 51% and 39% of the variability in the clinical frequency of FQ-R gyrA mutation events in Lineage 2 and Lineage 4 strains, respectively, can be attributed to how Mtb evolves FQ-R in vitro. As the Mtb genetic background strongly influenced the evolution of FQ-R in vitro, we conclude that the genetic background of Mtb also impacts the evolution of FQ-R in the clinic.

Resistance Sniffer: An online tool for prediction of drug resistance patterns of Mycobacterium tuberculosis isolates using next generation sequencing data

The effective control of multidrug resistant tuberculosis (MDR-TB) relies upon the timely diagnosis and correct treatment of all tuberculosis cases. Whole genome sequencing (WGS) has great potential as a method for the rapid diagnosis of drug resistant Mycobacterium tuberculosis (Mtb) isolates. This method overcomes most of the problems that are associated with current phenotypic drug susceptibility testing. However, the application of WGS in the clinical setting has been deterred by data complexities and skill requirements for implementing the technologies as well as clinical interpretation of the next generation sequencing (NGS) data. The proposed diagnostic application was drawn upon recent discoveries of patterns of Mtb clade-specific genetic polymorphisms associated with antibiotic resistance. A catalogue of genetic determinants of resistance to thirteen anti-TB drugs for each phylogenetic clade was created. A computational algorithm for the identification of states of diagnostic polymorphisms was implemented as an online software tool, Resistance Sniffer (http://resistance-sniffer.bi.up.ac.za/), and as a stand-alone software tool to predict drug resistance in Mtb isolates using complete or partial genome datasets in different file formats including raw Illumina fastq read files. The program was validated on sequenced Mtb isolates with data on antibiotic resistance trials available from GMTV database and from the TB Platform of South African Medical Research Council (SAMRC), Pretoria. The program proved to be suitable for probabilistic prediction of drug resistance profiles of individual strains and large sequence data sets.