Mycobacterium tuberculosis Whole Genome Sequences From Southern India Suggest Novel Resistance Mechanisms and the Need for Region-Specific Diagnostics

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.

Whole-genome sequencing of clinical isolates from tuberculosis patients in India: real-world data indicates a high proportion of pre-XDR cases

Treatment decisions for tuberculosis (TB) in the absence of full drug-susceptibility data can result in amplifying resistance and may compromise treatment outcomes. Genomics of Mycobacterium tuberculosis (M.tb) from clinical samples enables detection of drug resistance to multiple drugs. We performed whole-genome sequencing (WGS) for 600 clinical samples from patients with tuberculosis to identify the drug-resistance profile and mutation spectrum. We documented the reasons reported by clinicians for referral. WGS identified a high proportion (51%) of pre-extensively drug-resistant (pre-XDR) cases followed by multidrug-resistant tuberculosis (MDR-TB) (15.5%). This correlates with the primary reason for referral, as non-response to the first-line treatment (67%) and treatment failure or rifampicin resistance (14%). Multivariate analysis indicated that all young age groups (P < 0.05), male gender (P < 0.05), and Beijing strain (P < 0.01) were significant independent predictors of MDR-TB or MDR-TB+ [pre-extensively drug-resistant tuberculosis (XDR-TB) and XDR-TB]. Ser315Thr (72.5%) in the inhA gene and Ser450Leu in the rpoB gene (65.5%) were the most prevalent mutations, as were resistance-conferring mutations to pyrazinamide (41%) and streptomycin (61.33%). Mutations outside the rifampicin resistance-determining region (RRDR), Ile491Phe and Val170Phe, were seen in 1.3% of cases; disputed mutations in rpoB (Asp435Tyr, His445Asn, His445Leu, and Leu430Pro) were seen in 6% of cases, and mutations to newer drugs such as bedaquiline and linezolid in 1.0% and 7.5% of cases, respectively. This study on clinical samples highlights that there is a high proportion of pre-XDR cases and emerging resistance to newer drugs; ongoing transmission of these strains can cause serious threat to public health; and whole-genome sequencing can effectively identify and support precision medicine for TB.

IMPORTANCE

The current study is based on real-world data on the TB drug-resistance profile by whole-genome sequencing of 600 clinical samples from patients with TB in India. This study indicates the clinicians' reasons for sending samples for WGS, which is for difficult-to-treat cases and/or relapse and treatment failure. The study reports a significant proportion of cases with pre-XDR-TB strains that warrant policy makers' attention. It reflects the current iterative nature of the diagnostic tests under programmatic conditions that leads to delays in appropriate diagnosis and empirical treatment. India had an estimated burden of 2.95 million TB cases in 2020 and 135,000 multidrug-resistant cases. However, WGS profiles of M.tb from India remains disproportionately poorly represented. This study adds a significant body of data on the mutation profiles seen in M.tb isolated from patients with TB in India, mutations outside the RRDR, disputed mutations, and resistance-conferring mutations to newer drugs such as bedaquiline and linezolid.

Actionable mechanisms of drug tolerance and resistance in Mycobacterium tuberculosis

The emergence of antimicrobial resistance (AMR) across bacterial pathogens presents a serious threat to global health. This threat is further exacerbated in tuberculosis (TB), mainly due to a protracted treatment regimen involving a combination of drugs. A diversity of factors contributes to the emergence of drug resistance in TB, which is caused by the pathogen Mycobacterium tuberculosis (Mtb). While the traditional genetic mutation-driven drug resistance mechanisms operate in Mtb, there are also several additional unique features of drug resistance in this pathogen. Research in the past decade has enriched our understanding of such unconventional factors as efflux pumps, bacterial heterogeneity, metabolic states, and host microenvironment. Given that the discovery of new antibiotics is outpaced by the emergence of drug resistance patterns displayed by the pathogen, newer strategies for combating drug resistance are desperately needed. In the context of TB, such approaches include targeting the efflux capability of the pathogen, modulating the host environment to prevent bacterial drug tolerance, and activating the host anti-mycobacterial pathways. In this review, we discuss the traditional mechanisms of drug resistance in Mtb, newer understandings and the shaping of a set of unconventional approaches to target both the emergence and treatment of drug resistance in TB.

Role of the first WHO mutation catalogue in the diagnosis of antibiotic resistance in Mycobacterium tuberculosis in the Valencia Region, Spain: a retrospective genomic analysis

BACKGROUND

In June, 2021, WHO published the most complete catalogue to date of resistance-conferring mutations in Mycobacterium tuberculosis. Here, we aimed to assess the performance of genome-based antimicrobial resistance prediction using the catalogue and its potential for improving diagnostics in a real low-burden setting.

METHODS

In this retrospective population-based genomic study M tuberculosis isolates were collected from 25 clinical laboratories in the low-burden setting of the Valencia Region, Spain. Culture-positive tuberculosis cases reported by regional public health authorities between Jan 1, 2014, and Dec 31, 2016, were included. The drug resistance profiles of these isolates were predicted by the genomic identification, via whole-genome sequencing (WGS), of the high-confidence resistance-causing variants included in the catalogue and compared with the phenotype. We determined the minimum inhibitory concentration (MIC) of the isolates with discordant resistance profiles using the resazurin microtitre assay.

FINDINGS

WGS was performed on 785 M tuberculosis complex culture-positive isolates, and the WGS resistance prediction sensitivities were: 85·4% (95% CI 70·8-94·4) for isoniazid, 73·3% (44·9-92·2) for rifampicin, 50·0% (21·1-78·9) for ethambutol, and 57·1% (34·0-78·2) for pyrazinamide; all specificities were more than 99·6%. Sensitivity values were lower than previously reported, but the overall pan-susceptibility accuracy was 96·4%. Genotypic analysis revealed that four phenotypically susceptible isolates carried mutations (rpoB Leu430Pro and rpoB Ile491Phe for rifampicin and fabG1 Leu203Leu for isoniazid) known to give borderline resistance in standard phenotypic tests. Additionally, we identified three putative resistance-associated mutations (inhA Ser94Ala, katG Leu48Pro, and katG Gly273Arg for isoniazid) in samples with substantially higher MICs than those of susceptible isolates. Combining both genomic and phenotypic data, in accordance with the WHO diagnostic guidelines, we could detect two new multidrug-resistant cases. Additionally, we detected 11 (1·6%) of 706 isolates to be monoresistant to fluoroquinolone, which had been previously undetected.

INTERPRETATION

We showed that the WHO catalogue enables the detection of resistant cases missed in phenotypic testing in a low-burden region, thus allowing for better patient-tailored treatment. We also identified mutations not included in the catalogue, relevant at the local level. Evidence from this study, together with future updates of the catalogue, will probably lead in the future to the partial replacement of culture testing with WGS-based drug susceptibility testing in our setting.

FUNDING

European Research Council and the Spanish Ministerio de Ciencia.

The effect of M. tuberculosis lineage on clinical phenotype

