QKI shuttles internal m7G-modified transcripts into stress granules and modulates mRNA metabolism.

N7-methylguanosine (m7G) modification, routinely occurring at mRNA 5' cap or within tRNAs/rRNAs, also exists internally in messenger RNAs (mRNAs). Although m7G-cap is essential for pre-mRNA processing and protein synthesis, the exact role of mRNA internal m7G modification remains elusive. Here, we report that mRNA internal m7G is selectively recognized by Quaking proteins (QKIs). By transcriptome-wide profiling/mapping of internal m7G methylome and QKI-binding sites, we identified more than 1,000 high-confidence m7G-modified and QKI-bound mRNA targets with a conserved "GANGAN (N = A/C/U/G)" motif. Strikingly, QKI7 interacts (via C terminus) with the stress granule (SG) core protein G3BP1 and shuttles internal m7G-modified transcripts into SGs to regulate mRNA stability and translation under stress conditions. Specifically, QKI7 attenuates the translation efficiency of essential genes in Hippo signaling pathways to sensitize cancer cells to chemotherapy. Collectively, we characterized QKIs as mRNA internal m7G-binding proteins that modulate target mRNA metabolism and cellular drug resistance.

METTL1-mediated m7G modification of Arg-TCT tRNA drives oncogenic transformation.

The emerging "epitranscriptomics" field is providing insights into the biological and pathological roles of different RNA modifications. The RNA methyltransferase METTL1 catalyzes N7-methylguanosine (m7G) modification of tRNAs. Here we find METTL1 is frequently amplified and overexpressed in cancers and is associated with poor patient survival. METTL1 depletion causes decreased abundance of m7G-modified tRNAs and altered cell cycle and inhibits oncogenicity. Conversely, METTL1 overexpression induces oncogenic cell transformation and cancer. Mechanistically, we find increased abundance of m7G-modified tRNAs, in particular Arg-TCT-4-1, and increased translation of mRNAs, including cell cycle regulators that are enriched in the corresponding AGA codon. Accordingly, Arg-TCT expression is elevated in many tumor types and is associated with patient survival, and strikingly, overexpression of this individual tRNA induces oncogenic transformation. Thus, METTL1-mediated tRNA modification drives oncogenic transformation through a remodeling of the mRNA "translatome" to increase expression of growth-promoting proteins and represents a promising anti-cancer target.

NELF Regulates a Promoter-Proximal Step Distinct from RNA Pol II Pause-Release.

RNA polymerase II (RNA Pol II) is generally paused at promoter-proximal regions in most metazoans, and based on in vitro studies, this function has been attributed to the negative elongation factor (NELF). Here, we show that upon rapid depletion of NELF, RNA Pol II fails to be released into gene bodies, stopping instead around the +1 nucleosomal dyad-associated region. The transition to the 2nd pause region is independent of positive transcription elongation factor P-TEFb. During the heat shock response, RNA Pol II is rapidly released from pausing at heat shock-induced genes, while most genes are paused and transcriptionally downregulated. Both of these aspects of the heat shock response remain intact upon NELF loss. We find that NELF depletion results in global loss of cap-binding complex from chromatin without global reduction of nascent transcript 5' cap stability. Thus, our studies implicate NELF functioning in early elongation complexes distinct from RNA Pol II pause-release.

Transcriptome-wide Mapping of Internal N7-Methylguanosine Methylome in Mammalian mRNA.

N7-methylguanosine (m7G) is a positively charged, essential modification at the 5' cap of eukaryotic mRNA, regulating mRNA export, translation, and splicing. m7G also occurs internally within tRNA and rRNA, but its existence and distribution within eukaryotic mRNA remain to be investigated. Here, we show the presence of internal m7G sites within mammalian mRNA. We then performed transcriptome-wide profiling of internal m7G methylome using m7G-MeRIP sequencing (MeRIP-seq). To map this modification at base resolution, we developed a chemical-assisted sequencing approach that selectively converts internal m7G sites into abasic sites, inducing misincorporation at these sites during reverse transcription. This base-resolution m7G-seq enabled transcriptome-wide mapping of m7G in human tRNA and mRNA, revealing distribution features of the internal m7G methylome in human cells. We also identified METTL1 as a methyltransferase that installs a subset of m7G within mRNA and showed that internal m7G methylation could affect mRNA translation.

