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

The epidemiology, transmission, diagnosis, and management of drug-resistant tuberculosis-lessons from the South African experience

Drug-resistant tuberculosis (DR-TB) threatens to derail tuberculosis control efforts, particularly in Africa where the disease remains out of control. The dogma that DR-TB epidemics are fueled by unchecked rates of acquired resistance in inadequately treated or non-adherent individuals is no longer valid in most high DR-TB burden settings, where community transmission is now widespread. A large burden of DR-TB in Africa remains undiagnosed due to inadequate access to diagnostic tools that simultaneously detect tuberculosis and screen for resistance. Furthermore, acquisition of drug resistance to new and repurposed drugs, for which diagnostic solutions are not yet available, presents a major challenge for the implementation of novel, all-oral, shortened (6-9 months) treatment. Structural challenges including poverty, stigma, and social distress disrupt engagement in care, promote poor treatment outcomes, and reduce the quality of life for people with DR-TB. We reflect on the lessons learnt from the South African experience in implementing state-of-the-art advances in diagnostic solutions, deploying recent innovations in pharmacotherapeutic approaches for rapid cure, understanding local transmission dynamics and implementing interventions to curtail DR-TB transmission, and in mitigating the catastrophic socioeconomic costs of DR-TB. We also highlight globally relevant and locally responsive research priorities for achieving DR-TB control in South Africa.

Molecular determinants of multidrug-resistant tuberculosis in Sierra Leone

Multidrug-resistant tuberculosis (MDR-TB) management has become a serious global health challenge. Understanding its epidemic determinants on the regional level is crucial for developing effective control measures. We used whole genome sequencing data of 238 of Mycobacterium tuberculosis complex (MTBC) strains to determine drug resistance profiles, phylogeny, and transmission dynamics of MDR/rifampicin-resistant (RR) MTBC strains from Sierra Leone. Forty-two strains were classified as RR, 196 as MDR, 5 were resistant to bedaquiline (BDQ) and clofazimine (CFZ), but none was found to be resistant to fluoroquinolones. Sixty-one (26%) strains were resistant to all first-line drugs, three of which had additional resistance to BDQ/CFZ. The strains were classified into six major MTBC lineages (L), with strains of L4 being the most prevalent, 62% (n = 147), followed by L6 (Mycobacterium africanum) strains, (21%, n = 50). The overall clustering rate (using ≤d12 single-nucleotide polymorphism threshold) was 44%, stratified into 31 clusters ranging from 2 to 16 strains. The largest cluster (n = 16) was formed by sublineage 2.2.1 Beijing Ancestral 3 strains, which developed MDR several times. Meanwhile, 10 of the L6 strains had a primary MDR transmission. We observed a high diversity of drug resistance mutations, including borderline resistance mutations to isoniazid and rifampicin, and mutations were not detected by commercial assays. In conclusion, one in five strains investigated was resistant to all first-line drugs, three of which had evidence of BDQ/CFZ resistance. Implementation of interventions such as rapid diagnostics that prevent further resistance development and stop MDR-TB transmission chains in the country is urgently needed.

IMPORTANCE

A substantial proportion of MDR-TB strains in Sierra Leone were resistant against all first line drugs; however this makes the all-oral-six-month BPaLM regimen or other 6-9 months all oral regimens still viable, mainly because there was no FQ resistance.Resistance to BDQ was detected, as well as RR, due to mutations outside of the hotspot region. While the prevalence of those resistances was low, it is still cause for concern and needs to be closely monitored.

Clustering minimal inhibitory concentration data through Bayesian mixture models: An application to detect Mycobacterium tuberculosis resistance mutations

Antimicrobial resistance is becoming a major threat to public health throughout the world. Researchers are attempting to contrast it by developing both new antibiotics and patient-specific treatments. In the second case, whole-genome sequencing has had a huge impact in two ways: first, it is becoming cheaper and faster to perform whole-genome sequencing, and this makes it competitive with respect to standard phenotypic tests; second, it is possible to statistically associate the phenotypic patterns of resistance to specific mutations in the genome. Therefore, it is now possible to develop catalogues of genomic variants associated with resistance to specific antibiotics, in order to improve prediction of resistance and suggest treatments. It is essential to have robust methods for identifying mutations associated to resistance and continuously updating the available catalogues. This work proposes a general method to study minimal inhibitory concentration distributions and to identify clusters of strains showing different levels of resistance to antimicrobials. Once the clusters are identified and strains allocated to each of them, it is possible to perform regression method to identify with high statistical power the mutations associated with resistance. The method is applied to a new 96-well microtiter plate used for testing Mycobacterium tuberculosis.

