Mycobacterium tuberculosis Subculture Results in Loss of Potentially Clinically Relevant Heteroresistance

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.

Culture-Free Whole Genome Sequencing of Mycobacterium tuberculosis Using Ligand-Mediated Bead Enrichment Method

Background

Direct whole genome sequencing (WGS) of Mycobacterium tuberculosis (Mtb) can be used as a tool to study drug resistance, mixed infections, and within-host diversity. However, WGS is challenging to obtain from clinical samples due to low number of bacilli against a high background.

Methods

We prospectively collected 34 samples (sputum, n = 17; bronchoalveolar lavage, n = 13; and pus, n = 4) from patients with active tuberculosis (TB). Prior to DNA extraction, we used a ligand-mediated magnetic bead method to enrich Mtb from clinical samples and performed WGS on Illumina platform.

Results

Mtb was definitively identified based on WGS from 88.2% (30/34) of the samples, of which 35.3% (12/34) were smear negative. The overall median genome coverage was 15.2% (interquartile range [IQR], 7.7%-28.2%). There was a positive correlation between load of bacilli on smears and genome coverage (P < .001). We detected 58 genes listed in the World Health Organization mutation catalogue in each positive sample (median coverage, 85% [IQR, 61%-94%]), enabling the identification of mutations missed by routine diagnostics. Mutations causing resistance to rifampicin, isoniazid, streptomycin, and ethambutol were detected in 5 of 34 (14.7%) samples, including the rpoB S441A mutation that confers resistance to rifampicin, which is not covered by Xpert MTB/RIF.

Conclusions

We demonstrate the feasibility of magnetic bead-based enrichment for culture-free WGS of Mtb from clinical specimens, including smear-negative samples. This approach can also be integrated with low-cost sequencing workflows such as targeted sequencing for rapid detection of Mtb and drug resistance.

Exploring the potential of Oxford Nanopore Technologies sequencing for Mycobacterium tuberculosis sequencing: An assessment of R10 flowcells and V14 chemistry

Oxford Nanopore Technologies (ONT) sequencing is a promising technology. We assessed the performance of the new ONT R10 flowcells and V14 rapid sequencing chemistry for Mtb whole genome sequencing of Mycobacterium tuberculosis (Mtb) DNA extracted from clinical primary liquid cultures (CPLCs). Using the recommended protocols for MinION Mk1C, R10.4.1 MinION flowcells, and the ONT Rapid Sequencing Kit V14 on six CPLC samples, we obtained a pooled library yield of 10.9 ng/μl, generated 1.94 Gb of sequenced bases and 214k reads after 48h in a first sequencing run. Only half (49%) of all generated reads met the Phred Quality score threshold (>8). To assess if the low data output and sequence quality were due to impurities present in DNA extracted directly from CPLCs, we added a pre-library preparation bead-clean-up step and included purified DNA obtained from an Mtb subculture as a control sample in a second sequencing run. The library yield for DNA extracted from four CPLCs and one Mtb subculture (control) was similar (10.0 ng/μl), 2.38 Gb of bases and 822k reads were produced. The quality was slightly better with 66% of the produced reads having a Phred Quality >8. A third run of DNA from six CPLCs with bead clean-up pre-processing produced a low library yield (±1 Gb of bases, 166k reads) of low quality (51% of reads with a Phred Quality score >8). A median depth of coverage above 10× was only achieved for five of 17 (29%) sequenced libraries. Compared to Illumina WGS of the same samples, accurate lineage predictions and full drug resistance profiles from the generated ONT data could not be determined by TBProfiler. Further optimization of the V14 ONT rapid sequencing chemistry and library preparation protocol is needed for clinical Mtb WGS applications.

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%.

Co-resistance to isoniazid and second-line anti-tuberculosis drugs in isoniazid-resistant tuberculosis at a tertiary care hospital in Thailand

Isoniazid-resistant tuberculosis (Hr-TB) is an important drug-resistant tuberculosis (TB). In addition to rifampicin, resistance to other medications for Hr-TB can impact the course of treatment; however, there are currently limited data in the literature. In this study, the drug susceptibility profiles of Hr-TB treatment and resistance-conferring mutations were investigated for Hr-TB clinical isolates from Thailand. Phenotypic drug susceptibility testing (pDST) and genotypic drug susceptibility testing (gDST) were retrospectively and prospectively investigated using the Mycobacterium Growth Indicator Tube (MGIT), the broth microdilution (BMD) method, and whole-genome sequencing (WGS)-based gDST. The prevalence of Hr-TB cases was 11.2% among patients with TB. Most Hr-TB cases (89.5%) were newly diagnosed patients with TB. In the pDST analysis, approximately 55.6% (60/108) of the tested Hr-TB clinical isolates exhibited high-level isoniazid resistance. In addition, the Hr-TB clinical isolates presented co-resistance to ethambutol (3/161, 1.9%), levofloxacin (2/96, 2.1%), and pyrazinamide (24/118, 20.3%). In 56 Hr-TB clinical isolates, WGS-based gDST predicted resistance to isoniazid [katG S315T (48.2%) and fabG1 c-15t (26.8%)], rifampicin [rpoB L430P and rpoB L452P (5.4%)], and fluoroquinolones [gyrA D94G (1.8%)], but no mutation for ethambutol was detected. The categorical agreement for the detection of resistance to isoniazid, rifampicin, ethambutol, and levofloxacin between WGS-based gDST and the MGIT or the BMD method ranged from 80.4% to 98.2% or 82.1% to 100%, respectively. pDST and gDST demonstrated a low co-resistance rate between isoniazid and second-line TB drugs in Hr-TB clinical isolates.