RlmQ: a newly discovered rRNA modification enzyme bridging RNA modification and virulence traits in Staphylococcus aureus.

rRNA modifications play crucial roles in fine-tuning the delicate balance between translation speed and accuracy, yet the underlying mechanisms remain elusive. Comparative analyses of the rRNA modifications in taxonomically distant bacteria could help define their general, as well as species-specific, roles. In this study, we identified a new methyltransferase, RlmQ, in Staphylococcus aureus responsible for the Gram-positive specific m7G2601, which is not modified in Escherichia coli (G2574). We also demonstrate the absence of methylation on C1989, equivalent to E. coli C1962, which is methylated at position 5 by the Gram-negative specific RlmI methyltransferase, a paralog of RlmQ. Both modifications (S. aureus m7G2601 and E. coli m5C1962) are situated within the same tRNA accommodation corridor, hinting at a potential shared function in translation. Inactivation of S. aureus rlmQ causes the loss of methylation at G2601 and significantly impacts growth, cytotoxicity, and biofilm formation. These findings unravel the intricate connections between rRNA modifications, translation, and virulence in pathogenic Gram-positive bacteria.

The molecular language of RNA 5' ends: guardians of RNA identity and immunity.

RNA caps are deposited at the 5' end of RNA polymerase II transcripts. This modification regulates several steps of gene expression, in addition to marking transcripts as self to enable the innate immune system to distinguish them from uncapped foreign RNAs, including those derived from viruses. Specialized immune sensors, such as RIG-I and IFITs, trigger antiviral responses upon recognition of uncapped cytoplasmic transcripts. Interestingly, uncapped transcripts can also be produced by mammalian hosts. For instance, 5'-triphosphate RNAs are generated by RNA polymerase III transcription, including tRNAs, Alu RNAs, or vault RNAs. These RNAs have emerged as key players of innate immunity, as they can be recognized by the antiviral sensors. Mechanisms that regulate the presence of 5'-triphosphates, such as 5'-end dephosphorylation or RNA editing, prevent immune recognition of endogenous RNAs and excessive inflammation. Here, we provide a comprehensive overview of the complexity of RNA cap structures and 5'-triphosphate RNAs, highlighting their roles in transcript identity, immune surveillance, and disease.

Mettl1/Wdr4-Mediated m7G tRNA Methylome Is Required for Normal mRNA Translation and Embryonic Stem Cell Self-Renewal and Differentiation.

tRNAs are subject to numerous modifications, including methylation. Mutations in the human N7-methylguanosine (m7G) methyltransferase complex METTL1/WDR4 cause primordial dwarfism and brain malformation, yet the molecular and cellular function in mammals is not well understood. We developed m7G methylated tRNA immunoprecipitation sequencing (MeRIP-seq) and tRNA reduction and cleavage sequencing (TRAC-seq) to reveal the m7G tRNA methylome in mouse embryonic stem cells (mESCs). A subset of 22 tRNAs is modified at a "RAGGU" motif within the variable loop. We observe increased ribosome occupancy at the corresponding codons in Mettl1 knockout mESCs, implying widespread effects on tRNA function, ribosome pausing, and mRNA translation. Translation of cell cycle genes and those associated with brain abnormalities is particularly affected. Mettl1 or Wdr4 knockout mESCs display defective self-renewal and neural differentiation. Our study uncovers the complexity of the mammalian m7G tRNA methylome and highlights its essential role in ESCs with links to human disease.

HIV-1 hypermethylated guanosine cap licenses specialized translation unaffected by mTOR.