Designing molecular diagnostics for current tuberculosis drug regimens

Diagnostic development must occur in parallel with drug development to ensure the longevity of new treatment compounds. Despite an increasing number of novel and repurposed anti-tuberculosis compounds and regimens, there remains a large number of drugs for which no rapid and accurate molecular diagnostic option exists. The lack of rapid drug susceptibility testing for linezolid, bedaquiline, clofazimine, the nitroimidazoles (i.e pretomanid and delamanid) and pyrazinamide at any level of the healthcare system compromises the effectiveness of current tuberculosis and drug-resistant tuberculosis treatment regimens. In the context of current WHO tuberculosis treatment guidelines as well as promising new regimens, we identify the key diagnostic gaps for initial and follow-on tests to diagnose emerging drug resistance and aid in regimen selection. Additionally, we comment on potential gene targets for inclusion in rapid molecular drug susceptibility assays and sequencing assays for novel and repurposed drug compounds currently prioritized in current regimens, and evaluate the feasibility of mutation detection given the design of existing technologies. Based on current knowledge, we also propose design priorities for next generation molecular assays to support triage of tuberculosis patients to appropriate and effective treatment regimens. We encourage assay developers to prioritize development of these key molecular assays and support the continued evolution, uptake, and utility of sequencing to build knowledge of tuberculosis resistance mechanisms and further inform rapid treatment decisions in order to curb resistance to critical drugs in current regimens and achieve End TB targets.Trial registration: ClinicalTrials.gov identifier: NCT05117788..

Comparison of in silico predicted Mycobacterium tuberculosis spoligotypes and lineages from whole genome sequencing data

Bacterial strain-types in the Mycobacterium tuberculosis complex underlie tuberculosis disease, and have been associated with drug resistance, transmissibility, virulence, and host-pathogen interactions. Spoligotyping was developed as a molecular genotyping technique used to determine strain-types, though recent advances in whole genome sequencing (WGS) technology have led to their characterization using SNP-based sub-lineage nomenclature. Notwithstanding, spoligotyping remains an important tool and there is a need to study the congruence between spoligotyping-based and SNP-based sub-lineage assignation. To achieve this, an in silico spoligotype prediction method ("Spolpred2") was developed and integrated into TB-Profiler. Lineage and spoligotype predictions were generated for > 28 k isolates and the overlap between strain-types was characterized. Major spoligotype families detected were Beijing (25.6%), T (18.6%), LAM (13.1%), CAS (9.4%), and EAI (8.3%), and these broadly followed known geographic distributions. Most spoligotypes were perfectly correlated with the main MTBC lineages (L1-L7, plus animal). Conversely, at lower levels of the sub-lineage system, the relationship breaks down, with only 65% of spoligotypes being perfectly associated with a sub-lineage at the second or subsequent levels of the hierarchy. Our work supports the use of spoligotyping (membrane or WGS-based) for low-resolution surveillance, and WGS or SNP-based systems for higher-resolution studies.

Insight into pathogenomics and phylogeography of hypervirulent and highly-lethal Mycobacterium tuberculosis strain cluster

BACKGROUND

. The Mycobacterium tuberculosis Beijing genotype is globally spread lineage with important medical properties that however vary among its subtypes. M. tuberculosis Beijing 14717-15-cluster was recently discovered as both multidrug-resistant, hypervirulent, and highly-lethal strain circulating in the Far Eastern region of Russia. Here, we aimed to analyze its pathogenomic features and phylogeographic pattern.

RESULTS

