RNA modifications and cancer

m6A- and m5C- modified lncRNAs orchestrate the prognosis in cutaneous melanoma and m6A- modified LINC00893 regulates cutaneous melanoma cell metastasis

BACKGROUND

As the most important modifications on the RNA level, N6-methyladenosine (m6A-) and 5-methylcytosine (m5C-) modification could have a direct influence on the RNAs. Long non-coding RNAs (lncRNAs) could also be modified by methylcytosine modification. Compared with mRNAs, the function of lncRNAs could be more potent to some extent in biological processes like tumorigenesis. Until now, rare reports have been done associated with cutaneous melanoma. Herein, we wonder if the m6A- and m5C- modified lncRNAs could influence the immune landscape and prognosis in melanoma, and we also want to find some lncRNAs which could directly affect the malignant behaviors of melanoma.

METHODS

Systematically, we explored the expression pattern of m6A- and m5C- modified lncRNAs in melanoma from datasets including UCSC Xena and NCBI GEO, and the prognostic lncRNAs were selected. Then, according to the expression pattern of lncRNAs, melanoma samples from these datasets were divided into several subtypes. Prognostic model, nomogram survival model, drug sensitivity, GO, and KEGG pathway analysis were performed. Furthermore, among several selected lncRNAs, we identified one lncRNA named LINC00893 and investigated its expression pattern and its biological function in melanoma cell lines.

RESULTS

We identified 27 m6A- and m5C- related lncRNAs which were significantly associated with survival, and we made a subtype analysis of melanoma samples based on these 27 lncRNAs. Among the two subtypes, we found differences of immune cells infiltration between these two subtypes. Then, LASSO algorithm was used to screen the optimized lncRNAs combination including ZNF252P-AS1, MIAT, FAM13A-AS1, LINC-PINT, LINC00893, AGAP2-AS1, OIP5-AS1, and SEMA6A-AS1. We also found that there was a significant correlation between the different risk groups predicted based on RS model and the actual prognosis. The nomogram survival model based on independent survival prognostic factors was also constructed. Besides, sensitivity to chemotherapeutic agents, GO and KEGG analysis were performed. In different risk groups, a total of 14 drug molecules with different distributions were obtained, which included AZD6482, AZD7762, AZD8055, camptothecin, dasatinib, erlotinib, gefitinib, gemcitabine, GSK269962A, nilotinib, rapamycin, and sorafenib. A total of 55 significantly related biological processes and 17 KEGG signaling pathways were screened. At last, we noticed that LINC00893 had a relatively lower expression in melanoma tissue and cell lines compared with adjacent tissues and epidermal melanocyte, and down-regulation of LINC00893 could promote the malignant behavior of melanoma cells in A875 and MV3. In these two melanoma cell lines, down-regulation of m6A-related molecules like YTHDF3 and METTL3 could promote the expression of LINC00893.

CONCLUSION

We made an analysis of m6A- and m5C- related lncRNAs in melanoma samples and a prediction of these lncRNAs' role in prognosis, tumor microenvironment, immune infiltration, and clinicopathological features. We also found that LINC00893, which is potentially regulated by m6A modification, could serve as a tumor-suppressor in melanoma and play an inhibitory role in melanoma metastasis.

N6-methyladenosine methylation regulates the tumor microenvironment of Epstein-Barr virus-associated gastric cancer

BACKGROUND

N6-methyladenosine (m6A) methylation modification exists in Epstein-Barr virus (EBV) primary infection, latency, and lytic reactivation. It also modifies EBV latent genes and lytic genes. EBV-associated gastric cancer (EBVaGC) is a distinctive molecular subtype of GC. We hypothesized EBV and m6A methylation regulators interact with each other in EBVaGC to differentiate it from other types of GC.

AIM

To investigate the mechanisms of m6A methylation regulators in EBVaGC to determine the differentiating factors from other types of GC.

METHODS

First, The Cancer Gene Atlas and Gene Expression Omnibus databases were used to analyze the expression pattern of m6A methylation regulators between EBVaGC and EBV-negative GC (EBVnGC). Second, we identified Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment of m6A-related differentially expressed genes. We quantified the relative abundance of immune cells and inflammatory factors in the tumor microenvironment (TME). Finally, cell counting kit-8 cell proliferation test, transwell test, and flow cytometry were used to verify the effect of insulin-like growth factor binding protein 1 (IGFBP1) in EBVaGC cell lines.

RESULTS

m6A methylation regulators were involved in the occurrence and development of EBVaGC. Compared with EBVnGC, the expression levels of m6A methylation regulators Wilms tumor 1-associated protein, RNA binding motif protein 15B, CBL proto-oncogene like 1, leucine rich pentatricopeptide repeat containing, heterogeneous nuclear ribonucleoprotein A2B1, IGFBP1, and insulin-like growth factor 2 binding protein 1 were significantly downregulated in EBVaGC (P < 0.05). The overall survival rate of EBVaGC patients with a lower expression level of IGFBP1 was significantly higher (P = 0.046). GO and KEGG functional enrichment analyses showed that the immunity pathways were significantly activated and rich in immune cell infiltration in EBVaGC. Compared with EBVnGC, the infiltration of activated CD4+ T cells, activated CD8+ T cells, monocytes, activated dendritic cells, and plasmacytoid dendritic cells were significantly upregulated in EBVaGC (P < 0.001). In EBVaGC, the expression level of proinflammatory factors interleukin (IL)-17, IL-21, and interferon-γ and immunosuppressive factor IL-10 were significantly increased (P < 0.05). In vitro experiments demonstrated that the expression level of IGFBP1 was significantly lower in an EBVaGC cell line (SNU719) than in an EBVnGC cell line (AGS) (P < 0.05). IGFBP1 overexpression significantly attenuated proliferation and migration and promoted the apoptosis levels in SNU719. Interfering IGFBP1 significantly promoted proliferation and migration and attenuated the apoptosis levels in AGS.

CONCLUSION

m6A regulators could remodel the TME of EBVaGC, which is classified as an immune-inflamed phenotype and referred to as a “hot” tumor. Among these regulators, we demonstrated that IGFBP1 affected proliferation, migration, and apoptosis.

N6-methyladenosine methylation regulates the tumor microenvironment of Epstein-Barr virus-associated gastric cancer

BACKGROUND

N6-methyladenosine (m6A) methylation modification exists in Epstein-Barr virus (EBV) primary infection, latency, and lytic reactivation. It also modifies EBV latent genes and lytic genes. EBV-associated gastric cancer (EBVaGC) is a distinctive molecular subtype of GC. We hypothesized EBV and m6A methylation regulators interact with each other in EBVaGC to differentiate it from other types of GC.

AIM

To investigate the mechanisms of m6A methylation regulators in EBVaGC to determine the differentiating factors from other types of GC.

METHODS

First, The Cancer Gene Atlas and Gene Expression Omnibus databases were used to analyze the expression pattern of m6A methylation regulators between EBVaGC and EBV-negative GC (EBVnGC). Second, we identified Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment of m6A-related differentially expressed genes. We quantified the relative abundance of immune cells and inflammatory factors in the tumor microenvironment (TME). Finally, cell counting kit-8 cell proliferation test, transwell test, and flow cytometry were used to verify the effect of insulin-like growth factor binding protein 1 (IGFBP1) in EBVaGC cell lines.

RESULTS

m6A methylation regulators were involved in the occurrence and development of EBVaGC. Compared with EBVnGC, the expression levels of m6A methylation regulators Wilms tumor 1-associated protein, RNA binding motif protein 15B, CBL proto-oncogene like 1, leucine rich pentatricopeptide repeat containing, heterogeneous nuclear ribonucleoprotein A2B1, IGFBP1, and insulin-like growth factor 2 binding protein 1 were significantly downregulated in EBVaGC (P < 0.05). The overall survival rate of EBVaGC patients with a lower expression level of IGFBP1 was significantly higher (P = 0.046). GO and KEGG functional enrichment analyses showed that the immunity pathways were significantly activated and rich in immune cell infiltration in EBVaGC. Compared with EBVnGC, the infiltration of activated CD4+ T cells, activated CD8+ T cells, monocytes, activated dendritic cells, and plasmacytoid dendritic cells were significantly upregulated in EBVaGC (P < 0.001). In EBVaGC, the expression level of proinflammatory factors interleukin (IL)-17, IL-21, and interferon-γ and immunosuppressive factor IL-10 were significantly increased (P < 0.05). In vitro experiments demonstrated that the expression level of IGFBP1 was significantly lower in an EBVaGC cell line (SNU719) than in an EBVnGC cell line (AGS) (P < 0.05). IGFBP1 overexpression significantly attenuated proliferation and migration and promoted the apoptosis levels in SNU719. Interfering IGFBP1 significantly promoted proliferation and migration and attenuated the apoptosis levels in AGS.

CONCLUSION

m6A regulators could remodel the TME of EBVaGC, which is classified as an immune-inflamed phenotype and referred to as a “hot” tumor. Among these regulators, we demonstrated that IGFBP1 affected proliferation, migration, and apoptosis.

YTHDF1 promotes the osteolytic bone metastasis of breast cancer via inducing EZH2 and CDH11 translation

