Genetic variation in Mycobacterium tuberculosis isolates from a London outbreak associated with isoniazid resistance

A comparison of phenotypic and WGS drug susceptibility testing in Mycobacterium tuberculosis isolates from the Republic of Korea

Background

WGS has significant potential to help tackle the major public health problem of TB. The Republic of Korea has the third highest rates of TB of all Organisation for Economic Cooperation and Development countries but there has been very limited use of WGS in TB to date.

Objectives

A retrospective comparison of Mycobacterium tuberculosis (MTB) clinical isolates from 2015 to 2017 from two centres in the Republic of Korea using WGS to compare phenotypic drug susceptibility testing (pDST) and WGS drug susceptibility predictions (WGS-DSP).

Methods

Fifty-seven MTB isolates had DNA extracted and were sequenced using the Illumina HiSeq platform. The WGS analysis was performed using bwa mem, bcftools and IQ-Tree; resistance markers were identified using TB profiler. Phenotypic susceptibilities were carried out at the Supranational TB reference laboratory (Korean Institute of Tuberculosis).

Results

For first-line antituberculous drugs concordance for rifampicin, isoniazid, pyrazinamide and ethambutol was 98.25%, 92.98%, 87.72% and 85.96%, respectively. The sensitivity of WGS-DSP compared with pDST for rifampicin, isoniazid, pyrazinamide and ethambutol was 97.30%, 92.11%, 78.95% and 95.65%, respectively. The specificity for these first-line antituberculous drugs was 100%, 94.74%, 92.11% and 79.41%, respectively. The sensitivity and specificity for second-line drugs ranged from 66.67% to 100%, and from 82.98% to 100%, respectively.

Conclusions

This study confirms the potential role for WGS in drug susceptibility prediction, which would reduce turnaround times. However, further larger studies are needed to ensure current databases of drug resistance mutations are reflective of the TB present in the Republic of Korea.

Origin and Global Expansion of Mycobacterium tuberculosis Complex Lineage 3

Simple Summary

Tuberculosis still causes 1.5 million deaths annually and is mainly caused by Mycobacterium tuberculosis complex strains belonging to three evolutionary modern lineages (Lineages 2–4). While Lineage 2 and Lineage 4 virtually conquered the world, Lineage 3 is particularly successful in Northern and Eastern Africa, as well as in Southern Asia, the suspected evolutionary origin of these strains. Here, we sought to understand how Lineage 3 strains came to the African continent. To this end, we performed routine genotyping to characterize over 2500 clinical isolates from 38 countries. We then selected a representative collection of 373 isolates for a whole-genome analysis and a modeling approach to infer the geographic origin of different sublineages. In fact, the origin of Lineage 3 could be located in India, and we found evidence for independent introductions of four distinct sublineages into North/East Africa, in line with known ancient exchanges and migrations between both world regions. Our study illustrates that the evolutionary history of humans and their pathogens are closely connected and further provides a systematic understanding of the genomic diversity of Lineage 3, which could be important for the development of new tuberculosis vaccines or new therapeutics.

Abstract

Mycobacterium tuberculosis complex (MTBC) Lineage 3 (L3) strains are abundant in world regions with the highest tuberculosis burden. To investigate the population structure and the global diversity of this major lineage, we analyzed a dataset comprising 2682 L3 strains from 38 countries over 5 continents, by employing 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats genotyping (MIRU-VNTR) and drug susceptibility testing. We further combined whole-genome sequencing (WGS) and phylogeographic analysis for 373 strains representing the global L3 genetic diversity. Ancestral state reconstruction confirmed that the origin of L3 strains is located in Southern Asia and further revealed multiple independent introduction events into North-East and East Africa. This study provides a systematic understanding of the global diversity of L3 strains and reports phylogenetic variations that could inform clinical trials which evaluate the effectivity of new drugs/regimens or vaccine candidates.

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.

Role of Whole-Genome Sequencing in Characterizing the Mechanism of Action of para-Aminosalicylic Acid and Its Resistance

para-Aminosalicylic acid (PAS) remains one of the drugs of last resort for the treatment of tuberculosis, but its mechanism of action is still not completely understood. The main aim of this project was to identify new potential mechanisms of action and resistance to PAS by performing whole-genome sequencing on PAS-resistant laboratory mutants. A new variant in the folC gene was identified, as well as some other mutations that require further study.

Cloning and Protein Expression of eccB5 Gene in ESX-5 System from Mycobacterium tuberculosis

