Genome analysis of multi- and extensively-drug-resistant tuberculosis from KwaZulu-Natal, South Africa

Pan-genome de Bruijn graph using the bidirectional FM-index

BACKGROUND

Pan-genome graphs are gaining importance in the field of bioinformatics as data structures to represent and jointly analyze multiple genomes. Compacted de Bruijn graphs are inherently suited for this purpose, as their graph topology naturally reveals similarity and divergence within the pan-genome. Most state-of-the-art pan-genome graphs are represented explicitly in terms of nodes and edges. Recently, an alternative, implicit graph representation was proposed that builds directly upon the unidirectional FM-index. As such, a memory-efficient graph data structure is obtained that inherits the FM-index' backward search functionality. However, this representation suffers from a number of shortcomings in terms of functionality and algorithmic performance.

RESULTS

We present a data structure for a pan-genome, compacted de Bruijn graph that aims to address these shortcomings. It is built on the bidirectional FM-index, extending the ability of its unidirectional counterpart to navigate and search the graph in both directions. All basic graph navigation steps can be performed in constant time. Based on these features, we implement subgraph visualization as well as lossless approximate pattern matching to the graph using search schemes. We demonstrate that we can retrieve all occurrences corresponding to a read within a certain edit distance in a very efficient manner. Through a case study, we show the potential of exploiting the information embedded in the graph's topology through visualization and sequence alignment.

CONCLUSIONS

We propose a memory-efficient representation of the pan-genome graph that supports subgraph visualization and lossless approximate pattern matching of reads against the graph using search schemes. The C++ source code of our software, called Nexus, is available at https://github.com/biointec/nexus under AGPL-3.0 license.

Combining Radiological and Genomic TB Portals Data for Drug Resistance Analysis

Tuberculosis (TB) drug resistance is a worldwide public health problem. It decreases the likelihood of a positive outcome for the individual patient and increases the likelihood of disease spread. Therefore, early detection of TB drug resistance is crucial for improving outcomes and controlling disease transmission. While drug-sensitive tuberculosis cases are declining worldwide because of effective treatment, the threat of drug-resistant tuberculosis is growing, and the success rate of drug-resistant tuberculosis treatment is only around 60%. The TB Portals program provides a publicly accessible repository of TB case data with an emphasis on collecting drug-resistant cases. The dataset includes multi-modal information such as socioeconomic/geographic data, clinical characteristics, pathogen genomics, and radiological features. The program is an international collaboration whose participants are typically under a substantial burden of drug-resistant tuberculosis, with data collected from standard clinical care provided to the patients. Consequentially, the TB Portals dataset is heterogenous in nature, with data representing multiple treatment centers in different countries and containing cross-domain information. This study presents the challenges and methods used to address them when working with this real-world dataset. Our goal was to evaluate whether combining radiological features derived from a chest X-ray of the host and genomic features from the pathogen can potentially improve the identification of the drug susceptibility type, drug-sensitive (DS-TB) or drug-resistant (DR-TB), and the length of the first successful drug regimen. To perform these studies, significantly imbalanced data needed to be processed, which included a much larger number of DR-TB cases than DS-TB, many more cases with radiological findings than genomic ones, and the sparse high dimensional nature of the genomic information. Three evaluation studies were carried out. First, the DR-TB/DS-TB classification model achieved an average accuracy of 92.4% when using genomic features alone or when combining radiological and genomic features. Second, the regression model for the length of the first successful treatment had a relative error of 53.5% using radiological features, 25.6% using genomic features, and 22.0% using both radiological and genomic features. Finally, the relative error of the third regression model predicting the length of the first treatment using the most common drug combination varied depending on the feature type used. When using radiological features alone, the relative error was 17.8%. For genomic features alone, the relative error increased to 19.9%. The model had a relative error of 19.0% when both radiological and genomic features were combined. Although combining radiological and genomic features did not improve upon the use of genomic features when classifying DR-TB/DS-TB, the combination of the two feature types improved the relative error of the predictive model for the length of the first successful treatment. Furthermore, the regression model trained on radiological features achieved the best performance when predicting the treatment length of the most common drug combination.

Structure-directed identification of pyridine-2-methylamine derivatives as MmpL3 inhibitors for use as antitubercular agents

Mycobacterial membrane protein Large 3 (MmpL3), an inner membrane protein, plays a crucial role in the transport of mycolic acids that are essential for the viability of M. tuberculosis and has been a promising therapeutic target for new anti-TB agents. Herein, we report the discovery of pyridine-2-methylamine antitubercular compounds using a structure-based drug design strategy. Compound 62 stands out as the most potent compound with high activity against M. tb strain H37Rv (MIC = 0.016 μg/mL) as well as the clinically isolated strains of MDR/XDR-TB (MIC = 0.0039-0.0625 μg/mL), low Vero cell toxicity (IC50 ≥ 16 μg/mL), and moderate liver microsomal stability (CLint = 28 μL/min/mg). Furthermore, the resistant mutant of S288T due to single nucleotide polymorphism in mmpL3 was resistant to pyridine-2-methylamine 62, demonstrating compound 62 is likely target to MmpL3.

Whole-Genome Sequence-Based Characterization of Pre-XDR M. tuberculosis Clinical Isolates Collected in Kazakhstan

BACKGROUND

Kazakhstan has a high burden of multidrug-resistant tuberculosis in the Central Asian region. This study aimed to perform genomic characterization of Mycobacterium tuberculosis strains obtained from Kazakhstani patients with pre-extensively drug-resistant tuberculosis diagnosed in Kazakhstan.

METHODS

Whole-genome sequencing was performed on 10 pre-extensively drug-resistant M. tuberculosis strains from different regions of Kazakhstan. All strains had high-confidence resistance mutations according to the resistance grading system previously established by the World Health Organization. The genome analysis was performed using TB-Profiler, Mykrobe, CASTB, and ResFinder.

RESULTS

Valuable information for understanding the genetic diversity of tuberculosis in Kazakhstan can also be obtained from whole-genome sequencing. The results from the Phenotypic Drug Susceptibility Testing (DST) of bacterial strains were found to be consistent with the drug resistance information obtained from genomic data that characterized all isolates as pre-XDR. This information can help in developing targeted prevention and control strategies based on the local epidemiology of tuberculosis. Furthermore, the data obtained from whole-genome sequencing can help in tracing the transmission pathways of tuberculosis and facilitating early detection of outbreaks.

CONCLUSIONS

The results from whole-genome sequencing of tuberculosis clinical samples in Kazakhstan provide important insights into the drug resistance patterns and genetic diversity of tuberculosis in the country. These results can contribute to the improvement of tuberculosis control and management programs in Kazakhstan.

Antibiotic resistance acquisition versus primary transmission in the presentation of extensively drug-resistant tuberculosis

Mycobacterium tuberculosis is the leading cause of mortality worldwide due to a single bacterial pathogen. Of concern is the negative impact that the COVID-19 pandemic has had on the control of tuberculosis (TB) including drug-resistant forms of the disease. Antimicrobial resistance increases the likelihood of worsened outcomes in TB patients including treatment failure and death. Multidrug-resistant (MDR) strains, resistant to first-line drugs isoniazid and rifampin, and extensively drug-resistant (XDR) strains with further resistance to second-line drugs (SLD), threaten control programs designed to lower TB incidence and end the disease as a public health challenge by 2030, in accordance with UN Sustainable Development Goals. Tackling TB requires an understanding of the pathways through which drug resistance emerges. Here, the roles of acquired resistance mutation, and primary transmission, are examined with regard to XDR-TB. It is apparent that XDR-TB can emerge from MDR-TB through a small number of additional resistance mutations that occur in patients undergoing drug treatment. Rapid detection of resistance, to first-line drugs and SLD, at the initiation of and during treatment, and prompt adjustment of regimens are required to ensure treatment success in these patients. Primary transmission is predicted to make an increasing contribution to the XDR-TB caseload in the future. Much work is required to improve the implementation of the World Health Organization-recommended infection control practices and block onward transmission of XDR-TB patients to contacts including health-care workers. Finally, limiting background resistance to fluoroquinolones in pre-XDR strains of M. tuberculosis will necessitate better antimicrobial stewardship in the broader use of this drug class.

High clustering rate and genotypic drug-susceptibility screening for the newly recommended anti-tuberculosis drugs among global extensively drug-resistant Mycobacterium tuberculosis isolates