Bone metastasis is common in breast cancer and more effective therapies are required, however, its molecular mechanism is poorly understood. Additionally, the role of the m6A reader YTHDF1 in bone metastasis of breast cancer has not been reported. Here, we reveal that the increased expression of YTHDF1 is clinically correlated with breast cancer bone metastases. YTHDF1 promotes migration, invasion, and osteoblast adhesion and induces osteoclast differentiation of cancer cells in vitro and vivo. Mechanically, RNA-seq, MeRIP-seq and RIP-seq analysis, and molecular biology experiments demonstrate that YTHDF1 translationally enhances EZH2 and CDH11 expression by reading m6A-enriched sites of their transcripts. Moreover, adeno-associated virus (AAV) was used to deliver shYTHDF1 (shYTHDF1-AAV) in intratibial injection models, eliciting a significant suppressive effect on breast cancer bone metastatic formation and osteolytic destruction. Overall, we uncovered that YTHDF1 promotes osteolytic bone metastases of breast cancer by inducing EZH2 and CDH11 translation.

Prognostic signature and immune efficacy of m1A-, m5C-, m6A-, m7G-, and DNA methylation-related regulators in hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC) is the main type of primary liver cancer, and its related death ranks third worldwide. The curative methods and progress prediction markers of HCC are not sufficient enough. Nevertheless, little progress has been made in the signature of m1A-, m5C-, m6A-, m7G-, and DNA methylation of HCC. Results: We calibrated a risk gene signature model that can be used to categorize HCC patients based on univariate, multivariate, and LASSO Cox regression analysis. This gene signature classified the patients into high- and low-risk subgroups. Patients in the high-risk group showed significantly reduced overall survival (OS) compared with patients in the low-risk group. The gene set variation analysis (GSVA), immune infiltration, and immunotherapy response were analyzed. The results demonstrated that an immunosuppressive environment was exited and the high-risk group had higher sensitivity to 5-fluorouracil, cisplatin, sorafenib, tamoxifen, and epirubicin. These results indicated personalized therapy should be taken into consideration. Conclusions: Our findings enriched our understanding of the molecular heterogeneity, tumor microenvironment (TME), and drug susceptibility of HCC. m1A-, m5C-, m6A-, m7G-, and DNA methylation-related regulators may be promising biomarkers for future research.

Breast cancer clinical outcomes and tumor immune microenvironment: cross-dialogue of multiple epigenetic modification profiles

The discovery of RNA methylation alterations associated with cancer holds promise for their utilization as potential biomarkers in cancer diagnosis, prognosis, and prediction. RNA methylation has been found to impact the immunological microenvironment of tumors, but the specific role of methylation-related genes (MRGs), particularly in breast cancer (BC), the most common cancer among women globally, within the tumor microenvironment remains unknown. In this study, we obtained data from TCGA and GEO databases to investigate the expression patterns of MRGs in both genomic and transcriptional domains in BC. By analyzing the data, we identified two distinct genetic groupings that were correlated with clinicopathological characteristics, prognosis, degree of TME cell infiltration, and other abnormalities in MRGs among patients. Subsequently, an MRG model was developed to predict overall survival (OS) and its accuracy was evaluated in BC patients. Additionally, a highly precise nomogram was created to enhance the practical usability of the MRG model. In low-risk groups, we observed lower TBM values and higher TIDE scores. We further explored how MRGs influence a patient's prognosis, clinically significant characteristics, response to therapy, and the TME. These risk signatures have the potential to improve treatment strategies for BC patients and could be applied in future clinical settings. Moreover, they may also be utilized to determine prognosis and biological features in these patients.

Cigarette smoking combined with genetic variation regulates the m6A methylation of CRNKL1 and is associated with bladder cancer risk

