The implications of whole-genome sequencing in the control of tuberculosis

New insight in molecular detection of Mycobacterium tuberculosis

Mycobacterium tuberculosis, the causative agent of tuberculosis, is a pathogenic bacterium that has claimed millions of lives since the Middle Ages. According to the World Health Organization's report, tuberculosis ranks among the ten deadliest diseases worldwide. The presence of an extensive array of genes and diverse proteins within the cellular structure of this bacterium has provided us with a potent tool for diagnosis. While the culture method remains the gold standard for tuberculosis diagnosis, it is possible that molecular diagnostic methods, emphasis on the identification of mutation genes (e.g., rpoB and gyrA) and single nucleotide polymorphisms, could offer a safe and reliable alternative. Over the past few decades, as our understanding of molecular genetics has expanded, methods have been developed based on gene expansion and detection. These methods typically commence with DNA amplification through nucleic acid targeted techniques such as polymerase chain reaction. Various molecular compounds and diverse approaches have been employed in molecular assays. In this review, we endeavor to provide an overview of molecular assays for the diagnosis of tuberculosis with their properties (utilization, challenges, and functions). The ultimate goal is to explore the potential of replacing traditional bacterial methods with these advanced molecular diagnostic techniques.

Whole genome sequencing of drug resistance Mycobacterium tuberculosis from extra-pulmonary sites

This study aimed to determinate characteristics of drug resistance Mycobacterium tuberculosis from patients with extra-pulmonary tuberculosis (EPTB). Patients were retrospectively studied from January 2020 to December 2021. All the isolates were cultured, tested drug susceptibility, and detected the gene mutation using whole genome sequencing. The correlations of whole genome sequencing, pattern of DR, patients' distribution, and transmission were analyzed. 111 DR-EPTB isolates included pre-XDR-TB (53.2%), MDR-TB (29.7%), and poly-DR-TB (12.6%). The resistant drugs were INH followed by RFP and SM. The genotypes of 111 strains were lineage 2 and lineage 4. KatG_p.Ser315Thr was main gene mutation for resistance to INH; rpsL_p.Lys43Arg for SM, rpoB_p.Ser450Leu for rifampicin, embB_p.Met306Val for ethambutol, gyrA_p.Asp94Gly for FQs, and pncA_p.Thr76Pro for PZA. The residence was a significant risk factor for cluster transmission by patients and phenotypic DR types of strains for lineage 2 transmission. In the local area of southwest China INH, rifampicin and SM were main drugs in patients with DR-EPTB. KatG_p.Ser315, rpoB_p.Ser450Leu, and rpsL_p.Lys43Arg were main gene mutations. Phenotypic DR types and residence were main risk of transmission.

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.

Identification of a predominant genotype of Mycobacterium tuberculosis in Brazilian indigenous population

After nearly a century of vaccination and six decades of drug therapy, tuberculosis (TB) kills more people annually than any other infectious disease. Substantial challenges to disease eradication remain among vulnerable and underserved populations. The Guarani-Kaiowá people are an indigenous population in Paraguay and the Brazilian state of Mato Grosso do Sul. This community, marginalized in Brazilian society, experiences severe poverty. Like other South American indigenous populations, their TB prevalence is high, but the disease has remained largely unstudied in their communities. Herein, Mycobacterium tuberculosis isolates from local clinics were whole genome sequenced, and a population genetic framework was generated. Phylogenetics show M. tuberculosis isolates in the Guarani-Kaiowá people cluster away from selected reference strains, suggesting divergence. Most cluster in a single group, further characterized as M. tuberculosis sublineage 4.3.3. Closer analysis of SNPs showed numerous variants across the genome, including in drug resistance-associated genes, and with many unique changes fixed in each group. We report that local M. tuberculosis strains have acquired unique polymorphisms in the Guarani-Kaiowá people, and drug resistance characterization is urgently needed to inform public health to ensure proper care and avoid further evolution and spread of drug-resistant TB.