Appended to the 5' end of nascent RNA polymerase II transcripts is 7-methyl guanosine (m7G-cap) that engages nuclear cap-binding complex (CBC) to facilitate messenger RNA (mRNA) maturation. Mature mRNAs exchange CBC for eIF4E, the rate-limiting translation factor that is controlled through mTOR. Experiments in immune cells have now documented HIV-1 incompletely processed transcripts exhibited hypermethylated m7G-cap and that the down-regulation of the trimethylguanosine synthetase-1-reduced HIV-1 infectivity and virion protein synthesis by several orders of magnitude. HIV-1 cap hypermethylation required nuclear RNA helicase A (RHA)/DHX9 interaction with the shape of the 5' untranslated region (UTR) primer binding site (PBS) segment. Down-regulation of RHA or the anomalous shape of the PBS segment abrogated hypermethylated caps and derepressed eIF4E binding for virion protein translation during global down-regulation of host translation. mTOR inhibition was detrimental to HIV-1 proliferation and attenuated Tat, Rev, and Nef synthesis. This study identified mutually exclusive translation pathways and the calibration of virion structural/accessory protein synthesis with de novo synthesis of the viral regulatory proteins. The hypermethylation of select, viral mRNA resulted in CBC exchange to heterodimeric CBP80/NCBP3 that expanded the functional capacity of HIV-1 in immune cells.

Characterization of the m6A/m1A/m5C/m7G-related regulators on the prognosis and immune microenvironment of glioma by integrated analysis of scRNA-seq and bulk RNA-seq data.

BACKGROUND: Proliferation, metabolism, tumor occurrence and development in gliomas are greatly influenced by RNA modifications. However, no research has integrated the four RNA methylation regulators of m6A, m1A, m5C and m7G in gliomas to analyze their relationship with glioma prognosis and intratumoral heterogeneity. METHODS: Based on three in-house single-cell RNA-sequencing (scRNA-seq) data, the glioma heterogeneity and characteristics of m6A/m1A/m5C/m7G-related regulators were elucidated. Based on publicly available bulk RNA-sequencing (RNA-seq) data, a risk-score system for predicting the overall survival (OS) for gliomas was established by three machine learning methods and multivariate Cox regression analysis, and validated in an independent cohort. RESULTS: Seven cell types were identified in gliomas by three scRNA-seq data, and 22 m6A/m1A/m5C/m7G-related regulators among the marker genes of different cell subtypes were discovered. Three m6A/m1A/m5C/m7G-related regulators were selected to construct prognostic risk-score model, including EIFA, NSUN6 and TET1. The high-risk patients showed higher immune checkpoint expression, higher tumor microenvironment scores, as well as higher tumor mutation burden and poorer prognosis compared with low-risk patients. Additionally, the area under the curve values of the risk score and nomogram were 0.833 and 0.922 for 3 year survival and 0.759 and 0.885 for 5 year survival for gliomas. EIF3A was significantly highly expressed in glioma tissues in our in-house RNA-sequencing data (p < 0.05). CONCLUSION: These findings may contribute to further understanding of the role of m6A/m1A/m5C/m7G-related regulators in gliomas, and provide novel and reliable biomarkers for gliomas prognosis and treatment.

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.

Targeting N7-methylguanosine tRNA modification blocks hepatocellular carcinoma metastasis after insufficient radiofrequency ablation.

