Skip to main content

MDR M. tuberculosis outbreak clone in Eswatini missed by Xpert has elevated bedaquiline resistance dated to the pre-treatment era



Multidrug-resistant (MDR) Mycobacterium tuberculosis complex strains not detected by commercial molecular drug susceptibility testing (mDST) assays due to the RpoB I491F resistance mutation are threatening the control of MDR tuberculosis (MDR-TB) in Eswatini.


We investigate the evolution and spread of MDR strains in Eswatini with a focus on bedaquiline (BDQ) and clofazimine (CFZ) resistance using whole-genome sequencing in two collections ((1) national drug resistance survey, 2009–2010; (2) MDR strains from the Nhlangano region, 2014–2017).


MDR strains in collection 1 had a high cluster rate (95%, 117/123 MDR strains) with 55% grouped into the two largest clusters (gCL3, n = 28; gCL10, n = 40). All gCL10 isolates, which likely emerged around 1993 (95% highest posterior density 1987–1998), carried the mutation RpoB I491F that is missed by commercial mDST assays. In addition, 21 (53%) gCL10 isolates shared a Rv0678 M146T mutation that correlated with elevated minimum inhibitory concentrations (MICs) to BDQ and CFZ compared to wild type isolates. gCL10 isolates with the Rv0678 M146T mutation were also detected in collection 2.


The high clustering rate suggests that transmission has been driving the MDR-TB epidemic in Eswatini for three decades. The presence of MDR strains in Eswatini that are not detected by commercial mDST assays and have elevated MICs to BDQ and CFZ potentially jeopardizes the successful implementation of new MDR-TB treatment guidelines. Measures to limit the spread of these outbreak isolates need to be implemented urgently.


Tuberculosis (TB) is a significant global health problem that is exacerbated by the increasing incidence of drug-resistant TB (DR-TB) [1]. Delayed diagnosis and inappropriate treatment of DR-TB, together with suboptimal implementation of infection control measures, contribute to the emergence and transmission of multidrug-resistant and extensively drug-resistant (MDR, resistant to isoniazid [INH] and rifampicin [RMP]; XDR, MDR strains with additional resistance to a fluoroquinolone [FQ] and one second-line injectable drug [SLID]) Mycobacterium tuberculosis complex (Mtbc) strains [2, 3].

Eswatini, a small kingdom bordered by South Africa and Mozambique, has a high burden of TB, MDR-TB and HIV-TB, with an estimated TB incidence of 308 cases per 100,000 and an HIV prevalence of 27% among 15–45 year olds [2, 4]. Between 1997 and 2009, the burden of MDR-TB increased 7-fold from 0.9 to 7.7% among patients newly diagnosed with TB. Among retreatment cases, the proportions of MDR-TB quadrupled from 9.1% in 1997 to 33.9% in 2009 [4, 5]. To address the expanding MDR-TB epidemic, Eswatini has rolled-out the Xpert MTB/RIF, a rapid molecular TB diagnostic with the additional benefit of diagnosing RMP resistance; has expanded access to universal first- and second-line drug susceptibility testing (DST); and has decentralized DR-TB care [6].

A better understanding of the factors contributing and/or driving the MDR-TB epidemic in Eswatini is needed to inform these TB control strategies. This is particularly important because of two unusual findings in this setting. First, we previously showed that one in three MDR-Mtbc strains from the Eswatini 2009–2010 TB drug resistance survey (DRS) harbour the RpoB I491F mutation that is not interrogated by any of the molecular DST (mDST) assays endorsed by the World Health Organization (WHO) [4]. Thus, it is likely that a considerable number of patients with MDR-TB are incorrectly diagnosed as having drug-susceptible TB [4]. This mutation is estimated to account for 0.5% of RMP resistance globally but is frequent in Eswatini because of the clonal transmission of an outbreak strain that has also been described in neighbouring South Africa [7, 8]. By contrast, it is not clear to what extent clonal transmission is responsible for the proportion of RMP resistance that is detected by Xpert (i.e. strains with mutations in the RMP resistance determining region of rpoB).

Second, we have demonstrated that some of the RpoB I491F outbreak isolates harbour mutations in Rv0678, which encodes the repressor of the MmpS5-MmpL5 efflux pump [7]. If these mutations abolish or reduce the function of Rv0678, an increase in the minimum inhibitory concentrations (MICs) to both bedaquiline (BDQ) and clofazimine (CFZ) via increased efflux of both agents would be expected [9]. Both drugs have been prioritized for the longer MDR-TB regimens in the most recent WHO guidelines and are used in the standardized shorter MDR-TB regimen [10]. Moreover, CFZ is used throughout the entire duration of the recently WHO-endorsed shorter all-oral, BDQ-containing MDR-TB regimen [10, 11]. Clinically significant MIC increases to one or both drugs pre-dating the approval of BDQ would be a serious concern [7].

To investigate the resistance mechanisms, evolution, population structure, and transmission dynamics of MDR/XDR strains in Eswatini in detail, we first performed high-resolution genotyping, including whole-genome sequencing (WGS), for the isolates collected as part of the aforementioned national DRS from 2009 [4, 5]. Moreover, we conducted MIC testing for BDQ and CFZ to explore the potential effects of Rv0678 mutations on the resistance to both agents. Finally, we used WGS to assess the population structure of MDR strains and the frequency of Rv0678 mutations collected in the Nhlangano region between 2014 and 2017.


Strain collections

Collection 1

