Whole genome sequencing reveals complex evolution patterns of multidrug-resistant Mycobacterium tuberculosis Beijing strains in patients

A study on factors influencing delayed sputum conversion in newly diagnosed pulmonary tuberculosis based on bacteriology and genomics

Conversion of sputum from positive to negative is one of the indicators to evaluate the efficacy of anti-tuberculosis treatment (ATT). We investigate the factors associated with delayed sputum conversion after 2 or 5 months of ATT from the perspectives of bacteriology and genomics. A retrospective study of sputum conversion in sputum positive 1782 pulmonary tuberculosis (PTB) was conducted from 2021 to 2022 in Beijing, China. We also designed a case-matched study including 24 pairs of delayed-sputum-conversion patients (DSCPs) and timely-sputum-conversion patients (TSCPs), and collect clinical isolates from DSCPs before and after ATT and initial isolates of TSCPs who successfully achieved sputum conversion to negative after 2 months of ATT. A total of 75 strains were conducted drug sensitivity testing (DST) of 13 anti-TB drugs and whole-genome sequencing (WGS) to analyze the risk factors of delayed conversion and the dynamics changes of drug resistance and genomics of Mycobacterium tuberculosis (MTB) during ATT. We found TSCPs have better treatment outcomes and whose initial isolates show lower levels of drug resistance. Clinical isolates of DSCPs showed dynamically changing of resistance phenotypes and intra-host heterogeneity. Single nucleotide polymorphism (SNP) profiles showed large differences between groups. The study provided insight into the bacteriological and genomic variation of delayed sputum conversion. It would be helpful for early indication of sputum conversion and guidance on ATT.

SNPs in genes related to the repair of damage to DNA in clinical isolates of M. tuberculosis: A transversal and longitudinal approach

The presence of SNPs in genes related to DNA damage repair in M. tuberculosis can trigger hypermutagenic phenotypes with a higher probability of generating drug resistance. The aim of this research was to compare the presence of SNPs in genes related to DNA damage repair between sensitive and DR isolates, as well as to describe the dynamics in the presence of SNPs in M. tuberculosis isolated from recently diagnosed TB patients of the state of Veracruz, Mexico. The presence of SNPs in the coding regions of 65 genes related to DNA damage repair was analyzed. Eighty-six isolates from 67 patients from central Veracruz state, Mexico, were sequenced. The results showed several SNPs in 14 genes that were only present in drug-resistant genomes. In addition, by following of 15 patients, it was possible to describe three different dynamics of appearance and evolution of non-synonymous SNPs in genes related to DNA damage repair: 1) constant fixed SNPs, 2) population substitution, and 3) gain of fixed SNPs. Further research is required to discern the biological significance of each of these pathways and their utility as markers of DR or for treatment prognosis.

Association between mutations in a thyX-hsdS.1 region and para-aminosalicylic acid resistance in Mycobacterium tuberculosis clinical isolates

Although para-aminosalicylic acid (PAS) has been used to treat tuberculosis agent for decades, its mechanisms of resistance are still not thoroughly understood. Previously, sporadic studies showed that certain mutations in the thyX-hsdS.1 region caused PAS resistance in M. tuberculosis, but a comprehensive analysis is lacking. Recently, we found a G-10A mutation in thyX-hsdS.1 in a PAS-resistant clinical isolate, but it did not cause PAS resistance. SNPs in thyX-hsdS.1 in 6550 clinical isolates were analyzed, and 153 SNPs were identified. C-16 T was the most common SNP identified (54.25%, 83/153), followed by C-4T (7.19%, 11/153) and G-9A (6.54%, 10/153). Subsequently, the effects of those SNPs on the promoter activity of thyX were tested, and the results showed that mutations C-1T, G-3A, C-4T, C-4G, G-7A, G-9A, C-16T, G-18C, and C-19G led to increased promoter activity compared with the wild-type sequence, but other mutations did not. Then, thyX and wild-type thyX-hsdS.1, or thyX-hsdS.1 containing specific SNPs, were overexpressed in M. tuberculosis H37Ra. The results showed that mutations resulting in increased promoter activity also caused PAS resistance. Moreover, the results of an electrophoretic mobility shift assay showed that thyX-hsdS.1 containing the C-16T mutation had a higher binding capacity to RNA polymerase than did the wild-type sequence. Taken together, our data demonstrated that among the SNPs identified in thyX-hsdS.1 of M. tuberculosis clinical isolates, only those able to increase the promoter activity of thyX caused PAS resistance and therefore can be considered as molecular markers for PAS resistance.

Transmission Dynamics of a Mycobacterium tuberculosis Complex Outbreak in an Indigenous Population in the Colombian Amazon Region

Whole genome sequencing (WGS) has become the main tool for studying the transmission of Mycobacterium tuberculosis complex (MTBC) strains; however, the clonal expansion of one strain often limits its application in local MTBC outbreaks. The use of an alternative reference genome and the inclusion of repetitive regions in the analysis could potentially increase the resolution, but the added value has not yet been defined. Here, we leveraged short and long WGS read data of a previously reported MTBC outbreak in the Colombian Amazon Region to analyze possible transmission chains among 74 patients in the indigenous setting of Puerto Nariño (March to October 2016). In total, 90.5% (67/74) of the patients were infected with one distinct MTBC strain belonging to lineage 4.3.3. Employing a reference genome from an outbreak strain and highly confident single nucleotide polymorphisms (SNPs) in repetitive genomic regions, e.g., the proline-glutamic acid/proline-proline-glutamic-acid (PE/PPE) gene family, increased the phylogenetic resolution compared to a classical H37Rv reference mapping approach. Specifically, the number of differentiating SNPs increased from 890 to 1,094, which resulted in a more granular transmission network as judged by an increasing number of individual nodes in a maximum parsimony tree, i.e., 5 versus 9 nodes. We also found in 29.9% (20/67) of the outbreak isolates, heterogenous alleles at phylogenetically informative sites, suggesting that these patients are infected with more than one clone. In conclusion, customized SNP calling thresholds and employment of a local reference genome for a mapping approach can improve the phylogenetic resolution in highly clonal MTBC populations and help elucidate within-host MTBC diversity. IMPORTANCE The Colombian Amazon around Puerto Nariño has a high tuberculosis burden with a prevalence of 1,267/100,000 people in 2016. Recently, an outbreak of Mycobacterium tuberculosis complex (MTBC) bacteria among the indigenous populations was identified with classical MTBC genotyping methods. Here, we employed a whole-genome sequencing-based outbreak investigation in order to improve the phylogenetic resolution and gain new insights into the transmission dynamics in this remote Colombian Amazon Region. The inclusion of well-supported single nucleotide polymorphisms in repetitive regions and a de novo-assembled local reference genome provided a more granular picture of the circulating outbreak strain and revealed new transmission chains. Multiple patients from different settlements were possibly infected with at least two different clones in this high-incidence setting. Thus, our results have the potential to improve molecular surveillance studies in other high-burden settings, especially regions with few clonal multidrug-resistant (MDR) MTBC lineages/clades.

Mycobacterium tuberculosis Intra-Host Evolution Among Drug-Resistant Tuberculosis Patients Failing Treatment

Background

