Cryptic Microheteroresistance Explains Mycobacterium tuberculosis Phenotypic Resistance

Xpert MTB/RIF Ultra-resistant and MTBDRplus-susceptible rifampicin results in people with tuberculosis: utility of FluoroType MTBDR and deep sequencing

Xpert MTB/RIF Ultra (Ultra)-detected rifampicin-resistant tuberculosis (TB) is often programmatically confirmed using MTBDRplus. There are limited data on discordant results, including when re-tested using newer methods, like FluoroType MTBDR (FT-MTBDR) and targeted deep sequencing. MTBDRplus rifampicin-susceptible isolates from people with Ultra rifampicin-resistant sputum were identified from a South African programmatic laboratory. FT-MTBDR and single molecule-overlapping reads (SMOR; rpoB, inhA, katG) on isolate DNA were done (SMOR was used as a reference standard). Between 1 April 2021 and 30 September 2022, 8% (109/1347) of Ultra rifampicin-resistant specimens were MTBDRplus-susceptible. Of 89% (97/109) isolates with a sequenceable rpoB, SMOR resolved most in favor of Ultra (79% [77/97]). Sputum with lower mycobacterial load was associated with Ultra false-positive resistance (46% [11/24] of "very low" Ultra had false resistance vs 12% [9/73; P = 0.0004] of ≥"low"), as were Ultra heteroresistance calls (all wild-type probes, ≥1 mutant probe) (62% [23/37 vs 25% 15/60] for Ultra without heteroresistance calls; P = 0.0003). Of the 91% (88/97) of isolates successfully tested by FT-MTBDR, 55% (48/88) were FT-MTBDR rifampicin-resistant and 45% (40/88) susceptible, translating to 69% (47/68) sensitivity and 95% (19/20) specificity. In the 91% (99/109) of isolates with inhA and katG sequenced, 62% (61/99) were SMOR isoniazid-susceptible. When Ultra and MTBDRplus rifampicin results are discordant, Ultra is more likely to be correct, and FT-MTBDR agrees more with Ultra than MTBDRplus; however, lower load and the Ultra heteroresistance probe pattern were risk factors for Ultra false rifampicin-resistant results. Most people with Ultra-MTBDRplus discordant resistance results were isoniazid-susceptible. These data have implications for drug-resistant TB diagnosis.

Sequencing by binding rivals SMOR error-corrected sequencing by synthesis technology for accurate detection and quantification of minor (< 0.1%) subpopulation variants

BACKGROUND

Detecting very minor (< 1%) subpopulations using next-generation sequencing is a critical need for multiple applications, including the detection of drug resistant pathogens and somatic variant detection in oncology. A recently available sequencing approach termed 'sequencing by binding (SBB)' claims to have higher base calling accuracy data "out of the box." This paper evaluates the utility of using SBB for the detection of ultra-rare drug resistant subpopulations in Mycobacterium tuberculosis (Mtb) using a targeted amplicon assay and compares the performance of SBB to single molecule overlapping reads (SMOR) error corrected sequencing by synthesis (SBS) data.

RESULTS

SBS displayed an elevated error rate when compared to SMOR error-corrected SBS and SBB techniques. SMOR error-corrected SBS and SBB technologies performed similarly within the linear range studies and error rate studies.

CONCLUSIONS

With lower sequencing error rates within SBB sequencing, this technique looks promising for both targeted and unbiased whole genome sequencing, leading to the identification of minor (< 1%) subpopulations without the need for error correction methods.

Sequencing by Binding rivals error-corrected Sequencing by Synthesis technology for accurate detection and quantification of minor (<0.1%) subpopulation variants

Detecting very minor (< 1%) subpopulations using next-generation sequencing is a critical need for multiple applications including detection of drug resistant pathogens and somatic variant detection in oncology. To enable these applications, wet lab enhancements and bioinformatic error correction methods have been developed for 'sequencing by synthesis' technology to reduce its inherent sequencing error rate. A recently available sequencing approach termed 'sequencing by binding' claims to have higher base calling accuracy data "out of the box." This paper evaluates the utility of using 'sequencing by binding' for the detection of ultra-rare subpopulations down to 0.001%.

Genomic markers of drug resistance in Mycobacterium tuberculosis populations with minority variants

