Base-resolution m5C profiling across the mammalian transcriptome by bisulfite-free enzyme-assisted chemical labeling approach.

5-methylcytosine (m5C) is a prevalent RNA modification crucial for gene expression regulation. However, accurate and sensitive m5C sites identification remains challenging due to severe RNA degradation and reduced sequence complexity during bisulfite sequencing (BS-seq). Here, we report m5C-TAC-seq, a bisulfite-free approach combining TET-assisted m5C-to-f5C oxidation with selective chemical labeling, therefore enabling direct base-resolution m5C detection through pre-enrichment and C-to-T transitions at m5C sites. With m5C-TAC-seq, we comprehensively profiled the m5C methylomes in human and mouse cells, identifying a substantially larger number of confident m5C sites. Through perturbing potential m5C methyltransferases, we deciphered the responsible enzymes for most m5C sites, including the characterization of NSUN5's involvement in mRNA m5C deposition. Additionally, we characterized m5C dynamics during mESC differentiation. Notably, the mild reaction conditions and preservation of nucleotide composition in m5C-TAC-seq allow m5C detection in chromatin-associated RNAs. The accurate and robust m5C-TAC-seq will advance research into m5C methylation functional investigation.

SRSF2 plays an unexpected role as reader of m5C on mRNA, linking epitranscriptomics to cancer.

A common mRNA modification is 5-methylcytosine (m5C), whose role in gene-transcript processing and cancer remains unclear. Here, we identify serine/arginine-rich splicing factor 2 (SRSF2) as a reader of m5C and impaired SRSF2 m5C binding as a potential contributor to leukemogenesis. Structurally, we identify residues involved in m5C recognition and the impact of the prevalent leukemia-associated mutation SRSF2P95H. We show that SRSF2 binding and m5C colocalize within transcripts. Furthermore, knocking down the m5C writer NSUN2 decreases mRNA m5C, reduces SRSF2 binding, and alters RNA splicing. We also show that the SRSF2P95H mutation impairs the ability of the protein to read m5C-marked mRNA, notably reducing its binding to key leukemia-related transcripts in leukemic cells. In leukemia patients, low NSUN2 expression leads to mRNA m5C hypomethylation and, combined with SRSF2P95H, predicts poor outcomes. Altogether, we highlight an unrecognized mechanistic link between epitranscriptomics and a key oncogenesis driver.

Epitranscriptomic cytidine methylation of the hepatitis B viral RNA is essential for viral reverse transcription and particle production.

Epitranscriptomic RNA modifications have emerged as important regulators of the fate and function of viral RNAs. One prominent modification, the cytidine methylation 5-methylcytidine (m5C), is found on the RNA of HIV-1, where m5C enhances the translation of HIV-1 RNA. However, whether m5C functionally enhances the RNA of other pathogenic viruses remains elusive. Here, we surveyed a panel of commonly found RNA modifications on the RNA of hepatitis B virus (HBV) and found that HBV RNA is enriched with m5C as well as ten other modifications, at stoichiometries much higher than host messenger RNA (mRNA). Intriguingly, m5C is mostly found on the epsilon hairpin, an RNA element required for viral RNA encapsidation and reverse transcription, with these m5C mainly deposited by the cellular methyltransferase NSUN2. Loss of m5C from HBV RNA due to NSUN2 depletion resulted in a partial decrease in viral core protein (HBc) production, accompanied by a near-complete loss of the reverse transcribed viral DNA. Similarly, mutations introduced to remove the methylated cytidines resulted in a loss of HBc production and reverse transcription. Furthermore, pharmacological disruption of m5C deposition led to a significant decrease in HBV replication. Thus, our data indicate m5C methylations as a critical mediator of the epsilon elements' function in HBV virion production and reverse transcription, suggesting the therapeutic potential of targeting the m5C methyltransfer process on HBV epsilon as an antiviral strategy.

Construction of a prognostic model for lung adenocarcinoma based on m6A/m5C/m1A genes.

