Whole Genome Sequencing Investigation of a Tuberculosis Outbreak in Port-au-Prince, Haiti Caused by a Strain with a "Low-Level" rpoB Mutation L511P - Insights into a Mechanism of Resistance Escalation

Phenotypic and genotypic analysis of drug resistance in M. tuberculosis isolates in Gansu, China

Tuberculosis has posed a serious threat to human health. It is imperative to investigate the geographic prevalence of tuberculosis and medication resistance, as this information is essential for informing strategies for its prevention and treatment. Drug resistance was identified using a proportion method. Drug-resistant genes and pathways were predicted using whole genome sequencing. The drug resistance range of bedaquiline was identified using the microporous plate two-fold dilution method, and drug resistance genes were studied using sequencing. The study revealed that 19.99% of the tuberculosis cases had multidrug resistance. The genes of M. tuberculosis are predominantly involved in the synthesis of ABC transporters, two-component systems, and bacterial secretion systems, as well as in energy production and conversion, and lipid transport and metabolism. The genes encode for 82.45% of carbohydrate-related enzymes such as glycoside hydrolases, glycosyl transferases, and carbohydrate esterases. The minimum inhibitory concentration (MIC) of bedaquiline against clinical strains was approximately 0.06 μg/mL, with identified mutations in drug-resistant genes Rv0678, atpE, and pepQ, specifically V152A, P62A, and T222N, respectively. The multidrug resistance tuberculosis development was attributed to the strong medication resistance exhibited. It was concluded that tuberculosis had presented a high level of drug resistance. Phenotypic resistance was related to genes, existing potential genetic resistance in M. tuberculosis. Bedaquiline was found to possess effective antibacterial properties against M. tuberculosis.

Rapid detection of multidrug resistance in tuberculosis using nanopore-based targeted next-generation sequencing: a multicenter, double-blind study

Background

Resistance to anti-tuberculous drugs is a major challenge in the treatment of tuberculosis (TB). We aimed to evaluate the clinical availability of nanopore-based targeted next-generation sequencing (NanoTNGS) for the diagnosis of drug-resistant tuberculosis (DR-TB).

Methods

This study enrolled 253 patients with suspected DR-TB from six hospitals. The diagnostic efficacy of NanoTNGS for detecting Mycobacterium tuberculosis and its susceptibility or resistance to first- and second-line anti-tuberculosis drugs was assessed by comparing conventional phenotypic drug susceptibility testing (pDST) and Xpert MTB/RIF assays. NanoTNGS can be performed within 12 hours from DNA extraction to the result delivery.

Results

NanoTNGS showed a remarkable concordance rate of 99.44% (179/180) with the culture assay for identifying the Mycobacterium tuberculosis complex. The sensitivity of NanoTNGS for detecting drug resistance was 93.53% for rifampicin, 89.72% for isoniazid, 85.45% for ethambutol, 74.00% for streptomycin, and 88.89% for fluoroquinolones. Specificities ranged from 83.33% to 100% for all drugs tested. Sensitivity for rifampicin-resistant tuberculosis using NanoTNGS increased by 9.73% compared to Xpert MTB/RIF. The most common mutations were S531L (codon in E. coli) in the rpoB gene, S315T in the katG gene, and M306V in the embB gene, conferring resistance to rifampicin, isoniazid, and ethambutol, respectively. In addition, mutations in the pncA gene, potentially contributing to pyrazinamide resistance, were detected in 32 patients. Other prevalent variants, including D94G in the gyrA gene and K43R in the rpsL gene, conferred resistance to fluoroquinolones and streptomycin, respectively. Furthermore, the rv0678 R94Q mutation was detected in one sample, indicating potential resistance to bedaquiline.

Conclusion

NanoTNGS rapidly and accurately identifies resistance or susceptibility to anti-TB drugs, outperforming traditional methods. Clinical implementation of the technique can recognize DR-TB in time and provide guidance for choosing appropriate antituberculosis agents.