The isolates included in this analysis were collected as part of the TB-DRS in 2009–2010, which has been described in detail elsewhere [4]. Ethical approval was obtained from the Ministry of Health Scientific and Ethics Committee of Eswatini/Swaziland and the Ethics Review Board of Médecins Sans Frontières (MSF). Inclusion in the study was voluntary and after signing of an informed consent form. A total of 412 Mtbc strains were included in this study: 124 MDR, 1 XDR, 267 fully susceptible (DS), and 20 single-drug-resistant (SDR) or poly-drug-resistant (PDR) isolates but non-MDR. For general strain identification, genotyping by 24-loci mycobacterium interspersed repetitive unit-variable number of tandem repeat (MIRU-VNTR) typing and spoligotyping was performed. Additionally, WGS was undertaken for 273 isolates, including all DR isolates and a random selection of DS isolates (122 MDR [two MDR isolates could not be regrown to obtain sufficient DNA for WGS], 1 XDR, 20 SDR/PDR, and 130 DS isolates).

Four genomes from other sources (EMBL-EBI European Nucleotide Archive sequence read archive [study accession IDs: PRJNA393767 [12], PRJEB20942 [13], PRJEB5280 [14], and PRJNA395592 [15]]) were used to demonstrate that the Rv0678 N98D, G121R, and M146T mutations were homoplastic (Additional file 1: Table S1) [16, 17].

Collection 2

The genomes of 21 MDR isolates from a study evaluating thin-layer agar (TLA) DST conducted at the Nhlangano Health Centre microbiology laboratory were included to investigate cluster rates and genome characteristics of strains from a more recent period (i.e. between 2014 and 2017, with the majority of isolates from 2015 to 2016; Additional file 1: Table S2). The study protocol was approved by the Institutional Review Board of the Institute of Tropical Medicine (ITM); the Ethics Committee of the University Hospital of Antwerp, Belgium; and the Ministry of Health Scientific and Ethics Committee of Eswatini. Inclusion in the study was voluntary and after signing of an informed consent form. All consecutive patients investigated in Nhlangano (Shiselweni) for presumptive TB, older than 15 years, who had not received TB treatment in the previous month, and consented to be part of the study were included. Further details can be found in the supplementary methods (Additional file 2).

Laboratory procedures

Culture and drug susceptibility testing

Culture and phenotypic DST was performed as stated previously [4, 5]. CFZ and BDQ MICs were measured using the 1% proportion method with the MGIT960 system and the EpiCenter TBeXiST software according to the manufacturer’s instructions (Becton Dickinson, USA). The following concentrations were tested for both drugs: 0.125, 0.25, 0.5, 0.75, and 1 mg/L. Isolates susceptible at 0.125 or resistant at 1 mg/L tests were retested using an extended concentration range (i.e. 0.0312 and 0.0625, or 2 and 4 μg/mL, respectively). BDQ and CFZ MIC of Rv0678 mutant isolates (n = 25; 24 MDR and 1 DS), wild type isolates (n = 12; 7 MDR, 2 SDR/PDR, 4 DS; 4 of these were closely related to the Rv0678 mutant strains [3 gCL 10 strains, 1 S-type strain] and eight were randomly selected) were determined. H37Rv ATCC 27294 was included for quality control in each batch but excluded from the statistical analysis (see Additional file 2 and Additional file 1: Table S3).

Classical genotyping

Twenty-four-loci MIRU-VNTR typing and spoligotyping was done using standard approaches as described in supplemental methods previously (Additional file 2). Results were analysed using Bionumerics (version 7.6.3; Applied Maths [bioMérieux, Belgium]). Phylogenetic strain classification and MLVA-MTBC 15-9 nomenclature assignment was performed using the MIRU-VNTRplus database [18]. Clusters were defined as two or more Mtbc isolates sharing identical genotyping patterns according to both methods used.

Whole-genome sequencing and data analysis

Genomes were sequenced as described previously [19]. The analysis of WGS data was done with bioinformatic pipelines and parameters described before including those for evolutionary studies [19, 20]. A detailed description is provided in the supplement (Additional file 2). WGS data was submitted to the EMBL-EBI European Nucleotide Archive sequence read archive under the study IDs: PRJEB37777 [21], PRJEB6273 [22], PRJEB9680 [23], and PRJEB7281 [24] (Additional file 1: Tables S2 and S4).

Statistical analysis

Genotyping results were recorded at the molecular typing laboratory and added to the survey database.

Distributions of categorical variables between two comparison groups were compared with Fisher’s exact test. Comparisons of continuous variables were performed with a 2-sample t test or Wilcoxon rank-sum test. BDQ and CFZ MICs of mutant and wild type isolates were compared with the Mann-Whitney U test. We used an alpha level of 5% for all statistical tests.


Population structure, clustering rates, and evolution of Mycobacterium tuberculosis complex strains from national drug resistance survey from 2009 to 2010

Based on classical genotyping, the 412 Mtbc strains investigated were classified into nine previously defined phylogenetic lineages and sublineages, which showed no difference in distribution across the four study regions of HhoHho, Lubombo, Manzini, and Shiselweni (Table 1, Additional file 3: Fig. S1 and S2). Cluster (CL) analysis revealed that 278 of the 412 isolates investigated (67%) grouped into 60 MIRU-VNTR/spoligotype clusters (mCL) that comprised between 2 and 34 isolates (Additional file 1: Table S5). The mCL rate was significantly higher among MDR/XDR isolates (90%) compared to DS, SDR, or PDR isolates (58%) (p < 0.001). In fact, the two largest clusters mCL6 (n = 34, S-type) and mCL15 (n = 23, X-type) consisted of MDR isolates only (Additional file 1: Table S5).

Table 1 M. tuberculosis phylogenetic lineages identified in the isolates analysed

To further investigate the phylogeny and degree of clonality of MDR-Mtbc strains and to identify potential success markers on the genome level (e.g. potential compensatory mutations [19]), we performed WGS of a subset of 273 isolates (Additional file 3: Fig. S3). These included all drug-resistant isolates (122 MDR, 1 XDR, and 20 SDR/PDR) and a random selection of 130 DS isolates.