Minority variants of Mycobacterium tuberculosis harboring mutations conferring resistance can become dominant populations during tuberculosis (TB) treatment, leading to treatment failure. Our understanding of drug-resistant within-host subpopulations and the frequency of resistance-conferring mutations in minority variants remains limited. M. tuberculosis sequences recovered from liquid cultures of culture-confirmed TB cases notified between January 2017 and December 2021 in New South Wales, Australia were examined. Potential drug resistance-conferring minority variants were identified using LoFreq, and mixed populations of different M. tuberculosis strains (≥100 SNPs apart) were examined using QuantTB. A total of 1831 routinely sequenced M. tuberculosis strains were included in the analysis. Drug resistance-conferring minority variants were detected in 3.5% (65/1831) of sequenced cultures; 84.6% (55/65) had majority strains that were drug susceptible and 15.4% (10/65) had majority strains that were drug resistant. Minority variants with high-confidence drug resistance-conferring mutations were 1.5 times more common when the majority strains were drug resistant. Mixed M. tuberculosis strain populations were documented in 10.0% (183/1831) of specimens. Minority variants with high-confidence drug resistance-conferring mutations were more frequently detected in mixed M. tuberculosis strain populations (2.7%, 5/183) than in single strain populations (0.6%, 10/1648; P = 0.01). Drug-resistant minority variants require monitoring in settings that implement routine M. tuberculosis sequencing. The frequency with which drug-resistant minority variants are detected is likely influenced by pre-culture requirement. Culture-independent sequencing methods should provide a more accurate reflection of drug-resistant subpopulations.

Clinical implications of molecular drug resistance testing for Mycobacterium tuberculosis: a 2023 TBnet/RESIST-TB consensus statement

Drug-resistant tuberculosis is a substantial health-care concern worldwide. Despite culture-based methods being considered the gold standard for drug susceptibility testing, molecular methods provide rapid information about the Mycobacterium tuberculosis mutations associated with resistance to anti-tuberculosis drugs. This consensus document was developed on the basis of a comprehensive literature search, by the TBnet and RESIST-TB networks, about reporting standards for the clinical use of molecular drug susceptibility testing. Review and the search for evidence included hand-searching journals and searching electronic databases. The panel identified studies that linked mutations in genomic regions of M tuberculosis with treatment outcome data. Implementation of molecular testing for the prediction of drug resistance in M tuberculosis is key. Detection of mutations in clinical isolates has implications for the clinical management of patients with multidrug-resistant or rifampicin-resistant tuberculosis, especially in situations when phenotypic drug susceptibility testing is not available. A multidisciplinary team including clinicians, microbiologists, and laboratory scientists reached a consensus on key questions relevant to molecular prediction of drug susceptibility or resistance to M tuberculosis, and their implications for clinical practice. This consensus document should help clinicians in the management of patients with tuberculosis, providing guidance for the design of treatment regimens and optimising outcomes.

Detection of minor variants in Mycobacterium tuberculosis whole genome sequencing data

The study of genetic minority variants is fundamental to the understanding of complex processes such as evolution, fitness, transmission, virulence, heteroresistance and drug tolerance in Mycobacterium tuberculosis (Mtb). We evaluated the performance of the variant calling tool LoFreq to detect de novo as well as drug resistance conferring minor variants in both in silico and clinical Mtb next generation sequencing (NGS) data. The in silico simulations demonstrated that LoFreq is a conservative variant caller with very high precision (≥96.7%) over the entire range of depth of coverage tested (30x to1000x), independent of the type and frequency of the minor variant. Sensitivity increased with increasing depth of coverage and increasing frequency of the variant, and was higher for calling insertion and deletion (indel) variants than for single nucleotide polymorphisms (SNP). The variant frequency limit of detection was 0.5% and 3% for indel and SNP minor variants, respectively. For serial isolates from a patient with DR-TB; LoFreq successfully identified all minor Mtb variants in the Rv0678 gene (allele frequency as low as 3.22% according to targeted deep sequencing) in whole genome sequencing data (median coverage of 62X). In conclusion, LoFreq can successfully detect minor variant populations in Mtb NGS data, thus limiting the need for filtering of possible false positive variants due to sequencing error. The observed performance statistics can be used to determine the limit of detection in existing whole genome sequencing Mtb data and guide the required depth of future studies that aim to investigate the presence of minor variants.

Comparative Performance of Genomic Methods for the Detection of Pyrazinamide Resistance and Heteroresistance in Mycobacterium tuberculosis