Six lineages of Mycobacterium tuberculosis sensu stricto (which excludes M. africanum) are described. Single-country or small observational data suggest differences in clinical phenotype between lineages. We present strain lineage and clinical phenotype data from 12,246 patients from 3 low-incidence and 5 high-incidence countries. We used multivariable logistic regression to explore the effect of lineage on site of disease and on cavities on chest radiography, given pulmonary TB; multivariable multinomial logistic regression to investigate types of extra-pulmonary TB, given lineage; and accelerated failure time and Cox proportional-hazards models to explore the effect of lineage on time to smear and culture-conversion. Mediation analyses quantified the direct effects of lineage on outcomes. Pulmonary disease was more likely among patients with lineage(L) 2, L3 or L4, than L1 (adjusted odds ratio (aOR) 1.79, (95% confidence interval 1.49-2.15), p<0.001; aOR = 1.40(1.09-1.79), p = 0.007; aOR = 2.04(1.65-2.53), p<0.001, respectively). Among patients with pulmonary TB, those with L1 had greater risk of cavities on chest radiography versus those with L2 (aOR = 0.69(0.57-0.83), p<0.001) and L4 strains (aOR = 0.73(0.59-0.90), p = 0.002). L1 strains were more likely to cause osteomyelitis among patients with extra-pulmonary TB, versus L2-4 (p = 0.033, p = 0.008 and p = 0.049 respectively). Patients with L1 strains showed shorter time-to-sputum smear conversion than for L2. Causal mediation analysis showed the effect of lineage in each case was largely direct. The pattern of clinical phenotypes seen with L1 strains differed from modern lineages (L2-4). This has implications for clinical management and could influence clinical trial selection strategies.

Lineage classification and antitubercular drug resistance surveillance of Mycobacterium tuberculosis by whole-genome sequencing in Southern India

Whole-genome sequencing has created a revolution in tuberculosis management by providing a comprehensive picture of the various genetic polymorphisms with unprecedented accuracy. Studies mapping genomic heterogeneity in clinical isolates of Mycobacterium tuberculosis using a whole-genome sequencing approach from high tuberculosis burden countries are underrepresented. We report whole-genome sequencing results of 242 clinical isolates of culture-confirmed M. tuberculosis isolates from tuberculosis patients referred to a tertiary care hospital in Southern India. Phylogenetic analysis revealed that the isolates in our study belonged to five different lineages, with Indo-Oceanic (lineage 1, n = 122) and East-African Indian (lineage 3, n = 80) being the most prevalent. We report several mutations in genes conferring resistance to first and second line antitubercular drugs including the genes rpoB, katG, ahpC, inhA, fabG1, embB, pncA, rpsL, rrs, and gyrA. The majority of these mutations were identified in relatively high proportions in lineage 1. Our study highlights the utility of whole-genome sequencing as a potential supplemental tool to the existing genotypic and phenotypic methods, in providing expedited comprehensive surveillance of mutations that may be associated with antitubercular drug resistance as well as lineage characterization of M. tuberculosis isolates. Further larger-scale whole-genome datasets with linked minimum inhibition concentration testing are imperative for resolving the discrepancies between whole-genome sequencing and phenotypic drug sensitivity testing results and quantifying the level of the resistance associated with the mutations for optimization of antitubercular drug and precise dose selection in clinics. IMPORTANCE Studies mapping genetic heterogeneity of clinical isolates of M. tuberculosis for determining their strain lineage and drug resistance by whole-genome sequencing are limited in high tuberculosis burden settings. We carried out whole-genome sequencing of 242 M. tuberculosis isolates from drug-sensitive and drug-resistant tuberculosis patients, identified and collected as part of the TB Portals Program, to have a comprehensive insight into the genetic diversity of M. tuberculosis in Southern India. We report several genetic variations in M. tuberculosis that may confer resistance to antitubercular drugs. Further wide-scale efforts are required to fully characterize M. tuberculosis genetic diversity at a population level in high tuberculosis burden settings for providing precise tuberculosis treatment.

Whole-Genome Sequencing of Mycobacterium tuberculosis Isolates from Diabetic and Non-Diabetic Patients with Pulmonary Tuberculosis

The association of tuberculosis and type 2 diabetes mellitus has been a recognized re-emerging challenge in management of the convergence of the two epidemics. Though much of the literature has studied this association, there is less knowledge in the field of genetic diversities that might occur in strains infecting tuberculosis patients with and without diabetes. Our study focused on determining the extent of diversity of genotypes of Mycobacterium tuberculosis in both these categories of patients. We subjected 55 M. tuberculosis isolates from patients diagnosed with pulmonary TB with and without type 2 diabetes mellitus to whole-genome sequencing on Illumina Hi Seq platform. The most common lineage identified was lineage 1, the Indo-Oceanic lineage (n = 22%), followed by lineage 4, the Euro-American lineage (n = 18, 33%); lineage 3, the East-African Indian lineage (n = 13, 24%); and lineage 2, the East-Asian lineage (n = 1, 2%). There were no significant differences in the distribution of lineages in both diabetics and non-diabetics in the South Indian population, and further studies involving computational analysis and comparative transcriptomics are needed to provide deeper insights.

The effect of M. tuberculosis lineage on clinical phenotype

Eight lineages of Mycobacterium tuberculosis sensu stricto are described. Single-country or small observational data suggest differences in clinical phenotype between lineages. We present strain lineage and clinical phenotype data from 12,246 patients from 3 low-incidence and 5 high-incidence countries. We used multivariable logistic regression to explore the effect of lineage on site of disease and on cavities on chest radiography, given pulmonary TB; multivariable multinomial logistic regression to investigate types of extra-pulmonary TB, given lineage; and accelerated failure time and Cox proportional-hazards models to explore the effect of lineage on time to smear and culture-conversion. Mediation analyses quantified the direct effects of lineage on outcomes. Pulmonary disease was more likely among patients with lineage(L) 2, L3 or L4, than L1 (adjusted odds ratio (aOR) 1.79, (95% confidence interval 1.49-2.15), p<0.001; aOR=1.40(1.09-1.79), p=0.007; aOR=2.04(1.65-2.53), p<0.001, respectively). Among patients with pulmonary TB, those with L1 had greater risk of cavities on chest radiography versus those with L2 (aOR=0.69(0.57-0.83), p<0.001) and L4 strains (aOR=0.73(0.59-0.90), p=0.002). L1 strains were more likely to cause osteomyelitis among patients with extra-pulmonary TB, versus L2-4 (p=0.033, p=0.008 and p=0.049 respectively). Patients with L1 strains showed shorter time-to-sputum smear conversion than for L2. Causal mediation analysis showed the effect of lineage in each case was largely direct. The pattern of clinical phenotypes seen with L1 strains differed from modern lineages (L2-4). This has implications for clinical management and could influence clinical trial selection strategies.

Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia

Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.

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.

High fluoroquinolone resistance proportions among multidrug-resistant tuberculosis driven by dominant L2 Mycobacterium tuberculosis clones in the Mumbai Metropolitan Region

Background

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains are a serious health problem in India, also contributing to one-fourth of the global MDR tuberculosis (TB) burden. About 36% of the MDR MTBC strains are reported fluoroquinolone (FQ) resistant leading to high pre-extensively drug-resistant (pre-XDR) and XDR-TB (further resistance against bedaquiline and/or linezolid) rates. Still, factors driving the MDR/pre-XDR epidemic in India are not well defined.

Methods

In a retrospective study, we analyzed 1852 consecutive MTBC strains obtained from patients from a tertiary care hospital laboratory in Mumbai by whole genome sequencing (WGS). Univariate and multivariate statistics was used to investigate factors associated with pre-XDR. Core genome multi locus sequence typing, time scaled haplotypic density (THD) method and homoplasy analysis were used to analyze epidemiological success, and positive selection in different strain groups, respectively.

Results

