BACKGROUND: Radioresistance is the leading cause of death in advanced cervical cancer (CC). Dysregulation of RNA modification has recently emerged as a regulatory mechanism in radiation and drug resistance. We aimed to explore the biological function and clinical significance of 5-methylcytosine (m5C) in cervical cancer radiosensitivity.
METHODS: The abundance of RNA modification in radiotherapy-resistant and sensitive CC specimens was quantified by liquid chromatography-tandem mass spectrometry. The essential RNA modification-related genes involved in CC radiosensitivity were screened via RNA sequencing. The effect of NSUN6 on radiosensitivity was verified in CC cell lines, cell-derived xenograft (CDX), and 3D bioprinted patient-derived organoid (PDO). The mechanisms of NSUN6 in regulating CC radiosensitivity were investigated by integrative m5C sequencing, mRNA sequencing, and RNA immunoprecipitation.
RESULTS: We found a higher abundance of m5C modification in resistant CC samples, and NSUN6 was the essential m5C-regulating gene concerning radiosensitivity. NSUN6 overexpression was clinically correlated with radioresistance and poor prognosis in cervical cancer. Functionally, higher NSUN6 expression was associated with radioresistance in the 3D PDO model of cervical cancer. Moreover, silencing NSUN6 increased CC radiosensitivity in vivo and in vitro. Mechanistically, NDRG1 was one of the downstream target genes of NSUN6 identified by integrated m5C-seq, mRNA-seq, and functional validation. NSUN6 promoted the m5C modification of NDRG1 mRNA, and the m5C reader ALYREF bound explicitly to the m5C-labeled NDRG1 mRNA and enhanced NDRG1 mRNA stability. NDRG1 overexpression promoted homologous recombination-mediated DNA repair, which in turn led to radioresistance in cervical cancer.
CONCLUSIONS: Aberrant m5C hypermethylation and NSUN6 overexpression drive resistance to radiotherapy in cervical cancer. Elevated NSUN6 expression promotes radioresistance in cervical cancer by activating the NSUN6/ALYREF-m5C-NDRG1 pathway. The low expression of NSUN6 in cervical cancer indicates sensitivity to radiotherapy and a better prognosis.

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.
CONCLUSIONS: 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.

N5-methylcytosine (m5C) methylation modification plays a crucial role in the epigenetic mechanisms underlying tumorigenesis, aggressiveness, and malignancy in diffuse glioma. Our study aimed to develop a novel prognostic risk-scoring system to assess the impact of m5C modification in glioma patients. Initially, we identified two distinct m5C clusters based on the expression level of m5C regulators in The Cancer Genome Atlas glioblastoma (TCGA-GBM) dataset. Differentially expressed genes (DEGs) between the two m5C cluster groups were determined. Utilizing these m5C regulation-related DEGs, we classified glioma patients into three gene cluster groups: A, B, and C. Subsequently, an m5C scoring system was developed through a univariate Cox regression model, quantifying the m5C modification patterns utilizing six DEGs associated with disease prognosis. The resulting scoring system allowed us to categorize patients into high- or low-risk groups based on their m5C scores. In test (TCGA-GBM) and validation (Chinese Glioma Genome Atlas [CGGA]-1018 and CGGA-301) datasets, glioma patients with a higher m5C score consistently exhibited shorter survival durations, fewer isocitrate dehydrogenase (IDH) mutations, less 1p/19q codeletion and higher World Health Organization (WHO) grades. Additionally, distinct immune cell infiltration characteristics were observed among different m5C cluster groups and risk groups. Our study developed a novel prognostic scoring system based on m5C modification patterns for glioma patients, complementing existing molecular classifications and providing valuable insights into prognosis for glioma patients.