Tuberculosis Outbreak in an Educational Institution in Henan Province, China

On June 17, 2018, a case of pulmonary tuberculosis (TB) was reported among students at a senior high school in Luoning, China. The outbreak encompassed a total of 23 cases along with TB screening in the whole school by means of PPD and chest X-ray. By the end of September 2018, the entire 9 cases cultured positive had epidemiological association. All of the 9 Mycobacterium tuberculosis (Mtb) isolates available were sensitive to all drugs tested and had similar spoligotyping and 15 loci mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) profile. Whole-genome sequencing (WGS) of the Mtb isolates revealed 20 variable nucleotide positions within 8 cases, indicating a clonal outbreak. The index case, which was first identified and diagnosed, is separated from the cluster by a minimum number of 95 distinct SNPs. Minimum distance spanning tree (MST) indicted that the 8 cases were indeed part of a single transmission chain. It was concluded that this is an epidemic situation of TB outbreak exposed by the aggrieved index case at school, which was caused by the veiled infectious case wherein a student was suffering from TB and attending school simultaneously.

Rifampicin-Monoresistant Tuberculosis Is Not the Same as Multidrug-Resistant Tuberculosis: a Descriptive Study from Khayelitsha, South Africa

Rifampin monoresistance (RMR; rifampin resistance and isoniazid susceptibility) accounts for 38% of all rifampin-resistant tuberculosis (RR-TB) in South Africa and is increasing. We aimed to compare RMR-TB with multidrug-resistant TB (MDR-TB) in a setting with high TB, RR-TB, and HIV burdens. Patient-level clinical data and stored RR Mycobacterium tuberculosis isolates from 2008 to 2017 with available whole-genome sequencing (WGS) data were used to describe risk factors associated with RMR-TB and to compare RR-conferring mutations between RMR-TB and MDR-TB. A subset of isolates with particular RR-conferring mutations were subjected to semiquantitative rifampin phenotypic drug susceptibility testing. Among 2,041 routinely diagnosed RR-TB patients, 463 (22.7%) had RMR-TB. HIV-positive individuals (adjusted odds ratio [aOR], 1.4; 95% confidence interval [CI], 1.1 to 1.9) and diagnosis between 2013 and 2017 versus between 2008 and 2012 (aOR, 1.3; 95% CI, 1.1 to 1.7) were associated with RMR-TB. Among 1,119 (54.8%) patients with available WGS data showing RR-TB, significant differences in the distribution of rpoB RR-conferring mutations between RMR and MDR isolates were observed. Mutations associated with high-level RR were more commonly found among MDR isolates (811/889 [90.2%] versus 162/230 [70.4%] among RMR isolates; P < 0.0001). In particular, the rpoB L430P mutation, conferring low-level RR, was identified in 32/230 (13.9%) RMR isolates versus 10/889 (1.1%) in MDR isolates (P < 0.0001). Among 10 isolates with an rpoB L430P mutation, 7 were phenotypically susceptible using the critical concentration of 0.5 μg/ml (range, 0.125 to 1 μg/ml). The majority (215/230 [93.5%]) of RMR isolates showed susceptibility to all other TB drugs, highlighting the potential benefits of WGS for simplified treatment. These data suggest that the evolution of RMR-TB differs from MDR-TB with a potential contribution from HIV infection.

Low-Level Rifampin Resistance and rpoB Mutations in Mycobacterium tuberculosis: an Analysis of Whole-Genome Sequencing and Drug Susceptibility Test Data in New York