Multidrug-resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB) make TB difficult to control. Global susceptibility data for six newly recommended anti-TB drugs against M/XDR-TB are still limited. Using publicly available whole-genome sequences, we determined the proportion of 513 phenotypically XDR-TB isolates that carried mutations associated with resistance against these drugs (bedaquiline, clofazimine, linezolid, delamanid, pretomanid and cycloserine). Mutations of Rv0678 and Rv1979c were detected in 69/513 isolates (13.5%) for bedaquiline resistance and 79/513 isolates (15.4%) for clofazimine resistance with additional mmpL5 mutations. Mutations conferring resistance to delamanid were detected in fbiB and ddn genes for 11/513 isolates (2.1%). For pretomanid, a mutation was detected in the ddn gene for 3/513 isolates (0.6%). Nineteen mutations of pykA, cycA, ald, and alr genes, conferring resistance to cycloserine, were found in 153/513 isolates (29.8%). No known mutations associated with linezolid resistance were detected. Cluster analysis showed that 408/513 isolates fell within 99 clusters and that 354 of these isolates were possible primary drug-resistant TB (292 XDR-TB, 57 pre-XDR-TB and 5 MDR-TB). Clonal transmission of primary XDR isolates might contribute significantly to the high prevalence of DR-TB globally.

Clinical and pathological characteristics associated with the presence of the IS6110 Mycobacterim tuberculosis transposon in neoplastic cells from non-small cell lung cancer patients

Lung cancer (LC) and pulmonary tuberculosis (TB) are the deadliest neoplastic and bacterial infectious diseases worldwide, respectively. Clinicians and pathologists have long discussed the co-existence of LC and TB, and several epidemiologic studies have presented evidence indicating that TB could be associated with the development of LC, particularly adenocarcinoma. Nonetheless, this data remains controversial, and the mechanism which could underlie the association remains largely unexplored. Some bioinformatic studies have shown that human cancer biopsies have a very high frequency of bacterial DNA integration; since Mycobacterium Tuberculosis (MTb) is an intracellular pathogen, it could play an active role in the cellular transformation. Our group performed an exploratory study in a cohort of 88 LC patients treated at the Instituto Nacional de Cancelorogía (INCan) of Mexico City to evaluate the presence of MTb DNA in LC tissue specimens. For the first time, our results show the presence of the MTb IS6110 transposon in 40.9% (n = 36/88) of patients with lung adenocarcinomas. Additionally, through in-situ PCR we identified the presence of IS6110 in the nuclei of tumor cells. Furthermore, shotgun sequencing from two samples identified traces of MTb genomes present in tumor tissue, suggesting that similar Mtb strains could be infecting both patients.

Genomic Analysis of Multidrug-Resistant Mycobacterium tuberculosis Strains From Patients in Kazakhstan

Tuberculosis (TB) is an infectious disease that remains an essential public health problem in many countries. Despite decreasing numbers of new cases worldwide, the incidence of antibiotic-resistant forms (multidrug resistant and extensively drug-resistant) of TB is increasing. Next-generation sequencing technologies provide a high-throughput approach to identify known and novel potential genetic variants that are associated with drug resistance in Mycobacterium tuberculosis (Mtb). There are limited reports and data related to whole-genome characteristics of drug-resistant Mtb strains circulating in Kazakhstan. Here, we report whole-genome sequencing and analysis results of eight multidrug-resistant strains collected from TB patients in Kazakhstan. Genotyping and validation of all strains by MIRU-VNTR and spoligotyping methodologies revealed that these strains belong to the Beijing family. The spectrum of specific and potentially novel genomic variants (single-nucleotide polymorphisms, insertions, and deletions) related to drug resistance was identified and annotated. ResFinder, CARD, and CASTB antibiotic resistance databases were used for the characterization of genetic variants in genes associated with drug resistance. Our results provide reference data and genomic profiles of multidrug-resistant isolates for further comparative studies and investigations of genetic patterns in drug-resistant Mtb strains.

A fragment-based approach to assess the ligandability of ArgB, ArgC, ArgD and ArgF in the L-arginine biosynthetic pathway of Mycobacterium tuberculosis

The L-arginine biosynthesis pathway consists of eight enzymes that catalyse the conversion of L-glutamate to L-arginine. Arginine auxotrophs (argB/argF deletion mutants) of Mycobacterium tuberculosis are rapidly sterilised in mice, while inhibition of ArgJ with Pranlukast was found to clear chronic M. tuberculosis infection in a mouse model. Enzymes in the arginine biosynthetic pathway have therefore emerged as promising targets for anti-tuberculosis drug discovery. In this work, the ligandability of four enzymes of the pathway ArgB, ArgC, ArgD and ArgF is assessed using a fragment-based approach. We identify several hits against these enzymes validated with biochemical and biophysical assays, as well as X-ray crystallographic data, which in the case of ArgB were further confirmed to have on-target activity against M. tuberculosis. These results demonstrate the potential for more enzymes in this pathway to be targeted with dedicated drug discovery programmes.

Distribution and Clonality of drug-resistant tuberculosis in South Africa

Background

Studies have shown that drug-resistant tuberculosis (DR-TB) in South Africa (SA) is clonal and is caused mostly by transmission. Identifying transmission chains is important in controlling DR-TB. This study reports on the sentinel molecular surveillance data of Rifampicin-Resistant (RR) TB in SA, aiming to describe the RR-TB strain population and the estimated transmission of RR-TB cases.

Method

RR-TB isolates collected between 2014 and 2018 from eight provinces were genotyped using combination of spoligotyping and 24-loci mycobacterial interspersed repetitive-units-variable-number tandem repeats (MIRU-VNTR) typing.

Results

Of the 3007 isolates genotyped, 301 clusters were identified. Cluster size ranged between 2 and 270 cases. Most of the clusters (247/301; 82.0%) were small in size (< 5 cases), 12.0% (37/301) were medium sized (5–10 cases), 3.3% (10/301) were large (11–25 cases) and 2.3% (7/301) were very large with 26–270 cases. The Beijing genotype was responsible for majority of RR-TB cases in Western and Eastern Cape, while the East-African-Indian-Somalian (EAI1_SOM) genotype accounted for a third of RR-TB cases in Mpumalanga. The overall proportion of RR-TB cases estimated to be due to transmission was 42%, with the highest transmission-rate in Western Cape (64%) and the lowest in Northern Cape (9%).

Conclusion

Large clusters contribute to the burden of RR-TB in specific geographic areas such as Western Cape, Eastern Cape and Mpumalanga, highlighting the need for community-wide interventions. Most of the clusters identified in the study were small, suggesting close contact transmission events, emphasizing the importance of contact investigations and infection control as the primary interventions in SA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02232-z.

The promises and limitations of N-acetylcysteine as a potentiator of first-line and second-line tuberculosis drugs

N-acetylcysteine (NAC) is most commonly used for the treatment of acetaminophen overdose and acetaminophen-induced liver injury. In patients infected with Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), NAC is given to treat hepatotoxicity induced by TB drugs. We had previously shown that cysteine, a derivative of NAC, potentiated the activity of isoniazid, a first-line TB drug, by preventing the emergence of INH resistance and persistence in M. tuberculosis in vitro. Herein, we demonstrate that in vitro, NAC has the same boosting activity with various combinations of first- and second-line TB drugs against drug-susceptible and multidrug-resistant M. tuberculosis strains. Similar to cysteine, NAC increased M. tuberculosis respiration. However, in M. tuberculosis-infected mice, the addition of NAC did not augment the activity of first- or second-line TB drugs. A comparison of the activity of NAC combined with TB drugs in murine and human macrophage cell lines revealed that studies in mice might not be recapitulated during host infection in vivo.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Exact mapping of Illumina blind spots in the Mycobacterium tuberculosis genome reveals platform-wide and workflow-specific biases

Whole-genome sequencing (WGS) is fundamental to Mycobacterium tuberculosis basic research and many clinical applications. Coverage across Illumina-sequenced M. tuberculosis genomes is known to vary with sequence context, but this bias is poorly characterized. Here, through a novel application of phylogenomics that distinguishes genuine coverage bias from deletions, we discern Illumina ‘blind spots’ in the M. tuberculosis reference genome for seven sequencing workflows. We find blind spots to be widespread, affecting 529 genes, and provide their exact coordinates, enabling salvage of unaffected regions. Fifty-seven pe/ppe genes (the primary families assumed to exhibit Illumina bias) lack blind spots entirely, while the remaining pe/ppe genes account for 55.1 % of blind spots. Surprisingly, we find coverage bias persists in homopolymers as short as 6 bp, shorter tracts than previously reported. While G+C-rich regions challenge all Illumina sequencing workflows, a modified Nextera library preparation that amplifies DNA with a high-fidelity polymerase markedly attenuates coverage bias in G+C-rich and homopolymeric sequences, expanding the ‘Illumina-sequenceable’ genome. Through these findings, and by defining workflow-specific exclusion criteria, we spotlight effective strategies for handling bias in M. tuberculosis Illumina WGS. This empirical analysis framework may be used to systematically evaluate coverage bias in other species using existing sequencing data.