UNASSIGNED: Neuroblastoma (NB) is a prevalent pediatric tumor originating from primordial neural crest cells. As one of the latest epigenetics investigations focuses, RNA 5-methylcytosine (m5C) is closely related to cancer risk. TET methylcytosine dioxygenase 3 (TET3) is a demethylase for m5C modification. Whether there is an association between TET3 gene polymorphisms and neuroblastoma risk remains unclear.
UNASSIGNED: We conducted an epidemiological study in 402 patients and 473 controls to evaluate the relationship between TET3 gene SNPs (rs7560668 T > C, rs828867 G > A, and rs6546891 A > G) and NB susceptibility.
UNASSIGNED: Our results showed that rs828867 G > A significantly reduced NB risk in Chinese children [GA vs. GG, adjusted odds ratio (OR) = 0.72, 95% confidence interval (CI) = 0.52-0.98, P=0.040; GA/AA vs. GG, adjusted OR = 0.74, 95% CI = 0.55-0.998, P=0.048]. Individuals with 2-3 risk genotypes had a significantly higher NB risk than those with 0-1 risk genotypes (adjusted OR = 1.40, 95% CI = 1.04-1.88, P=0.027). The stratified analysis showed that the rs828867 G > A associated with decreased NB risk is remarkable among children aged >18 months (adjusted OR = 0.67, 95% CI = 0.46-0.96, P=0.029) and patients at clinical III + IV stages (adjusted OR = 0.67, 95% CI = 0.45-0.98, P=0.040). Compared with the 0-1 risk genotype, the concurrence of 2-3 risk genotypes significantly increased NB risk in the following subgroups: children aged >18 months and patients at clinical III + IV stages. GTEx analysis suggested that rs828867 G > A was significantly associated with RP11-287D1.4 and POLE4 mRNA expression.
UNASSIGNED: Overall, our results revealed that rs828867 G > A in the TET3 gene is significantly associated with predisposition to NB.

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors with a high mortality rate. The 5-methylcytosine (m5C), a type of RNA modification, plays crucial regulatory roles in HCC carcinogenesis, metastasis, and prognosis. However, a few studies have investigated the effect of genetic variants in m5C modification genes on survival of patients with hepatitis B virus (HBV)-related HCC. In the present study, we evaluated associations between 144 SNPs in 15 m5C modification genes and overall survival (OS) in 866 patients with the HBV-related HCC. Expression quantitative trait loci (eQTL) analysis and differential expression analysis were conducted to investigate biological mechanisms. As a result, we identified that two SNPs (NSUN7 rs2437325 A > G and TRDMT1 rs34434809 G > C) were significantly associated with HBV-related HCC OS with adjusted allelic hazards ratios of 1.25 (95% confidence interval = 1.05-1.48 and P = 0.011) and 1.19 (1.02-1.38 and P = 0.027), respectively, with a trend of combined risk genotypes (Ptrend < 0.001). Moreover, the results of eQTL analyses showed that both NSUN7 rs2437325 G and TRDMT1 rs34434809 C alleles were associated with a reduced mRNA expression level in 208 normal liver tissues (P = 0.007 and P < 0.001, respectively). Taken together, genetic variants in the m5C modification genes may be potential prognostic biomarkers of HBV-related HCC after hepatectomy, likely through mediating the mRNA expression of corresponding genes.

RNA 5-methylcytosine (m5C) is an abundant chemical modification in mammalian RNAs and plays crucial roles in regulating vital physiological and pathological processes, especially in cancer. However, the dysregulation of m5C and its underlying mechanisms in non-small cell lung cancer (NSCLC) remain unclear. Here we identified that NSUN2, a key RNA m5C methyltransferase, is highly expressed in NSCLC tumor tissue. We found elevated NSUN2 expression levels strongly correlate with tumor grade and size, predicting poor outcomes for NSCLC patients. Furthermore, RNA-seq and subsequent confirmation studies revealed the antioxidant-promoting transcription factor NRF2 is a target of NSUN2, and depleting NSUN2 decreases the expression of NRF2 and increases the sensitivity of NSCLC cells to ferroptosis activators both in vitro and in vivo. Intriguingly, the methylated-RIP-qPCR assay results indicated that NRF2 mRNA has a higher m5C level when NSUN2 is overexpressed in NSCLC cells but shows no significant changes in the NSUN2 methyltransferase-deficient group. Mechanistically, we confirmed that NSUN2 upregulates the expression of NRF2 by enhancing the stability of NRF2 mRNA through the m5C modification within its 5\'UTR region recognized by the specific m5C reader protein YBX1, rather than influencing its translation. In subsequent rescue experiments, we show knocking down NRF2 diminished the proliferation, migration, and ferroptosis tolerance mediated by NSUN2 overexpression. In conclusion, our study unveils a novel regulatory mechanism in which NSUN2 sustains NRF2 expression through an m5C-YBX1-axis, suggesting that targeting NSUN2 and its regulated ferroptosis pathway might offer promising therapeutic strategies for NSCLC patients.