An improved simple method for the identification of Mycobacteria by MALDI-TOF MS (Matrix-Assisted Laser Desorption- Ionization mass spectrometry)

The aim of this study was to establish a simple method for the rapid identification of Mycobacteria species by MALDI-TOF (Matrix-Assisted Laser Desorption/Ionization-Time of Flight Mass spectrometry) using the Bruker MALDI-TOF Biotyper system (Bruker Daltonik, Bremen, Germany). A multicentre, prospective, and single blind study was performed in three European Hospitals, two Spanish and one UK hospital from May to August 2018. The BD BACTEC MGIT (Becton Dickinson, Berks, UK) liquid culture system was used in all three centres for the growth of Mycobacteria. When signal positive, tubes were removed from the analyser and in addition to standard laboratory procedures were subcultured on blood agar plates for MALDI-TOF analysis. Plates were incubated aerobically for 1 to 7 days at 37 °C and inspected every day. Once any growth was visible, it was transferred to the steel target plate, overlaid with 1 μl of neat formic acid and 1 μl HCCA matrix (alpha hydroxyl 4 cinnamic acid), and analysed in a Bruker Biotyper MALDI-TOF. Results given by MALDI-TOF were compared with the reference methods used for identification in the different centres. At two Spanish hospitals, identification by MALDI-TOF was only attempted on presumptive non-tuberculosis mycobacteria (NTM) and the results were initially compared with the results obtained by a commercial reverse hybridisation assay, GenoType CM/AS (Hain Lifescience, Tübingen, Germany). At the UK Hospital, identification of any presumptive mycobacteria was attempted and compared with the results obtained by whole genome sequencing (WGS). Overall in 142/167 (85%) of cases the identifications obtained were concordant; all Mycobacterium tuberculosis (MTB) isolates 43/43 (100%), 57/76 (75%) of the rapid growing nontuberculous mycobacteria (NTM), and 42/48 (85%) slow growing NTM tested were identified correctly. We report a new, easy, cheap and quick method for isolation and identification of Mycobacterium spp. without the need for additional steps or equipment and this method is in routine used in all three centres.

Whole genome sequencing for drug-resistant tuberculosis management in South Africa: What gaps would this address and what are the challenges to implementation?

Global control of tuberculosis (TB) has been seriously impacted by the emergence and transmission of its drug-resistant forms. Delayed detection and incomplete characterisation of drug-resistant tuberculosis (DR-TB) contributes to morbidity and mortality, and to ongoing transmission of drug-resistant strains. Current culture-based and molecular diagnostic tools for TB present numerous disadvantages that could potentially lead to misdiagnosis, inappropriate treatment initiation and the amplification of drug resistance. The detection of drug-resistant tuberculosis (DR-TB) in South Africa relies on molecular diagnostic assays such as the Xpert MTB/RIF and line probe assays (MTBDRplus and MTBDRsl). However, these molecular assays are limited to detecting resistance to only a few first-line and second-line drugs. It is for this reason that next-generation sequencing (NGS) and bioinformatics pipelines have been developed for rapid detection of M. tuberculosis drug resistance, with the added advantage that sequence data could also have public health applications through understanding transmission patterns. This review highlights some of the challenges that are currently hampering the diagnosis and control of DR-TB in a high burden setting of the KwaZulu-Natal (KZN) province in South Africa. Shortfalls of current diagnostic techniques for DR-TB are discussed in detail and we also propose how these might be overcome with an accurate and rapid NGS system.

Deciphering drug resistance in Mycobacterium tuberculosis using whole-genome sequencing: progress, promise, and challenges