BACKGROUND: Developing a prognostic model for lung adenocarcinoma (LUAD) that utilizes m6A/m5C/m1A genes holds immense importance in providing precise prognosis predictions for individuals. METHODS: This study mined m6A/m5C/m1A-related differential genes in LUAD based on public databases, identified LUAD tumor subtypes based on these genes, and further built a risk prognostic model grounded in differential genes between subtypes. The immune status between high- and low-risk groups was investigated, and the distribution of feature genes in tumor immune cells was analyzed using single-cell analysis. Based on the expression levels of feature genes, a projection of chemotherapeutic and targeted drugs was made for individuals identified as high-risk. Ultimately, cell experiments were further verified. RESULTS: The 6-gene risk prognosis model based on differential genes between tumor subtypes had good predictive performance. Individuals classified as low-risk exhibited a higher (P < 0.05) abundance of infiltrating immune cells. Feature genes were mainly distributed in tumor immune cells like CD4+T cells, CD8+T cells, and regulatory T cells. Four drugs with relatively low IC50 values were found in the high-risk group: Elesclomol, Pyrimethamine, Saracatinib, and Temsirolimus. In addition, four drugs with significant positive correlation (P < 0.001) between IC50 values and feature gene expression were found, including Alectinib, Estramustine, Brigatinib, and Elesclomol. The low expression of key gene NTSR1 reduced the IC50 value of irinotecan. CONCLUSION: Based on the m6A/m5C/m1A-related genes in LUAD, LUAD patients were divided into 2 subtypes, and a m6A/m5C/m1A-related LUAD prognostic model was constructed to provide a reference for the prognosis prediction of LUAD.

5-Methylcytosine transferase NSUN2 drives NRF2-mediated ferroptosis resistance in non-small cell lung cancer.

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.

FMRP promotes transcription-coupled homologous recombination via facilitating TET1-mediated m5C RNA modification demethylation.

RNA modifications regulate a variety of cellular processes including DNA repair.The RNA methyltransferase TRDMT1 generates methyl-5-cytosine (m5C) on messen-ger RNA (mRNA) at DNA double-strand breaks (DSBs) in transcribed regions, pro-moting transcription-coupled homologous recombination (HR). Here, we identifiedthat Fragile X mental retardation protein (FMRP) promotes transcription-coupled HRvia its interaction with both the m5C writer TRDMT1 and the m5C eraser ten-eleventranslocation protein 1 (TET1). TRDMT1, FMRP, and TET1 function in a temporalorder at the transcriptionally active sites of DSBs. FMRP displays a higher affinity forDNA:RNA hybrids containing m5C-modified RNA than for hybrids without modifica-tion and facilitates demethylation of m5C by TET1 in vitro. Loss of either the chroma-tin- or RNA-binding domain of FMRP compromises demethylation of damage-inducedm5C in cells. Importantly, FMRP is required for R-loop resolving in cells. Due to unre-solved R-loop and m5C preventing completion of DSB repair, FMRP depletion or lowexpression leads to delayed repair of DSBs at transcriptionally active sites and sensitizescancer cells to radiation in a BRCA-independent manner. Together, ourfindings presentan m5C reader, FMRP, which acts as a coordinator between the m5C writer and eraserto promote mRNA-dependent repair and cell survival in cancer.

Development of m6A/m5C/m1A regulated lncRNA signature for prognostic prediction, personalized immune intervention and drug selection in LUAD.

