Whole genome sequencing of clinical strains of Mycobacterium tuberculosis from Mumbai, India: A potential tool for determining drug-resistance and strain lineage

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.

An integrated method for targeted Oxford Nanopore sequencing and automated bioinformatics for the simultaneous detection of bacteria, fungi, and ARG

AIMS

The use of metagenomics for pathogen identification in clinical practice has been limited. Here we describe a workflow to encourage the clinical utility and potential of NGS for the screening of bacteria, fungi, and antimicrobial resistance genes (ARGs).

METHODS AND RESULTS

The method includes target enrichment, long-read sequencing, and automated bioinformatics. Evaluation of several tools and databases was undertaken across standard organisms (n = 12), clinical isolates (n = 114), and blood samples from patients with suspected bloodstream infections (n = 33). The strategy used could offset the presence of host background DNA, error rates of long-read sequencing, and provide accurate and reproducible detection of pathogens. Eleven targets could be successfully tested in a single assay. Organisms could be confidently identified considering ≥60% of best hits of a BLAST-based threshold of e-value 0.001 and a percent identity of >80%. For ARGs, reads with percent identity of >90% and >60% overlap of the complete gene could be confidently annotated. A kappa of 0.83 was observed compared to standard diagnostic methods. Thus, a workflow for the direct-from-sample, on-site sequencing combined with automated genomics was demonstrated to be reproducible.

CONCLUSION

NGS-based technologies overcome several limitations of current day diagnostics. Highly sensitive and comprehensive methods of pathogen screening are the need of the hour. We developed a framework for reliable, on-site, screening of pathogens.

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

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.

Prevalence of Mycobacterium tuberculosis mutations associated with isoniazid and rifampicin resistance: A systematic review and meta-analysis

Tuberculosis (TB) is still one of the leading causes of worldwide death, especially following the emergence of strains resistant to isoniazid (INH) and rifampicin (RIF). This study aimed to systematically review published articles focusing on the prevalence of INH and/or RIF resistance-associated mutations of Mycobacterium tuberculosis isolates in recent years. Literature databases were searched using appropriate keywords. The data of the included studies were extracted and used for a random-effects model meta-analysis. Of the initial 1442 studies, 29 were finally eligible to be included in the review. The overall resistance to INH and RIF was about 17.2% and 7.3%, respectively. There was no difference between the frequency of INH and RIF resistance using different phenotypic or genotypic methods. The INH and/or RIF resistance was higher in Asia. The S315T mutation in KatG (23.7 %), C-15 T in InhA (10.7 %), and S531L in RpoB (13.5 %) were the most prevalent mutations. Altogether, the results showed that due to S531L in RpoB, S315T in KatG, and C-15 T in InhA mutations INH- and RIF-resistant M. tuberculosis isolates were widely distributed. Thus, it would be diagnostically and epidemiologically beneficial to track these gene mutations among resistant isolates.

Do private providers initiate anti-tuberculosis therapy on the basis of chest radiographs? A standardised patient study in urban India

Background

The initiation of anti-tuberculosis treatment (ATT) based on results of WHO-approved microbiological diagnostics is an important marker of quality tuberculosis (TB) care. Evidence suggests that other diagnostic processes leading to treatment initiation may be preferred in high TB incidence settings. This study examines whether private providers start anti-TB therapy on the basis of chest radiography (CXR) and clinical examinations.

Methods

This study uses the standardized patient (SP) methodology to generate accurate and unbiased estimates of private sector, primary care provider practice when a patient presents a standardized TB case scenario with an abnormal CXR. Using multivariate log-binomial and linear regressions with standard errors clustered at the provider level, we analyzed 795 SP visits conducted over three data collection waves from 2014 to 2020 in two Indian cities. Data were inverse-probability-weighted based on the study sampling strategy, resulting in city-wave-representative results.

Findings

Amongst SPs who presented to a provider with an abnormal CXR, 25% (95% CI: 21-28%) visits resulted in ideal management, defined as the provider prescribing a microbiological test and not offering a concurrent prescription for a corticosteroid or antibiotic (including anti-TB medications). In contrast, 23% (95% CI: 19-26%) of 795 visits were prescribed anti-TB medications. Of 795 visits, 13% (95% CI: 10-16%) resulted in anti-TB treatment prescriptions/dispensation and an order for confirmatory microbiological testing.

Interpretation

One in five SPs presenting with abnormal CXR were prescribed ATT by private providers. This study contributes novel insights to empiric treatment prevalence based on CXR abnormality. Further work is needed to understand how providers make trade-offs between existing diagnostic practices, new technologies, profits, clinical outcomes, and the market dynamics with laboratories.

Funding

