Whole genome sequencing, analyses of drug resistance-conferring mutations, and correlation with transmission of Mycobacterium tuberculosis carrying katG-S315T in Hanoi, Vietnam

Word-based GWAS harnesses the rich potential of genomic data for E. coli quinolone resistance

Quinolone resistance presents a growing global health threat. We employed word-based GWAS to explore genomic data, aiming to enhance our understanding of this phenomenon. Unlike traditional variant-based GWAS analyses, this approach simultaneously captures multiple genomic factors, including single and interacting resistance mutations and genes. Analyzing a dataset of 92 genomic E. coli samples from a wastewater treatment plant in Dresden, we identified 54 DNA unitigs significantly associated with quinolone resistance. Remarkably, our analysis not only validated known mutations in gyrA and parC genes and the results of our variant-based GWAS but also revealed new (mutated) genes such as mdfA, the AcrEF-TolC multidrug efflux system, ptrB, and hisI, implicated in antibiotic resistance. Furthermore, our study identified joint mutations in 14 genes including the known gyrA gene, providing insights into potential synergistic effects contributing to quinolone resistance. These findings showcase the exceptional capabilities of word-based GWAS in unraveling the intricate genomic foundations of quinolone resistance.

Rv3539 (PPE63) of Mycobacterium Tuberculosis Promotes Survival of Mycobacterium Smegmatis in Human Macrophages Cell Line via Cell Wall Modulation of Bacteria and Altering Host's Immune Response

The modulation of host's immune response plays an important role in the intracellular survival of Mycobacterium tuberculosis. The intracellular pathogen counteracts environmental stresses with help of the expression of several genes. The M. tuberculosis genome encodes several immune-modulatory proteins including PE (proline-glutamic acid)/PPE (proline-proline-glutamic acid) superfamily proteins. It is unclear how the unique PE/PPE proteins superfamily contributes to survival under different stress and pathophysiology conditions. Previously, we showed that PPE63 (Rv3539) has C-terminal esterase extension and was localized as a membrane attached and in extracellular compartment. Therefore, the probability of these proteins interacting with the host to modulate the host immune response cannot be ruled out. The physiological role of PPE63 was characterized by expressing the PPE63 in the M. smegmatis, a non-pathogenic strain intrinsically deficient of PPE63. The recombinant M. smegmatis expressing PPE63 altered the colony morphology, lipid composition, and integrity of the cell wall. It provided resistance to multiple hostile environmental stress conditions and several antibiotics. MS_Rv3539 demonstrated higher infection and intracellular survival in comparison to the MS_Vec in the PMA-differentiated THP-1 cells. The decreased intracellular level of ROS, NO, and expression of iNOS was observed in THP-1 cells upon infection with MS_Rv3539 in comparison to MS_Vec. Further, the decrease in expression of pro-inflammatory cytokines like IL-6, TNF-α, and IL-1β and enhanced anti-inflammatory cytokines like IL-10, pointed toward its role in immune modulation. Overall this study suggested the role of Rv3539 in enhanced intracellular survival of M. smegmatis via cell wall modulation and altered immune response of host.

Host-pathogen relationship in retreated tuberculosis with major rifampicin resistance-conferring mutations

Introduction

It is assumed that host defense systems eliminating the pathogen and regulating tissue damage make a strong impact on the outcome of tuberculosis (TB) disease and that these processes are affected by rifampicin (RIF) resistance-conferring mutations of Mycobacterium tuberculosis (Mtb). However, the host responses to the pathogen harboring different mutations have not been studied comprehensively in clinical settings. We analyzed clinico-epidemiological factors and blood transcriptomic signatures associated with major rpoB mutations conferring RIF resistance in a cohort study.

Methods