IMPORTANCE

The prevalence of isoniazid-resistant tuberculosis (Hr-TB) is the highest among other types of drug-resistant tuberculosis. Currently, the World Health Organization (WHO) guidelines recommend the treatment of Hr-TB with rifampicin, ethambutol, pyrazinamide, and levofloxacin for 6 months. The susceptibility profiles of Hr-TB clinical isolates, especially when they are co-resistant to second-line drugs, are critical in the selection of the appropriate treatment regimen to prevent treatment failure. This study highlights the susceptibility profiles of the WHO-recommended treatment regimen in Hr-TB clinical isolates from a tertiary care hospital in Thailand and the concordance and importance of using the phenotypic drug susceptibility testing or genotypic drug susceptibility testing for accurate and comprehensive interpretation of results.

Genotype-Phenotype Characterization of Serial Mycobacterium tuberculosis Isolates in Bedaquiline-Resistant Tuberculosis

BACKGROUND

Emerging resistance to bedaquiline (BDQ) threatens to undermine advances in the treatment of drug-resistant tuberculosis (DRTB). Characterizing serial Mycobacterium tuberculosis (Mtb) isolates collected during BDQ-based treatment can provide insights into the etiologies of BDQ resistance in this important group of DRTB patients.

METHODS

We measured mycobacteria growth indicator tube (MGIT)-based BDQ minimum inhibitory concentrations (MICs) of Mtb isolates collected from 195 individuals with no prior BDQ exposure who were receiving BDQ-based treatment for DRTB. We conducted whole-genome sequencing on serial Mtb isolates from all participants who had any isolate with a BDQ MIC >1 collected before or after starting treatment (95 total Mtb isolates from 24 participants).

RESULTS

Sixteen of 24 participants had BDQ-resistant TB (MGIT MIC ≥4 µg/mL) and 8 had BDQ-intermediate infections (MGIT MIC = 2 µg/mL). Participants with pre-existing resistance outnumbered those with resistance acquired during treatment, and 8 of 24 participants had polyclonal infections. BDQ resistance was observed across multiple Mtb strain types and involved a diverse catalog of mmpR5 (Rv0678) mutations, but no mutations in atpE or pepQ. Nine pairs of participants shared genetically similar isolates separated by <5 single nucleotide polymorphisms, concerning for potential transmitted BDQ resistance.

CONCLUSIONS

BDQ-resistant TB can arise via multiple, overlapping processes, including transmission of strains with pre-existing resistance. Capturing the within-host diversity of these infections could potentially improve clinical diagnosis, population-level surveillance, and molecular diagnostic test development.

Rapid Mycobacterium abscessus antimicrobial susceptibility testing based on antibiotic treatment response mapping via Raman Microspectroscopy

OBJECTIVES

Antimicrobial susceptibility tests (ASTs) are pivotal tools for detecting and combating infections caused by multidrug-resistant rapidly growing mycobacteria (RGM) but are time-consuming and labor-intensive.

DESIGN

We used a Mycobacterium abscessus-based RGM model to develop a rapid (24-h) AST from the beginning of the strain culture, the Clinical Antimicrobials Susceptibility Test Ramanometry for RGM (CAST-R-RGM). The ASTs obtained for 21 clarithromycin (CLA)-treated and 18 linezolid (LZD)-treated RGM isolates.

RESULTS

CAST-R-RGM employs D2O-probed Raman microspectroscopy to monitor RGM metabolic activity, while also revealing bacterial antimicrobial drug resistance mechanisms. The results of clarithromycin (CLA)-treated and linezolid (LZD)-treated RGM isolates exhibited 90% and 83% categorical agreement, respectively, with conventional AST results of the same isolates. Furthermore, comparisons of time- and concentration-dependent Raman results between CLA- and LZD-treated RGM strains revealed distinct metabolic profiles after 48-h and 72-h drug treatments, despite similar profiles obtained for both drugs after 24-h treatments.

CONCLUSIONS

Ultimately, the rapid, accurate, and low-cost CAST-R-RGM assay offers advantages over conventional culture-based ASTs that warrant its use as a tool for improving patient treatment outcomes and revealing bacterial drug resistance mechanisms.

Storage of Mycobacterium tuberculosis culture isolates in MicrobankTM beads at a South African laboratory

Background

Mycobacterium tuberculosis complex (MTBC) isolates are typically stored at -70 °C in cryovials containing 1 mL aliquots of a liquid medium, with or without 50% glycerol. Multiple uses of the culture stock may decrease the strain viability while increasing the risk of culture contamination. Small culture aliquots may be more practical; however, storage capacity remains challenging. MicrobankTM beads (25 beads/vial) for the long-term storage of fungal cultures is well documented, but their use for storing MTBC isolates is uninvestigated.

Objective

The study aimed to determine the feasibility of using MicrobankTM beads for long-term storage of MTBC isolates at a laboratory in South Africa.

Methods

In February 2020, 20 isolates in liquid culture were stored in MicrobankTM beads, following an in-house developed protocol, at -70 °C. At defined time points (16 months [15 June 2021] and 21 months [18 November 2021]), two beads were retrieved from each storage vial and assessed for viability and level of contamination.

Results

Stored liquid isolates demonstrated MTBC growth within an average time-to-detection of 18 days following retrieval, even at 21 months post storage. Contaminating organisms were detected in 2 of 80 (2.5%) culture isolates.

Conclusion

MicrobankTM beads will allow for the reculture of up to 25 culture isolates using a reduced culture volume compared to current storage methods. MicrobankTM beads represent a storage solution for the medium-term storage of MTBC isolates.

