Patterns of rpoC Mutations in Drug-Resistant Mycobacterium tuberculosis Isolated from Patients in South Korea

Large-scale statistical analysis of Mycobacterium tuberculosis genome sequences identifies compensatory mutations associated with multi-drug resistance

Tuberculosis (TB), caused by Mycobacterium tuberculosis, has a significant impact on global health worldwide. The development of multi-drug resistant strains that are resistant to the first-line drugs isoniazid and rifampicin threatens public health security. Rifampicin and isoniazid resistance are largely underpinned by mutations in rpoB and katG respectively and are associated with fitness costs. Compensatory mutations are considered to alleviate these fitness costs and have been observed in rpoC/rpoA (rifampicin) and oxyR'-ahpC (isoniazid). We developed a framework (CompMut-TB) to detect compensatory mutations from whole genome sequences from a large dataset comprised of 18,396 M. tuberculosis samples. We performed association analysis (Fisher's exact tests) to identify pairs of mutations that are associated with drug-resistance, followed by mediation analysis to identify complementary or full mediators of drug-resistance. The analyses revealed several potential mutations in rpoC (N = 47), rpoA (N = 4), and oxyR'-ahpC (N = 7) that were considered either 'highly likely' or 'likely' to confer compensatory effects on drug-resistance, including mutations that have previously been reported and validated. Overall, we have developed the CompMut-TB framework which can assist with identifying compensatory mutations which is important for more precise genome-based profiling of drug-resistant TB strains and to further understanding of the evolutionary mechanisms that underpin drug-resistance.

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.

Mutant RNA polymerase can reduce susceptibility to antibiotics via ppGpp-independent induction of a stringent-like response

Background

Mutations in RNA polymerase (RNAP) can reduce susceptibility to ciprofloxacin in Escherichia coli, but the mechanism of transcriptional reprogramming responsible is unknown. Strains carrying ciprofloxacin-resistant (Cip^R) rpoB mutations have reduced growth fitness and their impact on clinical resistance development is unclear.

Objectives

To assess the potential for Cip^RrpoB mutations to contribute to resistance development by estimating the number of distinct alleles. To identify fitness-compensatory mutations that ameliorate the fitness costs of Cip^RrpoB mutations. To understand how Cip^RrpoB mutations reprogramme RNAP.

Methods

E. coli strains carrying five different Cip^RrpoB alleles were evolved with selection for improved fitness and characterized for acquired mutations, relative fitness and MIC_Cip. The effects of dksA mutations and a ppGpp⁰ background on growth and susceptibility phenotypes associated with Cip^RrpoB alleles were determined.

Results

The number of distinct Cip^RrpoB mutations was estimated to be >100. Mutations in RNAP genes and in dksA can compensate for the fitness cost of Cip^RrpoB mutations. Deletion of dksA reduced the MIC_Cip for strains carrying Cip^RrpoB alleles. A ppGpp⁰ phenotype had no effect on drug susceptibility.

Conclusions

Cip^RrpoB mutations induce an ppGpp-independent stringent-like response. Approximately half of the reduction in ciprofloxacin susceptibility is caused by an increased affinity of RNAP to DksA while the other half is independent of DksA. Stringent-like response activating mutations might be the most diverse class of mutations reducing susceptibility to antibiotics.

Highly transmitted M. tuberculosis strains are more likely to evolve MDR/XDR and cause outbreaks, but what makes them highly transmitted?

Multi-Drug-Resistant strains of Mycobacterium tuberculosis (MDR-TB) are a serious obstacle to global TB eradication. While most MDR-TB strains are infrequently transmitted, a few cause large transmission clusters that contribute substantially to local MDR-TB burdens. Here we examine whether the known mutations in these strains can explain their success. Drug resistance mutations differ in fitness costs and strains can also acquire compensatory mutations (CM) to restore fitness, but some highly transmitted MDR strains have no CM. The acquisition of resistance mutations that maintain high transmissibility seems to occur by chance and are more likely in strains that are intrinsically highly transmitted and cause many cases. Modern Beijing lineage strains have caused several large outbreaks, but MDR outbreaks are also caused by ancient Beijing and lineage 4 strains, suggesting the lineage is less important than the characteristics of the individual strain. The development of fluoroquinolone resistance appears to represent another level of selection, in which strains must surmount unknown fitness costs of gyrA mutations. The genetic determinants of high transmission are poorly defined but may involve genes encoding proteins involved in molybdenum acquisition and the Esx systems. In addition, strains eliciting lower cytokine responses and producing more caseating granulomas may have advantages for transmission. Successful MDR/XDR strains generally evolve from highly transmitted drug sensitive parent strains due to selection pressures from deficiencies in local TB control programs. Until TB incidence is considerably reduced, there will likely be highly transmitted strains that develop resistance to any new antibiotic.