In total, 1016 MTBC strains were MDR, out of which 703 (69.2%) were pre-XDR and 45 (4.4%) were XDR. Cluster rates were high among MDR (57.8%) and pre-XDR/XDR (79%) strains with three dominant L2 (Beijing) strain clusters (Cl 1–3) representing half of the pre-XDR and 40% of the XDR-TB cases. L2 strains were associated with pre-XDR/XDR-TB (P < 0.001) and, particularly Cl 1–3 strains, had high first-line and FQ resistance rates (81.6–90.6%). Epidemic success analysis using THD showed that L2 strains outperformed L1, L3, and L4 strains in short- and long-term time scales. More importantly, L2 MDR and MDR + strains had higher THD success indices than their not-MDR counterparts. Overall, compensatory mutation rates were highest in L2 strains and positive selection was detected in genes of L2 strains associated with drug tolerance (prpB and ppsA) and virulence (Rv2828c). Compensatory mutations in L2 strains were associated with a threefold increase of THD indices, suggesting improved transmissibility.

Conclusions

Our data indicate a drastic increase of FQ resistance, as well as emerging bedaquiline resistance which endangers the success of newly endorsed MDR-TB treatment regimens. Rapid changes in treatment and control strategies are required to contain transmission of highly successful pre-XDR L2 strains in the Mumbai Metropolitan region but presumably also India-wide.

Mycobacterium tuberculosis Lineages Associated with Mutations and Drug Resistance in Isolates from India

Current knowledge on resistance-conferring determinants in Mycobacterium tuberculosis is biased toward globally dominant lineages 2 and 4. In contrast, lineages 1 and 3 are predominant in India. In this study, we performed whole-genome sequencing of 498 MDR M. tuberculosis isolates from India to determine the prevalence of drug resistance mutations and to understand the genomic diversity. A retrospective collection of 498 M. tuberculosis isolates submitted to the National Institute for Research in Tuberculosis for phenotypic susceptibility testing between 2014 to 2016 were sequenced. Genotypic resistance prediction was performed using known resistance-conferring determinants. Genotypic and phenotypic results for 12 antituberculosis drugs were compared, and sequence data were explored to characterize lineages and their association with drug resistance. Four lineages were identified although lineage 1 predominated (43%). The sensitivity of prediction for isoniazid and rifampicin was 92% and 98%, respectively. We observed lineage-specific variations in the proportion of isolates with resistance-conferring mutations, with drug resistance more common in lineages 2 and 3. Disputed mutations (codons 430, 435, 445, and 452) in the rpoB gene were more common in isolates other than lineage 2. Phylogenetic analysis and pairwise SNP difference revealed high genetic relatedness of lineage 2 isolates. WGS based resistance prediction has huge potential, but knowledge of regional and national diversity is essential to achieve high accuracy for resistance prediction. IMPORTANCE Current knowledge on resistance-conferring determinants in Mycobacterium tuberculosis is biased toward globally dominant lineages 2 and 4. In contrast, lineages 1 and 3 are predominant in India. We performed whole-genome sequencing of 498 MDR M. tuberculosis isolates from India to determine the prevalence of drug resistance mutations and to understand genomic diversity. Four lineages were identified although lineage 1 predominated (43%). The sensitivity of prediction for isoniazid and rifampicin was 92% and 98%, respectively. We observed lineage-specific variations in the proportion of isolates with resistance-conferring mutations, with drug resistance more common in lineages 2 and 3. Disputed mutations (codons 430, 435, 445, and 452) in the rpoB gene were more common in isolates other than lineage 2. Phylogenetic analysis and pairwise SNP difference revealed high genetic relatedness of lineage 2 isolates. WGS based resistance prediction has huge potential, but knowledge of regional and national diversity is essential to achieve high accuracy for resistance prediction.

Whole-genome sequencing of Mycobacterium tuberculosis from Cambodia

Cambodia has one of the highest tuberculosis (TB) incidence rates in the WHO Western Pacific region. Remarkably though, the prevalence of multidrug-resistant TB (MDR-TB) remains low. We explored the genetic diversity of Mycobacterium tuberculosis (MTB) circulating in this unique setting using whole-genome sequencing (WGS). From October 2017 until January 2018, we collected one hundred sputum specimens from consenting adults older than 21 years of age, newly diagnosed with bacteriologically confirmed TB in 3 districts of Phnom Penh and Takeo provinces of Cambodia before they commence on their TB treatment, where eighty MTB isolates were successfully cultured and sequenced. Majority of the isolates belonged to Lineage 1 (Indo-Oceanic) (69/80, 86.25%), followed by Lineage 2 (East Asian) (10/80, 12.5%) and Lineage 4 (Euro-American) (1/80, 1.25%). Phenotypic resistance to both streptomycin and isoniazid was found in 3 isolates (3/80, 3.75%), while mono-resistance to streptomycin and isoniazid was identical at 2.5% (N = 2 each). None of the isolates tested was resistant to either rifampicin or ethambutol. The specificities of genotypic prediction for resistance to all drugs tested were 100%, while the sensitivities of genotypic resistance predictions to isoniazid and streptomycin were lower at 40% (2/5) and 80% (4/5) respectively. We identified 8 clusters each comprising of two to five individuals all residing in the Takeo province, making up half (28/56, 50%) of all individuals sampled in the province, indicating the presence of multiple ongoing transmission events. All clustered isolates were of Lineage 1 and none are resistant to any of the drugs tested. This study while demonstrating the relevance and utility of WGS in predicting drug resistance and inference of disease transmission, highlights the need to increase the representation of genotype-phenotype TB data from low and middle income countries in Asia and Africa to improve the accuracies for prediction of drug resistance.

Distribution and identification of Mycobacterium tuberculosis lineage in Kashgar prefecture

Objectives

Kashgar prefecture is an important transportation and trade hub with a high incidence of tuberculosis. The following study analyzed the composition and differences in Mycobacterium tuberculosis (M.tb) lineage and specific tags to distinguish the lineage of the M.tb in Kashgar prefecture, thus providing a basis for the classification and diagnosis of tuberculosis in this area.

Methods

Whole-genome sequencing (WGS) of 161 M.tb clinical strains was performed. The phylogenetic tree was constructed using Maximum Likelihood (ML) based on single nucleotide polymorphisms (SNPs) and verified through principal component analysis (PCA). The composition structure of M.tb in different regions was analyzed by combining geographic information.

Results

M.tb clinical strains were composed of lineage 2 (73/161, 45.34%), lineage 3 (52/161, 32.30%) and lineage 4 (36/161, 22.36%). Moreover, the 3 lineages were subdivided into 11 sublineages, among which lineage 2 included lineage 2.2.2/Asia Ancestral 1 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 2 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 3 (18/73, 24.66%), and lineage 2.2.1-Modern Beijing (39/73, 53.42%). Lineage 3 included lineage 3.2 (14/52, 26.92%) and lineage 3.3 (38/52, 73.08%), while lineage 4 included lineage 4.1 (3/36, 8.33%), lineage 4.2 (2/36, 5.66%), lineage 4.4.2 (1/36, 2.78%), lineage 4.5 (28/36, 77.78%) and lineage 4.8 (2/36, 5.66%), all of which were consistent with the PCA results. One hundred thirty-six markers were proposed for discriminating known circulating strains. Reconstruction of a phylogenetic tree using the 136 SNPs resulted in a tree with the same number of delineated clades. Based on geographical location analysis, the composition of Lineage 2 in Kashgar prefecture (45.34%) was lower compared to other regions in China (54.35%-90.27%), while the composition of Lineage 3 (32.30%) was much higher than in other regions of China (0.92%-2.01%), but lower compared to the bordering Pakistan (70.40%).