What this study adds

This study evaluated the use of MicrobankTM beads as an alternate method for storing MTBC culture isolates at -70 °C and provided a suitable option for medium-term storage of MTBC.

High rifampicin peak plasma concentrations accelerate the slow phase of bacterial decline in tuberculosis patients: Evidence for heteroresistance

BACKGROUND

Antibiotic treatments are often associated with a late slowdown in bacterial killing. This separates the killing of bacteria into at least two distinct phases: a quick phase followed by a slower phase, the latter of which is linked to treatment success. Current mechanistic explanations for the in vitro slowdown are either antibiotic persistence or heteroresistance. Persistence is defined as the switching back and forth between susceptible and non-susceptible states, while heteroresistance is defined as the coexistence of bacteria with heterogeneous susceptibilities. Both are also thought to cause a slowdown in the decline of bacterial populations in patients and therefore complicate and prolong antibiotic treatments. Reduced bacterial death rates over time are also observed within tuberculosis patients, yet the mechanistic reasons for this are unknown and therefore the strategies to mitigate them are also unknown.

METHODS AND FINDINGS

We analyse a dose ranging trial for rifampicin in tuberculosis patients and show that there is a slowdown in the decline of bacteria. We show that the late phase of bacterial killing depends more on the peak drug concentrations than the total drug exposure. We compare these to pharmacokinetic-pharmacodynamic models of rifampicin heteroresistance and persistence. We find that the observation on the slow phase's dependence on pharmacokinetic measures, specifically peak concentrations are only compatible with models of heteroresistance and incompatible with models of persistence. The quantitative agreement between heteroresistance models and observations is very good ([Formula: see text]). To corroborate the importance of the slowdown, we validate our results by estimating the time to sputum culture conversion and compare the results to a different dose ranging trial.

CONCLUSIONS

Our findings indicate that higher doses, specifically higher peak concentrations may be used to optimize rifampicin treatments by accelerating bacterial killing in the slow phase. It adds to the growing body of literature supporting higher rifampicin doses for shortening tuberculosis treatments.

Strain diversity and gene mutations associated with presumptive multidrug-resistant Mycobacterium tuberculosis complex isolates in Northwest Ethiopia

OBJECTIVES

In this study, we assessed the genetic diversity and gene mutations that confer resistance to rifampicin (RIF), isoniazid (INH), fluoroquinolone (FQ), and second-line injectable (SLI) drugs in RIF-resistant (RR)/multidrug-resistant tuberculosis (MDR-TB) isolates in Northwest Ethiopia.

METHODS

Spoligotyping was used to assign isolates to TB lineages (Ls), and Hain line probe assays were used to detect resistance to RIF, INH, and FQs, and SLIs.

RESULTS

Among 130 analyzed strains, 68.5% were RR, and four major Mycobacterium tuberculosis complex lineages (L1, L3, L4, and L7) were identified with a predominance of the Euro-American L4 (72, 54.7%), while L7 genotypes were less common (3, 2.3%). Overall, the L4-T3-ETH (41, 32.0%), L3-CAS1-Delhi (29, 22.7%), and L3-CAS1-Killi (19, 14.8%) families were most common. Line probe analysis showed that among rpoB mutants, 65.2% were S450L, while 87.8% of katG mutants were S315T. Only three isolates showed mutation (c-15t) at the inhA gene, and no double mutation with katG and inhA genes was found. Six strains, two each of L1, L3, and L4, were resistant to FQs, having gyrA mutations (D94G, S91P), of which three isolates had additional resistance to SLI (rrs A1401G or C1402T mutations) including one isolate with low-level kanamycin (KAN) resistance.

CONCLUSIONS

This study showed a predominance of L4-T3-ETH, L3-CAS1-Delhi, and L3-CAS1-Killi families, with a high rate of rpoB_S450L and katG_S315T mutations and a low proportion of gyrA and rrs mutations. L7 was less frequently observed in this study. Further investigations are, therefore, needed to understand L7 and other lineages with undefined mutations.

Development and Assessment of a Novel Whole-Gene-Based Targeted Next-Generation Sequencing Assay for Detecting the Susceptibility of Mycobacterium tuberculosis to 14 Drugs

Targeted next-generation sequencing (tNGS) has emerged as an alternative method for detecting drug-resistant tuberculosis (DR-TB). To provide comprehensive drug susceptibility information and to address mutations missed by available commercial molecular diagnostics, we developed and evaluated a tNGS panel with 22 whole-gene targets using the Ion Torrent platform to predict drug resistance to 14 drugs, namely, rifampicin (RIF), isoniazid (INH), ethambutol (EMB), pyrazinamide (PZA), moxifloxacin (MFX), levofloxacin (LFX), amikacin (AMK), capreomycin (CM), kanamycin (KM), streptomycin (SM), bedaquiline (BDQ), clofazimine (CFZ), linezolid (LZD), and delamanid (DLM). We selected 50 and 35 Mycobacterium tuberculosis isolates with various DR profiles as the training set and the challenge set, respectively. Comparative variant analyses of the DR genes were performed using Sanger sequencing and whole-genome sequencing (WGS). Phenotypic drug susceptibility testing (pDST) results were used as gold standards. Regarding the limit of detection, the tNGS assay detected 2.9 to 3.8% minority variants in 4% mutant mixtures. The sensitivity and specificity of tNGS were 97.0% (95% confidence interval [CI] = 93.1 to 98.7%) and 99.1% (95% CI = 97.7 to 99.7%), respectively. The concordance of tNGS with pDST was 98.5% (95% CI = 97.2 to 99.2%), which was comparable to that of WGS (98.7%, 95% CI = 97.4 to 99.3%) and better than that of Sanger sequencing (96.9%, 95% CI = 95.3 to 98.0%). The agreement between tNGS and pDST was almost perfect for RIF, INH, EMB, MFX, LFX, AMK, CM, KM, SM, BDQ, and LZD (kappa value = 0.807 to 1.000) and substantial for PZA (kappa value = 0.791). Our customized novel whole-gene-based tNGS panel is highly consistent with pDST and WGS for comprehensive and accurate prediction of drug resistance in a strengthened and streamlined DR-TB laboratory program. IMPORTANCE We developed and validated a tNGS assay that was the first to target 22 whole genes instead of regions of drug resistance genes and comprehensively detected susceptibility to 14 anti-TB drugs, with great flexibility to include new or repurposed drugs. Notably, we demonstrated that our custom-designed Ion AmpliSeq TB research panel platform had high concordance with pDST and could significantly reduce turnaround time (by approximately 70%) to meet a clinically actionable time frame. Our tNGS assay is a promising DST solution for providing needed clinical information for precision medicine-guided therapies for DR-TB and allows the rollout of active pharmacovigilance.

