Comprehensive Whole-Genome Sequencing and Reporting of Drug Resistance Profiles on Clinical Cases of Mycobacterium tuberculosis in New York State

Evaluation of WHO catalog of mutations and five WGS analysis tools for drug resistance prediction of Mycobacterium tuberculosis isolates from China

The continuous advancement of molecular diagnostic techniques, particularly whole-genome sequencing (WGS), has greatly facilitated the early diagnosis of drug-resistant tuberculosis patients. Nonetheless, the interpretation of results from various types of mutations in drug-resistant-associated genes has become the primary challenge in the field of molecular drug-resistance diagnostics. In this study, our primary objective is to evaluate the diagnosis accuracy of the World Health Organization (WHO) catalog of mutations and five WGS analysis tools (PhyResSE, Mykrobe, TB Profiler, Gen-TB, and SAM-TB) in drug resistance to 10 anti-Mycobacterium tuberculosis (MTB) drugs. We utilized the data of WGS collected between 2014 and 2017 in Zhejiang Province, consisting of 110 MTB isolates as detailed in our previous study. Based on phenotypic drug susceptibility testing (DST) results using the proportion method on Löwenstein-Jensen medium with antibiotics, we evaluated the predictive accuracy of genotypic DST obtained by these tools. The results revealed that the WHO catalog of mutations and five WGS analysis tools exhibit robust predictive capabilities concerning resistance to isoniazid, rifampicin, ethambutol, streptomycin, amikacin, kanamycin, and capreomycin. Notably, Mykrobe, SAM-TB, and TB Profiler demonstrate the most accurate predictions for resistance to pyrazinamide, prothionamide, and para-aminosalicylic acid, respectively. These findings are poised to significantly guide and influence future clinical treatment strategies and resistance monitoring protocols.IMPORTANCEWhole-genome sequencing (WGS) has the potential for the early diagnosis of drug-resistant tuberculosis. However, the interpretation of mutations of drug-resistant-associated genes represents a significant challenge as the amount and complexity of WGS data. We evaluated the accuracy of the World Health Organization catalog of mutations and five WGS analysis tools in predicting drug resistance to first-line and second-line anti-TB drugs. Our results offer clinicians guidance on selecting appropriate WGS analysis tools for predicting resistance to specific anti-TB drugs.

New insight in molecular detection of Mycobacterium tuberculosis

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

Exploring the potential of Oxford Nanopore Technologies sequencing for Mycobacterium tuberculosis sequencing: An assessment of R10 flowcells and V14 chemistry

Oxford Nanopore Technologies (ONT) sequencing is a promising technology. We assessed the performance of the new ONT R10 flowcells and V14 rapid sequencing chemistry for Mtb whole genome sequencing of Mycobacterium tuberculosis (Mtb) DNA extracted from clinical primary liquid cultures (CPLCs). Using the recommended protocols for MinION Mk1C, R10.4.1 MinION flowcells, and the ONT Rapid Sequencing Kit V14 on six CPLC samples, we obtained a pooled library yield of 10.9 ng/μl, generated 1.94 Gb of sequenced bases and 214k reads after 48h in a first sequencing run. Only half (49%) of all generated reads met the Phred Quality score threshold (>8). To assess if the low data output and sequence quality were due to impurities present in DNA extracted directly from CPLCs, we added a pre-library preparation bead-clean-up step and included purified DNA obtained from an Mtb subculture as a control sample in a second sequencing run. The library yield for DNA extracted from four CPLCs and one Mtb subculture (control) was similar (10.0 ng/μl), 2.38 Gb of bases and 822k reads were produced. The quality was slightly better with 66% of the produced reads having a Phred Quality >8. A third run of DNA from six CPLCs with bead clean-up pre-processing produced a low library yield (±1 Gb of bases, 166k reads) of low quality (51% of reads with a Phred Quality score >8). A median depth of coverage above 10× was only achieved for five of 17 (29%) sequenced libraries. Compared to Illumina WGS of the same samples, accurate lineage predictions and full drug resistance profiles from the generated ONT data could not be determined by TBProfiler. Further optimization of the V14 ONT rapid sequencing chemistry and library preparation protocol is needed for clinical Mtb WGS applications.

Evaluation of Whole Genome Sequencing-Based Predictions of Antimicrobial Resistance to TB First Line Agents: A Lesson from 5 Years of Data

Phenotypic susceptibility testing of the Mycobacterium tuberculosis complex (MTBC) isolate requires culture growth, which can delay rapid detection of resistant cases. Whole genome sequencing (WGS) and data analysis pipelines can assist in predicting resistance to antimicrobials used in the treatment of tuberculosis (TB). This study compared phenotypic susceptibility testing results and WGS-based predictions of antimicrobial resistance (AMR) to four first-line antimicrobials-isoniazid, rifampin, ethambutol, and pyrazinamide-for MTBC isolates tested between the years 2018-2022. For this 5-year retrospective analysis, the WGS sensitivity for predicting resistance for isoniazid, rifampin, ethambutol, and pyrazinamide using Mykrobe was 86.7%, 100.0%, 100.0%, and 47.8%, respectively, and the specificity was 99.4%, 99.5%, 98.7%, and 99.9%, respectively. The predictive values improved slightly using Mykrobe corrections applied using TB Profiler, i.e., the WGS sensitivity for isoniazid, rifampin, ethambutol, and pyrazinamide was 92.31%, 100%, 100%, and 57.78%, respectively, and the specificity was 99.63%. 99.45%, 98.93%, and 99.93%, respectively. The utilization of WGS-based testing addresses concerns regarding test turnaround time and enables analysis for MTBC member identification, antimicrobial resistance prediction, detection of mixed cultures, and strain genotyping, all through a single laboratory test. WGS enables rapid resistance detection compared to traditional phenotypic susceptibility testing methods using the WHO TB mutation catalog, providing an insight into lesser-known mutations, which should be added to prediction databases as high-confidence mutations are recognized. The WGS-based methods can support TB elimination efforts in Canada and globally by ensuring the early start of appropriate treatment, rapidly limiting the spread of TB outbreaks.

Targeted next-generation sequencing to diagnose drug-resistant tuberculosis: a systematic review and meta-analysis

BACKGROUND

Targeted next-generation sequencing (NGS) can rapidly and simultaneously detect mutations associated with resistance to tuberculosis drugs across multiple gene targets. The use of targeted NGS to diagnose drug-resistant tuberculosis, as described in publicly available data, has not been comprehensively reviewed. We aimed to identify targeted NGS assays that diagnose drug-resistant tuberculosis, determine how widely this technology has been used, and assess the diagnostic accuracy of these assays.

METHODS

In this systematic review and meta-analysis, we searched MEDLINE, Embase, Cochrane Library, Web of Science Core Collection, Global Index Medicus, Google Scholar, ClinicalTrials.gov, and the WHO International Clinical Trials Registry Platform for published and unpublished reports on targeted NGS for drug-resistant tuberculosis from Jan 1, 2005, to Oct 14, 2022, with updates to our search in Embase and Google Scholar until Feb 13, 2024. Studies eligible for the systematic review described targeted NGS approaches to predict drug resistance in Mycobacterium tuberculosis infections using primary samples, reference strain collections, or cultured isolates from individuals with presumed or confirmed tuberculosis. Our search had no limitations on study type or language, although only reports in English, German, and French were screened for eligibility. For the meta-analysis, we included test accuracy studies that used any reference standard, and we assessed risk of bias using the Quality Assessment of Diagnostic Accuracy Studies-2 tool. The primary outcomes for the meta-analysis were sensitivity and specificity of targeted NGS to diagnose drug-resistant tuberculosis compared to phenotypic and genotypic drug susceptibility testing. We used a Bayesian bivariate model to generate summary receiver operating characteristic plots and diagnostic accuracy measures, overall and stratified by drug and sample type. This study is registered with PROSPERO, CRD42022368707.

FINDINGS

We identified and screened 2920 reports, of which 124 were eligible for our systematic review, including 37 review articles and 87 reports of studies collecting samples for targeted NGS. Sequencing was mainly done in the USA (14 [16%] of 87), western Europe (ten [11%]), India (ten [11%]), and China (nine [10%]). We included 24 test accuracy studies in the meta-analysis, in which 23 different tuberculosis drugs or drug groups were assessed, covering first-line drugs, injectable drugs, and fluoroquinolones and predominantly comparing targeted NGS with phenotypic drug susceptibility testing. The combined sensitivity of targeted NGS across all drugs was 94·1% (95% credible interval [CrI] 90·9-96·3) and specificity was 98·1% (97·0-98·9). Sensitivity for individual drugs ranged from 76·5% (52·5-92·3) for capreomycin to 99·1% (98·3-99·7) for rifampicin; specificity ranged from 93·1% (88·0-96·3) for ethambutol to 99·4% (98·3-99·8) for amikacin. Diagnostic accuracy was similar for primary clinical samples and culture isolates overall and for rifampicin, isoniazid, ethambutol, streptomycin, and fluoroquinolones, and similar after excluding studies at high risk of bias (overall sensitivity 95·2% [95% CrI 91·7-97·1] and specificity 98·6% [97·4-99·3]).

INTERPRETATION

