Aberrant translation regulated by METTL1/WDR4-mediated tRNA N7-methylguanosine modification drives head and neck squamous cell carcinoma progression

Insights into the role of the N6-methyladenosine reader IGF2BP3 in the progression of oral squamous cell carcinoma and its connection to cell-cycle control

The genome of oral squamous cell carcinoma (OSCC) has been extensively characterized via bulk sequencing, revealing a multitude of genetic changes. The gene IGF2BP3, which encodes for the insulin-like growth factor 2 mRNA-binding protein 3, has been observed to be highly expressed in several types of cancer. This finding suggests that IGF2BP3 may play a significant role in the initiation and advancement of cancer. Nevertheless, the mechanisms by which IGF2BP3 contribute to OSCC are yet to be fully understood. In this study, we have observed that IGF2BP3 exhibits overexpression in OSCC. Based on our findings from bulk sequencing analysis, we have concluded that IGF2BP3 could potentially serve as a biomarker for predicting poor prognosis in OSCC. Moreover, it has been demonstrated that IGF2BP3 exhibits a significant association with the initiation and advancement of tumors both in vivo and in vitro. The evaluation of IGF2BP3 expression levels in relation to the cell cycle stage was conducted using single-cell RNA sequencing data. Tumor cells characterized by elevated IGF2BP3 expression demonstrated a higher percentage of cells in the G2/M transition phase. This study presents new findings indicating that the molecular target IGF2BP3 can serve as a prognostic indicator for tumors and has an impact on the development and progression of OSCC by influencing the regulation of the cell cycle.

WDR4 promotes HCC pathogenesis through N7-methylguanosine by regulating and interacting with METTL1

BACKGROUND

The N7-methylguanosine (m7G), a modification at defined internal positions within tRNAs and rRNAs, is correlated with tumor progression. Methyltransferase like 1 (METTL1)/ WD repeat domain 4 (WDR4) mediated tRNA m7G modification, which could alter many oncogenic mRNAs translation to promote progress of multiple cancer types. However, whether and how the internal mRNA m7G modification is involved in tumorigenesis remains unclear.

METHODS

The immunohistochemistry assay was conducted to detect the expression of WDR4 and METTL1 in hepatocellular carcinoma (HCC) and the expression of both genes whether contributes to the prognosis of the survival rate of HCC patients. Then, CCK8, colony formation assays and tumor xenograft models were conducted to determine the effects of WDR4 on HCC cells in vitro and vivo. Besides, dot blot assay, m7G-MeRIP-seq and RNA-seq analysis were conducted to determine whether WDR4 contributes to m7G modification and underlying mechanism in HCC cells. Finally, rescue and CO-IP assay were conducted to explore whether WDR4 and METTL1 proteins form a complex in Huh7 cells.

RESULTS

WDR4 modulates m7G modification at the internal sites of tumor-promoting mRNAs by forming the WDR4-METTL1 complex. WDR4 knockdown downregulated the expression of mRNA and protein levels of METTL1 gene and thus further modulate the formation of WDR4-METTL1 complex indirectly. METTL1 expression was markedly correlated with WDR4 expression in HCC tissues. HCC patients with high expression of both genes had a poor prognosis.

CONCLUSIONS

WDR4 may contribute to HCC pathogenesis by interacting with and regulating the expression of METTL1 to synergistically modulate the m7G modification of target mRNAs in tumor cells.

METTL1 facilitates ameloblastoma invasive growth via MAPK signaling pathway

OBJECTIVES

Ameloblastoma (AM), a common odontogenic epithelial tumor, exhibits aggressive growth due to incomplete encapsulation within the jawbone. Postoperative recurrence is a significant concern, closely associated with its invasive nature. We investigate the role of tRNA N-7 methylguanosine (m7G) modification mediated by Methyltransferase-like 1 (METTL1) in AM's invasive growth and prognosis.

MATERIALS AND METHODS

METTL1 expression was analyzed in diverse cell lines and clinical AM tissues. Its association with postoperative AM recurrence was examined. Functional experiments included METTL1 gene silencing using shRNA in hTERT-AM cells, assessing cell proliferation, migration, and invasion. Xenograft tumor model was constructed to investigate tumor growth. Molecular mechanisms behind METTL1's role in AM invasiveness were elucidated using Ribosome nascent-chain complex-bound mRNA sequencing (RNC-seq) and experimental analysis.

RESULTS

High METTL1 expression was significantly associated with postoperative recurrence in AM. The inhibition of AM development following METTL1 knockdown has been corroborated by experiments conducted both in vitro and in vivo. Analysis of RNC-seq data revealed that downregulated genes were predominantly enriched in the mitogen-activated protein kinase (MAPK) signaling pathway, suggesting that METTL1 may promote AM's invasive growth through the MAPK signaling pathway.

CONCLUSION

Our study elucidates the functional role of METTL1 in AM's invasive development and prognosis. High METTL1 expression is linked to postoperative recurrence, and METTL1 appears to promote AM invasiveness through the MAPK signaling pathway. These findings contribute to a better understanding of AM pathogenesis and may guide future therapeutic strategies.

Identification of N7-methylguanosine-related miRNAs as potential biomarkers for prognosis and drug response in breast cancer

Objectives

The impact of N7-methylguanosine (m7G) on tumor progression and the regulatory role of microRNAs (miRNAs) in immune function significantly influence breast cancer (BC) prognosis. Investigating the interplay between m7G modification and miRNAs provides novel insights for assessing prognostics and drug responses in BC.

Materials and methods

RNA sequences (miRNA and mRNA profiles) and clinical data for BC were acquired from the Cancer Genome Atlas (TCGA) database. A miRNA signature associated with 15 m7G in this cohort was identified using Cox regression and LASSO. The risk score model was evaluated using Kaplan-Meier and time-dependent ROC analysis, categorizing patients into high-risk and low-risk groups. Functional enrichment analyses were conducted to explore potential pathways. The immune system, including scores, cell infiltration, function, and drug sensitivity, was examined and compared between high-risk and low-risk groups. A nomogram that combines risk scores and clinical factors was developed and validated. Single-sample gene set enrichment analysis (ssGSEA) was employed to explore m7G-related miRNA signatures and immune cell relationships in the tumor microenvironment. Additionally, drug susceptibility was compared between risk groups.

Results

Fifteen m7G-related miRNAs were independently correlated with overall survival (OS) in BC patients. Time-dependent ROC analysis yielded area under the curve (AUC) values of 0.742, 0.726, and 0.712 for predicting 3-, 5-, and 10-year survival rates, respectively. The Kaplan-Meier analysis revealed a significant disparity in OS between the high-risk and low-risk groups (p = 1.3e-6). Multiple regression identified the risk score as a significant independent prognostic factor. An excellent calibration nomogram with a C-index of 0.785 (95 % CI: 0.728-0.843) was constructed. In immune analysis, low-risk patients exhibited heightened immune function and increased responsiveness to immunotherapy and chemotherapy compared to high-risk patients.

Conclusion

This study systematically analyzed m7G-related miRNAs and revealed their regulatory mechanisms concerning the tumor microenvironment (TME), pathology, and the prognosis of BC patient. Based on these miRNAs, a prognostic model and nomogram were developed for BC patients, facilitating prognostic assessments. These findings can also assist in predicting treatment responses and guiding medication selection.

Mitochondrial outer membrane protein MTUS1/ATIP1 exerts antitumor effects through ROS-induced mitochondrial pyroptosis in head and neck squamous cell carcinoma

We showed that microtubule-associated tumor suppressor gene (MTUS1/ATIP) downregulation correlated with poor survival in head and neck squamous cell carcinoma (HNSCC) patients and that MTUS1/ATIP1 was the most abundant isoform in HNSCC tissue. However, the location and function of MTUS1/ATIP1 have remain unclear. In this study, we confirmed that MTUS1/ATIP1 inhibited proliferation, growth and metastasis in HNSCC in cell- and patient-derived xenograft models in vitro and in vivo. MTUS1/ATIP1 localized in the outer mitochondrial membrane, influence the morphology, movement and metabolism of mitochondria and stimulated oxidative stress in HNSCC cells by directly interacting with MFN2. MTUS1/ATIP1 activated ROS, recruiting Bax to mitochondria, facilitating cytochrome c release to the cytosol to activate caspase-3, and inducing GSDME-dependent pyroptotic death in HNSCC cells. Our findings showed that MTUS1/ATIP1 localized in the outer mitochondrial membrane in HNSCC cells and mediated anticancer effects through ROS-induced pyroptosis, which may provide a novel therapeutic strategy for HNSCC treatment.

SOX18 meditates the resistance of Bmi1-expressing cells to cetuximab in HNSCC

OBJECTIVE

