Whole genome sequence analysis of multidrug-resistant Mycobacterium tuberculosis Beijing isolates from an outbreak in Thailand

A large geno-spatial cluster of multi-drug resistant tuberculosis outbreak in a western district of Thailand

The growing issue of drug resistance, particularly multidrug-resistant TB (MDR-TB), has exacerbated this problem. The rise of drug resistance TB is a severe global health concern. In Thailand, a persistent community outbreak of primary MDR-TB has been confirmed in the Tha Maka district of Kanchanaburi province, with an increasing prevalence of MDR-TB among newly diagnosed pulmonary tuberculosis cases. It was the first site in Thailand where a cluster of MDR-TB, caused by the Asian African 3 Modern Beijing strain, and XDR-TB, caused by L2.1, outbreaks were reported. This study aims to assess the MDR-TB outbreak in detail by characterizing the genomic profiles of the prevalent MDR-TB strains and examining their geographical distribution within the affected district. Through whole-genome sequencing (WGS) and bioinformatic analysis of 188 MTB isolates, the study identified three major phylogenetic lineages: the East Asian lineage (L2, 92 %), the Indo-Oceanic lineage (L1, 5.9 %), and the Euro-American lineage (L4, 2.1 %). The detailed sub-lineage distribution offers valuable insights into the predominant genetic clusters of M. tuberculosis within the sampled population. Notably, Lineage 2, specifically the L2.2.M3 sub-lineage, stood out as the dominant strain of MDR-TB, accounting for 77.7 % of the isolates. This finding underscores the significant prevalence of the L2.2.M3 sub-lineage and its potential role in the local transmission dynamics of tuberculosis. The high proportion and genetic homogeneity of the L2.2.M3 cluster among MDR-TB patients may indicate the strain's adaptation for more effective transmission within the Thai population. The increasing prevalence of this pathogenic strain could significantly impact tuberculosis control programs. Early diagnosis and contact tracing with chemotherapeutic preventive therapy for MDR-TB will be essential in inhibiting the spread and reactivation of these strains. Additionally, further studies are needed to prospectively identify transmission routes through contact tracing and real-time genotypic methods. It will also be crucial to ensure that future vaccines and/or recommended chemoprophylaxis therapy for MDR-TB will provide protection against these emerging strains.

Diversification of gene content in the Mycobacterium tuberculosis complex is determined by phylogenetic and ecological signatures

We analyzed the pan-genome and gene content modulation of the most diverse genome data set of the Mycobacterium tuberculosis complex (MTBC) gathered to date. The closed pan-genome of the MTBC was characterized by reduced accessory and strain-specific genomes, compatible with its clonal nature. However, significantly fewer gene families were shared between MTBC genomes as their phylogenetic distance increased. This effect was only observed in inter-species comparisons, not within-species, which suggests that species-specific ecological characteristics are associated with changes in gene content. Gene loss, resulting from genomic deletions and pseudogenization, was found to drive the variation in gene content. This gene erosion differed among MTBC species and lineages, even within M. tuberculosis, where L2 showed more gene loss than L4. We also show that phylogenetic proximity is not always a good proxy for gene content relatedness in the MTBC, as the gene repertoire of Mycobacterium africanum L6 deviated from its expected phylogenetic niche conservatism. Gene disruptions of virulence factors, represented by pseudogene annotations, are mostly not conserved, being poor predictors of MTBC ecotypes. Each MTBC ecotype carries its own accessory genome, likely influenced by distinct selective pressures such as host and geography. It is important to investigate how gene loss confer new adaptive traits to MTBC strains; the detected heterogeneous gene loss poses a significant challenge in elucidating genetic factors responsible for the diverse phenotypes observed in the MTBC. By detailing specific gene losses, our study serves as a resource for researchers studying the MTBC phenotypes and their immune evasion strategies.IMPORTANCEIn this study, we analyzed the gene content of different ecotypes of the Mycobacterium tuberculosis complex (MTBC), the pathogens of tuberculosis. We found that changes in their gene content are associated with their ecological features, such as host preference. Gene loss was identified as the primary driver of these changes, which can vary even among different strains of the same ecotype. Our study also revealed that the gene content relatedness of these bacteria does not always mirror their evolutionary relationships. In addition, some genes of virulence can be variably lost among strains of the same MTBC ecotype, likely helping them to evade the immune system. Overall, our study highlights the importance of understanding how gene loss can lead to new adaptations in these bacteria and how different selective pressures may influence their genetic makeup.

