Droplet based whole genome amplification for sequencing minute amounts of purified Mycobacterium tuberculosis DNA.

Implementation of whole genome sequencing (WGS) for patient care is hindered by limited Mycobacterium tuberculosis (Mtb) in clinical specimens and slow Mtb growth. We evaluated droplet multiple displacement amplification (dMDA) for amplification of minute amounts of Mtb DNA to enable WGS as an alternative to other Mtb enrichment methods. Purified genomic Mtb-DNA (0.1, 0.5, 1, and 5 pg) was encapsulated and amplified using the Samplix Xdrop-instrument and sequenced alongside a control sample using standard Illumina protocols followed by MAGMA-analysis. The control and 5 pg input dMDA samples underwent nanopore sequencing followed by Nanoseq and TB-profiler analysis. dMDA generated 105-2400 ng DNA from the 0.1-5 pg input DNA, respectively. Followed by Illumina WGS, dMDA raised mean sequencing depth from 7 × for 0.1 pg input DNA to ≥ 60 × for 5 pg input and the control sample. Bioinformatic analysis revealed a high number of false positive and false negative variants when amplifying ≤ 0.5 pg input DNA. Nanopore sequencing of the 5 pg dMDA sample presented excellent coverage depth, breadth, and accurate strain characterization, albeit elevated false positive and false negative variants compared to Illumina-sequenced dMDA sample with identical Mtb DNA input. dMDA coupled with Illumina WGS for samples with ≥ 5 pg purified Mtb DNA, equating to approximately 1000 copies of the Mtb genome, offers precision for drug resistance, phylogeny, and transmission insights.

Bacterial pathogens in Xpert MTB/RIF Ultra-negative sputum samples of patients with presumptive tuberculosis in a high TB burden setting: a 16S rRNA analysis.

In patients with presumptive tuberculosis (TB) in whom the diagnosis of TB was excluded, understanding the bacterial etiology of lower respiratory tract infections (LRTIs) is important for optimal patient management. A secondary analysis was performed on a cohort of 250 hospitalized patients with symptoms of TB. Bacterial DNA was extracted from sputum samples for Illumina 16S rRNA sequencing to identify bacterial species based on amplicon sequence variant level. The bacterial pathogen most likely to be responsible for the patients' LRTI could only be identified in a minority (6.0%, 13/215) of cases based on 16S rRNA amplicon sequencing: Mycoplasma pneumoniae (n = 7), Bordetella pertussis (n = 2), Acinetobacter baumanii (n = 2), and Pseudomonas aeruginosa (n = 2). Other putative pathogens were present in similar proportions of Xpert Ultra-positive and Xpert Ultra-negative sputum samples. The presence of Streptococcus (pseudo)pneumoniae appeared to increase the odds of radiological abnormalities (aOR 2.5, 95% CI 1.12-6.16) and the presence of S. (pseudo)pneumoniae (aOR 5.31, 95% CI 1.29-26.6) and Moraxella catarrhalis/nonliquefaciens (aOR 12.1, 95% CI 2.67-72.8) increased the odds of 6-month mortality, suggesting that these pathogens might have clinical relevance. M. pneumoniae, B. pertussis, and A. baumanii appeared to be the possible causes of TB-like symptoms. S. (pseudo)pneumoniae and M. catarrhalis/nonliquefaciens also appeared of clinical relevance based on 16S rRNA amplicon sequencing. Further research using tools with higher discriminatory power than 16S rRNA sequencing is required to develop optimal diagnostic and treatment strategies for this population.IMPORTANCEThe objective of this study was to identify possible bacterial lower respiratory tract infection (LRTI) pathogens in hospitalized patients who were initially suspected to have TB but later tested negative using the Xpert Ultra test. Although 16S rRNA was able to identify some less common or difficult-to-culture pathogens such as Mycoplasma pneumoniae and Bordetella pertussis, one of the main findings of the study is that, in contrast to what we had hypothesized, 16S rRNA is not a method that can be used to assist in the management of patients with presumptive TB having a negative Xpert Ultra test. Even though this could be considered a negative finding, we believe it is an important finding to report as it highlights the need for further research using different approaches.

Bedaquiline resistance in patients with drug-resistant tuberculosis in Cape Town, South Africa: a retrospective longitudinal cohort study.