Head and neck squamous cell carcinoma (HNSCC) is the most common type of malignancy in the head and neck region worldwide. The therapeutic strategies for HNSCC remain unsatisfying and limited. Here, we found a population of resistant Bmi1-expressing cells in the presence of cetuximab treatment and reported a novel role of SRY-box transcription factor 18 (SOX18), a member of the SOX family, in promoting HNSCC resistance to cetuximab. This study aimed to investigate the regulatory mechanism of Sox18 in Bmi1-positive cells and to search for better therapeutic targets.

METHODS

We successfully obtained Bmi1CreER, RosatdTomato, and RosaDTA mice and identified Bmi1-expressing cells through lineage tracing. SOX18 expression in HNSCC and normal tissues was analyzed by immunohistochemistry, colocalization of Sox18, and Bmi1-expressing cells was analyzed by immunofluorescence, and SOX18 expression in SCC9 cell lines was quantified by western blotting and quantitative real-time PCR. The investigation of the mechanism of SOX18-mediated cetuximab resistance in Bmi1-positive cells was based on the analysis of single-cell RNA-seq data obtained from the Gene Expression Omnibus (GEO) database. Western blotting was performed to verify the results obtained from the single-cell RNA-seq analysis.

RESULTS

In our study, we demonstrated that Bmi1-expressing cells were resistant to cetuximab treatment and that depletion of Bmi1-expressing cells improved cetuximab efficacy in HNSCC. We then discovered that Sox18 mediated the stem cell-like properties of Bmi1-expressing cells and promoted cellular cetuximab resistance through an oxidative phosphorylation pathway. There was a significant downregulation of key genes in the oxidative phosphorylation pathway in Sox18 knockout cell lines.

CONCLUSIONS

Taken together, the findings of our study suggest that Sox18 mediates the resistance of Bmi1-expressing cells to cetuximab in HNSCC via the oxidative phosphorylation pathway.

The impact of tRNA modifications on translation in cancer: identifying novel therapeutic avenues

Recent advancements have illuminated the critical role of RNA modifications in post-transcriptional regulation, shaping the landscape of gene expression. This review explores how tRNA modifications emerge as critical players, fine-tuning functionalities that not only maintain the fidelity of protein synthesis but also dictate gene expression and translation profiles. Highlighting their dysregulation as a common denominator in various cancers, we systematically investigate the intersection of both cytosolic and mitochondrial tRNA modifications with cancer biology. These modifications impact key processes such as cell proliferation, tumorigenesis, migration, metastasis, bioenergetics and the modulation of the tumor immune microenvironment. The recurrence of altered tRNA modification patterns across different cancer types underscores their significance in cancer development, proposing them as potential biomarkers and as actionable targets to disrupt tumorigenic processes, offering new avenues for precision medicine in the battle against cancer.

Metabolic reprogramming driven by METTL1-mediated tRNA m7G modification promotes acquired anlotinib resistance in oral squamous cell carcinoma

Tyrosine kinase inhibitors (TKIs) are frequently utilized in the management of malignant tumors. Studies have indicated that anlotinib has a significant inhibitory effect on oral squamous cell carcinoma (OSCC). However, the mechanisms underlying the development of resistance with long-term anlotinib treatment remain obscure. Our research found that METTL1 expression was heightened in anlotinib-resistant OSCC cells. We observed that METTL1 played a role in fostering resistance to anlotinib in both transgenic mouse models and in vitro. Mechanistically, the elevated METTL1 levels in anlotinib-resistant OSCC cells contributed to enhanced global mRNA translation and stimulated oxidative phosphorylation (OXPHOS) through m7G tRNA modification. Bioenergetic profiling demonstrated that METTL1 drived a metabolic shift from glycolysis to OXPHOS in anlotinib-resistant OSCC cells. Additionally, inhibition of OXPHOS biochemically negated METTL1's impact on anlotinib resistance. Overall, this study underscores the pivotal role of METTL1-mediated m7G tRNA modification in anlotinib resistance and lays the groundwork for novel therapeutic interventions to counteract resistance in OSCC.

RNPS1 stabilizes NAT10 protein to facilitate translation in cancer via tRNA ac4C modification

Existing studies have underscored the pivotal role of N-acetyltransferase 10 (NAT10) in various cancers. However, the outcomes of protein-protein interactions between NAT10 and its protein partners in head and neck squamous cell carcinoma (HNSCC) remain unexplored. In this study, we identified a significant upregulation of RNA-binding protein with serine-rich domain 1 (RNPS1) in HNSCC, where RNPS1 inhibits the ubiquitination degradation of NAT10 by E3 ubiquitin ligase, zinc finger SWIM domain-containing protein 6 (ZSWIM6), through direct protein interaction, thereby promoting high NAT10 expression in HNSCC. This upregulated NAT10 stability mediates the enhancement of specific tRNA ac4C modifications, subsequently boosting the translation process of genes involved in pathways such as IL-6 signaling, IL-8 signaling, and PTEN signaling that play roles in regulating HNSCC malignant progression, ultimately influencing the survival and prognosis of HNSCC patients. Additionally, we pioneered the development of TRMC-seq, leading to the discovery of novel tRNA-ac4C modification sites, thereby providing a potent sequencing tool for tRNA-ac4C research. Our findings expand the repertoire of tRNA ac4C modifications and identify a role of tRNA ac4C in the regulation of mRNA translation in HNSCC.

HIF-1α/METTL1/m7G axis is involved in CRC response to hypoxia

BACKGROUND

In recent years, many studies have confirmed that hypoxia and hypoxia inducible factor (HIF)-1α drive the development of colorectal cancer (CRC). HIF-1α also modulates epitranscriptomic remodeling to regulate cancer development. However, the mechanism by which RNA methylation is altered under hypoxic conditions and the underlying regulatory mechanisms in CRC remain unclear.

METHODS

Here, seven common types of modifications of mRNA and tRNA were quantitated using liquid chromatography-tandem mass spectrometry. To validate the robustness of the profiling data, modifications that were consistently altered across the three CRC cell lines under hypoxia were validated via dot blot analysis. Then, 10 enzymes that could regulate the abundance of three RNA modifications in tRNA were measured in CRC cells after hypoxia treatment using quantitative real-time polymerase chain reaction. Furthermore, the regulatory role of HIF-1α in the expression of methyltransferase 1 (METTL1) under hypoxic conditions was confirmed using METTL1 promoter activity assays and HIF-1α small interfering RNA (siRNA). The binding capacity of HIF-1α to each hypoxia response element (HRE) in the promoter of METTL1 was investigated by performing Chromatin immunoprecipitation assay (ChIP).

RESULTS

Abundance of RNA modifications was altered more consistently and significantly in tRNA than in mRNA under hypoxic conditions. In addition, the abundance of N7-methyleguanosine (m7G) modification in tRNA decreased significantly under hypoxic conditions. As a methyltransferase of the m7G modification in tRNA, the expression of METTL1 mRNA was drastically downregulated under hypoxic conditions. Mechanistically, suppression of HIF-1α by siRNA upregulated the METTL1 promoter activity. Furthermore, ChIP showed that HIF-1α could bind with an HRE in the promoter region of METTL1, indicating that METTL1 is a direct target of HIF-1α in CRC cells under hypoxic conditions.

CONCLUSIONS

Our study revealed that the abundance of the m7G modification in tRNA was drastically reduced in CRC cells dependent on the HIF-1α-mediated inhibition of METTL1 transcription under hypoxic conditions.

Genetic variants of m7G modification genes influence neuroblastoma susceptibility

Objective

Neuroblastoma is a life-threatening pediatric solid tumor whose etiology remains unclear. N7-methylguanosine (m7G) is one of the most important epigenetic modifications of RNA, which plays a crucial role in tumorigenesis. The m7G-mediated genes METTL1 and WDR4 also have been reported to be dysregulated in various cancers. However, the implications of METTL1 and WDR4 in neuroblastoma have not been clarified.

Methods

Given the oncogenic potential of m7G modification, we performed a case-control study to assess the association of METTL1 and WDR4 genes polymorphisms with neuroblastoma risk in a Chinese population consisting of 402 cases and 473 controls. Odds ratios (ORs) and 95 % confidence intervals (CIs) were applied to evaluate the associations between studied polymorphisms and neuroblastoma risk. The adjusted odds ratio (AOR) was adjusted for age and gender.

Results

Overall, four polymorphisms were significantly associated with neuroblastoma risk, including METTL1 rs2291617 (recessive model: adjusted OR = 1.59, 95 % CI = 1.08-2.34, P = 0.019), WDR4 rs2156316 (dominant model: adjusted OR = 0.74, 95 % CI = 0.57-0.97, P = 0.028), WDR4 rs6586250 (dominant model: adjusted OR = 0.59, 95 % CI = 0.42-0.84, P = 0.004) and WDR4 rs15736 (dominant model: adjusted OR = 0.60, 95 % CI = 0.42-0.85, P = 0.004). Stratified analysis showed stronger correlations between significant polymorphisms and neuroblastoma risk among subgroups divided by age, gender, tumor origin, and clinical stage. Furthermore, expression quantitative trait loci (eQTL) analysis revealed that significant polymorphisms were associated with the expression of the adjacent genes.

