Whole genome sequencing of Mycobacterium tuberculosis for prediction of drug resistance

An overview of next generation sequencing strategies and genomics tools used for tuberculosis research

Tuberculosis (TB) is a grave public health concern and is considered the foremost contributor to human mortality resulting from infectious disease. Due to the stringent clonality and extremely restricted genomic diversity, conventional methods prove inefficient for in-depth exploration of minor genomic variations and the evolutionary dynamics operating in Mycobacterium tuberculosis (M.tb) populations. Until now, the majority of reviews have primarily focused on delineating the application of whole-genome sequencing (WGS) in predicting antibiotic resistant genes, surveillance of drug resistance strains, and M.tb lineage classifications. Despite the growing use of next generation sequencing (NGS) and WGS analysis in TB research, there are limited studies that provide a comprehensive summary of there role in studying macroevolution, minor genetic variations, assessing mixed TB infections, and tracking transmission networks at an individual level. This highlights the need for systematic effort to fully explore the potential of WGS and its associated tools in advancing our understanding of TB epidemiology and disease transmission. We delve into the recent bioinformatics pipelines and NGS strategies that leverage various genetic features and simultaneous exploration of host-pathogen protein expression profile to decipher the genetic heterogeneity and host-pathogen interaction dynamics of the M.tb infections. This review highlights the potential benefits and limitations of NGS and bioinformatics tools and discusses their role in TB detection and epidemiology. Overall, this review could be a valuable resource for researchers and clinicians interested in NGS-based approaches in TB research.

New insight in molecular detection of Mycobacterium tuberculosis

Mycobacterium tuberculosis, the causative agent of tuberculosis, is a pathogenic bacterium that has claimed millions of lives since the Middle Ages. According to the World Health Organization's report, tuberculosis ranks among the ten deadliest diseases worldwide. The presence of an extensive array of genes and diverse proteins within the cellular structure of this bacterium has provided us with a potent tool for diagnosis. While the culture method remains the gold standard for tuberculosis diagnosis, it is possible that molecular diagnostic methods, emphasis on the identification of mutation genes (e.g., rpoB and gyrA) and single nucleotide polymorphisms, could offer a safe and reliable alternative. Over the past few decades, as our understanding of molecular genetics has expanded, methods have been developed based on gene expansion and detection. These methods typically commence with DNA amplification through nucleic acid targeted techniques such as polymerase chain reaction. Various molecular compounds and diverse approaches have been employed in molecular assays. In this review, we endeavor to provide an overview of molecular assays for the diagnosis of tuberculosis with their properties (utilization, challenges, and functions). The ultimate goal is to explore the potential of replacing traditional bacterial methods with these advanced molecular diagnostic techniques.

Exploring diagnostic methods for drug-resistant tuberculosis: A comprehensive overview

Despite available global efforts and funding, Tuberculosis (TB) continues to affect a considerable number of patients worldwide. Policy makers and stakeholders set clear goals to reduce TB incidence and mortality, but the emergence of multidrug-resistant TB (MDR-TB) and extensively drug-resistant TB (XDR-TB) complicate the reach of these goals. Drug-resistance TB needs to be diagnosed rapidly and accurately to effectively treat patients, prevent the transmission of MDR-TB, minimise mortality, reduce treatment costs and avoid unnecessary hospitalisations. In this narrative review, we provide a comprehensive overview of laboratory methods for detecting drug resistance in MTB, focusing on phenotypic, molecular and other drug susceptibility testing (DST) techniques. We found a large variety of methods used, with the BACTEC MGIT 960 being the most common phenotypic DST and the Xpert MTB/RIF being the most common molecular DST. We emphasise the importance of integrating phenotypic and molecular DST to address issues like resistance to new drugs, heteroresistance, mixed infections and low-level resistance mutations. Notably, most of the analysed studies adhered to the outdated definition of XDR-TB and did not consider the pre-XDR definition, thus posing challenges in aligning diagnostic methods with the current landscape of TB resistance.

Identification of novel single nucleotide variants in the drug resistance mechanism of Mycobacterium tuberculosis isolates by whole-genome analysis

BACKGROUND

Tuberculosis (TB) represents a major global health challenge. Drug resistance in Mycobacterium tuberculosis (MTB) poses a substantial obstacle to effective TB treatment. Identifying genomic mutations in MTB isolates holds promise for unraveling the underlying mechanisms of drug resistance in this bacterium.

METHODS

In this study, we investigated the roles of single nucleotide variants (SNVs) in MTB isolates resistant to four antibiotics (moxifloxacin, ofloxacin, amikacin, and capreomycin) through whole-genome analysis. We identified the drug-resistance-associated SNVs by comparing the genomes of MTB isolates with reference genomes using the MuMmer4 tool.