Targeted NGS is highly sensitive and specific for detecting drug resistance across panels of tuberculosis drugs and can be performed directly on clinical samples. There is a paucity of data on performance for some currently recommended drugs. The barriers preventing the use of targeted NGS to diagnose drug-resistant tuberculosis in high-burden countries need to be addressed.

FUNDING

National Institutes of Allergy and Infectious Diseases and Swiss National Science Foundation.

Acquired rifamycin resistance among patients with tuberculosis and HIV in new York City, 2001-2023

Introduction

Acquired rifamycin resistance (ARR) in tuberculosis (TB) has been associated with HIV infection and can necessitate complicated TB treatment regimens, particularly in people living with HIV (PLWH). This work examines clinical characteristics and treatment outcomes of PLWH who developed ARR from 2001 to 2023 in New York City (NYC) to inform best practices for treating these patients.

Methods

PLWH who developed ARR 2001-2023 were identified from the NYC TB registry.

Results

Sixteen PLWH developed ARR; 15 were diagnosed 2001-2009 and the 16th was diagnosed in 2017. Median CD4 count was 48/mm3. On initial presentation, 14 had positive sputum cultures; of these, 12 culture-converted prior to developing ARR. Ten patients completed a course of TB treatment but subsequently relapsed; in six of these cases, ARR was discovered upon relapse, triggering treatment with a non-rifamycin-containing regimen, while in the other four, ARR was discovered during a second round of rifamycin-containing treatment. Three patients were lost to follow-up during their initial course of TB treatment and later returned to care; after being restarted on a rifamycin-containing regimen, ARR was discovered. Finally, three patients culture-converted during their first course of treatment but subsequently had cultures that grew rifamycin-resistant Mycobacterium tuberculosis prior to treatment completion, leading to changes in their treatment regimens. Among the 16 patients, eight died before being cured of TB, seven successfully completed treatment, and one was lost to follow-up.

Conclusions

PLWH should be monitored closely for the development of ARR during treatment for TB, and sputum culture conversion should be interpreted cautiously in this group. Collecting a final sputum sample may be especially important for PLWH, as treatment failure and relapse were common in this population. The decrease in the number of cases of ARR among PLWH during the study period may reflect the decrease in the total number of PLWH diagnosed with TB in NYC in recent years, improved immune status of PLWH due to increased uptake of antiretroviral drugs, and improvements in the way anti-TB regimens are designed for PLWH (such as recommending daily rather than intermittent rifamycin dosing).

Implementation of a high-throughput whole genome sequencing approach with the goal of maximizing efficiency and cost effectiveness to improve public health

This manuscript describes the development of a streamlined, cost-effective laboratory workflow to meet the demands of increased whole genome sequence (WGS) capacity while achieving mandated quality metrics. From 2020 to 2021, the Wadsworth Center Bacteriology Laboratory (WCBL) used a streamlined workflow to sequence 5,743 genomes that contributed sequence data to nine different projects. The combined use of the QIAcube HT, Illumina DNA Prep using quarter volume reactions, and the NextSeq allowed the WCBL to process all samples that required WGS while also achieving a median turn-around time of 7 days (range 4 to 10 days) and meeting minimum sequence quality requirements. Public Health Laboratories should consider implementing these methods to aid in meeting testing requirements within budgetary restrictions.

IMPORTANCE

Public Health Laboratories that implement whole genome sequencing (WGS) technologies may struggle to find the balance between sample volume and cost effectiveness. We present a method that allows for sequencing of a variety of bacterial isolates in a cost-effective manner. This report provides specific strategies to implement high-volume WGS, including an innovative, low-cost solution utilizing a novel quarter volume sequencing library preparation. The methods described support the use of high-throughput DNA extraction and WGS within budgetary constraints, strengthening public health responses to outbreaks and disease surveillance.

Multidrug-resistant Mycobacterium tuberculosis transmission in Shandong, China

Multidrug-resistant tuberculosis (MDR-TB) has imposed a significant economic and health burden worldwide, notably in China. Using whole genome sequence, we sought to understand the mutation and transmission of MDR-TB in Shandong. A retrospective study of patients diagnosed with pulmonary tuberculosis in Shandong from 2009 to 2018 was conducted. To explore transmission patterns, we performed whole genome sequencing on MDR-TB isolates, identified genomic clusters, and assessed the drug resistance of TB isolates. Our study analyzed 167 isolates of MDR-TB, finding that 100 were clustered. The predominant lineage among MDR-TB isolates was lineage 2, specifically with a notable 88.6% belonging to lineage 2.2.1. Lineage 4 constituted a smaller proportion, accounting for 4.2% of the isolates. We discovered that Shandong has a significant clustering percentage for MDR-TB, with Jining having the highest percentage among all Shandong cities. The clustering percentages of MDR-TB, pre-extensively drug-resistant tuberculosis, and extensively drug-resistant tuberculosis were 59.9%, 66.0%, and 71.4%, respectively, and the clustering percentages increased with the expansion of the anti-TB spectrum. Isolates from genomic clusters 1 and 3 belonged to lineage 2.2.1 and showed signs of cross-regional transmission. The distribution of rrs A1401G and katG S315T mutations in lineage 2.2.1 and 2.2.2 strains differed significantly (P < .05). MDR-TB isolates with rpoB I480V, embA-12C > T, and rrs A1401G mutations showed a higher likelihood of clustering (P < .05). Our findings indicate a significant problem of local transmission of MDR-TB in Shandong, China. Beijing lineage isolates and some drug-resistant mutations account for the MDR-TB transmission in Shandong.

Present and Future Non-Culture-Based Diagnostics: Stewardship Potentials and Considerations

The medical microbiologist plays a key role in the transition from culture-based to molecular test methods for diagnosis of infectious diseases. They must understand the scientific and technical bases underlying these tests along with their associated benefits and limitations and be able to educate administrators and patient providers on their proper use. Coordination of testing practices between clinical departments and the spectrum of public health and research laboratories is essential to optimize health care delivery.

Rifampin Mono-Resistant Tuberculosis in New York City, 2010-2021: A Retrospective Case Series

Background

Although relatively rare, rifampin mono-resistant tuberculosis (RMR TB) poses important challenges to effective TB treatment and control. Information on the burden of RMR TB and treatment outcomes is needed to inform diagnosis and management.

Methods

Standardized variables were collected from the New York City (NYC) tuberculosis surveillance system for patients treated for RMR TB in NYC during 2010-2021.

Results

Of 7097 TB cases reported in 2010-2021, 31 (<1%) were treated clinically as RMR TB. Five (16%) of these patients had HIV. Seventeen patients (55%) had TB that was rifampin-resistant by both molecular and phenotypic drug susceptibility testing; 2 (6%) had rifampin resistance by phenotypic tests, and molecular tests were not done; and 12 (39%) were identified based only on molecular tests. Among these 12, 7 were rifampin-sensitive by phenotypic tests, and phenotypic testing could not be done for the other 5. Ten of the 31 (32%) were diagnosed in 2010-2015; the other 21 (including 10/12 diagnosed by molecular tests alone) were diagnosed in 2016-2021. Of the 31 patients, 21 (68%) completed treatment (median treatment duration of 18 months). Although the interval between tuberculosis treatment initiation and change to a non-rifamycin-containing regimen decreased significantly during the study period, the overall duration of treatment did not decrease significantly between 2010 and 2021.

Conclusions

Molecular drug susceptibility tests identified cases of RMR TB that were not detected by phenotypic testing and helped enable timely adjustment of tuberculosis treatment regimens. Short-course regimens are needed to reduce duration of treatment for RMR TB.

Global status of phenotypic pyrazinamide resistance in Mycobacterium tuberculosis clinical isolates: an updated systematic review and meta-analysis

Pyrazinamide (PZA) is an essential first-line tuberculosis drug for its unique mechanism of action active against multidrug-resistant-TB (MDR-TB). Thus, the aim of updated meta-analysis was to estimate the PZA weighted pooled resistance (WPR) rate in M. tuberculosis isolates based on publication date and WHO regions. We systematically searched the related reports in PubMed, Scopus, and Embase (from January 2015 to July 2022). Statistical analyses were performed using STATA software. The 115 final reports in the analysis investigated phenotypic PZA resistance data. The WPR of PZA was 57% (95% CI 48-65%) in MDR-TB cases. According to the WHO regions, the higher WPRs of PZA were reported in the Western Pacific (32%; 95% CI 18-46%), South East Asian region (37%; 95% CI 31-43%), and the Eastern Mediterranean (78%; 95% CI 54-95%) among any-TB patients, high risk of MDR-TB patients, and MDR-TB patients, respectively. A negligible increase in the rate of PZA resistance were showed in MDR-TB cases (55% to 58%). The rate of PZA resistance has been rising in recent years among MDR-TB cases, underlines the essential for both standard and novel drug regimens development.

Clinical Whole-Genome Sequencing Assay for Rapid Mycobacterium tuberculosis Complex First-Line Drug Susceptibility Testing and Phylogenetic Relatedness Analysis

