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.

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 Detection of Mutations and Genotyping of Drug-Resistant Mycobacterium tuberculosis Strains Isolated from Patients in the Rural Eastern Cape Province.

Drug-resistant tuberculosis (DR-TB) is still a major public health concern in South Africa. Mutations in M. tuberculosis can cause varying levels of phenotypic resistance to anti-TB medications. There have been no prior studies on gene mutations and the genotyping of DR-TB in the rural Eastern Cape Province; hence, we aimed to identify DR-TB mutations, genetic diversity, and allocated lineages among patients in this area. Using Xpert® MTB/RIF, we assessed the rifampin resistance of sputum samples collected from 1157 patients suspected of having tuberculosis. GenoType MTBDR plus VER 2.0 was used for the detection of mutations causing resistance to anti-TB medications. The next step was to spoligotype 441 isolates. The most prevalent rifampin resistance-conferring mutations were in rpoB codon S531L in INH-resistant strains; the katG gene at codon S315TB and the inhA gene at codon C-15TB had the most mutations; 54.5% and 24.7%, respectively. In addition, 24.6% of strains showed mutations in both the rpoB and inhA genes, while 69.9% of strains showed mutations in both the katG and rpoB genes. Heteroresistance was seen in 17.9% of all cases in the study. According to spoligotyping analysis, Beijing families predominated. Investigation of the evolutionary lineages of M. tuberculosis isolates can be carried out using the information provided by the study's diversity of mutations. In locations wherein these mutations have been discovered, decision-making regarding the standardization of treatment regimens or individualized treatment may be aided by the detection frequency of rpoB, katG, and inhA mutations in various study areas.

The impact of genotype on the phenotype of Mycobacterium tuberculosis ΔsufR mutants.

Iron-sulphur (FeS) cluster biogenesis is a tightly regulated process in vivo. In Mycobacterium tuberculosis (Mtb), SufR functions as a transcriptional repressor of the operon encoding the primary FeS cluster biogenesis system. Previously, three independently isolated mutants (ΔRv1460stop_1.19, ΔRv1460stop _5.19 and ΔRv1460stop _5.20) harbouring the same deletion in sufR, displayed different growth kinetics in OADC supplemented 7H9 media. To investigate this discrepancy, we performed whole genome sequencing of the 3 mutants and the wild-type progenitor. Single nucleotide polymorphisms (SNPs) were identified in 3 genes in the ΔRv1460stop_1.19 mutant and one gene in the ΔRv1460stop_5.20 mutant. Phenotyping of the ΔRv1460stop_5.19 mutant, which had no additional SNPs, revealed increased susceptibility to clofazimine, DMNQ and menadione, while uptake and survival in THP-1 cells were not significantly different from the wild-type strain. Given that these results differ from those reported for other sufR deletion mutants (ΔSufRMTB and MtbΔSufR), they suggest that the position of the sufR deletion and the genotype of the progenitor strain impact the resulting phenotype.

Effect of compensatory evolution in the emergence and transmission of rifampicin-resistant Mycobacterium tuberculosis in Cape Town, South Africa: a genomic epidemiology study.