. The study collection included M. tuberculosis DNA collected between 1996 and 2020 in different world regions. The bacterial DNA was subjected to genotyping and whole genome sequencing followed by bioinformatics and phylogenetic analysis. The PCR-based assay to detect specific SNPs of the Beijing 14717-15-cluster was developed and used for its screening in the global collections. Phylogenomic and phylogeographic analysis confirmed endemic prevalence of the Beijing 14717-15-cluster in the Asian part of Russia, and distant common ancestor with isolates from Korea (> 115 SNPs). The Beijing 14717-15-cluster isolates had two common resistance mutations RpsL Lys88Arg and KatG Ser315Thr and belonged to spoligotype SIT269. The Russian isolates of this cluster were from the Asian Russia while 4 isolates were from the Netherlands and Spain. The cluster-specific SNPs that significantly affect the protein function were identified in silico in genes within different categories (lipid metabolism, regulatory proteins, intermediary metabolism and respiration, PE/PPE, cell wall and cell processes).

CONCLUSIONS

. We developed a simple method based on real-time PCR to detect clinically significant MDR and hypervirulent Beijing 14717-15-cluster. Most of the identified cluster-specific mutations were previously unreported and could potentially be associated with increased pathogenic properties of this hypervirulent M. tuberculosis strain. Further experimental study to assess the pathobiological role of these mutations is warranted.

Transcriptional regulation and drug resistance in Mycobacterium tuberculosis

Bacterial drug resistance is one of the major challenges to present and future human health, as the continuous selection of multidrug resistant bacteria poses at serious risk the possibility to treat infectious diseases in the near future. One of the infection at higher risk to become incurable is tuberculosis, due to the few drugs available in the market against Mycobacterium tuberculosis. Drug resistance in this species is usually due to point mutations in the drug target or in proteins required to activate prodrugs. However, another interesting and underexplored aspect of bacterial physiology with important impact on drug susceptibility is represented by the changes in transcriptional regulation following drug exposure. The main regulators involved in this phenomenon in M. tuberculosis are the sigma factors, and regulators belonging to the WhiB, GntR, XRE, Mar and TetR families. Better understanding the impact of these regulators in survival to drug treatment might contribute to identify new drug targets and/or to design new strategies of intervention.

Prediction of drug resistance by Sanger sequencing of Mycobacterium tuberculosis complex strains isolated from multidrug resistant tuberculosis suspect patients in Ethiopia

Background

Ethiopia is one of the high multidrug-resistant tuberculosis (MDR-TB) burden countries. However, phenotypic drug susceptibility testing can take several weeks due to the slow growth of Mycobacterium tuberculosis complex (MTBC) strains. In this study, we assessed the performance of a Sanger sequencing approach to predict resistance against five anti-tuberculosis drugs and the pattern of resistance mediating mutations.

Methods

We enrolled 226 MTBC culture-positive MDR-TB suspects and collected sputum specimens and socio-demographic and TB related data from each suspect between June 2015 and December 2016 in Addis Ababa, Ethiopia. Phenotypic drug susceptibility testing (pDST) for rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin using BACTEC MGIT 960 was compared with the results of a Sanger sequencing analysis of seven resistance determining regions in the genes rpoB, katG, fabG-inhA, pncA, embB, rpsL, and rrs.

Result

DNA isolation for Sanger sequencing was successfully extracted from 92.5% (209/226) of the MTBC positive cultures, and the remaining 7.5% (17/226) strains were excluded from the final analysis. Based on pDST results, drug resistance proportions were as follows: isoniazid: 109/209 (52.2%), streptomycin: 93/209 (44.5%), rifampicin: 88/209 (42.1%), ethambutol: 74/209 (35.4%), and pyrazinamide: 69/209 (33.0%). Resistance against isoniazid was mainly mediated by the mutation katG S315T (97/209, 46.4%) and resistance against rifampicin by rpoB S531L (58/209, 27.8%). The dominating resistance-conferring mutations for ethambutol, streptomycin, and pyrazinamide affected codon 306 in embB (48/209, 21.1%), codon 88 in rpsL (43/209, 20.6%), and codon 65 in pncA (19/209, 9.1%), respectively. We observed a high agreement between phenotypic and genotypic DST, such as 89.9% (at 95% confidence interval [CI], 84.2%–95.8%) for isoniazid, 95.5% (95% CI, 91.2%–99.8%) for rifampicin, 98.6% (95% CI, 95.9–100%) for ethambutol, 91.3% (95% CI, 84.6–98.1%) for pyrazinamide and 57.0% (95% CI, 46.9%–67.1%) for streptomycin.

Conclusion

We detected canonical mutations implicated in resistance to rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin. High agreement with phenotypic DST results for all drugs renders Sanger sequencing promising to be performed as a complementary measure to routine phenotypic DST in Ethiopia. Sanger sequencing directly from sputum may accelerate accurate clinical decision-making in the future.