Research indicates that there are links between m6A, m5C and m1A modifications and the development of different types of tumours. However, it is not yet clear if these modifications are involved in the prognosis of LUAD. The TCGA-LUAD dataset was used as for signature training, while the validation cohort was created by amalgamating publicly accessible GEO datasets including GSE29013, GSE30219, GSE31210, GSE37745 and GSE50081. The study focused on 33 genes that are regulated by m6A, m5C or m1A (mRG), which were used to form mRGs clusters and clusters of mRG differentially expressed genes clusters (mRG-DEG clusters). Our subsequent LASSO regression analysis trained the signature of m6A/m5C/m1A-related lncRNA (mRLncSig) using lncRNAs that exhibited differential expression among mRG-DEG clusters and had prognostic value. The model's accuracy underwent validation via Kaplan-Meier analysis, Cox regression, ROC analysis, tAUC evaluation, PCA examination and nomogram predictor validation. In evaluating the immunotherapeutic potential of the signature, we employed multiple bioinformatics algorithms and concepts through various analyses. These included seven newly developed immunoinformatic algorithms, as well as evaluations of TMB, TIDE and immune checkpoints. Additionally, we identified and validated promising agents that target the high-risk mRLncSig in LUAD. To validate the real-world expression pattern of mRLncSig, real-time PCR was carried out on human LUAD tissues. The signature's ability to perform in pan-cancer settings was also evaluated. The study created a 10-lncRNA signature, mRLncSig, which was validated to have prognostic power in the validation cohort. Real-time PCR was applied to verify the actual manifestation of each gene in the signature in the real world. Our immunotherapy analysis revealed an association between mRLncSig and immune status. mRLncSig was found to be closely linked to several checkpoints, such as IL10, IL2, CD40LG, SELP, BTLA and CD28, which could be appropriate immunotherapy targets for LUAD. Among the high-risk patients, our study identified 12 candidate drugs and verified gemcitabine as the most significant one that could target our signature and be effective in treating LUAD. Additionally, we discovered that some of the lncRNAs in mRLncSig could play a crucial role in certain cancer types, and thus, may require further attention in future studies. According to the findings of this study, the use of mRLncSig has the potential to aid in forecasting the prognosis of LUAD and could serve as a potential target for immunotherapy. Moreover, our signature may assist in identifying targets and therapeutic agents more effectively.

NSUN6-mediated 5-methylcytosine modification of NDRG1 mRNA promotes radioresistance in cervical cancer.

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.

The Emerging Roles of Cytosine-5 Methylation in mRNAs.

Recent studies have unequivocally confirmed the presence of 5-methylcytosine (m5C) in mammalian mRNAs while indicating significant functional roles for this internal base modification type. Here, a brief history of m5C epitranscriptome research and a discussion of the important ways in which the field may now progress is presented.

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

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

NSUN2 affects diabetic retinopathy progression by regulating MUC1 expression through RNA m5C methylation.

BACKGROUND: Diabetic retinopathy (DR) is the leading cause of blinding eye disease among working adults and is primarily attributed to the excessive proliferation of microvessels, which leads to vitreous hemorrhage and retinal traction, thereby significantly impairing patient vision. NSUN2-mediated RNA m5C methylation is implicated in various diseases, and in this investigation, we focused on elucidating the impact of NSUN2 on the regulation of the expression of the downstream gene MUC1, specifically through RNA m5C methylation, on the progression of DR. METHOD: Utilizing Microarray analysis, we examined patient vitreous fluid to pinpoint potential therapeutic targets for DR. Differential expression of NSUN2 was validated through qRT-PCR, Western blot, and immunofluorescence in human tissue, animal tissue, and cell model of DR. The relationship between NSUN2 and DR was explored in vitro and in vivo through gene knockdown and overexpression. Various techniques, such as MeRIP-qPCR and dot blot, were applied to reveal the downstream targets and mechanism of action of NSUN2. RESULTS: The levels of both NSUN2 and RNA m5C methylation were significantly elevated in the DR model. Knockdown of NSUN2 mitigated DR lesion formation both in vitro and in vivo. Mechanistically, NSUN2 promoted MUC1 expression by binding to the RNA m5C reader ALYREF. Knockdown of ALYREF resulted in DR lesion alterations similar to those observed with NSUN2 knockdown. Moreover, MUC1 overexpression successfully reversed a series of DR alterations induced by NSUN2 silencing. CONCLUSIONS: NSUN2 regulates the expression of MUC1 through ALYREF-mediated RNA m5C methylation, thereby regulating the progression of DR and providing a new option for the treatment of DR in the future.

Functional TET2 gene polymorphisms increase the risk of neuroblastoma in Chinese children.

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.