Understanding Mycobacterium tuberculosis (Mtb) intra-host evolution of drug resistance is important for successful drug-resistant tuberculosis (DR-TB) treatment and control strategies. This study aimed to characterise the acquisition of genetic mutations and low-frequency variants associated with treatment-emergent Mtb drug resistance in longitudinally profiled clinical isolates from patients who experienced DR-TB treatment failure.

Patients and Methods

We performed deep Whole Genome Sequencing on 23 clinical isolates obtained longitudinally across nine timepoints from five patients who experienced DR-TB treatment failure enrolled in the CAPRISA 020 InDEX study. The minimum inhibitory concentrations (MICs) were established on the BACTEC™ MGIT 960™ instrument on 15/23 longitudinal clinical isolates for eight anti-TB drugs (rifampicin, isoniazid, ethambutol, levofloxacin, moxifloxacin, linezolid, clofazimine, bedaquiline).

Results

In total, 22 resistance associated mutations/variants were detected. We observed four treatment-emergent mutations in two out of the five patients. Emerging resistance to the fluoroquinolones was associated with 16- and 64-fold elevated levofloxacin (2-8 mg/L) and moxifloxacin (1-2 mg/L) MICs, respectively, resulting from the D94G/N and A90V variants in the gyrA gene. We identified two novel mutations associated with elevated bedaquiline MICs (>66-fold): an emerging frameshift variant (D165) on the Rv0678 gene and R409Q variant on the Rv1979c gene present from baseline.

Conclusion

Genotypic and phenotypic resistance to the fluoroquinolones and bedaquiline was acquired in two out of five patients who experienced DR-TB treatment failure. Deep sequencing of multiple longitudinal clinical isolates for resistance-associated mutations coupled with phenotypic MIC testing confirmed intra-host Mtb evolution.

Clinical implications of molecular drug resistance testing for Mycobacterium tuberculosis: a 2023 TBnet/RESIST-TB consensus statement

Drug-resistant tuberculosis is a substantial health-care concern worldwide. Despite culture-based methods being considered the gold standard for drug susceptibility testing, molecular methods provide rapid information about the Mycobacterium tuberculosis mutations associated with resistance to anti-tuberculosis drugs. This consensus document was developed on the basis of a comprehensive literature search, by the TBnet and RESIST-TB networks, about reporting standards for the clinical use of molecular drug susceptibility testing. Review and the search for evidence included hand-searching journals and searching electronic databases. The panel identified studies that linked mutations in genomic regions of M tuberculosis with treatment outcome data. Implementation of molecular testing for the prediction of drug resistance in M tuberculosis is key. Detection of mutations in clinical isolates has implications for the clinical management of patients with multidrug-resistant or rifampicin-resistant tuberculosis, especially in situations when phenotypic drug susceptibility testing is not available. A multidisciplinary team including clinicians, microbiologists, and laboratory scientists reached a consensus on key questions relevant to molecular prediction of drug susceptibility or resistance to M tuberculosis, and their implications for clinical practice. This consensus document should help clinicians in the management of patients with tuberculosis, providing guidance for the design of treatment regimens and optimising outcomes.

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.

Genome-Wide Study of Drug Resistant Mycobacterium tuberculosis and Its Intra-Host Evolution during Treatment

The emergence of drug resistant Mycobacterium tuberculosis (MTB) strains has become a global public health problem, while, at the same time, there has been development of new antimicrobial agents. The main goals of this study were to determine new variants associated with drug resistance in MTB and to observe which polymorphisms emerge in MTB genomes after anti-tuberculosis treatment. We performed whole-genome sequencing of 152 MTB isolates including 70 isolates as 32 series of pre- and post-treatment MTB. Based on genotypes and phenotypic drug susceptibility, we conducted phylogenetic convergence-based genome-wide association study (GWAS) with streptomycin-, isoniazid-, rifampicin-, ethambutol-, fluoroquinolones-, and aminoglycosides-resistant MTB against susceptible ones. GWAS revealed statistically significant associations of SNPs within Rv2820c, cyp123 and indels in Rv1269c, Rv1907c, Rv1883c, Rv2407, Rv3785 genes with resistant MTB phenotypes. Comparisons of serial isolates showed that treatment induced different patterns of intra-host evolution. We found indels within Rv1435c and ppsA that were not lineage-specific. In addition, Beijing-specific polymorphisms within Rv0036c, Rv0678, Rv3433c, and dop genes were detected in post-treatment isolates. The appearance of Rv3785 frameshift insertion in 2 post-treatment strains compared to pre-treatment was also observed. We propose that the insertion within Rv3785, which was a GWAS hit, might affect cell wall biosynthesis and probably mediates a compensatory mechanism in response to treatment. These results may shed light on the mechanisms of MTB adaptation to chemotherapy and drug resistance formation.

IS6110 Copy Number in Multi-Host Mycobacterium bovis Strains Circulating in Bovine Tuberculosis Endemic French Regions

IS6110 is an insertion sequence found in the Mycobacterium tuberculosis complex, to which Mycobacterium bovis belongs, which can play a role in genome plasticity and in bacterial evolution. In this study, the abundance and location of IS6110 on M. bovis genomic data of French animal field strains were studied. A first analysis was performed on a panel of 81 strains that reflect the national M. bovis population’s genetic diversity. The results show that more than one-third of them are IS6110 multicopy and that 10% have IS6110 in a high copy number (more than 6 copies). Multicopy strains are those circulating in the regions where prevalence was above the national average. Further study of 93 such strains, with an IS6110 copy number of 10-12, showed stability of IS6110 copy number and genome location over time and between host species. The correlation between M. bovis multicopy strains and high bovine tuberculosis (bTB) prevalence leads us to consider whether their epidemiological success could be partly due to genetic changes originated by IS6110 transposition.

Prediction of Mycobacterium tuberculosis drug resistance by nucleotide MALDI-TOF-MS

OBJECTIVES

To evaluate the performance of nucleotide matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) in predicting the drug resistance of Mycobacterium tuberculosis.

METHODS

A total of 115 rifampin-resistant and 53 rifampin-susceptible tuberculosis (TB) clinical isolates were randomly selected from TB strains stored at -80℃ in the clinical laboratory of Shanghai Pulmonary Hospital. Nucleotide MALDI-TOF-MS was performed to predict the drug resistance using phenotypic susceptibility as the gold standard.

RESULTS

The overall assay sensitivities and specificities of nucleotide MALDI-TOF-MS were 92.2% and 100.0% for rifampin, 90.9% and 98.6% for isoniazid, 71.4% and 81.2% for ethambutol, 85.1% and 93.1% for streptomycin, 94.0% and 100.0% for amikacin, 77.8% and 99.3% for kanamycin, 75.0% and 93.3% for ofloxacin, and 75.0% and 93.3% for moxifloxacin. The concordances between nucleotide MALDI-TOF-MS antimicrobial susceptibility testing (AST) and phenotypic AST were 94.6% (rifampin), 90.1% (isoniazid), 79.2% (ethambutol), 89.9% (streptomycin), 99.4% (amikacin), 97.0% (kanamycin), 88.1% (ofloxacin), and 88.0% (moxifloxacin).

CONCLUSION