Demographic data were collected from 295 active pulmonary TB patients with treatment history in Hanoi, Vietnam. When recruited, drug resistance-conferring mutations and lineage-specific variations were identified using whole-genome sequencing of clinical Mtb isolates. Before starting retreatment, total RNA was extracted from the whole blood of HIV-negative patients infected with Mtb that carried either the rpoB H445Y or rpoB S450L mutation, and the total RNA was subjected to RNA sequencing after age-gender matching. The individual RNA expression levels in the blood sample set were also measured using real-time RT-PCR. Logistic and linear regression models were used to assess possible associations.

Results

In our cohort, rpoB S450L and rpoB H445Y were major RIF resistance-conferring mutations [32/87 (36.8%) and 15/87 (17.2%), respectively]. H445Y was enriched in the ancient Beijing genotype and was associated with nonsynonymous mutations of Rv1830 that has been reported to regulate antibiotic resilience. H445Y was also more frequently observed in genetically clustered strains and in samples from patients who had received more than one TB treatment episode. According to the RNA sequencing, gene sets involved in the interferon-γ and-α pathways were downregulated in H445Y compared with S450L. The qRT-PCR analysis also confirmed the low expression levels of interferon-inducible genes, including BATF2 and SERPING1, in the H445Y group, particularly in patients with extensive lesions on chest X-ray.

Discussion

Our study results showed that rpoB mutations as well as Mtb sublineage with additional genetic variants may have significant effects on host response. These findings strengthen the rationale for investigation of host-pathogen interactions to develop countermeasures against epidemics of drug-resistant TB.

Understanding drivers of phylogenetic clustering and terminal branch lengths distribution in epidemics of Mycobacterium tuberculosis

Detecting factors associated with transmission is important to understand disease epidemics, and to design effective public health measures. Clustering and terminal branch lengths (TBL) analyses are commonly applied to genomic data sets of Mycobacterium tuberculosis (MTB) to identify sub-populations with increased transmission. Here, I used a simulation-based approach to investigate what epidemiological processes influence the results of clustering and TBL analyses, and whether differences in transmission can be detected with these methods. I simulated MTB epidemics with different dynamics (latency, infectious period, transmission rate, basic reproductive number R0, sampling proportion, sampling period, and molecular clock), and found that all considered factors, except for the length of the infectious period, affect the results of clustering and TBL distributions. I show that standard interpretations of this type of analyses ignore two main caveats: (1) clustering results and TBL depend on many factors that have nothing to do with transmission, (2) clustering results and TBL do not tell anything about whether the epidemic is stable, growing, or shrinking, unless all the additional parameters that influence these metrics are known, or assumed identical between sub-populations. An important consequence is that the optimal SNP threshold for clustering depends on the epidemiological conditions, and that sub-populations with different epidemiological characteristics should not be analyzed with the same threshold. Finally, these results suggest that different clustering rates and TBL distributions, that are found consistently between different MTB lineages, are probably due to intrinsic bacterial factors, and do not indicate necessarily differences in transmission or evolutionary success.

Whole genome sequencing of Mycobacterium tuberculosis for prediction of drug resistance

Whole-genome sequencing (WGS) has shown tremendous potential in rapid diagnosis of drug-resistant tuberculosis (TB). In the current study, we performed WGS on drug-resistant Mycobacterium tuberculosis isolates obtained from Shanghai (n = 137) and Russia (n = 78). We aimed to characterise the underlying and high-frequency novel drug-resistance-conferring mutations, and also create valuable combinations of resistance mutations with high predictive sensitivity to predict multidrug- and extensively drug-resistant tuberculosis (MDR/XDR-TB) phenotype using a bootstrap method. Most strains belonged to L2.2, L4.2, L4.4, L4.5 and L4.8 lineages. We found that WGS could predict 82.07% of phenotypically drug-resistant domestic strains. The prediction sensitivity for rifampicin (RIF), isoniazid (INH), ethambutol (EMB), streptomycin (STR), ofloxacin (OFL), amikacin (AMK) and capreomycin (CAP) was 79.71%, 86.30%, 76.47%, 88.37%, 83.33%, 70.00% and 70.00%, respectively. The mutation combination with the highest sensitivity for MDR prediction was rpoB S450L + rpoB H445A/P + katG S315T + inhA I21T + inhA S94A, with a sensitivity of 92.17% (0.8615, 0.9646), and the mutation combination with highest sensitivity for XDR prediction was rpoB S450L + katG S315T + gyrA D94G + rrs A1401G, with a sensitivity of 92.86% (0.8158, 0.9796). The molecular information presented here will be of particular value for the rapid clinical detection of MDR- and XDR-TB isolates through laboratory diagnosis.