This study was funded by the Bill & Melinda Gates Foundation (grant OPP1091843), and the Knowledge for Change Program at The World Bank.

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.

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.

Whole-genome sequencing of presumptive MDR-TB isolates from a tertiary healthcare setting in Mumbai

OBJECTIVES

Whole-genome sequencing (WGS) of Mycobacterium tuberculosis (MTB), proven to be a better alternative when compared with the combined sensitivity and specificity of all other modalities for diagnosis of tuberculosis (TB), aids epidemiological surveillance investigations by combining the current research with diagnostics. This study was conducted to identify and resolve operational challenges in performing WGS-based drug resistance testing (DRT) for MTB in a TB culture and drug susceptibility testing (DST) laboratory. Three critical, non-redundant steps for WGS-based DRT were tested: viz. DNA extraction, high-throughput paired-end next-generation sequencing (NGS), and genomic analysis pipeline for automated reporting of WGS-based DRT.

METHODS

DNA was extracted from 100 liquid culture isolates on a mycobacterial growth indicator tube (MGIT) using DNEASY Ultraclean Microbial Kit (Qiagen, USA) as per the manufacturer's instructions. Illumina paired-end sequencing was performed. All analysis steps were automated using custom python scripts, requiring no intervention. Variant calling was performed as per the World Health Organization (WHO) technical guide.

RESULTS

The number of cultures resistant to rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin was 89, 88, 35, 67, and 73, respectively. Resistance to amikacin, kanamycin, and capreomycin was found in 15, 17, and 15 cultures, respectively. Seventy cultures were resistant to fluoroquinolones, four were resistant to ethionamide, and 12 were resistant to linezolid. Six cultures were resistant to only one of the 18 drugs tested. Seventy-five cultures were resistant to more than three anti-TB drugs. One culture was resistant to 13 of the 18 anti-TB drugs tested for this study. The maximum number of variants were observed in the rpoB gene (n = 93, 93%), wherein the Ser450Leu was the predominant mutation (n = 68, 73%). Ser315Thr was the most common variant (n = 86, 97%) that encoded resistance to isoniazid. The Lys43Arg variant encodes resistance to streptomycin and was the third most predominant variant (n = 65, 89%). In addition to the high levels of resistance observed in the dataset, we also observed a high proportion of Beijing strains (n = 63, 63%).

CONCLUSION

Compared with results from routine diagnostics based on the 'Guidelines on Programmatic Management of Drug-Resistant TB (PMDT) in India', none of the samples had DST available for all 18 drugs. This represents a gap in PMDT guidelines. The WGS-DRT must be considered as the primary DST method after a sample is flagged rifampicin-resistant by cartridge-based nucleic acid amplification testing (CBNAAT). With several research studies currently underway globally to identify novel variants associated with drug resistance and classifiy their minimum inhibitory coefficients, WGS-DRT presents a scalable technology that updates analytical pipelines, relegating the need for changing microbiological protocols.

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.

Differential Impact of the rpoB Mutant on Rifampin and Rifabutin Resistance Signatures of Mycobacterium tuberculosis Is Revealed Using a Whole-Genome Sequencing Assay

Drug resistance in Mycobacterium tuberculosis (MTB) has long been a serious health issue worldwide. Most drug-resistant MTB isolates were identified due to treatment failure or in clinical examinations 3~6 months postinfection. In this study, we propose a whole-genome sequencing (WGS) pipeline via the Nanopore MinION platform to facilitate the efficacy of phenotypic identification of clinical isolates. We used the Nanopore MinION platform to perform WGS of clinical MTB isolates, including susceptible (n = 30) and rifampin- (RIF) or rifabutin (RFB)-resistant isolates (n = 20) according to results of a susceptibility test. Nonsynonymous variants within the rpoB gene associated with RIF resistance were identified using the WGS analytical pipeline. In total, 131 variants within the rpoB gene in RIF-resistant isolates were identified. The presence of the emergent Asp531Gly or His445Gln was first identified to be associated with the rifampin and rifabutin resistance signatures of clinical isolates. The results of the minimum inhibitory concentration (MIC) test further indicated that the Ser450Leu or the mutant within the rifampin resistance-determining region (RRDR)-associated rifabutin-resistant signature was diminished in the presence of novel mutants, including Phe669Val, Leu206Ile, or Met148Leu, identified in this study.

IMPORTANCE Current approaches to diagnose drug-resistant MTB are time-consuming, consequently leading to inefficient intervention or further disease transmission. In this study, we curated lists of coding variants associated with differential rifampin and rifabutin resistant signatures using a single molecule real-time (SMRT) sequencing platform with a shorter hands-on time. Accordingly, the emerging WGS pipeline constitutes a potential platform for efficacious and accurate diagnosis of drug-resistant MTB isolates.

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.

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.

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.