Cigarette smoking was known to accelerate the occurrence and development of bladder cancer by regulating RNA modification. However, the association between the combination of cigarette smoking and RNA modification-related single nucleotide polymorphisms (RNAm-SNPs) and bladder cancer risk remains unclear. In this study, 1681 participants, including 580 cases and 1101 controls, were recruited for genetic association analysis. In total, 1 287 990 RNAm-SNPs involving nine RNA modifications (m6A, m1A, m6Am, 2'-O-Me, m5C, m7G, A-to-I, m5U, and pseudouridine modification) were obtained from the RMVar database. The interactive effect of cigarette smoking and RNAm-SNPs on bladder cancer risk was assessed through joint analysis. The susceptibility analysis revealed that 89 RNAm-SNPs involving m6A, m1A, and A-to-I modifications were associated with bladder cancer risk. Among them, m6A-related rs2273058 in CRNKL1 was associated with bladder cancer risk (odds ratios (OR) = 1.35, padj = 1.78 × 10-4), and CRNKL1 expression was increased in bladder cancer patients (p = 0.035). Cigarette smoking combined with the A allele of rs2273058 increased bladder cancer risk compared with nonsmokers with the G allele of rs2273058 (OR = 2.40, padj = 3.11 × 10-9). Mechanistically, the A allele of rs2273058 endowed CRNKL1 with an additional m6A motif, facilitating recognition by m6A reader IGF2BP1, thereby promoting CRNKL1 expression under cigarette smoking (r = 0.142, p = 0.017). Moreover, elevated CRNKL1 expression may accelerate cell cycle and proliferation, thereby increasing bladder cancer risk. In summary, our study demonstrated that cigarette smoking combined with RNAm-SNPs contributes to bladder cancer risk, which provides a potential target for bladder cancer prevention.

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

Dissecting the oncogenic properties of essential RNA-modifying enzymes: a focus on NAT10

Chemical modifications of ribonucleotides significantly alter the physicochemical properties and functions of RNA. Initially perceived as static and essential marks in ribosomal RNA (rRNA) and transfer RNA (tRNA), recent discoveries unveiled a dynamic landscape of RNA modifications in messenger RNA (mRNA) and other regulatory RNAs. These findings spurred extensive efforts to map the distribution and function of RNA modifications, aiming to elucidate their distribution and functional significance in normal cellular homeostasis and pathological states. Significant dysregulation of RNA modifications is extensively documented in cancers, accentuating the potential of RNA-modifying enzymes as therapeutic targets. However, the essential role of several RNA-modifying enzymes in normal physiological functions raises concerns about potential side effects. A notable example is N-acetyltransferase 10 (NAT10), which is responsible for acetylating cytidines in RNA. While emerging evidence positions NAT10 as an oncogenic factor and a potential target in various cancer types, its essential role in normal cellular processes complicates the development of targeted therapies. This review aims to comprehensively analyze the essential and oncogenic properties of NAT10. We discuss its crucial role in normal cell biology and aging alongside its contribution to cancer development and progression. We advocate for agnostic approaches to disentangling the intertwined essential and oncogenic functions of RNA-modifying enzymes. Such approaches are crucial for understanding the full spectrum of RNA-modifying enzymes and imperative for designing effective and safe therapeutic strategies.

METTL1 mediated tRNA m7G modification promotes leukaemogenesis of AML via tRNA regulated translational control

BACKGROUND

RNA modifications have been proven to play fundamental roles in regulating cellular biology process. Recently, maladjusted N7-methylguanosine (m7G) modification and its modifiers METTL1/WDR4 have been confirmed an oncogene role in multiple cancers. However, the functions and molecular mechanisms of METTL1/WDR4 in acute myeloid leukemia (AML) remain to be determined.

METHODS

METTL1/WDR4 expression levels were quantified using qRT-PCR, western blot analysis on AML clinical samples, and bioinformatics analysis on publicly available AML datasets. CCK-8 assays and cell count assays were performed to determine cell proliferation. Flow cytometry assays were conducted to assess cell cycle and apoptosis rates. Multiple techniques were used for mechanism studies in vitro assays, such as northern blotting, liquid chromatography-coupled mass spectrometry (LC-MS/MS), tRNA stability analysis, transcriptome sequencing, small non-coding RNA sequencing, quantitative proteomics, and protein synthesis measurements.

RESULTS

METTL1/WDR4 are significantly elevated in AML patients and associated with poor prognosis. METTL1 knockdown resulted in reduced cell proliferation and increased apoptosis in AML cells. Mechanically, METTL1 knockdown leads to significant decrease of m7G modification abundance on tRNA, which further destabilizes tRNAs and facilitates the biogenesis of tsRNAs in AML cells. In addition, profiling of nascent proteins revealed that METTL1 knockdown and transfection of total tRNAs that were isolated from METTL1 knockdown AML cells decreased global translation efficiency in AML cells.

CONCLUSIONS

Taken together, our study demonstrates the important role of METTL1/WDR4 in AML leukaemogenesis, which provides a promising target candidate for AML therapy.

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.

An ensemble learning-based feature selection algorithm for identification of biomarkers of renal cell carcinoma

Feature selection plays a crucial role in classification tasks as part of the data preprocessing process. Effective feature selection can improve the robustness and interpretability of learning algorithms, and accelerate model learning. However, traditional statistical methods for feature selection are no longer practical in the context of high-dimensional data due to the computationally complex. Ensemble learning, a prominent learning method in machine learning, has demonstrated exceptional performance, particularly in classification problems. To address the issue, we propose a three-stage feature selection algorithm framework for high-dimensional data based on ensemble learning (EFS-GINI). Firstly, highly linearly correlated features are eliminated using the Spearman coefficient. Then, a feature selector based on the F-test is employed for the first stage selection. For the second stage, four feature subsets are formed using mutual information (MI), ReliefF, SURF, and SURF* filters in parallel. The third stage involves feature selection using a combinator based on GINI coefficient. Finally, a soft voting approach is proposed to employ for classification, including decision tree, naive Bayes, support vector machine (SVM), k-nearest neighbors (KNN) and random forest classifiers. To demonstrate the effectiveness and efficiency of the proposed algorithm, eight high-dimensional datasets are used and five feature selection methods are employed to compare with our proposed algorithm. Experimental results show that our method effectively enhances the accuracy and speed of feature selection. Moreover, to explore the biological significance of the proposed algorithm, we apply it on the renal cell carcinoma dataset GSE40435 from the Gene Expression Omnibus database. Two feature genes, NOP2 and NSUN5, are selected by our proposed algorithm. They are directly involved in regulating m5c RNA modification, which reveals the biological importance of EFS-GINI. Through bioinformatics analysis, we shows that m5C-related genes play an important role in the occurrence and progression of renal cell carcinoma, and are expected to become an important marker to predict the prognosis of patients.

A review on the role of RNA methylation in aging-related diseases

Senescence is the underlying mechanism of organism aging and is robustly regulated at the post-transcriptional level. This regulation involves the chemical modifications, of which the RNA methylation is the most common. Recently, a rapidly growing number of studies have demonstrated that methylation is relevant to aging and aging-associated diseases. Owing to the rapid development of detection methods, the understanding on RNA methylation has gone deeper. In this review, we summarize the current understanding on the influence of RNA modification on cellular senescence, with a focus on mRNA methylation in aging-related diseases, and discuss the emerging potential of RNA modification in diagnosis and therapy.

LAMA3 Promotes Tumorigenesis of Oral Squamous Cell Carcinoma by METTL3-Mediated N6-Methyladenosine Modification

Laminin subunit alpha 3 (LAMA3) is a cancer regulator. However, its effects and regulatory pathways in oral squamous cell carcinoma (OSCC) progression remain unknown. This research aimed to determine the influence of LAMA3 regulation via methyltransferase-like 3 (METTL3) on OSCC progression. Using quantitative real-time polymerase chain reaction and bioinformatics analysis, the expression levels of LAMA3 and METTL3 in OSCC tissues were examined. The functional roles of LAMA3 and METTL3 were analyzed using cell functional experiments. Using methylated RNA immunoprecipitation and mRNA stability assays, LAMA3 and METTL3 regulation was investigated. In OSCC tissues, LAMA3 was upregulated. LAMA3 inhibition hampered OSCC cell proliferation, invasion, and migration while its overexpression facilitated OSCC cell progression. METTL3 serves as a crucial upstream regulator of LAMA3 in OSCC and upregulates LAMA3 expression via an m6A-dependent mechanism. The low METTL3 expression partially restored the enhanced malignant phenotype induced by LAMA3 overexpression. Our findings indicate that METTL3 and LAMA3 act as pro-oncogenic factors in OSCC, with METTL3 promoting OSCC malignancy via m6A modification-dependent stabilization of LAMA3 transcripts, representing a novel regulatory mechanism in OSCC.

An m1A/m6A/m5C-associated long non-coding RNA signature: Prognostic and immunotherapeutic insights into cervical cancer

BACKGROUND

Cervical cancer (CC) remains a significant clinical challenge, even though its fatality rate has been declining in recent years. Particularly in developing countries, the prognosis for CC patients continues to be suboptimal despite numerous therapeutic advances.

METHODS

Using The Cancer Genome Atlas database, we extracted CC-related data. From this, 52 methylation-related genes (MRGs) were identified, leading to the selection of a 10 long non-coding RNA (lncRNA) signature co-expressed with these MRGs. R programming was employed to filter out the methylation-associated lncRNAs. Through univariate, least absolute shrinkage and selection operator (i.e. LASSO) and multivariate Cox regression analysis, an MRG-associated lncRNA model was constructed. The established risk model was further assessed via the Kaplan-Meier method, principal component analysis, functional enrichment annotation and a nomogram. Furthermore, we explored the potential of this model with respect to guiding immune therapeutic interventions and predicting drug sensitivities.

RESULTS

The derived 10-lncRNA signature, linked with MRGs, emerged as an independent prognostic factor. Segmenting patients based on their immunotherapy responses allowed for enhanced differentiation between patient subsets. Lastly, we highlighted potential compounds for distinguishing CC subtypes.

CONCLUSIONS

The risk model, associated with MRG-linked lncRNA, holds promise in forecasting clinical outcomes and gauging the efficacy of immunotherapies for CC patients.

m5 C regulator-mediated methylation modification patterns and tumor microenvironment infiltration characteristics in acute myeloid leukemia

BACKGROUND

Recently, many studies have been conducted to examine immune response modification at epigenetic level, but the candidate effect of RNA 5-methylcytosine (m5 C) modification on tumor microenvironment (TME) of acute myeloid leukemia (AML) is still unknown at present.

METHODS

We assessed the patterns of m5 C modification among 417 AML cases by using nine m5 C regulators. Thereafter, we associated those identified modification patterns with TME cell infiltration features. Additionally, stepwise regression and LASSO Cox regression analyses were conducted for quantifying patterns of m5 C modification among AML cases to establish the m5 C-score. Meanwhile, we validated the expression of genes in the m5C-score model by qRT-PCR. Finally, the present work analyzed the association between m5 C-score and AML clinical characteristics and prognostic outcomes.

RESULTS

In total, three different patterns of m5 C modification (m5 C-clusters) were identified, and highly differentiated TME cell infiltration features were also identified. On this basis, evaluating patterns of m5 C modification in single cancer samples was important for evaluating the immune/stromal activities in TME and for predicting prognosis. In addition, the m5 C-score was established, which showed a close relation with the overall survival (OS) of test and training set samples. Moreover, multivariate Cox analysis suggested that our constructed m5 C-score served as the independent predicting factor for the prognosis of AML (hazard ratio = 1.57, 95% confidence interval = 1.38-1.79, p < 1e-5 ).

CONCLUSIONS

This study shows that m5 C modification may be one of the key roles in the formation of diversity and complexity of TME. Meanwhile, assessing the patterns of m5 C modification among individual cancer samples is of great importance, which provides insights into cell infiltration features within TME, thereby helping to develop relevant immunotherapy and predict patient prognostic outcomes.

Function of m5C RNA methyltransferase NOP2 in high-grade serous ovarian cancer

RNA methyltransferase nucleolar protein p120 (NOP2), commonly referred to as NOP2/Sun RNA methyltransferase family member 1 (NSUN1), is involved in cell proliferation and is highly expressed in various cancers. However, its role in high-grade serous ovarian cancer (HGSOC) remains unclear. Our study investigated the expression of NOP2 in HGSOC tissues and normal fimbria tissues, and found that NOP2 was significantly upregulated in HGSOC tissues. Our experiments showed that NOP2 overexpression promoted cell proliferation in vivo and in vitro and increased the migration and invasion ability of HGSOC cells in vitro. Furthermore, we identified Rap guanine nucleotide exchange factor 4 (RAPGEF4) as a potential downstream target of NOP2 in HGSOC. Finally, our findings suggest that the regulation of NOP2 and RAPGEF4 may depend on m5C methylation levels.

lncRNA POU6F2-AS1 Regulated by KIAA1429 Contributes to Colorectal Cancer Progression in an m6A Modification Manner

Long non-coding RNAs (lncRNAs) are participated in tumourigenesis, including colorectal cancer (CRC). However, the effects and mechanisms of lncRNA POU6F2-AS1 in CRC have not been investigated. KIAA1429 act as a member of N6-methyladenosine (m6A) modification, has been knew as an oncogenic factor in various cancer containing CRC. We focus to investigate the regulation effect of lncRNA POU6F2-AS1, and the mechanism among lncRNA POU6F2-AS1 and KIAA1429 in CRC. The lncRNA POU6F2-AS1 and KIAA1429 levels in CRC tissue samples as well as cells were clarified by qRT-PCR, and their relationship was predicted by bioinformatics, MeRIP and Pearson analysis. Cell survival, migration and invasion were analyzed via EdU, wound healing and Transwell assays after lncRNA POU6F2-AS1 was down-regulated and KIAA1429 was up-regulated. LncRNA POU6F2-AS1 and KIAA1429 were enriched in CRC tissue samples. LncRNA POU6F2-AS1 silencing suppressed CRC cell survival, migration, and invasion, and KIAA1429 overexpression facilitated CRC cell malignancy. KIAA1429 promoted lncRNA POU6F2-AS1 expression via m6A modification. Furthermore, KIAA1429 upregulation reversed the inhibitory effect of lncRNA POU6F2-AS1 interference on the malignant behavior of CRC cells. lncRNA POU6F2-AS1 was modulated by KIAA1429 in the form of m6A modification to regulate the malignant phenotype of CRC, which may provide new insights into the potential application of KIAA1429-m6A-lncRNA POU6F2-AS1-based CRC therapy.

The complex nature of lncRNA-mediated chromatin dynamics in multiple myeloma

Extensive genome-wide sequencing efforts have unveiled the intricate regulatory potential of long non-protein coding RNAs (lncRNAs) within the domain of haematological malignancies. Notably, lncRNAs have been found to directly modulate chromatin architecture, thereby impacting gene expression and disease progression by interacting with DNA, RNA, and proteins in a tissue- or condition-specific manner. Furthermore, recent studies have highlighted the intricate epigenetic control of lncRNAs in cancer. Consequently, this provides a rationale to explore the possibility of therapeutically targeting lncRNAs themselves or the epigenetic mechanisms that govern their activity. Within the scope of this review, we will assess the current state of knowledge regarding the epigenetic regulation of lncRNAs and how, in turn, lncRNAs contribute to chromatin remodelling in the context of multiple myeloma.

m6A-mediated lncRNA NEAT1 plays an oncogenic role in non-small cell lung cancer by upregulating the HMGA1 expression through binding miR-361-3p

BACKGROUND

Lung cancer is the most common primary malignant tumor of the lung, and 85% of lung cancer is non-small cell lung cancer (NSCLC). The N6-methyladenosine (m6A) and long noncoding RNAs (lncRNAs) have been widely reported to participate in the development of non-small cell lung cancer.

OBJECTIVE

However, the potential molecular mechanisms of m6A-regulated lncRNAs in NSCLC still need further investigation.

METHODS

The expression levels and the role of lncRNA NEAT1 in NSCLC tissues or cells were measured by RT-qPCR, Western blot, cell counting kit 8 (CCK-8), flow cytometry assay. RNA immunoprecipitation (RIP) was used to measure the levels of m6A modification of NEAT1. Bioinformatics analysis and dual-luciferase reporter gene assay were detected the relationship between miR-361-3p and NEAT1/HMGA1. Mouse xenograft tumor models were established to confirm the effects of lncRNA NEAT1 in vivo.

RESULTS

In this study, we verified whether m6A-modified lncRNA nuclear enriched abundant transcript 1 (NEAT1) is involved in NSCLC progression via miR-361-3p/HMGA1 axis. Firstly, we found that lncRNA NEAT1 was upregulated in NSCLC, and was associated with a poor survival in NSCLC patients. Methyltransferase like 3 (METTL3)-mediated m6A modification stabilized and upregulated NEAT1 expression. Next, function experiment indicated that depletion of METTL3 and NEAT1 induced cell apoptosis and inhibited cell proliferation, epithelial-mesenchymal transition (EMT). Likewise, in vivo experiments further supported the oncogenic role of NEAT1 in NSCLC. In addition, the molecular mechanism was uncovered in our study, and we found that lncRNA NEAT1 promoted the expression of high-mobility group AT-hook 1 (HMGA1) by sponging miR-361-3p and then promoted tumorigenesis of NSCLC.

CONCLUSION

In conclusion, our findings demonstrated that METTL3-mediated m6A modification accelerated NSCLC progression by regulating the NEAT1/miR-361-3p/HMGA1 axis, which provides important targets for the treatment of NSCLC.

CDK13 promotes lipid deposition and prostate cancer progression by stimulating NSUN5-mediated m5C modification of ACC1 mRNA

Cyclin-dependent kinases (CDKs) regulate cell cycle progression and the transcription of a number of genes, including lipid metabolism-related genes, and aberrant lipid metabolism is involved in prostate carcinogenesis. Previous studies have shown that CDK13 expression is upregulated and fatty acid synthesis is increased in prostate cancer (PCa). However, the molecular mechanisms linking CDK13 upregulation and aberrant lipid metabolism in PCa cells remain largely unknown. Here, we showed that upregulation of CDK13 in PCa cells increases the fatty acyl chains and lipid classes, leading to lipid deposition in the cells, which is positively correlated with the expression of acetyl-CoA carboxylase (ACC1), the first rate-limiting enzyme in fatty acid synthesis. Gain- and loss-of-function studies showed that ACC1 mediates CDK13-induced lipid accumulation and PCa progression by enhancing lipid synthesis. Mechanistically, CDK13 interacts with RNA-methyltransferase NSUN5 to promote its phosphorylation at Ser327. In turn, phosphorylated NSUN5 catalyzes the m5C modification of ACC1 mRNA, and then the m5C-modified ACC1 mRNA binds to ALYREF to enhance its stability and nuclear export, thereby contributing to an increase in ACC1 expression and lipid deposition in PCa cells. Overall, our results disclose a novel function of CDK13 in regulating the ACC1 expression and identify a previously unrecognized CDK13/NSUN5/ACC1 pathway that mediates fatty acid synthesis and lipid accumulation in PCa cells, and targeting this newly identified pathway may be a novel therapeutic option for the treatment of PCa.

Neo-Darwinian Principles Exemplified in Cancer Genomics

Within the last two decades, the advent of next-generation sequencing accompanied by single-cell technologies has enabled cancer researchers to study in detail mutations and other genetic aberrations that transpire during transformation of cells to a neoplastic state. This article covers the insights gained through these extensive studies where neo-Darwinian principles can be inferred to play roles throughout neoplastic transformation. The cells promoted during cancer development exhibit cancer hallmarks combined with the related enabling characteristics as outlined by Hanahan and Weinberg, analogous to natural selection and survival of the fittest. Selection of driver mutations that inactivate proteins encoded by tumor suppressor genes differs in profound ways from mutations that activate tumor promoter proteins. In most cases, the later stages of cancer development are characterized by sudden, extensive damage to chromosomes in a process that is not Darwinian in nature. Nevertheless, cells that survive these cataclysmic events remain subject to Darwinian selection promoting clones exhibiting the greatest rates of progression. Duplications of chromosomal segments containing oncogenes, deletions of segments harboring tumor suppressor genes, or distinctive chromosomal rearrangements are often found in cells progressing into later stages of cancer. In summary, the technological developments in genome sequencing since the start of the century have given us clear insights into genomic alterations promoting tumor progression where neo-Darwinian mechanisms of clonal selection can be inferred to play a primary role.

A Hitchhiker's guide to RNA-RNA structure and interaction prediction tools

RNA biology has risen to prominence after a remarkable discovery of diverse functions of noncoding RNA (ncRNA). Most untranslated transcripts often exert their regulatory functions into RNA-RNA complexes via base pairing with complementary sequences in other RNAs. An interplay between RNAs is essential, as it possesses various functional roles in human cells, including genetic translation, RNA splicing, editing, ribosomal RNA maturation, RNA degradation and the regulation of metabolic pathways/riboswitches. Moreover, the pervasive transcription of the human genome allows for the discovery of novel genomic functions via RNA interactome investigation. The advancement of experimental procedures has resulted in an explosion of documented data, necessitating the development of efficient and precise computational tools and algorithms. This review provides an extensive update on RNA-RNA interaction (RRI) analysis via thermodynamic- and comparative-based RNA secondary structure prediction (RSP) and RNA-RNA interaction prediction (RIP) tools and their general functions. We also highlighted the current knowledge of RRIs and the limitations of RNA interactome mapping via experimental data. Then, the gap between RSP and RIP, the importance of RNA homologues, the relationship between pseudoknots, and RNA folding thermodynamics are discussed. It is hoped that these emerging prediction tools will deepen the understanding of RNA-associated interactions in human diseases and hasten treatment processes.

Insights into the regulatory role of RNA methylation modifications in glioma

Epitranscriptomic abnormalities, which are highly prevalent in primary central nervous system malignancies, have been identified as crucial contributors to the development and progression of gliomas. RNA epitranscriptomic modifications, particularly the reversible modification methylation, have been observed throughout the RNA cycle. Epitranscriptomic modifications, which regulate RNA transcription and translation, have profound biological implications. These modifications are associated with the development of several cancer types. Notably, three main protein types-writers, erasers, and readers, in conjunction with other related proteins, mediate these epitranscriptomic changes. This review primarily focuses on the role of recently identified RNA methylation modifications in gliomas, such as N6-methyladenosine (m6A), 5-methylcytosine (m5C), N7-methylguanosine (m7G), and N1-methyladenosine (m1A). We delved into their corresponding writers, erasers, readers, and related binding proteins to propose new approaches and prognostic indicators for patients with glioma.

Regulation and signaling pathways in cancer stem cells: implications for targeted therapy for cancer

Cancer stem cells (CSCs), initially identified in leukemia in 1994, constitute a distinct subset of tumor cells characterized by surface markers such as CD133, CD44, and ALDH. Their behavior is regulated through a complex interplay of networks, including transcriptional, post-transcriptional, epigenetic, tumor microenvironment (TME), and epithelial-mesenchymal transition (EMT) factors. Numerous signaling pathways were found to be involved in the regulatory network of CSCs. The maintenance of CSC characteristics plays a pivotal role in driving CSC-associated tumor metastasis and conferring resistance to therapy. Consequently, CSCs have emerged as promising targets in cancer treatment. To date, researchers have developed several anticancer agents tailored to specifically target CSCs, with some of these treatment strategies currently undergoing preclinical or clinical trials. In this review, we outline the origin and biological characteristics of CSCs, explore the regulatory networks governing CSCs, discuss the signaling pathways implicated in these networks, and investigate the influential factors contributing to therapy resistance in CSCs. Finally, we offer insights into preclinical and clinical agents designed to eliminate CSCs.

A novel m6A/m5C/m1A score signature to evaluate prognosis and its immunotherapy value in colon cancer patients

BACKGROUND

Colon cancer features strong heterogeneity and invasiveness, with high incidence and mortality rates. Recently, RNA modifications involving m6A, m5C, and m1A play a vital part in tumorigenesis and immune cell infiltration. However, integrated analysis among various RNA modifications in colon cancer has not been performed.

METHODS

RNA-seq profiling, clinical data and mutation data were obtained from The Cancer Genome Atlas and Gene Expression Omnibus. We first explored the mutation status and expression levels of m6A/m5C/m1A regulators in colon cancer. Then, different m6A/m5C/m1A clusters and gene clusters were identified by consensus clustering analysis. We further constructed and validated a scoring system, which could be utilized to accurately assess the risk of individuals and guide personalized immunotherapy. Finally, m6A/m5C/m1A regulators were validated by immunohistochemical staining and RT-qPCR.

RESULTS

In our study, three m6A/m5C/m1A clusters and gene clusters were identified. Most importantly, we constructed a m6A/m5C/m1A scoring system to assess the clinical risk of the individuals. Besides, the prognostic value of the score was validated with three independent cohorts. Moreover, the level of the immunophenoscore of the low m6A/m5C/m1A score group increased significantly with CTLA-4/PD-1 immunotherapy. Finally, we validated that the mRNA and protein expression of VIRMA and DNMT3B increased in colon cancer tissues.

CONCLUSIONS

We constructed and validated a stable and powerful m6A/m5C/m1A score signature to assess the survival outcomes and immune infiltration characteristics of colon cancer patients, which further guides optimization of personalized treatment, making it valuable for clinical translation and implementation.

Comprehensive landscape of the miRNA-regulated prognostic marker LAYN with immune infiltration and stemness in pan-cancer

BACKGROUND

COVID-19 has created a significant risk to worldwide public health. According to recent research, C-type lectins may be SARS-CoV-2 receptors. Layilin (LAYN), a broadly expressed integral membrane hyaluronan receptor with a C-type lectin structural domain, is a gene related to cell senescence. There are a few studies on C-type lectins in pan-cancer, and no pan-cancer analysis has been conducted for LAYN.

METHODS

The genotype tissue expression (GTEx) portal and the cancer genome map (TCGA) database were used to collect samples from healthy and cancer patients. Bioinformatics methods are used to construct immune landscape, mutation landscape, and stemness landscape of LAYN. The single-cell sequencing data were used from the CancerSEA website to analyze the functions of LAYN. The prognosis potential of LAYN was discussed based on machine learning.

RESULTS

LAYN is differentially expressed among cancers. Survival analysis indicated that LAYN was related to a poor overall survival (OS) rate in cancers, like HNSC, MESO, and OV. Mutational landscapes of LAYN in SKCM and STAD were constructed. LAYN was negatively related to Tumor Mutation Burden (TMB) in THCA, PRAD, and UCEC, and with the Microsatellite Instability (MSI) in STAD, LUAD, and UCEC. The immune landscape in pan-cancer suggested that LAYN may be involved in tumor immune escape. LAYN plays a crucial role in the infiltration of immune cells in malignant tumors. LAYN participates in methylation modifications and affects tumor proliferation and metastasis by regulating stemness. Analysis of single-cell sequencing data suggests that LAYN may participate in several biological processes, like stemness, apoptosis, and DNA repair. LAYN transcript was predicted as a liquid-liquid phase separation (LLPS)-related RNA. The results of KIRC were verified in the GEO and ArrayExpress databases. Furthermore, prognostic models based on machine learning of LAYN-related genes were established. Hsa-miR-153-5p and hsa-miR-505-3p may be the upstream miRNAs of LAYN and have a high value for tumor prognosis.

CONCLUSION

This study elucidated the functional mechanisms of LAYN from a pan-cancer perspective and provided novel insights into cancer prognosis, metastasis, and immunotherapy. LAYN has the potential to become a new target of mRNA vaccines and molecular therapies in tumors.

Identification of m5C-related lncRNAs signature to predict prognosis and therapeutic responses in esophageal squamous cell carcinoma patients

Esophageal squamous cell carcinoma (ESCC) has a dismal prognosis because of atypical early symptoms and heterogeneous therapeutic responses. 5-methylcytosine (m5C) modification plays an important role in the onset and development of many tumors and is widespread in long non-coding RNA (lncRNA) transcripts. However, the functions of m5C and lncRNAs in ESCC have not been completely elucidated. Herein, this study aimed to explore the role of m5C-related lncRNAs in ESCC. The RNA-seq transcriptome profiles and clinical information were downloaded from the TCGA-ESCC database. Pearson analysis was used to identify m5C-related lncRNAs. Then we established the m5C-related lncRNAs prognostic signature (m5C-LPS) using univariate Cox and least absolute shrinkage and selection operator (LASSO) regression analysis. Then, the prognostic value of m5C-LPS was evaluated internally and externally using the TCGA-ESCC and GSE53622 databases through multiple methods. We also detected the expression of these lncRNAs in ESCC cell lines and patient tissues. Fluorescence in situ hybridization (FISH) was used to detect the prognostic value of specific lncRNA. In addition, clinical parameters, immune status, genomic variants, oncogenic pathways, enrichment pathways, and therapeutic response features associated with m5C-LPS were explored using bioinformatics methods. We constructed and validated a prognostic signature based on 9 m5C-related lncRNAs (AC002091.2, AC009275.1, CAHM, LINC02057.1, AC0006329.1, AC037459.3, AC064807.1, ATP2B1-AS1, and UBAC2-AS1). The quantitative real-time polymerase chain reaction (qRT-PCR) revealed that most lncRNAs were upregulated in ESCC cell lines and patient tissues. And AC002091.2 was validated to have significant prognostic value in ESCC patients. A composite nomogram was generated to facilitate clinical practice by integrating this signature with the N stage. Besides, patients in the low-risk group were characterized by good clinical outcomes, favorable immune status, and low oncogenic alteration. Function enrichment analysis indicated that the risk score was associated with mRNA splicing, ncRNA processing, and DNA damage repair response. At the same time, we found significant differences in the responses to chemoradiotherapy between the two groups, proving the value of m5C-LPS in treatment decision-making in ESCC. This study established a novel prognostic signature based on 9 m5C-related lncRNAs, which is a promising biomarker for predicting clinical outcomes and therapeutic response in ESCC.

RNA methylations in depression, from pathological mechanism to therapeutic potential

Depression is caused by a variety of factors such as genetic factors, biological factors, and psychosocial factors, and the pathogenesis is complex. RNA methylations and related downstream signaling pathways influence a variety of biological mechanisms, including cell differentiation, tumorigenesis, sex determination, and stress response. In this work, we searched the PubMed, Web of Science, National Library of Science and Technology (NSTL), and ScienceDirect Online (SDOL) databases to summarize the biological roles of RNA methylations and their impact on the pathological mechanisms of depression. RNA methylations play a key role in the development of many diseases, and current research shows that RNA methylations are also closely linked to depression. RNA methylations in depression mainly involve "writers" (mediating the methylation modification process of RNAs), "erasers" (mediating the demethylation modification process of RNA methylation). Fat Mass and Obesity Associated (FTO) influences the development of depression by increasing body mass index (BMI), decreases the dopamine level, inhibits the adrenoceptor beta 2 (ADRB2)-c-Myc-sirt1 pathway, results in the m6A/m6Am dysregulation in brain, and may be involved in the pathogenesis of depression. The study of RNA methylations in depression has further deepened our understanding of the pathogenesis and development process of depression, provides new perspectives for the study of the pathological mechanism of depression, and provides new targets for the prevention and treatment of this disease.

Four types of adenine-related RNA modification writers -mediated molecular subtypes contribute to predicting clinical outcomes and treatment options in bladder cancer

Rationale

RNA modifications, containing m6A, m1A, alternative polyadenylation and adenosine-to-inosine RNA editing, involve in critical cancerous immunity and cancerous processes. However, the functional roles of RNA modification writers in bladder cancer (BLCA) are largely unknown.

Methods

In this study, unsupervised clustering was used to identify novel RNA modification writers -mediated molecular subtypes in BLCA. A corresponding quantitative indicator called WriterScore was developed using univariate Cox and Least absolute shrinkage and selection operator (LASSO) analysis. Then, we systematically analyzed the correlation between RNA modification writer-related clusters (WriterScore) and immunological characteristics, classical molecular subtypes, clinicopathologic features and treatment options in BLCA. Finally, we validated the WriterScore in multiple other external BLCA datasets, clinical sample dataset in Shengjing Hospital and pancancer.

Results

Two RNA modification writer-related clusters and three DEGclusters were obtained. These RNA modification writer-related clusters (WriterScore) were strongly associated with immunological characteristics, classical molecular subtypes, clinicopathologic features of BLCA. Moreover, WriterScore can properly predict the clinical outcomes and immunotherapy of BLCA patients.

Conclusion

Our study systematically investigated the role of RNA modification writers and developed a significant WriterScore to guide several treatment options in BLCA, which might bring some potential benefits for BLCA patients.

Connecting the Dots: N6-Methyladenosine (m6 A) Modification in Spermatogenesis

N6-methyladenosine (m6 A) is the most common RNA modification found in eukaryotes and is involved in multiple biological processes, including neuronal development, tumorigenesis, and gametogenesis. It is well known that methylation-modifying enzymes (classified into writers, erasers, and readers) mediate catalysis, clearance, and recognition of m6 A. Recent studies suggest that these genes may be associated with spermatogenesis. Numerous studies have revealed the m6 A role during spermatogenesis. However, the expression patterns and relationships of these m6 A enzymes during various stages of spermatogenesis remain unknown. In this review, it is aimed to provide an overview of m6 A enzyme functions and elucidate their potential mechanisms and regulatory relationships at a specific phase during spermatogenesis, providing new insights into the m6 A modification underlying the spermatogenesis process.

m5C regulator-mediated methylation modification phenotypes characterized by distinct tumor microenvironment immune heterogenicity in colorectal cancer

The RNA 5-methylcytosine (m5C) modification has been demonstrated to be an important epigenetic regulator and to impact colorectal cancer (CRC) progression. However, the potential roles of m5C modification in immune cell infiltration in the CRC tumor microenvironment (TME) remain unknown. The m5C modification phenotypes were comprehensively evaluated based on 14 m5C regulators in a meta-CRC cohort of 1792 patients and systematically correlated with the m5C modification phenotypes, immune cell infiltration characteristics and known biological processes. The m5Cscore model was constructed by principal component analysis (PCA) algorithms to quantify the m5C modification phenotypes of individual CRC samples and was used to predict the immunotherapy response. We identified three m5C modification phenotypes associated with distinct clinical outcomes and biological processes among the 1792 meta-CRC patients. Three phenotypes with a highly consistent TME landscape and characteristics were revealed: immune excluded, immune desert and immune inflammation. The meta-CRC patients were divided into high and low m5Cscore subgroups based on the m5Cscore. The m5Cscore was confirmed to have a negative correlation with infiltrating immune cells and PD-L1 expression and a positive correlation with tumor mutation burden (TMB), mutation rate and microsatellite instability (MSI) score. Moreover, patients in the low m5Cscore group had better immunotherapy responses and significant durable survival benefits in independent anti-PD-1/L1 immunotherapy cohorts for the immune checkpoint inhibitor (ICI) therapeutic strategy. This study revealed that m5C modification plays a crucial role in TME composition and complexity. Comprehensive evaluation of the m5C modification phenotypes of individual patients will enhance our understanding of TME characteristics and promote the application of more appropriate and personalized treatment strategies.

CRISPR-based m6A modification and its potential applications in telomerase regulation

Telomerase determines cell lifespan by controlling chromosome stability and cell viability, m⁶A epigenetic modification plays an important role in the regulation of telomerase activity. Using CRISPR epigenome editing to analyze specific m⁶A modification sites in telomerase will provide an important tool for analyzing the molecular mechanism of m⁶A modification regulating telomerase activity. In this review, we clarified the relevant applications of CRISPR system, paid special attention to the regulation of m⁶A modification in stem cells and cancer cells based on CRISPR system, emphasized the regulation of m⁶A modification on telomerase activity, pointed out that m⁶A modification sites regulate telomerase activity, and discussed strategies based on telomerase activity and disease treatment, which are helpful to promote the research of anti-aging and tumor related diseases.

Crosstalk between m6A and coding/non-coding RNA in cancer and detection methods of m6A modification residues

N6-methyladenosine (m6A) is one of the most common and well-known internal RNA modifications that occur on mRNAs or ncRNAs. It affects various aspects of RNA metabolism, including splicing, stability, translocation, and translation. An abundance of evidence demonstrates that m6A plays a crucial role in various pathological and biological processes, especially in tumorigenesis and tumor progression. In this article, we introduce the potential functions of m6A regulators, including "writers" that install m6A marks, "erasers" that demethylate m6A, and "readers" that determine the fate of m6A-modified targets. We have conducted a review on the molecular functions of m6A, focusing on both coding and noncoding RNAs. Additionally, we have compiled an overview of the effects noncoding RNAs have on m6A regulators and explored the dual roles of m6A in the development and advancement of cancer. Our review also includes a detailed summary of the most advanced databases for m6A, state-of-the-art experimental and sequencing detection methods, and machine learning-based computational predictors for identifying m6A sites.

The prognostic value and immune landscaps of m6A/m5C-related lncRNAs signature in the low grade glioma

BACKGROUND

N6-methyladenosine (m6A) and 5-methylcytosine (m5C) are the main RNA methylation modifications involved in the oncogenesis of cancer. However, it remains obscure whether m6A/m5C-related long non-coding RNAs (lncRNAs) affect the development and progression of low grade gliomas (LGG).

METHODS

We summarized 926 LGG tumor samples with RNA-seq data and clinical information from The Cancer Genome Atlas and Chinese Glioma Genome Atlas. 105 normal brain samples with RNA-seq data from the Genotype Tissue Expression project were collected for control. We obtained a molecular classification cluster from the expression pattern of sreened lncRNAs. The least absolute shrinkage and selection operator Cox regression was employed to construct a m6A/m5C-related lncRNAs prognostic signature of LGG. In vitro experiments were employed to validate the biological functions of lncRNAs in our risk model.

RESULTS

The expression pattern of 14 sreened highly correlated lncRNAs could cluster samples into two groups, in which various clinicopathological features and the tumor immune microenvironment were significantly distinct. The survival time of cluster 1 was significantly reduced compared with cluster 2. This prognostic signature is based on 8 m6A/m5C-related lncRNAs (GDNF-AS1, HOXA-AS3, LINC00346, LINC00664, LINC00665, MIR155HG, NEAT1, RHPN1-AS1). Patients in the high-risk group harbored shorter survival times. Immunity microenvironment analysis showed B cells, CD4 + T cells, macrophages, and myeloid-derived DC cells were significantly increased in the high-risk group. Patients in high-risk group had the worse overall survival time regardless of followed TMZ therapy or radiotherapy. All observed results from the TCGA-LGG cohort could be validated in CGGA cohort. Afterwards, LINC00664 was found to promote cell viability, invasion and migration ability of glioma cells in vitro.

CONCLUSION

Our study elucidated a prognostic prediction model of LGG by 8 m6A/m5C methylated lncRNAs and a critical lncRNA regulation function involved in LGG progression. High-risk patients have shorter survival times and a pro-tumor immune microenvironment.

5-Methylcytosine-related lncRNAs: predicting prognosis and identifying hot and cold tumor subtypes in head and neck squamous cell carcinoma

BACKGROUND

5-Methylcytosine (m⁵C) methylation is recognized as an mRNA modification that participates in biological progression by regulating related lncRNAs. In this research, we explored the relationship between m⁵C-related lncRNAs (mrlncRNAs) and head and neck squamous cell carcinoma (HNSCC) to establish a predictive model.

METHODS

RNA sequencing and related information were obtained from the TCGA database, and patients were divided into two sets to establish and verify the risk model while identifying prognostic mrlncRNAs. Areas under the ROC curves were assessed to evaluate the predictive effectiveness, and a predictive nomogram was constructed for further prediction. Subsequently, the tumor mutation burden (TMB), stemness, functional enrichment analysis, tumor microenvironment, and immunotherapeutic and chemotherapeutic responses were also assessed based on this novel risk model. Moreover, patients were regrouped into subtypes according to the expression of model mrlncRNAs.

RESULTS

Assessed by the predictive risk model, patients were distinguished into the low-MLRS and high-MLRS groups, showing satisfactory predictive effects with AUCs of 0.673, 0.712, and 0.681 for the ROCs, respectively. Patients in the low-MLRS groups exhibited better survival status, lower mutated frequency, and lower stemness but were more sensitive to immunotherapeutic response, whereas the high-MLRS group appeared to have higher sensitivity to chemotherapy. Subsequently, patients were regrouped into two clusters: cluster 1 displayed immunosuppressive status, but cluster 2 behaved as a hot tumor with a better immunotherapeutic response.

CONCLUSIONS

Referring to the above results, we established a m⁵C-related lncRNA model to evaluate the prognosis, TME, TMB, and clinical treatments for HNSCC patients. This novel assessment system is able to precisely predict the patients' prognosis and identify hot and cold tumor subtypes clearly for HNSCC patients, providing ideas for clinical treatment.

Epitranscriptional Regulation: From the Perspectives of Cardiovascular Bioengineering

The central dogma of gene expression involves DNA transcription to RNA and RNA translation into protein. As key intermediaries and modifiers, RNAs undergo various forms of modifications such as methylation, pseudouridylation, deamination, and hydroxylation. These modifications, termed epitranscriptional regulations, lead to functional changes in RNAs. Recent studies have demonstrated crucial roles for RNA modifications in gene translation, DNA damage response, and cell fate regulation. Epitranscriptional modifications play an essential role in development, mechanosensing, atherogenesis, and regeneration in the cardiovascular (CV) system, and their elucidation is critically important to understanding the molecular mechanisms underlying CV physiology and pathophysiology. This review aims at providing biomedical engineers with an overview of the epitranscriptome landscape, related key concepts, recent findings in epitranscriptional regulations, and tools for epitranscriptome analysis. The potential applications of this important field in biomedical engineering research are discussed.

ac4C acetylation regulates mRNA stability and translation efficiency in osteosarcoma

Objective

N4-acetylcytidine (ac4C) acetylation can promote target gene expression through improved mRNA stability. To explore the role of ac4C acetylation in osteosarcoma, U2OS and MG63 cell lines were treated with the N-acetyltransferase 10 (NAT10) inhibitor Remodelin. Reverse transcription-polymerase chain reaction (RT-PCR) and Western blot were used to test the gene and protein expression efficiency.

Methods

The proliferation rate of osteosarcoma cells was measured by a cell counting kit-8 (CCK8) assay. The cell cycle and apoptosis were analyzed by flow cytometry. The invasiveness of osteosarcoma cells was detected by a transwell invasion assay. The ac4C acetylation of target genes was screened by acetylated RNA immunoprecipitation and sequencing (acRIP-seq).

Results

We found that when osteosarcoma cells were treated with Remodelin at the optimal concentration, their NAT10 expression and the cell proliferation was inhibited, the cells in the G1 phase increased (P < 0.05) but those in the S phase decreased, the apoptotic cells in the early and late stages increased, and the cells invasiveness decreased (P < 0.05).

Conclusions

The farnesyltransferase subunit beta gene (FNTB) was identified by acRIP-seq as one of the target genes of ac4C acetylation and was further verified by RT-PCR and Western blot analyses. Remodelin was demonstrated to reduce the stability and protein translation efficiency of target gene mRNA in osteosarcoma cells. In conclusion, inhibition of ac4C acetylation in osteosarcoma can block proliferation and metastasis as well as promote apoptosis and cell cycle arrest. Ac4C acetylation contributes to the stability and protein translation efficiency of the downstream target gene mRNA.

Disruption of the mouse liver epitranscriptome by long-term aroclor 1260 exposure

Chronic environmental exposure to polychlorinated biphenyls (PCBs) is associated with non-alcoholic fatty liver disease (NAFLD) and exacerbated by a high fat diet (HFD). Here, chronic (34 wks.) exposure of low fat diet (LFD)-fed male mice to Aroclor 1260 (Ar1260), a non-dioxin-like (NDL) mixture of PCBs, resulted in steatohepatitis and NAFLD. Twelve hepatic RNA modifications were altered with Ar1260 exposure including reduced abundance of 2'-O-methyladenosine (Am) and N(6)-methyladenosine (m6A), in contrast to increased Am in the livers of HFD-fed, Ar1260-exposed mice reported previously. Differences in 13 RNA modifications between LFD- and HFD- fed mice, suggest that diet regulates the liver epitranscriptome. Integrated network analysis of epitranscriptomic modifications identified a NRF2 (Nfe2l2) pathway in the chronic, LFD, Ar1260-exposed livers and an NFATC4 (Nfatc4) pathway for LFD- vs. HFD-fed mice. Changes in protein abundance were validated. The results demonstrate that diet and Ar1260 exposure alter the liver epitranscriptome in pathways associated with NAFLD.

Vital roles of m5C RNA modification in cancer and immune cell biology

RNA modification plays an important role in epigenetics at the posttranscriptional level, and 5-methylcytosine (m5C) has attracted increasing attention in recent years due to the improvement in RNA m5C site detection methods. By influencing transcription, transportation and translation, m5C modification of mRNA, tRNA, rRNA, lncRNA and other RNAs has been proven to affect gene expression and metabolism and is associated with a wide range of diseases, including malignant cancers. RNA m5C modifications also substantially impact the tumor microenvironment (TME) by targeting different groups of immune cells, including B cells, T cells, macrophages, granulocytes, NK cells, dendritic cells and mast cells. Alterations in immune cell expression, infiltration and activation are highly linked to tumor malignancy and patient prognosis. This review provides a novel and holistic examination of m5C-mediated cancer development by examining the exact mechanisms underlying the oncogenicity of m5C RNA modification and summarizing the biological effects of m5C RNA modification on tumor cells as well as immune cells. Understanding methylation-related tumorigenesis can provide useful insights for the diagnosis as well as the treatment of cancer.

Role of m6A modification in immune microenvironment of digestive system tumors

Digestive system tumors are huge health problem worldwide, largely attributable to poor dietary choices. The role of RNA modifications in cancer development is an emerging field of research. RNA modifications are associated with the growth and development of various immune cells, which, in turn, regulate the immune response. The majority of RNA modifications are methylation modifications, and the most common type is the N6-methyladenosine (m6A) modification. Here, we reviewed the molecular mechanism of m6A in the immune cells and the role of m6A in the digestive system tumors. However, further studies are required to better understand the role of RNA methylation in human cancers for designing diagnostic and treatment strategies and predicting the prognosis of patients.

Identification of m7G-Related miRNA Signatures Associated with Prognosis, Oxidative Stress, and Immune Landscape in Lung Adenocarcinoma

The role of N7-methylguanosine(m7G)-related miRNAs in lung adenocarcinoma (LUAD) remains unclear. We used LUAD data from The Cancer Genome Atlas (TCGA) to establish a risk model based on the m7G-related miRNAs, and divided patients into high-risk or low-risk subgroups. A nomogram for predicting overall survival (OS) was then constructed based on the independent risk factors. In addition, we performed a functional enrichment analysis and defined the oxidative stress-related genes, immune landscape as well as a drug response profile in the high-risk and low-risk subgroups. This study incorporated 28 m7G-related miRNAs into the risk model. The data showed a significant difference in the OS between the high-risk and low-risk subgroups. The receiver operating characteristic curve (ROC) predicted that the area under the curve (AUC) of one-year, three-year and five-year OS was 0.781, 0.804 and 0.853, respectively. The C-index of the prognostic nomogram for predicting OS was 0.739. We then analyzed the oxidative stress-related genes and immune landscape in the high-risk and low-risk subgroups. The data demonstrated significant differences in the expression of albumin (ALB), estimated score, immune score, stromal score, immune cell infiltration and functions between the high-risk and low-risk subgroups. In addition, the drug response analysis showed that low-risk subgroups may be more sensitive to tyrosine kinase inhibitor (TKI) and histone deacetylase (HDAC) inhibitors. We successfully developed a novel risk model based on m7G-related miRNAs in this study. The model can predict clinical prognosis and guide therapeutic regimens in patients with LUAD. Our data also provided new insights into the molecular mechanisms of m7G in LUAD.

The m6A/m5C/m1A regulator genes signature reveals the prognosis and is related with immune microenvironment for hepatocellular carcinoma

BACKGROUND

RNA methylation is a crucial in many biological functions, and its aberrant regulation is associated with cancer progression. N6-Methyladenosine (m6A), 5-Methylcytosine (m5C), N1-methyladenosine (m1A) are common modifications of RNA methylation. However, the effect of methylation of m6A/m5C/m1A in hepatocellular carcinoma (HCC) remains unclear.

METHOD

The transcriptome datasets, clinic information, and mutational data of 48 m6A/m5C/m1A regulator genes were acquired from the TCGA database, and the prognostic hazard model was established by univariate and Least absolute shrinkage and selection operator (Lasso) regression. The multivariate regression was performed to determine whether the risk score was an independent prognostic indicator. Kaplan-Meier survival analysis and ROC curve analysis were used to evaluate the predictive ability of the risk model. Decision curve analysis（DCA）analysis was conducted to estimate the clinical utility of the risk model. We further analyzed the association between risk score and functional enrichment, tumor immune microenvironment, and somatic mutation.

RESULT

The four-gene (YTHDF1, YBX1, TRMT10C, TRMT61A) risk signature was constructed. The high-risk group had shorter overall survival (OS) than the low-risk group. Univariate and multivariate regression analysis indicated that risk score was an independent prognostic indicator. Risk scores in male group, T3 + T4 group and Stage III + IV group were higher in female group, T1 + T2 group and stage I + II group. The AUC values for 1-, 2-, and 3-year OS in the TCGA dataset were 0.764, 0.693, and 0.689, respectively. DCA analysis showed that the risk score had a higher clinical net benefit in 1- and 2-year OS than other clinical features.The risk score was positively related to some immune cell infiltration and most immune checkpoints.

CONCLUSION

We developed a novel m6A/m5C/m1A regulator genes' prognostic model, which could be applied as a latent prognostic tool for HCC and might guide the choice of immunotherapies.

An innovative model based on N7-methylguanosine-related lncRNAs for forecasting prognosis and tumor immune landscape in bladder cancer

BACKGROUND

As a novel type of the prevalent post-transcriptional modifications, N7-methylguanosine (m7G) modification is essential in the tumorigenesis, progression, and invasion of many cancers, including bladder cancer (BCa). However, the integrated roles of m7G-related lncRNAs in BCa remain undiscovered. This study aims to develop a prognostic model based on the m7G-related lncRNAs and explore its predictive value of the prognosis and anti-cancer treatment sensitivity.

METHODS

We obtained RNA-seq data and corresponding clinicopathological information from the TCGA database and collected m7G-related genes from previous studies and GSEA. Based on LASSO and Cox regression analysis, we developed a m7G prognostic model. The Kaplan-Meier (K-M) survival analysis and ROC curves were performed to evaluate the predictive power of the model. Gene set enrichment analysis (GSEA) was conducted to explore the molecular mechanisms behind apparent discrepancies between the low- and high-risk groups. We also investigated immune cell infiltration, TIDE score, TMB, the sensitivity of common chemotherapy drugs, and the response to immunotherapy between the two risk groups. Finally, we validated the expression levels of these ten m7G-related lncRNAs in BCa cell lines by qRT-PCR.

RESULTS

We developed a m7G prognostic model (risk score) composed of 10 m7G-related lncRNAs that are significantly associated with the OS of BCa patients. The K-M survival curves revealed that the high-risk group patients had significantly worse OS than those in the low-risk group. The Cox regression analysis confirmed that the risk score was a significant independent prognostic factor for BCa patients. We found that the high-risk group had higher the immune scores and immune cell infiltration. Furthermore, the results of the sensitivity of common anti-BCa drugs showed that the high-risk group was more sensitive to neoadjuvant cisplatin-based chemotherapy and anti-PD1 immunotherapy. Finally, qRT-PCR revealed that AC006058.1, AC073133.2, LINC00677, and LINC01338 were significantly downregulated in BCa cell lines, while the expression levels of AC124312.2 and AL158209.1 were significantly upregulated in BCa cell lines compared with normal cell lines.

CONCLUSION

The m7G prognostic model can be applied to accurately predict the prognosis and provide robust directions for clinicians to develop better individual-based and precise treatment strategies for BCa patients.

The importance of pseudouridylation: human disorders related to the fifth nucleoside

Pseudouridylation is one of the most abundant RNA modifications in eukaryotes, making pseudouridine known as the "fifth nucleoside." This highly conserved alteration affects all non-coding and coding RNA types. Its role and importance have been increasingly widely researched, especially considering that its absence or damage leads to serious hereditary diseases. Here, we summarize the human genetic disorders described to date that are related to the participants of the pseudouridylation process.

The lncRNA epigenetics: The significance of m6A and m5C lncRNA modifications in cancer

Most of our transcribed RNAs are represented by non-coding sequences. Long non-coding RNAs (lncRNAs) are transcripts with no or very limited protein coding ability and a length &gt;200nt. They can be epigenetically modified. N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), 7-methylguanosine (m7G) and 2’-O-methylation (Nm) are some of the lncRNAs epigenetic modifications. The epigenetic modifications of RNA are controlled by three classes of enzymes, each playing a role in a specific phase of the modification. These enzymes are defined as “writers”, “readers” and “erasers”. m6A and m5C are the most studied epigenetic modifications in RNA. These modifications alter the structure and properties, thus modulating the functions and interactions of lncRNAs. The aberrant expression of several lncRNAs is linked to the development of a variety of cancers and the epigenetic signatures of m6A- or m5C-related lncRNAs are increasingly recognized as potential biomarkers of prognosis, predictors of disease stage and overall survival. In the present manuscript, the most up to date literature is reviewed with the focus on m6A and m5C modifications of lncRNAs and their significance in cancer.