Radiofrequency heat ablation is an ideal radical treatment for hepatocellular carcinoma (HCC). However, insufficient radiofrequency ablation (IRFA) could lead to high recurrence of HCC. N7-methylguanosine (m7G) on tRNAs, an evolutionally conservative modification in mammals and yeast, modulates heat stress responses and tumor progression, while its function in HCC recurrence after IRFA remains unknown. Here, we found that IRFA significantly upregulates the level of m7G tRNA modification and its methyltransferase complex components METTL1/WDR4 in multiple systems including HCC patient-derived xenograft (PDX) mouse, patients' HCC tissues, sublethal-heat-treated models of HCC cell lines, and organoids. Functionally, gain-/loss-of-function assays showed that METTL1-mediated m7G tRNA modification promotes HCC metastasis under sublethal heat exposure both in vitro and in vivo. Mechanistically, we found that METTL1 and m7G tRNA modification enhance the translation of SLUG/SNAIL in a codon frequency-dependent manner under sublethal heat stress. Overexpression of SLUG/SNAIL rescued the malignant potency of METTL1 knockdown HCC cells after sublethal heat exposure. Our study uncovers the key functions of m7G tRNA modification in heat stress responses and HCC recurrence after IRFA, providing molecular basis for targeting METTL1-m7G-SLUG/SNAIL axis to prevent HCC metastasis after radiofrequency heat ablation treatment.

Tumor-associated characteristics and immune dysregulation in nasopharyngeal carcinoma under the regulation of m7G-related tumor microenvironment cells.

BACKGROUND: Nasopharyngeal carcinoma (NPC) is a type of malignant tumor with high morbidity. Aberrant levels of N7-methylguanosine (m7G) are closely associated with tumor progression. However, the characteristics of the tumor microenvironment (TME) in NPC associated with m7G modification remain unclear. METHODS: A total of 68,795 single cells from single-cell RNA sequencing data derived from 11 NPC tumor samples and 3 nasopharyngeal lymphatic hyperplasia (NLH) samples were clustered using a nonnegative matrix factorization algorithm according to 61 m7G RNA modification regulators. RESULTS: The m7G regulators were found differential expression in the TME cells of NPC, and most m7G-related immune cell clusters in NPC tissues had a higher abundance compared to non-NPC tissues. Specifically, m7G scores in the CD4+ and CD8+ T cell clusters were significantly lower in NPC than in NLH. T cell clusters differentially expressed immune co-stimulators and co-inhibitors. Macrophage clusters differentially expressed EIF4A1, and high EIF4A1 expression was associated with poor survival in patients with head and neck squamous carcinoma. EIF4A1 was upregulated in NPC tissues compared to the non-NPC tissues and mainly expressed in CD86+ macrophages. Moreover, B cell clusters exhibited tumor biological characteristics under the regulation of m7G-related genes in NPC. The fibroblast clusters interacted with the above immune cell clusters and enriched tumor biological pathways, such as FGER2 signaling pathway. Importantly, there were correlations and interactions through various ligand-receptor links among epithelial cells and m7G-related TME cell clusters. CONCLUSION: Our study revealed tumor-associated characteristics and immune dysregulation in the NPC microenvironment under the regulation of m7G-related TME cells. These results demonstrated the underlying regulatory roles of m7G in NPC.

SPAAC-NAD-seq, a sensitive and accurate method to profile NAD+-capped transcripts.

Nicotinamide adenine diphosphate (NAD+) is a novel messenger RNA 5' cap in Escherichia coli, yeast, mammals, and Arabidopsis Transcriptome-wide identification of NAD+-capped RNAs (NAD-RNAs) was accomplished through NAD captureSeq, which combines chemoenzymatic RNA enrichment with high-throughput sequencing. NAD-RNAs are enzymatically converted to alkyne-RNAs that are then biotinylated using a copper-catalyzed azide-alkyne cycloaddition (CuAAC) reaction. Originally applied to E. coli RNA, which lacks the m7G cap, NAD captureSeq was then applied to eukaryotes without extensive verification of its specificity for NAD-RNAs vs. m7G-capped RNAs (m7G-RNAs). In addition, the Cu2+ ion in the CuAAC reaction causes RNA fragmentation, leading to greatly reduced yield and loss of full-length sequence information. We developed an NAD-RNA capture scheme utilizing the copper-free, strain-promoted azide-alkyne cycloaddition reaction (SPAAC). We examined the specificity of CuAAC and SPAAC reactions toward NAD-RNAs and m7G-RNAs and found that both prefer the former, but also act on the latter. We demonstrated that SPAAC-NAD sequencing (SPAAC-NAD-seq), when combined with immunodepletion of m7G-RNAs, enables NAD-RNA identification with accuracy and sensitivity, leading to the discovery of new NAD-RNA profiles in Arabidopsis Furthermore, SPAAC-NAD-seq retained full-length sequence information. Therefore, SPAAC-NAD-seq would enable specific and efficient discovery of NAD-RNAs in prokaryotes and, when combined with m7G-RNA depletion, in eukaryotes.

