The thiol methyltransferase activity of TMT1A (METTL7A) is conserved across species.

Although few resistance mechanisms for histone deacetylase inhibitors (HDACis) have been described, we recently demonstrated that TMT1A (formerly METTL7A) and TMT1B (formerly METTL7B) can mediate resistance to HDACis with a thiol as the zinc-binding group by methylating and inactivating the drug. TMT1A and TMT1B are poorly characterized, and their normal physiological role has yet to be determined. As animal model systems are often used to determine the physiological function of proteins, we investigated whether the ability of these methyltransferases to methylate thiol-based HDACis is conserved across different species. We found that TMT1A was conserved across rats, mice, chickens, and zebrafish, displaying 85.7%, 84.8%, 60.7%, and 51.0% amino acid sequence identity, respectively, with human TMT1A. Because TMT1B was not found in the chicken or zebrafish, we focused our studies on the TMT1A homologs. HEK-293 cells were transfected to express mouse, rat, chicken, or zebrafish homologs of TMT1A and all conferred resistance to the thiol-based HDACIs NCH-51, KD-5170, and romidepsin compared to empty vector-transfected cells. Additionally, all homologs blunted the downstream effects of HDACi treatment such as increased p21 expression, increased acetylated histone H3, and cell cycle arrest. Increased levels of dimethylated romidepsin were also found in the culture medium of cells transfected to express any of the TMT1A homologs after a 24 h incubation with romidepsin compared to empty-vector transfected cells. Our results indicate that the ability of TMT1A to methylate molecules is conserved across species. Animal models may therefore be useful in elucidating the role of these enzymes in humans.

The modification landscape of P. aeruginosa tRNAs.

RNA modifications have a substantial impact on tRNA function, with modifications in the anticodon loop contributing to translational fidelity and modifications in the tRNA core impacting structural stability. In bacteria, tRNA modifications are crucial for responding to stress and regulating the expression of virulence factors. Although tRNA modifications are well-characterized in a few model organisms, our knowledge of tRNA modifications in human pathogens, such as Pseudomonas aeruginosa, remains limited. Here we leveraged two orthogonal approaches to build a reference landscape of tRNA modifications in E. coli, which enabled us to identify similar modifications in P. aeruginosa. Our analysis revealed a substantial degree of conservation between the two organisms, while also uncovering potential sites of tRNA modification in P. aeruginosa tRNAs that are not present in E. coli. The mutational signature at one of these sites, position 46 of tRNAGln1(UUG) is dependent on the P. aeruginosa homolog of TapT, the enzyme responsible for the 3-(3-amino-3-carboxypropyl) uridine (acp3U) modification. Identifying which modifications are present on different tRNAs will uncover the pathways impacted by the different tRNA modifying enzymes, some of which play roles in determining virulence and pathogenicity.

The modification landscape of P. aeruginosa tRNAs.

RNA modifications have a substantial impact on tRNA function, with modifications in the anticodon loop contributing to translational fidelity and modifications in the tRNA core impacting structural stability. In bacteria, tRNA modifications are crucial for responding to stress and regulating the expression of virulence factors. Although tRNA modifications are well-characterized in a few model organisms, our knowledge of tRNA modifications in human pathogens, such as Pseudomonas aeruginosa, remains limited. Here we leveraged two orthogonal approaches to build a reference landscape of tRNA modifications in E. coli, which enabled us to identify similar modifications in P. aeruginosa. Our analysis revealed a substantial degree of conservation between the two organisms, while also uncovering potential sites of tRNA modification in P. aeruginosa tRNAs that are not present in E. coli. The mutational signature at one of these sites, position 46 of tRNAGln1(UUG) is dependent on the P. aeruginosa homolog of TapT, the enzyme responsible for the 3-(3-amino-3-carboxypropyl) uridine (acp3U) modification. Identifying which modifications are present on different tRNAs will uncover the pathways impacted by the different tRNA modifying enzymes, some of which play roles in determining virulence and pathogenicity.