Nucleotide MALDI-TOF-MS could be a promising tool for rapid detection of Mycobacterium tuberculosis drug sensitivity to rifampin, isoniazid, ethambutol, streptomycin, amikacin, kanamycin, ofloxacin, and moxifloxacin.

Beijing genotype of Mycobacterium tuberculosis is associated with extensively drug-resistant tuberculosis: A global analysis

We found that the frequency of Beijing genotype among XDR-TB strains was high. The data in this study would help guide the TB control program, and we however need further investigation to confirm the reliability of the present findings.

Single-nucleotide polymorphism-based epidemiological analysis of Korean Mycobacterium bovis isolates

Background

Bovine tuberculosis (TB) is caused by Mycobacterium bovis, a well-known cause of zoonotic tuberculosis in cattle and deer, and has been investigated in many physiological and molecular studies. However, detailed genome-level studies of M. bovis have not been performed in Korea.

Objectives

To survey whole genome-wide single-nucleotide polymorphism (SNP) variants in Korean M. bovis field isolates and to define M. bovis groups in Korea by comparing SNP typing with spoligotyping and variable number tandem repeat typing.

Methods

A total of 46 M. bovis field isolates, isolated from laryngopharyngeal lymph nodes and lungs of Korean cattle, wild boar, and Korean water deer, were used to identify SNPs by performing whole-genome sequencing. SNP sites were confirmed via polymerase chain reaction using 87 primer pairs.

Results

We identified 34 SNP sites with different frequencies across M. bovis isolates, and performed SNP typing and epidemiological analysis, which divided the 46 field isolates into 16 subtypes.

Conclusions

Through SNP analysis, detailed differences in samples with identical spoligotypes could be detected. SNP analysis is, therefore, a useful epidemiological tracing tool that could enable better management of bovine TB, thus preventing further outbreaks and reducing the impact of this disease.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

Selective multiplexed enrichment for the detection and quantitation of low-fraction DNA variants via low-depth sequencing

DNA sequence variants with allele fractions below 1% are difficult to detect and quantify by sequencing owing to intrinsic errors in sequencing-by-synthesis methods. Although molecular-identifier barcodes can detect mutations with a variant-allele frequency (VAF) as low as 0.1% using next-generation sequencing (NGS), sequencing depths of over 25,000× are required, thus hampering the detection of mutations at high sensitivity in patient samples and in most samples used in research. Here we show that low-frequency DNA variants can be detected via low-depth multiplexed NGS after their amplification, by a median of 300-fold, using polymerase chain reaction and rationally designed 'blocker' oligonucleotides that bind to the variants. Using an 80-plex NGS panel and a sequencing depth of 250×, we detected single nucleotide polymorphisms with a VAF of 0.019% and contamination in human cell lines at a VAF as low as 0.07%. With a 16-plex NGS panel covering 145 mutations across 9 genes involved in melanoma, we detected low-VAF mutations (0.2-5%) in 7 out of the 19 samples of freshly frozen tumour biopsies, suggesting that tumour heterogeneity could be notably higher than previously recognized.

Microevolution of Mycobacterium tuberculosis Subpopulations and Heteroresistance in a Patient Receiving 27 Years of Tuberculosis Treatment in Germany

Preexisting and newly emerging resistant pathogen subpopulations (heteroresistance) are potential risk factors for treatment failure of multi/extensively drug resistant (MDR/XDR) tuberculosis (TB). Intrapatient evolutionary dynamics of Mycobacterium tuberculosis complex (Mtbc) strains and their implications on treatment outcomes are still not completely understood. To elucidate how Mtbc strains escape therapy, we analyzed 13 serial isolates from a German patient by whole-genome sequencing. Sequencing data were compared with phenotypic drug susceptibility profiles and the patient's collective 27-year treatment history to further elucidate factors fostering intrapatient resistance evolution. The patient endured five distinct TB episodes, ending in resistance to 16 drugs and a nearly untreatable XDR-TB infection. The first isolate obtained, during the patient's 5th TB episode, presented fixed resistance mutations to 7 anti-TB drugs, including isoniazid, rifampin, streptomycin, pyrazinamide, prothionamide, para-aminosalicylic acid, and cycloserine-terizidone. Over the next 13 years, a dynamic evolution with coexisting, heterogeneous subpopulations was observed in 6 out of 13 sequential bacterial isolates. The emergence of drug-resistant subpopulations coincided with frequent changes in treatment regimens, which often included two or fewer active compounds. This evolutionary arms race between competing subpopulations ultimately resulted in the fixation of a single XDR variant. Our data demonstrate the complex intrapatient microevolution of Mtbc subpopulations during failing MDR/XDR-TB treatment. Designing effective treatment regimens based on rapid detection of (hetero) resistance is key to avoid resistance development and treatment failure.

Common Adaptive Strategies Underlie Within-Host Evolution of Bacterial Pathogens

Within-host adaptation is a hallmark of chronic bacterial infections, involving substantial genomic changes. Recent large-scale genomic data from prolonged infections allow the examination of adaptive strategies employed by different pathogens and open the door to investigate whether they converge toward similar strategies. Here, we compiled extensive data of whole-genome sequences of bacterial isolates belonging to miscellaneous species sampled at sequential time points during clinical infections. Analysis of these data revealed that different species share some common adaptive strategies, achieved by mutating various genes. Although the same genes were often mutated in several strains within a species, different genes related to the same pathway, structure, or function were changed in other species utilizing the same adaptive strategy (e.g., mutating flagellar genes). Strategies exploited by various bacterial species were often predicted to be driven by the host immune system, a powerful selective pressure that is not species specific. Remarkably, we find adaptive strategies identified previously within single species to be ubiquitous. Two striking examples are shifts from siderophore-based to heme-based iron scavenging (previously shown for Pseudomonas aeruginosa) and changes in glycerol-phosphate metabolism (previously shown to decrease sensitivity to antibiotics in Mycobacterium tuberculosis). Virulence factors were often adaptively affected in different species, indicating shifts from acute to chronic virulence and virulence attenuation during infection. Our study presents a global view on common within-host adaptive strategies employed by different bacterial species and provides a rich resource for further studying these processes.

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.

Clinical Interpretation of Drug Susceptibility Tests in Tuberculosis

:

Prompt and accurate diagnosis of drug resistance is essential for optimal treatment of drug-resistant tuberculosis. However, only 20% of the more than half a million patients eligible for the treatment of MDR-TB/RR-TB receive an appropriate drug regimen. Drug-resistant TB regimens must include a sufficient number of effective medications, a significant challenge for clinicians worldwide, as most are forced to make therapeutic decisions without any, or very little information on drug susceptibility testing. Although phenotypic DST is still commonly regarded as the gold standard for determining M. tuberculosis susceptibility to antituberculosis drugs, it has several limitations, mainly its prolonged turnaround time. Molecular methods based on M. tuberculosis genomic DNA sequencing have been developed during the past two decades, to identify the most common mutations involved in drug resistance. The Xpert ® MTB/RIF is a real-time polymerase chain reaction that offers results in less than two hours and has an overall sensitivity for rifampin resistance of 96% and 98% specificity. Line probe assays (LPAs) are commercial DNA strip-based tests for detecting the most frequent mutations responsible for resistance to rifampin, isoniazid, fluoroquinolones, and second-line injectable drugs.