Pyrazinamide is an important component of both drug-susceptible and drug-resistant tuberculosis treatment regimens. Although approximately 50% of rifampin-resistant isolates are also resistant to pyrazinamide, pyrazinamide susceptibility testing is not routinely performed due to the challenging nature of the assay. We investigated the diagnostic accuracy of genotypic and phenotypic methods and explored the occurrence of pyrazinamide heteroresistance. We assessed pyrazinamide susceptibility among 358 individuals enrolled in the South African EXIT-RIF cohort using Sanger and targeted deep sequencing (TDS) of the pncA gene, whole-genome sequencing (WGS), and phenotypic drug susceptibility testing. We calculated the diagnostic accuracy of the different methods and investigated the prevalence and clinical impact of pncA heteroresistance. True pyrazinamide susceptibility status was assigned to each isolate using the Köser classification and expert rules. We observed 100% agreement across genotypic methods for detection of pncA fixed mutations; only TDS confidently identified three isolates (0.8%) with minor variants. For the 355 (99.2%) isolates that could be assigned true pyrazinamide status with confidence, phenotypic DST had a sensitivity of 96.5% (95% confidence interval [CI], 93.8 to 99.3%) and specificity of 100% (95% CI, 100 to 100%), both Sanger sequencing and WGS had a sensitivity of 97.1% (95% CI, 94.6 to 99.6%) and specificity of 97.8% (95% CI, 95.7 to 99.9%), and TDS had sensitivity of 98.8% (95% CI, 97.2 to 100%) and specificity of 97.8% (95% CI, 95.7 to 99.9%). We demonstrate high sensitivity and specificity for pyrazinamide susceptibility testing among all assessed genotypic methods. The prevalence of pyrazinamide heteroresistance in Mycobacterium tuberculosis isolates was lower than that identified for other first-line drugs.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

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

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

Resistance-Conferring Mutations on Whole-Genome Sequencing of Fluoroquinolone-resistant and -Susceptible Mycobacterium tuberculosis Isolates: A Proposed Threshold for Identifying Resistance

BACKGROUND

Fluoroquinolone resistance in Mycobacterium tuberculosis (Mtb) is conferred by DNA gyrase mutations, but not all fluoroquinolone-resistant Mtb isolates have mutations detected. The optimal allele frequency threshold to identify resistance-conferring mutations by whole-genome sequencing is unknown.

METHODS

Phenotypically ofloxacin-resistant and lineage-matched ofloxacin-susceptible Mtb isolates underwent whole-genome sequencing at an average coverage depth of 868 reads. Polymorphisms within the quinolone-resistance-determining region (QRDR) of gyrA and gyrB were identified. The allele frequency threshold using the Genome Analysis Toolkit pipeline was ~8%; allele-level data identified the predominant variant allele frequency and mutational burden (ie, sum of all variant allele frequencies in the QRDR) in gyrA, gyrB, and gyrA + gyrB for each isolate. Receiver operating characteristic (ROC) curves assessed the optimal measure of allele frequency and potential thresholds for identifying phenotypically resistant isolates.

RESULTS

Of 42 ofloxacin-resistant Mtb isolates, area under the ROC curve (AUC) was highest for predominant variant allele frequency, so that measure was used to evaluate optimal mutation detection thresholds. AUCs for 8%, 2.5%, and 0.8% thresholds were 0.8452, 0.9286, and 0.9069, respectively. Sensitivity and specificity were 69% and 100% for 8%, 86% and 100% for 2.5%, 91% and 91% for 0.8%. The sensitivity of the 2.5% and 0.8% thresholds were significantly higher than the 8% threshold (P = .016 and .004, respectively) but not significantly different between one another (P = .5).

CONCLUSIONS

A predominant mutation allele frequency threshold of 2.5% had the highest AUC for detecting DNA gyrase mutations that confer ofloxacin resistance, and was therefore the optimal threshold.

Deep amplicon sequencing for culture-free prediction of susceptibility or resistance to 13 anti-tuberculous drugs

Conventional molecular tests for detecting Mycobacterium tuberculosis complex (MTBC) drug resistance on clinical samples cover a limited set of mutations. Whole-genome sequencing (WGS) typically requires culture.

Here, we evaluated the Deeplex Myc-TB targeted deep-sequencing assay for prediction of resistance to 13 anti-tuberculous drugs/drug classes, directly applicable on sputum.