Whole-genome sequencing and Mycobacterium tuberculosis: Challenges in sample preparation and sequencing data analysis

The numbers of patients with tuberculosis (TB) caused by resistant strains are still alarming. Therefore, it is necessary to determine resistance more quickly and precisely, than it is with the currently used phenotypic and genotypic methods. In recent years, technological advances have been made and the whole-genome sequencing (WGS) method has been introduced as a part of routine diagnostics in clinical laboratories. Comparing a wide range of mycobacterial genomic variations with a reference genome leads to a consistent evaluation of molecular-epidemiology and resistance of Mycobacterium tuberculosis (M. tuberculosis) to a wide range of anti-TB drugs. The quality of the obtained sequencing data is closely related to the type of sample and the method used for DNA extraction and sequencing library preparation. Moreover, the correct interpretation of results is also influenced by a bioinformatic data processing. A large number of bioinformatics pipelines are currently available, the sensitivity of which varies due to the different sizes of databases containing relevant mutations. This review focuses on the individual steps included in the sequencing workflow and factors that may affect the interpretation of final results.

Molecular characterisation of multidrug-resistant Mycobacterium tuberculosis isolates from a high-burden tuberculosis state in Brazil

Tuberculosis (TB) is the leading cause of death among infectious diseases worldwide. Among the estimated cases of drug-resistant TB, approximately 60% occur in the BRICS countries (Brazil, Russia, India, China and South Africa). Among Brazilian states, primary and acquired multidrug-resistant TB (MDR-TB) rates were the highest in Rio Grande do Sul (RS). This study aimed to perform molecular characterisation of MDR-TB in the State of RS, a high-burden Brazilian state. We performed molecular characterisation of MDR-TB cases in RS, defined by drug susceptibility testing, using 131 Mycobacterium tuberculosis (M.tb) DNA samples from the Central Laboratory. We carried out MIRU-VNTR 24loci, spoligotyping, sequencing of the katG, inhA and rpoB genes and RD^Rio sublineage identification. The most frequent families found were LAM (65.6%) and Haarlem (22.1%). RD^Rio deletion was observed in 42 (32%) of the M.tb isolates. Among MDR-TB cases, eight (6.1%) did not present mutations in the studied genes. In 116 (88.5%) M.tb isolates, we found mutations associated with rifampicin (RIF) resistance in rpoB gene, and in 112 isolates (85.5%), we observed mutations related to isoniazid resistance in katG and inhA genes. An insertion of 12 nucleotides (CCAGAACAACCC) at the 516 codon in the rpoB gene, possibly responsible for a decreased interaction of RIF and RNA polymerase, was found in 19/131 of the isolates, belonging mostly to LAM and Haarlem families. These results enable a better understanding of the dynamics of transmission and evolution of MDR-TB in the region.

Whole genome sequencing of Mycobacterium tuberculosis isolates and clinical outcomes of patients treated for multidrug-resistant tuberculosis in Tanzania

BACKGROUND

Tuberculosis (TB), particularly multi- and or extensive drug resistant TB, is still a global medical emergency. Whole genome sequencing (WGS) is a current alternative to the WHO-approved probe-based methods for TB diagnosis and detection of drug resistance, genetic diversity and transmission dynamics of Mycobacterium tuberculosis complex (MTBC). This study compared WGS and clinical data in participants with TB.

RESULTS

This cohort study performed WGS on 87 from MTBC DNA isolates, 57 (66%) and 30 (34%) patients with drug resistant and susceptible TB, respectively. Drug resistance was determined by Xpert® MTB/RIF assay and phenotypic culture-based drug-susceptibility-testing (DST). WGS and bioinformatics data that predict phenotypic resistance to anti-TB drugs were compared with participant's clinical outcomes. They were 47 female participants (54%) and the median age was 35 years (IQR): 29-44). Twenty (23%) and 26 (30%) of participants had TB/HIV co-infection BMI < 18 kg/m2 respectively. MDR-TB participants had MTBC with multiple mutant genes, compared to those with mono or polyresistant TB, and the majority belonged to lineage 3 Central Asian Strain (CAS). Also, MDR-TB was associated with delayed culture-conversion (median: IQR (83: 60-180 vs. 51:30-66) days). WGS had high concordance with both culture-based DST and Xpert® MTB/RIF assay in detecting drug resistance (kappa = 1.00).

CONCLUSION

This study offers comparison of mutations detected by Xpert and WGS with phenotypic DST of M. tuberculosis isolates in Tanzania. The high concordance between the different methods and further insights provided by WGS such as PZA-DST, which is not routinely performed in most resource-limited-settings, provides an avenue for inclusion of WGS into diagnostic matrix of TB including drug-resistant TB.

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.

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.