The genome-based phylogeny confirmed the strain classification based on classical genotyping (Fig. 1). WGS-based cluster analysis (gCL) based on a 12 SNP distance revealed that 173 of the 273 isolates were in 32 clusters ranging in size from 2 to 40 isolates (Fig. 1 and Additional file 1: Table S5). With 95%, the cluster rate among MDR/XDR was very high (117/123 analysed MDR/XDR isolates), whereas the rate of SDR/PDR/DS was only 37% (56/150 analysed SDR/PDR/DS isolates). Isolates of the two largest clusters (i.e. gCL3 [n = 28, X-type] and gCL10 [n = 40, S-type]) represented 55% of the M/XDR isolates investigated (Additional file 1: Table S5).

Fig. 1

Phylogenetic diversity of 273 M. tuberculosis complex isolates from Eswatini. Data are presented in a maximum likelihood tree (MLT) based on 12,062 SNP positions. Phylogenetic lineages are displayed. Cluster genome: isolates in one cluster share the same colour. Occurrence of resistance mediating and putative compensatory mutations in RpoB and RpoC are highlighted; identical mutations share the same colour

Bayesian coalescent analysis suggested that the most recent common ancestors of each of the 11 MDR genome clusters, respectively, emerged approximately 6.9 to 20.5 years prior to the Eswatini DRS in 2009–2010 (Table 2). The two largest MDR clusters, gCL10 and gCL3, likely emerged 17.2 years (12.4–22.6 years 95% HPD) and 20.5 years (13.6–24.5 95% HPD) prior to the DRS.

Table 2 Emergence of main MDR genome clusters and association with putative compensatory mutations

Resistance mediating and compensatory mutations

One hundred twenty-three isolates had variations in KatG, 129 in RpoB, 113 in EmbB, 120 in PncA, 67 in rrs (streptomycin [SM]), 8 in GyrA, and 65 had variations in RpsL (Fig. 1 and Additional file 1: Table S6). The most common mutations found in katG and rpoB were the low fitness cost mutations resulting in S315T amino acid substitution (84%, 117/139 INH-resistant isolates) and S450L (47%, 60/127 RMP-resistant isolates; S531L E. coli numbering), respectively. The second most frequent mutation in RMP-resistant isolates was RpoB I491F (31%, 40/127 RMP-resistant isolates; I572F E. coli numbering) (Additional file 1: Table S6). Twenty-five isolates were found to have mutations in Rv0678, a gene involved in resistance to BDQ and CFZ (one Beijing isolate, and 24 gCL10 MDR outbreak cluster isolates, see below and Additional file 1: Table S6).

The analysis of compensatory mutations revealed that 45 isolates had a variant in RpoB, 11 in RpoA, and 52 in RpoC (Additional file 1: Table S6). Among gCL and non-gCL RMP-resistant isolates, 73% (86/118) and 44% (4/9) harboured putative compensatory mutations, respectively.

We then correlated the detected resistance mutations with the gCL classification. This underlined the high clonality of the outbreak cluster isolates defined by the cluster analysis. For example, all 40 gCL10 isolates carried the same mutations in KatG (S315T), RpoB (I491F), EmbB (M306I), and PncA (H51D), including a subgroup of 24 isolates with additional mutations in Rv0678 (Additional file 1: Table S3). All 28 gCL3 isolates shared the same mutations in KatG (S315T), RpoB (S450L), EmbB (M306I), and PncA (R154G) (Additional file 1: Table S7).

Based on Bayesian coalescent analysis, we calculated a dated maximum clade credibility tree (MCCT) for the S-type lineage and mapped the identified drug resistance-associated mutations to the MCCT (Fig. 2). This analysis showed that isolates of gCL10 acquired resistance to all four first-line drugs and streptomycin (SM) before they started spreading in the community. Isolates of gCL10 acquired further variants in Rv0678 and compensatory mutations in RpoB and RpoC (Additional file 1: Table S6).

Fig. 2

Maximum clade credibility tree of the S-type lineage. gCL10, the largest MDR cluster in Eswatini, has been circulating for 17.2 years prior to the DRS conducted in 2009–2010. Three different Rv0678 mutations occurred within this cluster. The M146T mutation was the most common (n = 21) and likely arose 6.7 years prior to the DRS. The N98D (n = 2) and G121R (n = 1) must have evolved within 5.5 and 7.8 years prior to the DRS, respectively. Isolates, for which BDQ and CFZ MICs were measured, are highlighted (Fig. 3 and Additional file 1: Table S3)

Phylogenetic analysis of Rv0678 mutations and their impact on clofazimine and bedaquiline minimum inhibitory concentrations

We observed four different Rv0678 mutations in collection 1 (Additional file 1: Table S6). A single DS Beijing isolate harboured an A110V mutation, whereas the remaining three occurred in 24 of the 40 gCL10 MDR outbreak cluster isolates with the RpoB I491F mutation (21 isolates with Rv0678 M146T, 2 isolates with Rv0678 N98D, and 1 isolate with Rv0678 G121R; Fig. 2). The three mutations associated, i.e. Rv0678 G121R, N98D, and M146R, with lineage 4 (S-type) strains in the gCL10 MDR cluster evolved independently 5-8 years prior to the DRS (Additional file 1: Table S7). Further, all three mutations were found in other strain collections and importantly other Mtbc lineages such as lineage 4 (Haarlem sublineage) and lineage 2 (Beijing) (Additional file 1: Table S1). This homoplasy (i.e. independent acquisition of one trait/mutation in distinct phylogenetic backgrounds) is often a sign of positive selection [26].