Host cell transcriptomic response to the multidrug-resistant Mycobacterium tuberculosis clonal outbreak Beijing strain reveals its pathogenic features

The upsurge of multidrug-resistant infections has rendered tuberculosis the principal cause of death among infectious diseases. A clonal outbreak multidrug-resistant triggering strain of Mycobacterium tuberculosis was identified in Kanchanaburi Province, labelled "MKR superspreader," which was found to subsequently spread to other regions, as revealed by prior epidemiological reports in Thailand. Herein, we showed that the MKR displayed a higher growth rate upon infection into host macrophages in comparison with the H37Rv reference strain. To further elucidate MKR's biology, we utilized RNA-Seq and differential gene expression analyses to identify host factors involved in the intracellular viability of the MKR. A set of host genes function in the cellular response to lipid pathway was found to be uniquely up-regulated in host macrophages infected with the MKR, but not those infected with H37Rv. Within this set of genes, the IL-36 cytokines which regulate host cell cholesterol metabolism and resistance against mycobacteria attracted our interest, as our previous study revealed that the MKR elevated genes associated with cholesterol breakdown during its growth inside host macrophages. Indeed, when comparing macrophages infected with the MKR to H37Rv-infected cells, our RNA-Seq data showed that the expression ratio of IL-36RN, the negative regulator of the IL-36 pathway, to that of IL-36G was greater in macrophages infected with the MKR. Furthermore, the MKR's intracellular survival and increased intracellular cholesterol level in the MKR-infected macrophages were diminished with decreased IL-36RN expression. Overall, our results indicated that IL-36RN could serve as a new target against this emerging multidrug-resistant M. tuberculosis strain.

Rapid Identification of Drug-Resistant Tuberculosis Genes Using Direct PCR Amplification and Oxford Nanopore Technology Sequencing

Mycobacterium tuberculosis antimicrobial resistance has been continually reported and is a major public health issue worldwide. Rapid prediction of drug resistance is important for selecting appropriate antibiotic treatments, which significantly increases cure rates. Gene sequencing technology has proven to be a powerful strategy for identifying relevant drug resistance information. This study established a sequencing method and bioinformatics pipeline for resistance gene analysis using an Oxford Nanopore Technologies sequencer. The pipeline was validated by Sanger sequencing and exhibited 100% concordance with the identified variants. Turnaround time for the nanopore sequencing workflow was approximately 12 h, facilitating drug resistance prediction several weeks earlier than that of traditional phenotype drug susceptibility testing. This study produced a customized gene panel assay for rapid bacterial identification via nanopore sequencing, which improves the timeliness of tuberculosis diagnoses and provides a reliable method that may have clinical application.

Alterations in molecular response of Mycobacterium tuberculosis against anti-tuberculosis drugs

BACKGROUND

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis, has plagued humans since the early middle-ages. More than one million deaths are recorded annually due to TB, even in present times. These deaths are primarily attributed to the constant appearance of resistant TB strains. Even with the advent of new therapeutics and diagnostics techniques, tuberculosis remains challenging to control due to resistant M. tuberculosis strains. Aided by various molecular changes, these strains adapt to stress created by anti-tuberculosis drugs.

MATERIALS AND METHODS

The review thus is an overview of ongoing research in the genome and transcriptome of antibiotic-resistant TB. It explores omics-based research to identify mutation and utilization of differential gene expression.

CONCLUSIONS

This study shows several mutations distinctive in the first- and second-line drug-resistant M. tuberculosis strains. It also explores the expressional differences of genes involved in the fundamental process of the cells and how they help in drug resistance. With the development of transcriptomics-based studies, a new insight has developed to inquire about gene expression changes in drug resistance. This information on expressional pattern changes can be utilized to design the basic platform of anti-TB treatments and therapeutic approaches. These novel insights can be instrumental in disease diagnosis and global containment of resistant TB.