Comprehensive analysis of transcriptomic profiling of 5-methylcytosin modification in placentas from preeclampsia and normotensive pregnancies.

Increasing evidence suggests that RNA m5C modification and its regulators have been confirmed to be associated with the pathogenesis of many diseases. However, the distribution and biological functions of m5C in mRNAs of placental tissues remain unknown. we collected placentae from normotensive pregnancies (CTR) and preeclampsia patients (PE) to analyze the transcriptomic profiling of m5C RNA methylation through m5C RNA immunoprecipitation (UMI-MeRIP-Seq). we discovered that overall m5C methylation peaks were decreased in placental tissues from PE patients. And, 2844 aberrant m5C peaks were identified, of which respectively 1304 m5C peaks were upregulated and 1540 peaks were downregulated. The distribution of m5C peaks were mainly located in CDS (coding sequences) regions in placental tissues of both groups, but compared with the CTR group, the m5C peak in PE group before the stop code of CDS was significantly increased and even higher than the peak value after start code in CDS. Differentially methylated genes were mainly enriched in MAPK/cAMP signaling pathway. Moreover, the up-regulated genes with hypermethylated modification were enriched in the processes of hypoxia, inflammation/immune response. Finally, through analyzing the mRNA expression levels of m5C RNA methylation regulators, we found only DNMT3B and TET3 were significantly upregulated in PE samples than in control group. And they are not only negatively correlated with each other, but also closely related to those differentially expressed genes modified by differential methylation.Our findings provide new insights regarding alterations of m5C RNA modification into the pathogenic mechanisms of PE.

XGBoost framework with feature selection for the prediction of RNA N5-methylcytosine sites.

5-methylcytosine (m5C) is indeed a critical post-transcriptional alteration that is widely present in various kinds of RNAs and is crucial to the fundamental biological processes. By correctly identifying the m5C-methylation sites on RNA, clinicians can more clearly comprehend the precise function of these m5C-sites in different biological processes. Due to their effectiveness and affordability, computational methods have received greater attention over the last few years for the identification of methylation sites in various species. To precisely identify RNA m5C locations in five different species including Homo sapiens, Arabidopsis thaliana, Mus musculus, Drosophila melanogaster, and Danio rerio, we proposed a more effective and accurate model named m5C-pred. To create m5C-pred, five distinct feature encoding techniques were combined to extract features from the RNA sequence, and then we used SHapley Additive exPlanations to choose the best features among them, followed by XGBoost as a classifier. We applied the novel optimization method called Optuna to quickly and efficiently determine the best hyperparameters. Finally, the proposed model was evaluated using independent test datasets, and we compared the results with the previous methods. Our approach, m5C- pred, is anticipated to be useful for accurately identifying m5C sites, outperforming the currently available state-of-the-art techniques.

Epitranscriptomic orchestrations: Unveiling the regulatory paradigm of m6A, A-to-I editing, and m5C in breast cancer via long noncoding RNAs and microRNAs.

Breast cancer (BC) poses a persistent global health challenge, particularly in countries with elevated human development indices linked to factors such as increased life expectancy, education, and wealth. Despite therapeutic progress, challenges persist, and the role of epitranscriptomic RNA modifications in BC remains inadequately understood. The epitranscriptome, comprising diverse posttranscriptional modifications on RNA molecules, holds the potential to intricately modulate RNA function and regulation, implicating dysregulation in various diseases, including BC. Noncoding RNAs (ncRNAs), acting as posttranscriptional regulators, influence physiological and pathological processes, including cancer. RNA modifications in long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) add an extra layer to gene expression control. This review delves into recent insights into epitranscriptomic RNA modifications, such as N-6-methyladenosine (m6A), adenine-to-inosine (A-to-I) editing, and 5-methylcytosine (m5C), specifically in the context of lncRNA and miRNAs in BC, highlighting their potential implications in BC development and progression. Understanding this intricate regulatory landscape is vital for deciphering the molecular mechanisms underlying BC and identifying potential therapeutic targets.