BACKGROUND: Bedaquiline is a life-saving tuberculosis drug undergoing global scale-up. People at risk of weak tuberculosis drug regimens are a priority for novel drug access despite the potential source of Mycobacterium tuberculosis-resistant strains. We aimed to characterise bedaquiline resistance in individuals who had sustained culture positivity during bedaquiline-based treatment. METHODS: We did a retrospective longitudinal cohort study of adults (aged ≥18 years) with culture-positive pulmonary tuberculosis who received at least 4 months of a bedaquiline-containing regimen from 12 drug-resistant tuberculosis treatment facilities in Cape Town, South Africa, between Jan 20, 2016, and Nov 20, 2017. Sputum was programmatically collected at baseline (ie, before bedaquiline initiation) and each month to monitor treatment response per the national algorithm. The last available isolate from the sputum collected at or after 4 months of bedaquiline was designated the follow-up isolate. Phenotypic drug susceptibility testing for bedaquiline was done on baseline and follow-up isolates in MGIT960 media (WHO-recommended critical concentration of 1 µg/mL). Targeted deep sequencing for Rv0678, atpE, and pepQ, as well as whole-genome sequencing were also done. FINDINGS: In total, 40 (31%) of 129 patients from an estimated pool were eligible for this study. Overall, three (8%) of 38 patients assessable by phenotypic drug susceptibility testing for bedaquiline had primary resistance, 18 (47%) gained resistance (acquired or reinfection), and 17 (45%) were susceptible at both baseline and follow-up. Several Rv0678 and pepQ single-nucleotide polymorphisms and indels were associated with resistance. Although variants occurred in Rv0676c and Rv1979c, these variants were not associated with resistance. Targeted deep sequencing detected low-level variants undetected by whole-genome sequencing; however, none were in genes without variants already detected by whole-genome sequencing. Patients with baseline fluoroquinolone resistance, clofazimine exposure, and four or less effective drugs were more likely to have bedaquiline-resistant gain. Resistance gain was primarily due to acquisition; however, some reinfection by resistant strains occurred. INTERPRETATION: Bedaquiline-resistance gain, for which we identified risk factors, was common in these programmatically treated patients with sustained culture positivity. Our study highlights risks associated with implementing life-saving new drugs and shows evidence of bedaquiline-resistance transmission. Routine drug susceptibility testing should urgently accompany scale-up of new drugs; however, rapid drug susceptibility testing for bedaquiline remains challenging given the diversity of variants observed. FUNDING: Doris Duke Charitable Foundation, US National Institute of Allergy and Infectious Diseases, South African Medical Research Council, National Research Foundation, Research Foundation Flanders, Stellenbosch University Faculty of Medicine Health Sciences, South African National Research Foundation, Swiss National Science Foundation, and Wellcome Trust.

The MAGMA pipeline for comprehensive genomic analyses of clinical Mycobacterium tuberculosis samples.

BACKGROUND: Whole genome sequencing (WGS) holds great potential for the management and control of tuberculosis. Accurate analysis of samples with low mycobacterial burden, which are characterized by low (<20x) coverage and high (>40%) levels of contamination, is challenging. We created the MAGMA (Maximum Accessible Genome for Mtb Analysis) bioinformatics pipeline for analysis of clinical Mtb samples. METHODS AND RESULTS: High accuracy variant calling is achieved by using a long seedlength during read mapping to filter out contaminants, variant quality score recalibration with machine learning to identify genuine genomic variants, and joint variant calling for low Mtb coverage genomes. MAGMA automatically generates a standardized and comprehensive output of drug resistance information and resistance classification based on the WHO catalogue of Mtb mutations. MAGMA automatically generates phylogenetic trees with drug resistance annotations and trees that visualize the presence of clusters. Drug resistance and phylogeny outputs from sequencing data of 79 primary liquid cultures were compared between the MAGMA and MTBseq pipelines. The MTBseq pipeline reported only a proportion of the variants in candidate drug resistance genes that were reported by MAGMA. Notable differences were in structural variants, variants in highly conserved rrs and rrl genes, and variants in candidate resistance genes for bedaquiline, clofazmine, and delamanid. Phylogeny results were similar between pipelines but only MAGMA visualized clusters. CONCLUSION: The MAGMA pipeline could facilitate the integration of WGS into clinical care as it generates clinically relevant data on drug resistance and phylogeny in an automated, standardized, and reproducible manner.