BACKGROUND: Experimental data show that drug-resistance-conferring mutations are often associated with a decrease in the replicative fitness of bacteria in vitro, and that this fitness cost can be mitigated by compensatory mutations; however, the role of compensatory evolution in clinical settings is less clear. We assessed whether compensatory evolution was associated with increased transmission of rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. METHODS: We did a genomic epidemiological study by analysing available M tuberculosis isolates and their associated clinical data from individuals routinely diagnosed with rifampicin-resistant tuberculosis in primary care and hospitals in Khayelitsha, Cape Town, South Africa. Isolates were collected as part of a previous study. All individuals diagnosed with rifampicin-resistant tuberculosis and with linked biobanked specimens were included in this study. We applied whole-genome sequencing, Bayesian reconstruction of transmission trees, and phylogenetic multivariable regression analysis to identify individual and bacterial factors associated with the transmission of rifampicin-resistant M tuberculosis strains. FINDINGS: Between Jan 1, 2008, and Dec 31, 2017, 2161 individuals were diagnosed with multidrug-resistant or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. Whole-genome sequences were available for 1168 (54%) unique individual M tuberculosis isolates. Compensatory evolution was associated with smear-positive pulmonary disease (adjusted odds ratio 1·49, 95% CI 1·08-2·06) and a higher number of drug-resistance-conferring mutations (incidence rate ratio 1·38, 95% CI 1·28-1·48). Compensatory evolution was also associated with increased transmission of rifampicin-resistant disease between individuals (adjusted odds ratio 1·55; 95% CI 1·13-2·12), independent of other patient and bacterial factors. INTERPRETATION: Our findings suggest that compensatory evolution enhances the in vivo fitness of drug-resistant M tuberculosis genotypes, both within and between patients, and that the in vitro replicative fitness of rifampicin-resistant M tuberculosis measured in the laboratory correlates with the bacterial fitness measured in clinical settings. These results emphasise the importance of enhancing surveillance and monitoring efforts to prevent the emergence of highly transmissible clones capable of rapidly accumulating new drug resistance mutations. This concern becomes especially crucial at present, because treatment regimens incorporating novel drugs are being implemented. FUNDING: Funding for this study was provided by a Swiss and South Africa joint research award (grant numbers 310030_188888, CRSII5_177163, and IZLSZ3_170834), the European Research Council (grant number 883582), and a Wellcome Trust fellowship (to HC; reference number 099818/Z/12/Z). ZS-D was funded through a PhD scholarship from the South African National Research Foundation and RMW was funded through the South African Medical Research Council.

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.

Discovery and biological evaluation of an adamantyl-amide derivative with likely MmpL3 inhibitory activity.

A series of molecules containing bulky lipophilic scaffolds was screened for activity against Mycobacterium tuberculosis and a number of compounds with antimycobacterial activity were identified. The most active compound, (2E)-N-(adamantan-1-yl)-3-phenylprop-2-enamide (C1), has a low micromolar minimum inhibitory concentration, low cytotoxicity (therapeutic index = 32.26), low mutation frequency and is active against intracellular Mycobacterium tuberculosis. Whole genome sequencing of mutants resistant to C1 showed a mutation in mmpL3 which may point to the involvement of MmpL3 in the antimycobacterial activity of the compound. In silico mutagenesis and molecular modelling studies were performed to better understand the binding of C1 within MmpL3 and the role that the specific mutation may play in the interaction at protein level. These analyses revealed that the mutation increases the energy required for binding of C1 within the protein translocation channel of MmpL3. The mutation also decreases the solvation energy of the protein, suggesting that the mutant protein might be more solvent-accessible, thereby restricting its interaction with other molecules. The results reported here describe a new molecule that may interact with the MmpL3 protein, providing insights into the effect of mutations on protein-ligand interactions and enhancing our understanding of this essential protein as a priority drug target.

Spatial Distribution of Drug-Resistant Mycobacterium tuberculosis Infections in Rural Eastern Cape Province of South Africa.

Tuberculosis (TB), an infectious airborne disease caused by Mycobacterium tuberculosis (Mtb), is a serious public health threat reported as the leading cause of morbidity and mortality worldwide. South Africa is a high-TB-burden country with TB being the highest infectious disease killer. This study investigated the distribution of Mtb mutations and spoligotypes in rural Eastern Cape Province. The Mtb isolates included were 1157 from DR-TB patients and analysed by LPA followed by spoligotyping of 441 isolates. The distribution of mutations and spoligotypes was done by spatial analysis. The rpoB gene had the highest number of mutations. The distribution of rpoB and katG mutations was more prevalent in four healthcare facilities, inhA mutations were more prevalent in three healthcare facilities, and heteroresistant isolates were more prevalent in five healthcare facilities. The Mtb was genetically diverse with Beijing more prevalent and largely distributed. Spatial analysis and mapping of gene mutations and spoligotypes revealed a better picture of distribution.

New insights on tuberculosis transmission dynamics and drug susceptibility profiles among the prison population in Southern Brazil based on whole-genome sequencing.