Mycobacterium tuberculosis (M. tuberculosis) is the causative agent of tuberculosis in human. One of the major M. tuberculosis virulence factors is early secretory antigenic target of 6-kDa (ESAT-6), and EccB5 protein encoded by eccB5 is one of its components. EccB5 protein is a transmembrane protein in ESX-5 system. The aim of this study is to explore the characteristics of wild-type EccB5 and its mutant form N426I. We expressed the EccB5 protein by cloning the mutant and wild-type eccB5 gene in Escherichia coli (E. coli). We compared the protein structure of wild type and mutant form of EccB5 and found changes in structure around Asn426 (loop structure) in wild type and around Ile426 (β-strand) in the mutant. The truncated recombinant protein of EccB5 was successfully cloned and expressed using plasmid pCold I in E. coli DH5α and E. coli strain Rosetta-gami B (DE3) and purified as a 38.6 kDa protein by using the affinity column. There was no detectable adenosine triphosphatase activity in truncated forms of EccB5 and its mutant. In conclusion, our study reveals successful cloning and protein expression of truncated form of eccB5 gene of M. tuberculosis. EccB5 protein in ESX-5 system may be an important membrane component involved in the transport machinery of type VII secretion system, which is essential for growth and virulence.

Molecular and genomic typing for tuberculosis surveillance: A survey study in 26 European countries

BACKGROUND

Molecular typing and whole genome sequencing (WGS) information is used for (inter-) national outbreak investigations. To assist the implementation of these techniques for tuberculosis (TB) surveillance and outbreak investigations at European level there is a need for inter-country collaboration and standardization. This demands more information on molecular typing practices and capabilities of individual countries. We aimed to review the use of molecular/genomic typing for TB surveillance in European Union and European Economic Area countries in 2016; assess its public health value; and collect experiences on typing data use for cross-border cluster investigations.

METHOD

A web-based questionnaire was provided to all TB National Focal Points. The questionnaire consisted of three parts: i) Use and integration of molecular and genomic typing data into TB surveillance; ii) Cross-border cluster investigation and international collaboration, and iii) Perception and evaluation of public health benefits of molecular and genomic typing for TB surveillance.

RESULTS

Of 26 responding countries, 20 used molecular typing for TB surveillance, including nine applying WGS. The level of integration into the national surveillance was heterogeneous. Among six countries not using typing for TB surveillance, more than half planned its implementation soon. Overall, most countries perceived an added public health value of molecular typing for TB control. Concerning international cluster investigations, countries had little experience and did not have standard protocols to exchange typing data.

CONCLUSION

Our study shows a wide use of molecular and genomic typing data for TB surveillance in EU/EEA countries and reveals that transition to WGS-based typing is ongoing or is considered in most countries. However, our results also show a high heterogeneity in the use and integration of typing data for TB surveillance. Standardization of typing data use for TB surveillance is needed and formal procedures should be developed to facilitate international collaboration.

Drug resistance mechanisms and drug susceptibility testing for tuberculosis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (MTB) is the deadliest infectious disease and the associated global threat has worsened with the emergence of drug resistance, in particular multidrug-resistant TB (MDR-TB) and extensively drug-resistant TB (XDR-TB). Although the World Health Organization (WHO) End-TB Strategy advocates for universal access to antimicrobial susceptibility testing, this is not widely available and/or it is still underused. The majority of drug resistance in clinical MTB strains is attributed to chromosomal mutations. Resistance-related mutations could also exert certain fitness cost to the drug-resistant MTB strains and growth fitness could be restored by the presence of compensatory mutations. Understanding these underlying mechanisms could provide an important insight into TB pathogenesis and predict the future trend of MDR-TB global pandemic. This review covers the mechanisms of resistance in MTB and provides a comprehensive overview of current phenotypic and molecular approaches for drug susceptibility testing, with particular attention to the methods endorsed and recommended by the WHO.

Alternate splicing of transcripts upon Mycobacterium tuberculosis infection impacts the expression of functional protein domains

Previously, we reported that infection of human macrophages with Mycobacterium tuberculosis (Mtb) results in massive alterations in the pattern of RNA splicing in the host. The finding gained significance since alternate spliced variants of a same gene may have substantially different structure, function, stability, interaction partners, localization, and so forth, owing to inclusion or exclusion of specific exons. To establish a proof-of-concept; on how infection-induced RNA splicing could impact protein functions, here we used RNA-seq data from THP-1 macrophages that were infected with clinical isolate of Mtb. In addition to re-establishing the fact that Mtb infection may cause strain specific alterations in RNA splicing, we also developed a new analysis pipeline resulting in characterization of domain maps of the transcriptome postinfection. For the sake of simplicity, we restricted our analysis to all the kinases in the human genome and considered only pfam classified protein domains and checked their frequency of inclusion or exclusion due to alternate splicing across the conditions and time points. We report massive alterations in the domain architecture of most regulated proteins across the entire kinases highlighting the physiological importance of such an understanding. This study paves way for more detailed analysis of different functional classes of proteins and perturbations to their domain architecture as a consequence of mycobacterial infections. Such analysis would yield unprecedented depth to our understanding of host-pathogen interaction and allow in a more systematic manner targeting of host pathways for controlling the infections.

