Two billion people are infected with Mycobacterium tuberculosis, leading to 10 million new cases of active tuberculosis and 1.5 million deaths annually. Universal access to drug susceptibility testing (DST) has become a World Health Organization priority. We previously developed a software tool, Mykrobe predictor, which provided offline species identification and drug resistance predictions for M. tuberculosis from whole genome sequencing (WGS) data. Performance was insufficient to support the use of WGS as an alternative to conventional phenotype-based DST, due to mutation catalogue limitations. Here we present a new tool, Mykrobe, which provides the same functionality based on a new software implementation. Improvements include i) an updated mutation catalogue giving greater sensitivity to detect pyrazinamide resistance, ii) support for user-defined resistance catalogues, iii) improved identification of non-tuberculous mycobacterial species, and iv) an updated statistical model for Oxford Nanopore Technologies sequencing data. Mykrobe is released under MIT license at https://github.com/mykrobe-tools/mykrobe. We incorporate mutation catalogues from the CRyPTIC consortium et al. (2018) and from Walker et al. (2015), and make improvements based on performance on an initial set of 3206 and an independent set of 5845 M. tuberculosis Illumina sequences. To give estimates of error rates, we use a prospectively collected dataset of 4362 M. tuberculosis isolates. Using culture based DST as the reference, we estimate Mykrobe to be 100%, 95%, 82%, 99% sensitive and 99%, 100%, 99%, 99% specific for rifampicin, isoniazid, pyrazinamide and ethambutol resistance prediction respectively. We benchmark against four other tools on 10207 (=5845+4362) samples, and also show that Mykrobe gives concordant results with nanopore data. We measure the ability of Mykrobe-based DST to guide personalized therapeutic regimen design in the context of complex drug susceptibility profiles, showing 94% concordance of implied regimen with that driven by phenotypic DST, higher than all other benchmarked tools.

Pyrazinamide and fluoroquinolones are essential antituberculosis drugs in new rifampicin-sparing regimens. However, little information about the extent of resistance to these drugs at the population level is available. In a molecular epidemiology analysis, we used population-based surveys from Azerbaijan, Bangladesh, Belarus, Pakistan, and South Africa to investigate resistance to pyrazinamide and fluoroquinolones among patients with tuberculosis. Resistance to pyrazinamide was assessed by gene sequencing with the detection of resistance-conferring mutations in the pncA gene, and susceptibility testing to fluoroquinolones was conducted using the MGIT system. Pyrazinamide resistance was assessed in 4972 patients. Levels of resistance varied substantially in the surveyed settings (3·0-42·1%). In all settings, pyrazinamide resistance was significantly associated with rifampicin resistance. Among 5015 patients who underwent susceptibility testing to fluoroquinolones, proportions of resistance ranged from 1·0-16·6% for ofloxacin, to 0·5-12·4% for levofloxacin, and 0·9-14·6% for moxifloxacin when tested at 0·5 μg/mL. High levels of ofloxacin resistance were detected in Pakistan. Resistance to moxifloxacin and gatifloxacin when tested at 2 μg/mL was low in all countries. Although pyrazinamide resistance was significantly associated with rifampicin resistance, this drug may still be effective in 19-63% of patients with rifampicin-resistant tuberculosis. Even though the high level of resistance to ofloxacin found in Pakistan is worrisome because it might be the expression of extensive and unregulated use of fluoroquinolones in some parts of Asia, the negligible levels of resistance to fourth-generation fluoroquinolones documented in all survey sites is an encouraging finding. Rational use of this class of antibiotics should therefore be ensured to preserve its effectiveness. Bill & Melinda Gates Foundation, United States Agency for International Development, Global Alliance for Tuberculosis Drug Development.

Mycobacterium tuberculosis is a clonal pathogen proposed to have co-evolved with its human host for millennia, yet our understanding of its genomic diversity and biogeography remains incomplete. Here we use a combination of phylogenetics and dimensionality reduction to reevaluate the population structure of M. tuberculosis, providing an in-depth analysis of the ancient Indo-Oceanic Lineage 1 and the modern Central Asian Lineage 3, and expanding our understanding of Lineages 2 and 4. We assess sub-lineages using genomic sequences from 4939 pan-susceptible strains, and find 30 new genetically distinct clades that we validate in a dataset of 4645 independent isolates. We find a consistent geographically restricted or unrestricted pattern for 20 groups, including three groups of Lineage 1. The distribution of terminal branch lengths across the M. tuberculosis phylogeny supports the hypothesis of a higher transmissibility of Lineages 2 and 4, in comparison with Lineages 3 and 1, on a global scale. We define an expanded barcode of 95 single nucleotide substitutions that allows rapid identification of 69 M. tuberculosis sub-lineages and 26 additional internal groups. Our results paint a higher resolution picture of the M. tuberculosis phylogeny and biogeography.