Conclusion

Three lineages were identified in M.tb clinical strains from Kashgar prefecture, with 136 branch-specific SNP. Kashgar borders with countries that have a high incidence of tuberculosis, such as Pakistan and India, which results in a large difference between the M.tb lineage and sublineage distribution in this region and other provinces of China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07307-4.

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.

Population structure, biogeography and transmissibility of Mycobacterium tuberculosis

Mycobacterium tuberculosis is a clonal pathogen proposed to have co-evolved with its human host for millennia, yet our understanding of its genomic diversity and biogeography remains incomplete. Here we use a combination of phylogenetics and dimensionality reduction to reevaluate the population structure of M. tuberculosis, providing an in-depth analysis of the ancient Indo-Oceanic Lineage 1 and the modern Central Asian Lineage 3, and expanding our understanding of Lineages 2 and 4. We assess sub-lineages using genomic sequences from 4939 pan-susceptible strains, and find 30 new genetically distinct clades that we validate in a dataset of 4645 independent isolates. We find a consistent geographically restricted or unrestricted pattern for 20 groups, including three groups of Lineage 1. The distribution of terminal branch lengths across the M. tuberculosis phylogeny supports the hypothesis of a higher transmissibility of Lineages 2 and 4, in comparison with Lineages 3 and 1, on a global scale. We define an expanded barcode of 95 single nucleotide substitutions that allows rapid identification of 69 M. tuberculosis sub-lineages and 26 additional internal groups. Our results paint a higher resolution picture of the M. tuberculosis phylogeny and biogeography.

A Low-Prevalence Single-Nucleotide Polymorphism in the Sensor Kinase PhoR in Mycobacterium tuberculosis Suppresses Its Autophosphatase Activity and Reduces Pathogenic Fitness: Implications in Evolutionary Selection

The genome sequencing of Mycobacterium tuberculosis, the causative organism of tuberculosis, has significantly improved our understanding of the mechanisms that drive the establishment of infection and disease progression. Several clinical strains of M. tuberculosis exhibit single-nucleotide polymorphisms (SNPs), the implications of which are only beginning to be understood. Here, we examined the impact of a specific polymorphism in PhoR, the sensor kinase of the PhoPR two-component system. Biochemical analysis revealed reduced autophosphatase/ATPase activity, which led to enhanced downstream gene expression. We complemented M. tuberculosis H37Ra with the wild-type and mutant phoPR genes and characterized the strains in a cell line infection model. We provide an explanation for the low prevalence of the SNP in clinical strains (∼1%), as the mutation causes a survival disadvantage in the host cells. The study provides a rare example of selection of a signaling node under competing evolutionary forces, wherein a biochemically superior mutation aids bacterial adaptation within-host but has low fitness for infection and hence is not selected. Our study highlights the importance of accounting for such SNPs to test therapeutic and co-therapeutic methods to combat TB.

Characterization of genetic diversity and clonal complexes by whole genome sequencing of Mycobacterium tuberculosis isolates from Jalisco, Mexico

Whole genome sequencing (WGS) analysis in tuberculosis allows the prediction of drug-resistant phenotypes, identification of lineages, and to better understanding of the epidemiology and transmission chains. Nevertheless the procedure has been scarcely assessed in Mexico, in this work we analyze by WGS isolates of Mycobacterium tuberculosis circulating in Jalisco, Mexico. Lineage and phylogenetic characterization, drug resistant prediction, "in silico" spoligotyping determination, were provided by WGS in 32 M. tuberculosis clinical isolates. Lineage 4 (L4), with 28 isolates (87%) and eleven sublineages was dominant. Forty SNPs and INDELs were found in genes related to first-, and second-line drugs. Eleven isolates were sensitive, seven (22%) were predicted to be resistant to isoniazid, two resistant to rifampicin (6%) and two (6%) were multidrug-resistant tuberuclosis. Spoligotyping shows that SIT 53 (19%) and SIT 119 (16%) were dominant. Four clonal transmission complexes were found. This is the first molecular epidemiological description of TB isolates circulating in western Mexico, achieved through WGS. L4 was dominant and included a high diversity of sublineages. It was possible to track the transmission route of two clonal complexes. The WGS demonstrated to be of great utility and with further implications for clinical and epidemiological study of TB in the region.

Evaluation of the Roche cobas MTB and MTB-RIF/INH Assays in Samples from Germany and Sierra Leone

The Roche cobas MTB and MTB-RIF/INH assays allow for detection of Mycobacterium tuberculosis complex (MTBC) nucleic acid and rifampicin (RIF) and isoniazid (INH) resistance-associated mutations in an automated, high-throughput workflow. In this study, we evaluated the performance of these assays, employing samples from settings of low and high tuberculosis (TB) burdens. A total of 325 frozen, leftover respiratory samples collected from treatment-naive patients with presumptive TB in Germany (n = 280) and presumptive RIF-resistant TB in Sierra Leone (n = 45) were used in this study. cobas MTB results for detection of MTBC DNA from N-acetyl-l-cysteine-sodium hydroxide (NALC-NaOH)-treated samples were compared to culture results. Predictions of RIF and INH resistance by the cobas MTB-RIF/INH assay were compared to a composite reference standard (phenotypic drug susceptibility testing and line probe assay). Whole-genome sequencing was used to resolve discordances. The overall sensitivity of cobas MTB for detection of MTBC DNA in culture-positive samples (n = 102) was 89.2% (95% confidence interval [CI], 81.7 to 93.9%). The specificity of cobas MTB was 98.6% (95% CI, 96.1 to 99.5%). Sensitivity and specificity for detection of RIF and INH resistance were 88.4% (95% CI, 75.5 to 94.9%) and 97.6% (95% CI, 87.4 to 99.6%) and 76.6% (95% CI, 62.8 to 86.4%) and 100.0% (95% CI, 90.8 to 100.0%), respectively. Discordant results for RIF and INH resistance were mainly due to uncommon mutations in samples from Sierra Leone that were not covered by the cobas MTB-RIF/INH assay. In conclusion, cobas MTB and MTB-RIF/INH assays provide accurate detection of MTBC DNA and resistance-associated mutations in respiratory samples. The influence of regional variations in the prevalence of resistance-conferring mutations requires further investigation.

Estimation of the global burden of Mycobacterium tuberculosis lineage 1

Tuberculosis is still problematic as it affects large numbers of people globally. Mycobacterium tuberculosis Lineage 1 (L1) or Indo Oceanic Lineage, one of widespread major lineages, has a specific geographic distribution and high mortality. It is highly diverse and endemic in several high burden countries. However, studies on the global burden of L1 and its sublineages remain limited. This may lead to the underestimation of the importance of its variance in developing and applying tuberculosis control measures. This study aimed to estimate the number of patients infected with M. tuberculosis L1 and its sublineages worldwide. The proportion of L1 among tuberculosis patients was searched in published reports from countries around the world and the number of patients was calculated based on a WHO report on country incidences and populations. The numbers of patients infected with the five major sublineages, namely L1.1.1, L1.1.2, L1.1.3, L1.2.1, and L1.2.2 were estimated where information was available. It was found that L1 accounted for 28% of global tuberculosis cases in 2012 and 2018. Over 80% of the L1 global burden was in India, the Philippines, Indonesia and Bangladesh, which are also among the countries with highest absolute numbers of tuberculosis patients in the world. Globally, the estimated number of patients infected with M. tuberculosis L1.2.1 and L1.1.2 was over 1.1 million and of patients infected with L1.1.1 was about 200,000. This study demonstrated that L1 contributes significantly to the global burden of tuberculosis. To achieve the End TB Strategy, more attention needs to be paid to the responses of M. tuberculosis L1 to various control measures.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Cross-sectional Whole-genome Sequencing and Epidemiological Study of Multidrug-resistant Mycobacterium tuberculosis in China