Rapid Identification of Drug Resistance and Phylogeny in M. tuberculosis, Directly from Sputum Samples

Tuberculosis (TB) remains one of the most important infectious diseases globally. Establishing a resistance profile from the initial TB diagnosis is a priority. Rapid molecular tests evaluate only the most common genetic variants responsible for resistance to certain drugs, and Whole Genome Sequencing (WGS) needs culture prior to next-generation sequencing (NGS), limiting their clinical value. Targeted sequencing (TS) from clinical samples avoids these drawbacks, providing a signature of genetic markers that can be associated with drug resistance and phylogeny. In this study, a proof-of-concept protocol was developed for detecting genomic variants associated with drug resistance and for the phylogenetic classification of Mycobacterium Tuberculosis (Mtb) in sputum samples. Initially, a set of Mtb reference strains from the WHO were sequenced (WGS and TS). The results from the protocol agreed >95% with WHO reported data and phenotypic drug susceptibility testing (pDST). Lineage genetics results were 100% concordant with those derived from WGS. After that, the TS protocol was applied to sputum samples from TB patients to detect resistance to first- and second-line drugs and derive phylogeny. The accuracy was >90% for all evaluated drugs, except Eto/Pto (77.8%), and 100% were phylogenetically classified. The results indicate that the described protocol, which affords the complete drug resistance profile and phylogeny of Mtb from sputum, could be useful in the clinical area, advancing toward more personalized and more effective treatments in the near future. IMPORTANCE The COVID-19 pandemic negatively affected the progress in accessing essential Tuberculosis (TB) services and reducing the burden of TB disease, resulting in a decreased detection of new cases and increased deaths. Generating molecular diagnostic tests with faster results without losing reliability is considered a priority. Specifically, developing an antimicrobial resistance profile from the initial stages of TB diagnosis is essential to ensure appropriate treatment. Currently available rapid molecular tests evaluate only the most common genetic variants responsible for resistance to certain drugs, limiting their clinical value. In this work, targeted sequencing on sputum samples from TB patients was used to identify Mycobacterium tuberculosis mutations in genes associated with drug resistance and to derive a phylogeny of the infecting strain. This protocol constitutes a proof-of-concept toward the goal of helping clinicians select a timely and appropriate treatment by providing them with actionable information beyond current molecular approaches.

Evolution of Mycobacterium tuberculosis drug resistance in the genomic era

Mycobacterium tuberculosis has acquired drug resistance to all drugs that have been used against it, including those only recently introduced into clinical practice. Compared to other bacteria, it has a well conserved genome due to its role as an obligate human pathogen that has adapted to a niche over five to ten thousand years. These features facilitate reconstruction and dating of M. tuberculosis phylogenies, giving key insights into how resistance has been acquired and spread globally. Resistance to each new drug has occurred within five to ten years of clinical use and has occurred even more rapidly with recently introduced drugs. In most cases, resistance-conferring mutations come with a fitness cost, but this can be overcome by compensatory mutations which restore fitness to that of wild-type bacteria. It is likely that M. tuberculosis acquires drug resistance while maintaining limited genomic variability due the generation of low frequency within-host variation, combined with ongoing purifying selection causing loss of variants without a clear fitness advantage. However, variants that do confer an advantage, such as drug resistance, can increase in prevalence amongst all bacteria within a host and become the dominant clone. These resistant strains can then be transmitted leading to primary drug resistant infection in a new host. As many countries move towards genomic methods for diagnosis of M. tuberculosis infection and drug resistance, it is important to be aware of the implications for the evolution of resistance. Currently, understanding of resistance-conferring mutations is incomplete, and some targeted genetic diagnostics create their own selective pressures. We discuss an example where a rifampicin resistance-conferring mutation which was not routinely covered by standard testing became dominant. Finally, resistance to new drugs such as bedaquiline and delamanid is caused by individually rare mutations occurring across a large mutational genomic target that have been detected over a short time, and do not provide statistical power for genotype-phenotype correlation – in contrast to longer-established drugs that form the backbone of drug-sensitive antituberculosis therapy. Therefore, we need a different approach to identify resistance-conferring mutations of new drugs before their resistance becomes widespread, abrogating their usefulness.

Fluoroquinolone heteroresistance, antimicrobial tolerance, and lethality enhancement