Novel mutations detected from drug resistant Mycobacterium tuberculosis isolated from North East of Thailand

The emergence of drug-resistant tuberculosis is a major global public health threat. Thailand is one of the top 14 countries with high tuberculosis and multi-drug resistant tuberculosis rates. Immediate detection of drug-resistant tuberculosis is necessary to reduce mortality and morbidity by effectively providing treatment to ameliorate the formation of resistant strains. Limited data exist of mutation profiles in Northeastern Thailand. Here, 65 drug-resistant Mycobacterium tuberculosis isolates were used to detect mutations by polymerase chain reaction (PCR) and DNA sequencing. In the katG gene, mutations were occurred in 47 (79.7%) among 59 isoniazid resistant samples. For rpoB gene, 31 (96.9%) were observed as mutations in 32 rifampicin resistant isolates. Of 47 katG mutation samples, 45 (95.7%) had mutations in katG315 codon and 2 (4.3%) showed novel mutations at katG365 with amino acid substitution of CCG-CGG (Pro-Arg). Moreover, out of 31 rpoB mutation isolates, the codon positions rpoB516, rpoB526, rpoB531 and rpoB533 were 3 (9.7%), 8 (25.8%), 11 (35.5%) and 1 (3.2%), respectively. Seven isolates of double point mutation were found [rpoB516, 526; 1 (3.2%) and rpoB516, 531; 6 (19.4%)]. In addition, 1 (3.2%) sample had triple point mutation at codon positions rpoB516, 526 and 531. Common and novel mutation codons of the rpoB and katG genes were generated. Although DNA sequencing showed high accuracy, conventional PCR could be applied as an initial marker for screening drug-resistant Mycobacterium tuberculosis isolates in limit resources region. Mutations reported here should be considered when developing new molecular diagnostic methods for implementation in Northeastern Thailand.

Latent tuberculosis: interaction of virulence factors in Mycobacterium tuberculosis

Tuberculosis (TB) remains a prominent health concern worldwide. Besides extensive research and vaccinations available, attempts to control the pandemic are cumbersome due to the complex physiology of Mycobacterium tuberculosis (Mtb). Alongside the emergence of drug-resistant TB, latent TB has worsened the condition. The tubercle bacilli are unusually behaved and successful with its strategies to modulate genes to evade host immune system and persist within macrophages. Under latent/unfavorable conditions, Mtb conceals itself from immune system and modulates its genes. Among many intracellular modulated genes, important are those involved in cell entry, fatty acid degradation, mycolic acid synthesis, phagosome acidification inhibition, inhibition of phagosome-lysosome complex and chaperon protein modulation. Though the study on these genes date back to early times of TB, an insight on their inter-relation within and to newly evolved genes are still required. This review focuses on the findings and discussions on these genes, possible mechanism, credibility as target for novel drugs and repurposed drugs and their interaction that enables Mtb in survival, pathogenesis, resistance and latency.

Phenotypic and genotypic features of the Mycobacterium tuberculosis lineage 1 subgroup in central Vietnam