The global rise of drug resistant tuberculosis has highlighted the need for improved diagnostic technologies that provide rapid and reliable drug resistance results. Here, we develop and validate a whole genome sequencing (WGS)-based test for identification of mycobacterium tuberculosis complex (MTB) drug resistance to rifampin, isoniazid, pyrazinamide, ethambutol, and streptomycin. Through comparative analysis of drug resistance results from WGS-based testing and phenotypic drug susceptibility testing (DST) of 38 clinical MTB isolates from patients receiving care in Los Angeles, CA, we found an overall concordance between methods of 97.4% with equivalent performance across culture media. Critically, prospective analysis of 11 isolates showed that WGS-based testing provides results an average of 36 days faster than phenotypic culture-based methods. We showcase the additional benefits of WGS data by investigating a suspected laboratory contamination event and using phylogenetic analysis to search for cryptic local transmission, finding no evidence of community spread amongst our patient population in the past six years. WGS-based testing for MTB drug resistance has the potential to greatly improve diagnosis of drug resistant MTB by accelerating turnaround time while maintaining accuracy and providing additional benefits for infection control, lab safety, and public health applications.

Genomic revolution: Transforming tuberculosis diagnosis and treatment with the use of Whole Genome Sequencing - A consensus statement

Tuberculosis (TB) is a preventable, treatable, and curable disease. However, in 2020, 9∙9 million people were estimated to have developed tuberculosis, and 1.5 million people were estimated to have died from it. Whereas in India, 2.6 million were diagnosed with TB and 436,000 succumbed to TB in 2019. India (26%) is the major contributor to the global drop in TB cases. The COVID-19 pandemic has substantially reduced access to services for the diagnosis and treatment of TB, resulting in an increase in deaths and a reversal in global progress. [1] Presently, TB incidence is falling at a rate of 2% per year, obstructed mainly by the rearing pandemic of drug-resistant tuberculosis (DRTB). Particularly concerning is multi-drug resistant TB (MDRTB), defined as resistance towards isoniazid (INH) and rifampicin (RIF). [2] The World Health Organization (WHO) targeted to reduce worldwide TB incidence by 90% until 2035. (1) Early initiation of effective treatment based on susceptibility patterns of the Mycobacterium tuberculosis complex (MTBC) is considered key to successful TB control in countries with high DRTB incidence. Worldwide MDRTB treatment outcomes are poor, with cure rates less than 60% (2) due to the lack of comprehensive Drug Susceptibility Testing (DST) in most high MDRTB burden countries. This is leading to the inadequate anti-TB activity of the provided regimens (3-5), unlike regimens advised for DST assure optimal results. (6) In addition to resistances to the established regimens, the resistance to the newer DRTB drugs is increasing. On World TB Day 2022, Academy of Advanced Medical Education, Thyrocare Technologies Limited and HyastackAnalytics - IITB along with expert pulmonologist and renowned physicians from India convened for an advisory board meeting in Delhi on 20th March 2022 to discuss the role of Whole Genome Sequencing (WGS) in the diagnosis and management of TB. Objectives and specific topics relating to WGS in MDRTB were discussed, each expert shared their views, which led to a group discussion with a commitment to putting the patient first, and increasing their collective efforts, the organizations recognized that it is possible to make this goal a reality. The organizations involved in the discussion have declared their commitment to engaging in collaborative efforts to tackle DRTB detection efficiently. They advocate for strengthening access to WGS TB services, controlling and preventing TB, improving surveillance and drug resistance management, and investing in research and development. This Round Table serves as a framework to build on and ensure that the goal of ending TB is achievable with WGS services wherever needed. Post discussion, a uniform consensus was said to be arrived if more than 80% board members agreed to the statement. The present paper is the outcome of aspects presented and discussed in the advisory board meeting.

Analysis of Drug-Resistance Characteristics and Genetic Diversity of Multidrug-Resistant Tuberculosis Based on Whole-Genome Sequencing on the Hainan Island, China

Purpose

Given the high burden of Tuberculosis (TB) in China, the prevalence of multidrug-resistant tuberculosis (MDR-TB) is significant. Whole-genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) enables the identification of lineages, drug-resistant mutations, and transmission patterns, offering valuable insights for TB control, clinical diagnosis, and treatment.

Methods

We collected 202 MDR-MTB strains from 3519 suspected pulmonary TB patients treated at The Second Affiliated Hospital of Hainan Medical University between July 2019 and June 2021. Proportional drug-susceptibility testing was performed using 8 common anti-tuberculosis drugs. Subsequently, the genotypic drug resistance and genetic characteristics were analyzed by the WGS.

Results

Lineages are identified by TB-profiler revealed 202 MDR-MTB strains, showcasing three predominant lineages, with lineage 2 being the most prevalent. Close genomic relatedness analysis and evidence of MTB transmission led to the formation of 15 clusters comprising 42 isolates, resulting in a clustering rate of 20.8%. Novelty, lineage 2.1 (non-Beijing) accounted for 27.2% of the MDR-MTB strains, which is rare in China and Neighboring countries. Regarding first-line anti-TB drugs, genes associated with rifampicin resistance, primarily the rpoB gene, were detected in 200 strains (99.0%). Genes conferring resistance to isoniazid, ethambutol, and streptomycin were identified in 191 (94.5%), 125 (61.9%), and 100 (49.5%) strains, respectively. Among the second-line drugs, 97 (48.0%) strains exhibited genes encoding resistance to fluoroquinolones. Comparing the results to phenotypic drug susceptibility-based testing, the sensitivity of WGS for detecting resistance to each of the six drugs (rifampicin, isoniazid, ethambutol, ofloxacin, kanamycin, capreomycin) was 90% or higher. With the exception of ethambutol, the specificity of WGS prediction for the remaining drugs exceeded 88%.

Conclusion

Our study provides crucial insights into genetic mutation types, genetic diversity, and transmission of MDR-MTB on Hainan Island, serving as a significant reference for MDR-MTB surveillance and clinical decision-making.

Mycobacterium tuberculosis complex whole-genome sequencing in New York State: Implementation of a reduced phenotypic drug susceptibility testing algorithm

Whole-genome sequencing (WGS) can predict drug resistance and antimicrobial susceptibility in Mycobacterium tuberculosis complex (MTBC) and has shown promise in partially replacing culture-based phenotypic drug susceptibility testing (pDST). We performed a two-year side by side study comparing the prediction of drug resistance and antimicrobial susceptibility by WGS molecular DST (mDST) to pDST to determine resistance at the critical concentration by Mycobacterial Growth Indicator Tube (MGIT) and agar proportion testing. Negative predictive values of WGS results were consistently high for the first-line drugs: rifampin (99.9%), isoniazid (99.0%), pyrazinamide (98.5%), and ethambutol (99.8%); the rates of resistance to these drugs, among strains in our population, are 2.9%, 10.4%, 46.3%, and 2.3%, respectively. WGS results were available an average 8 days earlier than first-line MGIT pDST. Based on these findings, we implemented a new testing algorithm with an updated WGS workflow in which strains predicted pan-susceptible were no longer tested by pDST. This algorithm was applied to 1177 isolates between October 2018 and September 2020, eliminating pDST for 66.6% of samples and reducing pDST for an additional 22.0%. This algorithm change resulted in faster turnaround times and decreased cost while maintaining comprehensive antimicrobial susceptibility profiles of all culture-positive MTBC cases in New York.

Direct detection of drug-resistant Mycobacterium tuberculosis using targeted next generation sequencing

Mycobacterium tuberculosis complex (MTBC) infections are treated with combinations of antibiotics; however, these regimens are not as efficacious against multidrug and extensively drug resistant MTBC. Phenotypic (growth-based) drug susceptibility testing on slow growing bacteria like MTBC requires many weeks to months to complete, whereas sequencing-based approaches can predict drug resistance (DR) with reduced turnaround time. We sought to develop a multiplexed, targeted next generation sequencing (tNGS) assay that can predict DR and can be performed directly on clinical respiratory specimens. A multiplex PCR was designed to amplify a group of thirteen full-length genes and promoter regions with mutations known to be involved in resistance to first- and second-line MTBC drugs. Long-read amplicon libraries were sequenced with Oxford Nanopore Technologies platforms and high-confidence resistance mutations were identified in real-time using an in-house developed bioinformatics pipeline. Sensitivity, specificity, reproducibility, and accuracy of the tNGS assay was assessed as part of a clinical validation study. In total, tNGS was performed on 72 primary specimens and 55 MTBC-positive cultures and results were compared to clinical whole genome sequencing (WGS) performed on paired patient cultures. Complete or partial susceptibility profiles were generated from 82% of smear positive primary specimens and the resistance mutations identified by tNGS were 100% concordant with WGS. In addition to performing tNGS on primary clinical samples, this assay can be used to sequence MTBC cultures mixed with other mycobacterial species that would not yield WGS results. The assay can be effectively implemented in a clinical/diagnostic laboratory with a two to three day turnaround time and, even if batched weekly, tNGS results are available on average 15 days earlier than culture-derived WGS results. This study demonstrates that tNGS can reliably predict MTBC drug resistance directly from clinical specimens or cultures and provide critical information in a timely manner for the appropriate treatment of patients with DR tuberculosis.

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.

Analysis of Xpert MTB/RIF results in retested patients with very low initial bacterial loads: A retrospective study in China

BACKGROUND

The Xpert MTB/RIF (Xpert) assay has been widely used to diagnose suspected active tuberculosis (TB) and rifampicin-resistant TB cases. Despite its excellent performance record, false-positive Xpert rifampicin (RIF) resistance results are obtained for specimens with extremely low bacterial loads.