Kinesin family member C 1 overexpression exerts tumor-promoting properties in head and neck squamous cell carcinoma via the Rac1/Wnt/β-catenin pathway

Kinesin family member C 1 (KIFC1) is a kinesin-14 motor protein, and its abnormal upregulation promotes the malignant behavior of cancer cells. N6-methyladenosine (m6A) RNA methylation is a common modification of eukaryotic mRNA and affects RNA expression. Herein, we explored how KIFC1 regulated head and neck squamous cell carcinoma (HNSCC) tumorigenesis and how m6A modification affected KIFC1 expression. Bioinformatics analysis was performed to screen for genes of interest, and in vitro and in vivo studies were carried out to investigate the function and mechanism of KIFC1 in HNSCC. We observed that the expression of KIFC1 in HNSCC tissues was significantly higher than in normal or adjacent normal tissues. Cancer patients with higher KIFC1 expression have lower tumor differentiation status. Demethylase alkB homolog 5 (ALKBH5), a cancer-promoting factor in HNSCC, could interact with KIFC1 mRNA and post-transcriptionally activated KIFC1 through m6A modification. KIFC1 downregulation suppressed HNSCC cell growth and metastasis in vivo and in vitro. However, overexpression of KIFC1 promoted these malignant behaviors. We demonstrated that KIFC1 overexpression activated the oncogenic Wnt/β-catenin pathway. KIFC1 interacted with the small GTPase Ras-related C3 botulinum toxin substrate 1 (Rac1) at the protein level and increased activity. The Rho GTPase Rac1 was indicated to be an upstream activator of the Wnt/β-catenin signaling pathway, and its Rac1 inhibitor, NSC-23766, treatment reversed the effects caused by KIFC1 overexpression. Those observations demonstrate that abnormal expression of KIFC1 may be regulated by demethylase ALKBH5 in an m6A-dependent manner and promote HNSCC progression via the Rac1/Wnt/β-catenin pathway.