With tuberculosis, the emergence of fluoroquinolone resistance erodes the ability of treatment to interrupt the progression of MDR-TB to XDR-TB. One way to reduce the emergence of resistance is to identify heteroresistant infections in which subpopulations of resistant mutants are likely to expand and make the infections fully resistant: treatment modification can be instituted to suppress mutant enrichment. Rapid DNA-based detection methods exploit the finding that fluoroquinolone-resistant substitutions occur largely in a few codons of DNA gyrase. A second approach for restricting the emergence of resistance involves understanding fluoroquinolone lethality through studies of antimicrobial tolerance, a condition in which bacteria fail to be killed even though their growth is blocked by lethal agents. Studies with Escherichia coli guide work with Mycobacterium tuberculosis. Lethal action, which is mechanistically distinct from blocking growth, is associated with a surge in respiration and reactive oxygen species (ROS). Mutations in carbohydrate metabolism that attenuate ROS accumulation create pan-tolerance to antimicrobials, disinfectants, and environmental stressors. These observations indicate the existence of a general death pathway with respect to stressors. M. tuberculosis displays a variation on the death pathway idea, as stress-induced ROS is generated by NADH-mediated reductive stress rather than by respiration. A third approach, which emerges from lethality studies, uses a small molecule, N-acetyl cysteine, to artificially increase respiration and additional ROS accumulation. That enhances moxifloxacin lethality with M. tuberculosis in culture, during infection of cultured macrophages, and with infection of mice. Addition of ROS stimulators to fluoroquinolone treatment of tuberculosis constitutes a new direction for suppressing the transition of MDR-TB to XDR-TB.

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.

Antibiotic heteroresistance in Mycobacterium tuberculosis isolates: a systematic review and meta-analysis

BACKGROUND

Mycobacterium tuberculosis (MTB) is responsible for tuberculosis; that continues to be a public health threat across the globe. Furthermore, increasing heteroresistance (HR)-the presence of resistant and susceptible isolates among MTB strains- has been reported from around the world. This phenomenon can lead to full resistance development and treatment failure.

METHODS

We systematically searched the relevant studies in PubMed, Scopus, and Embase (Until October 21, 2020). The study outcomes revealed the weighted pooled prevalence of antibiotic HR in MTB isolates with subgroup analysis by year, quality of study, and heteroresistance detection method.

RESULTS

A total of 38 studies which had investigated MTB isolates were included in the meta-analysis. Geographically, the highest number of studies were reported from Asia (n  =  24), followed by Africa (n  =  5). Nineteen studies reported HR to isoniazid, with a weighted pooled prevalence of 5% (95% CI 0-12) among 11,761 MTB isolates. Also, there is no important trend for the subgroup analysis by the study period (2001-2014 vs 2015-2017 vs 2018-2020). HR to rifampin was reported in 17 studies, with a weighted pooled prevalence of 7% (95% CI 2-14) among 3782 MTB isolates. HR to fluoroquinolone and ethambutol were reported in 12 and 4 studies, respectively, with weighted pooled prevalence of 10% and 1% among 2153 and 1509 MTB isolates, correspondingly.

CONCLUSION

Based on our analysis, HR in MTB isolates with different frequency rate is present worldwide. Thus, the selection of appropriate and reliable methods for HR detection is crucial for TB eradication.

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.

Mortality from drug-resistant tuberculosis in high-burden countries comparing routine drug susceptibility testing with whole-genome sequencing: a multicentre cohort study

BACKGROUND

Drug resistance threatens global tuberculosis control. We aimed to examine mortality in patients with tuberculosis from high-burden countries, according to concordance or discordance of results from drug susceptibility testing done locally and whole-genome sequencing (WGS).

METHODS

In this multicentre cohort study, we collected pulmonary Mycobacterium tuberculosis isolates and clinical data from individuals with tuberculosis from antiretroviral therapy programmes and tuberculosis clinics in Côte d'Ivoire, Democratic Republic of the Congo, Kenya, Nigeria, Peru, South Africa, and Thailand, stratified by HIV status and drug resistance. Sites tested drug susceptibility using routinely available methods. WGS was done on Illumina HiSeq 2500 in the USA and Switzerland, and TBprofiler was used to analyse the genomes. We included individuals aged 16 years or older with pulmonary tuberculosis (bacteriologically confirmed or clinically diagnosed). We analysed mortality in multivariable logistic regression models adjusted for sex, age, HIV status, history of tuberculosis, and sputum positivity.

FINDINGS

Between Sept 1, 2014, and July 4, 2016, of 634 patients included in our previous analysis, we included 582 patients with tuberculosis (median age 33 years [IQR 27-43], 225 [39%] women, and 247 [42%] HIV-positive). Based on WGS, 339 (58%) isolates were pan-susceptible, 35 (6%) monoresistant, 146 (25%) multidrug-resistant, and 24 (4%) pre-extensively drug-resistant (pre-XDR) or XDR. The analysis of mortality was based on 530 patients; 63 (12%) died and 77 (15%) patients received inappropriate treatment. Mortality ranged from 6% (18 of 310) in patients with pan-susceptible tuberculosis to 39% (nine of 23) in patients with pre-XDR or XDR tuberculosis. The adjusted odds ratio for mortality was 4·92 (95% CI 2·47-9·78) among undertreated patients, compared with appropriately treated patients.

INTERPRETATION

In seven countries with a high burden of tuberculosis, we observed discrepancies between drug resistance patterns obtained locally and WGS. The underdiagnosis of drug resistance resulted in inappropriate treatment and higher mortality. WGS can provide accurate and detailed drug resistance information required to improve the outcomes of drug-resistant tuberculosis in high-burden settings. Our results support WHO's call for point-of-care tests based on WGS.

FUNDING

National Institutes of Allergy and Infectious Diseases, Swiss National Science Foundation, and Swiss National Center for Mycobacteria.