Rapid and reliable detection of rifampin (RIF) resistance is critical for the diagnosis and treatment of drug-resistant and multidrug-resistant (MDR) tuberculosis. Discordant RIF phenotype/genotype susceptibility results remain a challenge due to the presence of rpoB mutations that do not confer high levels of RIF resistance, as have been exhibited in strains with mutations such as Ser450Leu. These strains, termed low-level RIF resistant, exhibit elevated RIF MICs compared to fully susceptible strains but remain phenotypically susceptible by mycobacterial growth indicator tube (MGIT) testing and have been associated with poor patient outcomes. Here, we assess RIF resistance prediction by whole-genome sequencing (WGS) among a set of 1,779 prospectively tested strains by both prevalence of rpoB gene mutation and phenotype as part of routine clinical testing during a 2.5-year period. During this time, 139 strains were found to have nonsynonymous rpoB mutations, 53 of which were associated with RIF resistance, including both low-level and high-level resistance. Resistance to RIF (1.0 μg/ml in MGIT) was identified in 43 (81.1%) isolates. The remaining 10 (18.9%) strains were susceptible by MGIT but were confirmed to be low-level RIF resistant by MIC testing. Full rpoB gene sequencing overcame the limitations of critical concentration phenotyping, probe-based genotyping, and partial gene sequencing methods. Universal clinical WGS with concurrent phenotypic testing provided a more complete understanding of the prevalence and type of rpoB mutations and their association with RIF resistance in New York.

Drug-resistant TB prevalence study in 5 health institutions in Haiti

Objectives

Tuberculosis (TB) is the leading infectious cause of death in the world. Multi-drug resistant TB (MDR-TB) is a major public health problem as treatment is long, costly, and associated to poor outcomes. Here, we report epidemiological data on the prevalence of drug-resistant TB in Haiti.

Methods

This cross-sectional prevalence study was conducted in five health centers across Haiti. Adult, microbiologically confirmed pulmonary TB patients were included. Molecular genotyping (rpoB gene sequencing and spoligotyping) and phenotypic drug susceptibility testing were used to characterize rifampin-resistant MTB isolates detected by Xpert MTB/RIF.

Results

Between April 2016 and February 2018, 2,777 patients were diagnosed with pulmonary TB by Xpert MTB/RIF screening and positive MTB cultures. A total of 74 (2.7%) patients were infected by a drug-resistant (DR-TB) M. tuberculosis strain. Overall HIV prevalence was 14.1%. Patients with HIV infection were at a significantly higher risk for infection with DR-TB strains compared to pan-susceptible strains (28.4% vs. 13.7%, adjusted odds ratio 2.6, 95% confidence interval 1.5–4.4, P = 0.001). Among the detected DR-TB strains, T1 (29.3%), LAM9 (13.3%), and H3 (10.7%) were the most frequent clades. In comparison with previous spoligotypes studies with data collected in 2000–2002 and in 2008–2009 on both sensitive and resistant strains of TB in Haiti, we observed a significant increase in the prevalence of the drug-resistant MTB Spoligo-International-Types (SIT) 137 (X2 clade: 8.1% vs. 0.3% in 2000–02 and 0.9% in 2008–09, p<0.001), 5 (T1 clade: 6.8% vs 1.9 in 2000–02 and 1.7% in 2008–09, P = 0.034) and 455 (T1 clade: 5.4% vs 1.6% and 1.1%, P = 0.029). Newly detected spoligotypes (SIT 6, 7, 373, 909 and 1624) were also recorded.

Conclusion

This study describes the genotypic and phenotypic characteristics of DR-TB strains circulating in Haiti from April 2016 to February 2018. Newly detected MTB clades harboring multi-drug resistance patterns among the Haitian population as well as the higher risk of MDR-TB infection in HIV-positive people highlights the epidemiological relevance of these surveillance data. The importance of detecting RIF-resistant patients, as proxy for MDR-TB in peripheral sites via molecular techniques, is particularly important to provide adequate patient case management, prevent the transmission of resistant strains in the community and to contribute to the surveillance of resistant strains.

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.

Assessing Nanopore Sequencing for Clinical Diagnostics: a Comparison of Next-Generation Sequencing (NGS) Methods for Mycobacterium tuberculosis