OBJECTIVE

We aimed to study the feasibility of repeat Xpert testing as a strategy for reducing the odds of obtaining false-positive results when testing paucibacillary TB patients.

METHODS

We enrolled previously tested TB patients with very low initial bacterial loads from May 2016 to February 2022 for Xpert retesting. A total of 251 TB patients were retested using the Xpert assay.

RESULTS

RIF resistance was noted in 65 (25.9 %) patients when tested by Xpert at initial diagnosis. Only 107 (42.6 %) of 251 patients tested positive for MTB when retested via Xpert. The majority (98.6 %) of RIF-susceptible cases were still susceptible to RIF when retested. Initial Xpert testing yielded 35 positive results for MTB in the RIF-resistant group, of whom 25 (71.4 %) still exhibited RIF resistance when retested. All culture-positive MTB isolates in the RIF-susceptible group were also RIF-susceptible by phenotypic DST. In the RIF-resistant group, 10 of 14 culture-positive MTB isolates exhibited RIF resistance, of which 4 isolates were deemed RIF-susceptible by phenotypic DST. The proportion of double mutations within the MTB rpoB RRDR sequence, as detected by hybridization of Xpert D and E probes, was significantly higher in the RIF-susceptible group than in the RIF-susceptible group.

CONCLUSIONS

Our results demonstrated that initial RIF-susceptible results were more accurate than RIF-resistant results. Additionally, patients with double mutations that delayed probe D/E hybridization were more likely to have false-positive Xpert results. Our findings emphasize that repeat Xpert MTB/RIF testing is necessary for TB patients with extremely low bacterial loads who are at high risk for RIF-resistant TB.

Performances of bioinformatics tools for the analysis of sequencing data of Mycobacterium tuberculosis complex strains

Whole-genome sequencing of Mycobacterium tuberculosis complex (MTBC) strains is a rapidly growing tool to obtain results regarding the resistance and phylogeny of the strains. We evaluated the performances of two bioinformatics tools for the analysis of whole-genome sequences of MTBC strains. Two hundred and twenty-seven MTBC strains were isolated and whole-genome sequenced at the laboratory of Avicenne Hospital between 2015 and 2021. We investigated the resistance and susceptibility status of strains using two online tools, Mykrobe and PhyResSE. We compared the genotypic and phenotypic resistance results obtained by drug susceptibility testing. Unlike with the Mykrobe tool, sequencing quality data were obtained using PhyResSE: average coverage of 98% and average depth of 119X. We found a similar concordance between phenotypic and genotypic results when determining susceptibility to first-line anti-tuberculosis drugs (95%) with both tools. The sensitivity and specificity of each tool compared to the phenotypic method were respectively 72% [52-87] and 98% [96-99] for Mykrobe and 76% [57-90] and 97% [94-99] for PhyResSE. Mykrobe and PhyResSE were easy to use and efficient. These platforms are accessible to people not trained in bioinformatics and constitute a complementary approach to phenotypic methods for the study of MTBC strains.

A machine-learning based model for automated recommendation of individualized treatment of rifampicin-resistant tuberculosis

Background

Rifampicin resistant tuberculosis remains a global health problem with almost half a million new cases annually. In high-income countries patients empirically start a standardized treatment regimen, followed by an individualized regimen guided by drug susceptibility test (DST) results. In most settings, DST information is not available or is limited to isoniazid and fluoroquinolones. Whole genome sequencing could more accurately guide individualized treatment as the full drug resistance profile is obtained with a single test. Whole genome sequencing has not reached its full potential for patient care, in part due to the complexity of translating a resistance profile into the most effective individualized regimen.

Methods

We developed a treatment recommender clinical decision support system (CDSS) and an accompanying web application for user-friendly recommendation of the optimal individualized treatment regimen to a clinician.

Results

Following expert stakeholder meetings and literature review, nine drug features and 14 treatment regimen features were identified and quantified. Using machine learning, a model was developed to predict the optimal treatment regimen based on a training set of 3895 treatment regimen-expert feedback pairs. The acceptability of the treatment recommender CDSS was assessed as part of a clinical trial and in a routine care setting. Within the clinical trial setting, all patients received the CDSS recommended treatment. In 8 of 20 cases, the initial recommendation was recomputed because of stock out, clinical contra-indication or toxicity. In routine care setting, physicians rejected the treatment recommendation in 7 out of 15 cases because it deviated from the national TB treatment guidelines. A survey indicated that the treatment recommender CDSS is easy to use and useful in clinical practice but requires digital infrastructure support and training.

Conclusions

Our findings suggest that global implementation of the novel treatment recommender CDSS holds the potential to improve treatment outcomes of rifampicin resistant tuberculosis.

Comparative genomics of drug-resistant strains of Mycobacterium tuberculosis in Ecuador

BACKGROUND

Tuberculosis is a serious infectious disease affecting millions of people. In spite of efforts to reduce the disease, increasing antibiotic resistance has contributed to persist in the top 10 causes of death worldwide. In fact, the increased cases of multi (MDR) and extreme drug resistance (XDR) worldwide remains the main challenge for tuberculosis control. Whole genome sequencing is a powerful tool for predicting drug resistance-related variants, studying lineages, tracking transmission, and defining outbreaks. This study presents the identification and characterization of resistant clinical isolates of Mycobacterium tuberculosis including a phylogenetic and molecular resistance profile study by sequencing the complete genome of 24 strains from different provinces of Ecuador.

RESULTS

Genomic sequencing was used to identify the variants causing resistance. A total of 15/21 isolates were identified as MDR, 4/21 as pre-XDR and 2/21 as XDR, with three isolates discarded due to low quality; the main sub-lineage was LAM (61.9%) and Haarlem (19%) but clades X, T and S were identified. Of the six pre-XDR and XDR strains, it is noteworthy that five come from females; four come from the LAM sub-lineage and two correspond to the X-class sub-lineage. A core genome of 3,750 genes, distributed in 295 subsystems, was determined. Among these, 64 proteins related to virulence and implicated in the pathogenicity of M. tuberculosis and 66 possible pharmacological targets stand out. Most variants result in nonsynonymous amino acid changes and the most frequent genotypes were identified as conferring resistance to rifampicin, isoniazid, ethambutol, para-aminosalicylic acid and streptomycin. However, an increase in the resistance to fluoroquinolones was detected.

CONCLUSION

This work shows for the first time the variability of circulating resistant strains between men and women in Ecuador, highlighting the usefulness of genomic sequencing for the identification of emerging resistance. In this regard, we found an increase in fluoroquinolone resistance. Further sampling effort is needed to determine the total variability and associations with the metadata obtained to generate better health policies.

The Changing Paradigm of Drug-Resistant Tuberculosis Treatment: Successes, Pitfalls, and Future Perspectives

Drug-resistant tuberculosis (DR-TB) remains a global crisis due to the increasing incidence of drug-resistant forms of the disease, gaps in detection and prevention, models of care, and limited treatment options. The DR-TB treatment landscape has evolved over the last 10 years. Recent developments include the remarkable activity demonstrated by the newly approved anti-TB drugs bedaquiline and pretomanid against Mycobacterium tuberculosis. Hence, treatment of DR-TB has drastically evolved with the introduction of the short-course regimen for multidrug-resistant TB (MDR-TB), transitioning to injection-free regimens and the approval of the 6-month short regimens for rifampin-resistant TB and MDR-TB. Moreover, numerous clinical trials are under way with the aim to reduce pill burden and shorten the DR-TB treatment duration. While there have been apparent successes in the field, some challenges remain. These include the ongoing inclusion of high-dose isoniazid in DR-TB regimens despite a lack of evidence for its efficacy and the inclusion of ethambutol and pyrazinamide in the standard short regimen despite known high levels of background resistance to both drugs. Furthermore, antimicrobial heteroresistance, extensive cavitary disease and intracavitary gradients, the emergence of bedaquiline resistance, and the lack of biomarkers to monitor DR-TB treatment response remain serious challenges to the sustained successes. In this review, we outline the impact of the new drugs and regimens on patient treatment outcomes, explore evidence underpinning current practices on regimen selection and duration, reflect on the disappointments and pitfalls in the field, and highlight key areas that require continued efforts toward improving treatment approaches and rapid biomarkers for monitoring treatment response.

Bedaquiline resistance probability to guide treatment decision making for rifampicin-resistant tuberculosis: insights from a qualitative study

BACKGROUND

Bedaquiline (BDQ) is a core drug for rifampicin-resistant tuberculosis (RR-TB) treatment. Accurate prediction of a BDQ-resistant phenotype from genomic data is not yet possible. A Bayesian method to predict BDQ resistance probability from next-generation sequencing data has been proposed as an alternative.

METHODS

We performed a qualitative study to investigate the decision-making of physicians when facing different levels of BDQ resistance probability. Fourteen semi-structured interviews were conducted with physicians experienced in treating RR-TB, sampled purposefully from eight countries with varying income levels and burden of RR-TB. Five simulated patient scenarios were used as a trigger for discussion. Factors influencing the decision of physicians to prescribe BDQ at macro-, meso- and micro levels were explored using thematic analysis.

RESULTS

