Emergence of multidrug-resistant Mycobacterium tuberculosis of the Beijing lineage in Portugal and Guinea-Bissau: a snapshot of moving clones by whole-genome sequencing

The Beijing genotype comprises a highly disseminated strain type that is frequently associated with multidrug resistant (MDR) tuberculosis (TB) and increased transmissibility but, countries such as Portugal and Guinea-Bissau fall outside the regions phylogeographically associated with this specific genotype. Nevertheless, recent data shows that this genotype might be gradually emerging in these two countries as an underlying cause of primary MDR-TB. Here, we describe the emergence of Mycobacterium tuberculosis Beijing strains associated with MDR-TB in Portugal and Guinea-Bissau demonstrating the presence of the well described superclusters 100-32 and 94-32 in Portugal and Guinea-Bissau, respectively. Genome-wide analysis and comparison with a global genomic dataset of M. tuberculosis Beijing strains, revealed the presence of two genomic clusters encompassing isolates from Portugal and Guinea-Bissau, GC1 (n = 121) and GC2 (n = 39), both of which bore SNP signatures compatible with the 100-32/B0/W148 and 94-32/Central Asia Outbreak clades, respectively. Moreover, GC2 encompasses a cross-border cluster between Portugal, Guinea-Bissau and Brazil thus supporting migration-associated introduction of MDR-TB and subsequent clonal expansion at the community-level. The comparison with global Beijing datasets demonstrates the global reach of the disease and its complex dissemination across multiple countries while in parallel there are clear microevolutionary trajectories towards extensively drug resistant TB.

Machine Learning Predicts Accurately Mycobacterium tuberculosis Drug Resistance From Whole Genome Sequencing Data

Background: Tuberculosis disease, caused by Mycobacterium tuberculosis, is a major public health problem. The emergence of M. tuberculosis strains resistant to existing treatments threatens to derail control efforts. Resistance is mainly conferred by mutations in genes coding for drug targets or converting enzymes, but our knowledge of these mutations is incomplete. Whole genome sequencing (WGS) is an increasingly common approach to rapidly characterize isolates and identify mutations predicting antimicrobial resistance and thereby providing a diagnostic tool to assist clinical decision making. Methods: We applied machine learning approaches to 16,688 M. tuberculosis isolates that have undergone WGS and laboratory drug-susceptibility testing (DST) across 14 antituberculosis drugs, with 22.5% of samples being multidrug resistant and 2.1% being extensively drug resistant. We used non-parametric classification-tree and gradient-boosted-tree models to predict drug resistance and uncover any associated novel putative mutations. We fitted separate models for each drug, with and without "co-occurrent resistance" markers known to be causing resistance to drugs other than the one of interest. Predictive performance was measured using sensitivity, specificity, and the area under the receiver operating characteristic curve, assuming DST results as the gold standard. Results: The predictive performance was highest for resistance to first-line drugs, amikacin, kanamycin, ciprofloxacin, moxifloxacin, and multidrug-resistant tuberculosis (area under the receiver operating characteristic curve above 96%), and lowest for third-line drugs such as D-cycloserine and Para-aminosalisylic acid (area under the curve below 85%). The inclusion of co-occurrent resistance markers led to improved performance for some drugs and superior results when compared to similar models in other large-scale studies, which had smaller sample sizes. Overall, the gradient-boosted-tree models performed better than the classification-tree models. The mutation-rank analysis detected no new single nucleotide polymorphisms linked to drug resistance. Discordance between DST and genotypically inferred resistance may be explained by DST errors, novel rare mutations, hetero-resistance, and nongenomic drivers such as efflux-pump upregulation. Conclusion: Our work demonstrates the utility of machine learning as a flexible approach to drug resistance prediction that is able to accommodate a much larger number of predictors and to summarize their predictive ability, thus assisting clinical decision making and single nucleotide polymorphism detection in an era of increasing WGS data generation.

Spatial Network Mapping of Pulmonary Multidrug-Resistant Tuberculosis Cavities Using RNA Sequencing

Rationale: There is poor understanding about protective immunity and the pathogenesis of cavitation in patients with tuberculosis.Objectives: To map pathophysiological pathways at anatomically distinct positions within the human tuberculosis cavity.Methods: Biopsies were obtained from eight predetermined locations within lung cavities of patients with multidrug-resistant tuberculosis undergoing therapeutic surgical resection (n = 14) and healthy lung tissue from control subjects without tuberculosis (n = 10). RNA sequencing, immunohistochemistry, and bacterial load determination were performed at each cavity position. Differentially expressed genes were normalized to control subjects without tuberculosis, and ontologically mapped to identify a spatially compartmentalized pathophysiological map of the cavity. In silico perturbation using a novel distance-dependent dynamical sink model was used to investigate interactions between immune networks and bacterial burden, and to integrate these identified pathways.Measurements and Main Results: The median (range) lung cavity volume on positron emission tomography/computed tomography scans was 50 cm3 (15-389 cm3). RNA sequence reads (31% splice variants) mapped to 19,049 annotated human genes. Multiple proinflammatory pathways were upregulated in the cavity wall, whereas a downregulation "sink" in the central caseum-fluid interface characterized 53% of pathways including neuroendocrine signaling, calcium signaling, triggering receptor expressed on myeloid cells-1, reactive oxygen and nitrogen species production, retinoic acid-mediated apoptosis, and RIG-I-like receptor signaling. The mathematical model demonstrated that neuroendocrine, protein kinase C-θ, and triggering receptor expressed on myeloid cells-1 pathways, and macrophage and neutrophil numbers, had the highest correlation with bacterial burden (r > 0.6), whereas T-helper effector systems did not.Conclusions: These data provide novel insights into host immunity to Mycobacterium tuberculosis-related cavitation. The pathways defined may serve as useful targets for the design of host-directed therapies, and transmission prevention interventions.

Integrating informatics tools and portable sequencing technology for rapid detection of resistance to anti-tuberculous drugs

Background

Mycobacterium tuberculosis resistance to anti-tuberculosis drugs is a major threat to global public health. Whole genome sequencing (WGS) is rapidly gaining traction as a diagnostic tool for clinical tuberculosis settings. To support this informatically, previous work led to the development of the widely used TBProfiler webtool, which predicts resistance to 14 drugs from WGS data. However, for accurate and rapid high throughput of samples in clinical or epidemiological settings, there is a need for a stand-alone tool and the ability to analyse data across multiple WGS platforms, including Oxford Nanopore MinION.

Results

We present a new command line version of the TBProfiler webserver, which includes hetero-resistance calling and will facilitate the batch processing of samples. The TBProfiler database has been expanded to incorporate 178 new markers across 16 anti-tuberculosis drugs. The predictive performance of the mutation library has been assessed using > 17,000 clinical isolates with WGS and laboratory-based drug susceptibility testing (DST) data. An integrated MinION analysis pipeline was assessed by performing WGS on 34 replicates across 3 multi-drug resistant isolates with known resistance mutations. TBProfiler accuracy varied by individual drug. Assuming DST as the gold standard, sensitivities for detecting multi-drug-resistant TB (MDR-TB) and extensively drug-resistant TB (XDR-TB) were 94% (95%CI 93–95%) and 83% (95%CI 79–87%) with specificities of 98% (95%CI 98–99%) and 96% (95%CI 95–97%) respectively. Using MinION data, only one resistance mutation was missed by TBProfiler, involving an insertion in the tlyA gene coding for capreomycin resistance. When compared to alternative platforms (e.g. Mykrobe predictor TB, the CRyPTIC library), TBProfiler demonstrated superior predictive performance across first- and second-line drugs.

Conclusions

The new version of TBProfiler can rapidly and accurately predict anti-TB drug resistance profiles across large numbers of samples with WGS data. The computing architecture allows for the ability to modify the core bioinformatic pipelines and outputs, including the analysis of WGS data sourced from portable technologies. TBProfiler has the potential to be integrated into the point of care and WGS diagnostic environments, including in resource-poor settings.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0650-x) contains supplementary material, which is available to authorized users.

Integrating standardized whole genome sequence analysis with a global Mycobacterium tuberculosis antibiotic resistance knowledgebase