Higher genome mutation rates of Beijing lineage of Mycobacterium tuberculosis during human infection

Mycobacterium tuberculosis (Mtb) strains of Beijing lineage have caused great concern because of their rapid emergence of drug resistance and worldwide spread. DNA mutation rates that reflect evolutional adaptation to host responses and the appearance of drug resistance have not been elucidated in human-infected Beijing strains. We tracked and obtained an original Mtb isolate of Beijing lineage from the 1999 tuberculosis outbreak in Japan, as well as five other isolates that spread in humans, and two isolates from the patient caused recurrence. Three isolates were from patients who developed TB within one year after infection (rapid-progressor, RP), and the other three isolates were from those who developed TB more than one year after infection (slow-progressor, SP). We sequenced genomes of these isolates and analyzed the propensity and rate of genomic mutations. Generation time versus mutation rate curves were significantly higher for RP. The ratio of oxidative versus non-oxidation damages induced mutations was higher in SP than RP, suggesting that persistent Mtb are exposed to oxidative stress in the latent state. Our data thus demonstrates that higher mutation rates of Mtb Beijing strains during human infection is likely to account for the higher adaptability and an emergence ratio of drug resistance.

Molecular epidemiology of drug resistant Mycobacterium tuberculosis in Africa: a systematic review

Background

The burden of drug resistant tuberculosis in Africa is largely driven by the emergence and spread of multidrug resistant (MDR) and extensively drug resistant (XDR) Mycobacterium tuberculosis strains. MDR-TB is defined as resistance to isoniazid and rifampicin, while XDR-TB is defined as MDR-TB with added resistance to any of the second line injectable drugs and any fluoroquinolone.

The highest burden of drug resistant TB is seen in countries further experiencing an HIV epidemic. The molecular mechanisms of drug resistance as well as the evolution of drug resistant TB strains have been widely studied using various genotyping tools. The study aimed to analyse the drug resistant lineages in circulation and transmission dynamics of these lineages in Africa by describing outbreaks, nosocomial transmission and migration. Viewed as a whole, this can give a better insight into the transmission dynamics of drug resistant TB in Africa.

Methods

A systematic review was performed on peer reviewed original research extracted from PubMed reporting on the lineages associated with drug resistant TB from African countries, and their association with outbreaks, nosocomial transmission and migration. The search terms “Tuberculosis AND drug resistance AND Africa AND (spoligotyping OR molecular epidemiology OR IS6110 OR MIRU OR DNA fingerprinting OR RFLP OR VNTR OR WGS)” were used to identify relevant articles reporting the molecular epidemiology of drug resistant TB in Africa.

Results

Diverse genotypes are associated with drug resistant TB in Africa, with variations in strain predominance within the continent. Lineage 4 predominates across Africa demonstrating the ability of “modern strains” to adapt and spread easily. Most studies under review reported primary drug resistance as the predominant type of transmission. Drug resistant TB strains are associated with community and nosocomial outbreaks involving MDR- and XDR-TB strains. The under-use of molecular epidemiological tools is of concern, resulting in gaps in knowledge of the transmission dynamics of drug resistant TB on the continent.

Conclusions

Genetic diversity of M. tuberculosis strains has been demonstrated across Africa implying that diverse genotypes are driving the epidemiology of drug resistant TB across the continent.

Rapid Pyrazinamide Drug Susceptibility Testing using a Closed-Tube PCR Assay of the Entire pncA gene

The continued use of pyrazinamide in the treatment of tuberculosis in the absence of a rapid, accurate and standardized pyrazinamide drug susceptibility assays is of great concern. While whole genome sequencing holds promise, it is not yet feasible option in low resource settings as it requires expensive instruments and bioinformatic analysis. We investigated the diagnostic performance of a closed-tube Linear-After-The-Exponential (LATE)-PCR assay for pyrazinamide susceptibility in Mycobacterium tuberculosis. Based on a set of 654 clinical Mycobacterium tuberculosis culture isolates with known mutations throughout the pncA gene as determined by Sanger sequencing, the assay displays excellent sensitivity of 96.9% (95% CI: 95.2-98.6) and specificity of 97.9% (95% CI: 96.1-99.7). In a subset of 384 isolates with phenotypic drug susceptibility testing, we also observed high sensitivity of 98.9% (95% CI: 97.5-100) but lower specificity of 91.8% (95% CI: 87.9-95.8) when compared to phenotypic drug susceptibility testing. We conclude that the LATE PCR assay offers both a rapid and accurate prediction of pyrazinamide susceptibility.

Pre-detection history of extensively drug-resistant tuberculosis in KwaZulu-Natal, South Africa

Antimicrobial-resistant (AMR) infections pose a major threat to global public health. Similar to other AMR pathogens, both historical and ongoing drug-resistant tuberculosis (TB) epidemics are characterized by transmission of a limited number of predominant Mycobacterium tuberculosis (Mtb) strains. Understanding how these predominant strains achieve sustained transmission, particularly during the critical period before they are detected via clinical or public health surveillance, can inform strategies for prevention and containment. In this study, we employ whole-genome sequence (WGS) data from TB clinical isolates collected in KwaZulu-Natal, South Africa to examine the pre-detection history of a successful strain of extensively drug-resistant (XDR) TB known as LAM4/KZN, first identified in a widely reported cluster of cases in 2005. We identify marked expansion of this strain concurrent with the onset of the generalized HIV epidemic 12 y prior to 2005, localize its geographic origin to a location in northeastern KwaZulu-Natal ∼400 km away from the site of the 2005 outbreak, and use protein structural modeling to propose a mechanism for how strain-specific rpoB mutations offset fitness costs associated with rifampin resistance in LAM4/KZN. Our findings highlight the importance of HIV coinfection, high preexisting rates of drug-resistant TB, human migration, and pathoadaptive evolution in the emergence and dispersal of this critical public health threat. We propose that integrating whole-genome sequencing into routine public health surveillance can enable the early detection and local containment of AMR pathogens before they achieve widespread dispersal.

Bridging the TB data gap: in silico extraction of rifampicin-resistant tuberculosis diagnostic test results from whole genome sequence data

Background

Mycobacterium tuberculosis rapid diagnostic tests (RDTs) are widely employed in routine laboratories and national surveys for detection of rifampicin-resistant (RR)-TB. However, as next-generation sequencing technologies have become more commonplace in research and surveillance programs, RDTs are being increasingly complemented by whole genome sequencing (WGS). While comparison between RDTs is difficult, all RDT results can be derived from WGS data. This can facilitate continuous analysis of RR-TB burden regardless of the data generation technology employed. By converting WGS to RDT results, we enable comparison of data with different formats and sources particularly for low- and middle-income high TB-burden countries that employ different diagnostic algorithms for drug resistance surveys. This allows national TB control programs (NTPs) and epidemiologists to utilize all available data in the setting for improved RR-TB surveillance.

Methods

We developed the Python-based MycTB Genome to Test (MTBGT) tool that transforms WGS-derived data into laboratory-validated results of the primary RDTs-Xpert MTB/RIF, XpertMTB/RIF Ultra, GenoType MDRTBplus v2.0, and GenoscholarNTM+MDRTB II. The tool was validated through RDT results of RR-TB strains with diverse resistance patterns and geographic origins and applied on routine-derived WGS data.

Results

The MTBGT tool correctly transformed the single nucleotide polymorphism (SNP) data into the RDT results and generated tabulated frequencies of the RDT probes as well as rifampicin-susceptible cases. The tool supplemented the RDT probe reactions output with the RR-conferring mutation based on identified SNPs. The MTBGT tool facilitated continuous analysis of RR-TB and Xpert probe reactions from different platforms and collection periods in Rwanda.

Conclusion

Overall, the MTBGT tool allows low- and middle-income countries to make sense of the increasingly generated WGS in light of the readily available RDT results, and assess whether currently implemented RDTs adequately detect RR-TB in their setting. With its feature to transform WGS to RDT results and facilitate continuous RR-TB data analysis, the MTBGT tool may bridge the gap between and among data from periodic surveys, continuous surveillance, research, and routine tests, and may be integrated within the national information system for use by the NTP and epidemiologists to improve setting-specific RR-TB control. The MTBGT source code and accompanying documentation are available at https://github.com/KamelaNg/MTBGT.