Conclusions

Our study indicated that four polymorphisms in m7G-mediated genes contribute to neuroblastoma susceptibility in the eastern Chinese population. However, our findings should be verified further by large-scale and well-designed studies.

Global and single-nucleotide resolution detection of 7-methylguanosine in RNA

RNA modifications, including N-7-methylguanosine (m7G), are pivotal in governing RNA stability and gene expression regulation. The accurate detection of internal m7G modifications is of paramount significance, given recent associations between altered m7G deposition and elevated expression of the methyltransferase METTL1 in various human cancers. The development of robust m7G detection techniques has posed a significant challenge in the field of epitranscriptomics. In this study, we introduce two methodologies for the global and accurate identification of m7G modifications in human RNA. We introduce borohydride reduction sequencing (Bo-Seq), which provides base resolution mapping of m7G modifications. Bo-Seq achieves exceptional performance through the optimization of RNA depurination and scission, involving the strategic use of high concentrations of NaBH4, neutral pH and the addition of 7-methylguanosine monophosphate (m7GMP) during the reducing reaction. Notably, compared to NaBH4-based methods, Bo-Seq enhances the m7G detection performance, and simplifies the detection process, eliminating the necessity for intricate chemical steps and reducing the protocol duration. In addition, we present an antibody-based approach, which enables the assessment of m7G relative levels across RNA molecules and biological samples, however it should be used with caution due to limitations associated with variations in antibody quality between batches. In summary, our novel approaches address the pressing need for reliable and accessible methods to detect RNA m7G methylation in human cells. These advancements hold the potential to catalyse future investigations in the critical field of epitranscriptomics, shedding light on the complex regulatory roles of m7G in gene expression and its implications in cancer biology.

Epitranscriptomic modifications in mesenchymal stem cell differentiation: advances, mechanistic insights, and beyond

RNA modifications, known as the "epitranscriptome", represent a key layer of regulation that influences a wide array of biological processes in mesenchymal stem cells (MSCs). These modifications, catalyzed by specific enzymes, often termed "writers", "readers", and "erasers", can dynamically alter the MSCs' transcriptomic landscape, thereby modulating cell differentiation, proliferation, and responses to environmental cues. These enzymes include members of the classes METTL, IGF2BP, WTAP, YTHD, FTO, NAT, and others. Many of these RNA-modifying agents are active during MSC lineage differentiation. This review provides a comprehensive overview of the current understanding of different RNA modifications in MSCs, their roles in regulating stem cell behavior, and their implications in MSC-based therapies. It delves into how RNA modifications impact MSC biology, the functional significance of individual modifications, and the complex interplay among these modifications. We further discuss how these intricate regulatory mechanisms contribute to the functional diversity of MSCs, and how they might be harnessed for therapeutic applications. The review also highlights current challenges and potential future directions in the study of RNA modifications in MSCs, emphasizing the need for innovative tools to precisely map these modifications and decipher their context-specific effects. Collectively, this work paves the way for a deeper understanding of the role of the epitranscriptome in MSC biology, potentially advancing therapeutic strategies in regenerative medicine and MSC-based therapies.

Cancer metastasis under the magnifying glass of epigenetics and epitranscriptomics

Most of the cancer-associated mortality and morbidity can be attributed to metastasis. The role of epigenetic and epitranscriptomic alterations in cancer origin and progression has been extensively demonstrated during the last years. Both regulations share similar mechanisms driven by DNA or RNA modifiers, namely writers, readers, and erasers; enzymes responsible of respectively introducing, recognizing, or removing the epigenetic or epitranscriptomic modifications. Epigenetic regulation is achieved by DNA methylation, histone modifications, non-coding RNAs, chromatin accessibility, and enhancer reprogramming. In parallel, regulation at RNA level, named epitranscriptomic, is driven by a wide diversity of chemical modifications in mostly all RNA molecules. These two-layer regulatory mechanisms are finely controlled in normal tissue, and dysregulations are associated with every hallmark of human cancer. In this review, we provide an overview of the current state of knowledge regarding epigenetic and epitranscriptomic alterations governing tumor metastasis, and compare pathways regulated at DNA or RNA levels to shed light on a possible epi-crosstalk in cancer metastasis. A deeper understanding on these mechanisms could have important clinical implications for the prevention of advanced malignancies and the management of the disseminated diseases. Additionally, as these epi-alterations can potentially be reversed by small molecules or inhibitors against epi-modifiers, novel therapeutic alternatives could be envisioned.

A comprehensive model for tRNA methylation modification studies

The molecular mechanism of N7-methylguanosine (m7G)-tRNA modification mediated by the methyltransferase like 1-WD repeat domain 4 (METTL1-WDR4). Route A indicated that with WDR4 as the scaffold, METTL1 catalyzed G46 methylation modification after binding of the tRNA variable loop via its αC and α6 helices, as reported by Richard I. Gregory's group. Route B showed the METTL1-WDR4-mediated tRNA methylation modification model and the effect of binding of the S-adenosylmethionine on the N-terminus of METTL1, as reported by Yunsun Nam's group.

tRNA modifications: insights into their role in human cancers

Transfer RNA (tRNA) plays a central role in translation by functioning as a biological link between messenger RNA (mRNA) and proteins. One prominent feature of the tRNA molecule is its heavily modified status, which greatly affects its biogenesis and function. Modifications within the anticodon loop are crucial for translation efficiency and accuracy, whereas other modifications in the body region affect tRNA structure and stability. Recent research has revealed that these diverse modifications are critical regulators of gene expression. They are involved in many important physiological and pathological processes, including cancers. In this review we focus on six different tRNA modifications to delineate their functions and mechanisms in tumorigenesis and tumor progression, providing insights into their clinical potential as biomarkers and therapeutic targets.

A novel m7G regulator-based methylation patterns in head and neck squamous cell carcinoma

Abnormal RNA N7-methylguanosine (m7G) modification is known to contribute to effects on tumor occurrence and development. Nevertheless, the mechanisms of its function in immunoregulation, tumor microenvironment (TME) modulation, and tumor promotion remain largely unknown. A series of computer-aided bioinformatic analyses were conducted based on transcriptomic, single-cell sequence, and spatial transcriptomic data to determine the m7G modification patterns in head and neck squamous cell carcinoma (HNSCC). Consensus clustering approach was employed according to the expressions of 33 m7G regulators. ESTIMATE, CIBERSORT, and single sample gene set enrichment analysis algorithms were adopted to investigate the immune cell infiltration features. A prognostic model named m7Gscore was established. Seurat, SingleR, and Monocle2 were used to analyze the single-cell sequence profiling. STUtility was used to integrate multiple spatial transcriptomic datasets. Quantitative reverse transcription polymerase chain reaction, transwell, and wound-healing assay were performed to verify the oncogenes. Here, three different m7G modification patterns were highlighted in HNSCC patients, which were also related to various clinical manifestations and three representative immunophenotypes: immune-excluded, immune-desert, and inflamed, separately. Patients with lower m7Gscore were highlighted by higher immune cell infiltrations, better overall survival rates, lesser tumor mutation burden (TMB), lower sensitivities to target inhibitors therapies, and better immunotherapeutic response. Moreover, DCPS, EIF4E, EIF4E2, LSM1, NCBP2, NUDT1, and NUDT5 were identified to play critical roles in T-cell differentiation. Knockdown of LSM1/NUDT5 could restrain the malignancy of HNSCC cells. Collectively, quantitative assessment of m7G modification patterns in individual HNSCC patients could contribute to identifying more efficient immunotherapeutic approaches and improve the clinical outcome of HNSCC.

The m7G Reader NCBP2 Promotes Pancreatic Cancer Progression by Upregulating MAPK/ERK Signaling

PDAC is one of the most common malignant tumors worldwide. The difficulty of early diagnosis and lack of effective treatment are the main reasons for its poor prognosis. Therefore, it is urgent to identify novel diagnostic and therapeutic targets for PDAC patients. The m7G methylation is a common type of RNA modification that plays a pivotal role in regulating tumor development. However, the correlation between m7G regulatory genes and PDAC progression remains unclear. By integrating gene expression and related clinical information of PDAC patients from TCGA and GEO cohorts, m7G binding protein NCBP2 was found to be highly expressed in PDAC patients. More importantly, PDAC patients with high NCBP2 expression had a worse prognosis. Stable NCBP2-knockdown and overexpression PDAC cell lines were constructed to further perform in-vitro and in-vivo experiments. NCBP2-knockdown significantly inhibited PDAC cell proliferation, while overexpression of NCBP2 dramatically promoted PDAC cell growth. Mechanistically, NCBP2 enhanced the translation of c-JUN, which in turn activated MEK/ERK signaling to promote PDAC progression. In conclusion, our study reveals that m7G reader NCBP2 promotes PDAC progression by activating MEK/ERK pathway, which could serve as a novel therapeutic target for PDAC patients.