Rewiring of RNA methylation by the oncometabolite fumarate in renal cell carcinoma.

Metabolic reprogramming is a hallmark of cancer that facilitates changes in many adaptive biological processes. Mutations in the tricarboxylic acid cycle enzyme fumarate hydratase (FH) lead to fumarate accumulation and cause hereditary leiomyomatosis and renal cell cancer (HLRCC). HLRCC is a rare, inherited disease characterized by the development of non-cancerous smooth muscle tumors of the uterus and skin, and an increased risk of an aggressive form of kidney cancer. Fumarate has been shown to inhibit 2-oxoglutarate-dependent dioxygenases (2OGDDs) involved in the hydroxylation of HIF1α, as well as in DNA and histone demethylation. However, the link between fumarate accumulation and changes in RNA post-transcriptional modifications has not been defined. Here, we determine the consequences of fumarate accumulation on the activity of different members of the 2OGDD family targeting RNA modifications. By evaluating multiple RNA modifications in patient-derived HLRCC cell lines, we show that mutation of FH selectively affects the levels of N6-methyladenosine (m6A), while the levels of 5-formylcytosine (f5C) in mitochondrial tRNA are unaffected. This supports the hypothesis of a differential impact of fumarate accumulation on distinct RNA demethylases. The observation that metabolites modulate specific subsets of RNA-modifying enzymes offers new insights into the intersection between metabolism and the epitranscriptome.

The Methyltransferases METTL7A and METTL7B Confer Resistance to Thiol-Based Histone Deacetylase Inhibitors.

Histone deacetylase inhibitors (HDACi) are part of a growing class of epigenetic therapies used for the treatment of cancer. Although HDACis are effective in the treatment of T-cell lymphomas, treatment of solid tumors with this class of drugs has not been successful. Overexpression of the multidrug resistance protein P-glycoprotein (P-gp), encoded by ABCB1, is known to confer resistance to the HDACi romidepsin in vitro, yet increased ABCB1 expression has not been associated with resistance in patients, suggesting that other mechanisms of resistance arise in the clinic. To identify alternative mechanisms of resistance to romidepsin, we selected MCF-7 breast cancer cells with romidepsin in the presence of the P-gp inhibitor verapamil to reduce the likelihood of P-gp-mediated resistance. The resulting cell line, MCF-7 DpVp300, does not express P-gp and was found to be selectively resistant to romidepsin but not to other HDACis such as belinostat, panobinostat, or vorinostat. RNA-sequencing analysis revealed upregulation of the mRNA coding for the putative methyltransferase, METTL7A, whose paralog, METTL7B, was previously shown to methylate thiol groups on hydrogen sulfide and captopril. As romidepsin has a thiol as the zinc-binding moiety, we hypothesized that METTL7A could inactivate romidepsin and other thiol-based HDACis via methylation of the thiol group. We demonstrate that expression of METTL7A or METTL7B confers resistance to thiol-based HDACis and that both enzymes are capable of methylating thiol-containing HDACis. We thus propose that METTL7A and METTL7B confer resistance to thiol-based HDACis by methylating and inactivating the zinc-binding thiol.

The thiol methyltransferase activity of TMT1A (METTL7A) is conserved across species.