SINGLE NUCLEOTIDE POLYMORPHISM OF ECCB5 GENE OF MYCOBACTERIUM TUBERCULOSIS COMPLEX ISOLATES FROM SUSPECTED PULMONARY TB PATIENTS IN SURABAYA INDONESIA

Background:

Mycobacterium tuberculosis Complex (MTBC) is a group of Mycobacterium that causes tuberculosis (TB). TB is an infectious disease that remains a global health problem. Indonesia is one of the five countries in the world where TB is the most prevalent and became the country with tle second largest rate of TB in 2014 and 2015. MTBC has high pathogenicity that can cause infections in animals and humans. The most common route of transmission is via airborne droplet nuclei and contact with animals or humans infected with TB. MTBC has many virulence factors. One of these factors is EccB₅ that is encoded by eccB₅ gene. EccB5 is a transmembrane protein-conserved membrane protein and could play a role in inducing damage in host cells, macrophage infection, and may correlate with active disease. The characterization of eccB₅ gene needs to be studied to determine the nucleotide sequences, which may be associated with active disease. The aim of this research was to analyze the nuclotide sequences of eccB₅ gene of MTBC from suspected pulmonary tuberculosis patients, SNPs of eccB₅ gene and possible correlation with the disease, especially in Indonesia.

Materials and Methods:

Samples were collected from the Tuberculosis Laboratory, Clinical Microbiology of Dr. Soetomo Hospital Surabaya Indonesia. DNA extraction used boiling extraction method and continued nucleic acid amplification using PCR techniques. Primer pairs used eccB₅ SK.. The positivity of DNA specific revealed amplicon in 1592 bp. PCR product was sequenced by 1st Base (First BASE Laboratories Sdn Bhd, Selangor, Malaysia). The sequence analysis used Genetyx-Win version 10.0 (Genetyx Corporation, Tokyo, Japan).

Results:

Total isolates of Mycobacterium spp. were 28 and those that showed positive MTBC were 24 isolates and 4 nontuberculosis mycobacteria (NTM) using immunochromatographic test (ICT). The amount of homology from MTBC using blast NCBI was 99%-100%. Two SNPs were found in position c.1277 which revealed replacement of amino acid in 426 of codon position.

Conclusion:

The sequence of eccB₅ gene of MTBC showed high significant homology, while proposed non-synoymous single nucleotide polymorphisms (nsSNP) may associated with clinical outcomes.

Comparative Genomic Analysis of Two Clonally Related Multidrug Resistant Mycobacterium tuberculosis by Single Molecule Real Time Sequencing

Background: Multidrug-resistant tuberculosis (MDR-TB) is posing a major threat to global TB control. In this study, we focused on two consecutive MDR-TB isolated from the same patient before and after the initiation of anti-TB treatment. To better understand the genomic characteristics of MDR-TB, Single Molecule Real-Time (SMRT) Sequencing and comparative genomic analyses was performed to identify mutations that contributed to the stepwise development of drug resistance and growth fitness in MDR-TB under in vivo challenge of anti-TB drugs. Result: Both pre-treatment and post-treatment strain demonstrated concordant phenotypic and genotypic susceptibility profiles toward rifampicin, pyrazinamide, streptomycin, fluoroquinolones, aminoglycosides, cycloserine, ethionamide, and para-aminosalicylic acid. However, although both strains carried identical missense mutations at rpoB S531L, inhA C-15T, and embB M306V, MYCOTB Sensititre assay showed that the post-treatment strain had 16-, 8-, and 4-fold elevation in the minimum inhibitory concentrations (MICs) toward rifabutin, isoniazid, and ethambutol respectively. The results have indicated the presence of additional resistant-related mutations governing the stepwise development of MDR-TB. Further comparative genomic analyses have identified three additional polymorphisms between the clinical isolates. These include a single nucleotide deletion at nucleotide position 360 of rv0888 in pre-treatment strain, and a missense mutation at rv3303c (lpdA) V44I and a 6-bp inframe deletion at codon 67-68 in rv2071c (cobM) in the post-treatment strain. Multiple sequence alignment showed that these mutations were occurring at highly conserved regions among pathogenic mycobacteria. Using structural-based and sequence-based algorithms, we further predicted that the mutations potentially have deleterious effect on protein function. Conclusion: This is the first study that compared the full genomes of two clonally-related MDR-TB clinical isolates during the course of anti-TB treatment. Our work has demonstrated the robustness of SMRT Sequencing in identifying mutations among MDR-TB clinical isolates. Comparative genome analysis also suggested novel mutations at rv0888, lpdA, and cobM that might explain the difference in antibiotic resistance and growth pattern between the two MDR-TB strains.