Background

The increase in multidrug-resistant tuberculosis (MDR-TB) severely hampers tuberculosis prevention and control in China, a country with the second highest MDR-TB burden globally. The first nationwide drug-resistant tuberculosis surveillance program provides an opportunity to comprehensively investigate the epidemiological/drug-resistance characteristics, potential drug-resistance mutations, and effective population changes of Chinese MDR-TB.

Methods

We sequenced 357 MDR strains from 4600 representative tuberculosis-positive sputum samples collected during the survey (70 counties in 31 provinces). Drug-susceptibility testing was performed using 18 anti-tuberculosis drugs, representing the most comprehensive drug-resistance profile to date. We used 3 statistical and 1 machine-learning methods to identify drug-resistance genes/single-nucleotide polymorphisms (SNPs). We used Bayesian skyline analysis to investigate changes in effective population size.

Results

Epidemiological/drug-resistance characteristics showed different MDR profiles, co-resistance patterns, preferred drug combination/use, and recommended regimens among 7 Chinese administrative regions. These factors not only reflected the serious multidrug co-resistance and drug misuse but they were also potentially significant in facilitating the development of appropriate regimens for MDR-TB treatment in China. Further investigation identified 86 drug-resistance genes/intergenic regions/SNPs (58 new), providing potential targets for MDR-TB diagnosis and treatment. In addition, the effective population of Chinese MDR-TB displayed a strong expansion during 1993–2000, reflecting socioeconomic transition within the country. The phenomenon of expansion was restrained after 2000, likely attributable to the advances in diagnosis/treatment technologies and government support.

Conclusions

Our findings provide an important reference and improved understanding of MDR-TB in China, which are potentially significant in achieving the goal of precision medicine with respect to MDR-TB prevention and treatment.

Cryptic Resistance Mutations Associated With Misdiagnoses of Multidrug-Resistant Tuberculosis

Understanding why some multidrug-resistant tuberculosis cases are not detected by rapid phenotypic and genotypic routine clinical tests is essential to improve diagnostic assays and advance toward personalized tuberculosis treatment. Here, we combine whole-genome sequencing with single-colony phenotyping to identify a multidrug-resistant strain that had infected a patient for 9 years. Our investigation revealed the failure of rapid testing and genome-based prediction tools to identify the multidrug-resistant strain. The false-negative findings were caused by uncommon rifampicin and isoniazid resistance mutations. Although whole-genome sequencing data helped to personalize treatment, the patient developed extensively drug-resistant tuberculosis, highlighting the importance of coupling new diagnostic methods with appropriate treatment regimens.

The burden of pre-extensively and extensively drug-resistant tuberculosis among MDR-TB patients in the Amhara region, Ethiopia

BACKGROUND

The emergence of pre-extensively and extensively drug-resistant tuberculosis (Pre-XDR/XDR-TB) is the major hurdle for TB prevention and care programs especially in developing countries like Ethiopia. The less emphasis on universal access to laboratory techniques for the rapid diagnosis of TB and drug susceptibility testing (DST) makes the management of MDR-TB a challenge. Early detection of second line anti-TB drugs resistance is essential to reduce transmission of Pre-XDR/XDR-TB strains and adjusting the treatment regimen in MDR-TB.

OBJECTIVE

To determine the prevalence and resistance pattern of Pre-XDR- and XDR-TB among MDR-TB patients in the Amhara region, Ethiopia.

METHODS

A cross sectional study was carried out in nine MDR-TB treatment centers in the Amhara region. Sputum samples were collected from all pulmonary rifampicin resistant (RR) or MDR-TB patients prior to anti-TB treatment. Lӧwenstein-Jensen (LJ) culture, Ziehl Neelsen (ZN) smear, MTBDRplus and MTBDRsl assays were performed according to the standard procedures. Data were analyzed using SPSS 20 software. Chi-square and/or Fishers exact test was employed.

RESULTS

Overall, 6.3% of MDR-TB isolates were resistant to at least one second line drugs. Pre-XDR-TB and XDR-TB isolates accounted 5.7% and 0.6% respectively. Moreover, 3.4% were resistant to FQs and 3.4% were resistant to second line injectable drugs. All isolates were susceptible for low level kanamycin. Almost all pre-XDR-TB strains (90%) were previously treated with anti-TB drugs. Drug resistant Mycobacterium tuberculosis isolates were disproportionately distributed in districts of the Amhara region and the majorities were concentrated in urban areas.

CONCLUSIONS

The high proportion of MDR-TB patients resistant to at least one second line drug is alarming. Strengthening the laboratory facilities to monitor pre-XDR and XDR-TB patients is crucial. The TB programs need to give emphasis on the effective and rational use of second line drugs for newly diagnosed MDR-TB patients to prevent the emergence of pre-XDR/XDR-TB strains.

An ancestral genomic locus in Mycobacterium tuberculosis clinical isolates from India hints the genetic link with Mycobacterium canettii

BACKGROUND

Tuberculosis remains a worldwide public health emergency. To better understand M. tuberculosis and to identify genomic variations characteristic to the Indian clinical isolates by a low-cost method, a genomic subtractive hybridization between M. tuberculosis H37Rv and a clinical isolate from South India was performed.

RESULTS

This revealed a novel 0.4-kb subtractive fragment which was used as a handle to pull out a 4.5-kb genomic region characteristic to the clinical isolate and was absent in H37Rv. On further studies, this 4.5-kb region was found to be present in 91% of the M. tuberculosis clinical isolates screened from Kerala, a state in South India. Interestingly, this novel region has 99% identity (with 100% query coverage) with genomic regions of M. canettii.

DISCUSSION

The present study hypothesizes that this locus was present in the recent common environmental ancestor of mycobacteria, retained to the maximum extent in M. canettii and ancestral isolates of M. tuberculosis, and later deleted in other modern lineages of M. tuberculosis. Thus, this region may serve as one of the links between the pathogenic mycobacteria and the environmental species. We also propose that the Indian isolates of M. tuberculosis might be closely related to the putative progenitor M. prototuberculosis with respect to this locus. More studies on other genomic loci from different strains of M. tuberculosis are required to establish more links in this direction.

Tuberculosis in Liberia: high multidrug-resistance burden, transmission and diversity modelled by multiple importation events

Tuberculosis (TB) surveillance is scarce in most African countries, even though it is the continent with the greatest disease incidence according to the World Health Organization. Liberia is within the 30 countries with the highest TB burden, probably as a consequence of the long civil war and the recent Ebola outbreak, both crippling the health system and depreciating the TB prevention and control programmes. Due to difficulties working in the country, there is a lack of resistance surveys and bacillus characterization. Here, we use genome sequencing of Mycobacteriumtuberculosis clinical isolates to fill this gap. Our results highlight that the bacillus population structure is dominated by lineage 4 strains that harbour an outstanding genetic diversity, higher than in the rest of Africa as a whole. Coalescent analyses demonstrate that strains currently circulating in Liberia were introduced several times beginning in the early year 600 CE until very recently coinciding with migratory movements associated with the civil war and Ebola epidemics. A higher multidrug-resistant (MDR)-TB frequency (23.5 %) than current estimates was obtained together with non-catalogued drug-resistance mutations. Additionally, 39 % of strains were in genomic clusters revealing that ongoing transmission is a major contribution to the TB burden in the country. Our report emphasizes the importance of TB surveillance and control in African countries where bacillus diversity, MDR-TB prevalence and transmission are coalescing to jeopardize TB control programmes.