With MTBC DNA tests, the limit of detection was 100–1000 genome copies for fixed resistance mutations. Deeplex Myc-TB captured in silico 97.1–99.3% of resistance phenotypes correctly predicted by WGS from 3651 MTBC genomes. On 429 isolates, the assay predicted 92.2% of 2369 first- and second-line phenotypes, with a sensitivity of 95.3% and a specificity of 97.4%. 56 out of 69 (81.2%) residual discrepancies with phenotypic results involved pyrazinamide, ethambutol and ethionamide, and low-level rifampicin or isoniazid resistance mutations, all notoriously prone to phenotypic testing variability. Only two out of 91 (2.2%) resistance phenotypes undetected by Deeplex Myc-TB had known resistance-associated mutations by WGS analysis outside Deeplex Myc-TB targets. Phenotype predictions from Deeplex Myc-TB analysis directly on 109 sputa from a Djibouti survey matched those of MTBSeq/PhyResSE/Mykrobe, fed with WGS data from subsequent cultures, with a sensitivity of 93.5/98.5/93.1% and a specificity of 98.5/97.2/95.3%, respectively. Most residual discordances involved gene deletions/indels and 3–12% heteroresistant calls undetected by WGS analysis or natural pyrazinamide resistance of globally rare “Mycobacterium canettii” strains then unreported by Deeplex Myc-TB. On 1494 arduous sputa from a Democratic Republic of the Congo survey, 14 902 out of 19 422 (76.7%) possible susceptible or resistance phenotypes could be predicted culture-free.

Deeplex Myc-TB may enable fast, tailored tuberculosis treatment.

The novel Deeplex Myc-TB molecular assay shows a high degree of accuracy for extensive prediction of susceptibility and resistance to 13 anti-tuberculous drugs, directly achievable without culture, which may enable fast, tailored tuberculosis treatmenthttps://bit.ly/3bAvcAt

Linezolid resistance in patients with drug-resistant TB and treatment failure in South Africa

Objectives

Limited data exist on clinical associations and genotypic correlates of linezolid resistance in Mycobacterium tuberculosis. We aimed to describe mutations and clinical factors associated with phenotypic linezolid resistance from patients with drug-resistant TB at two public sector facilities in South Africa.

Methods

Adults and adolescents with treatment failure (culture positivity ≥4 months) on a linezolid-containing regimen were retrospectively identified. Phenotypic resistance, as defined by a linezolid MIC >1 mg/L, was assessed for retrieved isolates using broth microdilution. Targeted sequencing of rrl and rplC was performed, irrespective of growth on subculture.

Results

Thirty-nine patients with linezolid-based treatment failure were identified, 13 (33%) of whom had phenotypic or genotypic linezolid resistance after a median duration of 22 months (range = 7–32) of linezolid therapy. Paired MIC testing and genotyping was performed on 55 unique isolates. All isolates with phenotypic resistance (n = 16) were associated with known resistance mutations, most frequently due to the T460C substitution in rplC (n = 10); rrl mutations included G2814T, G2270C/T and A2810C. No mutations were detected in isolates with MICs at or below the critical concentration.

Conclusions

Linezolid resistance occurred in a third of patients with drug-resistant TB and treatment failure. Resistance occurred late and was predicted by a limited number of mutations in rrl and rplC. Screening for genotypic resistance should be considered for patients with a positive culture after 4 months of linezolid therapy in order to optimize treatment and avoid the toxicity of ineffective linezolid therapy.

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.

Genomic signatures of drug resistance in highly resistant Mycobacterium tuberculosis strains of the early ancient sublineage of Beijing genotype in Russia

The Mycobacterium tuberculosis Beijing genotype is a clinically and epidemiologically important lineage that is subdivided into ancient/ancestral and modern strains. In our previous study in western Siberia, we identified variable number of tandem repeats (VNTR)-based clusters within the early ancient sublineage of the Beijing genotype characterized by an unexpectedly high rate of extensive drug resistance (XDR). In the current study, next generation sequencing data were analysed to gain insight into genomic signatures underlying drug resistance of these strains. A total of 184 genomes of the Beijing early ancient sublineage from Russia (16), China (15), Japan (36), Korea (25), Vietnam (18), Thailand (73), and the USA (1) were used for phylogenetic analysis. The drug-resistant profile was deduced genotypically. The Russian isolates were distributed into two clusters and were all drug resistant, mainly pre-XDR and XDR. The largest of these clusters included only Russian isolates from remote locations in both Asian and European parts of the country. All its isolates had a quadruple drug resistance (to isoniazid, rifampin, ethambutol and streptomycin) due to the 6-mutation signature (KatG Ser315Thr, KatG Ile335Val, RpoB Ser450Leu, RpoC Asp485Asn, EmbB Gln497Arg, and RpsL Lys43Arg). In most samples, it was complemented with additional and different pncA, gyrA and rrs mutations leading to the pre-XDR/XDR genotype. Phylogenomic analysis indicates a distant origin of this Russian resistant cluster in the early 1970s but location and circumstances are yet to be clarified.

Detection of low-frequency resistance-mediating SNPs in next-generation sequencing data of Mycobacterium tuberculosis complex strains with binoSNP

