Phylogenomic analysis and Mycobacterium tuberculosis antibiotic resistance prediction by whole-genome sequencing from clinical isolates of Caldas, Colombia

A precision overview of genomic resistance screening in Ecuadorian isolates of Mycobacterium tuberculosis using web-based bioinformatics tools

INTRODUCTION

Tuberculosis (TB) is among the deadliest diseases worldwide, and its impact is mainly due to the continuous emergence of resistant isolates during treatment due to the laborious process of resistance diagnosis, nonadherence to treatment and circulation of previously resistant isolates of Mycobacterium tuberculosis. In this study, we evaluated the performance and functionalities of web-based tools, including Mykrobe, TB-profiler, PhyResSE, KvarQ, and SAM-TB, for detecting resistance in 88 Ecuadorian isolates of Mycobacterium tuberculosis drug susceptibility tested previously. Statistical analysis was used to determine the correlation between genomic and phenotypic analysis. Our results showed that with the exception of KvarQ, all tools had the highest correlation with the conventional drug susceptibility test (DST) for global resistance detection (98% agreement and 0.941 Cohen's kappa), while SAM-TB, PhyResSE, TB-profiler and Mykrobe had better correlations with DST for first-line drug analysis individually. We also identified that in our study, only 50% of mutations characterized by the web-based tools in the rpoB, katG, embB, pncA, gyrA and rrs regions were canonical and included in the World Health Organization (WHO) catalogue. Our findings suggest that SAM-TB, PhyResSE, TB-profiler and Mykrobe were efficient in determining canonical resistance-related mutations, but more analysis is needed to improve second-line detection. Improving surveillance programs using whole-genome sequencing tools for first-line drugs, MDR-TB and XDR-TB is essential to understand the molecular epidemiology of TB in Ecuador.

IMPORTANCE

Tuberculosis, an infectious disease caused by Mycobacterium tuberculosis, most commonly affects the lungs and is often spread through the air when infected people cough, sneeze, or spit. However, despite the existence of effective drug treatment, patient adherence, long duration of treatment, and late diagnosis have reduced the effectiveness of therapy and increased drug resistance. The increase in resistant cases, added to the impact of the COVID-19 pandemic, has highlighted the importance of implementing efficient and timely diagnostic methodologies worldwide. The significance of our research is in evaluating and identifying a more efficient and user-friendly web-based tool to characterize resistance in Mycobacterium tuberculosis by whole-genome sequencing, which will allow more routine application to improve TB strain surveillance programs locally.

Whole-Genome Sequencing of Mycobacterium tuberculosis Isolates from Diabetic and Non-Diabetic Patients with Pulmonary Tuberculosis

The association of tuberculosis and type 2 diabetes mellitus has been a recognized re-emerging challenge in management of the convergence of the two epidemics. Though much of the literature has studied this association, there is less knowledge in the field of genetic diversities that might occur in strains infecting tuberculosis patients with and without diabetes. Our study focused on determining the extent of diversity of genotypes of Mycobacterium tuberculosis in both these categories of patients. We subjected 55 M. tuberculosis isolates from patients diagnosed with pulmonary TB with and without type 2 diabetes mellitus to whole-genome sequencing on Illumina Hi Seq platform. The most common lineage identified was lineage 1, the Indo-Oceanic lineage (n = 22%), followed by lineage 4, the Euro-American lineage (n = 18, 33%); lineage 3, the East-African Indian lineage (n = 13, 24%); and lineage 2, the East-Asian lineage (n = 1, 2%). There were no significant differences in the distribution of lineages in both diabetics and non-diabetics in the South Indian population, and further studies involving computational analysis and comparative transcriptomics are needed to provide deeper insights.

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.

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.