Characterization of Polymorphisms Associated with Multidrug-Resistant Tuberculosis by Whole Genomic Sequencing: A Preliminary Report from Mexico

Whole genome sequencing (WGS) has been proposed as a tool for the diagnosis of drug resistance in tuberculosis (TB); however, there have been few studies on its effectiveness in countries with significantly high drug resistance rates. This study therefore aimed to evaluate the effectiveness of WGS to identify mutations related to drug resistance in TB isolates from an endemic region of Mexico. The results showed that, of 35 multidrug-resistant isolates analyzed, the values of congruence found between the phenotypic drug susceptibility testing and polymorphisms were 94% for isoniazid, 97% for rifampicin, 90% for ethambutol, and 82% for pyrazinamide. It was also possible to identify eight isolates as potential pre-extensive drug resistant (XDR) and one as XDR. Twenty nine isolates were classified within L4 and two transmission clusters were identified. The results show the potential utility of WGS for predicting resistance against first- and second-line drugs, as well as providing a phylogenetic characterization of TB drug-resistant isolates circulating in Mexico.

Isolation and comparative genomics of Mycobacterium tuberculosis isolates from cattle and their attendants in South India

The major human pathogen Mycobacterium tuberculosis is rarely reported to cause disease in other animals. Cases in livestock are thought to occur through contact with infected handlers, but previous studies evaluating putative livestock-human transmission used typing techniques with limited resolution. Here, we undertook cross-sectional surveillance for tuberculosis in 271 livestock handlers and 167 cattle on three farms in Chennai, India and defined the relatedness of cultured isolates using whole genome sequencing. Humans and livestock were screened for active mycobacterial infection, and opportunistic post-mortem examination was performed on comparative intradermal test-positive cattle that died. Four cattle and 6 handlers on two farms were culture-positive for M. tuberculosis; M. bovis was not isolated. All 10 isolates (one from each case) belonged to Lineage 1. Pairwise genome comparisons of single nucleotide polymorphism (SNP) differences ranged from 1 to 600 SNPs, but 3 isolate pairs were less than 5 SNPs different. Two pairs were from handlers and the third pair were from two cattle on the same farm. The minimum pairwise SNP difference between a cattle and human isolate was >250 SNPs. Our study confirms the presence of M. tuberculosis infection in cattle in India, sequencing of which characterised relatedness between human and cattle-derived isolates.

A review of published spoligotype data indicates the diversity of Mycobacterium tuberculosis from India is under-represented in global databases

BACKGROUND

Mycobacterium tuberculosis (MTBC) lineages differ in clinical presentation, virulence, transmission, drug resistance and immunological responses. Despite having the largest burden of tuberculosis (TB) in the world, strains from India are underrepresented in international databases. We reviewed published spoligotype data to determine the distribution and diversity of MTBC lineages in India.

METHODS

A Pubmed/MEDLINE search identified 34 M. tuberculosis spoligotyping studies from India. Spoligotype patterns were extracted and the Spoligotype International Type (SIT) number, sub-lineage and lineage determined. Minimum Spanning Trees were used to determine relationships between patterns.

RESULTS

We identified 1528 spoligotype patterns distributed across 8300 isolates; 6733 isolates belonged to 472 SITs, with 53% of all isolates belonging to 12 SITs with at least 100 isolates each. Lineage 1 and Lineage 3 made up 67% of all isolates, although a lineage could not be assigned for 16% of isolates. Lineage 1 isolates were most common in Southern, Western and Eastern India, and Lineage 3 was most common in Northern and Central India. The RULE, CBN and KBBN lineage prediction algorithms from the TB-lineage tools performed variably, with the correct lineage predicted correctly for only 64% of patterns with known lineage. Using a consensus definition, 64% of the 1359 isolates with unknown lineage were assigned to Lineage 1, and 14% each were assigned to Lineages 3 and 4. With these lineage assignments, 80% of all isolates belonged to either Lineage 1 or Lineage 3.

CONCLUSION

Our findings indicate significant M. tuberculosis diversity in India. The documentation of 1056 orphan and unreported patterns indicate that this diversity is under-represented in global databases.

Deciphering drug resistance in Mycobacterium tuberculosis using whole-genome sequencing: progress, promise, and challenges

Tuberculosis (TB) is a global infectious threat that is intensified by an increasing incidence of highly drug-resistant disease. Whole-genome sequencing (WGS) studies of Mycobacterium tuberculosis, the causative agent of TB, have greatly increased our understanding of this pathogen. Since the first M. tuberculosis genome was published in 1998, WGS has provided a more complete account of the genomic features that cause resistance in populations of M. tuberculosis, has helped to fill gaps in our knowledge of how both classical and new antitubercular drugs work, and has identified specific mutations that allow M. tuberculosis to escape the effects of these drugs. WGS studies have also revealed how resistance evolves both within an individual patient and within patient populations, including the important roles of de novo acquisition of resistance and clonal spread. These findings have informed decisions about which drug-resistance mutations should be included on extended diagnostic panels. From its origins as a basic science technique, WGS of M. tuberculosis is becoming part of the modern clinical microbiology laboratory, promising rapid and improved detection of drug resistance, and detailed and real-time epidemiology of TB outbreaks. We review the successes and highlight the challenges that remain in applying WGS to improve the control of drug-resistant TB through monitoring its evolution and spread, and to inform more rapid and effective diagnostic and therapeutic strategies.

Identification and Characterization of Genetic Determinants of Isoniazid and Rifampicin Resistance in Mycobacterium tuberculosis in Southern India

Drug-resistant tuberculosis (TB), one of the leading causes of death worldwide, arises mainly from spontaneous mutations in the genome of Mycobacterium tuberculosis. There is an urgent need to understand the mechanisms by which the mutations confer resistance in order to identify new drug targets and to design new drugs. Previous studies have reported numerous mutations that confer resistance to anti-TB drugs, but there has been little systematic analysis to understand their genetic background and the potential impacts on the drug target stability and/or interactions. Here, we report the analysis of whole-genome sequence data for 98 clinical M. tuberculosis isolates from a city in southern India. The collection was screened for phenotypic resistance and sequenced to mine the genetic mutations conferring resistance to isoniazid and rifampicin. The most frequent mutation among isoniazid and rifampicin isolates was S315T in katG and S450L in rpoB respectively. The impacts of mutations on protein stability, protein-protein interactions and protein-ligand interactions were analysed using both statistical and machine-learning approaches. Drug-resistant mutations were predicted not only to target active sites in an orthosteric manner, but also to act through allosteric mechanisms arising from distant sites, sometimes at the protein-protein interface.

Insight into genetic diversity of Mycobacterium tuberculosis in Kandy, Sri Lanka reveals predominance of the Euro-American lineage

OBJECTIVE

Sri Lanka is a country where the molecular epidemiology of Mycobacterium tuberculosis (MTB) is poorly explored. Therefore, this study was performed to identify circulating lineages/sub-lineages of MTB and their transmission patterns.

METHODS

DNA was extracted from 89 isolates of MTB collected during 2012 and 2013 from new pulmonary tuberculosis patients in Kandy, Sri Lanka and analyzed by spoligotyping, large sequence polymorphism (LSP), mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing, and drug resistance-associated gene sequencing.