Accurate drug resistance detection is key for guiding effective tuberculosis treatment. While genotypic resistance can be rapidly detected by molecular methods, their application is challenged by mixed mycobacterial populations comprising both susceptible and resistant cells (heteroresistance). For this, next-generation sequencing (NGS) based approaches promise the determination of variants even at low frequencies. However, accurate methods for a valid detection of low-frequency variants in NGS data are currently lacking. To tackle this problem, we developed the variant detection tool binoSNP which allows the determination of low-frequency single nucleotide polymorphisms (SNPs) in NGS datasets from Mycobacterium tuberculosis complex (MTBC) strains. By taking a reference-mapped file as input, binoSNP evaluates each genomic position of interest using a binomial test procedure. binoSNP was validated using in-silico, in-vitro, and serial patient isolates datasets comprising varying genomic coverage depths (100-500×) and SNP allele frequencies (1-30%). Overall, the detection limit for low-frequency SNPs depends on the combination of coverage depth and allele frequency of the resistance-associated mutation. binoSNP allows for valid detection of resistance associated SNPs at a 1% frequency with a coverage ≥400×. In conclusion, binoSNP provides a valid approach to detect low-frequency resistance-mediating SNPs in NGS data from clinical MTBC strains. It can be implemented in automated, end-user friendly analysis tools for NGS data and is a step forward towards individualized TB therapy.

Dynamics of within-host Mycobacterium tuberculosis diversity and heteroresistance during treatment

BACKGROUND

Studying within-host genetic diversity of Mycobacterium tuberculosis (Mtb) in patients during treatment may identify adaptations to antibiotic and immune pressure. Understanding the significance of genetic heteroresistance, and more specifically heterozygous resistance-associated variants (RAVs), is clinically important given increasing use of rapid molecular tests and whole genome sequencing (WGS).

METHODS

We analyse data from six studies in KwaZulu-Natal, South Africa. Most patients (>75%) had baseline rifampicin resistance. Sputum was collected for culture at baseline and at between two and nine intervals until month six. Positive cultures underwent WGS. Mixed infections and reinfections were excluded from analysis.

FINDINGS

Baseline Mtb overall genetic diversity (at treatment initiation or major change to regimen) was associated with cavitary disease, not taking antiretroviral therapy if HIV infected, infection with lineage 2 strains and absence of second-line drug resistance on univariate analyses. Baseline genetic diversity was not associated with six-month outcome. Genetic diversity increased from baseline to weeks one and two before returning to previous levels. Baseline genetic heteroresistance was most common for bedaquiline (6/10 [60%] of isolates with RAVs) and fluoroquinolones (9/62 [13%]). Most patients with heterozygous RAVs on WGS with sequential isolates available demonstrated RAV persistence or fixation (17/20, 85%). New RAVs emerged in 9/286 (3%) patients during treatment. We could detect low-frequency RAVs preceding emergent resistance in only one case, although validation of deep sequencing to detect rare variants is required.

INTERPRETATION

In this study of single-strain Mtb infections, baseline within-host bacterial genetic diversity did not predict outcome but may reveal adaptations to host and drug pressures. Predicting emergent resistance from low-frequency RAVs requires further work to separate transient from consequential mutations.

FUNDING

Wellcome Trust, NIH/NIAID.

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.

How Well Do Routine Molecular Diagnostics Detect Rifampin Heteroresistance in Mycobacterium tuberculosis?

Rifampin heteroresistance—where rifampin-resistant and -susceptible tuberculosis (TB) bacilli coexist—may result in failed standard TB treatment and potential spread of rifampin-resistant strains. The detection of rifampin heteroresistance in routine rapid diagnostic tests (RDTs) allows for patients to receive prompt and effective multidrug-resistant-TB treatment and may improve rifampin-resistant TB control.

