Search

Highlights • Unresolved previous infection as major cause of recurring tuberculosis (TB) in Ghana. • Genomic epidemiology identifies high relapse among recurrent TB cases in Ghana. • 15-locus MIRU-VNTR typing is sufficient to predict the cause of TB recurrence. • Evidence of recent household-related TB transmission in Ghana. • Need for increased education by national TB control program., Objective To retrospectively investigate the cause of recurring tuberculosis (rcTB) among participants with pulmonary TB recruited from a prospective population-based study conducted between July 2012 and December 2015. Methods Mycobacterium tuberculosis complex isolates obtained from rcTB cases were characterized by standard mycobacterial genotyping tools, whole-genome sequencing, and phylogenetic analysis carried out to assess strain relatedness. Results The majority (58.3%, 21/36) of study participants with rcTB episodes had TB recurrence within 12 months post treatment. TB strains with isoniazid (INH) resistance were found in 19.4% (7/36) of participants at the primary episode, of which 29% (2/7) were also rifampicin-resistant. On TB recurrence, an INH-resistant strain was found in a larger proportion of participants, 27.8% (10/36), of which 40% (4/10) were MDR-TB strains. rcTB was attributed to relapse (same strain) in 75.0% (27/36) of participants and 25.0% (9/36) to re-infection. Conclusion Our findings indicate that previous unresolved infectiondue to inadequate treatment, may be the major cause of rcTB.

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., Journal of Clinical Microbiology, 60 (3), ISSN:0095-1137, ISSN:1098-660X

International audience; Abstract Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

Summary Background South Africa has a high burden of rifampicin-resistant tuberculosis (including multidrug-resistant [MDR] tuberculosis), with increasing rifampicin-monoresistant (RMR) tuberculosis over time. Resistance acquisition during first-line tuberculosis treatment could be a key contributor to this burden, and HIV might increase the risk of acquiring rifampicin resistance. We assessed whether HIV during previous treatment was associated with RMR tuberculosis and resistance acquisition among a retrospective cohort of patients with MDR or rifampicin-resistant tuberculosis. Methods In this retrospective cohort study, we included all patients routinely diagnosed with MDR or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa, between Jan 1, 2008, and Dec 31, 2017. Patient-level data were obtained from a prospective database, complemented by data on previous tuberculosis treatment and HIV from a provincial health data exchange. Stored MDR or rifampicin-resistant tuberculosis isolates from patients underwent whole-genome sequencing (WGS). WGS data were used to infer resistance acquisition versus transmission, by identifying genomically unique isolates (single nucleotide polymorphism threshold of five). Logistic regression analyses were used to assess factors associated with RMR tuberculosis and genomic uniqueness. Findings The cohort included 2041 patients diagnosed with MDR or rifampicin-resistant tuberculosis between Jan 1, 2008, and Dec 31, 2017; of those, 463 (22·7%) with RMR tuberculosis and 1354 (66·3%) with previous tuberculosis treatment. In previously treated patients, HIV positivity during previous tuberculosis treatment versus HIV negativity (adjusted odds ratio [OR] 2·07, 95% CI 1·35–3·18), and three or more previous tuberculosis treatment episodes versus one (1·96, 1·21–3·17) were associated with RMR tuberculosis. WGS data showing MDR or rifampicin-resistant tuberculosis were available for 1169 patients; 360 (30·8%) isolates were identified as unique. In previously treated patients, RMR tuberculosis versus MDR tuberculosis (adjusted OR 4·96, 3·40–7·23), HIV positivity during previous tuberculosis treatment (1·71, 1·03–2·84), and diagnosis in 2013–17 (1·42, 1·02–1·99) versus 2008–12, were associated with uniqueness. In previously treated patients with RMR tuberculosis, HIV positivity during previous treatment (adjusted OR 5·13, 1·61–16·32) was associated with uniqueness as was female sex (2·50 [1·18–5·26]). Interpretation These data suggest that HIV contributes to rifampicin-resistance acquisition during first-line tuberculosis treatment and that this might be driving increasing RMR tuberculosis over time. Large-scale prospective cohort studies are required to further quantify this risk. Funding Swiss National Science Foundation, South African National Research Foundation, and Wellcome Trust.

Mycobacterial Interspersed Repetitive Units-Variable Tandem Repeats (MIRU-VNTR) typing has been widely used for molecular epidemiological studies of tuberculosis (TB). However, genotyping tools for Mycobacterium tuberculosis (Mtb) may be limiting in some settings due to high cost and workload. In this study developed a customized stepwise MIRU-VNTR typing that prioritizes high discriminatory loci and validated this method using penitentiary system cohort in the country of Georgia.We used a previously generated MIRU-VNTR dataset from recurrent TB cases (32 cases) in Georgia and a new dataset of TB cases from the penitentiary system (102 cases) recruited from 2014 to 2015. A Hunter-Gaston Discriminatory Index (HGDI) was calculated utilizing a 24 standard loci panel, to select high discriminatory power loci, subsequently defined as the customized Georgia-specific set of loci for initial typing. The remaining loci were scored and hierarchically grouped for second and third step typing of the cohort. We then compared the processing time and costs of the customized stepwise method to the standard 24-loci method.For the customized Georgia-specific set that was used for initial typing, 10 loci were selected with a minimum value of 0.32 to the highest HGDI score locus. Customized 10 loci (step 1) typing of 102 Mtb patient isolates revealed 35.7% clustered cases. This proportion was reduced to 19.5% after hierarchical application of 2nd and 3rd step typing with the corresponding groups of loci. Our customized stepwise MIRU-VNTR genotyping approach reduced the quantity of samples to be typed and therefore overall processing time and costs by 42.6% each.Our study shows that our customized stepwise MIRU-VNTR typing approach is a valid alternative of standard MIRI-VNTR typing panels for molecular epidemiological investigation in Georgia that saves time, workload and costs. Similar approaches could be developed for other settings.