Selection of M7G-related lncRNAs in kidney renal clear cell carcinoma and their putative diagnostic and prognostic role

BACKGROUND

Kidney renal clear cell carcinoma (KIRC) is a common malignant tumor of the urinary system. This study aims to develop new biomarkers for KIRC and explore the impact of biomarkers on the immunotherapeutic efficacy for KIRC, providing a theoretical basis for the treatment of KIRC patients.

METHODS

Transcriptome data for KIRC was obtained from the The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Weighted gene co-expression network analysis identified KIRC-related modules of long noncoding RNAs (lncRNAs). Intersection analysis was performed differentially expressed lncRNAs between KIRC and normal control samples, and lncRNAs associated with N(7)-methylguanosine (m7G), resulting in differentially expressed m7G-associated lncRNAs in KIRC patients (DE-m7G-lncRNAs). Machine Learning was employed to select biomarkers for KIRC. The prognostic value of biomarkers and clinical features was evaluated using Kaplan-Meier (K-M) survival analysis, univariate and multivariate Cox regression analysis. A nomogram was constructed based on biomarkers and clinical features, and its efficacy was evaluated using calibration curves and decision curves. Functional enrichment analysis was performed to investigate the functional enrichment of biomarkers. Correlation analysis was conducted to explore the relationship between biomarkers and immune cell infiltration levels and common immune checkpoint in KIRC samples.

RESULTS

By intersecting 575 KIRC-related module lncRNAs, 1773 differentially expressed lncRNAs, and 62 m7G-related lncRNAs, we identified 42 DE-m7G-lncRNAs. Using XGBoost and Boruta algorithms, 8 biomarkers for KIRC were selected. Kaplan-Meier survival analysis showed significant survival differences in KIRC patients with high and low expression of the PTCSC3 and RP11-321G12.1. Univariate and multivariate Cox regression analyses showed that AP000696.2, PTCSC3 and clinical characteristics were independent prognostic factors for patients with KIRC. A nomogram based on these prognostic factors accurately predicted the prognosis of KIRC patients. The biomarkers showed associations with clinical features of KIRC patients, mainly localized in the cytoplasm and related to cytokine-mediated immune response. Furthermore, immune feature analysis demonstrated a significant decrease in immune cell infiltration levels in KIRC samples compared to normal samples, with a negative correlation observed between the biomarkers and most differentially infiltrating immune cells and common immune checkpoints.

CONCLUSION

In summary, this study discovered eight prognostic biomarkers associated with KIRC patients. These biomarkers showed significant correlations with clinical features, immune cell infiltration, and immune checkpoint expression in KIRC patients, laying a theoretical foundation for the diagnosis and treatment of KIRC.

Transfer RNAs as dynamic and critical regulators of cancer progression

Transfer RNAs (tRNAs) have been historically viewed as non-dynamic adaptors that decode the genetic code into proteins. Recent work has uncovered dynamic regulatory roles for these fascinating molecules. Advances in tRNA detection methods have revealed that specific tRNAs can become modulated upon DNA copy number and chromatin alterations and can also be perturbed by oncogenic signalling and transcriptional regulators in cancer cells or the tumour microenvironment. Such alterations in the levels of specific tRNAs have been shown to causally impact cancer progression, including metastasis. Moreover, sequencing methods have identified tRNA-derived small RNAs that influence various aspects of cancer progression, such as cell proliferation and invasion, and could serve as diagnostic and prognostic biomarkers or putative therapeutic targets in various cancers. Finally, there is accumulating evidence, including from genetic models, that specific tRNA synthetases - the enzymes responsible for charging tRNAs with amino acids - can either promote or suppress tumour formation. In this Review, we provide an overview of how deregulation of tRNAs influences cancer formation and progression.

The CTBP2-PCIF1 complex regulates m6Am modification of mRNA in head and neck squamous cell carcinoma

PCIF1 can mediate the methylation of N6,2'-O-dimethyladenosine (m6Am) in mRNA. Yet, the detailed interplay between PCIF1 and the potential cofactors and its pathological significance remain elusive. Here, we demonstrated that PCIF1-mediated cap mRNA m6Am modification promoted head and neck squamous cell carcinoma progression both in vitro and in vivo. CTBP2 was identified as a cofactor of PCIF1 to catalyze m6Am deposition on mRNA. CLIP-Seq data demonstrated that CTBP2 bound to similar mRNAs as compared with PCIF1. We then used the m6Am-Seq method to profile the mRNA m6Am site at single-base resolution and found that mRNA of TET2, a well-known tumor suppressor, was a major target substrate of the PCIF1-CTBP2 complex. Mechanistically, knockout of CTBP2 reduced PCIF1 occupancy on TET2 mRNA, and the PCIF1-CTBP2 complex negatively regulated the translation of TET2 mRNA. Collectively, our study demonstrates the oncogenic function of the epitranscriptome regulator PCIF1-CTBP2 complex, highlighting the importance of the m6Am modification in tumor progression.

N7-methylguanosine methylation of tRNAs regulates survival to stress in cancer

Tumour progression and therapy tolerance are highly regulated and complex processes largely dependent on the plasticity of cancer cells and their capacity to respond to stress. The higher plasticity of cancer cells highlights the need for identifying targetable molecular pathways that challenge cancer cell survival. Here, we show that N7-guanosine methylation (m7G) of tRNAs, mediated by METTL1, regulates survival to stress conditions in cancer cells. Mechanistically, we find that m7G in tRNAs protects them from stress-induced cleavage and processing into 5' tRNA fragments. Our analyses reveal that the loss of tRNA m7G methylation activates stress response pathways, sensitising cancer cells to stress. Furthermore, we find that the loss of METTL1 reduces tumour growth and increases cytotoxic stress in vivo. Our study uncovers the role of m7G methylation of tRNAs in stress responses and highlights the potential of targeting METTL1 to sensitise cancer cells to chemotherapy.

Biological roles of RNA m7G modification and its implications in cancer

M7G modification, known as one of the common post-transcriptional modifications of RNA, is present in many different types of RNAs. With the accurate identification of m7G modifications within RNAs, their functional roles in the regulation of gene expression and different physiological functions have been revealed. In addition, there is growing evidence that m7G modifications are crucial in the emergence of cancer. Here, we review the most recent findings regarding the detection techniques, distribution, biological functions and Regulators of m7G. We also summarize the connections between m7G modifications and cancer development, drug resistance, and tumor microenvironment as well as we discuss the research's future directions and trends.

RNA modification: mechanisms and therapeutic targets

RNA modifications are dynamic and reversible chemical modifications on substrate RNA that are regulated by specific modifying enzymes. They play important roles in the regulation of many biological processes in various diseases, such as the development of cancer and other diseases. With the help of advanced sequencing technologies, the role of RNA modifications has caught increasing attention in human diseases in scientific research. In this review, we briefly summarized the basic mechanisms of several common RNA modifications, including m6A, m5C, m1A, m7G, Ψ, A-to-I editing and ac4C. Importantly, we discussed their potential functions in human diseases, including cancer, neurological disorders, cardiovascular diseases, metabolic diseases, genetic and developmental diseases, as well as immune disorders. Through the "writing-erasing-reading" mechanisms, RNA modifications regulate the stability, translation, and localization of pivotal disease-related mRNAs to manipulate disease development. Moreover, we also highlighted in this review all currently available RNA-modifier-targeting small molecular inhibitors or activators, most of which are designed against m6A-related enzymes, such as METTL3, FTO and ALKBH5. This review provides clues for potential clinical therapy as well as future study directions in the RNA modification field. More in-depth studies on RNA modifications, their roles in human diseases and further development of their inhibitors or activators are needed for a thorough understanding of epitranscriptomics as well as diagnosis, treatment, and prognosis of human diseases.

P300/SP1 complex mediating elevated METTL1 regulates CDK14 mRNA stability via internal m7G modification in CRPC

BACKGROUND

N7-methylguanosine (m7G) modification is, a more common epigenetic modification in addition to m6A modification, mainly found in mRNA capsids, mRNA interiors, transfer RNA (tRNA), pri-miRNA, and ribosomal RNA (rRNA). It has been found that m7G modifications play an important role in mRNA transcription, tRNA stability, rRNA processing maturation, and miRNA biosynthesis. However, the role of m7G modifications within mRNA and its "writer" methyltransferase 1(METTL1) in tumors, particularly prostate cancer (PCa), has not been revealed.

METHODS