Although few resistance mechanisms for histone deacetylase inhibitors (HDACis) have been described, we recently demonstrated that TMT1A (formerly METTL7A) and TMT1B (formerly METTL7B) can mediate resistance to HDACis with a thiol as the zinc-binding group by methylating and inactivating the drug. TMT1A and TMT1B are poorly characterized, and their normal physiological role has yet to be determined. As animal model systems are often used to determine the physiological function of proteins, we investigated whether the ability of these methyltransferases to methylate thiol-based HDACis is conserved across different species. We found that TMT1A was conserved across rats, mice, chickens, and zebrafish, displaying 85.7%, 84.8%, 60.7% and 51.0% amino acid sequence identity, respectively, with human TMT1A. Because TMT1B was not found in the chicken or zebrafish, we focused our studies on the TMT1A homologs. HEK-293 cells were transfected to express mouse, rat, chicken, or zebrafish homologs of TMT1A and all conferred resistance to the thiol-based HDACIs NCH-51, KD-5170 and romidepsin compared to empty vector-transfected cells. Additionally, all homologs blunted the downstream effects of HDACi treatment such as increased p21 expression, increased acetylated histone H3, and cell cycle arrest. Increased levels of dimethylated romidepsin were also found in the culture medium of cells transfected to express any of the TMT1A homologs after a 24 h incubation with romidepsin compared to empty-vector transfected cells. Our results indicate that the ability of TMT1A to methylate molecules is conserved across species. Animal models may therefore be useful in elucidating the role of these enzymes in humans.

Extensible benchmarking of methods that identify and quantify polyadenylation sites from RNA-seq data.

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, limitations, and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for continuous extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies, while the containers and reproducible workflows could easily be deployed and extended to evaluate new methods or data sets.

Cryptic splice mutation in the fumarate hydratase gene in patients with clinical manifestations of Hereditary Leiomyomatosis and Renal Cell Cancer.

Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is an autosomal dominant condition characterized by the development of cutaneous and uterine leiomyomas and risk for development of an aggressive form of papillary renal cell cancer. HLRCC is caused by germline inactivating pathogenic variants in the fumarate hydratase (FH) gene, which encodes the enzyme that catalyzes the interconversion of fumarate and L-malate. We utilized enzyme and protein mobility assays to evaluate the FH enzyme in a cohort of patients who showed clinical manifestations of HLRCC but were negative for known pathogenic FH gene variants. FH enzyme activity and protein levels were decreased by 50% or greater in three family members, despite normal FH mRNA expression levels as measured by quantitative PCR. Direct Nanopore RNA sequencing demonstrated 57 base pairs of retained intron sequence between exons 9 and 10 of polyadenylated FH mRNA in these patients, resulting in a truncated FH protein. Genomic sequencing revealed a heterozygous intronic alteration of the FH gene (chr1: 241498239 T/C) resulting in formation of a splice acceptor site near a polypyrimidine tract, and a uterine fibroid obtained from a patient showed loss of heterozygosity at this site. The same intronic FH variant was identified in an unrelated patient who also showed a clinical phenotype of HLRCC. These data demonstrate that careful clinical assessment as well as biochemical characterization of FH enzyme activity, protein expression, direct RNA sequencing, and genomic DNA sequencing of patient-derived cells can identify pathogenic variants outside of the protein coding regions of the FH gene.

Extensible benchmarking of methods that identify and quantify polyadenylation sites from RNA-seq data.

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, and limitations and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for seamless extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies. Furthermore, the containers and reproducible workflows generated in the course of this project can be seamlessly deployed and extended in the future to evaluate new methods or datasets.

Dual Inhibition of Histone Deacetylases and the Mechanistic Target of Rapamycin Promotes Apoptosis in Cell Line Models of Uveal Melanoma.