Rifampin heteroresistance—where rifampin-resistant and -susceptible tuberculosis (TB) bacilli coexist—may result in failed standard TB treatment and potential spread of rifampin-resistant strains. The detection of rifampin heteroresistance in routine rapid diagnostic tests (RDTs) allows for patients to receive prompt and effective multidrug-resistant-TB treatment and may improve rifampin-resistant TB control. The limit of detection (LOD) of rifampin heteroresistance for phenotypic drug susceptibility testing by the proportion method is 1% and, yet, is insufficiently documented for RDTs. We, therefore, aimed to determine, for the four RDTs (XpertMTB/RIF, XpertMTB/RIF Ultra, GenoTypeMTBDRplusv2.0, and GenoscholarNTM+MDRTBII), the LOD per probe and mutation, validated by CFU counting and targeted deep sequencing (Deeplex-MycTB). We selected one rifampin-susceptible and four rifampin-resistant strains, with mutations D435V, H445D, H445Y, and S450L, respectively, mixed them in various proportions in triplicate, tested them with each RDT, and determined the LODs per mutation type. Deeplex-MycTB revealed concordant proportions of the minority resistant variants in the mixtures. The Deeplex-MycTB-validated LODs ranged from 20% to 80% for XpertMTB/RIF, 20% to 70% for Xpert Ultra, 5% to 10% for GenoTypeMTBDRplusv2.0, and 1% to 10% for GenoscholarNTM+MDRTBII for the different mutations. Deeplex-MycTB, GenoTypeMTBDRplusv2.0, and GenoscholarNTM+MDRTBII provide explicit information on rifampin heteroresistance for the most frequently detected mutations. Classic Xpert and Ultra report rifampin heteroresistance as rifampin resistance, while Ultra may denote rifampin heteroresistance through “mixed patterns” of wild-type and mutant melt probe, melt peak temperatures. Overall, our findings inform end users that the threshold for reporting resistance in the case of rifampin heteroresistance is the highest for Classic Xpert and Ultra to resolve phenotypic and genotypic discordant rifampin-resistant TB results.

Management of drug-resistant tuberculosis

Drug-resistant tuberculosis is a major public health concern in many countries. Over the past decade, the number of patients infected with Mycobacterium tuberculosis resistant to the most effective drugs against tuberculosis (ie, rifampicin and isoniazid), which is called multidrug-resistant tuberculosis, has continued to increase. Globally, 4·6% of patients with tuberculosis have multidrug-resistant tuberculosis, but in some areas, like Kazakhstan, Kyrgyzstan, Moldova, and Ukraine, this proportion exceeds 25%. Treatment for patients with multidrug-resistant tuberculosis is prolonged (ie, 9-24 months) and patients with multidrug-resistant tuberculosis have less favourable outcomes than those treated for drug-susceptible tuberculosis. Individualised multidrug-resistant tuberculosis treatment with novel (eg, bedaquiline) and repurposed (eg, linezolid, clofazimine, or meropenem) drugs and guided by genotypic and phenotypic drug susceptibility testing can improve treatment outcomes. Some clinical trials are evaluating 6-month regimens to simplify management and improve outcomes of patients with multidrug-resistant tuberculosis. Here we review optimal diagnostic and treatment strategies for patients with drug-resistant tuberculosis and their contacts.

Utility of Targeted, Amplicon-Based Deep Sequencing To Detect Resistance to First-Line Tuberculosis Drugs in Botswana

Multidrug-resistant tuberculosis (TB) is an alarming threat, and targeted deep sequencing (DS) may be an effective method for rapid identification of drug-resistant profiles, including detection of heteroresistance. We evaluated the sensitivity and specificity of targeted DS versus phenotypic drug susceptibility testing (pDST) among patients starting first-line anti-TB therapy in Botswana. Overall, we found high concordance between DS and pDST. Lower sensitivity of DS, which targets established high-confidence resistance variants, was observed for detecting isoniazid resistance among HIV-infected patients.

Whole genome sequencing Mycobacterium tuberculosis directly from sputum identifies more genetic diversity than sequencing from culture

Background

Repeated culture reduces within-sample Mycobacterium tuberculosis genetic diversity due to selection of clones suited to growth in culture and/or random loss of lineages, but it is not known to what extent omitting the culture step altogether alters genetic diversity. We compared M. tuberculosis whole genome sequences generated from 33 paired clinical samples using two methods. In one method DNA was extracted directly from sputum then enriched with custom-designed SureSelect (Agilent) oligonucleotide baits and in the other it was extracted from mycobacterial growth indicator tube (MGIT) culture.

Results

DNA directly sequenced from sputum showed significantly more within-sample diversity than that from MGIT culture (median 5.0 vs 4.5 heterozygous alleles per sample, p = 0.04). Resistance associated variants present as HAs occurred in four patients, and in two cases may provide a genotypic explanation for phenotypic resistance.

Conclusions

Culture-free M. tuberculosis whole genome sequencing detects more within-sample diversity than a leading culture-based method and may allow detection of mycobacteria that are not actively replicating.

Electronic supplementary material

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

Transmission of drug-resistant tuberculosis in HIV-endemic settings