Mycobacterium tuberculosis (Mtb) has different features depending on different geographic areas. We collected Mtb strains from patients with smear-positive pulmonary tuberculosis in Da Nang, central Vietnam. Using a whole genome sequencing platform, including genome assembly complemented by long-read-sequencing data, genomic characteristics were studied. Of 181 Mtb isolates, predominant Vietnamese EAI4_VNM and EAI4-like spoligotypes (31.5%), ZERO strains (5.0%), and part of EAI5 (11.1%) were included in a lineage-1 (L1) sublineage, i.e., L1.1.1.1. These strains were found less often in younger people, and they genetically clustered less frequently than other modern strains. Patients infected with ZERO strains demonstrated less lung infiltration. A region in RD2bcg spanning six loci, i.e., PE_PGRS35, cfp21, Rv1985c, Rv1986, Rv1987, and erm(37), was deleted in EAI4_VNM, EAI4-like, and ZERO strains, whereas another 118 bp deletion in furA was specific only to ZERO strains. L1.1.1.1-sublineage-specific deletions in PE_PGRS4 and PE_PGRS22 were also identified. RD900, seen in ancestral lineages, was present in majority of the L1 members. All strains without IS6110 (5.0%) had the ZERO spoligo-pattern. Distinctive features of the ancestral L1 strains provide a basis for investigation of the modern versus ancestral Mtb lineages and allow consideration of countermeasures against this heterogeneous pathogen.

First description of the katG gene deletion in a Mycobacterium tuberculosis clinical isolate and its impact on the mycobacterial fitness

Isoniazid (INH) is the cornerstone of the anti-tuberculosis regimens and emergence of Mycobacterium tuberculosis (Mtb) resistant strains is a major threat to our ability to control tuberculosis (TB) at global level. Mutations in the gene coding the catalase KatG confer resistance to high level of INH. In this paper, we describe for the first time a complete deletion of the genomic region containing the katG gene in an Mtb clinical strain isolated in Italy in a patient with HIV infection that previously completed INH preventive therapy. We genotypically characterized the Mtb strain and showed that katG deletion confers high-level resistance to INH (MIC > 25.6 μg/mL). The katG deletion did not impact significantly on Mtb fitness as we did not detect enhanced susceptibility to H₂O₂ compared to the wild type Mtb strains nor impaired growth in in vitro infection models. These findings highlight the ability of Mtb to acquire resistance to INH while maintaining fitness and pathogenic potential.

Teleological cooption of Mycobacterium tuberculosis PE/PPE proteins as porins: Role in molecular immigration and emigration

Permeation through bacterial cells for exchange or uptake of biomolecules and ions invariably depend upon the existence of pore-forming proteins (porins) in their outer membrane. Mycobacterium tuberculosis (M. tb) harbours one of the most rigid cell envelopes across bacterial genera and is devoid of the classical porins for solute transport across the cell membrane. Though canonical porins are incompatible with the evolution of permeability barrier, porin like activity has been reported from membrane preparations of pathogenic mycobacteria. This suggests a sophisticated transport mechanism that has been elusive until now, along with the protein family responsible for it. Recent evidence suggests that these slow-growing mycobacteria have co-opted some of PE/PPE family proteins as molecular transport channels, in place of porins, to facilitate uptake of nutrients required to thrive in the restrictive host environment. These reports advocate that PE/PPE proteins, due to their structural ability, have a potential role in importing small molecules to the cell's interior. This mechanism unveils how a successful pathogen overcomes its restrictive membrane's transport limitations for selective uptake of nutrients. If extrapolated to have a role in drug transport, these channels could help understand the emergence of drug resistance. Further, as these proteins are associated with the export of virulence factors, they can be exploited as novel drug targets. There remains, however, an interesting question that as the PE/PPE proteins can allow the 'import' of molecules from outside the cell, is the reverse transport also possible across the M. tb membrane. In this review, we have discussed recent evidence supporting PE/PPE's role as a specific transport channel for selective uptake of small molecule nutrients and, as possible molecular export machinery of M. tb. This newly discovered role as transmembrane channels demands further research on this enigmatic family of proteins to comprehend the pathomechanism of this very smart pathogen.

Systematic Review of Mutations Associated with Isoniazid Resistance Points to Continuing Evolution and Subsequent Evasion of Molecular Detection, and Potential for Emergence of Multidrug Resistance in Clinical Strains of Mycobacterium tuberculosis