The differential expression level of METTL1 between hormone-sensitive prostate cancer (HSPC) and castrate-resistant prostate cancer (CRPC) was evaluated via RNA-seq and in vitro experiments. The effects of METTL1 on CRPC progression were investigated through in vitro and in vivo assays. The upstream molecular mechanism of METTL1 expression upregulation and the downstream mechanism of its action were explored via Chromatin Immunoprecipitation quantitative reverse transcription polymerase chain reaction (CHIP-qPCR), Co-immunoprecipitation (Co-IP), luciferase reporter assay, transcriptome-sequencing, m7G AlkAniline-Seq, and mRNA degradation experiments, etc. RESULTS AND CONCLUSION: Here, we found that METTL1 was elevated in CRPC and that patients with METTL1 elevation tended to have a poor prognosis. Functionally, the knockdown of METTL1 in CRPC cells significantly limited cell proliferation and invasive capacity. Mechanistically, we unveiled that P300 can form a complex with SP1 and bind to the promoter region of the METTL1 gene via SP1, thereby mediating METTL1 transcriptional upregulation in CRPC. Subsequently, our findings indicated that METTL1 leads to enhanced mRNA stability of CDK14 by adding m7G modifications inside its mRNA, ultimately promoting CRPC progression.

m7G-related genes predict prognosis and affect the immune microenvironment and drug sensitivity in osteosarcoma

Background: Osteosarcoma (OS), a primary malignant bone tumor, confronts therapeutic challenges rooted in multidrug resistance. Comprehensive understanding of disease occurrence and progression is imperative for advancing treatment strategies. m7G modification, an emerging post-transcriptional modification implicated in various diseases, may provide new insights to explore OS pathogenesis and progression. Methods: The m7G-related molecular landscape in OS was probed using diverse bioinformatics analyses, encompassing LASSO Cox regression, immune infiltration assessment, and drug sensitivity analysis. Furthermore, the therapeutic potential of AZD2014 for OS was investigated through cell apoptosis and cycle assays. Eventually, multivariate Cox analysis and experimental validations, were conducted to investigate the independent prognostic m7G-related genes. Results: A comprehensive m7G-related risk model incorporating eight signatures was established, with corresponding risk scores correlated with immune infiltration and drug sensitivity. Drug sensitivity analysis spotlighted AZD2014 as a potential therapeutic candidate for OS. Subsequent experiments corroborated AZD2014's capability to induce G1-phase cell cycle arrest and apoptosis in OS cells. Ultimately, multivariate Cox regression analysis unveiled the independent prognostic importance of CYFIP1 and EIF4A1, differential expressions of which were validated at histological and cytological levels. Conclusion: This study furnishes a profound understanding of the contribution of m7G-related genes to the pathogenesis of OS. The discerned therapeutic potential of AZD2014, in conjunction with the identification of CYFIP1 and EIF4A1 as independent risk factors, opens novel vistas for the treatment of OS.

Methyltransferase-like proteins in cancer biology and potential therapeutic targeting

RNA modification has recently become a significant process of gene regulation, and the methyltransferase-like (METTL) family of proteins plays a critical role in RNA modification, methylating various types of RNAs, including mRNA, tRNA, microRNA, rRNA, and mitochondrial RNAs. METTL proteins consist of a unique seven-beta-strand domain, which binds to the methyl donor SAM to catalyze methyl transfer. The most typical family member METTL3/METTL14 forms a methyltransferase complex involved in N6-methyladenosine (m6A) modification of RNA, regulating tumor proliferation, metastasis and invasion, immunotherapy resistance, and metabolic reprogramming of tumor cells. METTL1, METTL4, METTL5, and METTL16 have also been recently identified to have some regulatory ability in tumorigenesis, and the rest of the METTL family members rely on their methyltransferase activity for methylation of different nucleotides, proteins, and small molecules, which regulate translation and affect processes such as cell differentiation and development. Herein, we summarize the literature on METTLs in the last three years to elucidate their roles in human cancers and provide a theoretical basis for their future use as potential therapeutic targets.

Comprehensive analysis of circRNAs for N7-methylguanosine methylation modification in human oral squamous cell carcinoma

N7-methylguanosine (m7G) modification is closely related to the occurrence of tumors. However, the m7G modification of circRNAs in oral squamous cell carcinoma (OSCC) remains to be investigated. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) was used to measure the methylation levels of m7G and identify m7G sites in circRNAs in human OSCC and normal tissues. The host genes of differentially methylated and differentially expressed circRNAs were analyzed by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, and circRNA-miRNA-mRNA networks were predicted using the miRanda and miRDB databases. The analysis identified 2348 m7G peaks in 624 circRNAs in OSCC tissues. In addition, the source of m7G-methylated circRNAs in OSCC was mainly the sense overlap region compared with normal tissues. The most conserved m7G motif in OSCC tissues was CCUGU, whereas the most conserved motif in normal tissues was RCCUG (R = G/A). Importantly, GO enrichment and KEGG pathway analysis showed that the host genes of differentially methylated and differentially expressed circRNAs were involved in many cellular biological functions. Furthermore, the significantly differentially expressed circRNAs were analyzed to predict the circRNA-miRNA-mRNA networks. This study revealed the whole profile of circRNAs of differential m7G methylation in OSCC and suggests that m7G-modified circRNAs may impact the development of OSCC.

RNA modifications in cancer

Currently, more than 170 modifications have been identified on RNA. Among these RNA modifications, various methylations account for two-thirds of total cases and exist on almost all RNAs. Roles of RNA modifications in cancer are garnering increasing interest. The research on m6A RNA methylation in cancer is in full swing at present. However, there are still many other popular RNA modifications involved in the regulation of gene expression post-transcriptionally besides m6A RNA methylation. In this review, we focus on several important RNA modifications including m1A, m5C, m7G, 2'-O-Me, Ψ and A-to-I editing in cancer, which will provide a new perspective on tumourigenesis by peeking into the complex regulatory network of epigenetic RNA modifications, transcript processing, and protein translation.

M7G methylated core genes (METTL1 and WDR4) and associated RNA risk signatures are associated with prognosis and immune escape in HCC

N7 methylguanosine (m7G) has a crucial role the development of hepatocellular carcinoma (HCC). This study aimed to investigate the impact of the m7G methylation core genes (METTL1 and WDR4) and associated RNA risk signatures on HCC. we found m7G methylation core genes (METTL1 and WDR4) were upregulated in four HCC cell lines, and downregulation of METTL1 and WDR4 attenuated HCC cell proliferation, migration, and invasion. Moreover, METTL1 and WDR4 are upregulated in HCC tissues, and that there is a significant positive correlation between them. METTL1 and WDR4 were identified as independent prognostic markers for HCC by employing overall survival (OS), disease-specific survival (DSS), Progression Free Interval survival (PFI), and univariate/multivariate Cox analyses. We identified 1479 coding RNAs (mRNAs) and 232 long non-coding RNAs (lncRNAs) associated with METTL1 / WDR4 by using weighted coexpression network analysis (WGCNA) and co-clustering analysis. The least absolute shrinkage and selection operator (lasso) were used to constructing mRNA and lncRNA risk signatures associated with the METTL1 / WDR4. These risk were independent poor prognostic factors in HCC. Furthermore, we found that METTL1 / WDR4 expression and mRNA / lncRNA risk scores were closely associated with TP53 mutations. Clinicopathological features correlation results showed that METTL1 / WDR4 expression and mRNA / lncRNA risk score were associated with the stage and invasion depth (T) of HCC. To predict the overall survival of HCC individuals, we constructed a nomogram with METTL1/WDR4 expression, mRNA/lncRNA risk score, and clinicopathological features. In addition, we combined single-cell sequencing datasets and immune escape-related checkpoints to construct an immune escape-related protein-protein interaction(PPI) network. In conclusion, M7G methylated core genes (METTL1 and WDR4) and associated RNA risk signatures are associated with prognosis and immune escape in HCC.

Identification and validation of a novel prognostic model for gastric cancer based on m7G-related genes

Background

The role of N7-methyladenosine (m7G)-related genes in the progression and prognosis of gastric cancer (GC) remains unclear. This study aimed to explore prognostic biomarkers for GC based on m7G methylation regulators and to construct a prognostic risk model.

Methods

RNA sequencing profiles with corresponding clinicopathological information associated with GC of which the histological type was stomach adenocarcinoma (STAD) were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), respectively. A total of 29 m7G regulators were extracted from previous studies. According to the expression similarity of m7G regulators, the GC samples obtained from TCGA were further classified into 2 clusters demonstrating different overall survival (OS) rates and genetic heterogeneity, and the differentially expressed genes (DEGs) between these 2 clusters were defined as m7G-related genes. Univariate regression analysis and regression analysis were then used to obtain the prognostic m7G-related genes. The samples in TCGA and Genotype-Tissue Expression (GTEx) were used to verify the differential expression and prognostic value of these m7G-related genes contained in the prognostic model. Subsequently, the risk score was combined with other prognostic factors to develop a nomogram. The predictive ability of the nomogram was evaluated by the standard receiver operating characteristic (ROC) curve. Gene set enrichment analysis (GSEA) was used to identify activation pathways in both groups. Finally, the association between the prognostic model and the immune characteristics of GC were appraised.