Next-generation sequencing technologies are being rapidly adopted as a tool of choice for diagnostic and outbreak investigation in public health laboratories. However, costs of operation and the need for specialized staff remain major hurdles for laboratories with limited resources for implementing these technologies. This project aimed to assess the feasibility of using Oxford Nanopore MinION whole-genome sequencing data of Mycobacterium tuberculosis isolates for species identification, in silico spoligotyping, detection of mutations associated with antimicrobial resistance (AMR) to accurately predict drug susceptibility profiles, and phylogenetic analysis to detect transmission between cases. The results were compared prospectively in real time to those obtained with our current clinically validated Illumina MiSeq sequencing assay for M. tuberculosis and phenotypic drug susceptibility testing results when available. Our assessment of 431 sequenced samples over a 32-week period demonstrates that, when using the proper quality controls and thresholds, the MinION can achieve levels of genotyping analysis and phenotypic resistance predictions comparable to those of the Illumina MiSeq at a very competitive cost per sample. Our results indicate that nanopore sequencing can be a suitable alternative to, or complement, currently used sequencing platforms in a clinical setting and has the potential to be widely adopted in public health laboratories in the near future.

Role of Disputed Mutations in the rpoB Gene in Interpretation of Automated Liquid MGIT Culture Results for Rifampin Susceptibility Testing of Mycobacterium tuberculosis

Low-level rifampin resistance associated with specific rpoB mutations (referred as "disputed") in Mycobacterium tuberculosis is easily missed by some phenotypic methods. To understand the mechanism by which some mutations are systematically missed by MGIT phenotypic testing, we performed an in silico analysis of their effect on the structural interaction between the RpoB protein and rifampin. We also characterized 24 representative clinical isolates by determining MICs on 7H10 agar and testing them by an extended MGIT protocol. We analyzed 2,097 line probe assays, and 156 (7.4%) cases showed a hybridization pattern referred to here as "no wild type + no mutation." Isolates harboring "disputed" mutations (L430P, D435Y, H445C/L/N/S, and L452P) tested susceptible in MGIT, with prevalence ranging from 15 to 57% (overall, 16 out of 55 isolates [29%]). Our in silico analysis did not highlight any difference between "disputed" and "undisputed" substitutions, indicating that all rpoB missense mutations affect the rifampin binding site. MIC testing showed that "undisputed" mutations are associated with higher MIC values (≥20 mg/liter) compared to "disputed" mutations (4 to >20 mg/liter). Whereas "undisputed" mutations didn't show any delay (Δ) in time to positivity of the test tube compared to the control tube on extended MGIT protocol, "disputed" mutations showed a mean Δ of 7.2 days (95% confidence interval [CI], 4.2 to 10.2 days; P < 0.05), providing evidence that mutations conferring low-level resistance are associated with a delay in growth on MGIT. Considering the proved relevance of L430P, D435Y, H445C/L/N, and L452P mutations in determining clinical resistance, genotypic drug susceptibility testing (DST) should be used to replace phenotypic results (MGIT) when such mutations are found.

Revisiting the mutant prevention concentration to guide dosing in childhood tuberculosis

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

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.

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.

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

The availability of whole-genome sequencing (WGS) as a tool for the diagnosis and clinical management of tuberculosis (TB) offers considerable promise in the fight against this stubborn epidemic. However, like other new technologies, the best application of WGS remains to be determined, for both conceptual and technical reasons. In this review, we consider the potential value of WGS in the clinical laboratory for the detection of Mycobacterium tuberculosis and the prediction of antibiotic resistance. We also discuss issues pertaining to data generation, interpretation and dissemination, given that WGS has to date been generally performed in research labs where results are not necessarily packaged in a clinician-friendly format. Although WGS is far more accessible now than it was in the past, the transition from a research tool to study TB into a clinical test to manage this disease may require further fine-tuning. Improvements will likely come through iterative efforts that involve both the laboratories ready to move TB into the genomic era and the front-line clinical/public health staff who will be interpreting the results to inform management decisions.