The surveillance of drug resistance among tuberculosis (TB) patients is central to combatting the global TB epidemic and preventing the spread of antimicrobial resistance. Isoniazid and rifampicin are two of the most powerful first-line anti-TB medicines, and resistance to either of them increases the risk of treatment failure, relapse, or acquisition of resistance to other drugs. The global prevalence of rifampicin resistance is well documented, occurring in 3.4% (95% CI 2.5%–4.4%) of new TB patients and 18% (95% CI 7.6%–31%) of previously treated TB patients in 2018, whereas the prevalence of isoniazid resistance at global and regional levels is less understood. In 2018, the World Health Organization (WHO) recommended a modified 6-month treatment regimen for people with isoniazid-resistant, rifampicin-susceptible TB (Hr-TB), which includes rifampicin, pyrazinamide, ethambutol, and levofloxacin. We estimated the global prevalence of Hr-TB among TB patients and investigated associated phenotypic and genotypic drug resistance patterns. Aggregated drug resistance data reported to WHO from either routine continuous surveillance or nationally representative periodic surveys of TB patients for the period 2003–2017 were reviewed. Isoniazid data were available from 156 countries or territories for 211,753 patients. Among these, the global prevalence of Hr-TB was 7.4% (95% CI 6.5%–8.4%) among new TB patients and 11.4% (95% CI 9.4%–13.4%) among previously treated TB patients. Additional data on pyrazinamide and levofloxacin resistance were available from 6 countries (Azerbaijan, Bangladesh, Belarus, Pakistan, the Philippines, and South Africa). There were no cases of resistance to both pyrazinamide and levofloxacin among Hr-TB patients, except for the Philippines (1.8%, 95% CI 0.2–6.4) and Belarus (5.3%, 95% CI 0.1–26.0). Sequencing data for all genomic regions involved in isoniazid resistance were available for 4,563 patients. Among the 1,174 isolates that were resistant by either phenotypic testing or sequencing, 78.6% (95% CI 76.1%–80.9%) had resistance-conferring mutations in the katG gene and 14.6% (95% CI 12.7%–16.8%) in both katG and the inhA promoter region. For 6.8% (95% CI 5.4%–8.4%) of patients, mutations occurred in the inhA promoter alone, for whom an increased dose of isoniazid may be considered. The main limitations of this study are that most analyses were performed at the national rather than individual patient level and that the quality of laboratory testing may vary between countries. In this study, the prevalence of Hr-TB among TB patients was higher than the prevalence of rifampicin resistance globally. Many patients with Hr-TB would be missed by current diagnostic algorithms driven by rifampicin testing, highlighting the need for new rapid molecular technologies to ensure access to appropriate treatment and care. The low prevalence of resistance to pyrazinamide and fluoroquinolones among patients with Hr-TB provides further justification for the recommended modified treatment regimen.

In many countries, regular monitoring of the emergence of resistance to anti-tuberculosis drugs is hampered by the limitations of phenotypic testing for drug susceptibility. We therefore evaluated the use of genetic sequencing for surveillance of drug resistance in tuberculosis. Population-level surveys were done in hospitals and clinics in seven countries (Azerbaijan, Bangladesh, Belarus, Pakistan, Philippines, South Africa, and Ukraine) to evaluate the use of genetic sequencing to estimate the resistance of Mycobacterium tuberculosis isolates to rifampicin, isoniazid, ofloxacin, moxifloxacin, pyrazinamide, kanamycin, amikacin, and capreomycin. For each drug, we assessed the accuracy of genetic sequencing by a comparison of the adjusted prevalence of resistance, measured by genetic sequencing, with the true prevalence of resistance, determined by phenotypic testing. Isolates were taken from 7094 patients with tuberculosis who were enrolled in the study between November, 2009, and May, 2014. In all tuberculosis cases, the overall pooled sensitivity values for predicting resistance by genetic sequencing were 91% (95% CI 87-94) for rpoB (rifampicin resistance), 86% (74-93) for katG, inhA, and fabG promoter combined (isoniazid resistance), 54% (39-68) for pncA (pyrazinamide resistance), 85% (77-91) for gyrA and gyrB combined (ofloxacin resistance), and 88% (81-92) for gyrA and gyrB combined (moxifloxacin resistance). For nearly all drugs and in most settings, there was a large overlap in the estimated prevalence of drug resistance by genetic sequencing and the estimated prevalence by phenotypic testing. Genetic sequencing can be a valuable tool for surveillance of drug resistance, providing new opportunities to monitor drug resistance in tuberculosis in resource-poor countries. Before its widespread adoption for surveillance purposes, there is a need to standardise DNA extraction methods, recording and reporting nomenclature, and data interpretation. Bill & Melinda Gates Foundation, United States Agency for International Development, Global Alliance for Tuberculosis Drug Development.

Pakistan is fifth among high burden countries for tuberculosis. A steady increase is seen in extrapulmonary tuberculosis (EPTB), which now accounts for 20% of all notified TB cases. There is very limited information on the epidemiology of EPTB. This study was performed with the aim to describe the demographic characteristics, clinical manifestations and treatment outcomes of EPTB patients in Pakistan. We performed descriptive analysis on routinely collected data for cohorts of TB patients registered nationwide in 2016 at health facilities selected using stratified convenient sampling. Altogether 54092 TB including 15790 (29.2%) EPTB cases were registered in 2016 at 50 study sites. The median age was 24 years for EPTB and 30 years for PTB patients. The crude prevalence of EPTB in females was 30.5% (95%CI; 29.9–31.0) compared to 27.9% (95%CI; 27.3–28.4) in males. The likelihood of having EPTB (OR), was 1.1 times greater for females, 2.0 times for children, and 3.3 times for residents of provinces in the North-West. The most common forms of EPTB were pleural (29.6%), lymphatic (22.7%) and abdominal TB (21.0%). Pleural TB was the most common clinical manifestation in adults (34.2%) and abdominal TB in children (38.4%). An increase in the prevalence of pleural and osteoarticular and decline in lymphatic and abdominal TB was observed with advancing age. Diversity in demography and clinical manifestations were noted between provinces. The treatment success rate for all type EPTB was significantly high compared to bacteriology confirmed PTB with the exception of EPTB affecting CNS with a high mortality rate. The study provides an insight into demography, clinical manifestations and treatment outcomes of EPTB. Further studies are needed to explain significant diversities observed between provinces, specific risk factors and challenges concerning EPTB management.