Circular RNAs (circRNAs) are a class of stable non-coding RNAs that have emerged as key regulators in human diseases including cancer. This study investigates the role of circRNA_0102913 (circ_0102913) in malignant behavior of colorectal cancer (CRC) cells and the underpinning mechanisms. By analyzing CRC-related GSE197991, GSE159669, and GSE223001 datasets, we obtained circ_0102913 as an aberrantly upregulated circRNA in CRC. Increased circ_0102913 expression was detected in CRC tissues and cells. By querying multiple bioinformatics systems (circBank, Circular RNA Interactome, TargetScan, miRDIP, miRwalk, and miRDB), we identified microRNA-571 (miR-571) as a target of circ_0102913 and Rac family small GTPase 2 (RAC2) mRNA as a target of miR-571. Biotinylated-RNA pull-down and/or luciferase assays showed that circ_0102913 bound to miR-571 to restore the expression of RAC2 mRNA. Circ_0102913 silencing or miR-571 overexpression repressed proliferation, migration and invasion, and in vivo tumorigenesis abilities of CRC cells. However, the malignant properties of cells were restored by RAC2 overexpression. The increased circ_0102913 expression in CRC cells was attributed to increased 5-methylcytosine (m5C) modification levels. Silencing of NOP2/Sun RNA methyltransferase 5 reduced the m5C level and therefore reduced stability and expression of circ_0102913 expression in CRC cells. In conclusion, this study demonstrates that m5C-mediated upregulation of circ_0102913 augments malignant properties of CRC cells through a miR-571/RAC2 axis.

5-methylcytosine (m5C modification) plays an essential role in tumors, which affects different types of RNA, the expression of downstream target genes, and downstream pathways, thus participating in the tumor process. However, the effect of m5C modification on RNA in tumors and the exact mechanism have not been systematically reviewed. Therefore, we reviewed the status and sites of m5C modification, as well as the expression pattern and biological functions of m5C regulators in tumors, and further summarized the effects and regulation mechanism of m5C modification on messenger RNA (mRNA), ribosomal RNA (rRNA), transfer RNA (tRNA), long non-coding RNA (lncRNA) and other RNA in tumors. Finally, we summed up the interaction network, potential application, and value in clinical diagnosis and treatment of tumors. Taken together, this review benefits revealing the mechanism of m5C modification in tumor progression and provides new strategies for tumor diagnosis and treatment.

The 5-methylcytosine (m5C) is the key chemical modification in RNAs. As one of the demethylases in m5C, TET2 has been shown as a tumor suppressor. However, the impact of TET2 gene polymorphisms on neuroblastoma has not been elucidated. 402 neuroblastoma patients and 473 controls were genotyped for TET2 gene polymorphisms using the TaqMan method. The impact of TET2 gene polymorphisms on neuroblastoma susceptibility was determined using multivariate logistic regression analysis. We also adopted genotype-tissue expression database to explore the impact of TET2 gene polymorphisms on the expression of host and nearby genes. We used the R2 platform and Sangerbox tool to analyze the relationship between gene expression and neuroblastoma risk and prognosis through non-parametric testing and Kaplan-Meier analysis, respectively. We found the TET2 gene polymorphisms (rs10007915 G > C and rs7670522 A > C) and the combined 2-5 risk genotypes can significantly increase neuroblastoma risk. Stratification analysis showed that these significant associations were more prominent in certain subgroups. TET2 rs10007915 G > C and rs7670522 A > C are significantly associated with reduced expression of TET2 mRNA. Moreover, lower expression of TET2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastoma. The rs10007915 G > C and rs7670522 A > C are significantly related to the increased expression of inorganic pyrophosphatase 2 mRNA, and higher expression of PPA2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastomas. In summary, TET2 rs10007915 G > C and rs7670522 A > C significantly confer neuroblastoma susceptibility, and further research is needed to investigate the underlying mechanisms.

Lung cancer is the leading cause of cancer-related deaths worldwide. Long non-coding RNAs (lncRNAs) have emerged as key regulators of cancer development and progression, and as promising biomarkers for the diagnosis and prognosis of cancer. In this study, we identified a new lncRNA (LINC02159) that was upregulated in the tumor tissues and serum of non-small cell lung cancer (NSCLC) patients. We demonstrated that knockdown of LINC02159 inhibited NSCLC cell proliferation, migration, and invasion, but induced cell apoptosis and cell cycle arrest in vitro and retarded tumor growth in vivo, while overexpression of LINC02159 led to the opposite effect. We discovered that LINC02159 was highly correlated with cancer growth and metastasis-related pathways by using transcriptomic analysis and that YAP1 was a potential target gene of LINC02159. Mechanistically, LINC02159 bound to the Aly/REF export factor (ALYREF) to enhance the stability of YAP1 messenger RNA (mRNA) via m5C modification, which led to the overexpression of YAP1 and the activation of the Hippo and β-catenin signaling pathways in NSCLC cells. Rescue experiments showed that LINC01259 promoted NSCLC progression in a YAP1- and ALYREF-dependent manner. In conclusion, LINC02159 plays an oncogenic role in NSCLC progression by regulating ALYREF/YAP1 signaling, and it has the potential to be utilized as a diagnostic marker and therapeutic target for NSCLC.