Evolution and spread of a highly drug resistant strain of Mycobacterium tuberculosis in Papua New Guinea

BACKGROUND

Molecular mechanisms determining the transmission and prevalence of drug resistant tuberculosis (DR-TB) in Papua New Guinea (PNG) are poorly understood. We used genomic and drug susceptibility data to explore the evolutionary history, temporal acquisition of resistance and transmission dynamics of DR-TB across PNG.

METHODS

We performed whole genome sequencing on isolates from Central Public Health Laboratory, PNG, collected 2017-2019. Data analysis was done on a composite dataset that also included 100 genomes previously sequenced from Daru, PNG (2012-2015).

RESULTS

Sampled isolates represented 14 of the 22 PNG provinces, the majority (66/94; 70%) came from the National Capital District (NCD). In the composite dataset, 91% of strains were Beijing 2.2.1.1, identified in 13 provinces. Phylogenetic tree of Beijing strains revealed two clades, Daru dominant clade (A) and NCD dominant clade (B). Multi-drug resistance (MDR) was repeatedly and independently acquired, with the first MDR cases in both clades noted to have emerged in the early 1990s, while fluoroquinolone resistance emerged in 2009 (95% highest posterior density 2000-2016). We identified the presence of a frameshift mutation within Rv0678 (p.Asp47fs) which has been suggested to confer resistance to bedaquiline, despite no known exposure to the drug. Overall genomic clustering was significantly associated with rpoC compensatory and inhA promoter mutations (p < 0.001), with high percentage of most genomic clusters (12/14) identified in NCD, reflecting its role as a potential national amplifier.

CONCLUSIONS

The acquisition and evolution of drug resistance among the major clades of Beijing strain threaten the success of DR-TB treatment in PNG. With continued transmission of this strain in PNG, genotypic drug resistance surveillance using whole genome sequencing is essential for improved public health response to outbreaks. With occurrence of resistance to newer drugs such as bedaquiline, knowledge of full drug resistance profiles will be important for optimal treatment selection.

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.

Revised nomenclature and SNP barcode for Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis (Mtb) lineage 2 (L2) strains are present globally, contributing to a widespread tuberculosis (TB) burden, particularly in Asia where both prevalence of TB and numbers of drug resistant TB are highest. The increasing availability of whole-genome sequencing (WGS) data worldwide provides an opportunity to improve our understanding of the global genetic diversity of Mtb L2 and its association with the disease epidemiology and pathogenesis. However, existing L2 sublineage classification schemes leave >20 % of the Modern Beijing isolates unclassified. Here, we present a revised SNP-based classification scheme of L2 in a genomic framework based on phylogenetic analysis of >4000 L2 isolates from 34 countries in Asia, Eastern Europe, Oceania and Africa. Our scheme consists of over 30 genotypes, many of which have not been described before. In particular, we propose six main genotypes of Modern Beijing strains, denoted L2.2.M1–L2.2.M6. We also provide SNP markers for genotyping L2 strains from WGS data. This fine-scale genotyping scheme, which can classify >98 % of the studied isolates, serves as a basis for more effective monitoring and reporting of transmission and outbreaks, as well as improving genotype-phenotype associations such as disease severity and drug resistance. This article contains data hosted by Microreact.

Genetic variants and their association with phenotypic resistance to bedaquiline in Mycobacterium tuberculosis: a systematic review and individual isolate data analysis

Background

Methods

Findings

Interpretation

Funding

A Saudi Arabian Public Health Perspective of Tuberculosis

Tuberculosis is a global health challenge due to its spreading potential. The Kingdom of Saudi Arabia (KSA) faces a challenge in the spread of tuberculosis from migrant workers, but the foremost threat is the huge number of pilgrims who travel to visit sacred sites of the Islamic world located in the holy cities of Makkah and Al Madina. Pilgrims visit throughout the year but especially in the months of Ramadan and Zul-Hijah. The rise of resistance in Mycobacterium tuberculosis is an established global phenomenon that makes such large congregations likely hotspots in the dissemination and spread of disease at a global level. Although very stringent and effective measures exist, the threat remains due to the ever-changing dynamics of this highly pathogenic disease. This overview primarily highlights the current public health challenges posed by this disease to the Saudi health system, which needs to be highlighted not only to the concerned authorities of KSA, but also to the concerned global quarters since the pilgrims and migrants come from all parts of the world with a majority coming from high tuberculosis-burdened countries.