BACKGROUND: The rate of tuberculosis (TB) infection among the prison population (PP) in Brazil is 28 times higher than that in the general population, and prison environment favors the spread of TB. OBJECTIVE: To describe TB transmission dynamics and drug resistance profiles in PP using whole-genome sequencing (WGS). METHODS: This was a retrospective study of Mycobacterium tuberculosis cultivated from people incarcerated in 55 prisons between 2016 and 2019; only one isolate per prisoner was included. Information about movement from one prison to another was tracked. Clinical information was collected, and WGS was performed on isolates obtained at the time of TB diagnosis. RESULTS: Among 134 prisoners included in the study, we detected 16 clusters with a total of 58 (43%) cases of M. tuberculosis. Clusters ranged from two to seven isolates with five or fewer single nucleotide polymorphism (SNP) differences, suggesting a recent transmission. Six (4.4%) isolates were resistant to at least one anti-TB drug. Two of these clustered together and showed resistance to rifampicin, isoniazid, and fluoroquinolones, with 100% concordance between WGS and phenotypic drug-susceptibility testing. Prisoners with clustered isolates had a high amount of movement between prisons (two to eight moves) during the study period. CONCLUSIONS: WGS demonstrated the recent transmission of TB within prisons in Brazil. The high movement among prisoners seems to be related to the transmission of the same M. tuberculosis strain within the prison system. Screening for TB before and after the movement of prisoners using rapid molecular tests could play a role in reducing transmission.

Evaluation of Nanopore sequencing for Mycobacterium tuberculosis drug susceptibility testing and outbreak investigation: a genomic analysis.

BACKGROUND: Mycobacterium tuberculosis whole-genome sequencing (WGS) has been widely used for genotypic drug susceptibility testing (DST) and outbreak investigation. For both applications, Illumina technology is used by most public health laboratories; however, Nanopore technology developed by Oxford Nanopore Technologies has not been thoroughly evaluated. The aim of this study was to determine whether Nanopore sequencing data can provide equivalent information to Illumina for transmission clustering and genotypic DST for M tuberculosis. METHODS: In this genomic analysis, we analysed 151 M tuberculosis isolates from Madagascar, South Africa, and England, which were collected between 2011 and 2018, using phenotypic DST and matched Illumina and Nanopore data. Illumina sequencing was done with the MiSeq, HiSeq 2500, or NextSeq500 platforms and Nanopore sequencing was done on the MinION or GridION platforms. Using highly reliable PacBio sequencing assemblies and pairwise distance correlation between Nanopore and Illumina data, we optimise Nanopore variant filters for detecting single-nucleotide polymorphisms (SNPs; using BCFtools software). We then used those SNPs to compare transmission clusters identified by Nanopore with the currently used UK Health Security Agency Illumina pipeline (COMPASS). We compared Illumina and Nanopore WGS-based DST predictions using the Mykrobe software and mutation catalogue. FINDINGS: The Nanopore BCFtools pipeline identified SNPs with a median precision of 99·3% (IQR 99·1-99·6) and recall of 90·2% (88·1-94·2) compared with a precision of 99·6% (99·4-99·7) and recall of 91·9% (87·6-98·6) using the Illumina COMPASS pipeline. Using a threshold of 12 SNPs for putative transmission clusters, Illumina identified 98 isolates as unrelated and 53 as belonging to 19 distinct clusters (size range 2-7). Nanopore reproduced 15 out of 19 clusters perfectly; two clusters were merged into one cluster, one cluster had a single sample missing, and one cluster had an additional sample adjoined. Illumina-based clusters were also closely replicated using a five SNP threshold and clustering accuracy was maintained using mixed Illumina and Nanopore datasets. Genotyping resistance variants with Nanopore was highly concordant with Illumina, having zero discordant SNPs across more than 3000 SNPs and four insertions or deletions (indels), across 60 000 indels. INTERPRETATION: Illumina and Nanopore technologies can be used independently or together by public health laboratories performing M tuberculosis genotypic DST and outbreak investigations. As a result, clinical and public health institutions making decisions on which sequencing technology to adopt for tuberculosis can base the choice on cost (which varies by country), batching, and turnaround time. FUNDING: Academy for Medical Sciences, Oxford Wellcome Institutional Strategic Support Fund, and the Swiss South Africa Joint Research Award (Swiss National Science Foundation and South African National Research Foundation).