Results

A prognostic model consisting of 11 m7G-related genes was constructed. GC patients in the high-risk group were shown to have a poor prognosis and this result was further demonstrated in each group. The risk model can be applied for patients with different clinical features. The results of GSEA showed that cell adhesion, cell junction, and focal adhesion were highly enriched in the high-risk group. In addition, we found that the expression of programmed cell death ligand 1 (PD-L1) was significantly elevated in the low-risk group, whereas programmed cell death ligand 2 (PD-L2) and tumor necrosis factor receptor superfamily member 4 (TNFRSF4) were overexpressed in the high-risk group.

Conclusions

We successfully built and verified a m7G relevant prognostic model for predicting prognosis and providing a new train of thought for improving the treatment of GC.

Fast and Efficient Design of Deep Neural Networks for Predicting N7-Methylguanosine Sites Using autoBioSeqpy

N⁷-Methylguanosine (m⁷G) is a crucial post-transcriptional RNA modification that plays a pivotal role in regulating gene expression. Accurately identifying m⁷G sites is a fundamental step in understanding the biological functions and regulatory mechanisms associated with this modification. While whole-genome sequencing is the gold standard for RNA modification site detection, it is a time-consuming, expensive, and intricate process. Recently, computational approaches, especially deep learning (DL) techniques, have gained popularity in achieving this objective. Convolutional neural networks and recurrent neural networks are examples of DL algorithms that have emerged as versatile tools for modeling biological sequence data. However, developing an efficient network architecture with superior performance remains a challenging task, requiring significant expertise, time, and effort. To address this, we previously introduced a tool called autoBioSeqpy, which streamlines the design and implementation of DL networks for biological sequence classification. In this study, we utilized autoBioSeqpy to develop, train, evaluate, and fine-tune sequence-level DL models for predicting m⁷G sites. We provided detailed descriptions of these models, along with a step-by-step guide on their execution. The same methodology can be applied to other systems dealing with similar biological questions. The benchmark data and code utilized in this study can be accessed for free at http://github.com/jingry/autoBioSeeqpy/tree/2.0/examples/m7G.

Epigenetic regulation in the tumor microenvironment: molecular mechanisms and therapeutic targets

Over decades, researchers have focused on the epigenetic control of DNA-templated processes. Histone modification, DNA methylation, chromatin remodeling, RNA modification, and noncoding RNAs modulate many biological processes that are crucial to the development of cancers. Dysregulation of the epigenome drives aberrant transcriptional programs. A growing body of evidence suggests that the mechanisms of epigenetic modification are dysregulated in human cancers and might be excellent targets for tumor treatment. Epigenetics has also been shown to influence tumor immunogenicity and immune cells involved in antitumor responses. Thus, the development and application of epigenetic therapy and cancer immunotherapy and their combinations may have important implications for cancer treatment. Here, we present an up-to-date and thorough description of how epigenetic modifications in tumor cells influence immune cell responses in the tumor microenvironment (TME) and how epigenetics influence immune cells internally to modify the TME. Additionally, we highlight the therapeutic potential of targeting epigenetic regulators for cancer immunotherapy. Harnessing the complex interplay between epigenetics and cancer immunology to develop therapeutics that combine thereof is challenging but could yield significant benefits. The purpose of this review is to assist researchers in understanding how epigenetics impact immune responses in the TME, so that better cancer immunotherapies can be developed.

WDR4 gene polymorphisms and Wilms tumor susceptibility in Chinese children: A five-center case-control study

Wilms tumor is the most common embryonal renal malignancy in children. WDR4 is an indispensable noncatalytic subunit of the RNA N7-methylguanosine (m7G) methyltransferase complex and plays an essential role in tumorigenesis. However, the relationship between polymorphisms in the WDR4 gene and susceptibility to Wilms tumor remains to be fully investigated. We performed a large case-control study involving 414 patients and 1199 cancer-free controls to investigate whether single nucleotide polymorphisms (SNPs) in the WDR4 gene are associated with Wilms tumor susceptibility. WDR4 gene polymorphisms (rs2156315 C > T, rs2156316 C > G, rs6586250 C > T, rs15736 G > A, and rs2248490 C > G) were genotyped using the TaqMan assay. In addition, unconditioned logistic regression analysis was performed, odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the association between WDR4 gene SNPs and Wilms tumor susceptibility as well as the strength of the associations. We found that only the rs6586250 C>T polymorphism was significantly associated with an increased risk of Wilms tumor (adjusted OR=2.99, 95% CI = 1.28-6.97, P = 0.011 for the rs6586250 TT genotype; adjusted OR=3.08, 95% CI = 1.33-7.17, P = 0.009 for the rs6586250 CC/CT genotype). Furthermore, the stratification analysis revealed that patients with the rs6586250 TT genotype and carriers with 1-5 risk genotypes exhibited statistically significant associations with increased Wilms tumor risk in specific subgroups. However, the rs2156315 CT/TT genotype was identified as having a protective effect against Wilms tumor in the age >18 months subgroup compared with the rs2156315 CC genotype. In brief, our study demonstrated that the rs6586250 C > T polymorphism of the WDR4 gene was significantly associated with Wilms tumor. This finding may contribute to the understanding of the genetic mechanism of Wilms tumor.

When N7-methyladenosine modification meets cancer: Emerging frontiers and promising therapeutic opportunities

N⁷-methylguanosine (m⁷G) methylation, one of the most common RNA modifications in eukaryotes, has recently gained considerable attention. The biological functions of m⁷G modification in RNAs, including tRNA, rRNA, mRNA, and miRNA, remain largely unknown in human diseases. Owing to rapid advances in high-throughput technologies, increasing evidence suggests that m⁷G modification plays a critical role in cancer initiation and progression. As m⁷G modification and hallmarks of cancer are inextricably linked together, targeting m⁷G regulators may provide new possibilities for future cancer diagnoses and potential intervention targets. This review summarizes various detection methods for m⁷G modification, recent advances in m⁷G modification and tumor biology regarding their interplay and regulatory mechanisms. We conclude with an outlook on the future of diagnosing and treating m⁷G-related diseases.

N7-methylguanosin regulators-mediated methylation modification patterns and characterization of the immune microenvironment in lower-grade glioma

N7-methylguanosine (m7G) modification signature has recently emerged as a crucial regulator of tumor progression and treatment in cancer. However, there is limited information available on the genomic profile of lower-grade gliomas (LGGs) related to m7G methylation modification genes' function in tumorigenesis and progression. In this study, we employed bioinformatics methods to characterize m7G modifications in individuals with LGG from The Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA). We used gene set enrichment analysis (GSEA), single sample GSEA (ssGSEA), CIBERSORT algorithm, ESTIMATE algorithm, and TIDE to evaluate the association between m7G modification patterns, tumor microenvironment (TME) cell infiltration properties, and immune infiltration markers. The m7G scoring scheme using principal component analysis (PCA) was employed to investigate the m7G modification patterns quantitatively. We examined the m7G modification hub genes' expression levels in normal samples, refractory epilepsy samples, and LGG samples using immunohistochemistry, western-blotting, and qRT-PCR. Our findings revealed that individuals with LGG could be categorized into two groups based on m7G scores (high and low) according to the properties of m7G. Moreover, we observed that high m7G score was associated with significant clinical benefit and prolonged survival duration in the anti-PD-1 cohort, while low m7G score was associated with improved prognostic outcomes and increased likelihood of complete or partial response in the anti-PD-L1 cohort. Different m7G subtypes also showed varying Tumor Mutational Burden (TMB) and immune profiles and might have distinct responses to immunotherapy. Furthermore, we identified five potential genetic markers that were highly correlated with the m7G score signature index. These findings provide insight into the features and classification associated with m7G methylation modifications and may aid in improving the clinical outcome of LGG.

Single cell sequencing analysis constructed the N7-methylguanosine (m7G)-related prognostic signature in uveal melanoma

BACKGROUND

Uveal melanoma is a highly malignant tumor in the eye. Its recurrence and metastasis are common, and the prognosis is poor.

METHODS

The transcriptome data of UVM were downloaded from TCGA database, and the single cell sequencing dataset GSE139829 was downloaded from GEO database. Weighted co-expression network analysis was used to explore the modules associated with m7G. Lasso regression was used to construct M⁷G-related prognostic signature. Immune infiltration analysis was used to explore the significance of the model in the tumor immune microenvironment. Finally, cell assays were used to explore the function of key genes in the MUM-2B and OCM-1 cell lines of UVM.

RESULTS

The prognostic signature was constructed by Cox regression and Lasso regression. Patients could be divided into high-risk group and low-risk group by this signature, and the high-risk group had worse prognosis (P<0.05). Cell experiments showed that the proliferation, invasion and migration ability of UVM cell lines were significantly decreased after the knockdown of PAG1, a key gene in signature, which proved that PAG1 might be a potential target of UVM.