Thin-Layer-Agar-Based Direct Phenotypic Drug Susceptibility Testing on Sputum in Eswatini Rapidly Detects Mycobacterium tuberculosis Growth and Rifampicin Resistance Otherwise Missed by WHO-Endorsed Diagnostic Tests

Xpert MTB/RIF rapidly detects resistance to rifampicin (RR); however, this test misses I491F-RR conferring rpoB mutation, common in southern Africa. In addition, Xpert MTB/RIF does not distinguish between viable and dead Mycobacterium tuberculosis (MTB). We aimed to investigate the ability of thin-layer agar (TLA) direct drug-susceptibility testing (DST) to detect MTB and its drug-resistance profiles in field conditions in Eswatini. Consecutive samples were tested in parallel with Xpert MTB/RIF and TLA for rifampicin (1.0 μg/ml) and ofloxacin (2.0 μg/ml). TLA results were compared at the Reference Laboratory in Antwerp with indirect-DST on Löwenstein-Jensen or 7H11 solid media and additional phenotypic and genotypic testing to resolve discordance. TLA showed a positivity rate for MTB detection of 7.1% versus 10.0% for Xpert MTB/RIF. Of a total of 4,547 samples included in the study, 200 isolates were available for comparison to the composite reference. Within a median of 18.4 days, TLA detected RR with 93.0% sensitivity (95% confidence interval [CI], 77.4 to 98.0) and 99.4% specificity (95% CI, 96.7 to 99.9) versus 62.5% (95% CI, 42.7 to 78.8) and 99.3% (95% CI, 96.2 to 99.9) for Xpert MTB/RIF. Eight isolates, 28.6% of all RR-confirmed isolates, carried the I491F mutation, all detected by TLA. TLA also correctly identified 183 of the 184 ofloxacin-susceptible isolates (99.5% specificity; 95% CI, 97.0 to 99.9). In field conditions, TLA rapidly detects RR, and in this specific setting, it contributed to detection of additional RR patients over Xpert MTB/RIF, mainly but not exclusively due to I491F. TLA also accurately excluded fluoroquinolone resistance.

Genetic diversity of the Mycobacterium tuberculosis complex strains from newly diagnosed tuberculosis patients in Northwest Ethiopia reveals a predominance of East-African-Indian and Euro-American lineages

OBJECTIVES

This study described the population structure of M. tuberculosis complex (MTBc) strains among patients with pulmonary or lymph node tuberculosis (TB) in Northwest Ethiopia and tested the performance of culture isolation and MPT64-based speciation for Lineage 7 (L7).

METHODS

Patients were recruited between April 2017 and June 2019 in North Gondar, Ethiopia. The MPT64 assay was used to confirm MTBc, and spoligotyping was used to characterize mycobacterial lineages. Line probe assay (LPA) was used to detect resistance to rifampicin and isoniazid.

RESULTS

Among 274 MTBc genotyped isolates, there were five MTBc lineages: L1-L4 and L7 were identified, with predominant East-African-Indian (L3) (53.6%) and Euro-American (L4) (40.1%) strains, and low prevalence (2.6%) of Ethiopia L7. The genotypes were similarly distributed between pulmonary and lymph node TB, and all lineages were equally isolated by culture and recognized as MTBc by the MPT64 assay. Additionally, LPA showed that 259 (94.5%) MTBc were susceptible to both rifampicin and isoniazid, and one (0.4%) was multi-drug resistant (resistant to both rifampicin and isoniazid).

CONCLUSION

These findings show that TB in North Gondar, Ethiopia, is mainly caused by L3 and L4 strains, with low rates of L7, confirmed as MTBc by MPT64 assay and with limited resistance to rifampicin and isoniazid.

Review of automated DNA extraction systems for sequencing-based solutions for drug-resistant tuberculosis detection

Robust clinical specimen nucleic acid extraction instrumentation and methods are critical to the performance of downstream molecular diagnostics for the diagnosis of drug-resistant tuberculosis (DR-TB). Currently, there is a high level of interest in sequencing-based solutions for rapid and comprehensive DR-TB testing from primary specimens (i.e., sputum). However, there is no standardized or fully automated sputum extraction system that has been widely implemented for use with Mycobacterium tuberculosis complex-containing sputum specimens. For sequencing-based technologies to be widely adopted in clinical laboratory settings in low- and middle-income countries, automated extraction technologies will be important to enhance scalability and reliability and to standardize performance of the downstream assays. Additionally, the ease of automatic technologies allows for faster uptake in laboratories currently without the expertise or infrastructure to perform manual extractions at the same automated throughput. This work is intended to provide an initial specification comparison of available automated DNA extraction systems that could serve as front-end components for existing and future sequencing approaches and provide the framework for future evaluations.

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.

Prediction of the hidden genotype of mixed infection strains in Iranian tuberculosis patients

BACKGROUND

Patients with mixed-strain Mycobacterium tuberculosis infections may be at a high risk of poor treatment outcomes. However, the mechanisms through which mixed infections affect the clinical manifestations are not well recognized. Evidence suggests that failure to detect the pathogen diversity within the host can influence the clinical results. We aimed to investigate the effects of different genotypes in mixed infections and determine their relationship with heteroresistance in the treatment of Iranian tuberculosis patients.

METHODS

One of the genotypes was identified in the culture and another genotype pattern in the mixed infection was predicted by comparing the pattern of MIRU-VNTR between the clinical specimens and their respective cultures in each patient. For all patients, the drug susceptibility testing was carried out on three single colonies from each clinical sample. The follow-up of patients was carried out during six months of treatment.