Revised nomenclature and SNP barcode for Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis (Mtb) lineage 2 (L2) strains are present globally, contributing to a widespread tuberculosis (TB) burden, particularly in Asia where both prevalence of TB and numbers of drug resistant TB are highest. The increasing availability of whole-genome sequencing (WGS) data worldwide provides an opportunity to improve our understanding of the global genetic diversity of Mtb L2 and its association with the disease epidemiology and pathogenesis. However, existing L2 sublineage classification schemes leave >20 % of the Modern Beijing isolates unclassified. Here, we present a revised SNP-based classification scheme of L2 in a genomic framework based on phylogenetic analysis of >4000 L2 isolates from 34 countries in Asia, Eastern Europe, Oceania and Africa. Our scheme consists of over 30 genotypes, many of which have not been described before. In particular, we propose six main genotypes of Modern Beijing strains, denoted L2.2.M1–L2.2.M6. We also provide SNP markers for genotyping L2 strains from WGS data. This fine-scale genotyping scheme, which can classify >98 % of the studied isolates, serves as a basis for more effective monitoring and reporting of transmission and outbreaks, as well as improving genotype-phenotype associations such as disease severity and drug resistance. This article contains data hosted by Microreact.

Drug-resistant Mycobacterium tuberculosis and its genotypes isolated from an outbreak in western Thailand

BACKGROUND

Multidrug-resistant TB (MDR-TB) outbreaks have occurred in the Thamaka district, Kanchanaburi province in Thailand.

METHODS

Seventy-two isolates, which included 7% mono-, 30.6% MDR and extensively drug-resistant TB (XDR-TB), were genotyped by spoligotyping, mycobacterial interspersed repetitive unit-variable-number tandem repeat (MIRU-VNTR) and single nucleotide polymorphism genotyping, and their drug resistance was analysed.

RESULTS

The spoligotyping results showed that Beijing spoligo-international type (SIT)1 was predominant (n=38; 52.8%) while the remaining were non-Beijing sublineages (n=34). The MIRU-VNTR analysis showed that Beijing isolates, most of which belonged to the modern type (n=37), formed 5 clusters and 13 individual patterns. In katG, only mutation Ser315Thr was identified. In rpoB, Ser531Leu was predominant, except for His526Arg and Leu533Pro, which were found in two isolates. A cluster of 14 Beijing strains contained these common mutations and shared the MIRU-VNTR genotype with isolates in the Thamaka district that had spread previously. Two U SIT523 isolates contained the mutations A1400G in rrs and Asp94Gly in gyrA genes, indicating a spread of XDR-TB.

CONCLUSIONS

Most mutations were associated with drug resistance and the specific MDR Beijing and XDR-TB in U SIT523 isolates remain. This genotyping is a key tool for tracking TB transmission in the Thamaka district of Thailand.

Clusters of Drug-Resistant Mycobacterium tuberculosis Detected by Whole-Genome Sequence Analysis of Nationwide Sample, Thailand, 2014-2017

Multidrug-resistant tuberculosis (MDR TB), pre-extensively drug-resistant tuberculosis (pre-XDR TB), and extensively drug-resistant tuberculosis (XDR TB) complicate disease control. We analyzed whole-genome sequence data for 579 phenotypically drug-resistant M. tuberculosis isolates (28% of available MDR/pre-XDR and all culturable XDR TB isolates collected in Thailand during 2014–2017). Most isolates were from lineage 2 (n = 482; 83.2%). Cluster analysis revealed that 281/579 isolates (48.5%) formed 89 clusters, including 205 MDR TB, 46 pre-XDR TB, 19 XDR TB, and 11 poly–drug-resistant TB isolates based on genotypic drug resistance. Members of most clusters had the same subset of drug resistance-associated mutations, supporting potential primary resistance in MDR TB (n = 176/205; 85.9%), pre-XDR TB (n = 29/46; 63.0%), and XDR TB (n = 14/19; 73.7%). Thirteen major clades were significantly associated with geography (p<0.001). Clusters of clonal origin contribute greatly to the high prevalence of drug-resistant TB in Thailand.

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.

Lysosome repositioning as an autophagy escape mechanism by Mycobacterium tuberculosis Beijing strain