Molecular testing is rapidly becoming an integral component of global tuberculosis (TB) control. Uncommon mechanisms of resistance escape detection by these platforms and undermine our ability to contain outbreaks. This article is a systematic review of published articles that reported isoniazid (INH) resistance-conferring mutations between September 2013 and December 2019. The genes katG, inhA, and fabG1, and the intergenic region oxyR'-ahpC were considered in this review. Fifty-two articles were included that described 9,306 clinical isolates (5,804 INH resistant [INH^r] and 3,502 INH susceptible [INH^s]) from 31 countries. The three most frequently mutated loci continue to be locus 315 of katG (katG315; n = 4,271), locus -15 of inhA (inhA-15; n = 787), and locus -8 of inhA (inhA-8; 106). However, the diagnostic value of inhA-8 is far lower than previously thought, as it only appears in 25 (0.4%) of the INH^r isolates lacking the first two mutations. I catalogued 45 new loci (29 katG, nine inhA, and seven ahpC) associated with INH resistance and identified 59 loci (common to this and previous reviews) as a reliable basis for molecular diagnostics. Including all observed mutations provides a cumulative sensitivity of 85.6%. In 14.4% of resistant isolates, no mechanism of resistance was detected, making them likely to escape molecular detection, and in the case of INH monoresistance, likely to convert to multidrug-resistant TB (MDR-TB). Integrating the information cataloged in this study into current diagnostic tools is essential for combating the emergence of MDR-TB, and its exclusion can lead to an unintended selection against common mechanisms and to diversifying evolution. Observation of many low-frequency resistance-conferring mutations points to an advantage of whole-genome sequencing (WGS) for diagnostics. Finally, I provide five recommendations for future diagnostic platforms.

Higher genome mutation rates of Beijing lineage of Mycobacterium tuberculosis during human infection

Mycobacterium tuberculosis (Mtb) strains of Beijing lineage have caused great concern because of their rapid emergence of drug resistance and worldwide spread. DNA mutation rates that reflect evolutional adaptation to host responses and the appearance of drug resistance have not been elucidated in human-infected Beijing strains. We tracked and obtained an original Mtb isolate of Beijing lineage from the 1999 tuberculosis outbreak in Japan, as well as five other isolates that spread in humans, and two isolates from the patient caused recurrence. Three isolates were from patients who developed TB within one year after infection (rapid-progressor, RP), and the other three isolates were from those who developed TB more than one year after infection (slow-progressor, SP). We sequenced genomes of these isolates and analyzed the propensity and rate of genomic mutations. Generation time versus mutation rate curves were significantly higher for RP. The ratio of oxidative versus non-oxidation damages induced mutations was higher in SP than RP, suggesting that persistent Mtb are exposed to oxidative stress in the latent state. Our data thus demonstrates that higher mutation rates of Mtb Beijing strains during human infection is likely to account for the higher adaptability and an emergence ratio of drug resistance.

Chemical reactivity and bioactivity properties of pyrazinamide analogs of acetylsalicylic acid and salicylic acid using conceptual density functional theory

Conventional drugs used to treat Tuberculosis (TB) are becoming ineffective due to the occurrence of multiple drug resistant strains of tuberculosis (TB). This has made the TB disease a a serious global health dilemma. Hence, there is desperate necessity for the advancement of new drugs. In this work, the chemical reactivity and bioactivity of several analogs ofpyrazinamide (PZA) were investigated. PZA is one of the first-line of drugs used to treat tuberculosis and is a key contributor to shortening the treatment time for the disease. Chemical reactivity descriptors of pyrazinamide (PZA) and its analogs of acetylsalicyclic acid and salicyclic acid were investigated using conceptual density functional theory in water as a solvent at the MN12SX/Def2TZVP level of theory. Results have shown that all PZA analogs have improved their global and local reactivity indeces as compared to pyrazinamide based on its electronegativity, electrodonating power, electroaccepting power, eletrophilicity, global hardness and dual descriptor condensed fukui indexes. Moreover, their pKa values are slightly higher than PZA. In terms of its drug-likeness, all PZA analogs passed the Lipinski's Rule of Five criteria. Furthermore, their bioactivity scores are significantly better than pyrazinamide indicating good reaction to G-Protein Coupled Receptor (GPCR) ligands, kinase inhibitors, ion channel modulators, nuclear receptor ligands, protease inhibitors and other enzyme targets. Overall, the PZA analogs are found to be promising anti-tuberculosis drugs. Based on global and local reactivity descriptors, pKa and bioactivity scores, PZA analog of 5-n-Octanoylsalicylic acid is the most reactive among the PZA analogs tested.