To investigate their impact on resistance levels, we measured the BDQ and CFZ MICs for all Rv0678 mutants in comparison with genotypically wild type (gWT) control isolates that did not harbour any mutations in Rv0678, three of which were gCL10 strains without Rv0678 mutations (see the “Methods” section, Fig. 3A, B). The 21 gCL10 outbreak isolates with the Rv0678 M146T mutation had significantly elevated MICs for BDQ and CFZ relative to gWT strains (p < 0.001 for both drugs; Fig. 3A, B) as well as compared to gWT gCL10 strains (p < 0.05). This was also apparent when comparing the modes of the MIC distributions of the M146T mutation, which were three times higher than those of gWT isolates (i.e. 0.75 mg/L vs. 0.25 mg/L), although both distributions overlapped considerably. In fact, even if areas of technical uncertainty (ATUs, as outlined in Additional file 2) were introduced at 0.75–1 mg/L for both drugs, a proportion of these mutants would still be classified as susceptible because of the technical variation in MIC testing (i.e. the non-wild type cut-off value (NCOFF) for CFZ was 0.25 and 0.125 mg/L for BDQ, which would have to be tested to reliably detect this mutation, Fig. 3A–C).

Fig. 3

Effect of Rv0678 mutations on BDQ and CFZ MICs. The MICs for the DS Beijing mutant with the Rv0678 A110V mutation were compared with those of 24 gCL10-outbreak MDR isolates with either the Rv0678 N98D, G121R, or M146T mutations (Fig. 2). The NCOFFs for the Rv0678 M146T distributions for CFZ (part a) and BDQ (part b) were set using the eyeball method, and a corresponding idealized representation of the MIC distributions of the gWT and Rv0678 M146T populations for both drugs (and their area of overlap in grey) was prepared for illustrative purposes (part c) [25]. The MIC increases for M146T relative to gWT isolates were statistically significant for both CFZ and BDQ. The Rv0678 G121R and A110V mutations could not be evaluated statistically as the mutations only appeared once

The MIC increases for both drugs were also statistically significant for the two strains with the Rv0678 N98D mutation (BDQ: 100% of MICs ≥ 0.5 mg/L, p value 0.03297; CFZ: 100% of MICs ≥ 0.5 mg/L, p value 0.01099). The BDQ and CFZ MICs for the strain with Rv0678 mutation G121R were 2 mg/L and, thus, fell into the resistant range but could not be evaluated statistically as this mutation was observed only once. Similarly, statistical analyses could not be performed for the Rv0678 A110V mutation, which tested susceptible to CFZ and BDQ with MICs of 0.5 mg/L.

Population structure of multidrug-resistant Mycobacterium tuberculosis complex strains from the Nhlangano region between 2014 and 2017

The genome-based cluster analysis of MDR-Mtbc strains from the TLA study revealed a high cluster rate (20 out of 21) with isolates from four genome clusters already found in collection 1 (Additional file 1: Table S2). Isolates of the two dominant outbreaks now accounted for seven (gCL3, 33%) and nine (gCl10, 43%) of the 21 MDR strains investigated (because of the different sampling strategies, these frequencies could not be directly compared with those from the DRS). Six of the nine gCl10 isolates (67%) shared the aforementioned Rv0678 M146T mutation.


This study confirmed that clonal transmission of particular MDR-Mtbc outbreak clones has been the main driver of the MDR-TB epidemic in Eswatini. In fact, the two dominant gCL3 and gCL10 outbreak strains had already acquired resistance to all first-line drugs and SM when they began to circulate approximately 30 years ago. Our data show that in the DRS from 2009, the most frequent outbreak strain (gCL10) accounted for 30% of MDR isolates. Thereafter, it continued to spread, as shown by our analysis of collection 2, to reach 56% in the latest national DRS from 2017 to 2018 as outlined in the WHO Global tuberculosis report 2019 [1, 27]. Unfortunately, we did not have access to the genomes of the strains from this latest DRS, which would have provided an even clearer picture of the MDR-TB epidemiology in Eswatini. Nevertheless, our results supported a number of notable conclusions.

Although all WHO-endorsed mDST assays do not interrogate the RpoB I491F that is common to all gCL10 strains, we found that these strains all shared the KatG S315T mutation that is covered by all WHO-endorsed line probe assays for first-line drugs. Therefore, isolates with this KatG mutation could be prioritized for reflex testing with assays capable of detecting RpoB I491F to inform treatment decisions [28]. To this end, we have started to implement the targeted next-generation sequencing (NGS) assay Deeplex-MycTB in Eswatini, which not only covers all major resistance genes to traditional first- and second-line drugs as well as Rv0678 but also enables genotyping based on the spoligotype and phylogenetically informative mutations from primary specimens, e.g. sputum [29, 30]. Such a triage approach will, however, inevitably be slower than diagnosing RMP resistance testing with Xpert, which may provide a continued evolutionary advantage to gCL10 [6]. Accordingly, a prospective investigation of the current epidemiological scenario related to the spread of the gCL10 outbreak strains is planned.

In addition to being resistant to all first-line drugs, the gCL10 clone acquired multiple Rv0678 mutations that correlate with elevated MICs to BDQ and/or CFZ and have been circulating since 2009 (and likely several years earlier, as indicated in the dated phylogenies). The evolutionary pressures for this phenomenon are not clear. BDQ was not used in Eswatini at the time of the first DRS and was not introduced on a larger scale until 2015 in the first clinical trials [31]. Given that CFZ featured in an MDR-TB treatment guideline from 2008 [32], we cannot exclude that individual patients were treated with CFZ in Eswatini or neighbouring countries. Alternatively, antifungal azoles, which are commonly used in sub-Saharan Africa, may be responsible [33].