:

Discrepancies between phenotypic and genotyping methods are a problem that the clinician will face in everyday practice. However, any resistance result (with any type of test) in a person with risk factors for harboring resistant microorganisms should be considered appropriate while the results of complementary tests are available.

Lesion Heterogeneity Coincides With Long-Term Heteroresistance in MDR-TB

Tuberculosis (TB) heteroresistance, in which only a fraction of the bacteria in a TB patient contains drug-resistant mutations, has been a rising concern. However, its origins and prevalence remain elusive. Here, whole-genome sequencing was performed on 83 serial isolates from 31 MDR-TB patients and heteroresistance was detected in isolates from 21 (67.74%) patients. Heteroresistance persisted in the host for long periods, spanning months to years, and was associated with having multiple tubercular lesions. Our findings indicate that heteroresistance is common and persistent in MDR-TB patients and may affect the success of their treatment regimens.

Sources of Multidrug Resistance in Patients With Previous Isoniazid-Resistant Tuberculosis Identified Using Whole Genome Sequencing: A Longitudinal Cohort Study

BACKGROUND

Meta-analysis of patients with isoniazid-resistant tuberculosis (TB) given standard first-line anti-TB treatment indicated an increased risk of multidrug-resistant TB (MDR-TB) emerging (8%), compared to drug-sensitive TB (0.3%). Here we use whole genome sequencing (WGS) to investigate whether treatment of patients with preexisting isoniazid-resistant disease with first-line anti-TB therapy risks selecting for rifampicin resistance, and hence MDR-TB.

METHODS

Patients with isoniazid-resistant pulmonary TB were recruited and followed up for 24 months. Drug susceptibility testing was performed by microscopic observation drug susceptibility assay, mycobacterial growth indicator tube, and by WGS on isolates at first presentation and in the case of re-presentation. Where MDR-TB was diagnosed, WGS was used to determine the genomic relatedness between initial and subsequent isolates. De novo emergence of MDR-TB was assumed where the genomic distance was 5 or fewer single-nucleotide polymorphisms (SNPs), whereas reinfection with a different MDR-TB strain was assumed where the distance was 10 or more SNPs.

RESULTS

Two hundred thirty-nine patients with isoniazid-resistant pulmonary TB were recruited. Fourteen (14/239 [5.9%]) patients were diagnosed with a second episode of TB that was multidrug resistant. Six (6/239 [2.5%]) were identified as having evolved MDR-TB de novo and 6 as having been reinfected with a different strain. In 2 cases, the genomic distance was between 5 and 10 SNPs and therefore indeterminate.

CONCLUSIONS

In isoniazid-resistant TB, de novo emergence and reinfection of MDR-TB strains equally contributed to MDR development. Early diagnosis and optimal treatment of isoniazid-resistant TB are urgently needed to avert the de novo emergence of MDR-TB during treatment.

Integron gene cassettes harboring novel variants of D-alanine-D-alanine ligase confer high-level resistance to D-cycloserine

Antibiotic resistance poses an increasing threat to global health. To tackle this problem, the identification of principal reservoirs of antibiotic resistance genes (ARGs) plus an understanding of drivers for their evolutionary selection are important. During a PCR-based screen of ARGs associated with integrons in saliva-derived metagenomic DNA of healthy human volunteers, two novel variants of genes encoding a d-alanine-d-alanine ligase (ddl6 and ddl7) located within gene cassettes in the first position of a reverse integron were identified. Treponema denticola was identified as the likely host of the ddl cassettes. Both ddl6 and ddl7 conferred high level resistance to d-cycloserine when expressed in Escherichia coli with ddl7 conferring four-fold higher resistance to D-cycloserine compared to ddl6. A SNP was found to be responsible for this difference in resistance phenotype between the two ddl variants. Molecular dynamics simulations were used to explain the mechanism of this phenotypic change at the atomic scale. A hypothesis for the evolutionary selection of ddl containing integron gene cassettes is proposed, based on molecular docking of plant metabolites within the ATP and d-cycloserine binding pockets of Ddl.

Dynamics of Extensive Drug Resistance Evolution of Mycobacterium tuberculosis in a Single Patient During 9 Years of Disease and Treatment

Mycobacterium tuberculosis is one of the hardest to treat bacterial pathogens with a high capacity to develop antibiotic resistance by mutations. Here we have performed whole-genome sequencing of consecutive M. tuberculosis isolates obtained during 9 years from a patient with pulmonary tuberculosis. The infecting strain was isoniazid resistant and during treatment it stepwise accumulated resistance mutations to 8 additional antibiotics. Heteroresistance was common and subpopulations with up to 3 different resistance mutations to the same drug coexisted. Sweeps of different resistant clones dominated the population at different time points, always coupled to resistance mutations coinciding with changes in the treatment regimens. Resistance mutations were predominant and no hitch-hiking, compensatory, or virulence-increasing mutations were detected, showing that the dominant selection pressure was antibiotic treatment. The results highlight the dynamic nature of M. tuberculosis infection, population structure, and resistance evolution and the importance of rapid antibiotic susceptibility tests to battle this pathogen.

Advances in the diagnosis of tuberculosis- Journey from smear microscopy to whole genome sequencing

The laboratory plays an important role in diagnosing tuberculosis (TB) and the identification and drug sensitivity testing (DST) of Mycobacterium tuberculosis. With a timely diagnosis and treatment with appropriate anti-TB drugs, most people who develop TB can be cured and onward transmission of infection curtailed. For a long time, laboratories used only microscopy and conventional culture-based diagnosis, however these procedures are slow and may require 3-4 weeks to yield results. Given the increasing rate of drug resistance, it has been necessary to look for new and rapid diagnostic methods. Various molecular based diagnostic technologies became available in the beginning of early 90s, providing rapid detection, identification and DST of M. tuberculosis. Molecular technologies offer the greatest potential for laboratories because they have the highest sensitivity and specificity. The present article will review some of the new methodology that has been introduced in the clinical laboratory.

Metabolic Changes of Mycobacterium tuberculosis during the Anti-Tuberculosis Therapy

Tuberculosis, caused by Mycobacterium tuberculosis complex bacteria, remains one of the most pressing health problems. Despite the general trend towards reduction of the disease incidence rate, the situation remains extremely tense due to the distribution of the resistant forms. Most often, these strains emerge through the intra-host microevolution of the pathogen during treatment failure. In the present study, the focus was on three serial clinical isolates of Mycobacterium tuberculosis Beijing B0/W148 cluster from one patient with pulmonary tuberculosis, to evaluate their changes in metabolism during anti-tuberculosis therapy. Using whole genome sequencing (WGS), 9 polymorphisms were determined, which occurred in a stepwise or transient manner during treatment and were linked to the resistance (GyrA D94A; inhA t-8a) or virulence. The effect of the inhA t-8a mutation was confirmed on both proteomic and transcriptomic levels. Additionally, the amount of RpsL protein, which is a target of anti-tuberculosis drugs, was reduced. At the systemic level, profound changes in metabolism, linked to the evolution of the pathogen in the host and the effects of therapy, were documented. An overabundance of the FAS-II system proteins (HtdX, HtdY) and expression changes in the virulence factors have been observed at the RNA and protein levels.