The perception and interpretation of BDQ resistance probability values varied widely between physicians. The limited availability of other RR-TB drugs and the high cost of BDQ hindered physicians from altering the BDQ-containing regimen and incorporating BDQ resistance probability in their decision-making. The little experience with BDQ susceptibility testing and whole-genome sequencing results, and the discordance between phenotypic susceptibility and resistance probability were other barriers for physicians to interpret the resistance probability estimates. Especially for BDQ resistance probabilities between 25% and 70%, physicians interpreted the resistance probability value dynamically, and other factors such as clinical and bacteriological treatment response, history of exposure to BDQ, and resistance profile were often considered more important than the BDQ probability value for the decision to continue or stop BDQ. In this grey zone, some physicians opted to continue BDQ but added other drugs to strengthen the regimen.

CONCLUSIONS

This study highlights the complexity of physicians' decision-making regarding the use of BDQ in RR-TB regimens for different levels of BDQ resistance probability.. Ensuring sufficient access to BDQ and companion drugs, improving knowledge of the genotype-phenotype association for BDQ resistance, availability of a rapid molecular test, building next-generation sequencing capacity, and developing a clinical decision support system incorporating BDQ resistance probability will all be essential to facilitate the implementation of BDQ resistance probability in personalizing treatment for patients with RR-TB.

Whole genome sequencing for tuberculosis in Victoria, Australia: A genomic implementation study from 2017 to 2020

BACKGROUND

Whole genome sequencing (WGS) is increasingly used by tuberculosis (TB) programs to monitor Mycobacterium tuberculosis (Mtb) transmission. We aimed to characterise the molecular epidemiology of TB and Mtb transmission in the low-incidence setting of Victoria, Australia, and assess the utility of WGS.

METHODS

WGS was performed on all first Mtb isolates from TB cases from 2017 to 2020. Potential clusters (≤12 single nucleotide polymorphisms [SNPs]) were investigated for epidemiological links. Transmission events in highly-related (≤5 SNPs) clusters were classified as likely or possible, based on the presence or absence of an epidemiological link, respectively. Case characteristics and transmission settings (as defined by case relationship) were summarised. Poisson regression was used to examine associations with secondary case number.

FINDINGS

Of 1844 TB cases, 1276 (69.2%) had sequenced isolates, with 182 (14.2%) in 54 highly-related clusters, 2-40 cases in size. Following investigation, 140 cases (11.0% of sequenced) were classified as resulting from likely/possible local-transmission, including 82 (6.4%) for which transmission was likely. Common identified transmission settings were social/religious (26.4%), household (22.9%) and family living in different households (7.1%), but many were uncertain (41.4%). While household transmission featured in many clusters (n = 24), clusters were generally smaller (median = 3 cases) than the fewer that included transmission in social/religious settings (n = 12, median = 7.5 cases). Sputum-smear-positivity was associated with higher secondary case numbers.

INTERPRETATION

WGS results suggest Mtb transmission commonly occurs outside the household in our low-incidence setting. Further work is required to optimise the use of WGS in public health management of TB.

FUNDING

The Victorian Tuberculosis Program receives block funding for activities including case management and contact tracing from the Victorian Department of Health. No specific funding for this report was received by manuscript authors or the Victorian Tuberculosis Program, and the funders had no role in the study design, data collection, data analysis, interpretation or report writing.

Rapid Identification of Drug Resistance and Phylogeny in M. tuberculosis, Directly from Sputum Samples

Tuberculosis (TB) remains one of the most important infectious diseases globally. Establishing a resistance profile from the initial TB diagnosis is a priority. Rapid molecular tests evaluate only the most common genetic variants responsible for resistance to certain drugs, and Whole Genome Sequencing (WGS) needs culture prior to next-generation sequencing (NGS), limiting their clinical value. Targeted sequencing (TS) from clinical samples avoids these drawbacks, providing a signature of genetic markers that can be associated with drug resistance and phylogeny. In this study, a proof-of-concept protocol was developed for detecting genomic variants associated with drug resistance and for the phylogenetic classification of Mycobacterium Tuberculosis (Mtb) in sputum samples. Initially, a set of Mtb reference strains from the WHO were sequenced (WGS and TS). The results from the protocol agreed >95% with WHO reported data and phenotypic drug susceptibility testing (pDST). Lineage genetics results were 100% concordant with those derived from WGS. After that, the TS protocol was applied to sputum samples from TB patients to detect resistance to first- and second-line drugs and derive phylogeny. The accuracy was >90% for all evaluated drugs, except Eto/Pto (77.8%), and 100% were phylogenetically classified. The results indicate that the described protocol, which affords the complete drug resistance profile and phylogeny of Mtb from sputum, could be useful in the clinical area, advancing toward more personalized and more effective treatments in the near future. IMPORTANCE The COVID-19 pandemic negatively affected the progress in accessing essential Tuberculosis (TB) services and reducing the burden of TB disease, resulting in a decreased detection of new cases and increased deaths. Generating molecular diagnostic tests with faster results without losing reliability is considered a priority. Specifically, developing an antimicrobial resistance profile from the initial stages of TB diagnosis is essential to ensure appropriate treatment. Currently available rapid molecular tests evaluate only the most common genetic variants responsible for resistance to certain drugs, limiting their clinical value. In this work, targeted sequencing on sputum samples from TB patients was used to identify Mycobacterium tuberculosis mutations in genes associated with drug resistance and to derive a phylogeny of the infecting strain. This protocol constitutes a proof-of-concept toward the goal of helping clinicians select a timely and appropriate treatment by providing them with actionable information beyond current molecular approaches.

Characterization of Mycobacterium salfingeri sp. nov.: A novel nontuberculous mycobacteria isolated from a human wound infection

Nontuberculous mycobacteria (NTM) are environmental bacteria commonly found in soil and water in almost every part of the world. While usually non-pathogenic, they can cause acute respiratory and cutaneous infections under certain circumstances or in patients with underlying medical conditions. Contrary to members of the Mycobacterium tuberculosis complex, documented human-to-human transmissions of NTM have been rarely reported and most cases result from direct environmental exposure. Here we describe the identification of a new NTM species isolated from a hand laceration of a New York State patient after a fall. This new NTM forms rough, orange pigmented colonies and is naturally resistant to doxycycline and tobramycin. Whole genome analysis reveal no close relatives present in public databases, and our findings are in accordance with the recognition of a new taxonomic species of NTM. We propose the name Mycobacterium salfingeri sp. nov. for this new NTM representative. The type strain is 20-157661T (DSM = 113368T, BCCM = ITM 501207T).

Review and Updates on the Diagnosis of Tuberculosis

Diagnosis of tuberculosis, and especially the diagnosis of extrapulmonary tuberculosis, still faces challenges in clinical practice. There are several reasons for this. Methods based on the detection of Mycobacterium tuberculosis (Mtb) are insufficiently sensitive, methods based on the detection of Mtb-specific immune responses cannot always differentiate active disease from latent infection, and some of the serological markers of infection with Mtb are insufficiently specific to differentiate tuberculosis from other inflammatory diseases. New tools based on technologies such as flow cytometry, mass spectrometry, high-throughput sequencing, and artificial intelligence have the potential to solve this dilemma. The aim of this review was to provide an updated overview of current efforts to optimize classical diagnostic methods, as well as new molecular and other methodologies, for accurate diagnosis of patients with Mtb infection.

Use of Whole-Genome Sequencing to Explore Mycobacterium tuberculosis Complex Circulating in a Hotspot Department in France

The Seine-Saint-Denis is the French metropolitan department with the highest incidence of tuberculosis (TB). Our aim was to explore epidemiological and phylogenetic characteristics of TB strains in this hotspot department. We performed WGS on 227 strains of Mycobacterium tuberculosis complex isolated from patients at the Avicenne Hospital from 2016 to 2021 and randomly selected to represent the clinical diversity of French TB localization. Clinical and demographic data were recorded for each TB patient. The mean age of patients was 36 years old. They came from Africa (44%), Asia (27%), Europe (26%) and America (3%). Strains isolated from extrapulmonary samples were associated with Asian patients, whereas strains isolated from pulmonary samples were associated with European patients. We observed a high level of lineage diversity in line with the known worldwide diversity. Interestingly, lineage 3 was associated with lymph node TB. Additionally, the sensitivity of WGS for predicting resistance was 100% for rifampicin, isoniazid and ethambutol and 66.7% for pyrazinamide. The global concordance with drug-susceptibility testing using the phenotypic approach was 97%. In microbiology laboratories, WGS turns out to be an essential tool for better understanding local TB epidemiology, with direct access to circulating lineage identification and to drug susceptibilities to first- and second-line anti-TB drugs.

Novel Mycobacterium tuberculosis Complex Genotype Related to M. caprae

We report the unusual genotypic characterization of a bacterium isolated from a clinical sample of a patient who grew up in Bangladesh and lives in the United States. Using whole-genome sequencing, we identified the bacterium as a member of the Mycobacterium tuberculosis complex (MTBC). Phylogenetic placement of this strain suggests a new MTBC genotype. Even though it had the same spoligotype as M. caprae strains, single-nucleotide polymorphism–based phylogenetic analysis placed the isolate as a sister lineage distinct from M. caprae, most closely related to 5 previously sequenced genomes isolated from primates and elephants in Asia. We propose a new animal-associated lineage, La4, within MTBC.