Tuberculosis (TB) is a global infectious threat that is intensified by an increasing incidence of highly drug-resistant disease. Whole-genome sequencing (WGS) studies of Mycobacterium tuberculosis, the causative agent of TB, have greatly increased our understanding of this pathogen. Since the first M. tuberculosis genome was published in 1998, WGS has provided a more complete account of the genomic features that cause resistance in populations of M. tuberculosis, has helped to fill gaps in our knowledge of how both classical and new antitubercular drugs work, and has identified specific mutations that allow M. tuberculosis to escape the effects of these drugs. WGS studies have also revealed how resistance evolves both within an individual patient and within patient populations, including the important roles of de novo acquisition of resistance and clonal spread. These findings have informed decisions about which drug-resistance mutations should be included on extended diagnostic panels. From its origins as a basic science technique, WGS of M. tuberculosis is becoming part of the modern clinical microbiology laboratory, promising rapid and improved detection of drug resistance, and detailed and real-time epidemiology of TB outbreaks. We review the successes and highlight the challenges that remain in applying WGS to improve the control of drug-resistant TB through monitoring its evolution and spread, and to inform more rapid and effective diagnostic and therapeutic strategies.

Whole genome sequencing for drug resistance determination in Mycobacterium tuberculosis

South Africa remains challenged with a high tuberculosis burden accompanied by an increase in drug resistant cases. We assessed the use of the Illumina MiSeq, a next-generation sequencing platform for whole genome sequencing, followed by bioinformatic analysis using a commercial software package to determine resistance to selected drugs used for Mycobacterium tuberculosis treatment in our setting. Whole genome sequencing shows potential as a diagnostic platform for the detection of drug resistance in Mycobacterium tuberculosis with the provision of information for several drugs simultaneously.

Molecular epidemiology of cattle tuberculosis in Mexico through whole-genome sequencing and spoligotyping

Mycobacterium bovis infection in cattle persists in Mexico, posing a threat to human health. Control of bovine tuberculosis, through the National Program Against Bovine Tuberculosis, has led to the decrease of disease prevalence in most of the country, except for high dairy production regions. Genotyping of M. bovis has been performed mainly by spoligotyping and variable number tandem repeats (VNTR), but higher resolution power can be useful for a finer definition of the spread of the disease. Whole genome sequencing and spoligotyping was performed for a set of 322 M. bovis isolates from different sources in Mexico: Baja California, Coahuila, Estado de Mexico, Guanajuato, Hidalgo, Jalisco, Queretaro and Veracruz, from dairy and beef cattle, as well as humans. Twelve main genetic clades were obtained through WGS and genetic diversity analysis. A clear differentiation of the Baja California isolates was seen as they clustered together exclusively. However, isolates from the central states showed no specific clustering whatsoever. Although WGS proves to have higher resolving power than spoligotyping, and since there was concordance between WGS and spoligotyping results, we consider that the latter is still an efficient and practical method for monitoring bovine tuberculosis in developing countries, where resources for higher technology are scarce.

A Quantitative Evaluation of MIRU-VNTR Typing Against Whole-Genome Sequencing for Identifying Mycobacterium tuberculosis Transmission: A Prospective Observational Cohort Study

Background

Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) typing is widely used in high-income countries to determine Mycobacterium tuberculosis relatedness. Whole-genome sequencing (WGS) is known to deliver greater specificity, but no quantitative prospective comparison has yet been undertaken.

Methods

We studied isolates from the English Midlands, sampled consecutively between 1 January 2012 and 31 December 2015. In addition to routinely performed MIRU-VNTR typing, DNA was extracted from liquid cultures and sequenced using Illumina technology. Demographic and epidemiological data for the relevant patients were extracted from the Enhanced Tuberculosis Surveillance system run by Public Health England. Closely related samples, defined using a threshold of five single nucleotide variants (SNVs), were compared to samples with identical MIRU-VNTR profiles, to samples from individuals with shared epidemiological risk factors, and to those with both characteristics.

Findings

1999 patients were identified for whom at least one M. tuberculosis isolate had been MIRU-VNTR typed and sequenced. Comparing epidemiological risk factors with close genetic relatedness, only co-residence had a positive predictive value of over 5%. Excluding co-resident individuals, 18.6% of patients with identical MIRU-VNTR profiles were within 5 SNVs. Where patients also shared social risk factors and ethnic group, this rose to 48%. Only 8% of MIRU-VNTR linked pairs in lineage 1 were within 5 SNV, compared to 31% in lineage 4.