The emergence and expansion of the multidrug-resistant tuberculosis epidemic is a threat to the global control of tuberculosis. Multidrug-resistant tuberculosis is the result of the selection of resistance-conferring mutations during inadequate antituberculosis treatment. However, HIV has a profound effect on the natural history of tuberculosis, manifesting in an increased rate of disease progression, leading to increased transmission and amplification of multidrug-resistant tuberculosis. Interventions specific to HIV-endemic areas are urgently needed to block tuberculosis transmission. These interventions should include a combination of rapid molecular diagnostics and improved chemotherapy to shorten the duration of infectiousness, implementation of infection control measures, and active screening of multidrug-resistant tuberculosis contacts, with prophylactic regimens for individuals without evidence of disease. Development and improvement of the efficacy of interventions will require a greater understanding of the factors affecting the transmission of multidrug-resistant tuberculosis in HIV-endemic settings, including population-based molecular epidemiology studies. In this Series article, we review what we know about the transmission of multidrug-resistant tuberculosis in settings with high burdens of HIV and define the research priorities required to develop more effective interventions, to diminish ongoing transmission and the amplification of drug resistance.

Mixed Mycobacterium tuberculosis-Strain Infections Are Associated With Poor Treatment Outcomes Among Patients With Newly Diagnosed Tuberculosis, Independent of Pretreatment Heteroresistance

Background

Heteroresistant Mycobacterium tuberculosis infections (defined as concomitant infection with drug-resistant and drug-susceptible strains) may explain the higher risk of poor tuberculosis treatment outcomes observed among patients with mixed-strain M. tuberculosis infections. We investigated the clinical effect of mixed-strain infections while controlling for pretreatment heteroresistance in a population-based sample of patients with tuberculosis starting first-line tuberculosis therapy in Botswana.

Methods

We performed 24-locus mycobacterial interspersed repetitive unit-variable number tandem-repeat analysis and targeted deep sequencing on baseline primary cultured isolates to detect mixed infections and heteroresistance, respectively. Drug-sensitive, micro-heteroresistant, macro-heteroresistant, and fixed-resistant infections were defined as infections in which the frequency of resistance was <0.1%, 0.1%-4%, 5%-94%, and ≥95%, respectively, in resistance-conferring domains of the inhA promoter, the katG gene, and the rpoB gene.

Results

Of the 260 patients with tuberculosis included in the study, 25 (9.6%) had mixed infections and 30 (11.5%) had poor treatment outcomes. Micro-heteroresistance, macro-heteroresistance, and fixed resistance were found among 11 (4.2%), 2 (0.8%), and 11 (4.2%), respectively, for isoniazid and 21 (8.1%), 0 (0%), and 10 (3.8%), respectively, for rifampicin. In multivariable analysis, mixed infections but not heteroresistant infections independently predicted poor treatment outcomes.

Conclusions

Among patients starting first-line tuberculosis therapy in Botswana, mixed infections were associated with poor tuberculosis treatment outcomes, independent of heteroresistance.

Prevalence and Characterization of Heterogeneous Variable-Number Tandem-Repeat Clusters Comprising Drug-Susceptible and/or Variable Resistant Mycobacterium tuberculosis Complex Isolates in the Netherlands from 2004 to 2016

The variable-number tandem-repeat (VNTR) typing method is used to study tuberculosis (TB) transmission. Clustering of Mycobacterium tuberculosis isolates with identical VNTR patterns is assumed to reflect recent transmission. Hence, clusters are thought to be homogeneous regarding antibiotic resistance. In practice, however, heterogeneous clusters are also identified. This study investigates the prevalence and characteristics of heterogeneous VNTR clusters and assesses whether isolates in these clusters remain clustered when subjected to whole-genome sequencing (WGS). In the period from 2004 to 2016, 9,072 isolates were included. Demographic and epidemiological linkage data were obtained from the Netherlands Tuberculosis Register. VNTR clusters were defined as homogeneous when isolates shared identical resistance profiles or as heterogeneous if both susceptible and (variable) resistant isolates were found. Multivariate logistic regression analysis was performed to identify factors associated with heterogeneous clustering. Isolates from 2016 were subjected to WGS, and a genetic distance of 12 single nucleotide polymorphisms (SNPs) was used as the cutoff for WGS clustering. In total, 4,661/9,072 (51%) isolates were clustered into 985 different VNTR clusters, of which 217 (22%) were heterogeneous. Patient characteristics associated with heterogeneous clustering were non-Dutch ethnicity (odds ratio [OR], 1.46 [95% confidence interval {CI}, 1.22 to 1.75]), asylum seeker (OR, 1.51 [95% CI, 1.24 to 1.85]), extrapulmonary TB (OR, 1.26 [95% CI, 1.09 to 1.46]), previous TB diagnosis (OR, 1.38 [95% CI, 1.04 to 1.82]), and not being a contact of a TB patient (OR, 1.35 [95% CI, 1.08 to 1.69]). With WGS, 34% of heterogeneous and 78% of homogeneous isolates from 2016 remained clustered. Heterogeneous VNTR clusters are common but seem to be explained by a substantial degree of false clustering by VNTR typing compared to WGS.