Genomic transmission clusters and circulating lineages of Mycobacterium tuberculosis among refugees residing in refugee camps in Ethiopia.

BACKGROUND: Understanding the transmission dynamics of Mycobacterium tuberculosis (Mtb) could benefit the design of tuberculosis (TB) prevention and control strategies for refugee populations. Whole Genome Sequencing (WGS) has not yet been used to document the Mtb transmission dynamics among refugees in Ethiopia. We applied WGS to accurately identify transmission clusters and Mtb lineages among TB cases in refugee camps in Ethiopia. METHOD AND DESIGN: We conducted a cross-sectional study of 610 refugees in refugee camps in Ethiopia presenting with symptoms of TB. WGS data of 67 isolates was analyzed using the Maximum Accessible Genome for Mtb Analysis (MAGMA) pipeline; iTol and FigTree were used to visualize phylogenetic trees, lineages, and the presence of transmission clusters. RESULTS: Mtb culture-positive refugees originated from South Sudan (52/67, 77.6%), Somalia (9/67, 13.4%). Eritrea (4/67, 6%), and Sudan (2/67, 3%). The majority (52, 77.6%) of the isolates belonged to Mtb lineage (L) 3, and one L9 was identified from a Somalian refugee. The vast majority (82%) of the isolates were pan-susceptible Mtb, and none were multi-drug-resistant (MDR)-TB. Based on the 5-single nucleotide polymorphisms cutoff, we identified eight potential transmission clusters containing 23.9% of the isolates. Contact investigation confirmed epidemiological links with either family or social interaction within the refugee camps or with neighboring refugee camps. CONCLUSION: Four lineages (L1, L3, L4, and L9) were identified, with the majority of strains being L3, reflecting the Mtb L3 dominance in South Sudan, where the majority of refugees originated from. Recent transmission among refugees was relatively low (24%), likely due to the short study period. The improved understanding of the Mtb transmission dynamics using WGS in refugee camps could assist in designing effective TB control programs for refugees.

Correction: TBProfiler for automated calling of the association with drug resistance of variants in Mycobacterium tuberculosis.

[This corrects the article DOI: 10.1371/journal.pone.0279644.].

Integrating Molecular Diagnostics and GIS Mapping: A Multidisciplinary Approach to Understanding Tuberculosis Disease Dynamics in South Africa Using Xpert MTB/RIF.

An investigation was carried out to examine the use of national Xpert MTB/RIF data (2013-2017) and GIS technology for MTB/RIF surveillance in South Africa. The aim was to exhibit the potential of using molecular diagnostics for TB surveillance across the country. The variables analysed include Mycobacterium tuberculosis (Mtb) positivity, the mycobacterial proportion of rifampicin-resistant Mtb (RIF), and probe frequency. The summary statistics of these variables were generated and aggregated at the facility and municipal level. The spatial distribution patterns of the indicators across municipalities were determined using the Moran's I and Getis Ord (Gi) statistics. A case-control study was conducted to investigate factors associated with a high mycobacterial load. Logistic regression was used to analyse this study's results. There was striking spatial heterogeneity in the distribution of Mtb and RIF across South Africa. The median patient age, urban setting classification, and number of health care workers were found to be associated with the mycobacterial load. This study illustrates the potential of using data generated from molecular diagnostics in combination with GIS technology for Mtb surveillance in South Africa. Spatially targeted interventions can be implemented in areas where high-burden Mtb persists.

Estimating the effect of a rifampicin resistant tuberculosis diagnosis by the Xpert MTB/RIF assay on two-year mortality.