Interpretation

In the setting studied, this molecular epidemiological study shows MIRU-VNTR typing and epidemiological risk factors are poorly predictive of close genomic relatedness, assessed by SNV. MIRU-VNTR performance varies markedly by lineage.

Funding

Public Health England, Health Innovation Challenge Fund, NIHR Health Protection Research Unit Oxford, NIHR Oxford Biomedical Research Centre.

•

M. tuberculosis relatedness was measured by MIRU-VNTR typing and single nucleotide variants (SNV).

•

In consecutive isolates in England, only 20% of samples with identical MIRU-VNTR profiles were within 5 SNVs.

•

MIRU-VNTR is more predictive of close relatedness in lineage 4 (Euro-American) lineages than in other lineages.

Genotyping of Mycobacterium tuberculosis using whole genome sequencing

Tuberculosis (TB) is considered one of the most serious infectious diseases worldwide. Effective control of tuberculosis infection involves multiple steps, such as reliable detection, treatment, an epidemiological control as a part of case management, and further surveillance and monitoring of TB spread in the human population. Due to the accelerating advances in molecular biology, especially in DNA sequencing, in the past decade, the application of these methods has become crucial for TB evolution studies, differentiation of Mycobacterium tuberculosis genotypes, and their distribution. Currently, several molecular genetic methods are available. The oldest typing methods (e.g., IS6110-RFLP, spoligotyping, and MIRU-VNTR) can discover the chain of transmission to the patient. Currently, whole genome sequencing facilitates is furthermore able to identify the source of infection, the transmission trays among individuals sharing the same isolate, as well as determination of the TB evolution and its resistance to antituberculotic agents. It is obvious that this technique will become a new gold standard in genotyping methods in tuberculosis molecular epidemiological studies. In this article, molecular genetic typing methods with a special focus on whole genome sequencing and data management are reviewed.

Genotyping Multidrug-Resistant Mycobacterium tuberculosis from Primary Sputum and Decontaminated Sediment with an Integrated Microfluidic Amplification Microarray Test

There is a growing awareness that molecular diagnostics for detect-to-treat applications will soon need a highly multiplexed mutation detection and identification capability. In this study, we converted an open-amplicon microarray hybridization test for multidrug-resistant (MDR) Mycobacterium tuberculosis into an entirely closed-amplicon consumable (an amplification microarray) and evaluated its performance with matched sputum and sediment extracts. Reproducible genotyping (the limit of detection) was achieved with ∼25 M. tuberculosis genomes (100 fg of M. tuberculosis DNA) per reaction; the estimated shelf life of the test was at least 18 months when it was stored at 4°C. The test detected M. tuberculosis in 99.1% of sputum extracts and 100% of sediment extracts and showed 100% concordance with the results of real-time PCR. The levels of concordance between M. tuberculosis and resistance-associated gene detection were 99.1% and 98.4% for sputum and sediment extracts, respectively. Genotyping results were 100% concordant between sputum and sediment extracts. Relative to the results of culture-based drug susceptibility testing, the test was 97.1% specific and 75.0% sensitive for the detection of rifampin resistance in both sputum and sediment extracts. The specificity for the detection of isoniazid (INH) resistance was 98.4% and 96.8% for sputum and sediment extracts, respectively, and the sensitivity for the detection of INH resistance was 63.6%. The amplification microarray reported the correct genotype for all discordant phenotype/genotype results. On the basis of these data, primary sputum may be considered a preferred specimen for the test. The amplification microarray design, shelf life, and analytical performance metrics are well aligned with consensus product profiles for next-generation drug-resistant M. tuberculosis diagnostics and represent a significant ease-of-use advantage over other hybridization-based tests for diagnosing MDR tuberculosis.

Same-Day Diagnostic and Surveillance Data for Tuberculosis via Whole-Genome Sequencing of Direct Respiratory Samples