Toward the mechanism of eIF4F-mediated ribosomal attachment to mammalian capped mRNAs.

Ribosomal attachment to mammalian capped mRNAs is achieved through the cap-eukaryotic initiation factor 4E (eIF4E)-eIF4G-eIF3-40S chain of interactions, but the mechanism by which mRNA enters the mRNA-binding channel of the 40S subunit remains unknown. To investigate this process, we recapitulated initiation on capped mRNAs in vitro using a reconstituted translation system. Formation of initiation complexes at 5'-terminal AUGs was stimulated by the eIF4E-cap interaction and followed "the first AUG" rule, indicating that it did not occur by backward scanning. Initiation complexes formed even at the very 5' end of mRNA, implying that Met-tRNAi (Met) inspects mRNA from the first nucleotide and that initiation does not have a "blind spot." In assembled initiation complexes, the cap was no longer associated with eIF4E. Omission of eIF4A or disruption of eIF4E-eIF4G-eIF3 interactions converted eIF4E into a specific inhibitor of initiation on capped mRNAs. Taken together, these results are consistent with the model in which eIF4E-eIF4G-eIF3-40S interactions place eIF4E at the leading edge of the 40S subunit, and mRNA is threaded into the mRNA-binding channel such that Met-tRNAi (Met) can inspect it from the first nucleotide. Before entering, eIF4E likely dissociates from the cap to overcome steric hindrance. We also found that the m(7)G cap specifically interacts with eIF3l.

Novel m7G-related lncRNA signature for predicting overall survival in patients with gastric cancer.

Presenting with a poor prognosis, gastric cancer (GC) remains one of the leading causes of disease and death worldwide. Long non-coding RNAs (lncRNAs) regulate tumor formation and have been long used to predict tumor prognosis. N7-methylguanosine (m7G) is the most prevalent RNA modification. m7G-lncRNAs regulate GC onset and progression, but their precise mechanism in GC is unclear. The objective of this research was the development of a new m7G-related lncRNA signature as a biomarker for predicting GC survival rate and guiding treatment. The Cancer Genome Atlas database helped extract gene expression data and clinical information for GC. Pearson correlation analysis helped point out m7G-related lncRNAs. Univariate Cox analysis helped in identifying m7G-related lncRNA with predictive capability. The Lasso-Cox method helped point out seven lncRNAs for the purpose of establishing an m7G-related lncRNA prognostic signature (m7G-LPS), followed by the construction of a nomogram. Kaplan-Meier analysis, univariate and multivariate Cox regression analysis, calibration plot of the nomogram model, receiver operating characteristic curve and principal component analysis were utilized for the verification of the risk model's reliability. Furthermore, q-PCR helped verify the lncRNAs expression of m7G-LPS in-vitro. The study subjects were classified into high and low-risk groups based on the median value of the risk score. Gene enrichment analysis confirmed the constructed m7G-LPS' correlation with RNA transcription and translation and multiple immune-related pathways. Analysis of the clinicopathological features revealed more progressive features in the high-risk group. CIBERSORT analysis showed the involvement of m7G-LPS in immune cell infiltration. The risk score was correlated with immune checkpoint gene expression, immune cell and immune function score, immune cell infiltration, and chemotherapy drug sensitivity. Therefore, our study shows that m7G-LPS constructed using seven m7G-related lncRNAs can predict the survival time of GC patients and guide chemotherapy and immunotherapy regimens as biomarker.