Phenotypic and Genotypic Characterization of Streptomycin-Resistant Multidrug-Resistant Mycobacterium tuberculosis Clinical Isolates in Southern China

Streptomycin (STR) is the first antibiotic used in the treatment of tuberculosis (TB) and the earliest antituberculosis drug with acquired resistance developed by Mycobacterium tuberculosis. The high prevalence of such resistance in many parts of the world limits its use for treating multidrug-resistant (MDR) TB. The aims of this study are to characterize of mutations in rpsL, rrs, and gidB genes in MDR M. tuberculosis isolates originating from southern China and to investigate possible relationship between mutations and strain genotypes for precise diagnosis and treatment. Sequences of rpsL, rrs, and gidB genes and the resistance profiles were analyzed for 218 MDR M. tuberculosis isolates. Our study showed that 68.35% of MDR M. tuberculosis isolates were resistant to STR and 89.91% of STR-resistant (STR^R) isolates were Beijing lineage strains. Mutations were observed in STR^R MDR M. tuberculosis isolates at the following rates: 72.48% in rpsL, 36.91% in rrs, and 15.44% in gidB. Compared with the phenotypic data, the combination of mutations in rpsL, rrs, and gidB has sensitivity and specificity of 96.64% and 100.00%, respectively. The most common mutations in STR^R isolates were rpsL128,263 and rrs514,1401, of which rpsL128 showed association with Beijing lineage (p < 0.001). It is noteworthy that a1401g mutation was present in rrs, while MDR M. tuberculosis isolates were resistant to both STR and amikacin. Twenty two novel mutations were found in STR^R isolates. These findings could be helpful to develop rapid molecular diagnostic methods and understand STR resistance in China for developing TB precision medicine and disturbance of drug-resistant TB transmission.

Comparative fitness analysis of D-cycloserine resistant mutants reveals both fitness-neutral and high-fitness cost genotypes

Drug resistant infections represent one of the most challenging medical problems of our time. D-cycloserine is an antibiotic used for six decades without significant appearance and dissemination of antibiotic resistant strains, making it an ideal model compound to understand what drives resistance evasion. We therefore investigated why Mycobacterium tuberculosis fails to become resistant to D-cycloserine. To address this question, we employed a combination of bacterial genetics, genomics, biochemistry and fitness analysis in vitro, in macrophages and in mice. Altogether, our results suggest that the ultra-low rate of emergence of D-cycloserine resistance mutations is the dominant biological factor delaying the appearance of clinical resistance to this antibiotic. Furthermore, we also identified potential compensatory mechanisms able to minimize the severe fitness costs of primary D-cycloserine resistance conferring mutations.

D-cycloserine (DCS) has been used for decades to treat Mycobacterium tuberculosis (Mtb) but resistance is rarely observed in clinical isolates. Here, the authors report ultra-low rate of emergence of resistance mutations as the underlying mechanism of DCS resistance evasion in Mtb.

Whole-genome sequencing for rapid, reliable and routine investigation of Mycobacterium tuberculosis transmission in local communities

Contact investigations following the diagnosis of active tuberculosis (TB) are paramount for the control of the disease. Epidemiological data are very powerful for contact tracing but might be delayed and/or difficult to integrate, especially in the setting of multiple contact-tracing investigations. The aim of this study was to address the added-value of whole-genome sequencing (WGS) to routine local TB surveillance systems. From November 2016 to July 2017, the local TB surveillance system identified three clusters that could constitute a unique larger outbreak. Epidemiological and clinical information were integrated with WGS genotyping data of Mycobacterium tuberculosis strains obtained using a simple DNA extraction method coupled with sequencing using an Illumina MiSeq platform and an in-house bioinformatics pipeline for single nucleotide polymorphism (SNP) analysis. Epidemiological investigations identified three putative TB clusters potentially interrelated including eight patients with active TB. Seven M. tuberculosis isolates were available and analysed by WGS. Using a 5-SNP threshold to define recent transmission, WGS-based genotyping supported the occurrence of the three clusters as well as a link between clusters 1 and 2 (SNP ≤1), constituting a larger outbreak. This outbreak was clearly delineated by refuting a potential link with the third cluster (SNP >500). Genotyping data did not support the belonging of patient 7 to any studied cluster. This study illustrates the usefulness of WGS genotyping for routine TB surveillance systems in local communities to rapidly confirm or disprove epidemiological hypotheses and delineate TB clusters, especially in the context of multiple contact-tracing investigations.

Cross-Border Movement of Highly Drug-Resistant Mycobacterium tuberculosis from Papua New Guinea to Australia through Torres Strait Protected Zone, 2010-2015

In this retrospective study, we used whole-genome sequencing (WGS) to delineate transmission dynamics, characterize drug-resistance markers, and identify risk factors of transmission among Papua New Guinea residents of the Torres Strait Protected Zone (TSPZ) who had tuberculosis diagnoses during 2010-2015. Of 117 isolates collected, we could acquire WGS data for 100; 79 were Beijing sublineage 2.2.1.1, which was associated with active transmission (odds ratio 6.190, 95% CI 2.221-18.077). Strains were distributed widely throughout the TSPZ. Clustering occurred more often within than between villages (p = 0.0013). Including 4 multidrug-resistant tuberculosis isolates from Australia citizens epidemiologically linked to the TSPZ into the transmission network analysis revealed 2 probable cross-border transmission events. All multidrug-resistant isolates (33/104) belonged to Beijing sublineage 2.2.1.1 and had high-level isoniazid and ethionamide co-resistance; 2 isolates were extensively drug resistant. Including WGS in regional surveillance could improve tuberculosis transmission tracking and control strategies within the TSPZ.

Direct transmission of within-host Mycobacterium tuberculosis diversity to secondary cases can lead to variable between-host heterogeneity without de novo mutation: A genomic investigation

BACKGROUND

Whole genome sequencing (WGS) has enabled the development of new approaches to track Mycobacterium tuberculosis (Mtb) transmission between tuberculosis (TB) cases but its utility may be challenged by the discovery that Mtb diversifies within hosts. Nevertheless, there is limited data on the presence and degree of within-host evolution.

METHODS

We profiled a well-documented Mtb transmission cluster with three pulmonary TB cases to investigate within-host evolution and describe its impact on recent transmission estimates. We used deep sequencing to track minority allele frequencies (<50·0% abundance) during transmission and standard treatment.

FINDINGS

Pre-treatment (n = 3) and serial samples collected over 2 months of antibiotic treatment (n = 16) from all three cases were analysed. Consistent with the epidemiological data, zero fixed SNP separated all genomes. However, we identified six subclones between the three cases with an allele frequency ranging from 35·0% to 100·0% across sampling intervals. Five subclones were identified within the index case pre-treatment and shared with one secondary case, while only the dominant clone was observed in the other secondary case. By tracking the frequency of these heterogeneous alleles over the two-month therapy, we observed distinct signatures of drift and negative selection, but limited evidence for de novo mutations, even under drug pressure.

INTERPRETATION