RESULTS

The predominant lineage was lineage 4 (Euro-American, 45.9%), followed by lineage 1 (Indo-Oceanic, 29.4%), lineage 2 (East-Asian, 23.5%), and lineage 3 (Central-Asian, 1.2%). Among 26 spoligotype patterns, eight were undesignated or new types and seven of these belonged to lineage 4. Undesignated lineage 4/SIT124 (n=2/8) and SIT3234 (n=8/8) clustered together based on 24-locus MIRU-VNTR typing. The dominant sub-lineage was Beijing/SIT1 (n=19), with the isoniazid resistance katG G944C mutation (Ser315Thr) detected in two of them.

CONCLUSIONS

The population structure of MTB in Kandy, Sri Lanka was different from that in the South Asian region. The clonal expansion of locally evolved lineage 4/SIT3234 and detection of the pre-multidrug resistant Beijing isolates from new tuberculosis patients is alarming and will require continuous monitoring.

Whole genomic sequencing as a tool for diagnosis of drug and multidrug-resistance tuberculosis in an endemic region in Mexico

Background

Whole genome sequencing (WGS) has been proposed as a tool for diagnosing drug resistance in tuberculosis. However, reports of its effectiveness in endemic countries with important numbers of drug resistance are scarce. The goal of this study was to evaluate the effectiveness of this procedure in isolates from a tuberculosis endemic region in Mexico.

Methods

WGS analysis was performed in 81 tuberculosis positive clinical isolates with a known phenotypic profile of resistance against first-line drugs (isoniazid, rifampin, ethambutol, pyrazinamide and streptomycin). Mutations related to drug resistance were identified for each isolate; drug resistant genotypes were predicted and compared with the phenotypic profile. Genotypes and transmission clusters based on genetic distances were also characterized.

Findings

Prediction by WGS analysis of resistance against isoniazid, rifampicin, ethambutol, pyrazinamide and streptomycin showed sensitivity values of 84%, 96%, 71%, 75% and 29%, while specificity values were 100%, 94%, 90%, 90% and 98%, respectively. Prediction of multidrug resistance showed a sensitivity of 89% and specificity of 97%. Moreover, WGS analysis revealed polymorphisms related to second-line drug resistance, enabling classification of eight and two clinical isolates as pre- and extreme drug-resistant cases, respectively. Lastly, four lineages were identified in the population (L1, L2, L3 and L4). The most frequent of these was L4, which included 90% (77) of the isolates. Six transmission clusters were identified; the most frequent was TC6, which included 13 isolates with a L4.1.1 and a predominantly multidrug-resistant condition.

Conclusions

The results illustrate the utility of WGS for establishing the potential for prediction of resistance against first and second line drugs in isolates of tuberculosis from the region. They also demonstrate the feasibility of this procedure for use as a tool to support the epidemiological surveillance of drug- and multidrug-resistant tuberculosis.

Deciphering Within-Host Microevolution of Mycobacterium tuberculosis through Whole-Genome Sequencing: the Phenotypic Impact and Way Forward

The Mycobacterium tuberculosis genome is more heterogenous and less genetically stable within the host than previously thought. Currently, only limited data exist on the within-host microevolution, diversity, and genetic stability of M. tuberculosis As a direct consequence, our ability to infer M. tuberculosis transmission chains and to understand the full complexity of drug resistance profiles in individual patients is limited. Furthermore, apart from the acquisition of certain drug resistance-conferring mutations, our knowledge on the function of genetic variants that emerge within a host and their phenotypic impact remains scarce. We performed a systematic literature review of whole-genome sequencing studies of serial and parallel isolates to summarize the knowledge on genetic diversity and within-host microevolution of M. tuberculosis We identified genomic loci of within-host emerged variants found across multiple studies and determined their functional relevance. We discuss important remaining knowledge gaps and finally make suggestions on the way forward.

Genome Sequencing of Polydrug-, Multidrug-, and Extensively Drug-Resistant Mycobacterium tuberculosis Strains from South India

The genomes of 16 clinical Mycobacterium tuberculosis isolates were subjected to whole-genome sequencing to identify mutations related to resistance to one or more anti-Mycobacterium drugs. The sequence data will help in understanding the genomic characteristics of M. tuberculosis isolates and their resistance mutations prevalent in South India.

Whole Genome Sequencing of Mycobacterium tuberculosis Clinical Isolates From India Reveals Genetic Heterogeneity and Region-Specific Variations That Might Affect Drug Susceptibility

Whole genome sequencing (WGS) of Mycobacterium tuberculosis has been constructive in understanding its evolution, genetic diversity and the mechanisms involved in drug resistance. A large number of sequencing efforts from across the globe have revealed genetic diversity among clinical isolates and the genetic determinants for their resistance to anti-tubercular drugs. Considering the high TB burden in India, the availability of WGS studies is limited. Here we present, WGS results of 200 clinical isolates of M. tuberculosis from North India which are categorized as sensitive to first-line drugs, mono-resistant, multi-drug resistant and pre-extensively drug resistant isolates. WGS revealed that 20% of the isolates were co-infected with M. tuberculosis and non-tuberculous mycobacteria species. We identified 12,802 novel genetic variations in M. tuberculosis isolates including 343 novel SNVs in 38 genes which are known to be associated with drug resistance and are not currently used in the diagnostic kits for detection of drug resistant TB. We also identified M. tuberculosis lineage 3 to be predominant in the northern region of India. Additionally, several novel SNVs, which may potentially confer drug resistance were found to be enriched in the drug resistant isolates sampled. This study highlights the significance of employing WGS in diagnosis and for monitoring further development of MDR-TB strains.

The Future of TB Resistance Diagnosis: The Essentials on Whole Genome Sequencing and Rapid Testing Methods

Tuberculosis resistance diagnostics have vastly improved in recent years thanks to the development of standardised phenotypic and molecular testing methods. However, these methods are either slow or limited in the number of resistant genotypes they can detect. With the advent of next-generation sequencing (NGS) we can sidestep all those problems, as we can sequence whole tuberculosis genomes at increasingly smaller costs and requiring less and less DNA. In this review, we explain how accumulated knowledge in the field has allowed us to go from phenotypic testing to molecular methods to Whole Genome Sequencing (WGS) for resistance diagnostics. We compare current diagnostic methods with WGS as to their efficacy in detecting resistant cases, and show how forthcoming advances in NGS technologies will be crucial in widespread implementation of WGS as a diagnostic tool.

Evidence for Highly Variable, Region-Specific Patterns of T-Cell Epitope Mutations Accumulating in Mycobacterium tuberculosis Strains