Drug resistance, fitness and compensatory mutations in Mycobacterium tuberculosis

For tuberculosis to be eradicated, the transmission of Multi-Drug-Resistant and eXtensively Drug Resistant strains of Mycobacterium tuberculosis (MDR and XDR-TB) must be considerably reduced. Drug resistant strains were initially thought to have reduced fitness, and the majority of resistant strains may actually have compromised fitness because they are found in only one or a few patients. In contrast, some MDR/XDR-TB strains are highly transmitted and cause large outbreaks. Most antibiotics target essential bacterial functions and the mutations that confer resistance to anti-TB drugs can incur fitness costs manifested as slower growth and reduced viability. The fitness costs vary with different resistance mutations and the bacilli can also accumulate secondary mutations that compensate for the compromised functions and partially or fully restore lost fitness. The compensatory mutations (CM) are different for each antibiotic, as they mitigate the deleterious effects of the specific functions compromised by the resistance mutations. CM are generally more common in strains with resistance mutations incurring the greatest fitness costs, but for RIF resistance, CM are most frequent in strains with the mutation carrying the least fitness cost, Ser450Leu. Here, we review what is known about fitness costs, CM and mechanisms of resistance to the drugs that define a strain as MDR or XDR-TB. The relative fitness costs of the resistance mutations and the mitigating effects of CM largely explain why certain mutations are frequently found in highly transmitted clusters while others are less frequently, rarely or never found in clinical isolates. The CM illustrate how drug resistance affects bacteria and how bacteria evolve to overcome the effects of the antibiotics, and thus a paradigm for how mycobacteria can evolve in response to stress.

rpoB gene mutations in rifampin-resistant Mycobacterium tuberculosis isolates from rural areas of Zhejiang, China

Objective

The aim was to analyze genetic mutations in the rpoB gene of rifampin-resistant Mycobacterium tuberculosis isolates (RIF^R-MTB) from Zhejiang, China.

Methods

We prospectively analyzed RIF^R-associated mutations in 13 rural areas of Zhejiang. Isolates were subjected to species identification, phenotype drug susceptibility testing (DST), DNA extraction, and rpoB gene sequencing.

Results

A total of 103 RIF^R isolates were identified by DST (22 RIF^R only, 14 poly-drug resistant, 49 multidrug resistant, 13 pre-extensively drug resistant [pre-XDR], and 5 extensively drug resistant [XDR]) from 2152 culture-positive sputum specimens. Gene sequencing of rpoB showed that the most frequent mutation was S450L (37.86%, 39/103); mutations P280L, E521K, and D595Y were outside the rifampicin resistance-determining region (RRDR) but may be associated with RIF^R. Mutations associated with poly-drug resistant, pre-XDR, and XDR TB were mainly located at codon 445 or 450 in the RRDR.

Conclusions

The frequency of rpoB RRDR mutation in Zhejiang is high. Further studies are needed to clarify the relationships between RIF^R and the TTC insertion at codon 433 in the RRDR and the P280L and D595Y mutations outside the RRDR.

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.

Comparative genome analysis of multidrug-resistant and susceptible Glaesserella parasuis strains isolated from diseased pigs in China