We document within-host Mtb diversity in an index case, which led to transmission of minority alleles to a secondary case. Incorporating data on heterogeneous alleles may refine our understanding of Mtb transmission dynamics. However, more evidence is needed on the role of transmission bottleneck on observed heterogeneity between cases.

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 Mycobacterium tuberculosis directly from sputum identifies more genetic diversity than sequencing from culture

Background

Repeated culture reduces within-sample Mycobacterium tuberculosis genetic diversity due to selection of clones suited to growth in culture and/or random loss of lineages, but it is not known to what extent omitting the culture step altogether alters genetic diversity. We compared M. tuberculosis whole genome sequences generated from 33 paired clinical samples using two methods. In one method DNA was extracted directly from sputum then enriched with custom-designed SureSelect (Agilent) oligonucleotide baits and in the other it was extracted from mycobacterial growth indicator tube (MGIT) culture.

Results

DNA directly sequenced from sputum showed significantly more within-sample diversity than that from MGIT culture (median 5.0 vs 4.5 heterozygous alleles per sample, p = 0.04). Resistance associated variants present as HAs occurred in four patients, and in two cases may provide a genotypic explanation for phenotypic resistance.

Conclusions

Culture-free M. tuberculosis whole genome sequencing detects more within-sample diversity than a leading culture-based method and may allow detection of mycobacteria that are not actively replicating.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5782-2) contains supplementary material, which is available to authorized users.

Bacterial Genome-Wide Association Identifies Novel Factors That Contribute to Ethionamide and Prothionamide Susceptibility in Mycobacterium tuberculosis

In Mycobacterium tuberculosis, recent genome-wide association studies have identified a novel constellation of mutations that are correlated with high-level drug resistances. Interpreting the functional importance of the new resistance-associated mutations has been complicated, however, by a lack of experimental validation and a poor understanding of the epistatic factors influencing these correlations, including strain background and programmatic variation in treatment regimens. Here we perform a genome-wide association analysis in a panel of Mycobacterium tuberculosis strains from China to identify variants correlated with resistance to the second-line prodrug ethionamide (ETH). Mutations in a bacterial monooxygenase, Rv0565c, are significantly associated with ETH resistance. We demonstrate that Rv0565c is a novel activator of ETH, independent of the two known activators, EthA and MymA. Clinically prevalent mutations abrogate Rv0565c function, and deletion of Rv0565c confers a consistent fitness benefit on M. tuberculosis in the presence of partially inhibitory doses of ETH. Interestingly, Rv0565c activity affects susceptibility to prothionamide (PTH), the ETH analog used in China, to a greater degree. Further, clinical isolates vary in their susceptibility to both ETH and PTH, to an extent that correlates with the total expression of ETH/PTH activators (EthA, MymA, and Rv0565c). These results suggest that clinical strains considered susceptible to ETH/PTH are not equally fit during treatment due to both Rv0565c mutations and more global variation in the expression of the prodrug activators.IMPORTANCE Phenotypic antibiotic susceptibility testing in Mycobacterium tuberculosis is slow and cumbersome. Rapid molecular diagnostics promise to help guide therapy, but such assays rely on complete knowledge of the molecular determinants of altered antibiotic susceptibility. Recent genomic studies of antibiotic-resistant M. tuberculosis have identified several candidate loci beyond those already known to contribute to antibiotic resistance; however, efforts to provide experimental validation have lagged. Our study identifies a gene (Rv0565c) that is associated with resistance to the second-line antibiotic ethionamide at a population level. We then use bacterial genetics to show that the variants found in clinical strains of M. tuberculosis improve bacterial survival after ethionamide exposure.

Deciphering Within-Host Microevolution of Mycobacterium tuberculosis through Whole-Genome Sequencing: the Phenotypic Impact and Way Forward

The Mycobacterium tuberculosis genome is more heterogenous and less genetically stable within the host than previously thought. Currently, only limited data exist on the within-host microevolution, diversity, and genetic stability of M. tuberculosis As a direct consequence, our ability to infer M. tuberculosis transmission chains and to understand the full complexity of drug resistance profiles in individual patients is limited. Furthermore, apart from the acquisition of certain drug resistance-conferring mutations, our knowledge on the function of genetic variants that emerge within a host and their phenotypic impact remains scarce. We performed a systematic literature review of whole-genome sequencing studies of serial and parallel isolates to summarize the knowledge on genetic diversity and within-host microevolution of M. tuberculosis We identified genomic loci of within-host emerged variants found across multiple studies and determined their functional relevance. We discuss important remaining knowledge gaps and finally make suggestions on the way forward.

Evaluation of a gene-by-gene approach for prospective whole-genome sequencing-based surveillance of multidrug resistant Mycobacterium tuberculosis

Whole-genome sequencing (WGS) offers unprecedented resolution for tracking Mycobacterium tuberculosis transmission and antibiotic-resistance spread. Still, the establishment of standardized WGS-based pipelines and the definition of epidemiological clusters based on genetic relatedness are under discussion. We aimed to implement a dynamic gene-by-gene approach, fully relying on freely available software, for prospective WGS-based tuberculosis surveillance, demonstrating its application for detecting transmission chains by retrospectively analysing all M/XDR strains isolated in 2013-2017 in Portugal. We observed a good correlation between genetic relatedness and epidemiological links, with strongly epilinked clusters displaying mean pairwise allele differences (AD) always below 0.3% (ratio of mean AD over the total number of shared loci between same-cluster strains). This data parallels the genetic distances acquired by the core-SNV analysis, while providing higher resolution and epidemiological concordance than MIRU-VNTR genotyping. The dynamic analysis of strain sub-sets (i.e., increasing the number of shared loci within each sub-set) also strengthens the confidence in detecting epilinked clusters. This gene-by-gene strategy also offers several practical benefits (e.g., reliance on freely-available software, scalability and low computational requirements) that further consolidated its suitability for a timely and robust prospective WGS-based laboratory surveillance of M/XDR-TB cases.

Drug susceptibility testing and mortality in patients treated for tuberculosis in high-burden countries: a multicentre cohort study

BACKGROUND

Drug resistance is a challenge for the global control of tuberculosis. We examined mortality in patients with tuberculosis from high-burden countries, according to concordance or discordance of results from drug susceptibility testing done locally and in a reference laboratory.

METHODS

This multicentre cohort study was done in Côte d'Ivoire, Democratic Republic of the Congo, Kenya, Nigeria, South Africa, Peru, and Thailand. We collected Mycobacterium tuberculosis isolates and clinical data from adult patients aged 16 years or older. Patients were stratified by HIV status and tuberculosis drug resistance. Molecular or phenotypic drug susceptibility testing was done locally and at the Swiss National Center for Mycobacteria, Zurich, Switzerland. We examined mortality during treatment according to drug susceptibility test results and treatment adequacy in multivariable logistic regression models adjusting for sex, age, sputum microscopy, and HIV status.

FINDINGS