Vaccines that confer protection through induction of adaptive T-cell immunity rely on understanding T-cell epitope (TCE) evolution induced by immune escape. This is poorly understood in tuberculosis (TB), an ancient, chronic disease, where CD4 T-cell immunity is of recognized importance. We probed 905 functionally validated, curated human CD4 T cell epitopes in 79 Mycobacterium tuberculosis (Mtb) whole genomes from India. This screen resulted in identifying 64 mutated epitopes in these strains initially using a computational pipeline and subsequently verified by single nucleotide polymorphism (SNP) analysis. SNP based phylogeny revealed the 79 Mtb strains to cluster to East African Indian (EAI), Central Asian Strain (CAS), and Beijing (BEI) lineages. Eighty-nine percent of the mutated T-cell epitopes (mTCEs) identified in the 79 Mtb strains from India has not previously been reported. These mTCEs were encoded by genes with high nucleotide diversity scores including seven mTCEs encoded by six antigens in the top 10% of rapidly divergent Mtb genes encoded by these strains. Using a T cell functional assay readout, we demonstrate 62% of mTCEs tested to significantly alter CD4 T-cell IFNγ and/or IL2 secretion with associated changes in predicted HLA-DR binding affinity: the gain of function mutations displayed higher predicted HLA-DR binding affinity and conversely mutations resulting in loss of function displayed lower predicted HLA-DR binding affinity. Most mutated antigens belonged to the cell wall/cell processes, and, intermediary metabolism and respiration families though all known Mtb proteins encoded mutations. Analysis of the mTCEs in an SNP database of 5,310 global Mtb strains identified 82% mTCEs to be significantly more prevalent in Mtb strains isolated from India, including 36 mTCEs identified exclusively in strains from India. These epitopes had a significantly higher predicted binding affinity to HLA-DR alleles that were highly prevalent in India compared to HLA-DR alleles rare in India, highlighting HLA-DR maybe an important driver of these mutations. This first evidence of region-specific TCE mutations potentially employed by Mtb to escape host immunity has important implications for TB vaccine design.

Validation of Novel Mycobacterium tuberculosis Isoniazid Resistance Mutations Not Detectable by Common Molecular Tests

Resistance to the first-line antituberculosis (TB) drug isoniazid (INH) is widespread, and the mechanism of resistance is unknown in approximately 15% of INH-resistant (INH-R) strains. To improve molecular detection of INH-R TB, we used whole-genome sequencing (WGS) to analyze 52 phenotypically INH-R Mycobacterium tuberculosis complex (MTBC) clinical isolates that lacked the common katG S315T or inhA promoter mutations. Approximately 94% (49/52) of strains had mutations at known INH-associated loci that were likely to confer INH resistance. All such mutations would be detectable by sequencing more DNA adjacent to existing target regions. Use of WGS minimized the chances of missing infrequent INH resistance mutations outside commonly targeted hotspots. We used recombineering to generate 12 observed clinical katG mutations in the pansusceptible H37Rv reference strain and determined their impact on INH resistance. Our functional genetic experiments have confirmed the role of seven suspected INH resistance mutations and discovered five novel INH resistance mutations. All recombineered katG mutations conferred resistance to INH at a MIC of ≥0.25 μg/ml and should be added to the list of INH resistance determinants targeted by molecular diagnostic assays. We conclude that WGS is a useful tool for detecting uncommon INH resistance mutations that would otherwise be missed by current targeted molecular testing methods and suggest that its use (or use of expanded conventional or next-generation-based targeted sequencing) may provide earlier diagnosis of INH-R TB.

Use of whole genome sequencing in surveillance of drug resistant tuberculosis

INTRODUCTION

The threat of resistance to anti-tuberculosis drugs is of global concern. Current efforts to monitor resistance rely on phenotypic testing where cultured bacteria are exposed to critical concentrations of the drugs. Capacity for such testing is low in TB endemic countries. Drug resistance is caused by mutations in the Mycobacterium tuberculosis genome and whole genome sequencing to detect these mutations offers an alternative means of assessing resistance. Areas covered: The challenges of assessing TB drug resistance are discussed. Progress in elucidating the M. tuberculosis resistome and evidence of the accuracy of next generation sequencing for detecting resistance is reviewed. Expert Commentary: There are considerable advantages to using next generation sequencing for TB drug resistance surveillance. Accuracy is high for detecting resistance to the major first-line drugs but is currently lower for the second-line drugs due to our incomplete knowledge regarding resistance causing mutations. With the advances in sequencing technology and the opportunity to replace phenotypic drug susceptibility testing with safer and more cost effective methods it would appear that the question is when to implement. Current bottlenecks are sample extraction to allow whole genome sequencing directly from sputum and the lack of bioinformatics expertise in some TB endemic countries.

Evolution of drug resistance in Mycobacterium tuberculosis: a review on the molecular determinants of resistance and implications for personalized care

Drug-resistant TB (DR-TB) remains a significant challenge in TB treatment and control programmes worldwide. Advances in sequencing technology have significantly increased our understanding of the mechanisms of resistance to anti-TB drugs. This review provides an update on advances in our understanding of drug resistance mechanisms to new, existing drugs and repurposed agents. Recent advances in WGS technology hold promise as a tool for rapid diagnosis and clinical management of TB. Although the standard approach to WGS of Mycobacterium tuberculosis is slow due to the requirement for organism culture, recent attempts to sequence directly from clinical specimens have improved the potential to diagnose and detect resistance within days. The introduction of new databases may be helpful, such as the Relational Sequencing TB Data Platform, which contains a collection of whole-genome sequences highlighting key drug resistance mutations and clinical outcomes. Taken together, these advances will help devise better molecular diagnostics for more effective DR-TB management enabling personalized treatment, and will facilitate the development of new drugs aimed at improving outcomes of patients with this disease.

Rapid molecular diagnostics for multi-drug resistant tuberculosis in India

INTRODUCTION

Rapid molecular diagnostic methods help in the detection of TB and Rifampicin resistance. These methods detect TB early, are accurate and play a crucial role in reducing the burden of drug resistant tuberculosis. Areas covered: This review analyses rapid molecular diagnostic tools used in the diagnosis of MDR-TB in India, such as the Line Probe Assay and GeneXpert. We have discussed the burden of MDR-TB and the impact of recent diagnostic tools on case detection and treatment outcomes. This review also discusses the costs involved in establishing these new techniques in India. Expert commentary: Molecular methods have considerable advantages for the programmatic management of drug resistant TB. These include speed, standardization of testing, potentially high throughput and reduced laboratory biosafety requirements. There is a desperate need for India to adopt modern, rapid, molecular tools with point-of-care tests being currently evaluated. New molecular diagnostic tests appear to be cost effective and also help in detecting missing cases. There is enough evidence to support the scaling up of these new tools in India.

Whole-Genome Sequencing of Bacterial Pathogens: the Future of Nosocomial Outbreak Analysis

Outbreaks of multidrug-resistant bacteria present a frequent threat to vulnerable patient populations in hospitals around the world. Intensive care unit (ICU) patients are particularly susceptible to nosocomial infections due to indwelling devices such as intravascular catheters, drains, and intratracheal tubes for mechanical ventilation. The increased vulnerability of infected ICU patients demonstrates the importance of effective outbreak management protocols to be in place. Understanding the transmission of pathogens via genotyping methods is an important tool for outbreak management. Recently, whole-genome sequencing (WGS) of pathogens has become more accessible and affordable as a tool for genotyping. Analysis of the entire pathogen genome via WGS could provide unprecedented resolution in discriminating even highly related lineages of bacteria and revolutionize outbreak analysis in hospitals. Nevertheless, clinicians have long been hesitant to implement WGS in outbreak analyses due to the expensive and cumbersome nature of early sequencing platforms. Recent improvements in sequencing technologies and analysis tools have rapidly increased the output and analysis speed as well as reduced the overall costs of WGS. In this review, we assess the feasibility of WGS technologies and bioinformatics analysis tools for nosocomial outbreak analyses and provide a comparison to conventional outbreak analysis workflows. Moreover, we review advantages and limitations of sequencing technologies and analysis tools and present a real-world example of the implementation of WGS for antimicrobial resistance analysis. We aimed to provide health care professionals with a guide to WGS outbreak analysis that highlights its benefits for hospitals and assists in the transition from conventional to WGS-based outbreak analysis.