Glaesserella parasuis (G. parasuis) is a respiratory pathogen of swine and the etiological agent of Glässer's disease. Although the emergence of multidrug-resistant (MDR) G. parasuis is a critical problem in the swine industry, there are few publications on the genetic basis of antimicrobial resistance of G. parasuis. In this study, comparative genome analyses were used to identify genomic differences between two phenotypically distinct isolates, an MDR isolate (HPS-1) and a susceptible isolate (HPS-2), from diseased swines in China. These isolates were both serovar 4, which is predominant in cases of Glässer's disease and is the most prevalent serovar in China. Based on clusters of orthologous group (COG) annotations, genes assigned to the extracellular structure category were only detected in HPS-1 and genes related to cell motility were more abundant in HPS-1 than in HPS-2. A comparative genomic analysis showed that these two isolates are closely related, although there was a large-scale genomic rearrangement. Eighteen percent of the genome consisted of strain-specific accessory genes of HPS-1. Notably, only the two genes aac(6')-Ie-aph(2'')-Ia and bla_ROB-1 on a plasmid were specific to HPS-1. We also detected 30,599 single nucleotide polymorphisms (SNPs), including nonsynonymous SNPs in the aminoglycoside resistance gene aph(3'')-Ib, the fusidic acid resistance gene fusA, and the two rifampicin resistance genes rpoC and rpoB in HPS-1. These findings improve our understanding of the differences between MDR and susceptible isolates and will aid the development of treatment strategies to decrease the prevalence and disease burden caused by G. parasuis.

Genotypic characterization of 'inferred' rifampin mutations in GenoType MTBDRplus assay and its association with phenotypic susceptibility testing of Mycobacterium tuberculosis

In GenoType MTBDRplus assay [line probe assay (LPA)], when Mycobacterium tuberculosis (M. tuberculosis) sample DNA fails to hybridize to at least 1 rpoB wild-type probe and any mutation probe, it is inferred as rifampin (RIF)-resistant. In this study, we sought to identify such 'inferred' mutations in M. tuberculosis isolates (n = 203) by rpoB gene sequencing and determined their association with phenotypic resistance. D516Y, H526N, L511P mutations were associated with both phenotypically sensitive (59%, n = 38/64) and resistant (23.7%, n = 33/139) antimicrobial susceptibility testing (AST) results, whereas S531W mutation was associated with only RIF-resistant isolates (33%, n = 46/139). These results demonstrated that, at standard drug concentrations, some 'inferred' mutations may be missed by RIF-AST (phenotypically sensitive). The use of LPA permits identification of these RIF-resistant isolates, and incorporation of additional mutation probes (e.g., S531W) could further increase LPA specificity. Further studies are needed to establish the significance of the type of 'inferred' mutation with clinical/treatment outcomes.

Characteristics of compensatory mutations in the rpoC gene and their association with compensated transmission of Mycobacterium tuberculosis

The aim of this study was to characterize rpoC gene mutations in Mycobacterium tuberculosis (MTB) and investigate the factors associated with rpoC mutations and the relation between rpoC mutations and tuberculosis (TB) transmission. A total of 245 MTB clinical isolates from patients with TB in six provinces and two municipalities in China were characterized based on gene mutations through DNA sequencing of rpoC and rpoB genes, phenotyping via standard drug susceptibility testing, and genotypic profiling by mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing. Approximately 36.4% of the rifampin-resistant isolates harbored nonsynonymous mutations in the rpoC gene. Twenty-nine nonsynonymous single mutations and three double mutations were identified. The rpoC mutations at locus 483 (11.3%) were predominant, and the mutations at V483G, W484G, I491V, L516P, L566R, N698K, and A788E accounted for 54.5% of the total detected mutations. Fifteen new mutations in the rpoC gene were identified. Rifampin resistance and rpoB mutations at locus 531 were significantly associated with rpoC mutations. MIRU-VNTR genotype results indicated that 18.4% of the studied isolates were clustered, and the rpoC mutations were not significantly associated with MIRU-VNTR clusters. A large proportion of rpoC mutation was observed in the rifampicin-resistant MTB isolates. However, the findings of this study do not support the association of rpoC mutation with compensated transmissibility.

Phylogenetic Analysis for the Origin of Typhoid Fever Outbreak on Jeju Island, Korea, in 2017

In June 2017, on Jeju Island, six patients were diagnosed with typhoid fever. All six patients were admitted to our hospital with fever, all of them had been working at the same hotel. The gyrA and gyrB genes of four out of six Salmonella typhi cases were amplified by polymerase chain reaction. The phylogenetic tree of the genes showed that the S. Typhi isolates from Jeju Island in 2017 had the same genotype and were similar to isolates from the United Kingdom. The origin of the isolates of the outbreak was the same, and the disease strain may have been imported from outside of Korea.

Identification and Characterization of Genetic Determinants of Isoniazid and Rifampicin Resistance in Mycobacterium tuberculosis in Southern India

Drug-resistant tuberculosis (TB), one of the leading causes of death worldwide, arises mainly from spontaneous mutations in the genome of Mycobacterium tuberculosis. There is an urgent need to understand the mechanisms by which the mutations confer resistance in order to identify new drug targets and to design new drugs. Previous studies have reported numerous mutations that confer resistance to anti-TB drugs, but there has been little systematic analysis to understand their genetic background and the potential impacts on the drug target stability and/or interactions. Here, we report the analysis of whole-genome sequence data for 98 clinical M. tuberculosis isolates from a city in southern India. The collection was screened for phenotypic resistance and sequenced to mine the genetic mutations conferring resistance to isoniazid and rifampicin. The most frequent mutation among isoniazid and rifampicin isolates was S315T in katG and S450L in rpoB respectively. The impacts of mutations on protein stability, protein-protein interactions and protein-ligand interactions were analysed using both statistical and machine-learning approaches. Drug-resistant mutations were predicted not only to target active sites in an orthosteric manner, but also to act through allosteric mechanisms arising from distant sites, sometimes at the protein-protein interface.