Mycobacterium tuberculosis Next-Generation Whole Genome Sequencing: Opportunities and Challenges

Mycobacterium tuberculosis drug resistance is a threat to global tuberculosis (TB) control. Comprehensive and timely drug susceptibility determination is critical to inform appropriate treatment of drug-resistant tuberculosis (DR-TB). Phenotypic drug susceptibility testing (DST) is the gold standard for M. tuberculosis drug resistance determination. M. tuberculosis whole genome sequencing (WGS) has the potential to be a one-stop method for both comprehensive DST and epidemiological investigations. We discuss in this review the tremendous opportunities that next-generation WGS presents in terms of understanding the molecular epidemiology of tuberculosis and mechanisms of drug resistance. The potential clinical value and public health impact in the areas of DST for patient management and tracing of transmission chains for timely public health intervention are also discussed. We present the current challenges for the implementation of WGS in low and middle-income settings. WGS analysis has already been adapted routinely in laboratories to inform patient management and public health interventions in low burden high-income settings such as the United Kingdom. We predict that the technology will be adapted similarly in high burden settings where the impact on the epidemic will be greatest.

Approximate, simultaneous comparison of microbial genome architectures via syntenic anchoring of quiver representations

Motivation

A long-standing limitation in comparative genomic studies is the dependency on a reference genome, which hinders the spectrum of genetic diversity that can be identified across a population of organisms. This is especially true in the microbial world where genome architectures can significantly vary. There is therefore a need for computational methods that can simultaneously analyze the architectures of multiple genomes without introducing bias from a reference.

Results

In this article, we present Ptolemy: a novel method for studying the diversity of genome architectures-such as structural variation and pan-genomes-across a collection of microbial assemblies without the need of a reference. Ptolemy is a 'top-down' approach to compare whole genome assemblies. Genomes are represented as labeled multi-directed graphs-known as quivers-which are then merged into a single, canonical quiver by identifying 'gene anchors' via synteny analysis. The canonical quiver represents an approximate, structural alignment of all genomes in a given collection encoding structural variation across (sub-) populations within the collection. We highlight various applications of Ptolemy by analyzing structural variation and the pan-genomes of different datasets composing of Mycobacterium, Saccharomyces, Escherichia and Shigella species. Our results show that Ptolemy is flexible and can handle both conserved and highly dynamic genome architectures. Ptolemy is user-friendly-requires only FASTA-formatted assembly along with a corresponding GFF-formatted file-and resource-friendly-can align 24 genomes in ∼10 mins with four CPUs and <2 GB of RAM.

Availability and implementation

Github: https://github.com/AbeelLab/ptolemy.

Supplementary information

Supplementary data are available at Bioinformatics online.

Advances in the development of molecular genetic tools for Mycobacterium tuberculosis

Mycobacterium tuberculosis, a clinically relevant Gram-positive bacterium of great clinical relevance, is a lethal pathogen owing to its complex physiological characteristics and development of drug resistance. Several molecular genetic tools have been developed in the past few decades to study this microorganism. These tools have been instrumental in understanding how M. tuberculosis became a successful pathogen. Advanced molecular genetic tools have played a significant role in exploring the complex pathways involved in M. tuberculosis pathogenesis. Here, we review various molecular genetic tools used in the study of M. tuberculosis. Further, we discuss the applications of clustered regularly interspaced short palindromic repeat interference (CRISPRi), a novel technology recently applied in M. tuberculosis research to study target gene functions. Finally, prospective outcomes of the applications of molecular techniques in the field of M. tuberculosis genetic research are also discussed.

The occurrence and frequency of genomic mutations that mediate Isoniazid and Rifampicin resistance in Mycobacterium tuberculosis isolates from untreated pulmonary Tuberculosis cases in urban Blantyre, Malawi

Background

The emergence and spread of drug-resistant Tuberculosis (TB) is a major public health threat. TB resistance originates in the course of treatment due to genomic mutations in Mycobacterium tuberculosis (MTB). An increase in new cases with drug-resistant TB could be an indicator of high levels of circulating resistant strains. This study was conducted to determine the occurrence and frequency of genomic mutations that mediate Isoniazid (INH) and Rifampicin (RIF) resistance among isolates from untreated TB cases in urban Blantyre, Malawi.

Methods

A cross-sectional retrospective study was conducted on a panel of 141(n=141) MTB clinical isolates recovered between June 2010 and January 2012 from >2+ Ziehl-Neelsen smear positive new pulmonary-TB patients with no history of treatment. Frozen isolates were revived using the BACTEC MGIT detection system. DNA was extracted using GenoLyse DNA extraction kit and detection of genomic mutations was carried out using the GenoType MTBDRplus Ver 2.0 assay.

Results

Out of the 141 isolates studied, 3 (2.1%) were found carrying mutations in the katG gene that confer resistance to Isoniazid (INH). No mutations were detected in the inhA promoter region gene that confer weak INH resistance or in the rpoB gene that confer Rifampicin resistance. All katG mutant genes had a S315T1 single point mutation, a genomic alteration that mediates high INH resistance.

Conclusion

The katG mutant gene conferring resistance to INH was the only genomic mutation observed among the isolates studied and the frequency of occurrence was low. Our findings suggest low levels of circulating drug-resistant MTB strains in urban Blantyre, Malawi.

Fitness-compensatory mutations facilitate the spread of drug-resistant F15/LAM4/KZN and F28 Mycobacterium tuberculosis strains in KwaZulu-Natal, South Africa

While the acquisition of drug resistance is often accompanied by fitness costs, Mycobacterium tuberculosis has developed mechanisms to overcome these costs in the form of compensatory mutations. In an attempt to dissect strain-specific differences in biological fitness, 10 M. tuberculosis genomes, representing F15/LAM4/KZN, Beijing, F11 and F28 genotypes were sequenced on the Illumina MiSeq platform. Drug-susceptible F15/LAM4/KZN strains differed by 43 SNPs, demonstrating that heterogeneity exists even among closely-related strains. We found unique, nonsynonymous single-nucleotide polymorphisms (SNPs) in the sigA and grcC1 genes of multidrug resistant (MDR) and XDR F15/LAM4/KZN strains, respectively. The F28 MDR strain harboured a novel ubiA mutation in combination with its embB M306I mutation, which may be related to ethambutol resistance. In addition, it possessed a low-frequency rpoC mutation, suggesting that this strain was in the process of developing compensation. In contrast, no compensatory mutations were identified in Beijing and F11 MDR strains, corroborating its low in vitro fitness. Clinical strains also harboured unique SNPs in a number of important genes associated with virulence, highlighting the need for future studies which examine the correlation of genetic variations with phenotypic diversity. In summary, whole-genome sequencing revealed the presence of fitness-compensatory mutations in F15/LAM4/KZN and F28 genotypes which predominate in MDR and/or extensively drug resistant (XDR) forms in KwaZulu-Natal, South Africa.

Whole genome sequencing data of 1110 Mycobacterium tuberculosis isolates identifies insertions and deletions associated with drug resistance

Background

Drug resistance in Mycobacterium tuberculosis (MTB) is one of the major challenges in tuberculosis (TB) treatment. However, known mutations cannot explain all of the cases of resistance and little research has focused on the relationship between insertions / deletions (indels) and drug resistance.

Results

Here, we retrieved whole genome sequencing data of 743 drug-resistant MTB strains and 367 pan-susceptible strains from TB patients from the public domain to identify novel genomic markers of drug resistance. A total of 20 region markers containing genes and intergenic regions (IGRs) with significant statistical correlation with antibiotic resistance were revealed, four of which have been previously reported to be associated with drug resistance. In addition, 83 point markers containing frameshift (FS) mutations and IGR indels were also identified independently based on differences in their incidence rates between drug-sensitive and -resistant strains. Among the 83 point markers, eight indels were detected in known drug-associated genes or IGRs. Furthermore, the overlap between 20 region markers and 83 point markers further indicated their associations with drug resistance. The markers identified were involved in essential bacterial metabolic functions, including cell wall and transmembrane transporter functions. A strong correlation between FS mutations and mutations in DNA repair genes including I21V in alkA, R48G in mutT4 and P2R in nth was also found.

Conclusions

This study identified a set of novel genetic markers with FS mutations and IGR indels associated with MTB drug resistance, which greatly broadens the pool of mutations related to MTB drug resistance. This insight may be important in identifying novel mechanisms of drug resistance in MTB.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4734-6) contains supplementary material, which is available to authorized users.

Evolution of drug resistance in Mycobacterium tuberculosis: a review on the molecular determinants of resistance and implications for personalized care