The dominant mutation Rv0678 M146T shared by 21 (53%) gCL10 isolates correlated with MICs to both drugs that were approximately three times higher than for gWT isolates (62% of MICs ≥ 0.75 mg/L). To prove conclusively whether M146T is causative of this increase would require allelic exchange experiments, which were beyond the scope of this study. It is not clear whether these particular MIC increases are clinically significant given that no pharmacokinetic/pharmacodynamic targets for the efficacy of BDQ or CFZ have been set by the European Committee for Antimicrobial Susceptibility Testing (EUCAST) or WHO [9, 34]. Instead, conservative clinical breakpoints (CBs), as defined by EUCAST, have been set that correspond to the likely epidemiologic cut-off values (ECOFFs) of 1 mg/L [9, 35]. Consequently, it was assumed that any MIC increase above the ECOFF would result in worse treatment outcomes, until strong evidence to the contrary is presented. In such a case, the CB could be raised above the ECOFF and/or a second, higher CB could be introduced to define a range of susceptibility at increased exposure (previously known as “intermediate”) [35, 36]. This does not mean that isolates with MICs at or below the current CB are necessarily susceptible given that they could harbour a resistance mechanism that results in only modest MIC increases. To minimize the misclassification of such mutations, an ATU for both drugs could be set and the concentration below the CB would have to be tested, as opposed to just the CB, as currently recommended by WHO (Fig. 3) [36]. Indeed, Nimmo et al. have made a similar proposal for 7H11 DST [33].

In practice, it will be impossible to study the effect of all Rv0678 variants on resistance levels comprehensively, considering the large number of possible mutations with likely heterogeneous effects on MICs (i.e. depending on whether they abolish the function of the repressor completely or not, which is further complicated by the fact that mutations can only confer elevated MICs in strain backgrounds with a functional efflux pump) [37,38,39,40,41]. Moreover, simply correlating clinical outcomes with MICs is complicated by the overlap between MIC distributions. Compared with resistance mechanisms that confer marked MIC increases (e.g. delamanid), BDQ and CFZ will, therefore, pose significant challenges to regulatory agencies and diagnostic laboratories [9, 39, 42]. Nevertheless, the evidence is mounting that particular Rv0678 mutations emerge in patients on failing regimens, confer elevated MICs to both drugs, and adversely affect clinical outcomes [37,38,39,40, 43,44,45,46]. Indeed, Eswatini will likely provide a unique opportunity to study the effect of a single Rv0678 variant on clinical outcomes because the M146T mutation is so frequent in this setting, which would provide invaluable data regarding the efficacy of the all-oral shorter MDR-TB regimen.

A limitation of this study is selection bias, due to exclusion of some (mostly poorly growing) isolates from the genotyping analysis for the strains from the TB-DRS in 2009–2010. Because the excluded isolates were mainly DS, our findings about the MDR-TB clusters are unlikely to be biased significantly. But, any comparison between DS and MDR-TB results has to be interpreted with caution. Another limitation is the fact that the majority of the isolates were from 2009, and the overall number of MDR-Mtbc strains likely represents only about 30% of all MDR-Mtbc strains emerging in that period in Eswatini [4, 5, 47]. This is even more relevant for our second collection that sampled MDR-Mtbc strains from just one region that represents a small fraction (below 5%) of the MDR-Mtbc strains in that period [47]. Still, our analysis of collection 2 isolates and the recently presented data from the recent DRS survey performed in years 2017–2018 underline the current importance of RpoB I491F MDR outbreak cluster isolates in Eswatini [1, 27], but prospective more comprehensive data are needed to detail the ongoing MDR-TB transmission in the country.


From a public health perspective, our findings that transmission and increasing prevalence of highly resistant strains plays a key role in Eswatini’s expanding MDR-TB epidemic underline the urgent need to prioritize breaking transmission chains in the country. This, in turn, means focusing on interventions, such as improving diagnosis of MDR-TB, including RpoB I491F mutants, tracing of MDR-TB patient contacts, implementing effective infection control measures within healthcare facilities and communities, and careful epidemiological monitoring of the spread of MDR outbreak clones on a regional scale and beyond. Moreover, the effect of the Rv0678 M146T mutation on treatment outcomes with BDQ- and CFZ-containing regimens and enhanced transmission rates should be studied as a priority, including using the most recent DRS data.

Our study raises fundamental points regarding the development and introduction of novel anti-TB agents. First, just because an agent is novel, the assumption that isolates are equally susceptible to this agent or that isolates with elevated MICs are rare does not apply in all settings [26, 48]. In fact, it has been recently demonstrated that Rv0678 mutations also confer elevated MICs to OPC-167832, which is currently in phase 1 and 2 trials, illustrating the value of elucidating the genetic basis of resistance [49]. Second, companion DST assays with breakpoints that are set based on modern microbiological principles are needed at the time approval of novel agents, as opposed to years later [39, 45, 50]. Third, if new drugs are used as part of weak regimens, resistance will inevitably evolve and may transmit rapidly. Consequently, molecular and phenotypic DST for BDQ, CFZ, and other drugs need to be scaled up, and further research on resistance mechanisms and the rates of pre-existing resistance in different settings is required, ideally, before novel regimens are introduced. In this context, targeted NGS assays or WGS could play a key role for reflex testing to rapidly rule-in resistance to the vast majority of drugs and enabling real-time surveillance of DR isolates [28, 29, 51].

Availability of data and materials