Molecular Analysis of Streptomycin Resistance Genes in Clinical Strains of Mycobacterium tuberculosis and Biocomputational Analysis of the MtGidB L101F Variant

Globally, tuberculosis (TB) remains a prevalent threat to public health. In 2019, TB affected 10 million people and caused 1.4 million deaths. The major challenge for controlling this infectious disease is the emergence and spread of drug-resistant Mycobacterium tuberculosis, the causative agent of TB. The antibiotic streptomycin is not a current first-line anti-TB drug. However, WHO recommends its use in patients infected with a streptomycin-sensitive strain. Several mutations in the M. tuberculosisrpsL, rrs and gidB genes have proved association with streptomycin resistance. In this study, we performed a molecular analysis of these genes in clinical isolates to determine the prevalence of known or novel mutations. Here, we describe the genetic analysis outcome. Furthermore, a biocomputational analysis of the MtGidB L101F variant, the product of a novel mutation detected in gidB during molecular analysis, is also reported as a theoretical approach to study the apparent genotype-phenotype association.

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

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

Drug resistance, fitness and compensatory mutations in Mycobacterium tuberculosis

For tuberculosis to be eradicated, the transmission of Multi-Drug-Resistant and eXtensively Drug Resistant strains of Mycobacterium tuberculosis (MDR and XDR-TB) must be considerably reduced. Drug resistant strains were initially thought to have reduced fitness, and the majority of resistant strains may actually have compromised fitness because they are found in only one or a few patients. In contrast, some MDR/XDR-TB strains are highly transmitted and cause large outbreaks. Most antibiotics target essential bacterial functions and the mutations that confer resistance to anti-TB drugs can incur fitness costs manifested as slower growth and reduced viability. The fitness costs vary with different resistance mutations and the bacilli can also accumulate secondary mutations that compensate for the compromised functions and partially or fully restore lost fitness. The compensatory mutations (CM) are different for each antibiotic, as they mitigate the deleterious effects of the specific functions compromised by the resistance mutations. CM are generally more common in strains with resistance mutations incurring the greatest fitness costs, but for RIF resistance, CM are most frequent in strains with the mutation carrying the least fitness cost, Ser450Leu. Here, we review what is known about fitness costs, CM and mechanisms of resistance to the drugs that define a strain as MDR or XDR-TB. The relative fitness costs of the resistance mutations and the mitigating effects of CM largely explain why certain mutations are frequently found in highly transmitted clusters while others are less frequently, rarely or never found in clinical isolates. The CM illustrate how drug resistance affects bacteria and how bacteria evolve to overcome the effects of the antibiotics, and thus a paradigm for how mycobacteria can evolve in response to stress.

Five-year microevolution of a multidrug-resistant Mycobacterium tuberculosis strain within a patient with inadequate compliance to treatment

BACKGROUND

Whole-genome sequencing has shown that the Mycobacterium tuberculosis infection process can be more heterogeneous than previously thought. Compartmentalized infections, exogenous reinfections, and microevolution are manifestations of this clonal complexity. The analysis of the mechanisms causing the microevolution -the genetic variability of M. tuberculosis at short time scales- of a parental strain into clonal variants with a patient is a relevant issue that has not been yet completely addressed. To our knowledge, a whole genome sequence microevolution analysis in a single patient with inadequate adherence to treatment has not been previously reported.

CASE PRESENTATION

In this work, we applied whole genome sequencing analysis for a more in-depth analysis of the microevolution of a parental Mycobacterium tuberculosis strain into clonal variants within a patient with poor treatment compliance in Argentina. We analyzed the whole-genome sequence of 8 consecutive Mycobacterium tuberculosis isolates obtained from a patient within 57-months of intermittent therapy. Nineteen mutations (9 short-term, 10 fixed variants) emerged, most of them associated with drug resistance. The first isolate was already resistant to isoniazid, rifampicin, and streptomycin, thereafter the strain developed resistance to fluoroquinolones and pyrazinamide. Surprisingly, isolates remained susceptible to the pro-drug ethionamide after acquiring a frameshift mutation in ethA, a gene required for its activation. We also found a novel variant, (T-54G), in the 5' untranslated region of whiB7 (T-54G), a region allegedly related to kanamycin resistance. Notably, discrepancies between canonical and phage-based susceptibility testing to kanamycin were previously found for the isolate harboring this mutation. In our patient, microevolution was mainly driven by drug selective pressure. Rare short-term mutations fixed together with resistance-conferring mutations during therapy.