We obtained M tuberculosis isolates from 871 patients diagnosed between 2013 and 2016. After exclusion of 237 patients, 634 patients with tuberculosis were included in this analysis; the median age was 33·2 years (IQR 26·9-42·5), 239 (38%) were women, 272 (43%) were HIV-positive, and 69 (11%) patients died. Based on the reference laboratory drug susceptibility test, 394 (62%) strains were pan-susceptible, 45 (7%) monoresistant, 163 (26%) multidrug-resistant (MDR), and 30 (5%) had pre-extensively or extensively drug resistant (pre-XDR or XDR) tuberculosis. Results of reference and local laboratories were concordant for 513 (81%) of 634 patients and discordant for 121 (19%) of 634. Overall, sensitivity to detect any resistance was 90·8% (95% CI 86·5-94·2) and specificity 84·3% (80·3-87·7). Mortality ranged from 6% (20 of 336) in patients with pan-susceptible tuberculosis treated according to WHO guidelines to 57% (eight of 14) in patients with resistant strains who were under-treated. In logistic regression models, compared with concordant drug susceptibility test results, the adjusted odds ratio of death was 7·33 (95% CI 2·70-19·95) for patients with discordant results potentially leading to under-treatment.

INTERPRETATION

Inaccurate drug susceptibility testing by comparison with a reference standard leads to under-treatment of drug-resistant tuberculosis and increased mortality. Rapid molecular drug susceptibility test of first-line and second-line drugs at diagnosis is required to improve outcomes in patients with MDR tuberculosis and pre-XDR or XDR tuberculosis.

FUNDING

National Institutes of Allergy and Infectious Diseases, Swiss National Science Foundation, Swiss National Center for Mycobacteria.

Mycobacterial genomics and structural bioinformatics: opportunities and challenges in drug discovery

Of the more than 190 distinct species of Mycobacterium genus, many are economically and clinically important pathogens of humans or animals. Among those mycobacteria that infect humans, three species namely Mycobacterium tuberculosis (causative agent of tuberculosis), Mycobacterium leprae (causative agent of leprosy) and Mycobacterium abscessus (causative agent of chronic pulmonary infections) pose concern to global public health. Although antibiotics have been successfully developed to combat each of these, the emergence of drug-resistant strains is an increasing challenge for treatment and drug discovery. Here we describe the impact of the rapid expansion of genome sequencing and genome/pathway annotations that have greatly improved the progress of structure-guided drug discovery. We focus on the applications of comparative genomics, metabolomics, evolutionary bioinformatics and structural proteomics to identify potential drug targets. The opportunities and challenges for the design of drugs for M. tuberculosis, M. leprae and M. abscessus to combat resistance are discussed.

Insights into the processes that drive the evolution of drug resistance in Mycobacterium tuberculosis

At present, the successful transmission of drug-resistant Mycobacterium tuberculosis, including multidrug-resistant (MDR) and extensively drug-resistant (XDR) strains, in human populations, threatens tuberculosis control worldwide. Differently from many other bacteria, M. tuberculosis drug resistance is acquired mainly through mutations in specific drug resistance-associated genes. The panel of mutations is highly diverse, but depends on the affected gene and M. tuberculosis genetic background. The variety of genetic profiles observed in drug-resistant clinical isolates underlines different evolutionary trajectories towards multiple drug resistance, although some mutation patterns are prominent. This review discusses the intrinsic processes that may influence drug resistance evolution in M. tuberculosis, such as mutation rate, drug resistance-associated mutations, fitness cost, compensatory mutations and epistasis. This knowledge should help to better predict the risk of emergence of highly resistant M. tuberculosis strains and to develop new tools and strategies to limit the development and spread of MDR and XDR strains.

Insight into multidrug-resistant Beijing genotype Mycobacterium tuberculosis isolates in Myanmar

OBJECTIVES

Myanmar is a World Health Organization high tuberculosis (TB) burden country with a high multidrug-resistant (MDR)-TB burden. Of significance, a high prevalence of the Beijing genotype of Mycobacterium tuberculosis (MTB) among MDR-MTB has been reported previously. A detailed genetic characterization of TB clinical isolates was performed in order to explore whether there is an association between the prevalence of the Beijing MTB genotype and MDR-TB in Myanmar.

METHODS

A total of 265 MDR-MTB clinical isolates collected in 2010 and 2012 were subjected to spoligotyping, mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) analysis, single nucleotide polymorphism (SNP) typing, and drug resistance-associated gene sequencing, including rpoC to detect potential compensatory evolution.

RESULTS

Of the total MDR-MTB isolates, 79.2% (210/265) were of the Beijing genotype, the majority of which were the 'modern' subtype. Beijing genotype isolates were differentiated by 15-locus MIRU-VNTR and a high clustering rate (53.0%) was observed in the modern subtype. These MIRU-VNTR patterns were similar to Beijing genotype clones spreading across Russia and Central Asia. A high prevalence of katG Ser315Thr, and genetic evidence of extensive drug resistance (XDR) and pre-XDR and compensatory mutations in rpoC were observed among clustered isolates.

CONCLUSIONS

MDR-MTB strains of the Beijing genotype might be spreading in Myanmar and present a major challenge to TB control in this country.

Perspectives for personalized therapy for patients with multidrug-resistant tuberculosis

According to the World Health Organization (WHO), tuberculosis is the leading cause of death attributed to a single microbial pathogen worldwide. In addition to the large number of patients affected by tuberculosis, the emergence of Mycobacterium tuberculosis drug-resistance is complicating tuberculosis control in many high-burden countries. During the past 5 years, the global number of patients identified with multidrug-resistant tuberculosis (MDR-TB), defined as bacillary resistance at least against rifampicin and isoniazid, the two most active drugs in a treatment regimen, has increased by more than 20% annually. Today we experience a historical peak in the number of patients affected by MDR-TB. The management of MDR-TB is characterized by delayed diagnosis, uncertainty of the extent of bacillary drug-resistance, imprecise standardized drug regimens and dosages, very long duration of therapy and high frequency of adverse events which all translate into a poor prognosis for many of the affected patients. Major scientific and technological advances in recent years provide new perspectives through treatment regimens tailor-made to individual needs. Where available, such personalized treatment has major implications on the treatment outcomes of patients with MDR-TB. The challenge now is to bring these adances to those patients that need them most.

Unique Features and Anti-microbial Targeting of Folate- and Flavin-Dependent Methyltransferases Required for Accurate Maintenance of Genetic Information

Comparative genome analyses have led to the discovery and characterization of novel flavin- and folate-dependent methyltransferases that mainly function in DNA precursor synthesis and post-transcriptional RNA modification by forming (ribo) thymidylate and its derivatives. Here we discuss the recent literature on the novel mechanistic features of these enzymes sometimes referred to as “uracil methyltransferases,” albeit we prefer to refer to them as (ribo) thymidylate synthases. These enzyme families attest to the convergent evolution of nucleic acid methylation. Special focus is given to describing the unique characteristics of these flavin- and folate-dependent enzymes that have emerged as new models for studying the non-canonical roles of reduced flavin co-factors (FADH₂) in relaying carbon atoms between enzyme substrates. This ancient enzymatic methylation mechanism with a very wide phylogenetic distribution may be more commonly used for biological methylation reactions than previously anticipated. This notion is exemplified by the recent discovery of additional substrates for these enzymes. Moreover, similar reaction mechanisms can be reversed by demethylases, which remove methyl groups e.g., from human histones. Future work is now required to address whether the use of different methyl donors facilitates the regulation of distinct methylation reactions in the cell. It will also be of great interest to address whether the low activity flavin-dependent thymidylate synthases ThyX represent ancestral enzymes that were eventually replaced by the more active thymidylate synthases of the ThyA family to facilitate the maintenance of larger genomes in fast-growing microbes. Moreover, we discuss the recent efforts from several laboratories to identify selective anti-microbial compounds that target flavin-dependent thymidylate synthase ThyX. Altogether we underline how the discovery of the alternative flavoproteins required for methylation of DNA and/or RNA nucleotides, in addition to providing novel targets for antibiotics, has provided new insight into microbial physiology and virulence.