RESULTS

Based on MIRU-VNTR profiles of clinical samples, we showed that 55.6% (25/45) of the Iranian patients included in the study had mixed infections. Patients with mixed infections had a higher rate of treatment failure, compared to others (P=0.03). By comparing clinical sample profiles to profiles obtained after culture, we were able to distinguish between major and hidden genotypes. Among hidden genotypes, Haarlem (L4.1.2) and Beijing (L2) were associated to treatment failure (6/8 patients).

CONCLUSIONS

To conclude, we propose a procedure using the MIRU-VNTR method to identify the different genotypes in mixed infections. The present findings suggest that genotypes with potentially higher pathogenicity may not be detected when performing experimental culture in patients with mixed infections.

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.

Fluoroquinolone heteroresistance in Mycobacterium tuberculosis: detection by genotypic and phenotypic assays in experimentally mixed populations

Heteroresistance - the simultaneous presence of drug-susceptible and -resistant organisms - is common in Mycobacterium tuberculosis. In this study, we aimed to determine the limit of detection (LOD) of genotypic assays to detect gatifloxacin-resistant mutants in experimentally mixed populations. A fluoroquinolone-susceptible M. tuberculosis mother strain (S) and its in vitro selected resistant daughter strain harbouring the D94G mutation in gyrA (R) were mixed at different ratio’s. Minimum inhibitory concentrations (MICs) against gatifloxacin were determined, while PCR-based techniques included: line probe assays (Genotype MTBDRsl and GenoScholar-FQ + KM TB II), Sanger sequencing and targeted deep sequencing. Droplet digital PCR was used as molecular reference method. A breakpoint concentration of 0.25 mg/L allows the phenotypic detection of ≥1% resistant bacilli, whereas at 0.5 mg/L ≥ 5% resistant bacilli are detected. Line probe assays detected ≥5% mutants. Sanger sequencing required the presence of around 15% mutant bacilli to be detected as (hetero) resistant, while targeted deep sequencing detected ≤1% mutants. Deep sequencing and phenotypic testing are the most sensitive methods for detection of fluoroquinolone-resistant minority populations, followed by line probe assays (provided that the mutation is confirmed by a mutation band), while Sanger sequencing proved to be the least sensitive method.

Direct transmission of within-host Mycobacterium tuberculosis diversity to secondary cases can lead to variable between-host heterogeneity without de novo mutation: A genomic investigation

BACKGROUND

Whole genome sequencing (WGS) has enabled the development of new approaches to track Mycobacterium tuberculosis (Mtb) transmission between tuberculosis (TB) cases but its utility may be challenged by the discovery that Mtb diversifies within hosts. Nevertheless, there is limited data on the presence and degree of within-host evolution.

METHODS

We profiled a well-documented Mtb transmission cluster with three pulmonary TB cases to investigate within-host evolution and describe its impact on recent transmission estimates. We used deep sequencing to track minority allele frequencies (<50·0% abundance) during transmission and standard treatment.

FINDINGS

Pre-treatment (n = 3) and serial samples collected over 2 months of antibiotic treatment (n = 16) from all three cases were analysed. Consistent with the epidemiological data, zero fixed SNP separated all genomes. However, we identified six subclones between the three cases with an allele frequency ranging from 35·0% to 100·0% across sampling intervals. Five subclones were identified within the index case pre-treatment and shared with one secondary case, while only the dominant clone was observed in the other secondary case. By tracking the frequency of these heterogeneous alleles over the two-month therapy, we observed distinct signatures of drift and negative selection, but limited evidence for de novo mutations, even under drug pressure.

INTERPRETATION

We document within-host Mtb diversity in an index case, which led to transmission of minority alleles to a secondary case. Incorporating data on heterogeneous alleles may refine our understanding of Mtb transmission dynamics. However, more evidence is needed on the role of transmission bottleneck on observed heterogeneity between cases.

Effects of Host, Sample, and in vitro Culture on Genomic Diversity of Pathogenic Mycobacteria

Mycobacterium tuberculosis (M. tb), an obligate human pathogen and the etiological agent of tuberculosis (TB), remains a major threat to global public health. Comparative genomics has been invaluable for monitoring the emergence and spread of TB and for gaining insight into adaptation of M. tb. Most genomic studies of M. tb are based on single bacterial isolates that have been cultured for several weeks in vitro. However, in its natural human host, M. tb comprises complex, in some cases massive bacterial populations that diversify over the course of infection and cannot be wholly represented by a single genome. Recently, enrichment via hybridization capture has been used as a rapid diagnostic tool for TB, circumventing culturing protocols and enabling the recovery of M. tb genomes directly from sputum. This method has further applicability to the study of M. tb adaptation, as it enables a higher resolution and more direct analysis of M. tb genetic diversity within hosts with TB. Here we analyzed genomic material from M. tb and Mycobacterium bovis populations captured directly from sputum and from cultured samples using metagenomic and Pool-Seq approaches. We identified effects of sampling, patient, and sample type on bacterial genetic diversity. Bacterial genetic diversity was more variable and on average higher in sputum than in culture samples, suggesting that manipulation in the laboratory reshapes the bacterial population. Using outlier analyses, we identified candidate bacterial genetic loci mediating adaptation to these distinct environments. The study of M. tb in its natural human host is a powerful tool for illuminating host pathogen interactions and understanding the bacterial genetic underpinnings of virulence.

Feasibility of Direct Sputum Molecular Testing for Drug Resistance as Part of Tuberculosis Clinical Trials Eligibility Screening