RESULTS

We observed a strikingly high proportion (94.2%) of MTB isolates resistant to ofloxacin, underscoring the current prevalence of drug resistance in MTB. An average of 3529 SNVs were detected in a single ofloxacin-resistant isolate, indicating a mutation rate of approximately 0.08% under the selective pressure of ofloxacin exposure. We identified a set of 60 SNVs associated with extensively drug-resistant tuberculosis (XDR-TB), among which 42 SNVs were non-synonymous mutations located in the coding regions of nine key genes (ctpI, desA3, mce1R, moeB1, ndhA, PE_PGRS4, PPE18, rpsA, secF). Protein structure modeling revealed that SNVs of three genes (PE_PGRS4, desA3, secF) are close to the critical catalytic active sites in the three-dimensional structure of the coding proteins.

CONCLUSION

This comprehensive study elucidates novel resistance mechanisms in MTB against antibiotics, paving the way for future design and development of anti-tuberculosis drugs.

Rapid Diagnosis of Drug-Resistant Tuberculosis-Opportunities and Challenges

Global tuberculosis (TB) eradication is undermined by increasing prevalence of emerging resistance to available drugs, fuelling ongoing demand for more complex diagnostic and treatment strategies. Early detection of TB drug resistance coupled with therapeutic decision making guided by rapid characterisation of pre-treatment and treatment emergent resistance remains the most effective strategy for averting Drug-Resistant TB (DR-TB) transmission, reducing DR-TB associated mortality, and improving patient outcomes. Solid- and liquid-based mycobacterial culture methods remain the gold standard for Mycobacterium tuberculosis (MTB) detection and drug susceptibility testing (DST). Unfortunately, delays to result return, and associated technical challenges from requirements for specialised resource and capacity, have limited DST use and availability in many high TB burden resource-limited countries. There is increasing availability of a variety of rapid nucleic acid-based diagnostic assays with adequate sensitivity and specificity to detect gene mutations associated with resistance to one or more drugs. While a few of these assays produce comprehensive calls for resistance to several first- and second-line drugs, there is still no endorsed genotypic drug susceptibility test assay for bedaquiline, pretomanid, and delamanid. The global implementation of regimens comprising these novel drugs in the absence of rapid phenotypic drug resistance profiling has generated a new set of diagnostic challenges and heralded a return to culture-based phenotypic DST. In this review, we describe the available tools for rapid diagnosis of drug-resistant tuberculosis and discuss the associated opportunities and challenges.

A precision overview of genomic resistance screening in Ecuadorian isolates of Mycobacterium tuberculosis using web-based bioinformatics tools

INTRODUCTION