CONCLUSIONS

Our study explored the significance of m⁷G in UVM, provided biomarkers for its diagnosis and treatment.

Current and potential roles of RNA modification-mediated autophagy dysregulation in cancer

Autophagy, a cellular lysosomal degradation and survival pathway, supports nutrient recycling and adaptation to metabolic stress and participates in various stages of tumor development, including tumorigenesis, metastasis, and malignant state maintenance. Among the various factors contributing to the dysregulation of autophagy in cancer, RNA modification can regulate autophagy by directly affecting the expression of core autophagy proteins. We propose that autophagy disorder mediated by RNA modification is an important mechanism for cancer development. Therefore, this review mainly discusses the role of RNA modification-mediated autophagy regulation in tumorigenesis. We summarize the molecular basis of autophagy and the core proteins and complexes at different stages of autophagy, especially those involved in cancer development. Moreover, we describe the crosstalk of RNA modification and autophagy and review the recent advances and potential role of the RNA modification/autophagy axis in the development of multiple cancers. Furthermore, the dual role of the RNA modification/autophagy axis in cancer drug resistance is discussed. A comprehensive understanding and extensive exploration of the molecular crosstalk of RNA modifications with autophagy will provide important insights into tumor pathophysiology and provide more options for cancer therapeutic strategies.

Association of RNA m(7)G Modification Gene Polymorphisms with Pediatric Glioma Risk

Glioma stemming from glial cells of the central nervous system (CNS) is one of the leading causes of cancer death in childhood. The genetic predisposition of glioma is not fully understood. METTL1-WDR4 methyltransferase complex is implicated in tumorigenesis by catalyzing N7-methylguanosine (m⁷G) modification of RNA. This study is aimed at determining the association of glioma risk with three polymorphisms (rs2291617, rs10877013, and rs10877012) in METTL1 and five polymorphisms (rs2156315 rs2156316, rs6586250, rs15736, and rs2248490) in WDR4 gene in children of Chinese Han. We enrolled 314 cases and 380 controls from three independent hospitals. Genotypes of these polymorphisms were determined using the TaqMan assay. We found the WDR4 gene rs15736 was significantly associated with reduced glioma risk (GA/AA vs. GG: adjusted odds ratio = 0.63, 95%confidence interval = 0.42 - 0.94, P = 0.023) out of the eight studied polymorphisms. Stratified analyses showed that the association of rs15736 with the risk of glioma remained significant in children aged 60 months or older, girls, the subgroups with astrocytic tumors, or grade I + II glioma. We also found the combined effects of five WDR4 gene polymorphisms on glioma risk. Finally, expression quantitative trait locus (eQTL) analyses elucidated that the rs15736 polymorphism was related to the expression level of WDR4 and neighboring gene cystathionine-beta-synthase (CBS). Our finding provided evidence of a causal association between WDR4 gene polymorphisms and glioma susceptibility in Chinese Han children.

Structures and mechanisms of tRNA methylation by METTL1-WDR4

Specific, regulated modification of RNAs is important for proper gene expression^1,2. tRNAs are rich with various chemical modifications that affect their stability and function^3,4. 7-Methylguanosine (m⁷G) at tRNA position 46 is a conserved modification that modulates steady-state tRNA levels to affect cell growth^5,6. The METTL1-WDR4 complex generates m⁷G46 in humans, and dysregulation of METTL1-WDR4 has been linked to brain malformation and multiple cancers^7-22. Here we show how METTL1 and WDR4 cooperate to recognize RNA substrates and catalyse methylation. A crystal structure of METTL1-WDR4 and cryo-electron microscopy structures of METTL1-WDR4-tRNA show that the composite protein surface recognizes the tRNA elbow through shape complementarity. The cryo-electron microscopy structures of METTL1-WDR4-tRNA with S-adenosylmethionine or S-adenosylhomocysteine along with METTL1 crystal structures provide additional insights into the catalytic mechanism by revealing the active site in multiple states. The METTL1 N terminus couples cofactor binding with conformational changes in the tRNA, the catalytic loop and the WDR4 C terminus, acting as the switch to activate m⁷G methylation. Thus, our structural models explain how post-translational modifications of the METTL1 N terminus can regulate methylation. Together, our work elucidates the core and regulatory mechanisms underlying m⁷G modification by METTL1, providing the framework to understand its contribution to biology and disease.

N7-methylguanosine regulatory genes well represented by METTL1 define vastly different prognostic, immune and therapy landscapes in adrenocortical carcinoma

Although N7-methylguanosine (m7G) is one of the most frequent RNA modifications, it has received little attention. Adrenocortical carcinoma (ACC) is a highly malignant and easily metastatic tumor, eagerly needing for novel therapeutic strategy. Herein, a novel m7G risk signature (METTL1, NCBP1, NUDT1 and NUDT5) was constructed using the Lasso regression analysis. It possessed highly prognostic value and could improve the predictive accuracy and clinical making-decision benefit of traditional prognostic model. Its prognostic value was also successfully validated in GSE19750 cohort. Through CIBERSORT, ESTIMATE, ssGSEA and GSEA analyzes, high-m7G risk score was found to be closely associated with increased enrichment of glycolysis and suppression of anti-cancer immune response. Therapeutic correlation of m7G risk signature was also investigated using tumor mutation burden, the expressions of immune checkpoints, TIDE score, IMvigor 210 cohort and TCGA cohort. m7G risk score was a potential biomarker for predicting the efficacy of ICBs and mitotane. Furthermore, we explored the biofunctions of METTL1 in ACC cells through a series of experimentations. Overexpression of METTL1 stimulated the proliferation, migration and invasion of H295R and SW13 cells. Immunofluorescence assays revealed that the infiltrating levels of CD8+ T cells was lower and that of macrophages was higher in clinical ACC samples with high METTL1 expression compared to that in low expression ones. Silencing METTL1 could significantly inhibited tumor growth in mouse xenograft model. Western blot assays showed that METTL1 positively regulated the expression of glycolysis rate-limiting enzyme HK1. Finally, miR-885-5p and CEBPB were predicted as the upstream regulators of METTL1 through data mining of the public databases. In conclusions, m7G regulatory genes well represented by METTL1 profoundly affected the prognosis, tumor immune, therapeutic outcomes, and malignant progression of ACC.

Structural basis of regulated m7G tRNA modification by METTL1-WDR4

Chemical modifications of RNA have key roles in many biological processes1-3. N7-methylguanosine (m7G) is required for integrity and stability of a large subset of tRNAs4-7. The methyltransferase 1-WD repeat-containing protein 4 (METTL1-WDR4) complex is the methyltransferase that modifies G46 in the variable loop of certain tRNAs, and its dysregulation drives tumorigenesis in numerous cancer types8-14. Mutations in WDR4 cause human developmental phenotypes including microcephaly15-17. How METTL1-WDR4 modifies tRNA substrates and is regulated remains elusive18. Here we show,  through structural, biochemical and cellular studies of human METTL1-WDR4, that WDR4 serves as a scaffold for METTL1 and the tRNA T-arm. Upon tRNA binding, the αC region of METTL1 transforms into a helix, which together with the α6 helix secures both ends of the tRNA variable loop. Unexpectedly, we find that the predicted disordered N-terminal region of METTL1 is part of the catalytic pocket and essential for methyltransferase activity. Furthermore, we reveal that S27 phosphorylation in the METTL1 N-terminal region inhibits methyltransferase activity by locally disrupting the catalytic centre. Our results provide a molecular understanding of tRNA substrate recognition and phosphorylation-mediated regulation of METTL1-WDR4, and reveal the presumed disordered N-terminal region of METTL1 as a nexus of methyltransferase activity.

RNA Modifications in Cancer Metabolism and Tumor Microenvironment

RNA modifications have recently been recognized as essential posttranscriptional regulators of gene expression in eukaryotes. Investigations over the past decade have revealed that RNA chemical modifications have profound effects on tumor initiation, progression, refractory, and recurrence. Tumor cells are notorious for their robust plasticity in response to the stressful microenvironment and undergo metabolic adaptations to sustain rapid cell proliferation, which is termed as metabolic reprogramming. Meanwhile, cancer-associated metabolic reprogramming leads to substantial alterations of intracellular and extracellular metabolites, which further reshapes the tumor microenvironment (TME). Moreover, cancer cells compete with tumor-infiltrating immune cells for the limited nutrients to maintain their proliferation and function in the TME. In this chapter, we review recent interesting findings on the engagement of epitranscriptomic pathways, especially the ones associated with N6-methyladenosine (m6A), in the regulation of cancer metabolism and the surrounding microenvironment. We also discuss the promising therapeutic approaches targeting RNA modifications for anti-tumor therapy.