Whole-Genome Sequencing for Resistance Level Prediction in Multidrug-Resistant Tuberculosis

Defining the precise relationship between resistance mutations and quantitative phenotypic drug susceptibility testing will increase the value of whole-genome sequencing (WGS) for predicting tuberculosis drug resistance. However, a large number of WGS data sets currently lack corresponding quantitative phenotypic data—the MICs. Using MYCOTBI plates, we determined the MICs to nine antituberculosis drugs for 154 clinical multidrug-resistant tuberculosis isolates from the Shenzhen Center for Chronic Disease Control in Shenzhen, China. Comparing MICs with predicted drug-resistance profiles inferred by WGS showed that WGS could predict the levels of resistance to isoniazid, rifampicin, streptomycin, fluoroquinolones, and aminoglycosides. We also found some mutations that may not be associated with drug resistance, such as EmbB D328G, mutations in the gid gene, and C−12T in the eis promoter. However, some strains carrying the same mutations showed different levels of resistance to the corresponding drugs. The MICs of different strains with the RpsL K88R, fabG1 C−15T mutations and some with mutations in embB and rpoB, had MICs to the corresponding drugs that varied by 8-fold or more. This variation is unexplained but could be influenced by the bacterial genetic background. Additionally, we found that 32.3% of rifampicin-resistant isolates were rifabutin-susceptible, particularly those with rpoB mutations H445D, H445L, H445S, D435V, D435F, L452P, S441Q, and S441V. Studying the influence of bacterial genetic background on the MIC and the relationship between rifampicin-resistant mutations and rifabutin resistance levels should improve the ability of WGS to guide the selection of medical treatment regimens.

IMPORTANCE Whole-genome sequencing (WGS) has excellent potential in drug-resistance prediction. The MICs are essential indications of adding a particular antituberculosis drug dosage or changing the entire treatment regimen. However, the relationship between many known drug-resistant mutations and MICs is unclear, especially for rarer ones. The results showed that WGS could predict resistance levels to isoniazid, rifampicin, streptomycin, fluoroquinolones, and aminoglycosides. However, some mutations may not be associated with drug resistance, and some others may confer various MICs to strains carrying them. Also, 32.3% of rifampicin (RIF)-resistant strains were classified as sensitive to rifabutin (RFB), and some mutations in the rpoB gene may be associated with this phenotype. Our data on the MIC distribution of strains with some rarer mutations add to the accumulated data on the resistance level associated with such mutations to help guide further research and draw meaningful conclusions.

Prediction of drug resistance profile of multidrug-resistant Mycobacterium tuberculosis (MDR-MTB) isolates from newly diagnosed case by whole genome sequencing (WGS): a study from a high tuberculosis burden country

OBJECTIVES

Our aim was to assess the ability of the Whole-genome sequencing (WGS) in predicting drug resistance profile of multidrug-resistant mycobacterium tuberculosis (MDR-MTB) from newly diagnosed cases in China.

METHODS

We validated the Phenotypic drug Sensitivity Test (pDST) for 12 anti-tuberculosis drugs using the Bactec MGIT 960 system. We described the characteristics of the isolates enrolled and compared the pDST results with resistance profiles predicted by WGS.

RESULTS

The pDST showed that of the 43 isolates enrolled, 25.6% were sensitive to rifabutin (RFB); 97.7%、97.7%、93.0% and 93.0% were sensitive to cycloserine (Cs), amikacin/kanamycin (Ak/Km), para-aminosalicylic acid (Pas) and ethionamide Eto), respectively; 18.6% were resistant to fluoroquinolones (FQs) or second-line injections. Genotype DST determined by WGS of Ak/Km、Eto and RFP reached high consistency to 97.7% compared with pDST, followed by moxifloxacin (Mfx) 95.3%, levofloxaci (Lfx) and Pas 93%, streptomycin (Sm) 90.3%. The genotype DST of RFB and EMB showed low consistency with the pDST of 67.2 and 79.1%. WGS also detected 27.9% isolates of pyrazinamide(PZA)-related drug-resistant mutation. No mutations associated with linezolid (Lzd), bedaquiline (Bdq) and clofazimine (Cfz) were detectd.

CONCLUSIONS

WGS has the potential to infer resistance profiles without time-consuming phenotypic methods, which could be provide a basis to formulate reasonable treatment in high TB burden areas.

The impact of frequently neglected model violations on bacterial recombination rate estimation: a case study in Mycobacterium canettii and Mycobacterium tuberculosis

Mycobacterium canettii is a causative agent of tuberculosis in humans, along with the members of the Mycobacterium tuberculosis complex. Frequently used as an outgroup to the M. tuberculosis complex in phylogenetic analyses, M. canettii is thought to offer the best proxy for the progenitor species that gave rise to the complex. Here, we leverage whole-genome sequencing data and biologically relevant population genomic models to compare the evolutionary dynamics driving variation in the recombining M. canettii with that in the nonrecombining M. tuberculosis complex, and discuss differences in observed genomic diversity in the light of expected levels of Hill-Robertson interference. In doing so, we highlight the methodological challenges of estimating recombination rates through traditional population genetic approaches using sequences called from populations of microorganisms and evaluate the likely mis-inference that arises owing to a neglect of common model violations including purifying selection, background selection, progeny skew, and population size change. In addition, we compare performance when full within-host polymorphism data are utilized, versus the more common approach of basing analyses on within-host consensus sequences.

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.

Use of Whole-Genome Sequencing to Predict Mycobacterium tuberculosis Complex Drug Resistance from Early Positive Liquid Cultures

Our objective was to evaluate the performance of whole-genome sequencing (WGS) from early positive liquid cultures for predicting Mycobacterium tuberculosis complex (MTBC) drug resistance. Clinical isolates were obtained from tuberculosis patients at Shanghai Pulmonary Hospital (SPH). Antimicrobial susceptibility testing (AST) was performed, and WGS from early Bactec mycobacterial growth indicator tube (MGIT) 960-positive liquid cultures was performed to predict the drug resistance using the TB-Profiler informatics platform. A total of 182 clinical isolates were enrolled in this study. Using phenotypic AST as the gold standard, the overall sensitivity and specificity for WGS were, respectively, 97.1% (89.8 to 99.6%) and 90.4% (83.4 to 95.1%) for rifampin, 91.0% (82.4 to 96.3%) and 95.2% (89.1 to 98.4%) for isoniazid, 100.0% (89.4 to 100.0%) and 87.3% (80.8 to 92.1%) for ethambutol, 96.6% (88.3 to 99.6%) and 61.8% (52.6 to 70.4%) for streptomycin, 86.8% (71.9 to 95.6%) and 95.8% (91.2 to 98.5%) for moxifloxacin, 86.5% (71.2 to 91.5%) and 95.2% (90.3 to 98.0%) for ofloxacin, 100.0% (54.1 to 100.0%) and 67.6% (60.2 to 74.5%) for amikacin, 100.0% (63.1 to 100.0%) and 67.2% (59.7 to 74.2%) for kanamycin, 62.5% (24.5 to 91.5%) and 88.5% (82.8 to 92.8%) for ethionamide, 33.3% (4.3 to 77.7%) and 98.3% (95.1 to 99.7%) for para-aminosalicylic acid, and 0.0% (0.0 to 12.3%) and 100.0% (97.6 to 100.0%) for cycloserine. The concordances of WGS-based AST and phenotypic AST were as follows: rifampin (92.9%), isoniazid (93.4%), ethambutol (89.6%), streptomycin (73.1%), moxifloxacin (94.0%), ofloxacin (93.4%), amikacin (68.7%), kanamycin (68.7%), ethionamide (87.4%), para-aminosalicylic acid (96.2%) and cycloserine (84.6%). We conclude that WGS could be a promising approach to predict MTBC resistance from early positive liquid cultures.

IMPORTANCE In this study, we used whole-genome sequencing (WGS) from early positive liquid (MGIT) cultures instead of solid cultures to predict drug resistance of 182 Mycobacterium tuberculosis complex (MTBC) clinical isolates to predict drug resistance using the TB-Profiler informatics platform. Our study indicates that WGS may be a promising method for predicting MTBC resistance using early positive liquid cultures.

Rapid Identification of Mycobacterium tuberculosis Complex Using Mass Spectrometry: A Proof of Concept

Mycobacteria that form the Mycobacterium tuberculosis complex are responsible for deadly tuberculosis in animals and patients. Identification of these pathogens at the species level is of primary importance for treatment and source tracing and currently relies on DNA analysis, including whole genome sequencing (WGS), which requires a whole day. In this study, we report the unprecedented discrimination of M. tuberculosis complex species using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF-MS), with WGS as the comparative reference standard. In the first step, optimized peptide extraction applied to 36 isolates otherwise identified in five of the 11 M. tuberculosis complex variants by WGS yielded 139 MALDI-TOF spectra, which were used to identify biomarkers of interest that facilitate differentiation between variants. In a second step, 70/80 (88%) other isolates were correctly classified by an algorithm based on specific peaks. This study is the first to report a MALDI-TOF-MS method for discriminating M. tuberculosis complex mycobacteria that is easily implemented in clinical microbiology laboratories.

SAM-TB: a whole genome sequencing data analysis website for detection of Mycobacterium tuberculosis drug resistance and transmission