New insights on tuberculosis transmission dynamics and drug susceptibility profiles among the prison population in Southern Brazil based on whole-genome sequencing

ABSTRACT Background: The rate of tuberculosis (TB) infection among the prison population (PP) in Brazil is 28 times higher than that in the general population, and prison environment favors the spread of TB. Objective: To describe TB transmission dynamics and drug resistance profiles in PP using whole-genome sequencing (WGS). Methods: This was a retrospective study of Mycobacterium tuberculosis cultivated from people incarcerated in 55 prisons between 2016 and 2019; only one isolate per prisoner was included. Information about movement from one prison to another was tracked. Clinical information was collected, and WGS was performed on isolates obtained at the time of TB diagnosis. Results: Among 134 prisoners included in the study, we detected 16 clusters with a total of 58 (43%) cases of M. tuberculosis. Clusters ranged from two to seven isolates with five or fewer single nucleotide polymorphism (SNP) differences, suggesting a recent transmission. Six (4.4%) isolates were resistant to at least one anti-TB drug. Two of these clustered together and showed resistance to rifampicin, isoniazid, and fluoroquinolones, with 100% concordance between WGS and phenotypic drug-susceptibility testing. Prisoners with clustered isolates had a high amount of movement between prisons (two to eight moves) during the study period. Conclusions: WGS demonstrated the recent transmission of TB within prisons in Brazil. The high movement among prisoners seems to be related to the transmission of the same M. tuberculosis strain within the prison system. Screening for TB before and after the movement of prisoners using rapid molecular tests could play a role in reducing transmission.

Whole genome sequencing improves the discrimination between Mycobacterium bovis strains on the southern border of Kruger National Park, South Africa.

Background: Mycobacterium bovis forms part of the Mycobacterium tuberculosis complex and has an extensive host range and zoonotic potential. Various genotyping methods (e.g., spoligotyping) have been used to describe the molecular epidemiology of M. bovis. Advances in whole genome sequencing (WGS) have increased resolution to enable detection of genomic variants to the level of single nucleotide polymorphisms. This is especially relevant to One Health research on tuberculosis which benefits by being able to use WGS to identify epidemiologically linked cases, especially recent transmission. The use of WGS in molecular epidemiology has been extensively used in humans and cattle but is limited in wildlife. This approach appears to overcome the limitations of conventional genotyping methods due to lack of genetic diversity in M. bovis. Methods: This pilot study investigated the spoligotype and WGS of M. bovis isolates (n = 7) from wildlife in Marloth Park (MP) and compared these with WGS data from other South African M. bovis isolates. In addition, the greater resolution of WGS was used to explore the phylogenetic relatedness of M. bovis isolates in neighbouring wildlife populations. Results: The phylogenetic analyses showed the closest relatives to the seven isolates from MP were isolates from wildlife in Kruger National Park (KNP), which shares a border with MP. However, WGS data indicated that the KNP and MP isolates formed two distinct clades, even though they had similar spoligotypes and identical in silico genetic regions of difference profiles. Conclusions: Mycobacterium bovis isolates from MP were hypothesized to be directly linked to KNP wildlife, based on spoligotyping. However, WGS indicated more complex epidemiology. The presence of two distinct clades which were genetically distinct (SNP distance of 19-47) and suggested multiple transmission events. Therefore, WGS provided new insight into the molecular epidemiology of the M. bovis isolates from MP and their relationship to isolates from KNP. This approach will facilitate greater understanding of M. bovis transmission at wildlife-livestock-human interfaces and advances One Health research on tuberculosis, especially across different host species.

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.

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.

Ancient and recent differences in the intrinsic susceptibility of Mycobacterium tuberculosis complex to pretomanid.