Coscolla et al., Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto, as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum . Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum , and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally., M.C. is supported by the Ramón y Cajal programme from the Ministerio de Ciencia, Innovación y Universidades. This work was supported by the European Society of Clinical Microbiology and Infectious Diseases (ESCMID) (research award to M.C.), Ministerio de Ciencia, Innovación y Universidades (grant number RTI2018-094399-A-I00 to M.C.) and Consellería de Educació de la Generalitat Valenciana (grant number SEJI/2019/011 to M.C.), the Swiss National Science Foundation (grants 310030_188888, IZRJZ3_164171, IZLSZ3_170834 and CRSII5_177163 to S. G.), the European Research Council (883582-ECOEVODRTB to S. G.) and Wellcome (grant number 097134/Z/11/Z to D. Y.-M).

24 páginas, 3 figuras, 1 tabla. The sequence data generated by this study has been deposited on SRA (https://www.ncbi.nlm.nih.gov/sra) under the accession number PRJNA670836. Extended data is available here: https://github.com/fmenardo/MTBC_L1_L3. DOI: https://doi.org/10.5281/zenodo.4609804 (Menardo, 2021)., Background: Lineage 1 (L1) and 3 (L3) are two lineages of the Mycobacterium tuberculosis complex (MTBC) causing tuberculosis (TB) in humans. L1 and L3 are prevalent around the rim of the Indian Ocean, the region that accounts for most of the world's new TB cases. Despite their relevance for this region, L1 and L3 remain understudied. Methods: We analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We reconstructed the evolutionary history of these two lineages and identified genes under positive selection. Results: We found a strongly asymmetric pattern of migration from South Asia toward neighboring regions, highlighting the historical role of South Asia in the dispersion of L1 and L3. Moreover, we found that several genes were under positive selection, including genes involved in virulence and resistance to antibiotics. For L1 we identified signatures of local adaptation at the esxH locus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. Conclusions: Our study highlights the importance of genetic diversity in the MTBC, and sheds new light on two of the most important MTBC lineages affecting humans., This work was supported by the Swiss National Science Foundation (grants 310030_188888, CRSII5_177163, IZRJZ3_164171 and IZLSZ3_170834) and the European Research Council (309540 EVODRTB and 883582-ECOEVODRTB).

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 mu g/ml (range, 0.125 to 1 mu 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.

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients’ lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC., Heteroresistance is an important stepping-stone to how an initially monoclonal and drug-susceptible population of Mycobacterium tuberculosis becomes fully resistant to a given antimicrobial during the course of an infection, and this review discusses how bacterial mutation rates, bacterial population size, the number of mutations that can confer antimicrobial resistance (i.e. AMR target size) and the fitness of AMR mutations all modulate the magnitude of heteroresistance.

Mutations in the genes of the F420 signaling pathway of Mycobacterium tuberculosis complex, including dnn, fgd1, fbiA, fbiB, fbiC, and fbiD, can lead to delamanid resistance. We searched for such mutations among 129 M. tuberculosis strains from Asia, South America, and Africa using whole-genome sequencing; 70 (54%) strains had at least one mutation in one of the genes., Mutations in the genes of the F420 signaling pathway of Mycobacterium tuberculosis complex, including dnn, fgd1, fbiA, fbiB, fbiC, and fbiD, can lead to delamanid resistance. We searched for such mutations among 129 M. tuberculosis strains from Asia, South America, and Africa using whole-genome sequencing; 70 (54%) strains had at least one mutation in one of the genes. For 10 strains with mutations, we determined the MIC of delamanid. We found one strain from a delamanid-naive patient carrying the natural polymorphism Tyr29del (ddn) that was associated with a critical delamanid MIC.

Human tuberculosis is caused by members of the Mycobacterium tuberculosis complex (MTBC) that vary in virulence and transmissibility. While genome-wide association studies have uncovered several mutations conferring drug resistance, much less is known about the factors underlying other bacterial phenotypes. Variation in the outcome of tuberculosis infection and diseases has been attributed primarily to patient and environmental factors, but recent evidence indicates an additional role for the genetic diversity among MTBC clinical strains. Here, we used metabolomics to unravel the effect of genetic variation on the strain-specific metabolic adaptive capacity and vulnerability. To define the functionality of single-nucleotide polymorphisms (SNPs) systematically, we developed a constraint-based approach that integrates metabolomic and genomic data. Our model-based predictions correctly classify SNP effects in pyruvate kinase and suggest a genetic basis for strain-specific inherent baseline susceptibility to the antibiotic para-aminosalicylic acid. Our method is broadly applicable across microbial life, opening possibilities for the development of more selective treatment strategies.

Summary Background Drug resistance is a challenge for the global control of tuberculosis. We examined mortality in patients with tuberculosis from high-burden countries, according to concordance or discordance of results from drug susceptibility testing done locally and in a reference laboratory. Methods This multicentre cohort study was done in Cote d'Ivoire, Democratic Republic of the Congo, Kenya, Nigeria, South Africa, Peru, and Thailand. We collected Mycobacterium tuberculosis isolates and clinical data from adult patients aged 16 years or older. Patients were stratified by HIV status and tuberculosis drug resistance. Molecular or phenotypic drug susceptibility testing was done locally and at the Swiss National Center for Mycobacteria, Zurich, Switzerland. We examined mortality during treatment according to drug susceptibility test results and treatment adequacy in multivariable logistic regression models adjusting for sex, age, sputum microscopy, and HIV status. Findings We obtained M tuberculosis isolates from 871 patients diagnosed between 2013 and 2016. After exclusion of 237 patients, 634 patients with tuberculosis were included in this analysis; the median age was 33·2 years (IQR 26·9–42·5), 239 (38%) were women, 272 (43%) were HIV-positive, and 69 (11%) patients died. Based on the reference laboratory drug susceptibility test, 394 (62%) strains were pan-susceptible, 45 (7%) monoresistant, 163 (26%) multidrug-resistant (MDR), and 30 (5%) had pre-extensively or extensively drug resistant (pre-XDR or XDR) tuberculosis. Results of reference and local laboratories were concordant for 513 (81%) of 634 patients and discordant for 121 (19%) of 634. Overall, sensitivity to detect any resistance was 90·8% (95% CI 86·5–94·2) and specificity 84·3% (80·3–87·7). Mortality ranged from 6% (20 of 336) in patients with pan-susceptible tuberculosis treated according to WHO guidelines to 57% (eight of 14) in patients with resistant strains who were under-treated. In logistic regression models, compared with concordant drug susceptibility test results, the adjusted odds ratio of death was 7·33 (95% CI 2·70–19·95) for patients with discordant results potentially leading to under-treatment. Interpretation Inaccurate drug susceptibility testing by comparison with a reference standard leads to under-treatment of drug-resistant tuberculosis and increased mortality. Rapid molecular drug susceptibility test of first-line and second-line drugs at diagnosis is required to improve outcomes in patients with MDR tuberculosis and pre-XDR or XDR tuberculosis. Funding National Institutes of Allergy and Infectious Diseases, Swiss National Science Foundation, Swiss National Center for Mycobacteria.

TheMycobacterium tuberculosiscomplex (MTBC) causes tuberculosis (TB) in humans and various other mammals. The human-adapted members of the MTBC comprise seven phylogenetic lineages that differ in their geographical distribution. There is growing evidence that this phylogenetic diversity modulates the outcome of TB infection and disease. For decades, TB research and development has focused on the two canonical MTBC reference strains H37Rv and Erdman, both of which belong to Lineage 4. Relying on only a few laboratory-adapted strains can be misleading as study results might not be directly transferrable to clinical settings where patients are infected with a diverse array of strains, including drug-resistant variants. Here, we argue for the need to expand TB research and development by incorporating the phylogenetic diversity of the MTBC. To facilitate such work, we have assembled a group of 20 genetically well-characterized clinical strains representing the seven known human-adapted MTBC lineages. With the “MTBC clinical strains reference set” we aim to provide a standardized resource for the TB community. We hope it will enable more direct comparisons between studies that explore the physiology of MTBC beyond the lab strains used thus far. We anticipate that detailed phenotypic analyses of this reference strain set will increase our understanding of TB biology and assist in the development of new control tools that are universally effective.

Introduction Recurrent tuberculosis (TB) is one of the main challenges in TB control. Genotyping based on Mycobacterial Interspersed Repetitive Units–Variable Tandem Repeats (MIRU-VNTR) has been widely used to differentiate between relapse and reinfection, which are the two main causes of recurrent TB. There is a lack of data regarding the causes of TB recurrence in Georgia, and while differentiating between relapse and reinfection plays a key role in defining appropriate interventions, the required genotyping methodologies have not been implemented. The objective of this study was to implement MIRU-VNTR genotyping at the National Center for Tuberculosis and Lung Diseases (NCTBLD) and differentiate between relapse and reinfection in multidrug resistant (MDR-) TB patients from Tbilisi, Georgia. Methods Recurrent MDR tuberculosis cases from 2014–2016 diagnosed at NCTLD were included in the study when bacterial samples from both episodes were available. Genotyping based on the MIRU-VNTR 24 loci was implemented and used for differentiating between relapse and reinfection. Paired samples showing the same MIRU-VNTR pattern or one locus difference were classified as relapse, while two and more loci differences were treated as reinfection. Exact logistic regression was used to identify predictors of recurrence. Results Thirty two MDR-TB patients (64 samples) were included and MIRU-VNTR 24 typing was performed on the corresponding paired samples. Of the 32 patients, 25 (83.3%) were identified as relapse while 5 (16.7%) were due to re-infection. Patients with a history of incarceration were significantly associated with TB reinfection (p< 0.05). Conclusion Recurrent TB in MDR patients in Georgia are mainly caused by relapse, raising concerns on the efficacy of the TB control program. An association between incarceration and reinfection likely reflects high levels of ongoing TB transmission in prisons, indicating the need for better TB infection control measures in these settings. Our results add to the rationale for implementing genotypic surveillance of TB more broadly to support TB control in Georgia.

BackgroundHuman tuberculosis (TB) is caused by seven phylogenetic lineages of the Mycobacterium tuberculosis complex (MTBC), Lineage 1-7. Recent advances in rapid genotyping of MTBC based on single nucleotide polymorphisms (SNP), allow for phylogenetically robust strain classification, paving the way for defining genotype-phenotype relationships in clinical settings. Such studies have revealed that, in addition to host and environmental factors, strain variation in the MTBC influences the outcome of TB infection and disease. In Tanzania, such molecular epidemiological studies of TB however are scarce in spite of a high TB burden.Methods and findingsHere we used SNP-typing to characterize a nationwide collection of 2,039 MTBC clinical isolates representative of 1.6% of all new and retreatment TB cases notified in Tanzania during 2012 and 2013. Four lineages, namely Lineage 1-4 were identified within the study population. The distribution and frequency of these lineages varied across regions but overall, Lineage 4 was the most frequent (n = 866, 42.5%), followed by Lineage 3 (n = 681, 33.4%) and 1 (n = 336, 16.5%), with Lineage 2 being the least frequent (n = 92, 4.5%). We found Lineage 2 to be independently associated with female sex (adjusted odds ratio [aOR] 2.14; 95% confidence interval [95% CI] 1.31 - 3.50, p = 0.002) and retreatment cases (aOR 1.67; 95% CI 0.95 - 2.84, p = 0. 065) in the study population. We found no associations between MTBC lineage and patient age or HIV status. Our sublineage typing based on spacer oligotyping on a subset of Lineage 1, 3 and 4 strains revealed the presence of mainly EAI, CAS and LAM families. Finally, we detected low levels of multidrug resistant isolates among a subset of 144 retreatment cases.ConclusionsThis study provides novel insights into the MTBC lineages and the possible influence of pathogen-related factors on the TB epidemic in Tanzania.

Immigrants from regions with a high incidence of tuberculosis (TB) are a risk group for TB in low-incidence countries such as Switzerland. In a previous analysis of a nationwide collection of 520 Mycobacterium tuberculosis isolates from 2000 to 2008, we identified 35 clusters comprising 90 patients based on standard genotyping (24-locus mycobacterial interspersed repetitive-unit–variable-number tandem-repeat [MIRU-VNTR] typing and spoligotyping). Here, we used whole-genome sequencing (WGS) to revisit these transmission clusters. Genome-based transmission clusters were defined as isolate pairs separated by ≤12 single nucleotide polymorphisms (SNPs). WGS confirmed 17/35 (49%) MIRU-VNTR typing clusters; the other 18 clusters contained pairs separated by >12 SNPs. Most transmission clusters (3/4) of Swiss-born patients were confirmed by WGS, as opposed to 25% (4/16) of the clusters involving only foreign-born patients. The overall clustering proportion was 17% (90 patients; 95% confidence interval [CI], 14 to 21%) by standard genotyping but only 8% (43 patients; 95% CI, 6 to 11%) by WGS. The clustering proportion was 17% (67/401; 95% CI, 13 to 21%) by standard genotyping and 7% (26/401; 95% CI, 4 to 9%) by WGS among foreign-born patients and 19% (23/119; 95% CI, 13 to 28%) and 14% (17/119; 95% CI, 9 to 22%), respectively, among Swiss-born patients. Using weighted logistic regression, we found weak evidence of an association between birth origin and transmission (adjusted odds ratio of 2.2 and 95% CI of 0.9 to 5.5 comparing Swiss-born patients to others). In conclusion, standard genotyping overestimated recent TB transmission in Switzerland compared to WGS, particularly among immigrants from regions with a high TB incidence, where genetically closely related strains often predominate. We recommend the use of WGS to identify transmission clusters in settings with a low incidence of TB.

Whole-genome sequencing allows rapid detection of drug-resistant

Background Large sequence datasets are difficult to visualize and handle. Additionally, they often do not represent a random subset of the natural diversity, but the result of uncoordinated and convenience sampling. Consequently, they can suffer from redundancy and sampling biases. Results Here we present Treemmer, a simple tool to evaluate the redundancy of phylogenetic trees and reduce their complexity by eliminating leaves that contribute the least to the tree diversity. Conclusions Treemmer can reduce the size of datasets with different phylogenetic structures and levels of redundancy while maintaining a sub-sample that is representative of the original diversity. Additionally, it is possible to fine-tune the behavior of Treemmer including any kind of meta-information, making Treemmer particularly useful for empirical studies., BMC Bioinformatics, 19, ISSN:1471-2105

Drug resistance is best thought of as an ongoing biological process. Resistant bacteria must emerge, become established and ultimately transmit in order to be relevant to human health. In this context, genetic diversity can influence the rate and likelihood of resistance emerging; it can also modulate the net physiological impact of resistance and the propensity of an organism to improve any defects that arise from it. Combined, these effects can have an impact on a larger scale, with highly transmissible drug-resistant bacterial strains posing a formidable threat to global health. These considerations are pertinent to the future of tuberculosis control as well. In this chapter, we review our current understanding of the impact of genetic diversity in the broadest sense on the evolution of drug-resistant members of the Mycobacterium tuberculosis complex.

Background. Whole-genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single-nucleotide polymorphism (SNP) typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. Methods. On the basis of genome sequences of 3 historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1642 patient isolates and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. Results. We identified 68 patients associated with the outbreak strain. Most received a tuberculosis diagnosis in 1991-1995, but cases were observed until 2011. Two thirds were homeless and/or substance abusers. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into 3 subclusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. Conclusions. Strain-specific SNP genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real time and at high resolution

Antibiotic-resistant Mycobacterium tuberculosis strains are threatening progress in containing the global tuberculosis epidemic. Mycobacterium tuberculosis is intrinsically resistant to many antibiotics, limiting the number of compounds available for treatment. This intrinsic resistance is due to a number of mechanisms including a thick, waxy, hydrophobic cell envelope and the presence of drug degrading and modifying enzymes. Resistance to the drugs which are active against M. tuberculosis is, in the absence of horizontally transferred resistance determinants, conferred by chromosomal mutations. These chromosomal mutations may confer drug resistance via modification or overexpression of the drug target, as well as by prevention of prodrug activation. Drug resistance mutations may have pleiotropic effects leading to a reduction in the bacterium's fitness, quantifiable e.g. by a reduction in the in vitro growth rate. Secondary so-called compensatory mutations, not involved in conferring resistance, can ameliorate the fitness cost by interacting epistatically with the resistance mutation. Although the genetic diversity of M. tuberculosis is low compared to other pathogenic bacteria, the strain genetic background has been demonstrated to influence multiple aspects in the evolution of drug resistance. The rate of resistance evolution and the fitness costs of drug resistance mutations may vary as a function of the genetic background.

Recent surveillance data of multidrug-resistant tuberculosis (MDR-TB) reported the highest rates of resistance ever documented. As further amplification of resistance in MDR strains of Mycobacterium tuberculosis occurs, extensively drug-resistant (XDR) and totally drug resistant (TDR) TB are beginning to emerge. Whilst for the most part, the epidemiological factors involved in the spread of MDR-TB are understood, insights into the bacterial drivers of MDR-TB have been gained only recently, largely owing to novel technologies and research in other organisms. Herein, we review recent findings on how bacterial factors such as persistence, hypermutation, the complex interrelationship between drug resistance and fitness, compensatory evolution, and epistasis affect the evolution of multidrug resistance in M. tuberculosis. Improved knowledge of these factors will help better predict the future trajectory of MDR-TB, and contribute to the development of new tools and strategies to combat this growing public health threat.

Rifampin resistance in clinical isolates of Mycobacterium tuberculosis arises primarily through the selection of bacterial variants harboring mutations in the 81-bp rifampin resistance-determining region of the rpoB gene. While these mutations were shown to infer a fitness cost in the absence of antibiotic pressure, compensatory mutations in rpoA and rpoC were identified which restore the fitness of rifampin-resistant bacteria carrying mutations in rpoB . To investigate the epidemiological relevance of these compensatory mutations, we analyzed 286 drug-resistant and 54 drug-susceptible clinical M. tuberculosis isolates from the Western Cape, South Africa, a high-incidence setting of multidrug-resistant tuberculosis. Sequencing of a portion of the RpoA-RpoC interaction region of the rpoC gene revealed that 23.5% of all rifampin-resistant isolates tested carried a nonsynonymous mutation in this region. These putative compensatory mutations in rpoC were associated with transmission, as 30.8% of all rifampin-resistant isolates with an IS 6110 restriction fragment length polymorphism (RFLP) pattern belonging to a recognized RFLP cluster harbored putative rpoC mutations. Such mutations were present in only 9.4% of rifampin-resistant isolates with unique RFLP patterns ( P < 0.01). Moreover, these putative compensatory mutations were associated with specific strain genotypes and the rpoB S531L rifampin resistance mutation. Among isolates harboring this rpoB mutation, 44.1% also harbored rpoC mutations, while only 4.1% of the isolates with other rpoB mutations exhibited mutations in rpoC ( P < 0.001). Our study supports a role for rpoC mutations in the transmission of multidrug-resistant tuberculosis and illustrates how epistatic interactions between drug resistance-conferring mutations, compensatory mutations, and different strain genetic backgrounds might influence compensatory evolution in drug-resistant M. tuberculosis .

Generalist and specialist species differ in the breadth of their ecological niches. Little is known about the niche width of obligate human pathogens. Here we analyzed a global collection of Mycobacterium tuberculosis lineage 4 clinical isolates, the most geographically widespread cause of human tuberculosis. We show that lineage 4 comprises globally distributed and geographically restricted sublineages, suggesting a distinction between generalists and specialists. Population genomic analyses showed that, whereas the majority of human T cell epitopes were conserved in all sublineages, the proportion of variable epitopes was higher in generalists. Our data further support a European origin for the most common generalist sublineage. Hence, the global success of lineage 4 reflects distinct strategies adopted by different sublineages and the influence of human migration., We thank S. Lecher, S. Li and J. Zallet for technical support. Calculations were performed at the sciCORE scientific computing core facility at the University of Basel. This work was supported by the Swiss National Science Foundation (grants 310030_166687 (S.G.) and 320030_153442 (M.E.) and Swiss HIV Cohort Study grant 740 to L.F.), the European Research Council (309540-EVODRTB to S.G.), TB-PAN-NET (FP7-223681 to S.N.), PathoNgenTrace projects (FP7-278864-2 to S.N.), SystemsX.ch (S.G.), the German Center for Infection Research (DZIF; S.N.), the Novartis Foundation (S.G.), the Natural Science Foundation of China (91631301 to Q.G.), and the National Institute of Allergy and Infectious Diseases (5U01-AI069924-05) of the US National Institutes of Health (M.E.).

Background Mycobacterium africanum comprises two phylogenetic lineages within the M. tuberculosis complex (MTBC) and is an important cause of human tuberculosis (TB) in West Africa. The reasons for this geographic restriction of M. africanum remain unclear. Here, we performed a prospective study to explore associations between the characteristics of TB patients and the MTBC lineages circulating in Ghana. Method We genotyped 1,211 MTBC isolates recovered from pulmonary TB patients recruited between 2012 and 2014 using single nucleotide polymorphism typing and spoligotyping. Associations between patient and pathogen variables were assessed using univariate and multivariate logistic regression. Results Of the 1,211 MTBC isolates analysed, 71.9 % (871) belonged to Lineage 4; 12.6 % (152) to Lineage 5 (also known as M. africanum West-Africa 1), 9.2 % (112) to Lineage 6 (also known as M. africanum West-Africa 2) and 0.6 % (7) to Mycobacterium bovis. Univariate analysis revealed that Lineage 6 strains were less likely to be isoniazid resistant compared to other strains (odds ratio = 0.25, 95 % confidence interval (CI): 0.05–0.77, P

SETTING: Madang and surroundings, Papua New Guinea (PNG). OBJECTIVE: To characterise the genetic diversity and drug susceptibility of Mycobacterium tuberculosis isolates collected in Madang and surroundings. DESIGN: M. tuberculosis was isolated from sputum samples from active pulmonary tuberculosis cases. Drug resistance profiles were obtained by drug susceptibility testing. M. tuberculosis lineages were identified by single nucleotide polymorphisms and sub-typing was performed by spoligotyping. Spoligotyping and 24 locus mycobacterial interspersed repetitive units-variable number of tandem repeats were combined to identify clustered isolates. RESULTS: The 173 M. tuberculosis isolates collected belonged predominantly to the Euro-American lineage (Lineage 4) and the East-Asian lineage (Lineage 2). Multidrug-resistant M. tuberculosis were observed in 5.2% of isolates. Lineage 2 M. tuberculosis, which includes the 'Beijing' genotype, was significantly associated with any drug resistance (OR 5.2, 95%CI 1.8-15.1). Cluster analyses showed 44% molecularly clustered isolates, suggesting transmission of M. tuberculosis in the community, including transmission of primary drug-resistant M. tuberculosis. CONCLUSION: These data provide the first insight into the molecular characteristics of M. tuberculosis in the Madang area of PNG, and indicate substantial drug resistance with evidence of ongoing transmission

Mycobacterium tuberculosis harbours little DNA sequence diversity compared with other bacteria. However, there is mounting evidence that strain-to-strain variation in this organism has been underestimated. We review our current understanding of the genetic diversity among M. tuberculosis clinical strains and discuss the relevance of this diversity for the ongoing global epidemics of drug-resistant tuberculosis. Based on findings in other bacteria, we propose that epistatic interactions between pre-existing differences in strain genetic background, acquired drug-resistance-conferring mutations and compensatory changes could play a role in the emergence and spread of drug-resistant M. tuberculosis.

Summary In recent years, the percentage of tuberculosis (TB) cases in immigrants in Spain has increased markedly. In this context, discrimination between cases with potentially imported TB and cases likely to have acquired it by recent transmission after arrival is a basic issue. In this study, we evaluated molecular markers to obtain information on the geographic origin of Mycobacterium tuberculosis (MTB) strains and to determine whether they could help to identify imported TB cases. We analyzed the pks 15/1 gene, spoligotype, and MIRU locus 24 in 168 MTB isolates from immigrants of 32 nationalities and from a selection of autochthonous cases. In non-Asian immigrants we could not detect major differences with the autochthonous cases. For the Asian cases, we found some specific features for these markers. pks 15/1 was intact in 69% of the Asian patients, but this gene had a 7-bp deletion in all non-Asian and Spanish (non-Beijing) cases. The spoligotype-defined lineages EAI and CAS, and the allele with two repetitions in MIRU locus 24 was found exclusively among Asian immigrants. The analysis of certain bacterial molecular markers could help to discriminate between potentially imported TB cases and those more probably acquired in the host country.

BackgroundDifferent combinations of variable number of tandem repeat (VNTR) loci have been proposed for genotyping Mycobacterium tuberculosis complex (MTBC). Existing VNTR schemes show different discriminatory capacity among the six human MTBC lineages. Here, we evaluated the discriminatory power of a "customized MIRU12" loci format proposed previously by Comas et al. based on the standard 24 loci defined by Supply et al. for VNTR-typing of MTBC in Ghana.MethodOne hundred and fifty-eight MTBC isolates classified into Lineage 4 and Lineage 5 were used to compare a customized lineage-specific panel of 12 MIRU-VNTR loci ("customized MIRU-12") to the standard MIRU-15 genotyping scheme. The resolution power of each typing method was determined based on the Hunter-Gaston- Discriminatory Index (HGDI). A minimal set of customized MIRU-VNTR loci for typing Lineages 4 (Euro-American) and 5 (M. africanum West African 1) strains from Ghana was defined based on the cumulative HGDI.Results and conclusionAmong the 106 Lineage 4 strains, the customized MIRU-12 identified a total of 104 distinct genotypes consisting of 2 clusters of 2 isolates each (clustering rate 1.8%), and 102 unique strains while standard MIRU-15 yielded a total of 105 different genotypes, including 1 cluster of 2 isolates (clustering rate: 0.9%) and 104 singletons. Among, 52 Lineage 5 isolates, customized MIRU-12 genotyping defined 51 patterns with 1 cluster of 2 isolates (clustering rate: 0.9%) and 50 unique strains whereas MIRU-15 classified all 52 strains as unique. Cumulative HGDI values for customized MIRU-12 for Lineages 4 and 5 were 0.98 respectively whilst that of standard MIRU-15 was 0.99. A union of loci from the customised MIRU-12 and standard MIRU-15 revealed a set of customized eight highly discriminatory loci: 4052, 2163B, 40, 4165, 2165, 10,16 and 26 with a cumulative HGDI of 0.99 for genotyping Lineage 4 and 5 strains from Ghana.

Tuberculosis caused 20% of all human deaths in the Western world between the seventeenth and nineteenth centuries and remains a cause of high mortality in developing countries. In analogy to other crowd diseases, the origin of human tuberculosis has been associated with the Neolithic Demographic Transition, but recent studies point to a much earlier origin. We analyzed the whole genomes of 259 M. tuberculosis complex (MTBC) strains and used this data set to characterize global diversity and to reconstruct the evolutionary history of this pathogen. Coalescent analyses indicate that MTBC emerged about 70,000 years ago, accompanied migrations of anatomically modern humans out of Africa and expanded as a consequence of increases in human population density during the Neolithic period. This long coevolutionary history is consistent with MTBC displaying characteristics indicative of adaptation to both low and high host densities.

There is increasing evidence that strain variation in Mycobacterium tuberculosis complex (MTBC) might influence the outcome of tuberculosis infection and disease. To assess genotype-phenotype associations, phylogenetically robust molecular markers and appropriate genotyping tools are required. Most current genotyping methods for MTBC are based on mobile or repetitive DNA elements. Because these elements are prone to convergent evolution, the corresponding genotyping techniques are suboptimal for phylogenetic studies and strain classification. By contrast, single nucleotide polymorphisms (SNP) are ideal markers for classifying MTBC into phylogenetic lineages, as they exhibit very low degrees of homoplasy. In this study, we developed two complementary SNP-based genotyping methods to classify strains into the six main human-associated lineages of MTBC, the "Beijing" sublineage, and the clade comprising Mycobacterium bovis and Mycobacterium caprae. Phylogenetically informative SNPs were obtained from 22 MTBC whole-genome sequences. The first assay, referred to as MOL-PCR, is a ligation-dependent PCR with signal detection by fluorescent microspheres and a Luminex flow cytometer, which simultaneously interrogates eight SNPs. The second assay is based on six individual TaqMan real-time PCR assays for singleplex SNP-typing. We compared MOL-PCR and TaqMan results in two panels of clinical MTBC isolates. Both methods agreed fully when assigning 36 well-characterized strains into the main phylogenetic lineages. The sensitivity in allele-calling was 98.6% and 98.8% for MOL-PCR and TaqMan, respectively. Typing of an additional panel of 78 unknown clinical isolates revealed 99.2% and 100% sensitivity in allele-calling, respectively, and 100% agreement in lineage assignment between both methods. While MOL-PCR and TaqMan are both highly sensitive and specific, MOL-PCR is ideal for classification of isolates with no previous information, whereas TaqMan is faster for confirmation. Furthermore, both methods are rapid, flexible and comparably inexpensive.

1 Swiss Tropical and Public Health Institute, Basel, Switzerland, 2 University of Basel, Basel, Switzerland, 3 National Committee for Leprosy, Buruli Ulcer, Yaws and Leishmaniasis Control, Department of Disease Control, Ministry of Public Health, Yaounde, Cameroon, 4 Bankim District Hospital, Bankim, Cameroon, 5 Universitat Heidelberg, Heidelberg, Germany, 6 FAIRMED Africa Regional Office, Yaounde, Cameroon

In this multicenter study, the reliability of two nonradiometric, fully automated systems, the MB/BacT and BACTEC MGIT 960 systems, for testing the susceptibilities of 82 Mycobacterium tuberculosis strains to isoniazid, rifampin, ethambutol, and streptomycin was evaluated in comparison with the radiometric BACTEC 460TB system. The arbitration of discrepant results was done by the reanalysis of the strain, the determination of the MIC, and the molecular characterization of some resistance determinants. The overall level of agreement with BACTEC 460TB results was 96% with the MB/BacT test and 97.2% with the BACTEC MGIT 960 system. With both methods, the level of agreement with BACTEC 460TB results was 96.3% for isoniazid, 98.8% for rifampin, and 98.8% for ethambutol. The level of agreement for streptomycin was 90.2% with MB/BacT and 97.5% with BACTEC MGIT 960. Overall, there were 11 very major errors and 2 major errors with the MB/BacT method and 5 very major errors and 2 major errors with the BACTEC MGIT 960 system. In general, the MB/BacT and BACTEC MGIT 960 systems showed good performance for susceptibility testing with first-line antituberculosis drugs.

Mycobacterium africanum is a member of the Mycobacterium tuberculosis complex (MTBC) and an important cause of human tuberculosis in West Africa that is rarely observed elsewhere. Here we genotyped 613 MTBC clinical isolates from Ghana, and searched for associations between the different phylogenetic lineages of MTBC and patient variables. We found that 17.1% (105/613) of the MTBC isolates belonged to M. africanum, with the remaining belonging to M. tuberculosis sensu stricto. No M. bovis was identified in this sample. M. africanum was significantly more common in tuberculosis patients belonging to the Ewe ethnic group (adjusted odds ratio: 3.02; 95% confidence interval: 1.67–5.47, p, Author Summary Tuberculosis remains one of the main global public health problems. Human tuberculosis is caused by bacteria known as the Mycobacterium tuberculosis complex (MTBC). The MTBC includes a variant called Mycobacterium africanum, which causes up to half of all tuberculosis cases in West Africa. For reasons unknown, M. africanum does not occur in other parts of the world. To explore the possible reasons for this geographic restriction of M. africanum, we analysed a large collection of bacterial strains isolated from tuberculosis patients in Ghana. We genetically characterized these bacterial isolates and collected relevant socio-demographic and epidemiological data. We found tuberculosis patients infected with M. africanum were more likely to belong to the Ewe ethnic group, compared to patients carrying other MTBC bacteria. The Ewes are indigenous inhabitants of coastal regions in West Africa that have previously been shown to harbour a high prevalence of M. africanum. Our findings support the hypothesis that different variants of MTBC have adapted to different human populations, and offer a possible explanation for the geographical restriction of M. africanum to West Africa.

The phylogeographic population structure of Mycobacterium tuberculosis suggests local adaptation to sympatric human populations. We hypothesized that HIV infection, which induces immunodeficiency, will alter the sympatric relationship between M. tuberculosis and its human host. To test this hypothesis, we performed a nine-year nation-wide molecular-epidemiological study of HIV–infected and HIV–negative patients with tuberculosis (TB) between 2000 and 2008 in Switzerland. We analyzed 518 TB patients of whom 112 (21.6%) were HIV–infected and 233 (45.0%) were born in Europe. We found that among European-born TB patients, recent transmission was more likely to occur in sympatric compared to allopatric host–pathogen combinations (adjusted odds ratio [OR] 7.5, 95% confidence interval [95% CI] 1.21–infinity, p = 0.03). HIV infection was significantly associated with TB caused by an allopatric (as opposed to sympatric) M. tuberculosis lineage (OR 7.0, 95% CI 2.5–19.1, p, Author Summary Human tuberculosis (TB) caused by Mycobacterium tuberculosis kills 1.5 million people each year. M. tuberculosis has been affecting humans for millennia, suggesting that different strain lineages may be adapted to specific human populations. The combination of a particular strain lineage and its corresponding patient population can be classified as sympatric (e.g. Euro-American lineage in Europeans) or allopatric (e.g. East-Asian lineage in Europeans). We hypothesized that infection with the human immunodeficiency virus (HIV), which impairs the human immune system, will interfere with this host–pathogen relationship. We performed a nation-wide molecular-epidemiological study of HIV–infected and HIV–negative TB patients between 2000 and 2008 in Switzerland. We found that HIV infection was associated with the less adapted allopatric lineages among patients born in Europe, and this was not explained by social or other patient factors such as increased social mixing in HIV–infected individuals. Strikingly, the association between HIV infection and less adapted M. tuberculosis lineages was stronger in patients with more pronounced immunodeficiency. Our observation was replicated in a second independent panel of M. tuberculosis strains collected during a population-based study in the Canton of Bern. In summary, our study provides evidence that the sympatric host–pathogen relationship in TB is disrupted by HIV infection.

Background Tuberculosis (TB) is a major public health problem in Nepal. Strain variation in Mycobacterium tuberculosis may influence the outcome of TB infection and disease. To date, the phylogenetic diversity of M. tuberculosis in Nepal is unknown. Methods and Findings We analyzed 261 M. tuberculosis isolates recovered from pulmonary TB patients recruited between August 2009 and August 2010 in Nepal. M. tuberculosis lineages were determined by single nucleotide polymorphisms (SNP) typing and spoligotyping. Drug resistance was determined by sequencing the hot spot regions of the relevant target genes. Overall, 164 (62.8%) TB patients were new, and 97 (37.2%) were previously treated. Any drug resistance was detected in 50 (19.2%) isolates, and 16 (6.1%) were multidrug-resistant. The most frequent M. tuberculosis lineage was Lineage 3 (CAS/Delhi) with 106 isolates (40.6%), followed by Lineage 2 (East-Asian lineage, includes Beijing genotype) with 84 isolates (32.2%), Lineage 4 (Euro-American lineage) with 41 (15.7%) isolates, and Lineage 1 (Indo-Oceanic lineage) with 30 isolates (11.5%). Based on spoligotyping, we found 45 different spoligotyping patterns that were previously described. The Beijing (83 isolates, 31.8%) and CAS spoligotype (52, 19.9%) were the dominant spoligotypes. A total of 36 (13.8%) isolates could not be assigned to any known spoligotyping pattern. Lineage 2 was associated with female sex (adjusted odds ratio [aOR] 2.58, 95% confidence interval [95% CI] 1.42–4.67, p = 0.002), and any drug resistance (aOR 2.79; 95% CI 1.43–5.45; p = 0.002). We found no evidence for an association of Lineage 2 with age or BCG vaccination status. Conclusions We found a large genetic diversity of M. tuberculosis in Nepal with representation of all four major lineages. Lineages 3 and 2 were dominating. Lineage 2 was associated with clinical characteristics. This study fills an important gap on the map of the M. tuberculosis genetic diversity in the Asian region.

Drug-resistant tuberculosis is a growing threat to global public health. Recent efforts to understand the evolution of drug resistance have shown that changes in drug–target interactions are only the first step in a longer adaptive process. The emergence of transmissible drug-resistant Mycobacterium tuberculosis is the result of a multitude of additional genetic mutations, many of which interact, a phenomenon known as epistasis. The varied effects of these epistatic interactions include compensating for the reduction of the biological cost associated with the development of drug resistance, increasing the level of resistance, and possibly accommodating broader changes in the physiology of resistant bacteria. Knowledge of these processes and our ability to detect them as they happen informs the development of diagnostic tools and better control strategies. In particular, the use of whole genome sequencing combined with surveillance efforts in the field could provide a powerful instrument to prevent future epidemics of drug-resistant tuberculosis.

Multidrug-resistant tuberculosis (MDR-TB) accounts for one third of the annual deaths due to antimicrobial resistance(1). Drug resistance-conferring mutations frequently cause fitness costs in bacteria(2–5). Experimental work indicates that these drug resistance-related fitness costs might be mitigated by compensatory mutations(6–10). However, the clinical relevance of compensatory evolution remains poorly understood. Here we show that, in the country of Georgia, during a 6-year nationwide study, 63% of MDR-TB was due to patient-to-patient transmission. Compensatory mutations and patient incarceration were independently associated with transmission. Furthermore, compensatory mutations were overrepresented among isolates from incarcerated individuals that also frequently spilled over into the non-incarcerated population. As a result, up to 31% of MDR-TB in Georgia was directly or indirectly linked to prisons. We conclude that prisons fuel the epidemic of MDR-TB in Georgia by acting as ecological drivers of fitness-compensated strains with high transmission potential.

Summary Background Drug resistance threatens global tuberculosis control. We aimed to examine mortality in patients with tuberculosis from high-burden countries, according to concordance or discordance of results from drug susceptibility testing done locally and whole-genome sequencing (WGS). Methods In this multicentre cohort study, we collected pulmonary Mycobacterium tuberculosis isolates and clinical data from individuals with tuberculosis from antiretroviral therapy programmes and tuberculosis clinics in Cote d'Ivoire, Democratic Republic of the Congo, Kenya, Nigeria, Peru, South Africa, and Thailand, stratified by HIV status and drug resistance. Sites tested drug susceptibility using routinely available methods. WGS was done on Illumina HiSeq 2500 in the USA and Switzerland, and TBprofiler was used to analyse the genomes. We included individuals aged 16 years or older with pulmonary tuberculosis (bacteriologically confirmed or clinically diagnosed). We analysed mortality in multivariable logistic regression models adjusted for sex, age, HIV status, history of tuberculosis, and sputum positivity. Findings Between Sept 1, 2014, and July 4, 2016, of 634 patients included in our previous analysis, we included 582 patients with tuberculosis (median age 33 years [IQR 27–43], 225 [39%] women, and 247 [42%] HIV-positive). Based on WGS, 339 (58%) isolates were pan-susceptible, 35 (6%) monoresistant, 146 (25%) multidrug-resistant, and 24 (4%) pre-extensively drug-resistant (pre-XDR) or XDR. The analysis of mortality was based on 530 patients; 63 (12%) died and 77 (15%) patients received inappropriate treatment. Mortality ranged from 6% (18 of 310) in patients with pan-susceptible tuberculosis to 39% (nine of 23) in patients with pre-XDR or XDR tuberculosis. The adjusted odds ratio for mortality was 4·92 (95% CI 2·47–9·78) among undertreated patients, compared with appropriately treated patients. Interpretation In seven countries with a high burden of tuberculosis, we observed discrepancies between drug resistance patterns obtained locally and WGS. The underdiagnosis of drug resistance resulted in inappropriate treatment and higher mortality. WGS can provide accurate and detailed drug resistance information required to improve the outcomes of drug-resistant tuberculosis in high-burden settings. Our results support WHO's call for point-of-care tests based on WGS. Funding National Institutes of Allergy and Infectious Diseases, Swiss National Science Foundation, and Swiss National Center for Mycobacteria.