Whole genome sequencing (WGS) can provide insight into drug-resistance, transmission chains and the identification of outbreaks, but data analysis remains an obstacle to its routine clinical use. Although several drug-resistance prediction tools have appeared, until now no website integrates drug-resistance prediction with strain genetic relationships and species identification of nontuberculous mycobacteria (NTM). We have established a free, function-rich, user-friendly online platform for MTB WGS data analysis (SAM-TB, http://samtb.szmbzx.com) that integrates drug-resistance prediction for 17 antituberculosis drugs, detection of variants, analysis of genetic relationships and NTM species identification. The accuracy of SAM-TB in predicting drug-resistance was assessed using 3177 sequenced clinical isolates with results of phenotypic drug-susceptibility tests (pDST). Compared to pDST, the sensitivity of SAM-TB for detecting multidrug-resistant tuberculosis was 93.9% [95% confidence interval (CI) 92.6–95.1%] with specificity of 96.2% (95% CI 95.2–97.1%). SAM-TB also analyzes the genetic relationships between multiple strains by reconstructing phylogenetic trees and calculating pairwise single nucleotide polymorphism (SNP) distances to identify genomic clusters. The incorporated mlstverse software identifies NTM species with an accuracy of 98.2% and Kraken2 software can detect mixed MTB and NTM samples. SAM-TB also has the capacity to share both sequence data and analysis between users. SAM-TB is a multifunctional integrated website that uses WGS raw data to accurately predict antituberculosis drug-resistance profiles, analyze genetic relationships between multiple strains and identify NTM species and mixed samples containing both NTM and MTB. SAM-TB is a useful tool for guiding both treatment and epidemiological investigation.

Whole genome sequencing has the potential to improve treatment for rifampicin-resistant tuberculosis in high burden settings: a retrospective cohort study

Background Treatment of multidrug-resistant or rifampicin-resistant tuberculosis (MDR/RR-TB), although improved in recent years with shorter, more tolerable regimens, remains largely standardised and based on limited drug susceptibility testing (DST). More individualised treatment with expanded DST access is likely to improve patient outcomes. Methods To assess the potential of TB drug resistance prediction based on whole genome sequencing (WGS) to provide more effective treatment regimens, we applied current South African treatment recommendations to a retrospective cohort of MDR/RR-TB patients from Khayelitsha, Cape Town. Routine DST and clinical data were used to retrospectively categorise patients into a recommended regimen, either a standardised short regimen or a longer individualised regimen. Potential regimen changes were then described with the addition of WGS-derived DST. Findings WGS data were available for 1274 MDR/RR-TB patient treatment episodes across 2008-2017. Among 834 patients initially eligible for the shorter regimen, 385 (46%) may have benefited from reduced drug dosage or removing ineffective drugs when WGS data were considered. A further 187 (22%) may have benefited from more effective adjusted regimens. Among 440 patients initially eligible for a longer individualised regimen, 153 (35%) could have been switched to the short regimen. Overall, 305 (24%) patients had MDR/RR-TB with second-line TB drug resistance, where the availability of WGS-derived DST would have allowed more effective treatment individualisation. Interpretation These data suggest considerable benefits could accrue from routine access to WGS-derived resistance prediction. Advances in culture-free sequencing and expansion of the reference resistance mutation catalogue will increase the utility of WGS resistance prediction. Funding Swiss National Science Foundation, South African National Research Foundation, and Wellcome Trust.

A Recombinase Polymerase Amplification-Coupled Cas12a Mutant-Based Module for Efficient Detection of Streptomycin-Resistant Mutations in Mycobacterium tuberculosis

Drug-resistant tuberculosis (TB) is a serious public health problem and threat to global TB prevention and control. Streptomycin (STR) is the earliest and classical anti-TB drug, and it is the earliest drug that generated resistance to anti-TB treatment, which limits its use in treating TB and impedes TB control efforts. The rapid, economical, and highly sensitive detection of STR-resistant TB may help reduce disease transmission and morbimortality. CRISPR/CRISPR-associated protein (Cas) is a new-generation pathogen detection method that can detect single-nucleotide polymorphisms with high sensitivity and good specificity. In this study, a Cas12a RR detection system that can recognize more non-traditional protospacer-adjacent motif-targeting sequences was developed based on Cas12a combined with recombinase polymerase amplification technology. This system detects 0.1% of the target substance, and the entire detection process can be completed within 60 min. Its sensitivity and specificity for detecting clinical STR-resistant Mycobacterium tuberculosis were both 100%. Overall, the Cas12 RR detection system provides a novel alternative for the rapid, simple, sensitive, and specific detection of STR-resistant TB, which may contribute to the prompt treatment and prevention of disease transmission in STR-resistant TB.

Tuberculosis drug resistance profiling based on machine learning: A literature review

Tuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), is one of the top 10 causes of death worldwide. Drug-resistant tuberculosis (DR-TB) poses a major threat to the World Health Organization's “End TB” strategy which has defined its target as the year 2035. In 2019, there were close to 0.5 million cases of DRTB, of which 78% were resistant to multiple TB drugs. The traditional culture-based drug susceptibility test (DST - the current gold standard) often takes multiple weeks and the necessary laboratory facilities are not readily available in low-income countries. Whole genome sequencing (WGS) technology is rapidly becoming an important tool in clinical and research applications including transmission detection or prediction of DR-TB. For the latter, many tools have recently been developed using curated database(s) of known resistance conferring mutations. However, documenting all the mutations and their effect is a time-taking and a continuous process and therefore Machine Learning (ML) techniques can be useful for predicting the presence of DR-TB based on WGS data. This can pave the way to an earlier detection of drug resistance and consequently more efficient treatment when compared to the traditional DST.

The whole-genome sequencing in predicting Mycobacterium tuberculosis drug susceptibility and resistance in Papua, Indonesia

BACKGROUND

Tuberculosis is one of the deadliest disease caused by Mycobacterium tuberculosis. Its treatment still becomes a burden for many countries including Indonesia. Drug resistance is one of the problems in TB treatment. However, a development in the molecular field through Whole-genome sequencing (WGS) can be used as a solution in detecting mutations associated with TB- drugs. This investigation intended to implement this data for supporting the scientific community in deeply understanding any TB epidemiology and evolution in Papua along with detecting any mutations in genes associated with TB-Drugs.

RESULT

A whole-genome sequencing was performed on the random samples from TB Referral Laboratory in Papua utilizing MiSeq 600 cycle Reagent Kit (V3). Furthermore, TBProfiler was used for genome analysis, RAST Server was employed for annotation, while Gview server was applied for BLAST genome mapping and a Microscope server was implemented for Regions of Genomic Plasticity (RGP). The largest genome of M. tuberculosis obtained was at the size of 4,396,040 bp with subsystems number at 309 and the number of coding sequences at 4326. One sample (TB751) contained one RGP. The drug resistance analysis revealed that several mutations associated with TB-drug resistance existed. In details, mutations of rpoB gene which were identified as S450L, D435Y, H445Y, L430P, and Q432K had caused the reduced effectiveness of rifampicin; while the mutases in katG (S315T), kasA (312S), inhA (I21V), and Rv1482c-fabG1 (C-15 T) genes had contributed to the resistance in isoniazid. In streptomycin, the resistance was triggered by the mutations in rpsL (K43R) and rrs (A514C, A514T) genes, and, in Amikacin, its resistance was led by mutations in rrs (A514C) gene. Additionally, in Ethambutol and Pyrazinamide, their reduced effectiveness was provoked by embB gene mutases (M306L, M306V, D1024N) and pncA (W119R).

CONCLUSIONS

The results from whole-genome sequencing of TB clinical sample in Papua, Indonesia could contribute to the surveillance of TB-drug resistance. In the drug resistance profile, there were 15 Multi Drugs Resistance (MDR) samples. However, Extensively Drug-resistant (XDR) samples have not been found, but samples were resistant to only Amikacin, a second-line drug.

Utilization of whole genome sequencing for resolution of discrepant Mycobacterium tuberculosis drug susceptibility results: A case report

A 44-year-old woman undergoing therapy for acute promyelocytic leukemia (APL) developed disseminated tuberculosis. Mycobacterium tuberculosis (TB) was isolated from the blood and sputum. Initial drug susceptibility testing (DST) of the blood isolate revealed resistance to isoniazid and ethambutol but the sputum isolate showed no resistance. Due to drug resistance concerns, the patient was treated with multiple second and third-line drugs, and suffered from drug side effects. To further investigate the DST discrepancies, whole genome sequencing (WGS) was performed on both isolates. No known resistance mutations to first line or second line drugs were identified in either isolate, which was confirmed by additional susceptibility testing performed by a different reference laboratory and the California Department of Public Health (CDPH) laboratory. Treatment was reduced to a simpler and less toxic regimen due to these investigations. WGS is shown to be a valuable tool for resolving discordant phenotypic DST results of TB isolates and has the potential to provide accurate and timely results guiding appropriate therapy in the clinical setting.

Retrospective evaluation of routine whole genome sequencing of Mycobacterium tuberculosis at the Belgian National Reference Center, 2019

OBJECTIVES