Studies assessing patient-centred outcomes of novel rifampicin resistant tuberculosis (RR-TB) diagnostics are rare and mostly apply conventional methods which may not adequately address biases. Even though the Xpert MTB/RIF molecular assay was endorsed a decade ago for simultaneous diagnosis of tuberculosis and RR-TB, the impact of the assay on mortality among people with RR-TB has not yet been assessed. We analysed data of an observational prospective cohort study (EXIT-RIF) performed in South Africa. We applied a causal inference approach using inverse odds of sampling weights to rectify survivor bias and selection bias caused by differing screening guidelines. We also adjusted for confounding using a marginal structural model with inverse probability of treatment weights. We estimated the total effect of an RR-TB diagnosis made by the Xpert assay versus the pre-Xpert diagnostic algorithm (entailing a targeted Line Probe Assay (LPA) among TB-confirmed patients) on two-year mortality and we assessed mediation by RR-treatment initiation. Of the 749 patients diagnosed with RR-TB [247 (33%) by the pre-Xpert diagnostic algorithm and 502 (67%) by the Xpert assay], 42.7% died. Of these, 364 (48.6%) patients died in the pre-Xpert group and 200 (39.8%) in the Xpert group. People diagnosed with RR-TB by the Xpert assay had a higher odds of RR-TB treatment initiation compared to those diagnosed by the targeted LPA-based diagnostic process (OR 2.79; 95%CI 2.19-3.56). Receiving an RR-TB diagnosis by Xpert resulted in a 28% reduction in the odds of mortality within 2 years after presentation to the clinic (ORCI 0.72; 95%CI 0.53-0.99). Causal mediation analysis suggests that the higher rate of RR-TB treatment initiation in people diagnosed by the Xpert assay explains the effect of Xpert on 2-year mortality [natural indirect effect odds ratio 0.90 (95%CI 0.85-0.96). By using causal inference methods in combination with high quality observational data, we could demonstrate that the introduction of the Xpert assay caused a 28% reduction in 2-year odds of mortality of RR-TB. This finding highlights the need for advocacy for a worldwide roll-out of rapid molecular tests. Because the effect is mainly caused by increased RR-TB treatment initiation, health care systems should also ensure timely initiation of effective treatment upon an RR-TB diagnosis.

Site-directed mutagenesis of Mycobacterium tuberculosis and functional validation to investigate potential bedaquiline resistance-causing mutations.

Molecular detection of bedaquiline resistant tuberculosis is challenging as only a small proportion of mutations in candidate bedaquiline resistance genes have been statistically associated with phenotypic resistance. We introduced two mutations, atpE Ile66Val and Rv0678 Thr33Ala, in the Mycobacterium tuberculosis H37Rv reference strain using homologous recombineering or recombination to investigate the phenotypic effect of these mutations. The genotype of the resulting strains was confirmed by Sanger- and whole genome sequencing, and bedaquiline susceptibility was assessed by minimal inhibitory concentration (MIC) assays. The impact of the mutations on protein stability and interactions was predicted using mutation Cutoff Scanning Matrix (mCSM) tools. The atpE Ile66Val mutation did not elevate the MIC above the critical concentration (MIC 0.25-0.5 µg/ml), while the MIC of the Rv0678 Thr33Ala mutant strains (> 1.0 µg/ml) classifies the strain as resistant, confirming clinical findings. In silico analyses confirmed that the atpE Ile66Val mutation minimally disrupts the bedaquiline-ATP synthase interaction, while the Rv0678 Thr33Ala mutation substantially affects the DNA binding affinity of the MmpR transcriptional repressor. Based on a combination of wet-lab and computational methods, our results suggest that the Rv0678 Thr33Ala mutation confers resistance to BDQ, while the atpE Ile66Val mutation does not, but definite proof can only be provided by complementation studies given the presence of secondary mutations.

A Bayesian approach to estimate the probability of resistance to bedaquiline in the presence of a genomic variant.

BACKGROUND: Bedaquiline is a core drug for treatment of rifampicin-resistant tuberculosis. Few genomic variants have been statistically associated with bedaquiline resistance. Alternative approaches for determining the genotypic-phenotypic association are needed to guide clinical care. METHODS: Using published phenotype data for variants in Rv0678, atpE, pepQ and Rv1979c genes in 756 Mycobacterium tuberculosis isolates and survey data of the opinion of 33 experts, we applied Bayesian methods to estimate the posterior probability of bedaquiline resistance and corresponding 95% credible intervals. RESULTS: Experts agreed on the role of Rv0678, and atpE, were uncertain about the role of pepQ and Rv1979c variants and overestimated the probability of bedaquiline resistance for most variant types, resulting in lower posterior probabilities compared to prior estimates. The posterior median probability of bedaquiline resistance was low for synonymous mutations in atpE (0.1%) and Rv0678 (3.3%), high for missense mutations in atpE (60.8%) and nonsense mutations in Rv0678 (55.1%), relatively low for missense (31.5%) mutations and frameshift (30.0%) in Rv0678 and low for missense mutations in pepQ (2.6%) and Rv1979c (2.9%), but 95% credible intervals were wide. CONCLUSIONS: Bayesian probability estimates of bedaquiline resistance given the presence of a specific mutation could be useful for clinical decision-making as it presents interpretable probabilities compared to standard odds ratios. For a newly emerging variant, the probability of resistance for the variant type and gene can still be used to guide clinical decision-making. Future studies should investigate the feasibility of using Bayesian probabilities for bedaquiline resistance in clinical practice.

Microfluidic Capture of Mycobacterium tuberculosis from Clinical Samples for Culture-Free Whole-Genome Sequencing.

Mycobacterium tuberculosis whole-genome sequencing (WGS) is a powerful tool as it can provide data on population diversity, drug resistance, disease transmission, and mixed infections. Successful WGS is still reliant on high concentrations of DNA obtained through M. tuberculosis culture. Microfluidics technology plays a valuable role in single-cell research but has not yet been assessed as a bacterial enrichment strategy for culture-free WGS of M. tuberculosis. In a proof-of-principle study, we evaluated the use of Capture-XT, a microfluidic lab-on-chip cleanup and pathogen concentration platform to enrich M. tuberculosis bacilli from clinical sputum specimens for downstream DNA extraction and WGS. Three of the four (75%) samples processed by the microfluidics application passed the library preparation quality control, compared to only one of the four (25%) samples not enriched by the microfluidics M. tuberculosis capture application. WGS data were of sufficient quality, with mapping depth of ≥25× and 9 to 27% of reads mapping to the reference genome. These results suggest that microfluidics-based M. tuberculosis cell capture might be a promising method for M. tuberculosis enrichment in clinical sputum samples, which could facilitate culture-free M. tuberculosis WGS. IMPORTANCE Diagnosis of tuberculosis is effective using molecular methods; however, a comprehensive characterization of the resistance profile of Mycobacterium tuberculosis often requires culturing and phenotypic drug susceptibility testing or culturing followed by whole-genome sequencing (WGS). The phenotypic route can take anywhere from 1 to >3 months to result, by which point the patient may have acquired additional drug resistance. The WGS route is a very attractive option; however, culturing is the rate-limiting step. In this original article, we provide proof-of-principle evidence that microfluidics-based cell capture can be used on high-bacillary-load clinical samples for culture-free WGS.

A machine-learning based model for automated recommendation of individualized treatment of rifampicin-resistant tuberculosis.

Background: Rifampicin resistant tuberculosis remains a global health problem with almost half a million new cases annually. In high-income countries patients empirically start a standardized treatment regimen, followed by an individualized regimen guided by drug susceptibility test (DST) results. In most settings, DST information is not available or is limited to isoniazid and fluoroquinolones. Whole genome sequencing could more accurately guide individualized treatment as the full drug resistance profile is obtained with a single test. Whole genome sequencing has not reached its full potential for patient care, in part due to the complexity of translating a resistance profile into the most effective individualized regimen. Methods: We developed a treatment recommender clinical decision support system (CDSS) and an accompanying web application for user-friendly recommendation of the optimal individualized treatment regimen to a clinician. Results: Following expert stakeholder meetings and literature review, nine drug features and 14 treatment regimen features were identified and quantified. Using machine learning, a model was developed to predict the optimal treatment regimen based on a training set of 3895 treatment regimen-expert feedback pairs. The acceptability of the treatment recommender CDSS was assessed as part of a clinical trial and in a routine care setting. Within the clinical trial setting, all patients received the CDSS recommended treatment. In 8 of 20 cases, the initial recommendation was recomputed because of stock out, clinical contra-indication or toxicity. In routine care setting, physicians rejected the treatment recommendation in 7 out of 15 cases because it deviated from the national TB treatment guidelines. A survey indicated that the treatment recommender CDSS is easy to use and useful in clinical practice but requires digital infrastructure support and training. Conclusions: Our findings suggest that global implementation of the novel treatment recommender CDSS holds the potential to improve treatment outcomes of rifampicin resistant tuberculosis.

TBProfiler for automated calling of the association with drug resistance of variants in Mycobacterium tuberculosis.

Following a huge global effort, the first World Health Organization (WHO)-endorsed catalogue of 17,356 variants in the Mycobacterium tuberculosis complex along with their classification as associated with resistance (interim), not associated with resistance (interim) or uncertain significance was made public In June 2021. This marks a critical step towards the application of next generation sequencing (NGS) data for clinical care. Unfortunately, the variant format used makes it difficult to look up variants when NGS data is generated by other bioinformatics pipelines. Furthermore, the large number of variants of uncertain significance in the catalogue hamper its useability in clinical practice. We successfully converted 98.3% of variants from the WHO catalogue format to the standardized HGVS format. We also created TBProfiler version 4.4.0 to automate the calling of all variants located in the tier 1 and 2 candidate resistance genes along with their classification when listed in the WHO catalogue. Using a representative sample of 339 clinical isolates from South Africa containing 691 variants in a tier 1 or 2 gene, TBProfiler classified 105 (15%) variants as conferring resistance, 72 (10%) as not conferring resistance and 514 (74%) as unclassified, with an average of 29 unclassified variants per isolate. Using a second cohort of 56 clinical isolates from a TB outbreak in Spain containing 21 variants in the tier 1 and 2 genes, TBProfiler classified 13 (61.9%) as unclassified, 7 (33.3%) as not conferring resistance, and a single variant (4.8%) classified as conferring resistance. Continued global efforts using standardized methods for genotyping, phenotyping and bioinformatic analyses will be essential to ensure that knowledge on genomic variants translates into improved patient care.

Bedaquiline resistance probability to guide treatment decision making for rifampicin-resistant tuberculosis: insights from a qualitative study.

BACKGROUND: Bedaquiline (BDQ) is a core drug for rifampicin-resistant tuberculosis (RR-TB) treatment. Accurate prediction of a BDQ-resistant phenotype from genomic data is not yet possible. A Bayesian method to predict BDQ resistance probability from next-generation sequencing data has been proposed as an alternative. METHODS: We performed a qualitative study to investigate the decision-making of physicians when facing different levels of BDQ resistance probability. Fourteen semi-structured interviews were conducted with physicians experienced in treating RR-TB, sampled purposefully from eight countries with varying income levels and burden of RR-TB. Five simulated patient scenarios were used as a trigger for discussion. Factors influencing the decision of physicians to prescribe BDQ at macro-, meso- and micro levels were explored using thematic analysis. RESULTS: The perception and interpretation of BDQ resistance probability values varied widely between physicians. The limited availability of other RR-TB drugs and the high cost of BDQ hindered physicians from altering the BDQ-containing regimen and incorporating BDQ resistance probability in their decision-making. The little experience with BDQ susceptibility testing and whole-genome sequencing results, and the discordance between phenotypic susceptibility and resistance probability were other barriers for physicians to interpret the resistance probability estimates. Especially for BDQ resistance probabilities between 25% and 70%, physicians interpreted the resistance probability value dynamically, and other factors such as clinical and bacteriological treatment response, history of exposure to BDQ, and resistance profile were often considered more important than the BDQ probability value for the decision to continue or stop BDQ. In this grey zone, some physicians opted to continue BDQ but added other drugs to strengthen the regimen. CONCLUSIONS: This study highlights the complexity of physicians' decision-making regarding the use of BDQ in RR-TB regimens for different levels of BDQ resistance probability.. Ensuring sufficient access to BDQ and companion drugs, improving knowledge of the genotype-phenotype association for BDQ resistance, availability of a rapid molecular test, building next-generation sequencing capacity, and developing a clinical decision support system incorporating BDQ resistance probability will all be essential to facilitate the implementation of BDQ resistance probability in personalizing treatment for patients with RR-TB.

Sequencing Mycobacteria and Algorithm-determined Resistant Tuberculosis Treatment (SMARTT): a study protocol for a phase IV pragmatic randomized controlled patient management strategy trial.

BACKGROUND: Rifampicin-resistant tuberculosis (RR-TB) remains an important global health problem. Ideally, the complete drug-resistance profile guides individualized treatment for all RR-TB patients, but this is only practised in high-income countries. Implementation of whole genome sequencing (WGS) technologies into routine care in low and middle-income countries has not become a reality due to the expected implementation challenges, including translating WGS results into individualized treatment regimen composition. METHODS: This trial is a pragmatic, single-blinded, randomized controlled medical device trial of a WGS-guided automated treatment recommendation strategy for individualized treatment of RR-TB. Subjects are 18 years or older and diagnosed with pulmonary RR-TB in four of the five health districts of the Free State province in South Africa. Participants are randomized in a 1:1 ratio to either the intervention (a WGS-guided automated treatment recommendation strategy for individualized treatment of RR-TB) or control (RR-TB treatment according to the national South African guidelines). The primary effectiveness outcome is the bacteriological response to treatment measured as the rate of change in time to liquid culture positivity during the first 6 months of treatment. Secondary effectiveness outcomes include cure rate, relapse rate (recurrence of RR-TB disease) and TB free survival rate in the first 12 months following RR-TB treatment completion. Additional secondary outcomes of interest include safety, the feasibility of province-wide implementation of the strategy into routine care, and health economic assessment from a patient and health systems perspective. DISCUSSION: This trial will provide important real-life evidence regarding the feasibility, safety, cost, and effectiveness of a WGS-guided automated treatment recommendation strategy for individualized treatment of RR-TB. Given the pragmatic nature, the trial will assist policymakers in the decision-making regarding the integration of next-generation sequencing technologies into routine RR-TB care in high TB burden settings. TRIAL REGISTRATION: ClinicalTrials.gov NCT05017324. Registered on August 23, 2021.

Variants in Bedaquiline-Candidate-Resistance Genes: Prevalence in Bedaquiline-Naive Patients, Effect on MIC, and Association with Mycobacterium tuberculosis Lineage.

Studies have shown that variants in bedaquiline-resistance genes can occur in isolates from bedaquiline-naive patients. We assessed the prevalence of variants in all bedaquiline-candidate-resistance genes in bedaquiline-naive patients, investigated the association between these variants and lineage, and the effect on phenotype. We used whole-genome sequencing to identify variants in bedaquiline-resistance genes in isolates from 509 bedaquiline treatment naive South African tuberculosis patients. A phylogenetic tree was constructed to investigate the association with the isolate lineage background. Bedaquiline MIC was determined using the UKMYC6 microtiter assay. Variants were identified in 502 of 509 isolates (98.6%), with the highest (85%) prevalence of variants in the Rv0676c (mmpL5) gene. We identified 36 unique variants, including 19 variants not reported previously. Only four isolates had a bedaquiline MIC equal to or above the epidemiological cut-off value of 0.25 µg/mL. Phylogenetic analysis showed that 14 of the 15 variants observed more than once occurred monophyletically in one Mycobacterium tuberculosis (sub)lineage. The bedaquiline MIC differed between isolates belonging to lineage 2 and 4 (Fisher's exact test, P = 0.0004). The prevalence of variants in bedaquiline-resistance genes in isolates from bedaquiline-naive patients is high, but very few (<2%) isolates were phenotypically resistant. We found an association between variants in bedaquiline resistance genes and Mycobacterium tuberculosis (sub)lineage, resulting in a lineage-dependent difference in bedaquiline phenotype. Future studies should investigate the impact of the presence of variants on bedaquiline-resistance acquisition and treatment outcome.

Transmission, distribution and drug resistance-conferring mutations of extensively drug-resistant tuberculosis in the Western Cape Province, South Africa.

Extensively drug-resistant tuberculosis (XDR-TB), defined as resistance to at least isoniazid (INH), rifampicin (RIF), a fluoroquinolone (FQ) and a second-line injectable drug (SLID), is difficult to treat and poses a major threat to TB control. The transmission dynamics and distribution of XDR Mycobacterium tuberculosis (Mtb) strains have not been thoroughly investigated. Using whole genome sequencing data on 461 XDR-Mtb strains, we aimed to investigate the geographical distribution of XDR-Mtb strains in the Western Cape Province of South Africa over a 10 year period (2006-2017) and assess the association between Mtb sub-lineage, age, gender, geographical patient location and membership or size of XDR-TB clusters. First, we identified transmission clusters by excluding drug resistance-conferring mutations and using the 5 SNP cutoff, followed by merging clusters based on their most recent common ancestor. We then consecutively included variants conferring resistance to INH, RIF, ethambutol (EMB), pyrazinamide (PZA), SLIDs and FQs in the cluster definition. Cluster sizes were classified as small (2-4 isolates), medium (5-20 isolates), large (21-100 isolates) or very large (>100 isolates) to reflect the success of individual strains. We found that most XDR-TB strains were clustered and that including variants conferring resistance to INH, RIF, EMB, PZA and SLIDs in the cluster definition did not significantly reduce the proportion of clustered isolates (85.5-82.2 %) but increased the number of patients belonging to small clusters (4.3-12.4 %, P=0.56). Inclusion of FQ resistance-conferring variants had the greatest effect, with 11 clustered isolates reclassified as unique while the number of clusters increased from 17 to 37. Lineage 2 strains (lineage 2.2.1 typical Beijing or lineage 2.2.2 atypical Beijing) showed the large clusters which were spread across all health districts of the Western Cape Province. We identified a significant association between residence in the Cape Town metropole and cluster membership (P=0.016) but no association between gender, age and cluster membership or cluster size (P=0.39). Our data suggest that the XDR-TB epidemic in South Africa probably has its origin in the endemic spread of MDR Mtb and pre-XDR Mtb strains followed by acquisition of FQ resistance, with more limited transmission of XDR Mtb strains. This only became apparent with the inclusion of drug resistance-conferring variants in the definition of a cluster. In addition to the prevention of amplification of resistance, rapid diagnosis of MDR, pre-XDR and XDR-TB and timely initiation of appropriate treatment is needed to reduce transmission of difficult-to-treat TB.

A treatment recommender clinical decision support system for personalized medicine: method development and proof-of-concept for drug resistant tuberculosis.

BACKGROUND: Personalized medicine tailors care based on the patient's or pathogen's genotypic and phenotypic characteristics. An automated Clinical Decision Support System (CDSS) could help translate the genotypic and phenotypic characteristics into optimal treatment and thus facilitate implementation of individualized treatment by less experienced physicians. METHODS: We developed a hybrid knowledge- and data-driven treatment recommender CDSS. Stakeholders and experts first define the knowledge base by identifying and quantifying drug and regimen features for the prototype model input. In an iterative manner, feedback from experts is harvested to generate model training datasets, machine learning methods are applied to identify complex relations and patterns in the data, and model performance is assessed by estimating the precision at one, mean reciprocal rank and mean average precision. Once the model performance no longer iteratively increases, a validation dataset is used to assess model overfitting. RESULTS: We applied the novel methodology to develop a treatment recommender CDSS for individualized treatment of drug resistant tuberculosis as a proof of concept. Using input from stakeholders and three rounds of expert feedback on a dataset of 355 patients with 129 unique drug resistance profiles, the model had a 95% precision at 1 indicating that the highest ranked treatment regimen was considered appropriate by the experts in 95% of cases. Use of a validation data set however suggested substantial model overfitting, with a reduction in precision at 1 to 78%. CONCLUSION: Our novel and flexible hybrid knowledge- and data-driven treatment recommender CDSS is a first step towards the automation of individualized treatment for personalized medicine. Further research should assess its value in fields other than drug resistant tuberculosis, develop solid statistical approaches to assess model performance, and evaluate their accuracy in real-life clinical settings.

Nanopore Sequencing for Mycobacterium tuberculosis: a Critical Review of the Literature, New Developments, and Future Opportunities.

The next-generation, short-read sequencing technologies that generate comprehensive, whole-genome data with single nucleotide resolution have already advanced tuberculosis diagnosis, treatment, surveillance, and source investigation. Their high costs, tedious and lengthy processes, and large equipment remain major hurdles for research use in high tuberculosis burden countries and implementation into routine care. The portable next-generation sequencing devices developed by Oxford Nanopore Technologies (ONT) are attractive alternatives due to their long-read sequence capability, compact low-cost hardware, and continued improvements in accuracy and throughput. A systematic review of the published literature demonstrated limited uptake of ONT sequencing in tuberculosis research and clinical care. Of the 12 eligible articles presenting ONT sequencing data on at least one Mycobacterium tuberculosis sample, four addressed software development for long-read ONT sequencing data with potential applications for M. tuberculosis. Only eight studies presented results of ONT sequencing of M. tuberculosis, of which five performed whole-genome and three did targeted sequencing. Based on these findings, we summarize the standard processes, reflect on the current limitations of ONT sequencing technology, and the research needed to overcome the main hurdles. The low capital cost, portable nature and continued improvement in the performance of ONT sequencing make it an attractive option for sequencing for research and clinical care, but limited data are available on its application in the tuberculosis field. Important research investment is needed to unleash the full potential of ONT sequencing for tuberculosis research and care.