Functional interplay within the epitranscriptome: Reality or fiction?

RNA modifications have recently emerged as an important regulatory layer of gene expression. The most prevalent and reversible modification on messenger RNA (mRNA), N6-methyladenosine, regulates most steps of RNA metabolism and its dysregulation has been associated with numerous diseases. Other modifications such as 5-methylcytosine and N1-methyladenosine have also been detected on mRNA but their abundance is lower and still debated. Adenosine to inosine RNA editing is widespread on coding and non-coding RNA and can alter mRNA decoding as well as protect against autoimmune diseases. 2'-O-methylation of the ribose and pseudouridine are widespread on ribosomal and transfer RNA and contribute to proper RNA folding and stability. While the understanding of the individual role of RNA modifications has now reached an unprecedented stage, still little is known about their interplay in the control of gene expression. In this review we discuss the examples where such interplay has been observed and speculate that with the progress of mapping technologies more of those will rapidly accumulate.

Genetic variants in m5C modification genes are associated with survival of patients with HBV-related hepatocellular carcinoma.

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.

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.

An Overview of Current Detection Methods for RNA Methylation.

Epitranscriptomic mechanisms, which constitute an important layer in post-transcriptional gene regulation, are involved in numerous cellular processes under health and disease such as stem cell development or cancer. Among various such mechanisms, RNA methylation is considered to have vital roles in eukaryotes primarily due to its dynamic and reversible nature. There are numerous RNA methylations that include, but are not limited to, 2'-O-dimethyladenosine (m6Am), N7-methylguanosine (m7G), N6-methyladenosine (m6A) and N1-methyladenosine (m1A). These biochemical modifications modulate the fate of RNA by affecting the processes such as translation, target site determination, RNA processing, polyadenylation, splicing, structure, editing and stability. Thus, it is highly important to quantitatively measure the changes in RNA methylation marks to gain insight into cellular processes under health and disease. Although there are complicating challenges in identifying certain methylation marks genome wide, various methods have been developed recently to facilitate the quantitative measurement of methylated RNAs. To this end, the detection methods for RNA methylation can be classified in five categories such as antibody-based, digestion-based, ligation-based, hybridization-based or direct RNA-based methods. In this review, we have aimed to summarize our current understanding of the detection methods for RNA methylation, highlighting their advantages and disadvantages, along with the current challenges in the field.

Functional polymorphisms of the TET1 gene increase the risk of neuroblastoma in Chinese children.

Common genetic mutations are absent in neuroblastoma, one of the most common childhood tumours. As a demethylase of 5-methylcytosine (m5C) modification, TET1 plays an important role in tumourigenesis and differentiation. However, the association between TET1 gene polymorphisms and susceptibility to neuroblastoma has not been reported. Three TET1 gene polymorphisms (rs16925541 A > G, rs3998860 G > A and rs12781492 A > C) in 402 Chinese patients with neuroblastoma and 473 cancer-free controls were assessed using TaqMan. Multivariate logistic regression analysis was used to evaluate the association between TET1 gene polymorphisms and susceptibility to neuroblastoma. The GTEx database was used to analyse the impact of these polymorphisms on peripheral gene expression. The relationship between gene expression and prognosis was analysed using Kaplan-Meier analysis with the R2 platform. We found that both rs3998860 G > A and rs12781492 A > C were significantly associated with increased neuroblastoma risk. Stratified analysis further showed that rs3998860 G > A and rs12781492 A > C significantly increased neuroblastoma risk in certain subgroups. In the combined risk genotype model, 1-3 risk genotypes significantly increased risk of neuroblastoma compared with the 0 risk genotype. rs3998860 G > A and rs12781492 A > C were significantly associated with increased STOX1 mRNA expression in adrenal and whole blood, and high expression of STOX1 mRNA in adrenal and whole blood was significantly associated with worse prognosis. In summary, TET1 gene polymorphisms are significantly associated with increased neuroblastoma risk; further research is required for the potential mechanism and therapeutic prospects in neuroblastoma.