Whole-Genome Sequencing in Relation to Resistance of Mycobacterium Tuberculosis

Tuberculosis, a disease caused by Mycobacterium tuberculosis, represents one of the deadliest infections worldwide. The incidence of resistant forms is increasing year by year; therefore, it is necessary to involve new methods for rapid diagnostics and treatment. One of the possible solutions is the use of whole-genome sequencing (WGS).

The WGS provides an identification of complete genome of the microorganism, including all genes responsible for resistance, in comparison with other genotypic methods (eg. Xpert MTB / RIF or Hain line-probes) that are capable to detect only basic genes. WGS data are available in 1-9 days and several online software tools (TBProfiler, CASTB, Mykrobe PredictorTB) are used for their interpretation and analysis, compared to 3-8 weeks in the case of classic phenotypic evaluation.

Furthermore, WGS predicts resistance to the first-line antituberculotics with a sensitivity of 85-100% and a specificity of 85-100%.

This review elucidates the importance and summarizes the current knowledge about the possible use of WGS in diagnosis and treatment of resistant forms of tuberculosis elucidates.

WGS of M. tuberculosis brings new possibilities for rapid and accurate diagnostics of resistant forms of tuberculosis. Introducing WGS into routine practice can help to reduce the spread of resistant forms of tuberculosis as well as to increase the success rate of the treatment, especially through an appropriate combination of antituberculotics ATs. Introduction of WGS into routine diagnostics can, in spite of the financial difficulty, significantly improve patient care.

Beyond multidrug resistance: Leveraging rare variants with machine and statistical learning models in Mycobacterium tuberculosis resistance prediction

Background

Methods

Findings

Interpretation

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.

Meta-narrative review of molecular methods for diagnosis and monitoring of multidrug-resistant tuberculosis treatment in adults

Early and accurate diagnosis and rigorous clinical and microbiological monitoring of multidrug-resistant tuberculosis (MDR-TB) treatment can curb morbidity and mortality. While others are still under evaluation, the World Health Organization has recommended few novel molecular methods for MDR-TB diagnosis only. We present current molecular methods for diagnosis and monitoring of MDR-TB treatment in TB-endemic settings. A systematic meta-narrative review was conducted according to the RAMESES recommendations. Electronic databases were searched for relevant articles published in English language from January 2013 to June 2018. Based on predefined criteria, two independent reviewers extracted the key messages from relevant articles. Disagreement between them was resolved through discussion and the involvement of a third reviewer, if needed. Key messages were synthesized to create the meta-narratives for method's accuracy, drug-susceptibility capability, and laboratory infrastructure required. We included 33 articles out of 1213 records retrieved, of which 16 (48%) and 12 (36%) were conducted in high- and low-TB-endemic settings, respectively. Xpert® MTB/RIF, GenoType MTBDRplus, GenoType MTBDRsl, FlouroType™ MTBDR, TB TaqMan® array card, and DNA sequencers can accurately guide effective treatment regimens. Molecular bacterial load assay quantifies mycobactericidal impact of these regimens. Although they present inherent advantages compared to the current standard of care, they carry important limitations to implementation and/or scale-up. Therefore, considerable effort must now be directed to implementation and health systems research to maximize these forecasted benefits for individual patient's health outcomes.

Dissecting whole-genome sequencing-based online tools for predicting resistance in Mycobacterium tuberculosis: can we use them for clinical decision guidance?

Whole-genome sequencing (WGS)-based bioinformatics platforms for the rapid prediction of resistance will soon be implemented in the Tuberculosis (TB) laboratory, but their accuracy assessment still needs to be strengthened. Here, we fully-sequenced a total of 54 multidrug-resistant (MDR) and five susceptible TB strains and performed, for the first time, a simultaneous evaluation of the major four free online platforms (TB Profiler, PhyResSE, Mykrobe Predictor and TGS-TB). Overall, the sensitivity of resistance prediction ranged from 84.3% using Mykrobe predictor to 95.2% using TB profiler, while specificity was higher and homogeneous among platforms. TB profiler revealed the best performance robustness (sensitivity, specificity, PPV and NPV above 95%), followed by TGS-TB (all parameters above 90%). We also observed a few discrepancies between phenotype and genotype, where, in some cases, it was possible to pin-point some "candidate" mutations (e.g., in the rpsL promoter region) highlighting the need for their confirmation through mutagenesis assays and potential review of the anti-TB genetic databases. The rampant development of the bioinformatics algorithms and the tremendously reduced time-frame until the clinician may decide for a definitive and most effective treatment will certainly trigger the technological transition where WGS-based bioinformatics platforms could replace phenotypic drug susceptibility testing for TB.