Drug-resistant tuberculosis: challenges and opportunities for diagnosis and treatment

With an estimated incidence of 490000 cases in 2016, multidrug resistant tuberculosis (TB), against which key first-line anti-tuberculars are less efficacious, presents major challenges for global health. Poor treatment outcomes coupled with a yawning treatment gap between those in need of second-line therapy and those who receive it, underscore the urgent need for new approaches to tackle the scourge of drug-resistant TB. Against this background, significant progress has been made in understanding the complex biology of TB drug resistance and disease pathogenesis, and in establishing a pipeline for delivering new drugs and drug combinations. In this review, we highlight the challenges of drug-resistant TB and the ways in which new advances could be harnessed to improve treatment outcomes.

Drug resistance mechanisms and drug susceptibility testing for tuberculosis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (MTB) is the deadliest infectious disease and the associated global threat has worsened with the emergence of drug resistance, in particular multidrug-resistant TB (MDR-TB) and extensively drug-resistant TB (XDR-TB). Although the World Health Organization (WHO) End-TB Strategy advocates for universal access to antimicrobial susceptibility testing, this is not widely available and/or it is still underused. The majority of drug resistance in clinical MTB strains is attributed to chromosomal mutations. Resistance-related mutations could also exert certain fitness cost to the drug-resistant MTB strains and growth fitness could be restored by the presence of compensatory mutations. Understanding these underlying mechanisms could provide an important insight into TB pathogenesis and predict the future trend of MDR-TB global pandemic. This review covers the mechanisms of resistance in MTB and provides a comprehensive overview of current phenotypic and molecular approaches for drug susceptibility testing, with particular attention to the methods endorsed and recommended by the WHO.

Drug-resistant tuberculosis: An update on disease burden, diagnosis and treatment

The emergence of antimicrobial resistance against Mycobacterium tuberculosis, the leading cause of mortality due to a single microbial pathogen worldwide, represents a growing threat to public health and economic growth. The global burden of multidrug-resistant tuberculosis (MDR-TB) has recently increased by an annual rate of more than 20%. According to the World Health Organization approximately only half of all patients treated for MDR-TB achieved a successful outcome. For many years, patients with drug-resistant tuberculosis (TB) have received standardized treatment regimens, thereby accelerating the development of MDR-TB through drug-specific resistance amplification. Comprehensive drug susceptibility testing (phenotypic and/or genotypic) is necessary to inform physicians about the best drugs to treat individual patients with tailor-made treatment regimens. Phenotypic drug resistance can now often, but with variable sensitivity, be predicted by molecular drug susceptibility testing based on whole genome sequencing, which in the future could become an affordable method for the guidance of treatment decisions, especially in high-burden/resource-limited settings. More recently, MDR-TB treatment outcomes have dramatically improved with the use of bedaquiline-based regimens. Ongoing clinical trials with novel and repurposed drugs will potentially further improve cure-rates, and may substantially decrease the duration of MDR-TB treatment necessary to achieve relapse-free cure.

Mycobacterium tuberculosis Subculture Results in Loss of Potentially Clinically Relevant Heteroresistance

Multidrug-resistant tuberculosis (TB) presents a major public health dilemma. Heteroresistance, the coexistence of drug-resistant and drug-susceptible strains or of multiple drug-resistant strains with discrete haplotypes, may affect accurate diagnosis and the institution of effective treatment. Subculture, or passage of cells onto fresh growth medium, is utilized to preserve Mycobacterium tuberculosis cell lines and is universally employed in TB diagnostics. The impact of such passages, typically performed in the absence of drug, on drug-resistant subpopulations is hypothesized to vary according to the competitive costs of genotypic resistance-associated variants. We applied ultradeep next-generation sequencing to 61 phenotypically rifampin-monoresistant (n = 17) and preextensively (n = 41) and extensively (n = 3) drug-resistant isolates with presumptive heteroresistance at two time points in serial subculture. We found significant dynamic loss of minor-variant resistant subpopulations across all analyzed resistance-determining regions, including eight isolates (13%) whose antibiogram data would have transitioned from resistant to susceptible for at least one drug through subculture. Surprisingly, some resistance-associated variants appeared to be selected for in subculture.