Drug-resistant TB (DR-TB) remains a significant challenge in TB treatment and control programmes worldwide. Advances in sequencing technology have significantly increased our understanding of the mechanisms of resistance to anti-TB drugs. This review provides an update on advances in our understanding of drug resistance mechanisms to new, existing drugs and repurposed agents. Recent advances in WGS technology hold promise as a tool for rapid diagnosis and clinical management of TB. Although the standard approach to WGS of Mycobacterium tuberculosis is slow due to the requirement for organism culture, recent attempts to sequence directly from clinical specimens have improved the potential to diagnose and detect resistance within days. The introduction of new databases may be helpful, such as the Relational Sequencing TB Data Platform, which contains a collection of whole-genome sequences highlighting key drug resistance mutations and clinical outcomes. Taken together, these advances will help devise better molecular diagnostics for more effective DR-TB management enabling personalized treatment, and will facilitate the development of new drugs aimed at improving outcomes of patients with this disease.

Mutations inside rifampicin-resistance determining region of rpoB gene associated with rifampicin-resistance in Mycobacterium tuberculosis

BACKGROUND

Rifampicin (RIF) plays a pivotal role in the treatment of tuberculosis due to its bactericidal effects. Because the action of RIF is on rpoB gene encoding RNA polymerase β subunit, 95% of RIF resistant mutations are present in rpoB gene. The majority of the mutations in rpoB gene are found within an 81bp RIF-resistance determining region (RRDR).

METHODOLOGY

Literatures on RIF resistant mutations published between 2010 and 2016 were thoroughly reviewed.

RESULTS

The most commonly mutated codons in RRDR of rpoB gene are 531, 526 and 516. The possibilities of absence of mutation in RRDR of rpoB gene in MDR-TB isolates in few studies was due to existence of other rare rpoB mutations outside RRDR or different mechanism of rifampicin resistance.

CONCLUSION

Molecular methods which can identify extensive mutations associated with multiple anti-tuberculous drugs are in urgent need so that the research on drug resistant mutations should be extended.

Multi- and Extensively Drug Resistant Mycobacterium tuberculosis in South Africa: a Molecular Analysis of Historical Isolates

Modern advances in genomics provide an opportunity to reinterpret historical bacterial culture collections. In this study, genotypic antibiotic resistance profiles of Mycobacterium tuberculosis isolates from a historical 20-year-old multidrug-resistant tuberculosis (MDR-TB) culture collection in South Africa are described. DNA samples extracted from the phenotypically MDR-TB isolates (n = 240) were assayed by Hain line probe assay (LPA) for the confirmation of MDR-TB and by Illumina Miseq whole-genome sequencing (WGS) for the characterization of mutations in eight genes (rpoB, katG, inhA, rpsL, pncA, embB, gyrA, and rrs) that are known to code for resistance to commonly used anti-TB agents. LPA identified 71.3% of the TB isolates as MDR-TB, 18.3% as rifampin (RIF) monoresistant, 2% as isoniazid (INH) monoresistant, and 8.3% as susceptible to both RIF and INH (RIF+INH). In a subset of 42 randomly selected isolates designated as RIF+INH resistant by Löwenstein-Jensen (LJ) culture in 1993, LPA and WGS results confirmed MDR-TB. In all five INH-monoresistant isolates by LPA and in all but one (the wild type) of the 34 successfully sequenced RIF-monoresistant isolates, WGS revealed matching mutations. Only 26% of isolates designated as susceptible by LPA, however, were found to be wild type by WGS. Novel mutations were found in the rpoB (Thr480Ala, Gln253Arg, Val249Met, Val251Tyr, Val251Phe), katG (Trp477STOP, Gln88STOP, Trp198STOP, Trp412STOP), embB (Thr11Xaa, Gln59Pro), and pncA (Thr100Ile, Thr159Ala, Ala134Arg, Val163Ala, Thr153Ile, DelGpos7, Phe106Ser) genes. Three MDR-TB isolates showed mutations in both the gyrA and rrs genes, suggesting that extensively drug-resistant tuberculosis existed in South Africa well before its formal recognition in 2006.

Potential Application of Digitally Linked Tuberculosis Diagnostics for Real-Time Surveillance of Drug-Resistant Tuberculosis Transmission: Validation and Analysis of Test Results

BACKGROUND

Tuberculosis (TB) is the highest-mortality infectious disease in the world and the main cause of death related to antimicrobial resistance, yet its surveillance is still paper-based. Rifampicin-resistant TB (RR-TB) is an urgent public health crisis. The World Health Organization has, since 2010, endorsed a series of rapid diagnostic tests (RDTs) that enable rapid detection of drug-resistant strains and produce large volumes of data. In parallel, most high-burden countries have adopted connectivity solutions that allow linking of diagnostics, real-time capture, and shared repository of these test results. However, these connected diagnostics and readily available test results are not used to their full capacity, as we have yet to capitalize on fully understanding the relationship between test results and specific rpoB mutations to elucidate its potential application to real-time surveillance.

OBJECTIVE

We aimed to validate and analyze RDT data in detail, and propose the potential use of connected diagnostics and associated test results for real-time evaluation of RR-TB transmission.

METHODS

We selected 107 RR-TB strains harboring 34 unique rpoB mutations, including 30 within the rifampicin resistance-determining region (RRDR), from the Belgian Coordinated Collections of Microorganisms, Antwerp, Belgium. We subjected these strains to Xpert MTB/RIF, GenoType MTBDRplus v2.0, and Genoscholar NTM + MDRTB II, the results of which were validated against the strains' available rpoB gene sequences. We determined the reproducibility of the results, analyzed and visualized the probe reactions, and proposed these for potential use in evaluating transmission.

RESULTS

The RDT probe reactions detected most RRDR mutations tested, although we found a few critical discrepancies between observed results and manufacturers' claims. Based on published frequencies of probe reactions and RRDR mutations, we found specific probe reactions with high potential use in transmission studies: Xpert MTB/RIF probes A, Bdelayed, C, and Edelayed; Genotype MTBDRplus v2.0 WT2, WT5, and WT6; and Genoscholar NTM + MDRTB II S1 and S3. Inspection of probe reactions of disputed mutations may potentially resolve discordance between genotypic and phenotypic test results.

CONCLUSIONS

We propose a novel approach for potential real-time detection of RR-TB transmission through fully using digitally linked TB diagnostics and shared repository of test results. To our knowledge, this is the first pragmatic and scalable work in response to the consensus of world-renowned TB experts in 2016 on the potential of diagnostic connectivity to accelerate efforts to eliminate TB. This is evidenced by the ability of our proposed approach to facilitate comparison of probe reactions between different RDTs used in the same setting. Integrating this proposed approach as a plug-in module to a connectivity platform will increase usefulness of connected TB diagnostics for RR-TB outbreak detection through real-time investigation of suspected RR-TB transmission cases based on epidemiologic linking.

The TB Portals: an Open-Access, Web-Based Platform for Global Drug-Resistant-Tuberculosis Data Sharing and Analysis

The TB Portals program is an international consortium of physicians, radiologists, and microbiologists from countries with a heavy burden of drug-resistant tuberculosis working with data scientists and information technology professionals. Together, we have built the TB Portals, a repository of socioeconomic/geographic, clinical, laboratory, radiological, and genomic data from patient cases of drug-resistant tuberculosis backed by shareable, physical samples. Currently, there are 1,299 total cases from five country sites (Azerbaijan, Belarus, Moldova, Georgia, and Romania), 976 (75.1%) of which are multidrug or extensively drug resistant and 38.2%, 51.9%, and 36.3% of which contain X-ray, computed tomography (CT) scan, and genomic data, respectively. The top Mycobacterium tuberculosis lineages represented among collected samples are Beijing, T1, and H3, and single nucleotide polymorphisms (SNPs) that confer resistance to isoniazid, rifampin, ofloxacin, and moxifloxacin occur the most frequently. These data and samples have promoted drug discovery efforts and research into genomics and quantitative image analysis to improve diagnostics while also serving as a valuable resource for researchers and clinical providers. The TB Portals database and associated projects are continually growing, and we invite new partners and collaborations to our initiative. The TB Portals data and their associated analytical and statistical tools are freely available at https://tbportals.niaid.nih.gov/.

Evolution of Mycobacterium tuberculosis: New Insights into Pathogenicity and Drug Resistance

The tuberculosis agent Mycobacterium tuberculosis has undergone a long and selective evolution toward human infection and represents one of the most widely spread pathogens due to its efficient aerosol-mediated human-to-human transmission. With the availability of more and more genome sequences, the evolutionary trajectory of this obligate pathogen becomes visible, which provides us with new insights into the molecular events governing evolution of the bacterium and its ability to accumulate drug-resistance mutations. In this review, we summarize recent developments in mycobacterial research related to this matter that are important for a better understanding of the current situation and future trends and developments in the global epidemiology of tuberculosis, as well as for possible public health intervention possibilities.

