Genetic Diversity in Mycobacterium tuberculosis Clinical Isolates and Resulting Outcomes of Tuberculosis Infection and Disease.

Tuberculosis claims more human lives than any other bacterial infectious disease and represents a clear and present danger to global health as new tools for vaccination, treatment, and interruption of transmission have been slow to emerge. Additionally, tuberculosis presents with notable clinical heterogeneity, which complicates diagnosis, treatment, and the establishment of nonrelapsing cure. How this heterogeneity is driven by the diversity ofclinical isolates of the causative agent, Mycobacterium tuberculosis, has recently garnered attention. Herein, we review advances in the understanding of how naturally occurring variation in clinical isolates affects transmissibility, pathogenesis, immune modulation, and drug resistance. We also summarize how specific changes in transcriptional responses can modulate infection or disease outcome, together with strain-specific effects on gene essentiality. Further understanding of how this diversity of M. tuberculosis isolates affects disease and treatment outcomes will enable the development of more effective therapeutic options and vaccines for this dreaded disease.

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.

Capacity building for whole genome sequencing of Mycobacterium tuberculosis and bioinformatics in high TB burden countries.

BACKGROUND: Whole genome sequencing (WGS) is increasingly used for Mycobacterium tuberculosis (Mtb) research. Countries with the highest tuberculosis (TB) burden face important challenges to integrate WGS into surveillance and research. METHODS: We assessed the global status of Mtb WGS and developed a 3-week training course coupled with long-term mentoring and WGS infrastructure building. Training focused on genome sequencing, bioinformatics and development of a locally relevant WGS research project. The aim of the long-term mentoring was to support trainees in project implementation and funding acquisition. The focus of WGS infrastructure building was on the DNA extraction process and bioinformatics. FINDINGS: Compared to their TB burden, Asia and Africa are grossly underrepresented in Mtb WGS research. Challenges faced resulted in adaptations to the training, mentoring and infrastructure building. Out-of-date laptop hardware and operating systems were overcome by using online tools and a Galaxy WGS analysis pipeline. A case studies approach created a safe atmosphere for students to formulate and defend opinions. Because quality DNA extraction is paramount for WGS, a biosafety level 3 and general laboratory skill training session were added, use of commercial DNA extraction kits was introduced and a 2-week training in a highly equipped laboratory was combined with a 1-week training in the local setting. INTERPRETATION: By developing and sharing the components of and experiences with a sequencing and bioinformatics training program, we hope to stimulate capacity building programs for Mtb WGS and empower high-burden countries to play an important role in WGS-based TB surveillance and research.

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.

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.

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.

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.

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.

Rifampicin-Monoresistant Tuberculosis Is Not the Same as Multidrug-Resistant Tuberculosis: a Descriptive Study from Khayelitsha, South Africa.

Rifampin monoresistance (RMR; rifampin resistance and isoniazid susceptibility) accounts for 38% of all rifampin-resistant tuberculosis (RR-TB) in South Africa and is increasing. We aimed to compare RMR-TB with multidrug-resistant TB (MDR-TB) in a setting with high TB, RR-TB, and HIV burdens. Patient-level clinical data and stored RR Mycobacterium tuberculosis isolates from 2008 to 2017 with available whole-genome sequencing (WGS) data were used to describe risk factors associated with RMR-TB and to compare RR-conferring mutations between RMR-TB and MDR-TB. A subset of isolates with particular RR-conferring mutations were subjected to semiquantitative rifampin phenotypic drug susceptibility testing. Among 2,041 routinely diagnosed RR-TB patients, 463 (22.7%) had RMR-TB. HIV-positive individuals (adjusted odds ratio [aOR], 1.4; 95% confidence interval [CI], 1.1 to 1.9) and diagnosis between 2013 and 2017 versus between 2008 and 2012 (aOR, 1.3; 95% CI, 1.1 to 1.7) were associated with RMR-TB. Among 1,119 (54.8%) patients with available WGS data showing RR-TB, significant differences in the distribution of rpoB RR-conferring mutations between RMR and MDR isolates were observed. Mutations associated with high-level RR were more commonly found among MDR isolates (811/889 [90.2%] versus 162/230 [70.4%] among RMR isolates; P < 0.0001). In particular, the rpoB L430P mutation, conferring low-level RR, was identified in 32/230 (13.9%) RMR isolates versus 10/889 (1.1%) in MDR isolates (P < 0.0001). Among 10 isolates with an rpoB L430P mutation, 7 were phenotypically susceptible using the critical concentration of 0.5 µg/ml (range, 0.125 to 1 µg/ml). The majority (215/230 [93.5%]) of RMR isolates showed susceptibility to all other TB drugs, highlighting the potential benefits of WGS for simplified treatment. These data suggest that the evolution of RMR-TB differs from MDR-TB with a potential contribution from HIV infection.