A rapid diagnosis of drug-resistant tuberculosis (TB) is critical for early initiation of effective therapy. Molecular testing with line probe assays (MTBDRplus and MTBDRsl) on culture isolates has been available for some time and significantly reduces the time to diagnosis of drug resistance. However, routine use of this test directly on sputum is less common. As part of enrollment screening procedures for tuberculosis clinical trials conducted in Hanoi, Vietnam, we evaluated the feasibility and performance of line probe assay (LPA) testing directly on sputum samples from 315 participants with no prior history of TB treatment. Test performance characteristics for the detection of rifampin (RIF) and isoniazid (INH) drug resistance as compared to culture-based drug susceptibility testing (DST) reference standard were calculated. LPA demonstrated high sensitivity and specificity for the diagnosis of drug resistance. Scaling up molecular testing on sputum as part of time-sensitive clinical trial screening procedures in high TB burden settings is feasible and will reduce both time to initiation of appropriate therapy and the risk of late exclusions due to microbiologic ineligibility.

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.

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

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

Recent developments in the diagnosis of drug-resistant tuberculosis

Urgent steps are required to control the drug-resistant tuberculosis (TB) epidemic worldwide. Individualised treatment, using detailed drug-susceptibility test results to guide choice of antibiotics, improves patient outcomes and minimises adverse effects. Recent years have seen substantial advances in our ability to provide rapid, detailed drug-resistance profiles using genotypic methods for detection of mutations conferring drug-resistance. Rapid testing using real-time PCR to target the most important drug-resistance mutations allows the diagnosis of drug resistance to be made with the first diagnostic test, even in low resource settings. The use of whole genome sequencing to infer resistance to a range of different drugs facilitates earlier tailoring of therapy and detection of resistant subpopulations in mixed infections. Low burden countries, such as Australia are well positioned to lead the development and refinement of these new methods, to accelerate the incorporation of these new tools into TB control programs in high burden countries.

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.

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.

Cryptic Microheteroresistance Explains Mycobacterium tuberculosis Phenotypic Resistance

RATIONALE

Minority drug-resistant Mycobacterium tuberculosis subpopulations can be associated with phenotypic resistance but are poorly detected by Sanger sequencing or commercial molecular diagnostic assays.

OBJECTIVES

To determine the role of targeted next-generation sequencing in resolving these minor variant subpopulations.

METHODS

We used single molecule overlapping reads (SMOR), a targeted next-generation sequencing approach that dramatically reduces sequencing error, to analyze primary cultured isolates phenotypically resistant to rifampin, fluoroquinolones, or aminoglycosides, but for which Sanger sequencing found no resistance-associated variants (RAVs) within respective resistance-determining regions (study group). Isolates also underwent single-colony selection on antibiotic-containing agar, blinded to sequencing results. As a positive control, isolates with multiple colocalizing chromatogram peaks were also analyzed (control group).

MEASUREMENTS AND MAIN RESULTS

Among 61 primary culture isolates (25 study group and 36 control group), SMOR described 66 (49%) and 45 (33%) of 135 total heteroresistant RAVs at frequencies less than 5% and less than 1% of the total mycobacterial population, respectively. In the study group, SMOR detected minor resistant variant subpopulations in 80% (n = 20/25) of isolates with no Sanger-identified RAVs (median subpopulation size, 1.0%; interquartile range, 0.2-3.9%). Single-colony selection on drug-containing media corroborated SMOR results for 90% (n = 18/20) of RAV-containing specimens, and the absence of RAVs in 60% (n = 3/5) of isolates. Overall, Sanger sequencing was concordant with SMOR for 77% (n = 53/69) of macroheteroresistant (5-95% total population), but only 5% of microheteroresistant (<5%) subpopulations (n = 3/66) across both groups.

CONCLUSIONS

Cryptic minor variant mycobacterial subpopulations exist below the resolving capability of current drug susceptibility testing methodologies, and may explain an important proportion of false-negative resistance determinations.

Mycobacterial genomic DNA from used Xpert MTB/RIF cartridges can be utilised for accurate second-line genotypic drug susceptibility testing and spoligotyping

Xpert MTB/RIF (Xpert) is a widely-used test for tuberculosis (TB) and rifampicin-resistance. Second-line drug susceptibility testing (DST), which is recommended by policymakers, typically requires additional specimen collection that delays effective treatment initiation. We examined whether cartridge extract (CE) from used Xpert TB-positive cartridges was, without downstream DNA extraction or purification, suitable for both genotypic DST (MTBDRplus, MTBDRsl), which may permit patients to rapidly receive a XDR-TB diagnosis from a single specimen, and spoligotyping, which could facilitate routine genotyping. To determine the limit-of-detection and diagnostic accuracy, CEs from dilution series of drug-susceptible and -resistant bacilli were tested (MTBDRplus, MTBDRsl). Xpert TB-positive patient sputa CEs (n = 85) were tested (56 Xpert-rifampicin-susceptible, MTBDRplus and MTBDRsl; 29 Xpert-rifampicin-resistant, MTBDRsl). Spoligotyping was done on CEs from dilution series and patient sputa (n = 10). MTBDRplus had high non-valid result rates. MTBDRsl on CEs from dilutions ≥10³CFU/ml (C_T ≤ 24, >“low” Xpert semiquantitation category) was accurate, had low indeterminate rates and, on CE from sputa, highly concordant with MTBDRsl isolate results. CE spoligotyping results from dilutions ≥10³CFU/ml and sputa were correct. MTBDRsl and spoligotyping on CE are thus highly feasible. These findings reduce the need for additional specimen collection and culture, for which capacity is limited in high-burden countries, and have implications for diagnostic laboratories and TB molecular epidemiology.