Mycobacterium tuberculosis and whole-genome sequencing: how close are we to unleashing its full potential?

BACKGROUND

Nearly two decades after completion of the genome sequence of Mycobacterium tuberculosis (MTB), and with the advent of next generation sequencing technologies (NGS), whole-genome sequencing (WGS) has been applied to a wide range of clinical scenarios. Starting in 2017, England is the first country in the world to pioneer its use on a national scale for the diagnosis of tuberculosis, detection of drug resistance, and typing of MTB.

AIMS

This narrative review critically analyses the current applications of WGS for MTB and explains how close we are to realizing its full potential as a diagnostic, epidemiologic, and research tool.

SOURCES

We searched for reports (both original articles and reviews) published in English up to 31 May 2017, with combinations of the following keywords: whole-genome sequencing, Mycobacterium, and tuberculosis. MEDLINE, Embase, and Scopus were used as search engines. We included articles that covered different aspects of whole-genome sequencing in relation to MTB.

CONTENT

This review focuses on three main themes: the role of WGS for the prediction of drug susceptibility, MTB outbreak investigation and genetic diversity, and research applications of NGS.

IMPLICATIONS

Many of the original expectations have been accomplished, and we believe that with its unprecedented sensitivity and power, WGS has the potential to address many unanswered questions in the near future. However, caution is still needed when interpreting WGS data as there are some important limitations to be aware of, from correct interpretation of drug susceptibilities to the bioinformatic support needed.

Twenty years and counting: epidemiology of an outbreak of isoniazid-resistant tuberculosis in England and Wales, 1995 to 2014

An outbreak of isoniazid-resistant tuberculosis first identified in London has now been ongoing for 20 years, making it the largest drug-resistant outbreak of tuberculosis documented to date worldwide. We identified culture-confirmed cases with indistinguishable molecular strain types and extracted demographic, clinical, microbiological and social risk factor data from surveillance systems. We summarised changes over time and used kernel-density estimation and k-function analysis to assess geographic clustering. From 1995 to 2014, 508 cases were reported, with a declining trend in recent years. Overall, 70% were male (n = 360), 60% born in the United Kingdom (n = 306), 39% white (n = 199), and 26% black Caribbean (n = 134). Median age increased from 25 years in the first 5 years to 42 in the last 5. Approximately two thirds of cases reported social risk factors: 45% drug use (n = 227), 37% prison link (n = 189), 25% homelessness (n = 125) and 13% alcohol dependence (n = 64). Treatment was completed at 12 months by 52% of cases (n = 206), and was significantly lower for those with social risk factors (p < 0.05), but increased over time for all patients (p < 0.05). The outbreak remained focused in north London throughout. Control of this outbreak requires continued efforts to prevent and treat further active cases through targeted screening and enhanced case management.

Use of whole-genome sequencing to distinguish relapse from reinfection in a completed tuberculosis clinical trial

BACKGROUND

RIFAQUIN was a tuberculosis chemotherapy trial in southern Africa including regimens with high-dose rifapentine with moxifloxacin. Here, the application of whole-genome sequencing (WGS) is evaluated within RIFAQUIN for identifying new infections in treated patients as either relapses or reinfections. WGS is further compared with mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR) typing. This is the first report of WGS being used to evaluate new infections in a completed clinical trial for which all treatment and epidemiological data are available for analysis.

METHODS

DNA from 36 paired samples of Mycobacterium tuberculosis cultured from patients before and after treatment was typed using 24-loci MIRU-VNTR, in silico spoligotyping and WGS. Following WGS, the sequences were mapped against the reference strain H37Rv, the single-nucleotide polymorphism (SNP) differences between pairs were identified, and a phylogenetic reconstruction was performed.

RESULTS

WGS indicated that 32 of the paired samples had a very low number of SNP differences (0-5; likely relapses). One pair had an intermediate number of SNP differences, and was likely the result of a mixed infection with a pre-treatment minor genotype that was highly related to the post-treatment genotype; this was reclassified as a relapse, in contrast to the MIRU-VNTR result. The remaining three pairs had very high SNP differences (>750; likely reinfections).

CONCLUSIONS

WGS and MIRU-VNTR both similarly differentiated relapses and reinfections, but WGS provided significant extra information. The low proportion of reinfections seen suggests that in standard chemotherapy trials with up to 24 months of follow-up, typing the strains brings little benefit to an analysis of the trial outcome in terms of differentiating relapse and reinfection. However, there is a benefit to using WGS as compared to MIRU-VNTR in terms of the additional genotype information obtained, in particular for defining the presence of mixed infections and the potential to identify known and novel drug-resistance markers.