OBJECTIVES: To develop a robust phenotypic antimicrobial susceptibility testing (AST) method with a correctly set breakpoint for pretomanid (Pa), the most recently approved anti-tuberculosis drug. METHODS: The Becton Dickinson Mycobacterial Growth Indicator Tube™ (MGIT) system was used at six laboratories to determine the MICs of a phylogenetically diverse collection of 356 Mycobacterium tuberculosis complex (MTBC) strains to establish the epidemiological cut-off value for pretomanid. MICs were correlated with WGS data to study the genetic basis of differences in the susceptibility to pretomanid. RESULTS: We observed ancient differences in the susceptibility to pretomanid among various members of MTBC. Most notably, lineage 1 of M. tuberculosis, which is estimated to account for 28% of tuberculosis cases globally, was less susceptible than lineages 2, 3, 4 and 7 of M. tuberculosis, resulting in a 99th percentile of 2 mg/L for lineage 1 compared with 0.5 mg/L for the remaining M. tuberculosis lineages. Moreover, we observed that higher MICs (≥8 mg/L), which probably confer resistance, had recently evolved independently in six different M. tuberculosis strains. Unlike the aforementioned ancient differences in susceptibility, these recent differences were likely caused by mutations in the known pretomanid resistance genes. CONCLUSIONS: In light of these findings, the provisional critical concentration of 1 mg/L for MGIT set by EMA must be re-evaluated. More broadly, these findings underline the importance of considering the global diversity of MTBC during clinical development of drugs and when defining breakpoints for AST.

Whole-Genome Sequencing Has the Potential To Improve Treatment for Rifampicin-Resistant Tuberculosis in High-Burden Settings: a Retrospective Cohort Study.

Treatment of multidrug-resistant or rifampicin-resistant tuberculosis (MDR/RR-TB), although improved in recent years with shorter, more tolerable regimens, remains largely standardized and based on limited drug susceptibility testing (DST). More individualized treatment with expanded DST access is likely to improve patient outcomes. To assess the potential of TB drug resistance prediction based on whole-genome sequencing (WGS) to provide more effective treatment regimens, we applied current South African treatment recommendations to a retrospective cohort of MDR/RR-TB patients from Khayelitsha, Cape Town. Routine DST and clinical data were used to retrospectively categorize patients into a recommended regimen, either a standardized short regimen or a longer individualized regimen. Potential regimen changes were then described with the addition of WGS-derived DST. WGS data were available for 1274 MDR/RR-TB patient treatment episodes across 2008 to 2017. Among 834 patients initially eligible for the shorter regimen, 385 (46%) may have benefited from reduced drug dosage or removing ineffective drugs when WGS data were considered. A further 187 (22%) patients may have benefited from more effective adjusted regimens. Among 440 patients initially eligible for a longer individualized regimen, 153 (35%) could have been switched to the short regimen. Overall, 305 (24%) patients had MDR/RR-TB with second-line TB drug resistance, where the availability of WGS-derived DST would have allowed more effective treatment individualization. These data suggest considerable benefits could accrue from routine access to WGS-derived resistance prediction. Advances in culture-free sequencing and expansion of the reference resistance mutation catalogue will increase the utility of WGS resistance prediction.

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.

Optimizing liquefaction and decontamination of sputum for DNA extraction from Mycobacterium tuberculosis.

Whole genome sequencing (WGS) can investigate the entire Mycobacterium tuberculosis (Mtb) genome but currently requires large amounts of mycobacterial DNA, necessitating culture. Culture-free Mtb WGS could revolutionize the clinical use of WGS but is hampered by the high viscosity, low mycobacterial load, and high contamination with bacterial and human DNA in sputum samples. To improve the sputum liquefaction and decontamination step prior to DNA extraction, we assessed the efficiency of Myco-TB, MycoPrep, and Sputolysin with/without TiKa-Kic in liquefying and decontaminating sputum and aimed to evaluate the effect of these approaches on mycobacterial viability, and Mtb DNA quality and quantity. Experiments using spiked sputum samples showed that Myco-TB and BD MycoPrep with standard (15 min) or increased (30 min) incubation time, but not reduced (7,5 min) incubation time performed well in liquefying and decontaminating sputum. No difference in DNA quality or quantity, contamination, or the amount of human DNA present was observed. In comparison, Sputolysin with/without TiKa-Kic was less effective for liquefaction and decontamination of sputum. PCR amplification of the human GAPDH gene after sputum treatment, showed the presence of human DNA in all samples, regardless of sputum treatment. Focused efforts are needed to deplete contaminating DNA for culture-free Mtb WGS.

Detection of minor variants in Mycobacterium tuberculosis whole genome sequencing data.

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