Breast cancer prognosis and immunological characteristics are predicted using the m6A/m5C/m1A/m7G-related long noncoding RNA signature.

According to statistics, breast cancer (BC) has replaced lung cancer as the most common cancer in the world. Therefore, specific detection markers and therapeutic targets need to be explored as a way to improve the survival rate of BC patients. We first identified m6A/m5C/m1A/m7G-related long noncoding RNAs (MRlncRNAs) and developed a model of 16 MRlncRNAs. Kaplan-Meier survival analysis was applied to assess the prognostic power of the model, while univariate Cox analysis and multivariate Cox analysis were used to assess the prognostic value of the constructed model. Then, we constructed a nomogram to illustrate whether the predicted results were in good agreement with the actual outcomes. We tried to use the model to distinguish the difference in sensitivity to immunotherapy between the two groups and performed some analyses such as immune infiltration analysis, ssGSEA and IC50 prediction. To explore the novel anti-tumor drug response, we reclassified the patients into two clusters. Next, we assessed their response to clinical treatment by the R package pRRophetic, which is determined by the IC50 of each BC patient. We finally identified 11 MRlncRNAs and based on them, a risk model was constructed. In this model, we found good agreement between calibration plots and prognosis prediction. The AUC of ROC curves was 0.751, 0.734, and 0.769 for 1-year, 2-year, and 3-year overall survival (OS), respectively. The results showed that the IC50 was significantly different between the risk groups, suggesting that the risk groups can be used as a guide for systemic treatment. We regrouped patients into two clusters based on 11 MRlncRNAs expression. Next, we conducted immune scores for 2 clusters, which showed that cluster 1 had higher stromal scores, immune scores and higher estimated (microenvironment) scores, demonstrating that TME of cluster 1 was different from cluster 2. The results of this study support that MRlncRNAs can predict tumor prognosis and help differentiate patients with different sensitivities to immunotherapy as a basis for individualized treatment for BC patients.

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.

m7G-related genes-NCBP2 and EIF4E3 determine immune contexture in head and neck squamous cell carcinoma by regulating CCL4/CCL5 expression.

Aberrant N7 -methylguanosine (m7G) levels closely correlate with tumor genesis and progression. NCBP2 and EIF4E3 are two important m7G-related cap-binding genes. This study aimed to identify the relationship between the EIF4E3/NCBP2 function and immunological characteristics of head and neck squamous cell carcinoma (HNSCC). Hierarchical clustering was employed in classifying HNSCC patients into two groups based on the expressions of NCBP2 and EIF4E3. The differentially expressed genes were identified between the two groups, and GO functional enrichment was subsequently performed. Weighted gene co-expression network analysis was conducted to identify the hub genes related to EIF4E3/NCBP2 expression and immunity. The differential infiltration of immune cells and the response to immunotherapy were compared between the two groups. Single-cell sequence and trajectory analyses were performed to predict cell differentiation and display the expression of EIF4E3/NCBP2 in each state. In addition, quantitative real-time PCR, spatial transcriptome analysis, transwell assay, and western blotting were conducted to verify the biological function of EIF4E3/NCBP2. Here, group A showed a higher EIF4E3 expression and a lower NCBP2 expression, which had higher immune scores, proportion of most immune cells, immune activities, expression of immunomodulatory targets, and a better response to cancer immunotherapy. Besides, 56 hub molecules with notable immune regulation significance were identified. A risk model containing 17 hub genes and a prognostic nomogram was successfully established. Moreover, HNSCC tissues had a lower EIF4E3 expression and a higher NCBP2 expression than normal tissues. NCBP2 and EIF4E3 played a vital role in the differentiation of monocytes. Furthermore, the expression of CCL4/CCL5 can be regulated via EIF4E3 overexpression and NCBP2 knockdown. Collectively, NCBP2 and EIF4E3 can affect downstream gene expression, as well as immune contexture and response to immunotherapy, which could induce "cold-to-hot" tumor transformation in HNSCC patients.