Revisiting the mutant prevention concentration to guide dosing in childhood tuberculosis

The mutant prevention concentration (MPC) is a well-known concept in the chemotherapy of many bacterial infections, but is seldom considered in relation to tuberculosis (TB) treatment, as the required concentrations are generally viewed as unachievable without undue toxicity. Early studies revealed single mutations conferring high MICs of first- and second-line anti-TB agents; however, the growing application of genomics and quantitative drug susceptibility testing in TB suggests a wide range of MICs often determined by specific mutations and strain type. In paediatric TB, pharmacokinetic studies indicate that despite increasing dose recommendations, a proportion of children still do not achieve adult-derived targets. When considering the next stage in anti-TB drug dosing and the introduction of novel therapies for children, we suggest consideration of MPC and its incorporation into pharmacokinetic studies to more accurately determine appropriate concentration targets in children, to restrict the growth of resistant mutants and better manage drug-resistant TB.

Direct Detection of Rifampin and Isoniazid Resistance in Sputum Samples from Tuberculosis Patients by High-Resolution Melt Curve Analysis

Drug-resistant tuberculosis (TB) is a major threat to TB control worldwide. Globally, only 40% of the 340,000 notified TB patients estimated to have multidrug-resistant-TB (MDR-TB) were detected in 2015. This study was carried out to evaluate the utility of high-resolution melt curve analysis (HRM) for the rapid and direct detection of MDR-TB in Mycobacterium tuberculosis in sputum samples. A reference plasmid library was first generated of the most frequently observed mutations in the resistance-determining regions of rpoB, katG, and an inhA promoter and used as positive controls in HRM. The assay was first validated in 25 MDR M. tuberculosis clinical isolates. The assay was evaluated on DNA isolated from 99 M. tuberculosis culture-positive sputum samples that included 84 smear-negative sputum samples, using DNA sequencing as gold standard. Mutants were discriminated from the wild type by comparing melting-curve patterns with those of control plasmids using HRM software. Rifampin (RIF) and isoniazid (INH) monoresistance were detected in 11 and 21 specimens, respectively, by HRM. Six samples were classified as MDR-TB by sequencing, one of which was missed by HRM. The HRM-RIF, INH-katG, and INH-inhA assays had 89% (95% confidence interval [CI], 52, 100%), 85% (95% CI, 62, 97%), and 100% (95% CI, 74, 100%) sensitivity, respectively, in smear-negative samples, while all assays had 100% sensitivity in smear-positive samples. All assays had 100% specificity. Concordance of 97% to 100% (κ value, 0.9 to 1) was noted between sequencing and HRM. Heteroresistance was observed in 5 of 99 samples by sequencing. In conclusion, the HRM assay was a cost-effective (Indian rupee [INR]400/US$6), rapid, and closed-tube method for the direct detection of MDR-TB in sputum, especially for direct smear-negative cases.

Transcriptomics as a tool to discover new antibacterial targets

The emergence of antibiotic-resistant pathogens, multiple drug-resistance, and extremely drug-resistant strains demonstrates the need for improved strategies to discover new drug-based compounds. The development of transcriptomics, proteomics, and metabolomics has provided new tools for global studies of living organisms. However, the compendium of expression profiles produced by these methods has introduced new scientific challenges into antimicrobial research. In this review, we discuss the practical value of transcriptomic techniques as well as their difficulties and pitfalls. We advocate the construction of new databases of transcriptomic data, using standardized formats in addition to standardized models of bacterial and yeast similar to those used in systems biology. The inclusion of proteomic and metabolomic data is also essential, as the resulting networks can provide a landscape to rationally predict and exploit new drug targets and to understand drug synergies.

Genomic epidemiology of Lineage 4 Mycobacterium tuberculosis subpopulations in New York City and New Jersey, 1999-2009

Background

Whole genome sequencing (WGS) has rapidly become an important research tool in tuberculosis epidemiology and is likely to replace many existing methods in public health microbiology in the near future. WGS-based methods may be particularly useful in areas with less diverse Mycobacterium tuberculosis populations, such as New York City, where conventional genotyping is often uninformative and field epidemiology often difficult. This study applies four candidate strategies for WGS-based identification of emerging M. tuberculosis subpopulations, employing both phylogenomic and population genetics methods.

Results

M. tuberculosis subpopulations in New York City and New Jersey can be distinguished via phylogenomic reconstruction, evidence of demographic expansion and subpopulation-specific signatures of selection, and by determination of subgroup-defining nucleotide substitutions. These methods identified known historical outbreak clusters and previously unidentified subpopulations within relatively monomorphic M. tuberculosis endemic clone groups. Neutrality statistics based on the site frequency spectrum were less useful for identifying M. tuberculosis subpopulations, likely due to the low levels of informative genetic variation in recently diverged isolate groups. In addition, we observed that isolates from New York City endemic clone groups have acquired multiple non-synonymous SNPs in virulence- and growth-associated pathways, and relatively few mutations in drug resistance-associated genes, suggesting that overall pathoadaptive fitness, rather than the acquisition of drug resistance mutations, has played a central role in the evolutionary history and epidemiology of M. tuberculosis subpopulations in New York City.

Conclusions

Our results demonstrate that some but not all WGS-based methods are useful for detection of emerging M. tuberculosis clone groups, and support the use of phylogenomic reconstruction in routine tuberculosis laboratory surveillance, particularly in areas with relatively less diverse M. tuberculosis populations. Our study also supports the use of wider-reaching phylogenomic and population genomic methods in tuberculosis public health practice, which can support tuberculosis control activities by identifying genetic polymorphisms contributing to epidemiological success in local M. tuberculosis populations and possibly explain why certain isolate groups are apparently more successful in specific host populations.

Electronic supplementary material

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

Mechanisms of first-line antimicrobial resistance in multi-drug and extensively drug resistant strains of Mycobacterium tuberculosis in KwaZulu-Natal, South Africa

BACKGROUND

In South Africa, drug resistant tuberculosis is a major public health crisis in the face of the colossal HIV pandemic.

METHODS

In an attempt to understand the distribution of drug resistance in our setting, we analysed the rpoB, katG, inhA, pncA and embB genes associated with resistance to key drugs used in the treatment of tuberculosis in clinical isolates of Mycobacterium tuberculosis in the KwaZulu-Natal province.

RESULTS

Classical mutations were detected in the katG, inhA and embB genes associated with resistance to isoniazid and ethambutol. Diverse mutations were recorded in the multidrug resistant (MDR) and extensively drug resistant (XDR) isolates for the rpoB and pncA gene associated with resistance to rifampicin and pyrazinamide.

CONCLUSIONS

M.tuberculosis strains circulating in our setting display a combination of previously observed mutations, each mediating resistance to a different drug. The MDR and XDR TB isolates analysed in this study displayed classical mutations linked to INH and EMB resistance, whilst diverse mutations were linked to RIF and PZA resistance. The similarity of the XDR strains confirms reports of the clonality of the XDR epidemic. The successful dissemination of the drug resistant strains in the province underscores the need for rapid diagnostics to effectively diagnose drug resistance and guide treatment.

Whole genome sequencing of Mycobacterium tuberculosis for detection of drug resistance: a systematic review

OBJECTIVES

We conducted a systematic review to determine the diagnostic accuracy of whole genome sequencing (WGS) of Mycobacterium tuberculosis for the detection of resistance to first- and second-line anti-tuberculosis (TB) drugs.

METHODS

The study was conducted according to the criteria of the Preferred Reporting Items for Systematic Reviews group. A total of 20 publications were included. The sensitivity, specificity, positive-predictive value and negative-predictive value of WGS using phenotypic drug susceptibility testing methods as a reference standard were determined.

RESULTS

Anti-TB agents tested included all first-line drugs, a variety of reserve drugs, as well as new drugs. Polymorphisms in a total of 53 genes were tested for associations with drug resistance. Pooled sensitivity and specificity values for detection of resistance to selected first-line drugs were 0.98 (95% CI 0.93-0.98) and 0.98 (95% CI 0.98-1.00) for rifampicin and 0.97 (95% CI 0.94-0.99) and 0.93 (95% CI 0.91-0.96) for isoniazid, respectively. Due to high heterogeneity in study designs, lack of data, knowledge of resistance mechanisms and clarity on exclusion of phylogenetic markers, there was a significant variation in analytical performance of WGS for the remaining first-line, reserved drugs and new drugs.