m5C-dependent cross-regulation between nuclear reader ALYREF and writer NSUN2 promotes urothelial bladder cancer malignancy through facilitating RABL6/TK1 mRNAs splicing and stabilization

The significance of 5-methylcytosine (m⁵C) methylation in human malignancies has become an increasing focus of investigation. Here, we show that m⁵C regulators including writers, readers and erasers, are predominantly upregulated in urothelial carcinoma of the bladder (UCB) derived from Sun Yat-sen University Cancer Center and The Cancer Genome Atlas cohort. In addition, NOP2/Sun RNA methyltransferase family member 2 (NSUN2) as a methyltransferase and Aly/REF export factor (ALYREF) as a nuclear m⁵C reader, are frequently coexpressed in UCB. By applying patient-derived organoids model and orthotopic xenograft mice model, we demonstrate that ALYREF enhances proliferation and invasion of UCB cells in an m⁵C-dependent manner. Integration of tanscriptome-wide RNA bisulphite sequencing (BisSeq), RNA-sequencing (RNA-seq) and RNA Immunoprecipitation (RIP)-seq analysis revealed that ALYREF specifically binds to hypermethylated m⁵C site in RAB, member RAS oncogene family like 6 (RABL6) and thymidine kinase 1 (TK1) mRNA via its K171 domain. ALYREF controls UCB malignancies through promoting hypermethylated RABL6 and TK1 mRNA for splicing and stabilization. Moreover, ALYREF recognizes hypermethylated m⁵C site of NSUN2, resulting in NSUN2 upregulation in UCB. Clinically, the patients with high coexpression of ALYREF/RABL6/TK1 axis had the poorest overall survival. Our study unveils an m⁵C dependent cross-regulation between nuclear reader ALYREF and m⁵C writer NSUN2 in activation of hypermethylated m⁵C oncogenic RNA through promoting splicing and maintaining stabilization, consequently leading to tumor progression, which provides profound insights into therapeutic strategy for UCB.