Induction of host cell autophagy by starvation was shown to enhance lysosomal delivery to mycobacterial phagosomes, resulting in the restriction of Mycobacterium tuberculosis reference strain H37Rv. Our previous study showed that strains belonging to M. tuberculosis Beijing genotype resisted starvation-induced autophagic elimination but the factors involved remained unclear. Here, we conducted RNA-Seq of macrophages infected with the autophagy-resistant Beijing strain (BJN) compared to macrophages infected with H37Rv upon autophagy induction by starvation. Results identified several genes uniquely upregulated in BJN-infected macrophages but not in H37Rv-infected cells, including those encoding Kxd1 and Plekhm2, which function in lysosome positioning towards the cell periphery. Unlike H37Rv, BJN suppressed enhanced lysosome positioning towards the perinuclear region and lysosomal delivery to its phagosome upon autophagy induction by starvation, while depletion of Kxd1 and Plekhm2 reverted such effects, resulting in restriction of BJN intracellular survival upon autophagy induction by starvation. Taken together, these data indicated that Kxd1 and Plekhm2 are important for the BJN strain to suppress lysosome positioning towards the perinuclear region and lysosomal delivery into its phagosome during autophagy induction by starvation to evade starvation-induced autophagic restriction.

Transcriptional response to the host cell environment of a multidrug-resistant Mycobacterium tuberculosis clonal outbreak Beijing strain reveals its pathogenic features

Tuberculosis is a global public health problem with emergence of multidrug-resistant infections. Previous epidemiological studies of tuberculosis in Thailand have identified a clonal outbreak multidrug-resistant strain of Mycobacterium tuberculosis in the Kanchanaburi province, designated “MKR superspreader”, and this particular strain later was found to also spread to other regions. In this study, we elucidated its biology through RNA-Seq analyses and identified a set of genes involved in cholesterol degradation to be up-regulated in the MKR during the macrophage cell infection, but not in the H37Rv reference strain. We also found that the bacterium up-regulated genes associated with the ESX-1 secretion system during its intracellular growth phase, while the H37Rv did not. All results were confirmed by qRT-PCR. Moreover, we showed that compounds previously shown to inhibit the mycobacterial ESX-1 secretion system and cholesterol utilisation, and FDA-approved drugs known to interfere with the host cholesterol transportation were able to decrease the intracellular survival of the MKR when compared to the untreated control, while not that of the H37Rv. Altogether, our findings suggested that such pathways are important for the MKR’s intracellular growth, and potentially could be targets for the discovery of new drugs against this emerging multidrug-resistant strain of M. tuberculosis.

Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes.IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

Genomic evidence supporting the clonal expansion of extensively drug-resistant tuberculosis bacteria belonging to a rare protoBeijing genotype

Tuberculosis disease (TB), caused by Mycobacterium tuberculosis, is a major public health issue in Thailand. The high prevalence of modern Beijing (Lineage 2.2.1) strains has been associated with multi- and extensively drug-resistant infections (MDR-, XDR-TB), complicating disease control. The impact of rarer proto-Beijing (L2.1) strains is less clear. In our study of thirty-seven L2.1 clinical isolates spanning thirteen years, we found a high prevalence of XDR-TB cases (32.4%). With ≤ 12 pairwise SNP distances, 43.2% of L2.1 patients belong to MDR-TB or XDR-TB transmission clusters suggesting a high level of clonal expansion across four Thai provinces. All XDR-TB (100%) were likely due to transmission rather than inadequate treatment. We found a 47 mutation signature and a partial deletion of the fadD14 gene in the circulating XDR-TB cluster, which can be used for surveillance of this rare and resilient M. tuberculosis strain-type that is causing increasing health burden. We also detected three novel deletion positions, a deletion of 1285 bp within desA3 (Rv3230c), large deletions in the plcB, plcA, and ppe38 gene which may play a role in the virulence, pathogenesis or evolution of the L2.1 strain-type.

Geno-Spatial Distribution of Mycobacterium Tuberculosis and Drug Resistance Profiles in Myanmar-Thai Border Area