Drug-resistant tuberculosis poses a persistent public health threat. The ReSeqTB platform is a collaborative, curated knowledgebase, designed to standardize and aggregate global Mycobacterium tuberculosis complex (MTBC) variant data from whole genome sequencing (WGS) with phenotypic drug susceptibility testing (DST) and clinical data. We developed a unified analysis variant pipeline (UVP) ( https://github.com/CPTR-ReSeqTB/UVP ) to identify variants and assign lineage from MTBC sequence data. Stringent thresholds and quality control measures were incorporated in this open source tool. The pipeline was validated using a well-characterized dataset of 90 diverse MTBC isolates with conventional DST and DNA Sanger sequencing data. The UVP exhibited 98.9% agreement with the variants identified using Sanger sequencing and was 100% concordant with conventional methods of assigning lineage. We analyzed 4636 publicly available MTBC isolates in the ReSeqTB platform representing all seven major MTBC lineages. The variants detected have an above 94% accuracy of predicting drug based on the accompanying DST results in the platform. The aggregation of variants over time in the platform will establish confidence-graded mutations statistically associated with phenotypic drug resistance. These tools serve as critical reference standards for future molecular diagnostic assay developers, researchers, public health agencies and clinicians working towards the control of drug-resistant tuberculosis.

Identifying mixed Mycobacterium tuberculosis infections from whole genome sequence data

BACKGROUND

Mixed, polyclonal Mycobacterium tuberculosis infection occurs in natural populations. Developing an effective method for detecting such cases is important in measuring the success of treatment and reconstruction of transmission between patients. Using whole genome sequence (WGS) data, we assess two methods for detecting mixed infection: (i) a combination of the number of heterozygous sites and the proportion of heterozygous sites to total SNPs, and (ii) Bayesian model-based clustering of allele frequencies from sequencing reads at heterozygous sites.

RESULTS

In silico and in vitro artificially mixed and known pure M. tuberculosis samples were analysed to determine the specificity and sensitivity of each method. We found that both approaches were effective in distinguishing between pure strains and mixed infection where there was relatively high (> 10%) proportion of a minor strain in the mixture. A large dataset of clinical isolates (n = 1963) from the Karonga Prevention Study in Northern Malawi was tested to examine correlations with patient characteristics and outcomes with mixed infection. The frequency of mixed infection in the population was found to be around 10%, with an association with year of diagnosis, but no association with age, sex, HIV status or previous tuberculosis.

CONCLUSIONS

Mixed Mycobacterium tuberculosis infection was identified in silico using whole genome sequence data. The methods presented here can be applied to population-wide analyses of tuberculosis to estimate the frequency of mixed infection, and to identify individual cases of mixed infections. These cases are important when considering the evolution and transmission of the disease, and in patient treatment.

Use of whole genome sequencing in surveillance of drug resistant tuberculosis

INTRODUCTION

The threat of resistance to anti-tuberculosis drugs is of global concern. Current efforts to monitor resistance rely on phenotypic testing where cultured bacteria are exposed to critical concentrations of the drugs. Capacity for such testing is low in TB endemic countries. Drug resistance is caused by mutations in the Mycobacterium tuberculosis genome and whole genome sequencing to detect these mutations offers an alternative means of assessing resistance. Areas covered: The challenges of assessing TB drug resistance are discussed. Progress in elucidating the M. tuberculosis resistome and evidence of the accuracy of next generation sequencing for detecting resistance is reviewed. Expert Commentary: There are considerable advantages to using next generation sequencing for TB drug resistance surveillance. Accuracy is high for detecting resistance to the major first-line drugs but is currently lower for the second-line drugs due to our incomplete knowledge regarding resistance causing mutations. With the advances in sequencing technology and the opportunity to replace phenotypic drug susceptibility testing with safer and more cost effective methods it would appear that the question is when to implement. Current bottlenecks are sample extraction to allow whole genome sequencing directly from sputum and the lack of bioinformatics expertise in some TB endemic countries.

Author Correction: Genome-wide analysis of multi- and extensively drug-resistant Mycobacterium tuberculosis

In the version of this article initially published, the URL listed for TubercuList was incorrect. The correct URL is https://mycobrowser.epfl.ch/. The error has been corrected in the HTML and PDF versions of the article.

Clonal expansion across the seas as seen through CPLP-TB database: A joint effort in cataloguing Mycobacterium tuberculosis genetic diversity in Portuguese-speaking countries

Tuberculosis (TB) remains a major health problem within the Community of Portuguese Language Speaking Countries (CPLP). Despite the marked variation in TB incidence across its member-states and continued human migratory flux between countries, a considerable gap in the knowledge on the Mycobacterium tuberculosis population structure and strain circulation between the countries still exists. To address this, we have assembled and analysed the largest CPLP M. tuberculosis molecular and drug susceptibility dataset, comprised by a total of 1447 clinical isolates, including 423 multidrug-resistant isolates, from five CPLP countries. The data herein presented reinforces Latin American and Mediterranean (LAM) strains as the hallmark of M. tuberculosis populational structure in the CPLP coupled with country-specific differential prevalence of minor clades. Moreover, using high-resolution typing by 24-loci MIRU-VNTR, six cross-border genetic clusters were detected, thus supporting recent clonal expansion across the Lusophone space.

To make this data available to the scientific community and public health authorities we developed CPLP-TB (available at http://cplp-tb.ff.ulisboa.pt), an online database coupled with web-based tools for exploratory data analysis. As a public health tool, it is expected to contribute to improved knowledge on the M. tuberculosis population structure and strain circulation within the CPLP, thus supporting the risk assessment of strain-specific trends.

•

The Community of Portuguese Speaking Countries (CPLP) occupies a vast geographical area.

•

Three CPLP countries are shortlisted in the WHO's list of Top 30 high-burden countries.

•

Common Mycobacterium tuberculosis population structure denote historical strain flow.

•

Cross-border clusters suggest recent intercontinental tuberculosis transmission.

•

CPLP-TB: a novel strain database and framework for collaborative studies and strain tracing.

Genome-wide analysis of multi- and extensively drug-resistant Mycobacterium tuberculosis

To characterize the genetic determinants of resistance to antituberculosis drugs, we performed a genome-wide association study (GWAS) of 6,465 Mycobacterium tuberculosis clinical isolates from more than 30 countries. A GWAS approach within a mixed-regression framework was followed by a phylogenetics-based test for independent mutations. In addition to mutations in established and recently described resistance-associated genes, novel mutations were discovered for resistance to cycloserine, ethionamide and para-aminosalicylic acid. The capacity to detect mutations associated with resistance to ethionamide, pyrazinamide, capreomycin, cycloserine and para-aminosalicylic acid was enhanced by inclusion of insertions and deletions. Odds ratios for mutations within candidate genes were found to reflect levels of resistance. New epistatic relationships between candidate drug-resistance-associated genes were identified. Findings also suggest the involvement of efflux pumps (drrA and Rv2688c) in the emergence of resistance. This study will inform the design of new diagnostic tests and expedite the investigation of resistance and compensatory epistatic mechanisms.

Methylation in Mycobacterium tuberculosis is lineage specific with associated mutations present globally

DNA methylation is an epigenetic modification of the genome involved in regulating crucial cellular processes, including transcription and chromosome stability. Advances in PacBio sequencing technologies can be used to robustly reveal methylation sites. The methylome of the Mycobacterium tuberculosis complex is poorly understood but may be involved in virulence, hypoxic survival and the emergence of drug resistance. In the most extensive study to date, we characterise the methylome across the 4 major lineages of M. tuberculosis and 2 lineages of M. africanum, the leading causes of tuberculosis disease in humans. We reveal lineage-specific methylated motifs and strain-specific mutations that are abundant globally and likely to explain loss of function in the respective methyltransferases. Our work provides a set of sixteen new complete reference genomes for the Mycobacterium tuberculosis complex, including complete lineage 5 genomes. Insights into lineage-specific methylomes will further elucidate underlying biological mechanisms and other important phenotypes of the epi-genome.

A standardised method for interpreting the association between mutations and phenotypic drug resistance in Mycobacterium tuberculosis

A clear understanding of the genetic basis of antibiotic resistance in Mycobacterium tuberculosis is required to accelerate the development of rapid drug susceptibility testing methods based on genetic sequence.Raw genotype-phenotype correlation data were extracted as part of a comprehensive systematic review to develop a standardised analytical approach for interpreting resistance associated mutations for rifampicin, isoniazid, ofloxacin/levofloxacin, moxifloxacin, amikacin, kanamycin, capreomycin, streptomycin, ethionamide/prothionamide and pyrazinamide. Mutation frequencies in resistant and susceptible isolates were calculated, together with novel statistical measures to classify mutations as high, moderate, minimal or indeterminate confidence for predicting resistance.We identified 286 confidence-graded mutations associated with resistance. Compared to phenotypic methods, sensitivity (95% CI) for rifampicin was 90.3% (89.6-90.9%), while for isoniazid it was 78.2% (77.4-79.0%) and their specificities were 96.3% (95.7-96.8%) and 94.4% (93.1-95.5%), respectively. For second-line drugs, sensitivity varied from 67.4% (64.1-70.6%) for capreomycin to 88.2% (85.1-90.9%) for moxifloxacin, with specificity ranging from 90.0% (87.1-92.5%) for moxifloxacin to 99.5% (99.0-99.8%) for amikacin.This study provides a standardised and comprehensive approach for the interpretation of mutations as predictors of M. tuberculosis drug-resistant phenotypes. These data have implications for the clinical interpretation of molecular diagnostics and next-generation sequencing as well as efficient individualised therapy for patients with drug-resistant tuberculosis.

A multiple genome analysis of Mycobacterium tuberculosis reveals specific novel genes and mutations associated with pyrazinamide resistance

BACKGROUND

Tuberculosis (TB) is a major global health problem and drug resistance compromises the efforts to control this disease. Pyrazinamide (PZA) is an important drug used in both first and second line treatment regimes. However, its complete mechanism of action and resistance remains unclear.

RESULTS

We genotyped and sequenced the complete genomes of 68 M. tuberculosis strains isolated from unrelated TB patients in Peru. No clustering pattern of the strains was verified based on spoligotyping. We analyzed the association between PZA resistance with non-synonymous mutations and specific genes. We found mutations in pncA and novel genes significantly associated with PZA resistance in strains without pncA mutations. These included genes related to transportation of metal ions, pH regulation and immune system evasion.

CONCLUSIONS

These results suggest potential alternate mechanisms of PZA resistance that have not been found in other populations, supporting that the antibacterial activity of PZA may hit multiple targets.

The epidemiology, pathogenesis, transmission, diagnosis, and management of multidrug-resistant, extensively drug-resistant, and incurable tuberculosis

Global tuberculosis incidence has declined marginally over the past decade, and tuberculosis remains out of control in several parts of the world including Africa and Asia. Although tuberculosis control has been effective in some regions of the world, these gains are threatened by the increasing burden of multidrug-resistant (MDR) and extensively drug-resistant (XDR) tuberculosis. XDR tuberculosis has evolved in several tuberculosis-endemic countries to drug-incurable or programmatically incurable tuberculosis (totally drug-resistant tuberculosis). This poses several challenges similar to those encountered in the pre-chemotherapy era, including the inability to cure tuberculosis, high mortality, and the need for alternative methods to prevent disease transmission. This phenomenon mirrors the worldwide increase in antimicrobial resistance and the emergence of other MDR pathogens, such as malaria, HIV, and Gram-negative bacteria. MDR and XDR tuberculosis are associated with high morbidity and substantial mortality, are a threat to health-care workers, prohibitively expensive to treat, and are therefore a serious public health problem. In this Commission, we examine several aspects of drug-resistant tuberculosis. The traditional view that acquired resistance to antituberculous drugs is driven by poor compliance and programmatic failure is now being questioned, and several lines of evidence suggest that alternative mechanisms-including pharmacokinetic variability, induction of efflux pumps that transport the drug out of cells, and suboptimal drug penetration into tuberculosis lesions-are likely crucial to the pathogenesis of drug-resistant tuberculosis. These factors have implications for the design of new interventions, drug delivery and dosing mechanisms, and public health policy. We discuss epidemiology and transmission dynamics, including new insights into the fundamental biology of transmission, and we review the utility of newer diagnostic tools, including molecular tests and next-generation whole-genome sequencing, and their potential for clinical effectiveness. Relevant research priorities are highlighted, including optimal medical and surgical management, the role of newer and repurposed drugs (including bedaquiline, delamanid, and linezolid), pharmacokinetic and pharmacodynamic considerations, preventive strategies (such as prophylaxis in MDR and XDR contacts), palliative and patient-orientated care aspects, and medicolegal and ethical issues.

Removing the bottleneck in whole genome sequencing of Mycobacterium tuberculosis for rapid drug resistance analysis: a call to action

Whole genome sequencing (WGS) can provide a comprehensive analysis of Mycobacterium tuberculosis mutations that cause resistance to anti-tuberculosis drugs. With the deployment of bench-top sequencers and rapid analytical software, WGS is poised to become a useful tool to guide treatment. However, direct sequencing from clinical specimens to provide a full drug resistance profile remains a serious challenge. This article reviews current practices for extracting M. tuberculosis DNA and possible solutions for sampling sputum. Techniques under consideration include enzymatic digestion, physical disruption, chemical degradation, detergent solubilization, solvent extraction, ligand-coated magnetic beads, silica columns, and oligonucleotide pull-down baits. Selective amplification of genomic bacterial DNA in sputum prior to WGS may provide a solution, and differential lysis to reduce the levels of contaminating human DNA is also being explored. To remove this bottleneck and accelerate access to WGS for patients with suspected drug-resistant tuberculosis, it is suggested that a coordinated and collaborative approach be taken to more rapidly optimize, compare, and validate methodologies for sequencing from patient samples.

Effect of new tuberculosis diagnostic technologies on community-based intensified case finding: a multicentre randomised controlled trial

BACKGROUND

Inadequate case detection results in high levels of undiagnosed tuberculosis in sub-Saharan Africa. Data for the effect of new diagnostic tools when used for community-based intensified case finding are not available, so we investigated whether the use of sputum Xpert-MTB/RIF and the Determine TB LAM urine test in two African communities could be effective.

METHODS

In a pragmatic, randomised, parallel-group trial with individual randomisation stratified by country, we compared sputum Xpert-MTB/RIF, and if HIV-infected, the Determine TB LAM urine test (novel diagnostic group), with laboratory-based sputum smear microscopy (routine diagnostic group) for intensified case finding in communities with high tuberculosis and HIV prevalence in Cape Town, South Africa, and Harare, Zimbabwe. Participants were randomly assigned (1:1) to these groups with computer-generated allocation lists, using culture as the reference standard. In Cape Town, participants were randomised and tested at an Xpert-equipped mobile van, while in Harare, participants were driven to a local clinic where the same diagnostic tests were done. The primary endpoint was the proportion of culture-positive tuberculosis cases initiating tuberculosis treatment in each study group at 60 days. This trial is registered at ClinicalTrials.gov, number NCT01990274.

FINDINGS

Between Oct 18, 2013, and March 31, 2015, 2261 individuals were screened and 875 (39%) of these met the criteria for diagnostic testing. 439 participants were randomly assigned to the novel group and 436 to the routine group. 74 (9%) of 875 participants had confirmed tuberculosis. If late culture-based treatment initiation was excluded, more patients with culture-positive tuberculosis were initiated on treatment in the novel group at 60 days (36 [86%] of 42 in the novel group vs 18 [56%] of 32 in the routine group). Thus the difference in the proportion initiating treatment between groups was 29% (95% CI 9-50, p=0·0047) and 53% more patients initiated therapy in the novel diagnostic group than in the routine diagnostic group. One culture-positive patient was treated based only on a positive LAM test.

INTERPRETATION

Compared with traditional tools, Xpert-MTB/RIF for community-based intensified case finding in HIV and tuberculosis-endemic settings increased the proportion of patients initiating treatment. By contrast, urine LAM testing was not found to be useful for intensive case finding in this setting.

FUNDING

European and Developing Countries Clinical Trials Partnership and South African Medical Research Council.

The variability and reproducibility of whole genome sequencing technology for detecting resistance to anti-tuberculous drugs

BACKGROUND

The emergence of resistance to anti-tuberculosis drugs is a serious and growing threat to public health. Next-generation sequencing is rapidly gaining traction as a diagnostic tool for investigating drug resistance in Mycobacterium tuberculosis to aid treatment decisions. However, there are few little data regarding the precision of such sequencing for assigning resistance profiles.

METHODS

We investigated two sequencing platforms (Illumina MiSeq, Ion Torrent PGM™) and two rapid analytic pipelines (TBProfiler, Mykrobe predictor) using a well characterised reference strain (H37Rv) and clinical isolates from patients with tuberculosis resistant to up to 13 drugs. Results were compared to phenotypic drug susceptibility testing. To assess analytical robustness individual DNA samples were subjected to repeated sequencing.

RESULTS

The MiSeq and Ion PGM systems accurately predicted drug-resistance profiles and there was high reproducibility between biological and technical sample replicates. Estimated variant error rates were low (MiSeq 1 per 77 kbp, Ion PGM 1 per 41 kbp) and genomic coverage high (MiSeq 51-fold, Ion PGM 53-fold). MiSeq provided superior coverage in GC-rich regions, which translated into incremental detection of putative genotypic drug-specific resistance, including for resistance to para-aminosalicylic acid and pyrazinamide. The TBProfiler bioinformatics pipeline was concordant with reported phenotypic susceptibility for all drugs tested except pyrazinamide and para-aminosalicylic acid, with an overall concordance of 95.3%. When using the Mykrobe predictor concordance with phenotypic testing was 73.6%.

CONCLUSIONS

We have demonstrated high comparative reproducibility of two sequencing platforms, and high predictive ability of the TBProfiler mutation library and analytical pipeline, when profiling resistance to first- and second-line anti-tuberculosis drugs. However, platform-specific variability in coverage of some genome regions may have implications for predicting resistance to specific drugs. These findings may have implications for future clinical practice and thus deserve further scrutiny, set within larger studies and using updated mutation libraries.

Alternate efflux pump mechanism may contribute to drug resistance in extensively drug-resistant isolates of Mycobacterium tuberculosis

INTRODUCTION

Extensively drug-resistant tuberculosis (XDR-TB) has emerged as one of the biggest threats to public health and TB control programs worldwide. XDR-TB is caused by Mycobacterium tuberculosis (MTB) strains resistant to rifampin and isoniazid, as well as to a fluoroquinolone and to at least one injectable aminoglycoside. Drug resistance in MTB has primarily been associated with single nucleotide polymorphisms (SNPs) in particular genes. However, it has also been shown that efflux pumps may play a role in resistance of MTB. Upregulation of drug efflux pumps can decrease the intracellular concentration of drugs and reduce their efficacy.

METHODS

Whole genome sequencing was performed on 32 XDR-TB clinical isolates. Sequence data were used to investigate SNPs in efflux pump genes as compared with the H37Rv reference genome.

RESULTS

Of the XDR MTB strains, eight (21.62%) were wild type for rpsL, rrs (500 region), and gidB genes, but had non-synonymous (ns) SNPs (aspartic acid to histidine) in the drrA efflux pump gene at position 3273138. Three of eight (37.5%) XDR MTB strains, wild type for rpsL, rrs (500 region), gidB, and gyrB genes were phenotypically streptomycin sensitive and five (62.5%) XDR MTB strains were streptomycin resistant, while all XDR MTB strains, wild type for rpsL, rrs, gidB, and gyrB genes were resistant to fluoroquinolone (ofloxacin) and ethambutol. In addition, three XDR MTB strains wild type for rpsL, rrs, gidB, and drrA genes showed nsSNPs (isoleucine to valine) in the major facilitator superfamily, Rv1634 efflux pump gene at position 1839306.

CONCLUSION

Our data show an nsSNP in the drrA efflux pump gene that may result in upregulation of drug efflux mechanisms in MTB strains. It is therefore imperative to understand the mechanism of efflux and its role in drug resistance, which will enable the identification of new drug targets and development of new drug regimens to counteract the drug efflux mechanism of MTB.

Turkish and Japanese Mycobacterium tuberculosis sublineages share a remote common ancestor

Two geographically distant M. tuberculosis sublineages, Tur from Turkey and T3-Osaka from Japan, exhibit partially identical genotypic signatures (identical 12-loci MIRU-VNTR profiles, distinct spoligotyping patterns). We investigated T3-Osaka and Tur sublineages characteristics and potential genetic relatedness, first using MIRU-VNTR locus analysis on 21 and 25 samples of each sublineage respectively, and second comparing Whole Genome Sequences of 8 new samples to public data from 45 samples uncovering human tuberculosis diversity. We then tried to date their Most Recent Common Ancestor (MRCA) using three calibrations of SNP accumulation rate (long-term=0.03SNP/genome/year, derived from a tuberculosis ancestor of around 70,000years old; intermediate=0.2SNP/genome/year derived from a Peruvian mummy; short-term=0.5SNP/genome/year). To disentangle between these scenarios, we confronted the corresponding divergence times with major human history events and knowledge on human genetic divergence. We identified relatively high intrasublineage diversity for both T3-Osaka and Tur. We definitively proved their monophyly; the corresponding super-sublineage (referred to as "T3-Osa-Tur") shares a common ancestor with T3-Ethiopia and Ural sublineages but is only remotely related to other Euro-American sublineages such as X, LAM, Haarlem and S. The evolutionary scenario based on long-term evolution rate being valid until T3-Osa-Tur MRCA was not supported by Japanese fossil data. The evolutionary scenario relying on short-term evolution rate since T3-Osa-Tur MRCA was contradicted by human history and potential traces of past epidemics. T3-Osaka and Tur sublineages were found likely to have diverged between 800y and 2000years ago, potentially at the time of Mongol Empire. Altogether, this study definitively proves a strong genetic link between Turkish and Japanese tuberculosis. It provides a first hypothesis for calibrating TB Euro-American lineage molecular clock; additional studies are needed to reliably date events corresponding to intermediate depths in tuberculosis phylogeny.

Perspectives on Advances in Tuberculosis Diagnostics, Drugs, and Vaccines

Despite concerted efforts over the past 2 decades at developing new diagnostics, drugs, and vaccines with expanding pipelines, tuberculosis remains a global emergency. Several novel diagnostic technologies show promise of better point-of-care rapid tests for tuberculosis including nucleic acid-based amplification tests, imaging, and breath analysis of volatile organic compounds. Advances in new and repurposed drugs for use in multidrug-resistant (MDR) or extensively drug-resistant (XDR) tuberculosis have focused on development of several new drug regimens and their evaluation in clinical trials and now influence World Health Organization guidelines. Since the failure of the MVA85A vaccine 2 years ago, there have been no new tuberculosis vaccine candidates entering clinical testing. The current status quo of the lengthy treatment duration and poor treatment outcomes associated with MDR/XDR tuberculosis and with comorbidity of tuberculosis with human immunodeficiency virus and noncommunicable diseases is unacceptable. New innovations and political and funder commitment for early rapid diagnosis, shortening duration of therapy, improving treatment outcomes, and prevention are urgently required.

Regulatory In Vitro Diagnostics Landscape in Africa: Update on Regional Activities

Improved diagnostic tests for tuberculosis case detection are urgently needed that are affordable, robust, and easy to use so that they can be implemented widely. The mandate of national regulatory authorities is to ensure the safety and effectiveness of diagnostics, protecting the population against unsafe products while expediting access to beneficial new devices. However, regulatory approval processes in the developing world are often complex, lengthy, and not transparent. Recent progress in building regulatory capacity using harmonized approaches will reduce duplication in clinical performance studies and manufacturing audits, facilitate information sharing through trust and mutual confidence building, and ultimately improve efficiency. These savings can be passed onto the consumers in the form of more affordable pricing and allowing new high-quality tests for tuberculosis to be introduced more quickly and without delay.

Mycobacterium tuberculosis whole genome sequencing and protein structure modelling provides insights into anti-tuberculosis drug resistance

BACKGROUND

Combating the spread of drug resistant tuberculosis is a global health priority. Whole genome association studies are being applied to identify genetic determinants of resistance to anti-tuberculosis drugs. Protein structure and interaction modelling are used to understand the functional effects of putative mutations and provide insight into the molecular mechanisms leading to resistance.

METHODS

To investigate the potential utility of these approaches, we analysed the genomes of 144 Mycobacterium tuberculosis clinical isolates from The Special Programme for Research and Training in Tropical Diseases (TDR) collection sourced from 20 countries in four continents. A genome-wide approach was applied to 127 isolates to identify polymorphisms associated with minimum inhibitory concentrations for first-line anti-tuberculosis drugs. In addition, the effect of identified candidate mutations on protein stability and interactions was assessed quantitatively with well-established computational methods.

RESULTS

The analysis revealed that mutations in the genes rpoB (rifampicin), katG (isoniazid), inhA-promoter (isoniazid), rpsL (streptomycin) and embB (ethambutol) were responsible for the majority of resistance observed. A subset of the mutations identified in rpoB and katG were predicted to affect protein stability. Further, a strong direct correlation was observed between the minimum inhibitory concentration values and the distance of the mutated residues in the three-dimensional structures of rpoB and katG to their respective drugs binding sites.

CONCLUSIONS

Using the TDR resource, we demonstrate the usefulness of whole genome association and convergent evolution approaches to detect known and potentially novel mutations associated with drug resistance. Further, protein structural modelling could provide a means of predicting the impact of polymorphisms on drug efficacy in the absence of phenotypic data. These approaches could ultimately lead to novel resistance mutations to improve the design of tuberculosis control measures, such as diagnostics, and inform patient management.

Recombination in pe/ppe genes contributes to genetic variation in Mycobacterium tuberculosis lineages

Background

Approximately 10 % of the Mycobacterium tuberculosis genome is made up of two families of genes that are poorly characterized due to their high GC content and highly repetitive nature. The PE and PPE families are typified by their highly conserved N-terminal domains that incorporate proline-glutamate (PE) and proline-proline-glutamate (PPE) signature motifs. They are hypothesised to be important virulence factors involved with host-pathogen interactions, but their high genetic variability and complexity of analysis means they are typically disregarded in genome studies.

Results

To elucidate the structure of these genes, 518 genomes from a diverse international collection of clinical isolates were de novo assembled. A further 21 reference M. tuberculosis complex genomes and long read sequence data were used to validate the approach. SNP analysis revealed that variation in the majority of the 168 pe/ppe genes studied was consistent with lineage. Several recombination hotspots were identified, notably pe_pgrs3 and pe_pgrs17. Evidence of positive selection was revealed in 65 pe/ppe genes, including epitopes potentially binding to major histocompatibility complex molecules.

Conclusions

This, the first comprehensive study of the pe and ppe genes, provides important insight into M. tuberculosis diversity and has significant implications for vaccine development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2467-y) contains supplementary material, which is available to authorized users.

Whole Genome Sequencing Shows a Low Proportion of Tuberculosis Disease Is Attributable to Known Close Contacts in Rural Malawi

Background

The proportion of tuberculosis attributable to transmission from close contacts is not well known. Comparison of the genome of strains from index patients and prior contacts allows transmission to be confirmed or excluded.

Methods

In Karonga District, Malawi, all tuberculosis patients are asked about prior contact with others with tuberculosis. All available strains from culture-positive patients were sequenced. Up to 10 single nucleotide polymorphisms between index patients and their prior contacts were allowed for confirmation, and ≥ 100 for exclusion. The population attributable fraction was estimated from the proportion of confirmed transmissions and the proportion of patients with contacts.

Results

From 1997–2010 there were 1907 new culture-confirmed tuberculosis patients, of whom 32% reported at least one family contact and an additional 11% had at least one other contact; 60% of contacts had smear-positive disease. Among case-contact pairs with sequences available, transmission was confirmed from 38% (62/163) smear-positive prior contacts and 0/17 smear-negative prior contacts. Confirmed transmission was more common in those related to the prior contact (42.4%, 56/132) than in non-relatives (19.4%, 6/31, p = 0.02), and in those with more intense contact, to younger index cases, and in more recent years. The proportion of tuberculosis attributable to known contacts was estimated to be 9.4% overall.

Conclusions

In this population known contacts only explained a small proportion of tuberculosis cases. Even those with a prior family contact with smear positive tuberculosis were more likely to have acquired their infection elsewhere.

The draft genome of Mycobacterium aurum, a potential model organism for investigating drugs against Mycobacterium tuberculosis and Mycobacterium leprae

Mycobacterium aurum (M. aurum) is an environmental mycobacteria that has previously been used in studies of anti-mycobacterial drugs due to its fast growth rate and low pathogenicity. The M. aurum genome has been sequenced and assembled into 46 contigs, with a total length of 6.02Mb containing 5684 annotated protein-coding genes. A phylogenetic analysis using whole genome alignments positioned M. aurum close to Mycobacterium vaccae and Mycobacterium vanbaalenii, within a clade related to fast-growing mycobacteria. Large-scale genomic rearrangements were identified by comparing the M. aurum genome to those of Mycobacterium tuberculosis and Mycobacterium leprae. M. aurum orthologous genes implicated in resistance to anti-tuberculosis drugs in M. tuberculosis were observed. The sequence identity at the DNA level varied from 68.6% for pncA (pyrazinamide drug-related) to 96.2% for rrs (streptomycin, capreomycin). We observed two homologous genes encoding the catalase-peroxidase enzyme (katG) that is associated with resistance to isoniazid. Similarly, two embB homologues were identified in the M. aurum genome. In addition to describing for the first time the genome of M. aurum, this work provides a resource to aid the use of M. aurum in studies to develop improved drugs for the pathogenic mycobacteria M. tuberculosis and M. leprae.

Rapid determination of anti-tuberculosis drug resistance from whole-genome sequences

Mycobacterium tuberculosis drug resistance (DR) challenges effective tuberculosis disease control. Current molecular tests examine limited numbers of mutations, and although whole genome sequencing approaches could fully characterise DR, data complexity has restricted their clinical application. A library (1,325 mutations) predictive of DR for 15 anti-tuberculosis drugs was compiled and validated for 11 of them using genomic-phenotypic data from 792 strains. A rapid online ‘TB-Profiler’ tool was developed to report DR and strain-type profiles directly from raw sequences. Using our DR mutation library, in silico diagnostic accuracy was superior to some commercial diagnostics and alternative databases. The library will facilitate sequence-based drug-susceptibility testing.

Diagnostics for Developing Countries

Improving the availability of high quality diagnostic tests for infectious diseases is a global priority. Lack of access by people living in low income countries may deprive them of life saving treatment and reduces opportunities to prevent onward transmission and spread of the disease. Diagnostic laboratories are often poorly resourced in developing countries, and sparsely distributed. Improved access may be achieved by using tests that do not require laboratory support, including rapid tests for use at the point-of-care. Despite increased interest, few new in vitro diagnostic (IVD) products reach the majority populations in low income countries. Barriers to uptake include cost and lack of robustness, with reduced test performances due to environmental pressures such as high ambient temperatures or dust. In addition to environmental factors test developers must consider the local epidemiology. Confounding conditions such as immunosuppression or variations in antigen presentation or genotype can affect test performance. Barriers to product development include access to finance to establish manufacturing capacity and cover the costs of market entry for new devices. Costs and delays may be inflated by current regulatory preregistration processes to ensure product safety and quality, and more harmonized approaches are needed.

Lifestyle, attitudes and needs of uncured XDR-TB patients living in the communities of South Africa: a qualitative study

OBJECTIVE

Patient-level data are required to inform strategies interrupting transmission and default in patients with extensively drug-resistant TB (XDR-TB) to improve models of care and identify potential routes of transmission. We therefore explored the experiences, lifestyle, attitudes and needs of patients with uncured XDR-TB, who failed or interrupted therapy, living without treatment in the community.

METHODS

We conducted in-depth interviews with 12 community-based patients from South Africa. Family members were interviewed when patients were unavailable. Interviews were analysed using inductive thematic analysis.

RESULTS

The thematic experiences identified from the interviews were as follows: (i) living with but not being cured of XDR-TB, (ii) altered lifestyle in the community, (iii) experiences with community health care, (iv) local community members, and (v) wants and needs. Patients identified mistrust in health care, futility of treatment regimens, a need for a purpose in life and subsistence as major concerns. Restriction of living in the community for patients whose treatment had failed resulted in self-imposed isolation. Defaulters focused more on the never-ending drug regimen and bad experiences with health care contributing to non-adherence. Family members emphasised an under-recognised experience of unforeseen burden, obligation, worry and discomfort. Lack of knowledge and lack of concern about transmission was evident.

CONCLUSION

Current models of care are not adequately meeting the needs of patients with uncured XDR-TB and relatives. These data inform the need for community-based palliative care, vocational facilities to improve economic opportunities, home-based infection control and improved psychosocial support to increase patient adherence, reduce transmission, provide income and relieve the burden on family members.

PhyTB: Phylogenetic tree visualisation and sample positioning for M. tuberculosis

Background

Phylogenetic-based classification of M. tuberculosis and other bacterial genomes is a core analysis for studying evolutionary hypotheses, disease outbreaks and transmission events. Whole genome sequencing is providing new insights into the genomic variation underlying intra- and inter-strain diversity, thereby assisting with the classification and molecular barcoding of the bacteria. One roadblock to strain investigation is the lack of user-interactive solutions to interrogate and visualise variation within a phylogenetic tree setting.

Results

We have developed a web-based tool called PhyTB (http://pathogenseq.lshtm.ac.uk/phytblive/index.php) to assist phylogenetic tree visualisation and identification of M. tuberculosis clade-informative polymorphism. Variant Call Format files can be uploaded to determine a sample position within the tree. A map view summarises the geographical distribution of alleles and strain-types. The utility of the PhyTB is demonstrated on sequence data from 1,601 M. tuberculosis isolates.

Conclusion

PhyTB contextualises M. tuberculosis genomic variation within epidemiological, geographical and phylogenic settings. Further tool utility is possible by incorporating large variants and phenotypic data (e.g. drug-resistance profiles), and an assessment of genotype-phenotype associations. Source code is available to develop similar websites for other organisms (http://sourceforge.net/projects/phylotrack).

Development of sample clean up methods for the analysis of Mycobacterium tuberculosis methyl mycocerosate biomarkers in sputum extracts by gas chromatography-mass spectrometry

•

We developed a sample clean-up method to detect tuberculosis from sputum by GC–MS.

•

Biomarkers recovered: 64–70% (standards solution), and 36–68% (sputum extracts).

•

Cholesterol removed: 93–98% (standards solution) and 62–92% (sputum extracts).

•

Less cholesterol in the filtered extracts avoids overloading of the analytical system.

•

Analyzing large sample batches will need fewer interruptions for system cleaning.

A proof of principle gas chromatography–mass spectrometry method is presented, in combination with clean up assays, aiming to improve the analysis of methyl mycocerosate tuberculosis biomarkers from sputum. Methyl mycocerosates are generated from the transesterification of phthiocerol dimycocerosates (PDIMs), extracted in petroleum ether from sputum of tuberculosis suspect patients. When a high matrix background is present in the sputum extracts, the identification of the chromatographic peaks corresponding to the methyl derivatives of PDIMs analytes may be hindered by the closely eluting methyl ether of cholesterol, usually an abundant matrix constituent frequently present in sputum samples. The purification procedures involving solid phase extraction (SPE) based methods with both commercial Isolute-Florisil cartridges, and purpose designed molecularly imprinted polymeric materials (MIPs), resulted in cleaner chromatograms, while the mycocerosates are still present. The clean-up performed on solutions of PDIMs and cholesterol standards in petroleum ether show that, depending on the solvent mix and on the type of SPE used, the recovery of PDIMs is between 64 and 70%, whilst most of the cholesterol is removed from the system. When applied to petroleum ether extracts from representative sputum samples, the clean-up procedures resulted in recoveries of 36–68% for PDIMs, allowing some superior detection of the target analytes.

Large-scale whole genome sequencing of M. tuberculosis provides insights into transmission in a high prevalence area

To improve understanding of the factors influencing tuberculosis transmission and the role of pathogen variation, we sequenced all available specimens from patients diagnosed over 15 years in a whole district in Malawi. Mycobacterium tuberculosis lineages were assigned and transmission networks constructed, allowing ≤10 single nucleotide polymorphisms (SNPs) difference. We defined disease as due to recent infection if the network-determined source was within 5 years, and assessed transmissibility from forward transmissions resulting in disease. High-quality sequences were available for 1687 disease episodes (72% of all culture-positive episodes): 66% of patients linked to at least one other patient. The between-patient mutation rate was 0.26 SNPs/year (95% CI 0.21–0.31). We showed striking differences by lineage in the proportion of disease due to recent transmission and in transmissibility (highest for lineage-2 and lowest for lineage-1) that were not confounded by immigration, HIV status or drug resistance. Transmissions resulting in disease decreased markedly over time.

DOI:http://dx.doi.org/10.7554/eLife.05166.001

Tuberculosis is an important public health threat around the globe and is particularly common in developing countries. It is difficult to control the spread of the disease because the bacteria that cause it can spread when an infected individual coughs or sneezes. It may take years for an infected individual to develop symptoms of tuberculosis so it can be hard to trace the source of an outbreak, and people infected with HIV are particularly susceptible to the disease.

The bacterium that causes the majority of cases of tuberculosis is called Mycobacterium tuberculosis. There are several different varieties or ‘lineages’ of M. tuberculosis, and it is thought that they may vary in their ability to spread and cause disease. However, the results of previous studies have been inconsistent and there also seems to be a lot of variation between strains within the same lineage.

In this study, Guerra-Assunção et al. used an approach called whole genome sequencing alongside more traditional methods to study the spread of tuberculosis in Malawi. They sequenced the genomes of every available sample of M. tuberculosis collected from patients in the Karonga district of Malawi over a 15-year period. This produced high-quality DNA sequence data about the bacteria responsible for almost 1700 cases of disease.

Using this massive amount of data, Guerra-Assunção et al. constructed networks that showed how the bacteria had spread in the community. This revealed that there were differences between the ability of the various M. tuberculosis lineages to cause disease and to spread in communities. For example, lineage 1 was less likely than the other lineages to cause disease soon after infecting an individual and was less able to spread.

The data also show that the proportion of cases of disease due to recent infection declined substantially during the 15-year period. This indicates that the tuberculosis and HIV control programmes in the area have been successful.

Guerra-Assunção et al.'s findings show that it is possible to understand how tuberculosis is transmitted on a large scale. The next challenge is to understand why the lineages differ in their ability to cause disease and spread between individuals.

DOI:http://dx.doi.org/10.7554/eLife.05166.002

Characterization of genomic variations in SNPs of PE_PGRS genes reveals deletions and insertions in extensively drug resistant (XDR) M. tuberculosis strains from Pakistan

BACKGROUND

Mycobacterium tuberculosis (MTB) PE_PGRS genes belong to the PE multigene family. Although the function of PE_PGRS genes is unknown, it is hypothesized that the PE_PGRS genes may be associated with antigenic variability in MTB.

MATERIAL AND METHODS

Whole genome sequencing analysis was performed on (n=37) extensively drug-resistant (XDR) MTB strains from Pakistan, which included Lineage 1 (East African Indian, n=2); Other lineage 1 (n=3); Lineage 3 (Central Asian, n=24); Other lineage 3 (n=4); Lineage 4 (X3, n=1) and T group (n=3) MTB strains.

RESULTS

There were 107 SNPs identified from the analysis of 42 PE_PGRS genes; of these, 13 were non-synonymous SNPs (nsSNPs). The nsSNPs identified in PE_PGRS genes - 6, 9 and 10 - were common in all EAI, CAS, Other lineages (1 and 3), T1 and X3. Deletions (DELs) in PE_PGRS genes - 3 and 19 - were observed in 17 (80.9%) CAS1 and 6 (85.7%) in Other lineages (1 and 3) XDR MTB strains, while DELs in the PE_PGRS49 were observed in all CAS1, CAS, CAS2 and Other lineages (1 and 3) XDR MTB strains. All CAS, EAI and Other lineages (1 and 3) strains showed insertions (INS) in PE_PGRS6 gene, while INS in the PE_PGRS genes 19 and 33 were observed in 20 (95.2%) CAS1, all CAS, CAS2, EAI and Other lineages (1 and 3) XDR MTB strains.

CONCLUSION

Genetic diversity in PE_PGRS genes contributes to antigenic variability and may result in increased immunogenicity of strains. This is the first study identifying variations in nsSNPs and INDELs in the PE_PGRS genes of XDR-TB strains from Pakistan. It highlights common genetic variations which may contribute to persistence.

A proportion of mutations fixed in the genomes of in vitro selected isogenic drug-resistant Mycobacterium tuberculosis mutants can be detected as minority variants in the parent culture

We studied genomic variation in a previously selected collection of isogenic Mycobacterium tuberculosis laboratory strains subjected to one or two rounds of antibiotic selection. Whole genome sequencing analysis identified eleven single, unique mutations (four synonymous, six non-synonymous, one intergenic), in addition to drug resistance-conferring mutations, that were fixed in the genomes of six monoresistant strains. Eight loci, present as minority variants (five non-synonymous, three synonymous) in the genome of the susceptible parent strain, became fixed in the genomes of multiple daughter strains. None of these mutations are known to be involved with drug resistance. Our results confirm previously observed genomic stability for M. tuberculosis, although the parent strain had accumulated allelic variants at multiple locations in an antibiotic-free in vitro environment. It is therefore likely to assume that these so-called hitchhiking mutations were co-selected and fixed in multiple daughter strains during antibiotic selection. The presence of multiple allelic variations, accumulated under non-selective conditions, which become fixed during subsequent selective steps, deserves attention. The wider availability of 'deep' sequencing methods could help to detect multiple bacterial (sub)populations within patients with high resolution and would therefore be useful in assisting in the detailed investigation of transmission chains.

Unraveling Mycobacterium tuberculosis genomic diversity and evolution in Lisbon, Portugal, a highly drug resistant setting

Background

Multidrug- (MDR) and extensively drug resistant (XDR) tuberculosis (TB) presents a challenge to disease control and elimination goals. In Lisbon, Portugal, specific and successful XDR-TB strains have been found in circulation for almost two decades.

Results

In the present study we have genotyped and sequenced the genomes of 56 Mycobacterium tuberculosis isolates recovered mostly from Lisbon. The genotyping data revealed three major clusters associated with MDR-TB, two of which are associated with XDR-TB. Whilst the genomic data contributed to elucidate the phylogenetic positioning of circulating MDR-TB strains, showing a high predominance of a single SNP cluster group 5. Furthermore, a genome-wide phylogeny analysis from these strains, together with 19 publicly available genomes of Mycobacterium tuberculosis clinical isolates, revealed two major clades responsible for M/XDR-TB in the region: Lisboa3 and Q1 (LAM).

The data presented by this study yielded insights on microevolution and identification of novel compensatory mutations associated with rifampicin resistance in rpoB and rpoC. The screening for other structural variations revealed putative clade-defining variants. One deletion in PPE41, found among Lisboa3 isolates, is proposed to contribute to immune evasion and as a selective advantage. Insertion sequence (IS) mapping has also demonstrated the role of IS6110 as a major driver in mycobacterial evolution by affecting gene integrity and regulation.

Conclusions

Globally, this study contributes with novel genome-wide phylogenetic data and has led to the identification of new genomic variants that support the notion of a growing genomic diversity facing both setting and host adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-991) contains supplementary material, which is available to authorized users.

Regulation of medical diagnostics and medical devices in the East African community partner states

Background

Medical devices and in vitro diagnostic tests (IVD) are vital components of health delivery systems but access to these important tools is often limited in Africa. The regulation of health commodities by National Regulatory Authorities is intended to ensure their safety and quality whilst ensuring timely access to beneficial new products. Streamlining and harmonizing regulatory processes may reduce delays and unnecessary expense and improve access to new products. Whereas pharmaceutical products are widely regulated less attention has been placed on the regulation of other health products. A study was undertaken to assess regulation of medical diagnostics and medical devices across Partner States of the East African Community (EAC).

Methods

Data was collected during October 2012 through desk based review of documents and field research, including face to face interviews with the assistance of a structured questionnaire with closed and open ended questions. Key areas addressed were (i) existence and role of National Regulatory Authorities; (ii) policy and legal framework for regulation; (iii) premarket control; (iv) marketing controls; (v) post-marketing control and vigilance; (vi) country capacity for regulation; (vii) country capacity for evaluation studies for IVD and (viii) priorities and capacity building for harmonization in EAC Partner States.

Results

Control of medical devices and IVDs in EAC Partner States is largely confined to national disease programmes such as tuberculosis, HIV and malaria. National Regulatory Authorities for pharmaceutical products do not have the capacity to regulate medical devices and in some countries laboratory based organisations are mandated to ensure quality of products used. Some activities to evaluate IVDs are performed in research laboratories but post market surveillance is rare. Training in key areas is considered essential to strengthening regulatory capacity for IVDs and other medical devices.

Conclusions

Regulation of medical devices and in vitro diagnostics has been neglected in EAC Partner States. Regulation is weak across the region, and although the majority of States have a legal mandate to regulate medical devices there is limited capacity to do so. Streamlining regulation in the EAC is seen as a positive aspiration with diagnostic tests considered a priority area for harmonisation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12913-014-0524-2) contains supplementary material, which is available to authorized users.

Recurrence due to relapse or reinfection with Mycobacterium tuberculosis: a whole-genome sequencing approach in a large, population-based cohort with a high HIV infection prevalence and active follow-up

BACKGROUND

Recurrent tuberculosis is a major health burden and may be due to relapse with the original strain or reinfection with a new strain.

METHODS

In a population-based study in northern Malawi, patients with tuberculosis diagnosed from 1996 to 2010 were actively followed after the end of treatment. Whole-genome sequencing with approximately 100-fold coverage was performed on all available cultures. Results of IS6110 restriction fragment-length polymorphism analyses were available for cultures performed up to 2008.

RESULTS

Based on our data, a difference of ≤10 single-nucleotide polymorphisms (SNPs) was used to define relapse, and a difference of >100 SNPs was used to define reinfection. There was no evidence of mixed infections among those classified as reinfections. Of 1471 patients, 139 had laboratory-confirmed recurrences: 55 had relapse, and 20 had reinfection; for 64 type of recurrence was unclassified. Almost all relapses occurred in the first 2 years. Human immunodeficiency virus infection was associated with reinfection but not relapse. Relapses were associated with isoniazid resistance, treatment before 2007, and lineage-3 strains. We identified several gene variants associated with relapse. Lineage-2 (Beijing) was overrepresented and lineage-1 underrepresented among the reinfecting strains (P = .004).

CONCLUSIONS

While some of the factors determining recurrence depend on the patient and their treatment, differences in the Mycobacterium tuberculosis genome appear to have a role in both relapse and reinfection.

PolyTB: a genomic variation map for Mycobacterium tuberculosis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) is the second major cause of death from an infectious disease worldwide. Recent advances in DNA sequencing are leading to the ability to generate whole genome information in clinical isolates of M. tuberculosis complex (MTBC). The identification of informative genetic variants such as phylogenetic markers and those associated with drug resistance or virulence will help barcode Mtb in the context of epidemiological, diagnostic and clinical studies. Mtb genomic datasets are increasingly available as raw sequences, which are potentially difficult and computer intensive to process, and compare across studies. Here we have processed the raw sequence data (>1500 isolates, eight studies) to compile a catalogue of SNPs (n = 74,039, 63% non-synonymous, 51.1% in more than one isolate, i.e. non-private), small indels (n = 4810) and larger structural variants (n = 800). We have developed the PolyTB web-based tool (http://pathogenseq.lshtm.ac.uk/polytb) to visualise the resulting variation and important meta-data (e.g. in silico inferred strain-types, location) within geographical map and phylogenetic views. This resource will allow researchers to identify polymorphisms within candidate genes of interest, as well as examine the genomic diversity and distribution of strains. PolyTB source code is freely available to researchers wishing to develop similar tools for their pathogen of interest.

Improving access to new diagnostics through harmonised regulation: priorities for action

A new generation of diagnostic tests is being developed for use at the point of care that could save lives and reduce the spread of infectious diseases through early detection and treatment. It is important that patients in developing countries have access to these products at affordable prices and without delay. Regulation of medical products is intended to ensure safety and quality whilst balancing the need for timely access to beneficial new products. Current regulatory oversight of diagnostic tests in developing countries is highly variable and weak regulation allows poor-quality tests to enter the market. However, inefficient or overzealous regulation results in unnecessary delays, increases costs and acts as a barrier to innovation and market entry. Setting international standards and streamlining the regulatory process could reduce these barriers. Four priority activities have been identified where convergence of standards and protocols or joint review of data would be advantageous: (1) adoption of a common registration file for pre-market approval; (2) convergence of quality standards for manufacturing site inspections; (3) use of common evaluation protocols, as well as joint review of data, to reduce unnecessary duplication of lengthy and costly clinical performance studies; and (4) use of networks of laboratories for post-market surveillance in order to monitor ongoing quality of diagnostic devices. The adoption and implementation of such measures in developing countries could accelerate access to new diagnostic tests that are safe and affordable.

Treatment outcomes of new tuberculosis patients hospitalized in Kampala, Uganda: a prospective cohort study

BACKGROUND

In most resource limited settings, new tuberculosis (TB) patients are usually treated as outpatients. We sought to investigate the reasons for hospitalisation and the predictors of poor treatment outcomes and mortality in a cohort of hospitalized new TB patients in Kampala, Uganda.

METHODS AND FINDINGS

Ninety-six new TB patients hospitalised between 2003 and 2006 were enrolled and followed for two years. Thirty two were HIV-uninfected and 64 were HIV-infected. Among the HIV-uninfected, the commonest reasons for hospitalization were low Karnofsky score (47%) and need for diagnostic evaluation (25%). HIV-infected patients were commonly hospitalized due to low Karnofsky score (72%), concurrent illness (16%) and diagnostic evaluation (14%). Eleven HIV uninfected patients died (mortality rate 19.7 per 100 person-years) while 41 deaths occurred among the HIV-infected patients (mortality rate 46.9 per 100 person years). In all patients an unsuccessful treatment outcome (treatment failure, death during the treatment period or an unknown outcome) was associated with duration of TB symptoms, with the odds of an unsuccessful outcome decreasing with increasing duration. Among HIV-infected patients, an unsuccessful treatment outcome was also associated with male sex (P = 0.004) and age (P = 0.034). Low Karnofsky score (aHR = 8.93, 95% CI 1.88 - 42.40, P = 0.001) was the only factor significantly associated with mortality among the HIV-uninfected. Mortality among the HIV-infected was associated with the composite variable of CD4 and ART use, with patients with baseline CD4 below 200 cells/µL who were not on ART at a greater risk of death than those who were on ART, and low Karnofsky score (aHR = 2.02, 95% CI 1.02 - 4.01, P = 0.045).

CONCLUSION

Poor health status is a common cause of hospitalisation for new TB patients. Mortality in this study was very high and associated with advanced HIV Disease and no use of ART.

Elucidating emergence and transmission of multidrug-resistant tuberculosis in treatment experienced patients by whole genome sequencing

BACKGROUND

Understanding the emergence and spread of multidrug-resistant tuberculosis (MDR-TB) is crucial for its control. MDR-TB in previously treated patients is generally attributed to the selection of drug resistant mutants during inadequate therapy rather than transmission of a resistant strain. Traditional genotyping methods are not sufficient to distinguish strains in populations with a high burden of tuberculosis and it has previously been difficult to assess the degree of transmission in these settings. We have used whole genome analysis to investigate M. tuberculosis strains isolated from treatment experienced patients with MDR-TB in Uganda over a period of four years.

METHODS AND FINDINGS

We used high throughput genome sequencing technology to investigate small polymorphisms and large deletions in 51 Mycobacterium tuberculosis samples from 41 treatment-experienced TB patients attending a TB referral and treatment clinic in Kampala. This was a convenience sample representing 69% of MDR-TB cases identified over the four year period. Low polymorphism was observed in longitudinal samples from individual patients (2-15 SNPs). Clusters of samples with less than 50 SNPs variation were examined. Three clusters comprising a total of 8 patients were found with almost identical genetic profiles, including mutations predictive for resistance to rifampicin and isoniazid, suggesting transmission of MDR-TB. Two patients with previous drug susceptible disease were found to have acquired MDR strains, one of which shared its genotype with an isolate from another patient in the cohort.

CONCLUSIONS

Whole genome sequence analysis identified MDR-TB strains that were shared by more than one patient. The transmission of multidrug-resistant disease in this cohort of retreatment patients emphasises the importance of early detection and need for infection control. Consideration should be given to rapid testing for drug resistance in patients undergoing treatment to monitor the emergence of resistance and permit early intervention to avoid onward transmission.

Development and optimization of a gas chromatography/mass spectrometry method for the analysis of thermochemolytic degradation products of phthiocerol dimycocerosate waxes found in Mycobacterium tuberculosis

RATIONALE

The phthiocerol dimycocerosates (PDIMs) are certain stable and hydrophobic waxes found in the cell membrane of Mycobacterium tuberculosis, bacteria that cause an infectious disease of growing concern worldwide. Previous studies report the analysis of derivatives of the hydrolysed PDIMs from biological samples, following complex extraction and offline derivatization of PDIMs biomarkers, prior to their analysis by gas chromatography/mass spectrometry (GC/MS).

METHODS

We developed and optimized a GC/MS method based on selected ion monitoring (SIM) to detect the derivatives produced via the thermally assisted hydrolysis and methylation (THM) of the PDIMs from the cell membrane of M. tuberculosis. The extraction of PDIMs from culture is simple, and their thermochemolysis is carried out automatically online, thus avoiding the time-consuming derivatization steps of hydrolysis and esterification, usually performed offline.

RESULTS

For standard PDIMs in petroleum ether, our optimized method gave an excellent linearity (R(2) = 0.99) at concentrations between 0.172 and 27.5 ng/mL, a good precision (RSD = 11.42%), and a limit of detection (LOD) of 100 pg/mL. For the PDIMs extracted from dilutions of M. tuberculosis culture, the method gave good linearity (R(2) = 0.9685) and an estimated LOD of 400 CFU/mL (CFU = colony forming units) in sterile distilled water.

CONCLUSIONS

A GC/MS(SIM) method is presented for the rapid and quantitative detection of M. tuberculosis, based on the online thermochemolysis of lipidic biomarkers extracted from the bacterial culture. The method has the potential to be applied in human and veterinary clinical laboratories for the rapid diagnosis of tuberculosis in infected biological samples.

Advances in tuberculosis diagnostics: the Xpert MTB/RIF assay and future prospects for a point-of-care test

Rapid progress has been made in the development of new diagnostic assays for tuberculosis in recent years. New technologies have been developed and assessed, and are now being implemented. The Xpert MTB/RIF assay, which enables simultaneous detection of Mycobacterium tuberculosis (MTB) and rifampicin (RIF) resistance, was endorsed by WHO in December, 2010. This assay was specifically recommended for use as the initial diagnostic test for suspected drug-resistant or HIV-associated pulmonary tuberculosis. By June, 2012, two-thirds of countries with a high tuberculosis burden and half of countries with a high multidrug-resistant tuberculosis burden had incorporated the assay into their national tuberculosis programme guidelines. Although the development of the Xpert MTB/RIF assay is undoubtedly a landmark event, clinical and programmatic effects and cost-effectiveness remain to be defined. We review the rapidly growing body of scientific literature and discuss the advantages and challenges of using the Xpert MTB/RIF assay in areas where tuberculosis is endemic. We also review other prospects within the developmental pipeline. A rapid, accurate point-of-care diagnostic test that is affordable and can be readily implemented is urgently needed. Investment in the tuberculosis diagnostics pipeline should remain a major priority for funders and researchers.