Since January 2019, the Belgian National Reference Center for Mycobacteria (NRC) has switched from conventional typing to prospective whole-genome sequencing (WGS) of all submitted Mycobacterium tuberculosis complex (MTB) isolates. The ISO17025 validated procedure starts with semi-automated extraction and purification of gDNA directly from the submitted MGIT tubes, without preceding subculturing. All samples are then sequenced on an Illumina MiSeq sequencer and analyzed using an in-house developed and validated bioinformatics workflow to determine the species and antimicrobial resistance. In this study, we retrospectively compare results obtained via WGS to conventional phenotypic and genotypic testing, for all Belgian MTB strains analyzed in 2019 (n = 306).

RESULTS

In all cases, the WGS-based procedure was able to identify correctly the MTB species. Compared to MGIT drug susceptibility testing (DST), the sensitivity and specificity of genetic prediction of resistance to first-line antibiotics were respectively 100 and 99% (rifampicin, RIF), 90.5 and 100% (isoniazid, INH), 100 and 98% (ethambutol, EMB) and 61.1 and 100% (pyrazinamide, PZA). The negative predictive value was above 95% for these four first-line drugs. A positive predictive value of 100% was calculated for INH and PZA, 80% for RIF and 45% for EMB.

CONCLUSIONS

Our study confirms the effectiveness of WGS for the rapid detection of M. tuberculosis complex and its drug resistance profiles for first-line drugs even when working directly on MGIT tubes, and supports the introduction of this test into the routine workflow of laboratories performing tuberculosis diagnosis.

TB or not TB? Definitive determination of species within the Mycobacterium tuberculosis complex in unprocessed sputum from adults with presumed multidrug-resistant tuberculosis

Objectives

Differences among Mycobacterium tuberculosis complex (MTC) species may predict drug resistance or treatment success. Thus, we optimised and deployed the genotype MTBC assay (gMTBC) to identify MTC to the species level, and then performed comparative genotypic drug‐susceptibility testing to anti‐tuberculosis drugs from direct sputum of patients with presumed multidrug‐resistant tuberculosis (MDR‐TB) by the MTBDRplus/sl reference method.

Methods

Patients with positive Xpert^® MTB/RIF (Xpert) results were consented to provide early‐morning‐sputum for testing by the gMTBC and the reference MTBDRplus/sl. Chi‐square or Fisher’s exact test compared proportions. Modified Poisson regression modelled detection of MTC by gMTBC.

Results

Among 73 patients, 53 (73%) were male and had a mean age of 43 (95% CI; 40–45) years. In total, 34 (47%), 36 (49%) and 38 (55%) had positive gMTBC, culture and MTBDR respectively. Forty patients (55%) had low quantity MTC by Xpert, including 31 (78%) with a negative culture. gMTBC was more likely to be positive in patients with chest cavity 4.18 (1.31–13.32, P = 0.016), high‐quantity MTC by Xpert 3.03 (1.35–6.82, P = 0.007) and sputum smear positivity 1.93 (1.19–3.14, P = 0.008). The accuracy of gMTBC in detecting MTC was 95% (95% CI; 86–98; κ = 0.89) compared to MTBDRplus/sl. All M. tuberculosis/canettii identified by gMTB were susceptible to fluoroquinolone and aminoglycosides/capreomycin.

Conclusions

The concordance between the gMTBC assay and MTBDRplus/sl in detecting MTC was high but lagged behind the yield of Xpert MTB/RIF. All M. tuberculosis/canettii were susceptible to fluoroquinolones, a core drug in MDR‐TB treatment regimens.

Molecular epidemiology and drug-resistance of tuberculosis in Luodian revealed by whole genome sequencing

In this study, we aimed to investigate the molecular epidemiology and drug-resistance profiles of tuberculosis (TB) in Luodian, an area with highest TB incidence and limited healthcare resources in Guizhou, China. The passive case finding strategy was used to identify suspected pulmonary TB with symptoms, and individuals with positive Mycobacterium tuberculosis (MTB) culture were enrolled from May 22, 2018 to April 21, 2019. All the 107 cases except three came from nine towns, including 55.1% from Longping and Bianyang. The phylogeny tree showed that 53.3% of strains were Lineage 2 (Beijing genotype), while 46.7% were Lineage 4 (Euro-American genotype). Among Lineage 2 strains, 66.7% were of "modern" Beijing type. Seven clusters with genomic distance within 12 SNPs were identified. The clusters included 14 strains, accounting for a clustering rate of 13.1%. The distance separating the clustered cases was between 2.1 and 71.0 km (Km), with an average paired distance of 21.8 Km (interquartile range, 2.8-38.0 Km). Based on the gene mutations associated with drug-resistance, we predicted that 4.8% of strains were resistant to isoniazid, 3.7% to rifampicin, and 3.7% to streptomycin; only one strain (0.9%) had multidrug resistance (MDR). This study found low drug-resistance rates in Luodian, and the sub-lineage of the "modern" Beijing branch has recent expansion in Luodian. This work may also serve as a genomic baseline to assess the evolution and spread of MTB in Guizhou.

Genomic Prediction of Antimicrobial Resistance: Ready or Not, Here It Comes!

BACKGROUND

Next-generation sequencing (NGS) technologies are being used to predict antimicrobial resistance. The field is evolving rapidly and transitioning out of the research setting into clinical use. Clinical laboratories are evaluating the accuracy and utility of genomic resistance prediction, including methods for NGS, downstream bioinformatic pipeline components, and the clinical settings in which this type of testing should be offered.

CONTENT

We describe genomic sequencing as it pertains to predicting antimicrobial resistance in clinical isolates and samples. We elaborate on current methodologies and workflows to perform this testing and summarize the current state of genomic resistance prediction in clinical settings. To highlight this aspect, we include 3 medically relevant microorganism exemplars: Mycobacterium tuberculosis, Staphylococcus aureus, and Neisseria gonorrhoeae. Last, we discuss the future of genomic-based resistance detection in clinical microbiology laboratories.

SUMMARY

Antimicrobial resistance prediction by genomic approaches is in its infancy for routine patient care. Genomic approaches have already added value to the current diagnostic testing landscape in specific circumstances and will play an increasingly important role in diagnostic microbiology. Future advancements will shorten turnaround time, reduce costs, and improve our analysis and interpretation of clinically actionable results.

A Bioinformatics Whole-Genome Sequencing Workflow for Clinical Mycobacterium tuberculosis Complex Isolate Analysis, Validated Using a Reference Collection Extensively Characterized with Conventional Methods and In Silico Approaches

The use of whole-genome sequencing (WGS) for routine typing of bacterial isolates has increased substantially in recent years. For Mycobacterium tuberculosis (MTB), in particular, WGS has the benefit of drastically reducing the time required to generate results compared to most conventional phenotypic methods. Consequently, a multitude of solutions for analyzing WGS MTB data have been developed, but their successful integration in clinical and national reference laboratories is hindered by the requirement for their validation, for which a consensus framework is still largely absent. We developed a bioinformatics workflow for (Illumina) WGS-based routine typing of MTB complex (MTBC) member isolates allowing complete characterization, including (sub)species confirmation and identification (16S, csb/RD, hsp65), single nucleotide polymorphism (SNP)-based antimicrobial resistance (AMR) prediction, and pathogen typing (spoligotyping, SNP barcoding, and core genome multilocus sequence typing). Workflow performance was validated on a per-assay basis using a collection of 238 in-house-sequenced MTBC isolates, extensively characterized with conventional molecular biology-based approaches supplemented with public data. For SNP-based AMR prediction, results from molecular genotyping methods were supplemented with in silico modified data sets, allowing us to greatly increase the set of evaluated mutations. The workflow demonstrated very high performance with performance metrics of >99% for all assays, except for spoligotyping, where sensitivity dropped to ∼90%. The validation framework for our WGS-based bioinformatics workflow can aid in the standardization of bioinformatics tools by the MTB community and other SNP-based applications regardless of the targeted pathogen(s). The bioinformatics workflow is available for academic and nonprofit use through the Galaxy instance of our institute at https://galaxy.sciensano.be.

Value of routine whole genome sequencing for Mycobacterium tuberculosis drug resistance detection

Routine whole genome sequencing (WGS) of pathogens is becoming more feasible as sequencing costs decrease and access to benchtop sequencing equipment and bioinformatics pipelines increases. This study examined the added value gained from implementing routine WGS of all Mycobacterium tuberculosis isolates in New South Wales, Australia. Drug resistance markers inferred from WGS data were compared to commercial genotypic drug susceptibility testing (DST) assays and conventional phenotypic DST in all isolates sequenced between 2016 and 2019. Of the 1107 clinical M. tuberculosis isolates sequenced, 29 (2.6%) were multi-drug resistant (MDR); most belonged to Beijing (336; 30.4%) or East-African Indian (332; 30%) lineages. Compared with conventional phenotypic DST, WGS identified an additional 1% of isolates which were likely drug resistant, explained by mutations previously associated with treatment failure and mixed bacterial populations. However, WGS provided a 20% increase in drug resistance detection in comparison with commercial genotypic assays by identifying mutations outside of the classic resistance determining regions in rpoB, inhA, katG, pncA and embB genes. Gains in drug resistance detection were significant (p = 0.0137, paired t-test), but varied substantially for different phylogenetic lineages. In low incidence settings, routine WGS of M. tuberculosis provides better guidance for person-centered management of drug resistant tuberculosis than commercial genotypic assays.

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.