Worldwide, studies investigating the relationship between the lineage of Mycobacterium tuberculosis (MTB) across geographic areas has empowered the "End TB" program and understand transmission across national boundaries. Genomic diversity of MTB varies with geographical locations and ethnicity. Genomic diversity can also affect the emergence of drug resistance. In Myanmar, we still have limited genetic information about geographical, ethnicity, and drug resistance linkage to MTB genetic information. This study aimed to describe the geno-spatial distribution of MTB and drug resistance profiles in Myanmar-Thailand border areas. A cross-sectional study was conducted with a total of 109 sequenced isolates. The lineages of MTB and the potential associated socio-demographic, geographic and clinical factors were analyzed using Fisher's exact tests. p value of statistically significance was set at < 0.05. We found that 67% of the isolates were lineage 1 (L1)/East-African-Indian (EAI) (n = 73), followed by lineage 2 (L2)/Beijing (n = 26), lineage 4 (L4)/European American (n = 6) and lineage 3 (L3)/Delhi/Central Asian (n = 4). "Gender", "type of TB patient", "sputum smear grading" and "streptomycin resistance" were significantly different with the lineages of MTB. Sublineages of L1, which had never been reported elsewhere in Myanmar, were detected in this study area. Moreover, both ethnicity and lineage of MTB significantly differed in distribution by patient location. Diversity of the lineage of MTB and detection of new sublineages suggested that this small area had been resided by a heterogeneous population group who actively transmitted the disease. This information on distribution of lineage of MTB can be linked in the future with those on the other side of the border to evaluate cross-border transmission.

A brief primer on genomic epidemiology: lessons learned from Mycobacterium tuberculosis

Genomics is now firmly established as a technique for the investigation and reconstruction of communicable disease outbreaks, with many genomic epidemiology studies focusing on revealing transmission routes of Mycobacterium tuberculosis. In this primer, we introduce the basic techniques underlying transmission inference from genomic data, using illustrative examples from M. tuberculosis and other pathogens routinely sequenced by public health agencies. We describe the laboratory and epidemiological scenarios under which genomics may or may not be used, provide an introduction to sequencing technologies and bioinformatics approaches to identifying transmission-informative variation and resistance-associated mutations, and discuss how variation must be considered in the light of available clinical and epidemiological information to infer transmission.

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.

Bioinformatics tools and databases for whole genome sequence analysis of Mycobacterium tuberculosis

Tuberculosis (TB) is an infectious disease of global public health importance caused by Mycobacterium tuberculosis complex (MTC) in which M. tuberculosis (Mtb) is the major causative agent. Recent advancements in genomic technologies such as next generation sequencing have enabled high throughput cost-effective generation of whole genome sequence information from Mtb clinical isolates, providing new insights into the evolution, genomic diversity and transmission of the Mtb bacteria, including molecular mechanisms of antibiotic resistance. The large volume of sequencing data generated however necessitated effective and efficient management, storage, analysis and visualization of the data and results through development of novel and customized bioinformatics software tools and databases. In this review, we aim to provide a comprehensive survey of the current freely available bioinformatics software tools and publicly accessible databases for genomic analysis of Mtb for identifying disease transmission in molecular epidemiology and in rapid determination of the antibiotic profiles of clinical isolates for prompt and optimal patient treatment.

Whole-Genome Sequencing Analysis of Serially Isolated Multi-Drug and Extensively Drug Resistant Mycobacterium tuberculosis from Thai Patients

Multi-drug and extensively drug-resistant tuberculosis (MDR and XDR-TB) are problems that threaten public health worldwide. Only some genetic markers associated with drug-resistant TB are known. Whole-genome sequencing (WGS) is a promising tool for distinguishing between re-infection and persistent infection in isolates taken at different times from a single patient, but has not yet been applied in MDR and XDR-TB. We aim to detect genetic markers associated with drug resistance and distinguish between reinfection and persistent infection from MDR and XDR-TB patients based on WGS analysis. Samples of Mycobacterium tuberculosis (n = 7), serially isolated from 2 MDR cases and 1 XDR-TB case, were retrieved from Siriraj Hospital, Bangkok. The WGS analysis used an Illumina Miseq sequencer. In cases of persistent infection, MDR-TB isolates differed at an average of 2 SNPs across the span of 2–9 months whereas in the case of reinfection, isolates differed at 61 SNPs across 2 years. Known genetic markers associated with resistance were detected from strains susceptible to streptomycin (2/7 isolates), p-aminosalicylic acid (3/7 isolates) and fluoroquinolone drugs. Among fluoroquinolone drugs, ofloxacin had the highest phenotype-genotype concordance (6/7 isolates), whereas gatifloxcain had the lowest (3/7 isolates). A putative candidate SNP in Rv2477c associated with kanamycin and amikacin resistance was suggested for further validation. WGS provided comprehensive results regarding molecular epidemiology, distinguishing between persistent infection and reinfection in M/XDR-TB and potentially can be used for detection of novel mutations associated with drug resistance.