tRNA Modifications and Modifying Enzymes in Disease, the Potential Therapeutic Targets

tRNA is one of the most conserved and abundant RNA species, which plays a key role during protein translation. tRNA molecules are post-transcriptionally modified by tRNA modifying enzymes. Since high-throughput sequencing technology has developed rapidly, tRNA modification types have been discovered in many research fields. In tRNA, numerous types of tRNA modifications and modifying enzymes have been implicated in biological functions and human diseases. In our review, we talk about the relevant biological functions of tRNA modifications, including tRNA stability, protein translation, cell cycle, oxidative stress, and immunity. We also explore how tRNA modifications contribute to the progression of human diseases. Based on previous studies, we discuss some emerging techniques for assessing tRNA modifications to aid in discovering different types of tRNA modifications.

RNA Modifications in Cancer Stem Cell Biology

Post-transcriptional regulation of gene expression shapes the cell state both in health and disease. RNA modifications-especially N6-methyladenosine (m6A)-have recently emerged as key players in RNA processing that depends on a sophisticated interplay between proteins of the RNA modification machinery. Importantly, the RNA epitranscriptome becomes dysregulated in cancer and promotes cancer-associated gene expression programs as well as cancer cell adaptation to the tumor microenvironment. At the top of the tumor hierarchy, cancer stem cells (CSCs) are master regulators of tumorigenesis and resistance to therapeutic intervention. Therefore, defining how RNA modifications influence the CSC state is of great importance for cancer drug development. In this chapter, we summarize the current knowledge of the roles of RNA modifications in shaping the CSC state and driving gene expression programs that confer stem-like properties to CSCs, promote CSC adaptation to the local microenvironment, and endow CSCs with metastatic potential and drug resistance.

N7-methylguanosine methylation-related regulator genes as biological markers in predicting prognosis for melanoma

The aim of this study is to find those N7-methylguanosine (m7G) methylation-related regulator genes (m7GMRRGs) which were associated with melanoma prognosis and use them to develop a prognostic prediction model. Clinical information was retrieved online from The Cancer Gene Atlas (TCGA) and the Gene Expression Omnibus (GEO). R software was used to extract m⁷GMRRGs by differential expression analysis. To create a prognostic risk model, univariate and multivariate Cox regression analyses were employed for the evaluation of the prognostic significance of m⁷G methylation modifiers. Internal validation using cohort from TCGA (training set) and external validation using cohort from GEO (validation set) of the model were carried out. The model's predictive performance was confirmed by using the Kaplan-Meier, univariate, and multivariate Cox regression, and receiver operating characteristic curve (ROC) by constructing column line plots incorporating clinical factor characteristics. Immune infiltration analyses were performed to assess the immune function of m⁷GMRRGs. Drug sensitivity analysis was conducted to study chemotherapeutic drug treatment cues. Prognostic models using four m⁷GMRRGs (EIF4E3, LARP1, NCBP3, and IFIT5) showed good prognostic power in training and validation sets. The area under the curve (AUC) at 1, 3, and 5 years for GEO-melanoma were 0.689, 0.704, and 0.726, respectively. The prediction model could distinctly classify patients with melanoma into different risk subgroups (P < 0.001 for TCGA-melanoma and P < 0.05 for GEO-melanoma). Clinical characteristics were taken into account in Cox regression and AUC analysis, which highlighted that the risk score served as an independent risk factor determining the prognosis of patients with melanoma. Immuno-infiltration analysis showed that m⁷GMRRGs could potentially regulate CD8⁺ T cells as well as regulatory T cells (Treg cells). Results of our study indicate a association between m⁷GMRRGs and melanoma prognosis, and the prognostic prediction model using m⁷GMRRGs may predict the prognosis of patients with melanoma well. Nevertheless, these results may provide a clue for potential better options of melanoma treatment but need further validation in futural studies.

Construction of m7G subtype classification on heterogeneity of sepsis

Sepsis is a highly heterogeneous disease and a major factor in increasing mortality from infection. N7-Methylguanosine (m7G) is a widely RNA modification in eukaryotes, which involved in regulation of different biological processes. Researchers have found that m7G methylation contributes to a variety of human diseases, but its research in sepsis is still limited. Here, we aim to establish the molecular classification of m7G gene-related sepsis, reveal its heterogeneity and explore the underlying mechanism. We first identified eight m7G related prognostic genes, and identified two different molecular subtypes of sepsis through Consensus Clustering. Among them, the prognosis of C2 subtype is worse than that of C1 subtype. The signal pathways enriched by the two subtypes were analyzed by ssGSEA, and the results showed that the amino acid metabolism activity of C2 subtype was more active than that of C1 subtype. In addition, the difference of immune microenvironment among different subtypes was explored through CIBERSORT algorithm, and the results showed that the contents of macrophages M0 and NK cells activated were significantly increased in C2 subtype, while the content of NK cells resting decreased significantly in C2 subtype. We further explored the relationship between immune regulatory genes and inflammation related genes between C2 subtype and C1 subtype, and found that C2 subtype showed higher expression of immune regulatory genes and inflammation related genes. Finally, we screened the key genes in sepsis by WGCNA analysis, namely NUDT4 and PARN, and verified their expression patterns in sepsis in the datasets GSE131761 and GSE65682. The RT-PCR test further confirmed the increased expression of NUDTA4 in sepsis patients. In conclusion, sepsis clustering based on eight m7G-related genes can well distinguish the heterogeneity of sepsis patients and help guide the personalized treatment of sepsis patients.

Identification of m7G Methylation-Related miRNA Signature Associated with Survival and Immune Microenvironment Regulation in Uterine Corpus Endometrial Carcinoma

Background. N7-methylguanosine (m7G) has been implicated in the development of cancer. The role of m7G-related miRNAs in the survival prediction of UCEC patients has not been investigated. Current research was the first to construct an m7G-related miRNA model to accurately predict the survival of patients with uterine corpus endometrial carcinoma (UCEC) and to explore immune cell infiltration and immune activity in the tumor microenvironment. Methods. RNA-seq data and clinical information of UCEC patients were derived from The Cancer Genome Atlas (TCGA) database. Using the TargetScan online database, we predicted miRNAs linked to the m7G-related genes and identified miRNAs which were significantly associated with the survival in UCEC patients and constructed a risk scoring model. The TCGA-UCEC cases were scored according to the risk model, and the high- and low-risk groups were divided by the median risk value. Gene enrichment analysis and immune cell infiltration and immune function analysis were performed using “clusterProfiler” and “GSVA” packages in R. Results. The survival prediction model consisted of 9 miRNAs, namely, hsa-miR-1301, hsa-miR-940, hsa-miR-592, hsa-miR-3170, hsa-miR-876, hsa-miR-215, hsa-miR-934, hsa-miR-3920, and hsa-miR-216b. Survival of UCEC patients in the high-risk group was worse than that in the low-risk group ( $p<0.001$ ). The receiver operating characteristic (ROC) curve showed that the model had good predictive performance, and the area under the curve was 0.800, 0.690, and 0.705 for 1-, 3-, and 5-year survival predictions, respectively. There were differences in the degree of immune cell infiltration and immune activity between the low-risk and high-risk groups. The expression levels of the identified differentially expressed genes correlated with the susceptibility to multiple anticancer drugs. Conclusions. The survival prediction model constructed based on 9 m7G-related miRNAs had good predictive performance.

M7G-related LncRNAs: A comprehensive analysis of the prognosis and immunity in glioma

Today, numerous international researchers have demonstrated that N7-methylguanosine (m7G) related long non-coding RNAs (m7G-related lncRNAs) are closely linked to the happenings and developments of various human beings’ cancers. However, the connection between m7G-related lncRNAs and glioma prognosis has not been investigated. We did this study to look for new potential biomarkers and construct an m7G-related lncRNA prognostic signature for glioma. We identified those lncRNAs associated with DEGs from glioma tissue sequences as m7G-related lncRNAs. First, we used Pearson’s correlation analysis to identify 28 DEGs by glioma and normal brain tissue gene sequences and predicated 657 m7G-related lncRNAs. Then, eight lncRNAs associated with prognosis were obtained and used to construct the m7G risk score model by lasso and Cox regression analysis methods. Furthermore, we used Kaplan-Meier analysis, time-dependent ROC, principal component analysis, clinical variables, independent prognostic analysis, nomograms, calibration curves, and expression levels of lncRNAs to determine the model’s accuracy. Importantly, we validated the model with external and internal validation methods and found it has strong predictive power. Finally, we performed functional enrichment analysis (GSEA, aaGSEA enrichment analyses) and analyzed immune checkpoints, associated pathways, and drug sensitivity based on predictors. In conclusion, we successfully constructed the formula of m7G-related lncRNAs with powerful predictive functions. Our study provides instructional value for analyzing glioma pathogenesis and offers potential research targets for glioma treatment and scientific research.