Comprehensive Pan-Cancer Analysis of KIF18A as a Marker for Prognosis and Immunity

KIF18A belongs to the Kinesin family, which participates in the occurrence and progression of tumors. However, few pan-cancer analyses have been performed on KIF18A to date. We used multiple public databases such as TIMER, The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Human Protein Atlas (HPA) to explore KIF18A mRNA expression in 33 tumors. We performed immunohistochemistry on liver cancer and pancreatic cancer tissues and corresponding normal tissues to examine the expression of KIF18A protein. Univariate Cox regression and Kaplan-Meier survival analysis were applied to detect the effect of KIF18A on overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI) of patients with these tumors. Subsequently, we explored KIF18A gene alterations in different tumor tissues using cBioPortal. The relationship between KIF18A and clinical characteristics, tumor microenvironment (TME), immune regulatory genes, immune checkpoints, tumor mutational burden (TMB), microsatellite instability (MSI), mismatch repairs (MMRs), DNA methylation, RNA methylation, and drug sensitivity was applied for further study using the R language. Gene Set Enrichment Analysis (GSEA) was utilized to explore the molecular mechanism of KIF18A. Bioinformatic analysis and immunohistochemical experiments confirmed that KIF18A was up-regulated in 27 tumors and was correlated with the T stage, N stage, pathological stage, histological grade, and Ki-67 index in many cancers. The overexpression of KIF18A had poor OS, DSS, and PFI in adrenocortical carcinoma (ACC), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), brain lower-grade glioma (LGG), liver cancer (LIHC), lung adenocarcinoma (LUAD), and pancreatic cancer (PAAD). Univariate and multivariate regression analysis confirmed KIF18A as an independent prognostic factor for LIHC and PAAD. The mutation frequency of KIF18A is the highest in endometrial cancer. KIF18A expression levels were positively associated with immunocyte infiltration, immune regulatory genes, immune checkpoints, TMB, MSI, MMRs, DNA methylation, RNA methylation, and drug sensitivity in certain cancers. In addition, we discovered that KIF18A participated in the cell cycle at the single-cell level and GSEA analysis for most cancers. These findings suggested that KIF18A could be regarded as a latent prognostic marker and a new target for cancer immunological therapy.