CONCLUSIONS

Whole genome sequencing could be considered a promising alternative to existing phenotypic and molecular drug susceptibility testing methods for rifampicin and isoniazid pending standardization of analytical pipelines. To ensure clinical relevance of WGS for detection of M. tuberculosis complex drug resistance, future studies should include information on clinical outcomes.

Recent developments in genomics, bioinformatics and drug discovery to combat emerging drug-resistant tuberculosis

Emergence of drug-resistant tuberculosis (DR-TB) is a big challenge in TB control. The delay in diagnosis of DR-TB leads to its increased transmission, and therefore prevalence. Recent developments in genomics have enabled whole genome sequencing (WGS) of Mycobacterium tuberculosis (M. tuberculosis) from 3-day-old liquid culture and directly from uncultured sputa, while new bioinformatics tools facilitate to determine DR mutations rapidly from the resulting sequences. The present drug discovery and development pipeline is filled with candidate drugs which have shown efficacy against DR-TB. Furthermore, some of the FDA-approved drugs are being evaluated for repurposing, and this approach appears promising as several drugs are reported to enhance efficacy of the standard TB drugs, reduce drug tolerance, or modulate the host immune response to control the growth of intracellular M. tuberculosis. Recent developments in genomics and bioinformatics along with new drug discovery collectively have the potential to result in synergistic impact leading to the development of a rapid protocol to determine the drug resistance profile of the infecting strain so as to provide personalized medicine. Hence, in this review, we discuss recent developments in WGS, bioinformatics and drug discovery to perceive how they would transform the management of tuberculosis in a timely manner.

Rifampicin resistance mutations in the 81 bp RRDR of rpoB gene in Mycobacterium tuberculosis clinical isolates using Xpert MTB/RIF in Khyber Pakhtunkhwa, Pakistan: a retrospective study

BACKGROUND

Multi-drug resistant tuberculosis (MDR-TB) is a major public health problem especially in developing countries. World Health Organization (WHO) recommends use of Xpert MTB/RIF assay to simultaneously detecting Mycobacterium tuberculosis (MTB) and rifampicin (RIF) resistance. The primary objective of this study was to determine the frequency of MDR-TB in patients suspected to have drug resistance in Khyber Pakhtunkhwa. The frequency of probes for various rpoB gene mutations using Xpert MTB/RIF assay within 81 bp RRDR (Rifampicin Resistance Determining Region) was the secondary objective.

METHODS

A total of 2391 specimens, received at Programmatic Management of Drug Resistant TB (PMDT) Unit, Lady Reading Hospital (LRH) Peshawar, Pakistan, between October 2011 and December 2014, were analyzed by Xpert MTB/RIF test. MTB positive with rifampicin resistance were further analyzed to first line anti-mycobacterial drug susceptibility testing (DST) using middle brook 7H10 medium. The data was analyzed using statistical software; SPSS version 18.

RESULTS

Out of 2391 specimens, 1408 (59 %) were found positive for MTB and among them, 408 (29 %) showed rifampicin-resistance with four different rpoB gene mutations within 81 bp RRDR. The frequency of various probes among RIF-resistant isolates was observed as: probe E, in 314 out of 408 isolates; B, 44 out of 408; A, 5 out of 408; D, 34 out of 408; and probe C was observed among 6 out of 408 RIF-resistant isolates. The probe A&B and E&D mutation combination was found in only 1 isolate in each case, while B&D mutation combination was detected among 3 out of 408 RIF-resistant isolates.

CONCLUSIONS

Hence, it is concluded from our study on a selected population, 29 % of patients had MDR-TB. Probe E related mutations (also known as codon 531and 533) were the most common rpoB genetic mutation [314 (77 %)], acknowledged by Xpert MTB/RIF assay. Least mutation was detected within the sequence 511 (1.2 %).

Paradoxical Hypersusceptibility of Drug-resistant Mycobacteriumtuberculosis to β-lactam Antibiotics

Mycobacterium tuberculosis (M. tuberculosis) is considered innately resistant to β-lactam antibiotics. However, there is evidence that susceptibility to β-lactam antibiotics in combination with β–lactamase inhibitors is variable among clinical isolates, and these may present therapeutic options for drug-resistant cases. Here we report our investigation of susceptibility to β-lactam/β–lactamase inhibitor combinations among clinical isolates of M. tuberculosis, and the use of comparative genomics to understand the observed heterogeneity in susceptibility. Eighty-nine South African clinical isolates of varying first and second-line drug susceptibility patterns and two reference strains of M. tuberculosis underwent minimum inhibitory concentration (MIC) determination to two β-lactams: amoxicillin and meropenem, both alone and in combination with clavulanate, a β–lactamase inhibitor. 41/91 (45%) of tested isolates were found to be hypersusceptible to amoxicillin/clavulanate relative to reference strains, including 14/24 (58%) of multiple drug-resistant (MDR) and 22/38 (58%) of extensively drug-resistant (XDR) isolates. Genome-wide polymorphisms identified using whole-genome sequencing were used in a phylogenetically-aware linear mixed model to identify polymorphisms associated with amoxicillin/clavulanate susceptibility. Susceptibility to amoxicillin/clavulanate was over-represented among isolates within a specific clade (LAM4), in particular among XDR strains. Twelve sets of polymorphisms were identified as putative markers of amoxicillin/clavulanate susceptibility, five of which were confined solely to LAM4. Within the LAM4 clade, ‘paradoxical hypersusceptibility’ to amoxicillin/clavulanate has evolved in parallel to first and second-line drug resistance. Given the high prevalence of LAM4 among XDR TB in South Africa, our data support an expanded role for β-lactam/β-lactamase inhibitor combinations for treatment of drug-resistant M. tuberculosis.

•

Paradoxical hypersusceptibility is observed drug susceptibility despite innate resistance in the wild type state.

•

Many MDR and XDR M. tuberculosis strains are susceptible to amoxicillin/clavulanate.

•

Whole-genome sequencing identified mutations associated with paradoxical hypersusceptibility.

•

An expanded role for β-lactams in drug-resistant M. tuberculosis is supported.

The global increase in drug-resistant tuberculosis has prompted a search for alternative therapies, including repurposing existing antibiotics. β-lactam antibiotics are safe drugs, however, they have previously been thought to be of limited use for tuberculosis due to innate resistance to this drug class. In this study, the authors found many drug-resistant tuberculosis isolates from South Africa to be susceptible to a β-lactam and β-lactamase combination, amoxicillin/clavulanate. With the use of comparative genomics, multiple genetic mutations were identified to be associated with this hypersusceptible phenotype. These findings support an expanded role of β-lactam/β-lactamase inhibitor combinations for treatment of drug-resistant TB.

Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogen's phylogeography

Currently, Mycobacterium tuberculosis isolates of Latin-American Mediterranean (LAM) family may be detected far beyond the geographic areas that coined its name 15years ago. Here, we established the framework phylogeny of this geographically intriguing and pathobiologically important mycobacterial lineage and hypothesized how human demographics and migration influenced its phylogeography. Phylogenetic analysis of LAM isolates from all continents based on 24 variable number of tandem repeats (VNTR) loci and other markers identified three global sublineages with certain geographic affinities and defined by large deletions RD115, RD174, and by spoligotype SIT33. One minor sublineage (spoligotype SIT388) appears endemic in Japan. One-locus VNTR signatures were established for sublineages and served for their search in published literature and geographic mapping. We suggest that the LAM family originated in the Western Mediterranean region. The most widespread RD115 sublineage seems the most ancient and encompasses genetically and geographically distant branches, including extremely drug resistant KZN in South Africa and LAM-RUS recently widespread across Northern Eurasia. The RD174 sublineage likely started its active spread in Brazil; its earlier branch is relatively dominated by isolates from South America and the derived one is dominated by Portuguese and South/Southeastern African isolates. The relatively most recent SIT33-sublineage is marked with enigmatic gaps and peaks across the Americas and includes South African clade F11/RD761, which likely emerged within the SIT33 subpopulation after its arrival to Africa. In addition to SIT388-sublineage, other deeply rooted, endemic LAM sublineages may exist that remain to be discovered. As a general conclusion, human mass migration appears to be the major factor that shaped the M. tuberculosis phylogeography over large time-spans.

A Novel Molecular Strategy for Surveillance of Multidrug Resistant Tuberculosis in High Burden Settings

BACKGROUND