Interpreting whole genome sequencing for investigating tuberculosis transmission: a systematic review

BACKGROUND

Whole genome sequencing (WGS) is becoming an important part of epidemiological investigations of infectious diseases due to greater resolution and cost reductions compared to traditional typing approaches. Many public health and clinical teams will increasingly use WGS to investigate clusters of potential pathogen transmission, making it crucial to understand the benefits and assumptions of the analytical methods for investigating the data. We aimed to understand how different approaches affect inferences of transmission dynamics and outline limitations of the methods.

METHODS

We comprehensively searched electronic databases for studies that presented methods used to interpret WGS data for investigating tuberculosis (TB) transmission. Two authors independently selected studies for inclusion and extracted data. Due to considerable methodological heterogeneity between studies, we present summary data with accompanying narrative synthesis rather than pooled analyses.

RESULTS

Twenty-five studies met our inclusion criteria. Despite the range of interpretation tools, the usefulness of WGS data in understanding TB transmission often depends on the amount of genetic diversity in the setting. Where diversity is small, distinguishing re-infections from relapses may be impossible; interpretation may be aided by the use of epidemiological data, examining minor variants and deep sequencing. Conversely, when within-host diversity is large, due to genetic hitchhiking or co-infection of two dissimilar strains, it is critical to understand how it arose. Greater understanding of microevolution and mixed infection will enhance interpretation of WGS data.

CONCLUSIONS

As sequencing studies have sampled more intensely and integrated multiple sources of information, the understanding of TB transmission and diversity has grown, but there is still much to be learnt about the origins of diversity that will affect inferences from these data. Public health teams and researchers should combine epidemiological, clinical and WGS data to strengthen investigations of transmission.

Evidence for the critical role of a secondary site rpoB mutation in the compensatory evolution and successful transmission of an MDR tuberculosis outbreak strain

BACKGROUND

MDR Mycobacterium tuberculosis clinical strains that cause large outbreaks, particularly among HIV-negative patients, are likely to have undergone the most successful compensatory evolution. Hence, mutations secondary to the acquisition of drug resistance are worthy of consideration in these highly transmissible strains. Here, we assessed the role of a mutation within rpoB, rpoB V615M, secondary to the rifampicin resistance-conferring mutation rpoB S531L, which is associated with a major MDR tuberculosis outbreak strain that evolved in an HIV-negative context in northern Tunisia.

METHODS

Using BCG as a model organism, we engineered strains harbouring either the rpoB S531L mutation alone or the double mutation rpoB S531L, V615M. Individual and competitive in vitro growth assays were performed in order to assess the relative fitness of each BCG mutant.

RESULTS

The rpoB V615M mutation was found to be invariably associated with rpoB S531L. Structural analysis mapped rpoB V615M to the same bridge helix region as rpoB compensatory mutations previously described in Salmonella. Compared with the rpoB single-mutant BCG, the double mutant displayed improved growth characteristics and fitness rates equivalent to WT BCG. Strikingly, the rpoB double mutation conferred high-level resistance to rifampicin.

CONCLUSIONS

Here, we demonstrated the fitness compensatory role of a mutation within rpoB, secondary to the rifampicin resistance mutation rpoB S531L, which is characteristic of an MDR M. tuberculosis major outbreak strain. The finding that this secondary mutation concomitantly increased the resistance level to rifampicin argues for its significant contribution to the successful transmission of the MDR-TB strain.