Comparative genomic analysis of multidrug-resistant Streptococcus pneumoniae isolates

Introduction

Multidrug resistance in Streptococcus pneumoniae has emerged as a serious problem to public health. A further understanding of the genetic diversity in antibiotic-resistant S. pneumoniae isolates is needed.

Methods

We conducted whole-genome resequencing for 25 pneumococcal strains isolated from children with different antimicrobial resistance profiles. Comparative analysis focus on detection of single-nucleotide polymorphisms (SNPs) and insertions and deletions (indels) was conducted. Moreover, phylogenetic analysis was applied to investigate the genetic relationship among these strains.

Results

The genome size of the isolates was ~2.1 Mbp, covering >90% of the total estimated size of the reference genome. The overall G+C% content was ~39.5%, and there were 2,200–2,400 open reading frames. All isolates with different drug resistance profiles harbored many indels (range 131–171) and SNPs (range 16,103–28,128). Genetic diversity analysis showed that the variation of different genes were associated with specific antibiotic resistance. Known antibiotic resistance genes (pbps, murMN, ciaH, rplD, sulA, and dpr) were identified, and new genes (regR, argH, trkH, and PTS-EII) closely related with antibiotic resistance were found, although these genes were primarily annotated with functions in virulence as well as carbohydrate and amino acid transport and metabolism. Phylogenetic analysis unambiguously indicated that isolates with different antibiotic resistance profiles harbored similar genetic backgrounds. One isolate, 14-LC.ER1025, showed a much weaker phylogenetic relationship with the other isolates, possibly caused by genomic variation.

Conclusion

In this study, although pneumococcal isolates had similar genetic backgrounds, strains were diverse at the genomic level. These strains exhibited distinct variations in their indel and SNP compositions associated with drug resistance.

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.

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

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

Mycobacterium tuberculosis Acquires Limited Genetic Diversity in Prolonged Infections, Reactivations and Transmissions Involving Multiple Hosts

Background:Mycobacterium tuberculosis (MTB) has limited ability to acquire variability. Analysis of its microevolution might help us to evaluate the pathways followed to acquire greater infective success. Whole-genome sequencing (WGS) in the analysis of the transmission of MTB has elucidated the magnitude of variability in MTB. Analysis of transmission currently depends on the identification of clusters, according to the threshold of variability (<5 SNPs) between isolates. Objective: We evaluated whether the acquisition of variability in MTB, was more frequent in situations which could favor it, namely intrapatient, prolonged infections or reactivations and interpatient transmissions involving multiple sequential hosts. Methods: We used WGS to analyze the accumulation of variability in sequential isolates from prolonged infections or translations from latency to reactivation. We then measured microevolution in transmission clusters with prolonged transmission time, high number of involved cases, simultaneous involvement of latency and active transmission. Results: Intrapatient and interpatient acquisition of variability was limited, within the ranges expected according to the thresholds of variability proposed, even though bursts of variability were observed. Conclusions: The thresholds of variability proposed for MTB seem to be valid in most circumstances, including those theoretically favoring acquisition of variability. Our data point to multifactorial modulation of microevolution, although further studies are necessary to elucidate the factors underlying this modulation.

Multi-clonal evolution of multi-drug-resistant/extensively drug-resistant Mycobacterium tuberculosis in a high-prevalence setting of Papua New Guinea for over three decades

An outbreak of multi-drug resistant (MDR) tuberculosis (TB) has been reported on Daru Island, Papua New Guinea. Mycobacterium tuberculosis strains driving this outbreak and the temporal accrual of drug resistance mutations have not been described. Whole genome sequencing of 100 of 165 clinical isolates referred from Daru General Hospital to the Supranational reference laboratory, Brisbane, during 2012–2015 revealed that 95 belonged to a single modern Beijing sub-lineage strain. Molecular dating suggested acquisition of streptomycin and isoniazid resistance in the 1960s, with potentially enhanced virulence mediated by an mycP1 mutation. The Beijing sub-lineage strain demonstrated a high degree of co-resistance between isoniazid and ethionamide (80/95; 84.2 %) attributed to an inhA promoter mutation combined with inhA and ndh coding mutations. Multi-drug resistance, observed in 78/95 samples, emerged with the acquisition of a typical rpoB mutation together with a compensatory rpoC mutation in the 1980s. There was independent acquisition of fluoroquinolone and aminoglycoside resistance, and evidence of local transmission of extensively drug resistant (XDR) strains from 2009. These findings underline the importance of whole genome sequencing in informing an effective public health response to MDR/XDR TB.

Role of Alanine Racemase Mutations in Mycobacterium tuberculosis d-Cycloserine Resistance

A screening of more than 1,500 drug-resistant strains of Mycobacterium tuberculosis revealed evolutionary patterns characteristic of positive selection for three alanine racemase (Alr) mutations. We investigated these mutations using molecular modeling, in vitro MIC testing, as well as direct measurements of enzymatic activity, which demonstrated that these mutations likely confer resistance to d-cycloserine.

Diversity and evolution of drug resistance mechanisms in Mycobacterium tuberculosis

Despite the efficacy of antibiotics to protect humankind against many deadly pathogens, such as Mycobacterium tuberculosis, nothing can prevent the emergence of drug-resistant strains. Several mechanisms facilitate drug resistance in M. tuberculosis including compensatory evolution, epistasis, clonal interference, cell wall integrity, efflux pumps, and target mimicry. In this study, we present recent findings relevant to these mechanisms, which can enable the discovery of new drug targets and subsequent development of novel drugs for treatment of drug-resistant M. tuberculosis.

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.

Early diagnosis and effective treatment regimens are the keys to tackle antimicrobial resistance in tuberculosis (TB): A report from Euroscicon's international TB Summit 2016

To say that tuberculosis (TB) has regained a strong foothold in the global human health and wellbeing scenario would be an understatement. Ranking alongside HIV/AIDS as the top reason for mortality due to a single infectious disease, the impact of TB extends far into socio-economic context worldwide. As global efforts led by experts and political bodies converge to mitigate the predicted outcome of growing antimicrobial resistance, the academic community of students, practitioners and researchers have mobilised to develop integrated, inter-disciplinary programmes to bring the plans of the former to fruition. Enabling this crucial requirement for unimpeded dissemination of scientific discovery was the TB Summit 2016, held in London, United Kingdom. This report critically discusses the recent breakthroughs made in diagnostics and treatment while bringing to light the major hurdles in the control of the disease as discussed in the course of the 3-day international event. Conferences and symposia such as these are the breeding grounds for successful local and global collaborations and therefore must be supported to expand the understanding and outreach of basic science research.