Routine full characterization of Mycobacterium tuberculosis is culture based, taking many weeks. Whole-genome sequencing (WGS) can generate antibiotic susceptibility profiles to inform treatment, augmented with strain information for global surveillance; such data could be transformative if provided at or near the point of care. We demonstrate a low-cost method of DNA extraction directly from patient samples for M. tuberculosis WGS. We initially evaluated the method by using the Illumina MiSeq sequencer (40 smear-positive respiratory samples obtained after routine clinical testing and 27 matched liquid cultures). M. tuberculosis was identified in all 39 samples from which DNA was successfully extracted. Sufficient data for antibiotic susceptibility prediction were obtained from 24 (62%) samples; all results were concordant with reference laboratory phenotypes. Phylogenetic placement was concordant between direct and cultured samples. With Illumina MiSeq/MiniSeq, the workflow from patient sample to results can be completed in 44/16 h at a reagent cost of £96/£198 per sample. We then employed a nonspecific PCR-based library preparation method for sequencing on an Oxford Nanopore Technologies MinION sequencer. We applied this to cultured Mycobacterium bovis strain BCG DNA and to combined culture-negative sputum DNA and BCG DNA. For flow cell version R9.4, the estimated turnaround time from patient to identification of BCG, detection of pyrazinamide resistance, and phylogenetic placement was 7.5 h, with full susceptibility results 5 h later. Antibiotic susceptibility predictions were fully concordant. A critical advantage of MinION is the ability to continue sequencing until sufficient coverage is obtained, providing a potential solution to the problem of variable amounts of M. tuberculosis DNA in direct samples.

Accurate and effective multidrug-resistant Mycobacterium tuberculosis detection method using gap-filling ligation coupled with high-resolution capillary electrophoresis-based single strand conformation polymorphism

Tuberculosis (TB) has severely threatened public health via emerging multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mycobacterium tuberculosis (MTB) strains. For effective TB treatment, rapid, accurate, and multiplex detection of drug resistance is extremely important. However, conventional methods for TB diagnosis are time consuming and have a limited effect on treatment. Nucleic acid-based molecular detection methods have been developed as an effective MDR/XDR-TB diagnosis technology. Among the nucleic acid-based methods, ligation-dependent methods are attractive as MDR/XDR-MTB detection technologies, but multiplex analysis is limited by the detection method. Although an electrophoresis-based method is considered for multiple target detection because it is free from the errors pertaining to hybridization-based systems, the procedure of multiplex analysis is quite complicated owing to the DNA size-based separation system. In this study, we report an accurate, rapid, and simple multiple MDR/XDR-MTB detection technology using gap-filling ligation reaction coupled with high-resolution capillary electrophoresis-based single-strand conformation polymorphism. Using this system, rapid and accurate MDR/XDR-MTB detection is feasible via similar length probes without the complicated step of probe design. We found that this method could accurately and effectively detect highly polymorphic regions in specific codons associated with drug resistance.

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.

Real-Time Sequencing of Mycobacterium tuberculosis: Are We There Yet?

Whole-genome sequencing has taken a leading role in epidemiologic studies of tuberculosis, but thus far, its real-time clinical utility has been low, in part because of the requirement for culture. In their report in this issue, Votintseva et al. (A. A. Votintseva, P. Bradley, L. Pankhurst, C. del Ojo Elias, M. Loose, K. Nilgiriwala, A. Chatterjee, E. G. Smith, N. Sanderson, T. M. Walker, M. R. Morgan, D. H. Wyllie, A. S. Walker, T. E. A. Peto, D. W. Crook, and Z. Iqbal, J Clin Microbiol 55:1285-1298, 2017, https://doi.org/10.1128/JCM.02483-16) present a new method for extracting Mycobacterium tuberculosis DNA directly from smear-positive respiratory samples, making it feasible to generate drug resistance predictions and phylogenetic trees in 44 h with the Illumina MiSeq. They also illustrate the potential for a <24-h turnaround time from DNA extraction to clinically relevant results with Illumina MiniSeq and Oxford Nanopore Technologies MinION. We comment on the promise and limitations of these approaches.