Differential DNA methylomes of clinical MDR, XDR and XXDR Mycobacterium tuberculosis isolates revealed by using single-molecule real-time sequencing

Roles of the m6A methyltransferases METTL3, METTL14, and WTAP in pulmonary tuberculosis

Objective

The aim of the current study was to investigate the contributing role of gene variation and transcription levels among the m6A methyltransferases METTL3, METTL14, and WTAP in pulmonary tuberculosis (PTB).

Methods

A case-control study including 461 PTB patients and 467 normal controls was designed for genotyping. Three SNPs in METTL3 (rs1061027, rs1139130, rs1061026), three SNPs in METTL14 (rs62328061, rs4834698, rs1064034), and two SNPs in WTAP (rs1853259, rs11752345) were genotyped via the SNPscan™ technique. METTL3, METTL14, and WTAP transcription levels were determined in 78 PTB patients and 86 controls via quantitative real-time reverse-transcription PCR.

Results

Frequencies of the METTL14 rs62328061 GG genotype, WTAP rs11752345 CT genotype, and T allele were significantly increased in PTB patients compared to controls. An increased risk of rs62328061 was detected in a recessive model, and a decreased risk of rs11752345 was detected in a dominant model in the PTB group. METTL3 gene variation was not associated with PTB risk. The METTL3 rs1139130 GG genotype was significantly increased with drug resistance, and the G allele was significantly decreased with drug-induced liver injury in PTB patients. A reduced frequency of the METTL14 rs62328061 G allele was associated with leukopenia, a reduced frequency of the WTAP rs11752345 T allele was associated with sputum smear positivity, and a higher frequency of the METTL14 rs4834698 TC genotype was evident in PTB patients with hypoproteinemia. Compared to controls, METTL3, METTL14, and WTAP transcription levels in PTB patients were significantly decreased, and the level of WTAP was increased in PTB patients with drug resistance. METTL3 level was negatively associated with erythrocyte sedimentation rate and aspartate aminotransferase, and METTL14 level was negatively correlated with alanine aminotransferase and aspartate aminotransferase.

Conclusion

METTL14 rs62328061 and WTAP rs11752345 variants were associated with the genetic background of PTB, and METTL3, METTL14, and WTAP levels were abnormally decreased, suggesting that these m6A methyltransferases may play important roles in PTB.

mbtD and celA1 association with ethambutol resistance in Mycobacterium tuberculosis: A multiomics analysis

Ethambutol (EMB) is a first-line antituberculosis drug currently being used clinically to treat tuberculosis. Mutations in the embCAB operon are responsible for EMB resistance. However, the discrepancies between genotypic and phenotypic EMB resistance have attracted much attention. We induced EMB resistance in Mycobacterium tuberculosis in vitro and used an integrated genome–methylome–transcriptome–proteome approach to study the microevolutionary mechanism of EMB resistance. We identified 509 aberrantly methylated genes (313 hypermethylated genes and 196 hypomethylated genes). Moreover, some hypermethylated and hypomethylated genes were identified using RNA-seq profiling. Correlation analysis revealed that the differential methylation of genes was negatively correlated with transcription levels in EMB-resistant strains. Additionally, two hypermethylated candidate genes (mbtD and celA1) were screened by iTRAQ-based quantitative proteomics analysis, verified by qPCR, and corresponded with DNA methylation differences. This is the first report that identifies EMB resistance-related genes in laboratory-induced mono-EMB-resistant M. tuberculosis using multi-omics profiling. Understanding the epigenetic features associated with EMB resistance may provide new insights into the underlying molecular mechanisms.

Diversity and novel mutations in membrane transporters of Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB), the causative agent of tuberculosis (TB), encodes a family of membrane proteins belonging to Resistance-Nodulation-Cell Division (RND) permeases also called multidrug resistance pumps. Mycobacterial membrane protein Large (MmpL) transporters represent a subclass of RND transporters known to participate in exporting of lipid components across the cell envelope. These proteins perform an essential role in MTB survival; however, there are no data regarding mutations in MmpL, polyketide synthase (PKS) and acyl-CoA dehydrogenase FadE proteins from Khyber Pakhtunkhwa, Pakistan. This study aimed to screen mutations in transmembrane transporter proteins including MmpL, PKS and Fad through whole-genome sequencing (WGS) in local isolates of Khyber Pakhtunkhwa province, Pakistan. Fourteen samples were collected from TB patients and drug susceptibility testing was performed. However, only three samples were completely sequenced. Moreover, 209 whole-genome sequences of the same geography were also retrieved from NCBI GenBank to analyze the diversity of mutations in MmpL, PKS and Fad proteins. Among the 212 WGS (Accession ID: PRJNA629298, PRJNA629388, and ERR2510337-ERR2510345, ERR2510546-ERR2510645), numerous mutations in Fad (n = 756), PKS (n = 479), and MmpL (n = 306) have been detected. Some novel mutations were also detected in MmpL, PKS and acyl-CoA dehydrogenase Fad. Novel mutations including Asn576Ser in MmpL8, Val943Gly in MmpL9 and Asn145Asp have been detected in MmpL3. The presence of a large number of mutations in the MTB membrane may have functional consequences on proteins. However, further experimental studies are needed to elucidate the variants' effect on MmpL, PKS and FadE functions.

RNA Expression Analysis of Mycobacterial Methyltransferases Genes in Different Resistant Strains of Mycobacterium tuberculosis

Background

Tuberculosis infection still represents a global health issue affecting patients worldwide. Strategies for its control may be not as effective as it should be, specifically in case of resistant strains of Mycobacterium tuberculosis (M.tb.) In this regard, the role of mycobacterial methyltransferases (MTases) in TB infection can be fundamental, though it has not been broadly deciphered.

Methods

Five resistant isolates of M.tb were obtained. M.tb H37Rv (ATCC 27249) was used as a reference strain. Seven putative mycobacterial MTase genes (Rv0645c, Rv2966c, Rv1988, Rv1694, Rv3919c, Rv2756c, and Rv3263) and Rv1392 as SAM synthase were selected for analysis. PCR-sequencing and qRT-PCR were performed to compare mutations and expression levels of MTases in different strains. The 2-ΔΔCt method was employed to calculate the relative expression levels of these genes.

Results

Only two mutations were found in isoniazid resistance (INHR) strain for Rv3919c (T to G in codon 341) and Rv1392 (G to A in codon 97) genes. Overexpression of Rv0645c, Rv2756c, Rv3263, and Rv2966c was detected in all sensitive and resistant isolates. However, Rv1988 and Rv3919c decreased and Rv1694 increased in the sensitive strains. The Rv1392 expression level also decreased in INHR isolate.

Conclusion

We found a correlation between mycobacterial MTases expression and resistance to antibiotics in M.tb strains. Some MTases undeniably are virulence factors that specifically hijack the host defense mechanism. Further evaluations are needed to explore the complete impact of mycobacterial MTases within specific strains of M.tb to introduce novel diagnosis and treatment strategies.