Tuberculosis (TB) is among the deadliest diseases worldwide, and its impact is mainly due to the continuous emergence of resistant isolates during treatment due to the laborious process of resistance diagnosis, nonadherence to treatment and circulation of previously resistant isolates of Mycobacterium tuberculosis. In this study, we evaluated the performance and functionalities of web-based tools, including Mykrobe, TB-profiler, PhyResSE, KvarQ, and SAM-TB, for detecting resistance in 88 Ecuadorian isolates of Mycobacterium tuberculosis drug susceptibility tested previously. Statistical analysis was used to determine the correlation between genomic and phenotypic analysis. Our results showed that with the exception of KvarQ, all tools had the highest correlation with the conventional drug susceptibility test (DST) for global resistance detection (98% agreement and 0.941 Cohen's kappa), while SAM-TB, PhyResSE, TB-profiler and Mykrobe had better correlations with DST for first-line drug analysis individually. We also identified that in our study, only 50% of mutations characterized by the web-based tools in the rpoB, katG, embB, pncA, gyrA and rrs regions were canonical and included in the World Health Organization (WHO) catalogue. Our findings suggest that SAM-TB, PhyResSE, TB-profiler and Mykrobe were efficient in determining canonical resistance-related mutations, but more analysis is needed to improve second-line detection. Improving surveillance programs using whole-genome sequencing tools for first-line drugs, MDR-TB and XDR-TB is essential to understand the molecular epidemiology of TB in Ecuador.

IMPORTANCE

Tuberculosis, an infectious disease caused by Mycobacterium tuberculosis, most commonly affects the lungs and is often spread through the air when infected people cough, sneeze, or spit. However, despite the existence of effective drug treatment, patient adherence, long duration of treatment, and late diagnosis have reduced the effectiveness of therapy and increased drug resistance. The increase in resistant cases, added to the impact of the COVID-19 pandemic, has highlighted the importance of implementing efficient and timely diagnostic methodologies worldwide. The significance of our research is in evaluating and identifying a more efficient and user-friendly web-based tool to characterize resistance in Mycobacterium tuberculosis by whole-genome sequencing, which will allow more routine application to improve TB strain surveillance programs locally.

Pediatric Drug-Resistant Tuberculosis: The Current and Future Prospects for Management and Prevention

In the continued battle against one of the oldest enemies known to mankind, Mycobacterium tuberculosis (MTB), the emergence of drug resistance to antituberculosis drugs among children poses multiple challenges for early detection and treatment. Molecular diagnostics and newer drugs like bedaquiline and delamanid have strengthened the armamentarium and helped design convenient, safe, and child-friendly therapeutic regimens against drug-resistant tuberculosis (TB). Preventive strategies like treatment of TB infection among children living in close contact with patients with drug-resistant TB and effective vaccines against TB are currently in the investigative stages of development and implementation. In addition to the implementation of recent novel diagnostics and treatment modalities, effective psychosocial and nutritional support, as well as dedicated monitoring for compliance and adverse effects, are crucial determinants for successful treatment outcomes in these children.

Resistance patterns and transmission of mono- and polyresistant TB: clinical impact of WGS

Objectives

Rapidly diagnosing drug-resistant TB is crucial for improving treatment and transmission control. WGS is becoming increasingly accessible and has added value to the diagnosis and treatment of TB. The aim of the study was to perform WGS to determine the rate of false-positive results of phenotypic drug susceptibility testing (pDST) and characterize the molecular mechanisms of resistance and transmission of mono- and polyresistant Mycobacterium (M.) tuberculosis.

Methods

WGS was performed on 53 monoresistant and 25 polyresistant M. tuberculosis isolates characterized by pDST. Sequencing data were bioinformatically processed to infer mutations encoding resistance and determine the origin of resistance and phylogenetic relationship between isolates studied.

Results

The data showed the variable sensitivity and specificity of WGS in comparison with pDST as the gold standard: isoniazid 92.7% and 92.3%; streptomycin 41.9% and 100.0%; pyrazinamide 15% and 94.8%; and ethambutol 75.0% and 98.6%, respectively. We found novel mutations encoding resistance to streptomycin (in gidB) and pyrazinamide (in kefB). Most isolates belonged to lineage 4 (80.1%) and the overall clustering rate was 11.5%. We observed lineage-specific gene variations encoding resistance to streptomycin and pyrazinamide.

Conclusions

This study highlights the clinical potential of WGS in ruling out false-positive drug resistance following phenotypic or genetic drug testing, and recommend this technology together with the WHO catalogue in designing an optimal individualized treatment regimen and preventing the development of MDR TB. Our results suggest that resistance is primarily developed through spontaneous mutations or selective pressure.

Nanomotion technology in combination with machine learning: a new approach for a rapid antibiotic susceptibility test for Mycobacterium tuberculosis

Nanomotion technology is a growth-independent approach that can be used to detect and record the vibrations of bacteria attached to cantilevers. We have developed a nanomotion-based antibiotic susceptibility test (AST) protocol for Mycobacterium tuberculosis (MTB). The protocol was used to predict strain phenotype towards isoniazid (INH) and rifampicin (RIF) using a leave-one-out cross-validation (LOOCV) and machine learning techniques. This MTB-nanomotion protocol takes 21 h, including cell suspension preparation, optimized bacterial attachment to functionalized cantilever, and nanomotion recording before and after antibiotic exposure. We applied this protocol to MTB isolates (n = 40) and were able to discriminate between susceptible and resistant strains for INH and RIF with a maximum sensitivity of 97.4% and 100%, respectively, and a maximum specificity of 100% for both antibiotics when considering each nanomotion recording to be a distinct experiment. Grouping recordings as triplicates based on source isolate improved sensitivity and specificity to 100% for both antibiotics. Nanomotion technology can potentially reduce time-to-result significantly compared to the days and weeks currently needed for current phenotypic ASTs for MTB. It can further be extended to other anti-TB drugs to help guide more effective TB 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.

Mutations and insights into the molecular mechanisms of resistance of Mycobacterium tuberculosis to first-line

Genetically antimicrobial resistance in Mycobacterium tuberculosis is currently one of the most important aspects of tuberculosis, considering that there are emerging resistant strains for almost every known drug used for its treatment. There are multiple antimicrobials used for tuberculosis treatment, and the most effective ones are the first-line drugs, which include isoniazid, pyrazinamide, rifampicin, and ethambutol. In this context, understanding the mechanisms of action and resistance of these molecules is essential for proposing new therapies and strategies of treatment. Additionally, understanding how and where mutations arise conferring a resistance profile to the bacteria and their effect on bacterial metabolism is an important requisite to be taken in producing safer and less susceptible drugs to the emergence of resistance. In this review, we summarize the most recent literature regarding novel mutations reported between 2017 and 2022 and the advances in the molecular mechanisms of action and resistance against first-line drugs used in tuberculosis treatment, highlighting recent findings in pyrazinamide resistance involving PanD and, additionally, resistance-conferring mutations for novel drugs such as bedaquiline, pretomanid, delamanid and linezolid.

Comparative genomics of drug-resistant strains of Mycobacterium tuberculosis in Ecuador

BACKGROUND

Tuberculosis is a serious infectious disease affecting millions of people. In spite of efforts to reduce the disease, increasing antibiotic resistance has contributed to persist in the top 10 causes of death worldwide. In fact, the increased cases of multi (MDR) and extreme drug resistance (XDR) worldwide remains the main challenge for tuberculosis control. Whole genome sequencing is a powerful tool for predicting drug resistance-related variants, studying lineages, tracking transmission, and defining outbreaks. This study presents the identification and characterization of resistant clinical isolates of Mycobacterium tuberculosis including a phylogenetic and molecular resistance profile study by sequencing the complete genome of 24 strains from different provinces of Ecuador.

RESULTS

Genomic sequencing was used to identify the variants causing resistance. A total of 15/21 isolates were identified as MDR, 4/21 as pre-XDR and 2/21 as XDR, with three isolates discarded due to low quality; the main sub-lineage was LAM (61.9%) and Haarlem (19%) but clades X, T and S were identified. Of the six pre-XDR and XDR strains, it is noteworthy that five come from females; four come from the LAM sub-lineage and two correspond to the X-class sub-lineage. A core genome of 3,750 genes, distributed in 295 subsystems, was determined. Among these, 64 proteins related to virulence and implicated in the pathogenicity of M. tuberculosis and 66 possible pharmacological targets stand out. Most variants result in nonsynonymous amino acid changes and the most frequent genotypes were identified as conferring resistance to rifampicin, isoniazid, ethambutol, para-aminosalicylic acid and streptomycin. However, an increase in the resistance to fluoroquinolones was detected.

CONCLUSION

This work shows for the first time the variability of circulating resistant strains between men and women in Ecuador, highlighting the usefulness of genomic sequencing for the identification of emerging resistance. In this regard, we found an increase in fluoroquinolone resistance. Further sampling effort is needed to determine the total variability and associations with the metadata obtained to generate better health policies.

Association between Mycobacterium tuberculosis genotype and diabetes mellitus/hypertension: a molecular study

Background

A paucity of studies focused on the genetic association that tuberculosis (TB) patients with non-communicable diseases (NCDs) are more likely to be infected with Mycobacterium tuberculosis (MTB) with more potent virulence on anti-TB drug resistance than those without NCDs. The study aimed to document the predominant genotype, determine the association between MTB genotypes and NCD status and drug resistance.

Methods

We conducted a molecular study in 105 TB patients based on a cross-sectional study focused on the comorbid relationship between chronic conditions and TB among 1773 subjects from September 1, 2019 to August 30, 2020 in Guizhou, China. The participants were investigated through face-to-face interviews, followed by NCDs screening. The DNA of MTB isolates was extracted prior to genotyping using 24 loci MIRU-VNTR. The subsequent evaluations were performed by phylogenetic trees, combined with tests of statistical power, Chi-square or Fisher and multivariate logistic regression analysis.

Results

The Beijing family of Lineage 2 (East Asia) was the predominant genotype accounting for 43.8% (46/105), followed by Lineage 4 (Euro-America) strains, including Uganda I (34.3%, 36/105), and the NEW-1 (9.5%, 10/105). The proportion of Beijing strain in patients with and without NCDS was 28.6% (8/28) and 49.4% (38/77), respectively, with a statistical power test value of 24.3%. No significant association was detected between MTB genotype and NCD status. A low clustering rate (2.9%) was identified, consisting of two clusters. The rates of global, mono-, poly- and multi-drug resistance were 16.2% (17/105), 14.3% (15/105), 1.0% (1/105) and 4.8% (5/105), respectively. The drug-resistant rates of rifampicin, isoniazid, and streptomycin, were 6.7% (7/105), 11.4% (12/105) and 5.7% (6/105), respectively. Isoniazid resistance was significantly associated with the Beijing genotype of Lineage 2 (19.6% versus 5.1%).

Conclusions

The Lineage 2 East Asia/Beijing genotype is the dominant genotype of the local MTB with endogenous infection preponderating. Not enough evidence is detected to support the association between the MTB genotype and diabetes/hypertension. Isoniazid resistance is associated with the Lineage 2 East Asia/Beijing strain.

Supplementary information

The online version contains supplementary material available at 10.1186/s12879-022-07344-z.