Purpose: Over 90% of uveal melanomas harbor pathogenic variants of the GNAQ or GNA11 genes that activate survival pathways. As previous studies found that Ras-mutated cell lines were vulnerable to a combination of survival pathway inhibitors and the histone-deacetylase inhibitor romidepsin, we investigated whether this combination would be effective in models of uveal melanoma. Methods: A small-scale screen of inhibitors of bromodomain-containing protein 4 (BRD4; OTX-015), extracellular signal-related kinase (ERK; ulixertinib), mechanistic target of rapamycin (mTOR; AZD-8055), or phosphoinositide 3-kinase (PI3K; GDC-0941) combined with a clinically relevant administration of romidepsin was performed on a panel of uveal melanoma cell lines (92.1, Mel202, MP38, and MP41) and apoptosis was quantified by flow cytometry after 48 hours. RNA sequencing analysis was performed on Mel202 cells treated with romidepsin alone, AZD-8055 alone, or the combination, and protein changes were validated by immunoblot. Results: AZD-8055 with romidepsin was the most effective combination in inducing apoptosis in the cell lines. Increased caspase-3 and PARP cleavage were noted in the cell lines when they were treated with romidepsin and mTOR inhibitors. RNA sequencing analysis of Mel202 cells revealed that apoptosis was the most affected pathway in the romidepsin/AZD-8055-treated cells. Increases in pro-apoptotic BCL2L11 and decreases in anti-apoptotic BIRC5 and BCL2L1 transcripts noted in the sequencing analysis were confirmed at the protein level in Mel202 cells. Conclusions: Our data suggest that romidepsin in combination with mTOR inhibition could be an effective treatment strategy against uveal melanoma due in part to changes in apoptotic proteins.

Controlling CRISPR-Cas9 with ligand-activated and ligand-deactivated sgRNAs.

The CRISPR-Cas9 system provides the ability to edit, repress, activate, or mark any gene (or DNA element) by pairing of a programmable single guide RNA (sgRNA) with a complementary sequence on the DNA target. Here we present a new method for small-molecule control of CRISPR-Cas9 function through insertion of RNA aptamers into the sgRNA. We show that CRISPR-Cas9-based gene repression (CRISPRi) can be either activated or deactivated in a dose-dependent fashion over a >10-fold dynamic range in response to two different small-molecule ligands. Since our system acts directly on each target-specific sgRNA, it enables new applications that require differential and opposing temporal control of multiple genes.

It's complicated m6A-dependent regulation of gene expression in cancer.

Cellular function relies on multiple pathways that are coordinated to ensure the proper execution of gene expression networks. Failure to coordinate the multiple programs active in the cell can have catastrophic consequences and lead to diseases such as cancer. At the post-transcriptional level, RNA modifications play important roles in the regulation of gene expression. N6-methyladenosine (m6A) is the most abundant internal messenger RNA (mRNA) modification and has gained increasing interest in the last few years as a dynamic regulator of RNA metabolism. Modifications regulate all stages of the RNA life cycle, from transcription to decay. Recent studies have pointed to the role of RNA methylation in cancer initiation and progression, and aberrant modification has served as a biomarker of early-stage diagnosis in several cancers. Here, we review the regulation of m6A, disruptions to methylation-dependent pathways that influence carcinogenesis, and potential avenues for m6A-related therapeutic strategies.

Determinants of tRNA Recognition by the Radical SAM Enzyme RlmN.

RlmN, a bacterial radical SAM methylating enzyme, has the unusual ability to modify two distinct types of RNA: 23S rRNA and tRNA. In rRNA, RlmN installs a methyl group at the C2 position of A2503 of 23S rRNA, while in tRNA the modification occurs at nucleotide A37, immediately adjacent to the anticodon triplet. Intriguingly, only a subset of tRNAs that contain an adenosine at position 37 are substrates for RlmN, suggesting that the enzyme carefully probes the highly conserved tRNA fold and sequence features to identify its targets. Over the past several years, multiple studies have addressed rRNA modification by RlmN, while relatively few investigations have focused on the ability of this enzyme to modify tRNAs. In this study, we utilized in vitro transcribed tRNAs as model substrates to interrogate RNA recognition by RlmN. Using chimeras and point mutations, we probed how the structure and sequence of RNA influences methylation, identifying position 38 of tRNAs as a critical determinant of substrate recognition. We further demonstrate that, analogous to previous mechanistic studies with fragments of 23S rRNA, tRNA methylation requirements are consistent with radical SAM reactivity. Together, our findings provide detailed insight into tRNA recognition by a radical SAM methylating enzyme.