CONCLUSIONS

This report highlights the relevance of whole-genome sequencing analysis in the clinic for characterization of pre-XDR and MDR resistance profile, particularly in patients with incomplete and/or intermittent treatment.

Population-level emergence of bedaquiline and clofazimine resistance-associated variants among patients with drug-resistant tuberculosis in southern Africa: a phenotypic and phylogenetic analysis

BACKGROUND

Bedaquiline and clofazimine are important drugs in the treatment of drug-resistant tuberculosis and are commonly used across southern Africa, although drug susceptibility testing is not routinely performed. In this study, we did a genotypic and phenotypic analysis of drug-resistant Mycobacterium tuberculosis isolates from cohort studies in hospitals in KwaZulu-Natal, South Africa, to identify resistance-associated variants (RAVs) and assess the extent of clofazimine and bedaquiline cross-resistance. We also used a comprehensive dataset of whole-genome sequences to investigate the phylogenetic and geographical distribution of bedaquiline and clofazimine RAVs in southern Africa.

METHODS

In this study, we included M tuberculosis isolates reported from the PRAXIS study of patients with drug-resistant tuberculosis treated with bedaquiline (King Dinuzulu Hospital, Durban) and three other cohort studies of drug-resistant tuberculosis in other KwaZulu-Natal hospitals, and sequential isolates from six persistently culture-positive patients with extensively drug-resistant tuberculosis at the KwaZulu-Natal provincial referral laboratory. Samples were collected between 2013 and 2019. Microbiological cultures were done as part of all parent studies. We sequenced whole genomes of included isolates and measured bedaquiline and clofazimine minimum inhibitory concentrations (MICs) for isolates identified as carrying any Rv0678 variant or previously published atpE, pepQ, and Rv1979c RAVs, which were the subject of the phenotypic study. We combined all whole-genome sequences of M tuberculosis obtained in this study with publicly available sequence data from other tuberculosis studies in southern Africa (defined as the countries of the Southern African Development Community), including isolates with Rv0678 variants identified by screening public genomic databases. We used this extended dataset to reconstruct phylogenetic relationships across lineage 2 and 4 M tuberculosis isolates.

FINDINGS

We sequenced the whole genome of 648 isolates from 385 patients with drug-resistant tuberculosis recruited into cohort studies in KwaZulu-Natal, and 28 isolates from six patients from the KwaZulu-Natal referral laboratory. We identified 30 isolates with Rv0678 RAVs from 16 (4%) of 391 patients. We did not identify any atpE, pepQ, or Rv1979c RAVs. MICs were measured for 21 isolates with Rv0678 RAVs. MICs were above the critical concentration for bedaquiline resistance in nine (43%) of 21 isolates, in the intermediate category in nine (43%) isolates, and within the wild-type range in three (14%) isolates. Clofazimine MICs in genetically wild-type isolates ranged from 0·12-0·5 μg/mL, and in isolates with RAVs from 0·25-4·0 μg/mL. Phylogenetic analysis of the extended dataset including M tuberculosis isolates from southern Africa resolved multiple emergences of Rv0678 variants in lineages 2 and 4, documented two likely nosocomial transmission events, and identified the spread of a possibly bedaquiline and clofazimine cross-resistant clone in eSwatini. We also identified four patients with pepQ frameshift mutations that may confer resistance.

INTERPRETATION

Bedaquiline and clofazimine cross-resistance in southern Africa is emerging repeatedly, with evidence of onward transmission largely due to Rv0678 mutations in M tuberculosis. Roll-out of bedaquiline and clofazimine treatment in the setting of limited drug susceptibility testing could allow further spread of resistance. Designing strong regimens would help reduce the emergence of resistance. Drug susceptibility testing is required to identify where resistance does emerge.

FUNDING

Wellcome Trust, National Institute of Allergy and Infectious Diseases and National Center for Advancing Translational Sciences of the National Institutes of Health.