Genomic signatures of drug resistance in highly resistant Mycobacterium tuberculosis strains of the early ancient sublineage of Beijing genotype in Russia

The Mycobacterium tuberculosis Beijing genotype is a clinically and epidemiologically important lineage that is subdivided into ancient/ancestral and modern strains. In our previous study in western Siberia, we identified variable number of tandem repeats (VNTR)-based clusters within the early ancient sublineage of the Beijing genotype characterized by an unexpectedly high rate of extensive drug resistance (XDR). In the current study, next generation sequencing data were analysed to gain insight into genomic signatures underlying drug resistance of these strains. A total of 184 genomes of the Beijing early ancient sublineage from Russia (16), China (15), Japan (36), Korea (25), Vietnam (18), Thailand (73), and the USA (1) were used for phylogenetic analysis. The drug-resistant profile was deduced genotypically. The Russian isolates were distributed into two clusters and were all drug resistant, mainly pre-XDR and XDR. The largest of these clusters included only Russian isolates from remote locations in both Asian and European parts of the country. All its isolates had a quadruple drug resistance (to isoniazid, rifampin, ethambutol and streptomycin) due to the 6-mutation signature (KatG Ser315Thr, KatG Ile335Val, RpoB Ser450Leu, RpoC Asp485Asn, EmbB Gln497Arg, and RpsL Lys43Arg). In most samples, it was complemented with additional and different pncA, gyrA and rrs mutations leading to the pre-XDR/XDR genotype. Phylogenomic analysis indicates a distant origin of this Russian resistant cluster in the early 1970s but location and circumstances are yet to be clarified.

Association Between the Phenotype and Genotype of Isoniazid Resistance Among Mycobacterium tuberculosis Isolates in Thailand

Purpose

The emergence of isoniazid-resistant tuberculosis (HR-TB) is a global public health problem, causing treatment failure and high mortality rates. This study aimed to determine the minimal inhibitory concentration (MIC) of isoniazid and detect the gene mutation in HR-TB and any association between the level of isoniazid resistance and gene mutation.

Methods

We collected 74 clinical HR-TB isolates from two tertiary-care centers in Thailand. MICs were established using broth macrodilution. A line probe assay (LPA) was used to detect gene mutations that confer resistance to isoniazid, rifampicin, aminoglycosides, and fluoroquinolones.

Results

Sixty-one (82.4%) isolates were monoresistant to isoniazid and 44 (72.1%) were highly resistant to isoniazid. From the clinical isolates, the range of isoniazid MICs was 0.4–16 μg/mL. The katG S315T gene mutation was the prominent mutation in both isoniazid-monoresistant TB (70.5%) and multidrug-resistant TB (72.7%) isolates. The positive predictive value (PPV) of katG was 100% in detecting high levels of isoniazid resistance. The PPV of the inhA mutation was 93.8% in detecting low levels of isoniazid resistance. Five isolates (6.8%) exhibited low-level phenotypic resistance, whereas an LPA failed to detect an isoniazid gene mutation. Our study found one HR-TB isolate with a gyrA fluoroquinolone-resistant gene mutation.

Conclusion

Most HR-TB isolates had high isoniazid-resistance levels associated with the katG gene mutation. High-dose isoniazid should be used with caution in patients with HR-TB. Early detection of drug resistance by genotypic assay can help determine an appropriate regimen.