Identification of m7G-related hub biomarkers and m7G regulator expression pattern in immune landscape during the progression of osteoarthritis.

BACKGROUND: Accumulating evidence has shown that aberrant N7-methylguanosine (m7G) RNA methylation played an important role in the occurrence and development of cancer. However, knowledge of m7G modifications in inflammatory diseases is limited. Osteoarthritis (OA) is the most common arthritic disease with poor prognosis. Our research aimed to identify m7G-related hub biomarkers and investigate m7G regulator expression pattern in immune landscape of OA patients. METHODS: Gene expression profiles and their clinical information were obtained from the Gene Expression Omnibus (GEO) database, and differential analysis of 14 m7G-related regulators between elective OA and normal samples was performed. M7G-related hub genes for OA were mined based on single-sample gene set enrichment analysis (ssGSEA) and the random forest (RF) algorithm, and qRT-PCR was performed to confirm the abnormal expression of hub genes. Enrichment, protein-protein interaction (PPI), transcription factor (TF)-gene interaction and microRNA (miRNA)-gene coregulatory analysis based on m7G hub genes were performed. Then we predicted several candidate drugs related to m7G hub genes using DSigDB database. Moreover, we comprehensively evaluated m7G methylation patterns in OA samples and systematically correlated these modification patterns with the characteristics of immune cell infiltration. The m7G score was generated to quantify m7G methylation patterns for individual OA patients by the application of principal component analysis (PCA) algorithm. RESULTS: We constructed an OA predictive model based on 4 m7G hub genes (SNUPN, METTL1, EIF4E2 and CYFIP1). Two m7G methylation patterns in OA were discovered to show distinct biological characteristics, and an m7G score were generated. M7G cluster A and a higher m7G score were found to be related to an inflamed phenotype. CONCLUSIONS: Our study was the first to comprehensively investigate the m7G methylation dysregulations in immune landscape during the progression of OA. These 4 m7G gene-related signatures can be used as novel OA biomarkers to predict the occurrence of OA. Evaluating the m7G methylation patterns of OA individuals will contribute to enhancing our cognition of immune infiltration characterization and guiding more effective immunotherapy strategies.

Identification of methyltransferase modification genes associated with prognosis and immune features of pancreatic adenocarcinoma.

BACKGROUND: Pancreatic adenocarcinoma (PAAD) is a malignant tumor with a high mortality rate. Methylation modifications acted a crucial role to affect cancer progression. The current study aimed to explore the potential role of methylase regulators in PAAD prognosis and immune microenvironment. METHODS: PubMed and TCGA databases were used to systematically analyze methylase regulators in PAAD. We identified three methylase clusters based on RNA methylase transcriptome data and obtained three gene clusters based on methylase modification-related differently expressed genes using principal component analysis (PCA) analysis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene Ontology (GO) biological processes were performed to explore the processes enriched in the different subgroups and single sample gene-set enrichment analysis (ssGSEA) was used to analyze the relationship between subgroups and immune infiltration in PAAD. RESULTS: We systematically screened 43 methylase regulators in PAAD samples and identified three methylase clusters with different clinical outcomes, as well as detected a significant relationship between methylase clusters and tumor immune infiltration. The top ten mutated genes include TP53, Kirsten rat sarcoma viral oncogene homolog (KRAS), titin gene (TTN), mucin 16 (MUC16), SMAD4, cyclin-dependent kinase inhibitor 2a (CDKN2A), Ryanodine receptor isoform-1 (RYR1), ring finger 43 (RNF43), protocadherin-15 (PCDH15), and AT-rich interacting domain-containing protein 1 A gene (ARID1A). CONCLUSION: The current study constructed an m6A/m5C/m1A/m7G modulator genes and explored methylase modification-related genes, which were related to the prognosis of PAAD patients and the immune checkpoint point cytotoxic T-lymphocyte associated protein 4 (CTLA4). These findings may provide prognostic predictors and direction for immunotherapy strategies for the treatment of PAAD.