A landscape of genomic alterations at the root of a near-untreatable tuberculosis epidemic.

BACKGROUND: Atypical Beijing genotype Mycobacterium tuberculosis strains are widespread in South Africa and have acquired resistance to up to 13 drugs on multiple occasions. It is puzzling that these strains have retained fitness and transmissibility despite the potential fitness cost associated with drug resistance mutations. METHODS: We conducted Illumina sequencing of 211 Beijing genotype M. tuberculosis isolates to facilitate the detection of genomic features that may promote acquisition of drug resistance and restore fitness in highly resistant atypical Beijing forms. Phylogenetic and comparative genomic analysis was done to determine changes that are unique to the resistant strains that also transmit well. Minimum inhibitory concentration (MIC) determination for streptomycin and bedaquiline was done for a limited number of isolates to demonstrate a difference in MIC between isolates with and without certain variants. RESULTS: Phylogenetic analysis confirmed that two clades of atypical Beijing strains have independently developed resistance to virtually all the potent drugs included in standard (pre-bedaquiline) drug-resistant TB treatment regimens. We show that undetected drug resistance in a progenitor strain was likely instrumental in this resistance acquisition. In this cohort, ethionamide (ethA A381P) resistance would be missed in first-line drug-susceptible isolates, and streptomycin (gidB L79S) resistance may be missed due to an MIC close to the critical concentration. Subsequent inadequate treatment historically led to amplification of resistance and facilitated spread of the strains. Bedaquiline resistance was found in a small number of isolates, despite lack of exposure to the drug. The highly resistant clades also carry inhA promoter mutations, which arose after ethA and katG mutations. In these isolates, inhA promoter mutations do not alter drug resistance, suggesting a possible alternative role. CONCLUSION: The presence of the ethA mutation in otherwise susceptible isolates from ethionamide-naïve patients demonstrates that known exposure is not an adequate indicator of drug susceptibility. Similarly, it is demonstrated that bedaquiline resistance can occur without exposure to the drug. Inappropriate treatment regimens, due to missed resistance, leads to amplification of resistance, and transmission. We put these results into the context of current WHO treatment regimens, underscoring the risks of treatment without knowledge of the full drug resistance profile.

Proteogenomic Investigation of Strain Variation in Clinical Mycobacterium tuberculosis Isolates.

Mycobacterium tuberculosis consists of a large number of different strains that display unique virulence characteristics. Whole-genome sequencing has revealed substantial genetic diversity among clinical M. tuberculosis isolates, and elucidating the phenotypic variation encoded by this genetic diversity will be of the utmost importance to fully understand M. tuberculosis biology and pathogenicity. In this study, we integrated whole-genome sequencing and mass spectrometry (GeLC-MS/MS) to reveal strain-specific characteristics in the proteomes of two clinical M. tuberculosis Latin American-Mediterranean isolates. Using this approach, we identified 59 peptides containing single amino acid variants, which covered ∼9% of all coding nonsynonymous single nucleotide variants detected by whole-genome sequencing. Furthermore, we identified 29 distinct peptides that mapped to a hypothetical protein not present in the M. tuberculosis H37Rv reference proteome. Here, we provide evidence for the expression of this protein in the clinical M. tuberculosis SAWC3651 isolate. The strain-specific databases enabled confirmation of genomic differences (i.e., large genomic regions of difference and nonsynonymous single nucleotide variants) in these two clinical M. tuberculosis isolates and allowed strain differentiation at the proteome level. Our results contribute to the growing field of clinical microbial proteogenomics and can improve our understanding of phenotypic variation in clinical M. tuberculosis isolates.

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.

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.

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.

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.

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.

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

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.