Inactivation of KDM5A suppresses growth and enhances chemosensitivity in liver cancer by modulating ROCK1/PTEN/AKT pathway

Liver cancer is a kind of malignant tumor with poor sensitivity to chemotherapy. It is urgent to investigate approaches to improve the outcome of chemotherapy. KDM5A has been reported to be an oncogene in various cancers and is associated with drug resistance. However, the functions of KDM5A in chemotherapeutic sensitivity of liver cancer not been well illustrated. In this study, we found that KDM5A was upregulated in liver cancer tissue and cell lines. KDM5A knockdown using a gene interference strategy suppressed the growth of liver cancer in vitro and in vivo. CPI-455, a pharmacological inactivation of KDM5A enhanced the cytotoxicity of cisplatin (CDDP) in liver cells. CPI-455 and CDDP cotreatment resulted in apoptosis and mitochondrial dysfunction. We also found that knockdown or inactivation of KDM5A resulted in the downregulation of ROCK1, an oncogene regulating the activation of the PTEN/AKT signaling pathway. In particular, overexpression of ROCK1 or SF1670, a pharmacological inhibitor of PTEN, alleviated the cytotoxicity of CPI-455 and CDDP cotreatment. In HCCLM3 xenografts, CPI-455 and CDDP cotreatment dramatically inhibited the growth of xenograft tumor compared to CPI-455 or CDDP treatment alone. In conclusion, this study suggested that targeting the inactivation of KDM5A is an efficient strategy to enhance the chemosensitivity of liver cancer cells to CDDP by modulating the ROCK1/PTEN/AKT signaling pathway.