In South Africa and other high prevalence countries, transmission is a significant contributor to rising rates of multidrug resistant tuberculosis (MDR-TB). Thus, there is a need to develop an early detection system for transmission clusters suitable for high burden settings. We have evaluated the discriminatory power and clustering concordance of a novel and simple genotyping approach, combining spoligotyping with pncA sequencing (SpoNC), against two well-established methods: IS6110-RFLP and 24-loci MIRU-VNTR.

METHODS

A total of 216 MDR-TB isolates collected from January to June 2010 from the NHLS Central TB referral laboratory in Braamfontein, Johannesburg, representing a diversity of strains from South Africa, were included. The isolates were submitted for genotyping, pncA sequencing and analysis to the Centre for Tuberculosis in South Africa and the Public Health Research Institute Tuberculosis Center at Rutgers University in the United States. Clustering rates, Hunter-Gaston Discriminatory Indexes (HGI) and Wallace coefficients were compared between the methods.

RESULTS

Overall clustering rates were high by both IS6110-RFLP (52.8%) and MIRU-VNTR (45.8%), indicative of on-going transmission. Both 24-loci MIRU-VNTR and IS6110-RFLP had similar HGI (0.972 and 0.973, respectively), with close numbers of unique profiles (87 vs. 70), clustered isolates (129 vs. 146), and cluster sizes (2 to 26 vs. 2 to 25 isolates). Spoligotyping alone was the least discriminatory (80.1% clustering, HGI 0.903), with 28 unique types. However, the discriminatory power of spoligotyping was improved when combined with pncA sequencing using the SpoNC approach (61.8% clustering, HGI 0.958). A high proportion of MDR-TB isolates had mutations in pncA (68%, n = 145), and pncA mutations were significantly associated with clustering (p = 0.007 and p = 0.0013 by 24-loci MIRU-VNTR and IS6110-RFLP, respectively), suggesting high rates of resistance to pyrazinamide among all MDR-TB cases and particularly among clustered cases.

CONCLUSION

We conclude that SpoNC provides good discrimination for MDR-TB surveillance and early identification of outbreaks in South Africa, with 24-loci MIRU-VNTR applied for pncA wild-type strains as needed.

Mycobacterium tuberculosis Type II Toxin-Antitoxin Systems: Genetic Polymorphisms and Functional Properties and the Possibility of Their Use for Genotyping

Various genetic markers such as IS-elements, DR-elements, variable number tandem repeats (VNTR), single nucleotide polymorphisms (SNPs) in housekeeping genes and other groups of genes are being used for genotyping. We propose a different approach. We suggest the type II toxin-antitoxin (TA) systems, which play a significant role in the formation of pathogenicity, tolerance and persistence phenotypes, and thus in the survival of Mycobacterium tuberculosis in the host organism at various developmental stages (colonization, infection of macrophages, etc.), as the marker genes. Most genes of TA systems function together, forming a single network: an antitoxin from one pair may interact with toxins from other pairs and even from other families. In this work a bioinformatics analysis of genes of the type II TA systems from 173 sequenced genomes of M. tuberculosis was performed. A number of genes of type II TA systems were found to carry SNPs that correlate with specific genotypes. We propose a minimally sufficient set of genes of TA systems for separation of M. tuberculosis strains at nine basic genotype and for further division into subtypes. Using this set of genes, we genotyped a collection consisting of 62 clinical isolates of M. tuberculosis. The possibility of using our set of genes for genotyping using PCR is also demonstrated.

Evolution of Extensively Drug-Resistant Tuberculosis over Four Decades: Whole Genome Sequencing and Dating Analysis of Mycobacterium tuberculosis Isolates from KwaZulu-Natal

BACKGROUND

The continued advance of antibiotic resistance threatens the treatment and control of many infectious diseases. This is exemplified by the largest global outbreak of extensively drug-resistant (XDR) tuberculosis (TB) identified in Tugela Ferry, KwaZulu-Natal, South Africa, in 2005 that continues today. It is unclear whether the emergence of XDR-TB in KwaZulu-Natal was due to recent inadequacies in TB control in conjunction with HIV or other factors. Understanding the origins of drug resistance in this fatal outbreak of XDR will inform the control and prevention of drug-resistant TB in other settings. In this study, we used whole genome sequencing and dating analysis to determine if XDR-TB had emerged recently or had ancient antecedents.

METHODS AND FINDINGS

We performed whole genome sequencing and drug susceptibility testing on 337 clinical isolates of Mycobacterium tuberculosis collected in KwaZulu-Natal from 2008 to 2013, in addition to three historical isolates, collected from patients in the same province and including an isolate from the 2005 Tugela Ferry XDR outbreak, a multidrug-resistant (MDR) isolate from 1994, and a pansusceptible isolate from 1995. We utilized an array of whole genome comparative techniques to assess the relatedness among strains, to establish the order of acquisition of drug resistance mutations, including the timing of acquisitions leading to XDR-TB in the LAM4 spoligotype, and to calculate the number of independent evolutionary emergences of MDR and XDR. Our sequencing and analysis revealed a 50-member clone of XDR M. tuberculosis that was highly related to the Tugela Ferry XDR outbreak strain. We estimated that mutations conferring isoniazid and streptomycin resistance in this clone were acquired 50 y prior to the Tugela Ferry outbreak (katG S315T [isoniazid]; gidB 130 bp deletion [streptomycin]; 1957 [95% highest posterior density (HPD): 1937-1971]), with the subsequent emergence of MDR and XDR occurring 20 y (rpoB L452P [rifampicin]; pncA 1 bp insertion [pyrazinamide]; 1984 [95% HPD: 1974-1992]) and 10 y (rpoB D435G [rifampicin]; rrs 1400 [kanamycin]; gyrA A90V [ofloxacin]; 1995 [95% HPD: 1988-1999]) prior to the outbreak, respectively. We observed frequent de novo evolution of MDR and XDR, with 56 and nine independent evolutionary events, respectively. Isoniazid resistance evolved before rifampicin resistance 46 times, whereas rifampicin resistance evolved prior to isoniazid only twice. We identified additional putative compensatory mutations to rifampicin in this dataset. One major limitation of this study is that the conclusions with respect to ordering and timing of acquisition of mutations may not represent universal patterns of drug resistance emergence in other areas of the globe.

CONCLUSIONS

In the first whole genome-based analysis of the emergence of drug resistance among clinical isolates of M. tuberculosis, we show that the ancestral precursor of the LAM4 XDR outbreak strain in Tugela Ferry gained mutations to first-line drugs at the beginning of the antibiotic era. Subsequent accumulation of stepwise resistance mutations, occurring over decades and prior to the explosion of HIV in this region, yielded MDR and XDR, permitting the emergence of compensatory mutations. Our results suggest that drug-resistant strains circulating today reflect not only vulnerabilities of current TB control efforts but also those that date back 50 y. In drug-resistant TB, isoniazid resistance was overwhelmingly the initial resistance mutation to be acquired, which would not be detected by current rapid molecular diagnostics employed in South Africa that assess only rifampicin resistance.

Genome Analysis of the First Extensively Drug-Resistant (XDR) Mycobacterium tuberculosis in Malaysia Provides Insights into the Genetic Basis of Its Biology and Drug Resistance

The outbreak of extensively drug-resistant tuberculosis (XDR-TB) has become an increasing problem in many TB-burdened countries. The underlying drug resistance mechanisms, including the genetic variation favored by selective pressure in the resistant population, are partially understood. Recently, the first case of XDR-TB was reported in Malaysia. However, the detailed genotype family and mechanisms of the formation of multiple drugs resistance are unknown. We sequenced the whole genome of the UM 1072388579 strain with a 2-kb insert-size library and combined with that from previously sequenced 500-bp-insert paired-end reads to produce an improved sequence with maximal sequencing coverage across the genome. In silico spoligotyping and phylogenetic analyses demonstrated that UM 1072388579 strain belongs to an ancestral-like, non-Beijing clade of East Asia lineage. This is supported by the presence of a number of lineage-specific markers, including fadD28, embA, nuoD and pks7. Polymorphism analysis showed that the drug-susceptibility profile is correlated with the pattern of resistance mutations. Mutations in drug-efflux pumps and the cell wall biogenesis pathway such as mmpL, pks and fadD genes may play an important role in survival and adaptation of this strain to its surrounding environment. In this work, fifty-seven putative promoter SNPs were identified. Among them, we identified a novel SNP located at -4 T allele of TetR/acrR promoter as an informative marker to recognize strains of East Asian lineage. Our work indicates that the UM 1072388579 harbors both classical and uncommon SNPs that allow it to escape from inhibition by many antibiotics. This study provides a strong foundation to dissect the biology and underlying resistance mechanisms of the first reported XDR M. tuberculosis in Malaysia.