Fastq files (raw sequencing data) for all strains analysed in this study are available from the EMBL-EBI European Nucleotide Archive sequence read archive ( [52] under the study IDs: PRJNA393767 [12], PRJEB20942 [13], PRJEB5280 [14], PRJNA395592 [15], PRJEB37777 [21], PRJEB6273 [22], PRJEB9680 [23], and PRJEB7281 [24], and details can be found in Additional file 1: Tables S1, S2 and S4.


  1. 1.

    World Health Organization. Global tuberculosis report 2019. 2019.

    Google Scholar 

  2. 2.

    World Health Organisation. Global tuberculosis report 2018. World Health Organization; 2018. Available from: Accessed 20 Nov 2018.

  3. 3.

    O’Donnell MR, Schluger NW. Long walk to treatment for XDR tuberculosis in South Africa. Lancet. 2014;383:1193–4 Available from: [cited 2014 Mar 25].

    Article  Google Scholar 

  4. 4.

    Sanchez-Padilla E, Merker M, Beckert P, Jochims F, Dlamini T, Kahn P, et al. Detection of drug-resistant tuberculosis by Xpert MTB/RIF in Swaziland. N Engl J Med. 2015;372:1181–2.

    CAS  Article  Google Scholar 

  5. 5.

    Sanchez-Padilla E, Dlamini T, Ascorra A, Rüsch-Gerdes S, Tefera ZD, Calain P, et al. High prevalence of multidrug-resistant tuberculosis, Swaziland, 2009-2010. Emerg Infect Dis. 2012;18:29–37.

    Article  Google Scholar 

  6. 6.

    Sikhondze W, Dlamini T, Khumalo D, Maphalala G, Dlamini S, Zikalala T, et al. Countrywide roll-out of Xpert(®) MTB/RIF in Swaziland: the first three years of implementation. Public Health Action. 2015;5:140–6.

    CAS  Article  Google Scholar 

  7. 7.

    Makhado NA, Matabane E, Faccin M, Pinçon C, Jouet A, Boutachkourt F, et al. Outbreak of multidrug-resistant tuberculosis in South Africa undetected by WHO-endorsed commercial tests: an observational study. Lancet Infect Dis. 2018;18:1350–9.

    Article  Google Scholar 

  8. 8.

    Zignol M, Cabibbe AM, Dean AS, Glaziou P, Alikhanova N, Ama C, et al. Genetic sequencing for surveillance of drug resistance in tuberculosis in highly endemic countries: a multi-country population-based surveillance study. Lancet Infect Dis. 2018;18:675–83 Available from: [cited 2020 Mar 19].

    Article  Google Scholar 

  9. 9.

    World Health Organisation. Technical Report on critical concentrations for drug susceptibility testing of medicines used in the treatment of drug-resistant tuberculosis. Geneva; 2018. Available from: Accessed 8 Aug 2019.

  10. 10.

    World Health Organization. WHO consolidated guidelines on drug-resistant tuberculosis treatment. 2019. Available from: [cited 2019 Aug 5].

    Google Scholar 

  11. 11.

    World Health Organization. Rapid communication: key changes to treatment of drug-resistant tuberculosis.2019. Available from:

    Google Scholar 

  12. 12.

    BioProject PRJNA393767. BioProject PRJNA393767. Available from: Accessed 26 Oct 2020.

  13. 13.

    BioProject PRJEB20942. BioProject PRJEB20942. Available from: Accessed 26 Oct 2020.

  14. 14.

    BioProject PRJEB5280. BioProject PRJEB5280. Available from: Accessed 26 Oct 2020.

  15. 15.

    BioProject PRJNA395592. BioProject PRJNA395592. Available from: Accessed 26 Oct 2020.

  16. 16.

    Grandjean L, Gilman RH, Iwamoto T, Köser CU, Coronel J, Zimic M, et al. Convergent evolution and topologically disruptive polymorphisms among multidrug-resistant tuberculosis in Peru. Plos One. 2017;12:e0189838.

    Article  Google Scholar 

  17. 17.

    Xu J, Wang B, Hu M, Huo F, Guo S, Jing W, et al. Primary clofazimine and bedaquiline resistance among isolates from patients with multidrug-resistant tuberculosis. Antimicrob Agents Chemother. 2017;61 Available from: [cited 2019 Aug 1].

  18. 18.

    Weniger T, Krawczyk J, Supply P, Harmsen D, Niemann S. Online tools for polyphasic analysis of Mycobacterium tuberculosis complex genotyping data: now and next. Infect Genet Evol. 2012;12:748–54 Available from: [cited 2013 Aug 7].

    Article  Google Scholar 

  19. 19.

    Merker M, Barbier M, Cox H, Rasigade J-P, Feuerriegel S, Kohl TA, et al. Compensatory evolution drives multidrug-resistant tuberculosis in Central Asia. eLife. 2018;7 Available from: [cited 2019 Mar 27].

  20. 20.

    Roetzer A, Diel R, Kohl TA, Rückert C, Nübel U, Blom J, et al. Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study. Neyrolles O, editor. Plos Med 2013;10:e1001387. Available from: [cited 2013 Aug 7].

  21. 21.

    BioProject PRJEB37777. BioProject PRJEB37777. Available from: Accessed 26 Oct 2020.

  22. 22.

    BioProject PRJEB6273. BioProject PRJEB6273. Available from: Accessed 26 Oct 2020.

  23. 23.

    BioProject PRJEB9680. BioProject PRJEB9680. Available from: Accessed 26 Oct 2020.

  24. 24.

    BioProject PRJEB7281. BioProject PRJEB7281. Available from: Accessed 26 Oct 2020.

  25. 25.

    Valsesia G, Hombach M, Maurer FP, Courvalin P, Roos M, Böttger EC. The resistant-population cutoff (RCOFF): a new concept for improved characterization of antimicrobial susceptibility patterns of non-wild-type bacterial populations. Diekema DJ, editor. J Clin Microbiol 2015;53:1806–1811. Available from: [cited 2019 Aug 9].

  26. 26.

    Merker M, Kohl TA, Barilar I, Andres S, Fowler PW, Chryssanthou E, et al. Phylogenetically informative mutations in genes implicated in antibiotic resistance in Mycobacterium tuberculosis complex. Genome Med. 2020;12:27.

    CAS  Article  Google Scholar 

  27. 27.

    Sikhondze, Welile, Dlamini, Themba, Joloba, Moses, Cirillo, Daniela, Ershova, Julia, Tosas, Olga, et al. XPERT MTB/RIF misses more than 50% of rifampicin cases in Eswatini: results of the 2nd national anti-TB drug resistance survey (2017/2018). Int J Tuberc Lung Dis. 2019;2 3:S585.

  28. 28.

    Gröschel MI, Walker TM, van der Werf TS, Lange C, Niemann S, Merker M. Pathogen-based precision medicine for drug-resistant tuberculosis. PLoS Pathog. 2018;14:e1007297.

    Article  Google Scholar 

  29. 29.

    Feuerriegel S, Kohl TA, Utpatel C, Andres S, Maurer FP, Heyckendorf J, et al. Rapid genomic first- and second-line drug resistance prediction from clinical Mycobacterium tuberculosis specimens using Deeplex®-MycTB. Eur Respir J. 2020:2001796 Available from: [cited 2020 Aug 7].

  30. 30.

    El Achkar S, Demanche C, Osman M, Rafei R, Ismail MB, Yaacoub H, et al. Drug-resistant tuberculosis, Lebanon, 2016–2017. Emerg Infect Dis. 2019;25:564–8 Available from: [cited 2019 Aug 15].

    Article  Google Scholar 

  31. 31.

    Vambe D, Dlamini T, Furin J, Gracia-Edwards C, Keus K, Kunene K, et al. Operational aspects of bedaquiline implementation in Swaziland: report from the field. Public Health Action. 2017;7:240–2 Available from: Accessed 9 Aug 2019.

  32. 32.

    Ministry of Health and Social Welfare, Eswatini. Drug resistant tuberculosis management guidelines and manual. Available from:

  33. 33.

    Nimmo C, Millard J, van Dorp L, Brien K, Moodley S, Wolf A, et al. Population-level emergence of bedaquiline and clofazimine resistance-associated variants among patients with drug-resistant tuberculosis in southern Africa: a phenotypic and phylogenetic analysis. Lancet Microbe. 2020;1:e165–74 Available from: [cited 2020 Aug 10].

    CAS  Article  Google Scholar 

  34. 34.

    European Committee for Antimicrobial Susceptibility Testing. Rationale for EUCAST clinical breakpoints bedaquiline. 2019. Available from:

    Google Scholar 

  35. 35.

    Kahlmeter G. The 2014 Garrod Lecture: EUCAST – are we heading towards international agreement? J Antimicrob Chemother. 2015;70:2427–39 Available from: [cited 2019 Aug 9].

    Article  Google Scholar 

  36. 36.

    Kahlmeter G, Giske CG, Kirn TJ, Sharp SE. Point-counterpoint: differences between the European Committee on Antimicrobial Susceptibility Testing and the Clinical Laboratory Standards Institute recommendations for reporting antimicrobial susceptibility results. J Clin Microbiol. 2019; Available from: [cited 2019 Aug 9].

  37. 37.

    Villellas C, Coeck N, Meehan CJ, Lounis N, de Jong B, Rigouts L, et al. Unexpected high prevalence of resistance-associated Rv0678 variants in MDR-TB patients without documented prior use of clofazimine or bedaquiline. J Antimicrob Chemother. 2017;72:684–90.

    CAS  PubMed  Google Scholar 

  38. 38.

    Nguyen TVA, Anthony RM, Bañuls A-L, Nguyen TVA, Vu DH, Alffenaar J-WC. Bedaquiline resistance: its emergence, mechanism, and prevention. Clin Infect Dis. 2018;66:1625–30.

    Article  Google Scholar 

  39. 39.

    Kranzer K, Kalsdorf B, Heyckendorf J, Andres S, Merker M, Hofmann-Thiel S, et al. New WHO treatment recommendations for multidrug-resistant tuberculosis: are we well enough prepared? Am J Respir Crit Care Med. 2019;200:514–5.

    Article  Google Scholar 

  40. 40.

    Ismail NA, Omar SV, Joseph L, Govender N, Blows L, Ismail F, et al. Defining bedaquiline susceptibility, resistance, cross-resistance and associated genetic determinants: a retrospective cohort study. EBioMedicine. 2018;28:136–42.

    Article  Google Scholar 

  41. 41.

    Kadura S, King N, Nakhoul M, Zhu H, Theron G, Köser CU, et al. Systematic review of mutations associated with resistance to the new and repurposed Mycobacterium tuberculosis drugs bedaquiline, clofazimine, linezolid, delamanid and pretomanid. J Antimicrob Chemother. 2020:dkaa136 Available from: [cited 2020 Aug 10].

  42. 42.

    Köser CU, Maurer FP, Kranzer K. ‘Those who cannot remember the past are condemned to repeat it’: drug-susceptibility testing for bedaquiline and delamanid. Int J Infect Dis. 2019;80:S32–5 Available from: [cited 2020 Aug 11].

    Article  Google Scholar 

  43. 43.

    de Vos M, Ley SD, Wiggins KB, Derendinger B, Dippenaar A, Grobbelaar M, et al. Bedaquiline microheteroresistance after cessation of tuberculosis treatment. N Engl J Med. 2019;380:2178–80.

    Article  Google Scholar 

  44. 44.

    Polsfuss S, Hofmann-Thiel S, Merker M, Krieger D, Niemann S, Rüssmann H, et al. Emergence of low-level delamanid and bedaquiline resistance during extremely drug-resistant tuberculosis treatment. Clin Infect Dis. 2019;69:1229–31.

    CAS  Article  Google Scholar 

  45. 45.

    Andres S, Merker M, Heyckendorf J, Kalsdorf B, Rumetshofer R, Indra A, et al. Bedaquiline-resistant tuberculosis: dark clouds on the horizon. Am J Respir Crit Care Med. 2020;201:1564–8.

    Article  Google Scholar 

  46. 46.

    Nimmo C, Millard J, Brien K, Moodley S, van Dorp L, Lutchminarain K, et al. Bedaquiline resistance in drug-resistant tuberculosis HIV co-infected patients. Eur Respir J. 2020;55:1902383.

    Article  Google Scholar 

  47. 47.

    World Health Organization. Global health observatory data repository 2020. Available from: MDRTB?lang=en. [cited 2020 Oct 8].

    Google Scholar 

  48. 48.

    Köser CU, Feuerriegel S, Summers DK, Archer JAC, Niemann S. Importance of the genetic diversity within the Mycobacterium tuberculosis complex for the development of novel antibiotics and diagnostic tests of drug resistance. Antimicrob Agents Chemother. 2012;56:6080–7.

    Article  Google Scholar 

  49. 49.

    Hariguchi N, Chen X, Hayashi Y, Kawano Y, Fujiwara M, Matsuba M, et al. OPC-167832, a novel carbostyril derivative with potent anti-tuberculosis activity as a DprE1 inhibitor. Antimicrob Agents Chemother 2020;AAC.02020–19, aac;AAC.02020–19v1. Available from: [cited 2020 Apr 2].

  50. 50.

    Schön T, Köser CU, Werngren J, Viveiros M, Georghiou S, Kahlmeter G, et al. Commentary: what is the role of the EUCAST reference method for MIC testing of the Mycobacterium tuberculosis complex? Clin Microbiol Infect. 2020;S1198743X20304456. Available from: [cited 2020 Aug 10].

  51. 51.

    Jouet A, Gaudin C, Badalato N, Allix-Béguec C, Duthoy S, Ferré A, et al. Deep amplicon sequencing for culture-free prediction of susceptibility or resistance to 13 anti-tuberculous drugs. Eur Respir J. 2020:2002338 Available from: [cited 2020 Oct 8].

  52. 52.

    Leinonen R, Akhtar R, Birney E, Bower L, Cerdeno-Tarraga A, Cheng Y, et al. The European nucleotide archive. Nucleic Acids Res. 2011;39:D28–31 Available from: [cited 2020 Oct 23].

    CAS  Article  Google Scholar 

Download references


We thank J. Zallet, I. Razio, T. Struve-Sonnenschein, V. Mohr, and T. Niemann, from Borstel, Germany, for the technical assistance. We would also like to thank all partners in the National Tuberculosis Control Programme and Médecins Sans Frontières-Switzerland in Eswatini.


Parts of this work have been supported by the European Union TB-PAN-NET (FP7-223681) project, by Médecins Sans Frontières-Switzerland, and by German Center for Infection Research, Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germanys Excellence Strategy – EXC 2167, and Leibniz Science Campus Evolutionary Medicine of the LUNG (EvoLUNG). The funders had no role in the study design; in the collection, analysis, and interpretation of the data; in the writing of the report; and in the decision to submit the paper for publication.

Author information




PB, ESP, MM, VD, TAK, CU, CUK, IB, NI, SHO, MK, RMW, HH, GM, EA, BDJ, BK, BS, SA, KK, FPM, MB, and SN conceived the idea and designed the study and analysed and interpreted the data. NI, SVO, MK, EA, BSJ, and RMW provided study materials or patients. All authors contributed to obtaining and assembling the data. PB, ESP, MM, MB, and SN wrote the initial draft of the paper. All authors contributed to the data interpretation and final draft of the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Stefan Niemann.

Ethics declarations

Ethics approval and consent to participate

Collection 1: Ethical approval was obtained from the Ministry of Health Scientific and Ethics Committee of Eswatini/Swaziland and the Ethics Review Board of Médecins Sans Frontières. Inclusion in the study was voluntary and after signing of an informed consent form. Collection 2: The study protocol was approved by the Institutional Review Board of the Institute of Tropical Medicine; the Ethics Committee of the University Hospital of Antwerp, Belgium; and the Ministry of Health Scientific and Ethics Committee of Eswatini. Inclusion in the study was voluntary and after signing of an informed consent form. The research performed conformed to the principles expressed in the Helsinki Declaration.

Consent for publication

Not applicable

Competing interests

Dr. Niemann reports grants from German Center for Infection Research, grants from European Union TB-PAN-NET, grants from Excellenz Cluster Precision Medicine in Chronic Inflammation EXC 2167, and grants from Leibniz Science Campus Evolutionary Medicine of the LUNG (EvoLUNG), during the conduct of the study; Dr. Köser reports personal fees from the World Health Organization Regional Office for Europe, personal fees from Becton Dickinson, personal fees from QuantuMDx Group Ldt, personal fees from Foundation for Innovative New Diagnostics, others from Hain Lifescience, others from Global Alliance for TB Drug Development Inc. and Otsuka Novel Products GmbH, and others from YD Diagnostics, outside the submitted work, and is an unpaid advisor to GenoScreen.

The remaining authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Origin of Rv0678 mutant isolates included in homoplastic analysis. Table S2. Drug resistance mutations and cluster analysis of collection 2 strains and EMBL-EBI ENA sequence read archive numbers. Table S3. BDQ and CFZ Minimum inhibitory concentration of 25 Rv0678 mutant and 12 wild type strains. Table S4. WGS of 273 isolates submitted to the EMBL-EBI ENA sequence read archive. Table S5. Epidemiological and genotyping data of 412 Mtbc-strains. Table S6. Phenotypic DST data of 412 Mtbc-strains, resistance mediating mutations for 273 whole genome sequenced isolates. Table S7. Drug resistance mediating mutations of all MDR clusters.

Additional file 2.

Supplemental methods and description of resistance mutations.

Additional file 3: Fig. S1.

Phylogenetic diversity and drug susceptibility of the 412 M. tuberculosis complex isolates from Eswatini. Fig. S2. Lineage distribution across Eswatini. Fig. S3. Isolates selected for further whole-genome sequencing analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Beckert, P., Sanchez-Padilla, E., Merker, M. et al. MDR M. tuberculosis outbreak clone in Eswatini missed by Xpert has elevated bedaquiline resistance dated to the pre-treatment era. Genome Med 12, 104 (2020).

Download citation


  • Tuberculosis
  • Multidrug resistance
  • Resistance evolution
  • MDR outbreak strains
  • Diagnostice escape
  • Treatment escape
  • Treatment failure