The role of m6A RNA methylation in human cancer

The RNA methyltransferase NSUN6 suppresses pancreatic cancer development by regulating cell proliferation

BACKGROUND

Pancreatic cancer (PC) is one of the most lethal solid malignancies in the world due to its excessive cell proliferation and aggressive metastatic features. Emerging evidences revealed the importance of posttranscriptional modifications of RNAs in PC progression. However, knowledge about the 5-methylcytosine (m5C) RNA modification in PC is still extremely limited. In this study, we attempted to explore the expression changes and clinical significances of 12 known m5C-related genes among PC patients.

METHODS

A total of 362 normal and 382 tumor specimens from PC patients were examined for candidate m5C-related gene and protein expression by using quantitative PCR (qPCR) and immunohistochemistry (IHC). The proliferation rate of PC cells was detected by MTS assay. Xenograft mouse models were used to assess the role of NSUN6 in PC tumor formation.

FINDINGS

Through analyzing the four Gene Expression Omnibus (GEO) databases, six m5C-related genes shown significant and consistent alterations were selected for further examination in our 3 independent PC cohorts. Finally, we identified the reduction of NSUN6 as a common feature of all PC sample sets examined. NSUN6 expression correlated with clinicopathologic parameters including T stage, and Ki67⁺ cell rate. Further assessing the transcriptional profiles of 50 PC tissues, we found biological processes associated with cell proliferation like cell cycle and G2M checkpoint were enriched in NSUN6 lower expression group. Helped by in vitro PC cell lines and in vivo xenograft mouse models, we confirmed the role of NSUN6 in regulating cell proliferation and PC tumor growth. Last but also importantly, we also show the good performance of NSUN6 in evaluating tumor recurrence and survival among PC patients.

INTERPRETATION

Our data suggested that NSUN6 is an important factor involved in regulating cell proliferation of PC, and highlights the potential of novel m5C-based clinical modalities as a therapeutic approach in PC patients.

FUNDING

This study was supported by the National Natural Science Foundation of China (Grant Nos. 81803014, 81802424, and 81802911).

Identification and Verification of Molecular Subtypes with Enhanced Immune Infiltration Based on m6A Regulators in Cutaneous Melanoma

Background

As the most aggressive type of skin cancer, cutaneous melanoma (CM) is experiencing a rapidly rising mortality in recent years. Exploring potential prognostic biomarkers or mechanisms of disease progression therefore has a great significance for CM. The purpose of this study was to identify genetic markers and prognostic performance of N6-methyladenosine (m6A) regulators in CM.

Method

Gene expression profiles, copy number variation (CNV), and single nucleotide polymorphism (SNP) data of patients were obtained from The Cancer Genome Atlas (TCGA) database.

Results

Genomic variation and association analysis of gene expressions revealed a high degree of genomic variation in the presence of m6A-regulated genes. m6A patients with high-frequency genomic variants in the regulatory gene tended to develop a worse prognosis (p < 0.01). Unsupervised cluster analysis of the expression profiles of m6A-regulated genes identified three clinically distinct molecular subtypes, including degradation-enhanced subgroup and immune-enhanced subgroup, with significant prognostic differences (p = 0.046). A novel prognostic signature, which was established according to m6A-related characteristic genes identified through genome-wide expression spectrum, could effectively identify samples with poor prognosis and enhanced immune infiltration, and the effectiveness was also verified in the dataset of the chip.

Conclusion

We identified genetic changes in the m6A regulatory gene in CM and related survival outcomes. The findings of this study provide new insights into the epigenetic understanding of m6A in CM.

The pan-cancer analysis of the two types of uterine cancer uncovered clinical and prognostic associations with m6A RNA methylation regulators

The role of m6A RNA methylation modification in uterine cancer has not been studied until now. We explored the relationship between m6A regulators and clinical characteristics and prognosis in uterine corpus endometrial carcinoma (UCEC) and uterine carcinosarcoma (UCS) with the data from the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx). We found that several regulators were up-regulated or down-regulated in the two types of cancer, and identified two cluster subgroups with statistically significant differences in pathological grade, age and survival rate. Multivariate Cox regression analysis showed that methyltransferase-like 16 (METTL16) had a low hazard ratio in UCEC. We used several regulators to construct a risk signature and divided tumor patients into high-risk and low-risk groups, and found that the high-risk group had significantly lower survival rates. Independent prognostic analysis showed that the insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) was a pan-prognostic regulator of uterine cancer. This result was further verified in the Gene Expression Omnibus (GEO) database. Based on above results, we conducted gene-ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to further reveal a potential mechanism for m6A RNA methylation regulators. We found that IGF2BP1 was enriched in gene expression (GO:0010467), poly(A) RNA binding (GO:0044822) and RNA binding (GO:0003723) pathways. KEGG analysis showed that IGF2BP1 was enriched in microRNAs in the cancer (hsa05206) pathway. Our study systematically elucidated the relationship between m6A RNA methylation regulators and uterine cancer and constructed the risk signature that can predict the prognosis and clinicopathological characteristics of uterine cancer.

ALKBH5 Gene Polymorphisms and Hepatoblastoma Susceptibility in Chinese Children

Incidence of hepatoblastoma has been increasing, but the causes of this disease remain unclear. Some studies have suggested that abnormal expressions of ALKBH5 gene are associated with multiple cancers. This study aims to test the hypothesis that hepatoblastoma risk may be modulated by genetic polymorphisms in ALKBH5 gene based on genotyped data from samples of 328 cases and 1476 controls enrolled from eight hospitals in China. We used TaqMan assay to genotype ALKBH5 gene single nucleotide polymorphisms (SNPs) rs1378602G > A and rs8400G > A. We calculated the odds ratios (ORs) and P values using logistic regression models to estimate the association between hepatoblastoma risk and ALKBH5 gene SNPs. We found the rs1378602G > A and rs8400G > A could not impact hepatoblastoma risk in single or combined analysis. Stratified analysis revealed that subjects with the rs8400 AA genotype are prone to getting hepatoblastoma in the clinical stage III + IV subgroup (adjusted OR = 1.93, 95% CI = 1.20-3.10, P=0.007), when compared to those with GG/GA genotype. False-positive report probability validated the reliability of the significant results. Preliminary functional annotations revealed that rs8400 A is correlated with increased expression of ALKBH5 gene in the expression quantitative trait locus (eQTL) analysis. In all, our investigation presents evidence of a weak impact of ALKBH5 gene polymorphisms on hepatoblastoma risk, using the largest hepatoblastoma sample size. These findings shed some light on the genetic basis of hepatoblastoma, implicating the role of ALKBH5 gene polymorphisms in the etiology of hepatoblastoma.

N6-Methyladenosine RNA Methylation Regulator-Related Alternative Splicing (AS) Gene Signature Predicts Non-Small Cell Lung Cancer Prognosis

Aberrant N6-methyladenosine (m6A) RNA methylation regulatory genes and related gene alternative splicing (AS) could be used to predict the prognosis of non-small cell lung carcinoma. This study focused on 13 m6A regulatory genes (METTL3, METTL14, WTAP, KIAA1429, RBM15, ZC3H13, YTHDC1, YTHDC2, YTHDF1, YTHDF2, HNRNPC, FTO, and ALKBH5) and expression profiles in TCGA-LUAD (n = 504) and TCGA-LUSC (n = 479) datasets from the Cancer Genome Atlas database. The data were downloaded and bioinformatically and statistically analyzed, including the gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. There were 43,948 mRNA splicing events in lung adenocarcinoma (LUAD) and 46,020 in lung squamous cell carcinoma (LUSC), and the data suggested that m6A regulators could regulate mRNA splicing. Differential HNRNPC and RBM15 expression was associated with overall survival (OS) of LUAD and HNRNPC and METTL3 expression with the OS of LUSC patients. Furthermore, the non-small cell lung cancer prognosis-related AS events signature was constructed and divided patients into high- vs. low-risk groups using seven and 14 AS genes in LUAD and LUSC, respectively. The LUAD risk signature was associated with gender and T, N, and TNM stages, but the LUSC risk signature was not associated with any clinical features. In addition, the risk signature and TNM stage were independent prognostic predictors in LUAD and the risk signature and T stage were independent prognostic predictors in LUSC after the multivariate Cox regression and receiver operating characteristic analyses. In conclusion, this study revealed the AS prognostic signature in the prediction of LUAD and LUSC prognosis.

m6A Regulators in Human Adipose Tissue - Depot-Specificity and Correlation With Obesity

BACKGROUND

N⁶-methyladenosine (m6A) is one of the most abundant post-transcriptional modifications on mRNA influencing mRNA metabolism. There is emerging evidence for its implication in metabolic disease. No comprehensive analyses on gene expression of m6A regulators in human adipose tissue, especially in paired adipose tissue depots, and its correlation with clinical variables were reported so far. We hypothesized that inter-depot specific gene expression of m6A regulators may differentially correlate with clinical variables related to obesity and fat distribution.

METHODS

We extracted intra-individually paired gene expression data (omental visceral adipose tissue (OVAT) N=48; subcutaneous adipose tissue (SAT) N=56) of m6A regulators from an existing microarray dataset. We also measured gene expression in another sample set of paired OVAT and SAT (N=46) using RT-qPCR. Finally, we extracted existing gene expression data from peripheral mononuclear blood cells (PBMCs) and single nucleotide polymorphisms (SNPs) in METTL3 and YTHDF3 from genome wide data from the Sorbs population (N=1049). The data were analysed for differential gene expression between OVAT and SAT; and for association with obesity and clinical variables. We further tested for association of SNP markers with gene expression and clinical traits.

RESULTS

In adipose tissue we observed that several m6A regulators (WTAP, VIRMA, YTHDC1 and ALKBH5) correlate with obesity and clinical variables. Moreover, we found adipose tissue depot specific gene expression for METTL3, WTAP, VIRMA, FTO and YTHDC1. In PBMCs, we identified ALKBH5 and YTHDF3 correlated with obesity. Genetic markers in METTL3 associate with BMI whilst SNPs in YTHDF3 are associated with its gene expression.

CONCLUSIONS

Our data show that expression of m6A regulators correlates with obesity, is adipose tissue depot-specific and related to clinical traits. Genetic variation in m6A regulators adds an additional layer of variability to the functional consequences.

The contribution of YTHDF2 gene rs3738067 A>G to the Wilms tumor susceptibility

YTHDF2 is responsible for maintaining the dynamic N⁶-methyladenosine (m⁶A) modification balance and influences a variety of cancers. We tested whether YTHDF2 gene rs3738067 A>G polymorphism is related to Wilms tumor by genotyping samples of Chinese children (450 cases and 1317 controls). However, the rs3738067 A>G polymorphism showed no statistical significance with Wilms tumor susceptibility. Stratification analysis also revealed that there was no remarkable association of rs3738067 variant AG/GG genotype with Wilms tumor risk in every subgroup (age, gender, and clinical stages). In all, the results indicated YTHDF2 gene rs3738067 A>G polymorphism could not alter Wilms tumor risk significantly.

Clinicopathological and immunological characterization of RNA m6 A methylation regulators in ovarian cancer

BACKGROUND

N⁶ -methyladenosine (m⁶ A) modification is one of the critical gene regulatory mechanisms implicated in cancer biology. However, the roles of m⁶ A regulators in ovarian cancer are still poorly understood.

METHODS

We integrated multiple databases including Gene Expression Omnibus (GEO), ROC Plotter, Kaplan-Meier Plotter, and Tumor Immune Estimation Resource (TIMER) to explore clinicopathological significance of m⁶ A regulators in ovarian cancer.

RESULTS

We showed that alterations in the expression of m⁶ A regulators were related to the malignancy and poor prognosis of ovarian cancer. We found decreased YTHDC1 and increased RBM15 expressions were associated with ovarian cancer cell metastases and HNRNPC was a predictor of paclitaxel resistance. Moreover, dysregulated m⁶ A regulators were enriched in the activation of cancer-related pathways. Our results further demonstrated that the level of immune cell infiltration and the expression of various immune gene markers were closely associated with the expressions of specific m⁶ A regulators (RBM15B, ZC3H13, YTHDF1, and IGF2BP1).

CONCLUSIONS

Our study establishes a new prognostic profile of ovarian cancer patients based on m⁶ A regulators, and highlights the potential roles of m⁶ A regulators in ovarian cancer development.

MODOMICS: An Operational Guide to the Use of the RNA Modification Pathways Database

MODOMICS is an established database of RNA modifications that provides comprehensive information concerning chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. This chapter covers the resources available on MODOMICS web server and the basic steps that can be undertaken by the user to explore them. MODOMICS is available at http://www.genesilico.pl/modomics .

MiR-9-1 Suppresses Cell Proliferation and Promotes Apoptosis by Targeting UHRF1 in Lung Cancer

Lung cancer is listed as the most common reason for cancer-related death all over the world despite diagnostic improvements and the development of chemotherapy and targeted therapies. MicroRNAs control both physiological and pathological processes including development and cancer. A microRNA-9 to 1 (miR-9 to 1) overexpression model in lung cancer cell lines was established and miR-9 to 1 was found to significantly suppress the proliferation rate in lung cancer cell lines, colony formation in vitro, and tumorigenicity in nude mice of A549 cells. Ubiquitin-like containing PHD and RING finger domains 1 (UHRF1) was then identified to direct target of miR-9 to 1. The inhibition of UHRF1 by miR-9 to 1 causes G1 arrest and p15, p16, and p21 were re-expressed in miR-9 to 1 group in mRNA level and protein level. Silence of UHRF1 expression in A549 cells resulted in the similar re-expression of p15, p16, p21 which is similar with miR-9 to 1 infection. Therefore, we concluded that UHRF1 is a new target for miR-9 to 1 to suppress cell proliferation by re-expression of tumor suppressors p15, p16, and p21 mediated by UHRF1.

m6A regulator-based methylation modification patterns characterized by distinct tumor microenvironment immune profiles in colon cancer

Recent studies have highlighted the biological significance of RNA N⁶-methyladenosine (m⁶A) modification in tumorigenicity and progression. However, it remains unclear whether m⁶A modifications also have potential roles in immune regulation and tumor microenvironment (TME) formation. Methods: In this study, we curated 23 m⁶A regulators and performed consensus molecular subtyping with NMF algorithm to determine m⁶A modification patterns and the m⁶A-related gene signature in colon cancer (CC). The ssGSEA and CIBERSORT algorithms were employed to quantify the relative infiltration levels of various immune cell subsets. An PCA algorithm based m⁶Sig scoring scheme was used to evaluate the m⁶A modification patterns of individual tumors with an immune response. Results: Three distinct m6A modification patterns were identified among 1307 CC samples, which were also associated with different clinical outcomes and biological pathways. The TME characterization revealed that the identified m⁶A patterns were highly consistent with three known immune profiles: immune-inflamed, immune-excluded, and immune-desert, respectively. Based on the m⁶Sig score, which was extracted from the m⁶A-related signature genes, CC patients can be divided into high and low score subgroups. Patients with lower m⁶Sig score was characterized by prolonged survival time and enhanced immune infiltration. Further analysis indicated that lower m⁶Sig score also correlated with greater tumor mutation loads, PD-L1 expression, and higher mutation rates in SMGs (e.g., PIK3CA and SMAD4). In addition, patients with lower m⁶Sig scores showed a better immune responses and durable clinical benefits in three independent immunotherapy cohorts. Conclusions: This study highlights that m⁶A modification is significantly associated with TME diversity and complexity. Quantitatively evaluating the m⁶A modification patterns of individual tumors will strengthen our understanding of TME characteristics and promote more effective immunotherapy strategies.

Blocking the IGF2BP1-promoted glucose metabolism of colon cancer cells via direct de-stabilizing mRNA of the LDHA enhances anticancer effects

Colorectal cancer (CRC) is a commonly diagnosed cancer with poor prognosis and high mortality rate. Hyperthermia (HT) is an adjunctive therapy to enhance the antitumor effects of traditional chemo- or radio- therapy. Here, we report that a cluster of essential regulator genes and speed-limit enzymes of glucose metabolism were significantly elevated under HT from a glucose metabolism PCR array analysis. Under low glucose supply or glucose metabolism inhibition, CRC cells displayed increased sensitivity to HT treatments. By transcript sequencing from the established HT resistant (HTR) colon cancer cell line LoVo HTR, we observed that IGF2BP1, an RNA-binding protein, was significantly upregulated in HTR cells compared with parental cells. Furthermore, LDHA mRNA was identified as an IGF2BP1 direct target. An RNA immunoprecipitation assay and RNA pull-down assay consistently illustrated IGF2BP1 specifically bonds to the 3′ UTR of LDHA mRNA, leading to enhanced stability of LDHA mRNA. Finally, we demonstrated that inhibiting the IGF2BP1-promoted glycolysis sensitized colon cancer cells to HT treatment via both in vitro and in vivo experiments. Our findings suggest that targeting the IGF2BP1-LDHA-glycolysis pathway might be a promising therapeutic approach to enhance the anti-cancer effects of HT treatment.

Long non‑coding RNA CASC15 facilitates esophageal squamous cell carcinoma tumorigenesis via decreasing SIM2 stability via FTO‑mediated demethylation

Long non-coding RNAs (lncRNAs) are involved in the regulation of esophageal squamous cell carcinoma (ESCC) progression. However, the function and mechanism of lncRNA cancer susceptibility candidate 15 (CASC15) are poorly defined. In the present study, tumor and normal adjacent tissues were collected from 45 patients with ESCC. Expression levels of CASC15, fat mass and obesity-associated (FTO) protein and single-minded 2 (SIM2) were examined via reverse transcription-quantitative PCR and western blot assays. Cell proliferation and apoptosis were evaluated via MTT, flow cytometry and caspase-3 activity assays, respectively. Additionally, an ESCC mouse xenograft model was used to assess the function of CASC15 in vivo. The interaction between FTO and CASC15/SIM2 was analyzed via RNA immunoprecipitation and RNA pull-down assays. The results revealed that CASC15 expression was elevated in ESCC tissues, and patients with ESCC exhibiting high CASC15 expression had a poor prognosis. CASC15-knockdown inhibited ESCC cell proliferation and facilitated apoptosis. Additionally, CASC15-knockdown decreased the growth of ESCC xenograft tumors. CASC15 decreased SIM2 stability via FTO-mediated demethylation. Additionally, FTO loss markedly weakened CASC15-mediated pro-proliferative and anti-apoptotic effects in ESCC cells. SIM2 downregulation weakened the effect of CASC15-knockdown on cell proliferation and inhibited the increase of the apoptotic rate and caspase-3 activity induced by CASC15 depletion in ESCC cells. In conclusion, CASC15 promoted ESCC tumorigenesis by decreasing SIM2 stability via FTO-mediated demethylation.

Head and Neck Squamous Cell Carcinoma: Epigenetic Landscape

Head and neck squamous carcinoma (HNSCC) constitutes the sixth most prevalent cancer worldwide. The molecular pathogenesis of HNSCC includes disorders in cell cycle, intercellular signaling, proliferation, squamous cell differentiation and apoptosis. In addition to the genetic mutations, changes in HNSCC are also characterized by the accumulation of epigenetic alterations such as DNA methylation, histone modifications, non-coding RNA activity and RNA methylation. In fact, some of them may promote cancer formation and progression by controlling the gene expression machinery, hence, they could be used as biomarkers in the clinical surveillance of HNSCC or as targets for therapeutic strategies. In this review, we focus on the current knowledge regarding epigenetic modifications observed in HNSCC and its predictive value for cancer development.

The m6A RNA Demethylase ALKBH5 Promotes Radioresistance and Invasion Capability of Glioma Stem Cells

Simple Summary

Glioblastoma stem cells (GBMSCs), which are particularly radio-resistant and invasive, are responsible for the high recurrence of glioblastoma (GBM). Therefore, there is a real need for a better understanding of the mechanisms involved in these processes and to identify new factors that might be targeted to radiosensitize GBMSC and decrease their invasive capability. Here, we report that the m6A RNA demethylase ALKBH5, which is overexpressed in GBMSCs, promotes their radioresistance by controlling the homologous repair. ALKBH5 was also involved in GBMSC invasion. These data suggest that ALKBH5 inhibition might be a novel approach to radiosensitize GBMSCs and to overcome their invasiveness.

Abstract

Recurrence of GBM is thought to be due to GBMSCs, which are particularly chemo-radioresistant and characterized by a high capacity to invade normal brain. Evidence is emerging that modulation of m6A RNA methylation plays an important role in tumor progression. However, the impact of this mRNA modification in GBM is poorly studied. We used patient-derived GBMSCs to demonstrate that high expression of the RNA demethylase, ALKBH5, increases radioresistance by regulating homologous recombination (HR). In cells downregulated for ALKBH5, we observed a decrease in GBMSC survival after irradiation likely due to a defect in DNA-damage repair. Indeed, we observed a decrease in the expression of several genes involved in the HR, including CHK1 and RAD51, as well as a persistence of γ-H2AX staining after IR. We also demonstrated in this study that ALKBH5 contributes to the aggressiveness of GBM by favoring the invasion of GBMSCs. Indeed, GBMSCs deficient for ALKBH5 exhibited a significant reduced invasion capability relative to control cells. Our data suggest that ALKBH5 is an attractive therapeutic target to overcome radioresistance and invasiveness of GBMSCs.

Gene Alterations of N6-Methyladenosine (m6A) Regulators in Colorectal Cancer: A TCGA Database Study

N6-methyladenosine (m⁶A) plays an important role in many cancers. However, few studies have examined the role of m6A in colorectal CRC. To examine the effect of m6A on CRC, we studied the genome of 591 CRC cases from The Cancer Genome Atlas (TCGA). The relationship between the messenger RNA (mRNA) expression, copy number variation (CNVs), and mutations of m6A “Writers,” “Readers,” and “Erasers,” prognosis, immune cell infiltration, and genetic mutations in CRC cases were analyzed. CNVs and mutations were found in thirteen m6A regulators. As expected, gain and amplification of m6A regulators increased the mRNA expression of these regulators, while deletion led to reduction in the mRNA expression. Moreover, CNVs and mutation of these regulators were significantly associated with APC, TP53, and microsatellite instability (MSI) status (p < 0.001, p < 0.001, and p = 0.029, respectively). CNVs of m6A regulators also correlated with inferred immune cell infiltration in CRC tissues, especially in colon tissues. Additionally, alterations of RBM15, YTHDF2, YTHDC1, YTHDC2, and METTL14 genes were related to the worse overall survival and disease-free survival (DFS) of CRC patients. Specifically, the deletion status of “Writers” was also correlated to the DFS of CRC patients (p = 0.02). Gene set enrichment analysis found that FTO was involved in mRNA 3′ end processing, polyubiquitin binding, and RNA polymerase promoter elongation, while YTHDC1 was related to interferon-alpha and gamma response. In conclusion, a novel relationship was identified between CNVs and mutations of m6A regulators with prognosis and inferred immune function of CRC. These findings will improve the understanding of the relationship of m6A in CRC.

Emerging roles of RNA methylation in gastrointestinal cancers

RNA methylation has emerged as a fundamental process in epigenetic regulation. Accumulating evidences indicate that RNA methylation is essential for many biological functions, and its dysregulation is associated with human cancer progression, particularly in gastrointestinal cancers. RNA methylation has a variety of biological properties, including N6-methyladenosine (m6A), 2-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), 5-methylcytosine (m5C) and 7-methyl guanosine (m7G). Dynamic and reversible methylation on RNA is mediated by RNA modifying proteins called "writers" (methyltransferases) and "erasers" (demethylases). "Readers" (modified RNA binding proteins) recognize and bind to RNA methylation sites, which influence the splicing, stability or translation of modified RNAs. Herein, we summarize the biological functions and mechanisms of these well-known RNA methylations, especially focusing on the roles of m6A in gastrointestinal cancer development.

DNA adenine methylation in eukaryotes: Enzymatic mark or a form of DNA damage?

6-methyladenine (6mA) is fairly abundant in nuclear DNA of basal fungi, ciliates and green algae. In these organisms, 6mA is maintained near transcription start sites in ApT context by a parental-strand instruction dependent maintenance methyltransferase and is positively associated with transcription. In animals and plants, 6mA levels are high only in organellar DNA. The 6mA levels in nuclear DNA are very low. They are attributable to nucleotide salvage and the activity of otherwise mitochondrial METTL4, and may be considered as a price that cells pay for adenine methylation in RNA and/or organellar DNA. Cells minimize this price by sanitizing dNTP pools to limit 6mA incorporation, and by converting 6mA that has been incorporated into DNA back to adenine. Hence, 6mA in nuclear DNA should be described as an epigenetic mark only in basal fungi, ciliates and green algae, but not in animals and plants.

Exploring diagnostic m6A regulators in endometriosis

Endometriosis is an estrogen-dependent inflammatory disorder, usually causing infertility, pelvic pain, and ovarian masses. This study intended to investigate the implication of N6-methyladenosine (m6A) regulators in endometriosis. We acquired 34 normal, 127 eutopic, and 46 ectopic, samples of endometrium from the Gene Expression Omnibus (GSE7305, GSE7307, GSE51981) database and the Array-express (E-MTAB-694) database. These samples were then used to profile the expression of 20 m6A regulators in endometriosis. The results indicated that most dysregulated (19/20) m6A regulators were significantly downregulated in eutopic vs. normal endometrium and also significantly downregulated in ectopic vs. eutopic endometrium. Several dysregulated m6A regulators were common to both contrast matrices: METTL3, YTHDF2, YTHDF3, HNRNPA2B1, HNRNPC, and FTO. Both HNRNPA2B1 and HNRNPC were associated with the severity of endometriosis in eutopic samples, and also exhibited diagnostic potential for endometriosis. HNRNPA2B1 and HNRNPC may influence immune pathways and the infiltration of immune cells in endometriosis. Abnormalities in the gene transcription factors network associated with endometriosis might affect the expression of HNRNPA2B1 and HNRNPC. In conclusion, we observed significant dysregulation of m6A regulators in endometriosis, and found that HNRNPA2B1 and HNRNPC might correlate with the immune response and serve as useful diagnostic biomarkers for endometriosis.

Clinical and prognostic pan-cancer analysis of m6A RNA methylation regulators in four types of endocrine system tumors

N6-methyladenosine (m6A), internal modification of mRNA, has recently been reported to be an important regulatory mechanism affecting tumor proliferation. However, its role in endocrine system tumors is poorly understood. We obtained datasets for four types tumors from the TCGA database, analyzed the GTEx database as a supplement to the control group, and used “Perl” and “R” software to analyze the datasets. Then we differentiated the expression level, used it to cluster consensus. Besides, we established lasso regression model to screen variables, used univariate and multivariate cox analyses to explore the independent risk factors associated with cancer prognosis. The results indicated that except for WTAP, the expression level of METTL3 was negatively correlated with other genes. The expression level of WTAP and METTL16 was positively correlated with overall survival (OS). Moreover, we found that different clinical subtypes of adrenal cortical carcinoma had significant differences in survival status, histologic grading, pathological T grade, and OS. Furthermore, different clinical subtypes of thyroid carcinoma had significant differences in histologic grading and pathological T grade. The differential expression of m6A regulatory genes is closely associated with the presence of endocrine-system-related tumors, and risk scores can be used to assess prognosis.

The m6A methylation landscape stratifies hepatocellular carcinoma into 3 subtypes with distinct metabolic characteristics

Objective

Epigenetic aberration plays an important role in the development and progression of hepatocellular carcinoma (HCC). However, the alteration of RNA N6-methyladenosine (m6A) modifications and its role in HCC progression remain unclear. We therefore aimed to provide evidence using bioinformatics analysis.

Methods

We comprehensively analyzed the m6A regulator modification patterns of 605 HCC samples and correlated them with metabolic alteration characteristics. We elucidated 390 gene-based m6A-related signatures and defined an m6Ascore to quantify m6A modifications. We then assessed their values for predicting prognoses and therapeutic responses in HCC patients.

Results

We identified 3 distinct m6A modification patterns in HCC, and each pattern had distinct metabolic characteristics. The evaluation of m6A modification patterns using m6Ascores could predict the prognoses, tumor stages, and responses to sorafenib treatments of HCC patients. A nomogram based on m6Ascores showed high accuracy in predicting the overall survival of patients. The area under the receiver operating characteristic curve of predictions of 1, 3, and 5-year overall survivals were 0.71, 0.69, and 0.70 in the training cohort, and in the test cohort it was 0.74, 0.75, and 0.71, respectively. M6Acluster C1, which corresponded to hypoactive mRNA methylation, lower expression of m6A regulators, and a lower m6Ascore, was characterized by metabolic hyperactivity, lower tumor stage, better prognosis, and lower response to sorafenib treatment. In contrast, m6Acluster C3 was distinct in its hyperactive mRNA methylations, higher expression of m6A regulators, and higher m6Ascores, and was characterized by hypoactive metabolism, advanced tumor stage, poorer prognosis, and a better response to sorafenib. The m6Acluster, C2, was intermediate between C1 and C3.

Conclusions

HCCs harbored distinct m6A regulator modification patterns that contributed to the metabolic heterogeneity and diversity of HCC. Development of m6A gene signatures and the m6Ascore provides a more comprehensive understanding of m6A modifications in HCC, and helps predict the prognosis and treatment response.

Systematic Investigation of mRNA N 6-Methyladenosine Machinery in Primary Prostate Cancer

Background

Appreciable findings have pointed out pivotal roles of N⁶-methyladenosine (m⁶A) machinery in cancer onset and progression. However, limited efforts have been directed towards relevant research in the prostate cancer area.

Methods

A PubMed search was conducted to acquire components of the mRNA m⁶A machinery. Multiomics integration was performed to systematically investigate the mRNA m⁶A machinery in primary prostate cancer. Furthermore, RNA interference assays of two prognostic m⁶A readers EIF3D and HNRNPA2B1 were conducted to explore m⁶A dependence of their functions in prostate cancer cell proliferation and migration.

Results

A total of 41 mRNA m⁶A regulators have been identified to date. A small degree of copy number aberrations and an extremely low frequency of somatic mutations were observed in the regulators across prostate tumors. Enrichment of CpG sites and extensive changes of DNA methylation in the m⁶A machinery were also found. Impact of copy number variation on m⁶A regulator expression was stronger than that of DNA methylation disturbance. Furthermore, our study identified a set of m⁶A regulators related to clinical features and/or survival which were largely m⁶A-binding proteins. The translation initiation factor subunit EIF3D and the splicing factor HNRNPA2B1 can be independent prognostic factors which may contribute to retardation and promotion of cancer progression, respectively, through affecting cancer-related processes such as cell cycle. Moreover, in vitro assays demonstrated that m⁶A impacted the EIF3D and HNRNPA2B1 roles in proliferation and migration of prostate cancer cells.

Conclusions

Our report systematically described molecular features of the mRNA m⁶A machinery and their potential roles in primary prostate cancer. Knowledge gained from this work may pave the way for further studies on the m⁶A system in prostate cancer.

N6-methyladenine modification in noncoding RNAs and its function in cancer

Non-coding RNAs are the main component of the extensive transcription results of the mammalian genome. They are not transcribed into proteins but play critical roles in regulating multiple biological processes and affecting cancer progression. m6A modification is one of the most abundant internal RNA modification of mammalian cells, and it involves almost all aspects of RNA metabolism. Recent research revealed tight correlations between m6A modification and ncRNAs and indicated the interaction between m6A and ncRNAs act a pivotal part in the development of cancer. The correlation between m6A modification and ncRNAs provides a new perspective for exploring the potential regulatory mechanism of tumor gene expression, and suggest that m6A modification and ncRNAs may be important prognostic markers and therapeutic targets for multiple cancers. In this review, we summarize the potential regulatory mechanisms between m6A methylation and ncRNAs, highlighting how their relationship affects biological functions in cancer.

When the Smoke Clears m6A from a Y Chromosome-Linked lncRNA, Men Get an Increased Risk of Cancer

Long noncoding RNAs (lncRNA) have been implicated in many diseases, including cancer. Although these disease-associated effects have been mostly attributed to the ability of lncRNAs to function as regulatory noncoding transcripts, there is growing evidence that lncRNAs may also encode functional micropeptides. In the current issue of Cancer Research, Wu and colleagues report a micropeptide encoded by a Y chromosome-linked lncRNA that may explain the higher incidence of esophageal cancer in male smokers. Furthermore, this report provides broader insights related to the molecular epidemiology of male-dominant and smoking-driven cancers and may also help explain some cancer-related associations with mosaic Y chromosome loss.See related article by Wu et al., p. 2790.

FTO demethylates m6A modifications in HOXB13 mRNA and promotes endometrial cancer metastasis by activating the WNT signalling pathway

Although many studies have confirmed the relationship between obesity and endometrial cancer (EC), the molecular mechanism between obesity and EC progression has not been elucidated. Overexpression of fat mass and the obesity associated protein FTO leads to weight gain, although recently it has been discovered that FTO can serve as a demethylase which erases N6-methyladenosine (m6A) modification and regulates the metabolization of mRNAs. In this study, we found high expression of FTO in metastatic EC and that this action promote both metastasis and invasion in vivo and in vitro. Mechanistically, FTO can catalyse demethylation modification in 3'UTR region of HOXB13 mRNA, thereby abolishing m6A modification recognition with the YTHDF2 protein. Decreasing HOXB13 mRNA decay and increasing HOXB13 protein expression was accompanied by WNT signalling pathway activation and the expression of downstream proteins, leading to tumour metastasis and invasion. We also found the WNT signalling pathway inhibitor ICG-001 can block HOXB13 gene-induced tumour metastasis, therefore ICG-001 may be a promising molecular intervention. This study provides insight into the relationship between obesity and the pathogenesis of endometrial cancer while highlighting future areas of research.

The N6-Methyladenosine Features of mRNA and Aberrant Expression of m6A Modified Genes in Gastric Cancer and Their Potential Impact on the Risk and Prognosis

Although N6-methyladenosine (m6A) mRNA methylation is known to be closely related to tumor events, its role in carcinogenesis and the development of gastric cancer (GC) is not yet clear. The aim of this study was to identify common m6A features and novel aberrant expression of m6A modified genes in GC and to further explore their potential impact on risk and prognosis. Three paired GC and paracancerous (PCa) tissues were collected to perform an m6A sequencing by MeRIP-seq and microarray assays. The expression profile of m6A and mRNA were determined. Gene function note and enrichment analysis were performed, and protein–protein interaction networks of differentially m6A methylated genes (DMGs) were generated using the DAVID and STRING databases, respectively. Validation of the m6A related differentially expressed genes by matching TCGA and GTEx data and human tissues. Clinical and pathological correlation and survival analysis were performed by TCGA data. The m6A motif sequence GGACAR (R = U or A) C was the consensus in both GC and PCa tissues. m6A peaks were significantly related to different coordinates, however, for most samples, the end of the coding sequence (CDS) was more prominent than the start of CDS. The genes with higher levels of m6A in their mRNAs were mainly enriched in transcriptional misregulation in carcinogenesis pathways, whereas the genes with decreased methylation mainly regulated digestion and absorption of protein. There are genes with differential m6A modifications in GC and paired PCa tissues, and these genes are mainly enriched in transcriptional misregulation and digestion/absorption pathways. m6A-GC with the down- and up-regulated genes may play an important role in gastric carcinogenesis, which can affect the risk and prognosis in GC.

Prognostic Value of YTHDF2 in Clear Cell Renal Cell Carcinoma

m6A, the main form of mRNA modification, participates in regulating multiple normal and pathological biological events, especially in tumorigenesis. However, there is little known about the association of m6A-related genes with prognosis of clear cell renal cell cancer (ccRCC). Therefore, the prognostic value of m6A-related genes was investigated using Kaplan–Meier curves of overall survival (OS) with the log-rank test and Cox regression analysis. The differential expression of YTHDF2 mRNA in ccRCC and tumor-adjacent normal tissues and associated with clinicopathological characteristics was also analyzed. The alteration of cancer signaling pathways was screened by Gene Set Enrichment Analysis (GSEA). Univariate analysis showed that 15 m6A-related genes (including YTHDF2) were closely related to prognosis. Multivariate analysis further confirmed that YTHDF2 could serve as an independent prognostic factor for the OS of ccRCC patients (P < 0.001). Low-level expression of YTHDF2 had poor prognosis in ccRCC patients with lower tumor–node–metastasis (TNM) stage, age > 61, non-distant metastasis, non-lymph node metastasis, female gender, and higher histological grade (P < 0.05). Moreover, YTHDF2 expression in ccRCC tissues (N = 529) is significantly lower than that of tumor-adjacent normal tissues (N = 72, P = 0.0086). Furthermore, GSEA demonstrated that AKT/mTOR/GSK3 pathway, EIF4 pathway, CHREBP2 pathway, MET pathway, NFAT pathway, FAS pathway, EDG1 pathway, and CTCF pathway are altered in tumors with high YTHDF2 expression. Taken together, our results demonstrated that YTHDF2 (an m6A-related gene) could serve as a potential prognostic biomarker of ccRCC, and targeting epigenetic modification may be a novel therapeutic strategy for the treatment of ccRCC.

Epigenetic crosstalk between hypoxia and tumor driven by HIF regulation

Hypoxia is the major influence factor in physiological and pathological courses which are mainly mediated by hypoxia-inducible factors (HIFs) in response to low oxygen tensions within solid tumors. Under normoxia, HIF signaling pathway is inhibited due to HIF-α subunits degradation. However, in hypoxic conditions, HIF-α is activated and stabilized, and HIF target genes are successively activated, resulting in a series of tumour-specific activities. The activation of HIFs, including HIF-1α, HIF-2α and HIF-3α, subsequently induce downstream target genes which leads to series of responses, the resulting abnormal processes or metabolites in turn affect HIFs stability. Given its functions in tumors progression, HIFs have been regarded as therapeutic targets for improved treatment efficacy. Epigenetics refers to alterations in gene expression that are stable between cell divisions, and sometimes between generations, but do not involve changes in the underlying DNA sequence of the organism. And with the development of research, epigenetic regulation has been found to play an important role in the development of tumors, which providing accumulating basic or clinical evidences for tumor treatments. Here, given how little has been reported about the overall association between hypoxic tumors and epigenetics, we made a more systematic review from epigenetic perspective in hope of helping others better understand hypoxia or HIF pathway, and providing more established and potential therapeutic strategies in tumors to facilitate epigenetic studies of tumors.

N6-methyladenosine METTL3 promotes cervical cancer tumorigenesis and Warburg effect through YTHDF1/HK2 modification

N6-methyladenosine (m6A) serves as the most common and conserved internal transcriptional modification. However, the roles of m6A on cervical cancer (CC) tumorigenesis are still unclear. Here, results indicated that METTL3 was significantly upregulated in CC tissue and cells, which was closely correlated with the lymph node metastasis and poor prognosis of CC patients. MeRIP-Seq analysis revealed the m6A profiles in CC cells. Functionally, METTL3 promoted the proliferation and Warburg effect (aerobic glycolysis) of CC cells. Mechanistically, METTL3 targeted the 3'-Untranslated Region (3'-UTR) of hexokinase 2 (HK2) mRNA. Moreover, METTL3 recruited YTHDF1, a m6A reader, to enhance HK2 stability. These findings demonstrated that METTL3 enhanced the HK2 stability through YTHDF1-mediated m6A modification, thereby promoting the Warburg effect of CC, which might promote a novel insight for the CC treatment.

Multiomics integrative analysis for gene signatures and prognostic values of m6A regulators in pancreatic adenocarcinoma: a retrospective study in The Cancer Genome Atlas project

N6-methyladenosine(m⁶A) is the most abundant post-transcriptional RNA modification in eukaryotes. However, little is known about its role in pancreatic adenocarcinoma (PAAD). The aim of our study was to identify gene signatures and prognostic values of m⁶A regulators in PAAD. Patients from 3 different datasets with complete genomic and transcriptomic sequencing data were enrolled. Survival analysis for different gene alterations was performed using log-rank tests and Cox regression model. The association between alteration of m⁶A regulators and clinicopathological characteristics was examined using chi-square test. Results showed a high frequency of copy number alterations (CNAs) of m⁶A regulatory genes in PAAD patients, but somatic mutations were rarely happened. CNAs and mutations of m⁶A regulatory genes was associated with patient's gender, pathologic stage and resected tumor size. Patients with "gain of function" for m⁶A "reader" genes combined with copy number loss of "writers" or "erasers" had worse overall survival (OS) compared with other patterns. Moreover, copy number gain of m⁶A "reader" gene insulin growth factor 2 binding protein 2 (IGF2BP2) was an independent risk factor for OS (HR = 2.392, 95%CI: 1.392-4.112, p<0.001) and disease-free survival (DFS) (HR = 2.400, 95%CI: 1.236-4.659, p=0.010). Gene Set Enrichment Analysis (GSEA) indicated that IGF2BP2 was correlated with multiple biological processes associated with cancer, of which the most significant processes were relevant to cancer cell cycle, cell immortalization and tumor immunity. To sum up, a significant relationship was found between m⁶A genomic alterations and worse clinical outcomes. These innovative findings are expected to guide further research on the mechanism of m⁶A in PAAD.

Connecting METTL3 and intratumoural CD33+ MDSCs in predicting clinical outcome in cervical cancer

Background

Methyltransferase-like 3 (METTL3) is a member of the m⁶A methyltransferase family and acts as an oncogene in cancers. Recent studies suggest that host innate immunity is regulated by the enzymes controlling m⁶A epitranscriptomic changes. Here, we aim to explore the associations between the levels of METTL3 and CD33⁺ myeloid-derived suppressor cells (MDSCs) in tumour tissues and the survival of patients with cervical cancer (CC).

Methods

Specimens of paraffin embedded tumour from 197 CC patients were collected. The expression levels of METTL3 and CD33 were measured by immunohistochemical (IHC) staining. The clinical associations of the IHC variants were analysed by Pearson’s or Spearman’s chi-square tests. Overall survival (OS) and disease-free survival (DFS) were estimated by the Kaplan–Meier method and log-rank test. Hazard ratios (HRs) and independent significance were obtained via Cox proportional hazards models for multivariate analyses. METTL3 in CD33⁺ cells or CC-derived cells was knocked down by METTL3-specific siRNA, and MDSC induction in vitro was performed in a co-culture system in the presence of METTL3-siRNA and METTL3-knockdown-CC-derived cells compared with that of the corresponding controls.

Results

We found that tumour tissues displayed increased levels of METTL3 and CD33⁺ MDSCs compared with tumour-adjacent tissues from the same CC patients. Importantly, METTL3 expression was positively related to the density of CD33⁺ cells in tumour tissues (P = 0.011). We further found that the direct CD33⁺CD11b⁺HLA-DR⁻ MDSC induction and tumour-derived MDSC induction in vitro were decreased in the absence of METTL3. The level of METTL3 in tumour microenvironments was significantly related to advanced tumour stage. The levels of METTL3 and CD33⁺ MDSCs in tumour tissues were notably associated with reduced DFS or OS. Cox model analysis revealed that the level of METTL3 in tumour cells was an independent factor for patient survival, specifically for DFS (HR = 3.157, P = 0.022) and OS (HR = 3.271, P = 0.012), while the CD33⁺ MDSC number was an independent predictor for DFS (HR: 3.958, P = 0.031). Interestingly, in patients with advanced-disease stages (II–IV), METTL3 in tumour cells was an independent factor for DFS (HR = 6.725, P = 0.010) and OS (HR = 5.140, P = 0.021), while CD33⁺ MDSC density was an independent factor for OS (HR = 8.802, P = 0.037).

Conclusion

Our findings suggest that CD33⁺ MDSC expansion is linked to high levels of METTL3 and that METTL3 and CD33⁺ MDSCs are independent prognostic factors in CC.

N6 -methyladenosine (m6 A) RNA modification in human cancer

N⁶ -methyladenosine (m⁶ A) RNA modification, first discovered in 1974, is the most prevalent, abundant and penetrating messenger RNA (mRNA) modification in eukaryotes. This governs the fate of modified transcripts, regulates RNA metabolism and biological processes, and participates in pathogenesis of numerous human diseases, especially in cancer through the reciprocal regulation of m⁶ A methyltransferases ("writers") and demethylases ("erasers") and the binding proteins decoding m⁶ A methylation ("readers"). Accumulating evidence indicates a complicated regulation network of m⁶ A modification involving multiple m⁶ A-associated regulatory proteins whose biological functions have been further analysed. This review aimed to summarize the current knowledge on the potential significance and molecular mechanisms of m⁶ A RNA modification in the initiation and progression of cancer.

Functional Versatility of the CDK Inhibitor p57Kip2

The cyclin/CDK inhibitor p57^Kip2 belongs to the Cip/Kip family, with p21^Cip1 and p27^Kip1, and is the least studied member of the family. Unlike the other family members, p57^Kip2 has a unique role during embryogenesis and is the only CDK inhibitor required for embryonic development. p57^Kip2 is encoded by the imprinted gene CDKN1C, which is the gene most frequently silenced or mutated in the genetic disorder Beckwith-Wiedemann syndrome (BWS), characterized by multiple developmental anomalies. Although initially identified as a cell cycle inhibitor based on its homology to other Cip/Kip family proteins, multiple novel functions have been ascribed to p57^Kip2 in recent years that participate in the control of various cellular processes, including apoptosis, migration and transcription. Here, we will review our current knowledge on p57^Kip2 structure, regulation, and its diverse functions during development and homeostasis, as well as its potential implication in the development of various pathologies, including cancer.

Role of m6A RNA methylation in cardiovascular disease (Review)

N⁶-methyladenosine (m⁶A) is the most prevalent and abundant type of internal post-transcriptional RNA modification in eukaryotic cells. Multiple types of RNA, including mRNAs, rRNAs, tRNAs, long non-coding RNAs and microRNAs, are involved in m⁶A methylation. The biological function of m⁶A modification is dynamically and reversibly mediated by methyltransferases (writers), demethylases (erasers) and m⁶A binding proteins (readers). The methyltransferase complex is responsible for the catalyzation of m⁶A modification and is typically made up of methyltransferase-like (METTL)3, METTL14 and Wilms tumor 1-associated protein. Erasers remove methylation by fat mass and obesity-associated protein and ALKB homolog 5. Readers play a role through the recognition of m⁶A-modified targeted RNA. The YT521-B homology domain family, heterogeneous nuclear ribonucleoprotein and insulin-like growth factor 2 mRNA-binding protein serve as m⁶A readers. The m⁶A methylation on transcripts plays a pivotal role in the regulation of downstream molecular events and biological functions, such as RNA splicing, transport, stability and translatability at the post-transcriptional level. The dysregulation of m⁶A modification is associated with cancer, drug resistance, virus replication and the pluripotency of embryonic stem cells. Recently, a number of studies have identified aberrant m⁶A methylation in cardiovascular diseases (CVDs), including cardiac hypertrophy, heart failure, arterial aneurysm, vascular calcification and pulmonary hypertension. The aim of the present review article was to summarize the recent research progress on the role of m⁶A modification in CVD and give a brief perspective on its prospective applications in CVD.

N6-methyladenosine RNA modification in cancer therapeutic resistance: Current status and perspectives

Several strategies, including chemotherapy and radiotherapy, have improved therapeutic outcomes among cancer patients in clinical practice. However, due to their heterogeneity, cancer cells frequently display primary or acquired therapeutic resistance, thereby resulting in treatment failure. The mechanisms underlying cancer therapeutic resistance are complex and varied. Among them, N6-methyladenosine (m6A) RNA modification has gained increasing attention as a potential determinant of therapy resistance within various cancers. In this review, we primarily describe evidence for the effect of the m6A epitranscriptome on RNA homeostasis modulation, which has been shown to alter multiple cellular pathways in cancer research and treatment. Additionally, we discuss the profiles and biological implications of m6A RNA methylation, which is undergoing intensive investigation for its effect on the control of therapeutic resistance.

m6A-YTHDF1-mediated TRIM29 upregulation facilitates the stem cell-like phenotype of cisplatin-resistant ovarian cancer cells

Ovarian cancer is the deadliest gynaecologic malignancy, and the five-year survival rate of patients is less than 35% worldwide. Cancer stem cells (CSCs) are a population of cells with stem-like characteristics that are thought to cause chemoresistance and recurrence. TRIM29 is aberrantly expressed in various cancers and associated with cancer development and progression. Previous studies showed that the upregulation of TRIM29 expression in pancreatic cancer is related to stem-like characteristics. However, the role of TRIM29 in ovarian cancer is poorly understood. In this study, we found that TRIM29 expression was increased at the translational level in both the cisplatin-resistant ovarian cancer cells and clinical tissues. Increased TRIM29 expression was associated with a poor prognosis of patients with ovarian cancer. In addition, TRIM29 could enhance the CSC-like characteristics of the cisplatin-resistant ovarian cancer cells. Recruitment of YTHDF1 to m6A-modified TRIM29 was involved in promoting TRIM29 translation in the cisplatin-resistant ovarian cancer cells. Knockdown of YTHDF1 suppressed the CSC-like characteristics of the cisplatin-resistant ovarian cancer cells, which could be rescued by ectopic expression of TRIM29. This study suggests TRIM29 may act as an oncogene to promote the CSC-like features of cisplatin-resistant ovarian cancer in an m6A-YTHDF1-dependent manner. Due to the roles of TRIM29 and YTHDF1 in the promotion of CSC-like features, they may become potential therapeutic targets to combat the recurrence of ovarian cancer.

Molecular characterization, biological function, tumor microenvironment association and clinical significance of m6A regulators in lung adenocarcinoma

N6-methyladenosine (m6A) modification can regulate a variety of biological processes. However, the implications of m6A modification in lung adenocarcinoma (LUAD) remain largely unknown. Here, we systematically evaluated the m6A modification features in more than 2400 LUAD samples by analyzing the multi-omics features of 23 m6A regulators. We depicted the genetic variation features of m6A regulators, and found mutations of FTO and YTHDF3 were linked to worse overall survival. Many m6A regulators were aberrantly expressed in tumors, among which FTO, IGF2BP3, YTHDF1 and RBM15 showed consistent alteration features across 11 independent cohorts. Besides, the regulator-pathway interaction network demonstrated that m6A modification was associated with various biological pathways, including immune-related pathways. The correlation between m6A regulators and tumor microenvironment was also assessed. We found that LRPPRC was negatively correlated with most tumor-infiltrating immune cells. On the other hand, we established a scoring tool named m6Sig, which was positively correlated with PD-L1 expression and could reflect both the tumor microenvironment characterization and prognosis of LUAD patients. Comparison of CNV between high and low m6Sig groups revealed differences on chromosome 7. Application of m6Sig on an anti-PD-L1 immunotherapy cohort confirmed that the high m6Sig group demonstrated therapeutic advantages and clinical benefits. Our study indicated that m6A modification is involved in many aspects of LUAD and contributes to tumor microenvironment formation. A better understanding of m6A modification will provide more insights into the molecular mechanisms of LUAD and facilitate developing more effective personalized treatment strategies. A web application was built along with this study (http://www.bioinfo-zs.com/luadexpress/).

m6A Reader HNRNPA2B1 Promotes Esophageal Cancer Progression via Up-Regulation of ACLY and ACC1

N6-methyladenosine (m⁶A) modification is the most abundant modification on eukaryotic RNA. In recent years, lots of studies have reported that m⁶A modification and m⁶A RNA methylation regulators were involved in cancer progression. However, the m⁶A level and its regulators in esophageal cancer (ESCA) remain poorly understood. In this study, we analyzed the expression of m⁶A regulators using The Cancer Genome Atlas data and found 14 of 19 m⁶A regulators are significantly increased in ESCA samples. Then we performed a univariate Cox regression analysis and LASSO (least absolute shrinkage and selection operator) Cox regression model to investigate the prognostic role of m⁶A regulators in ESCA, and the results indicated that a two-gene prognostic signature including ALKBH5 and HNRNPA2B1 could predict overall survival of ESCA patients. Moreover, HNRNPA2B1 is higher expressed in high-risk scores subtype of ESCA, indicating that HNRNPA2B1 may be involved in ESCA development. Subsequently, we confirmed that the level of m⁶A and HNRNPA2B1 was significantly increased in ESCA. We also found that HNRNPA2B1 expression positively correlated with tumor diameter and lymphatic metastasis of ESCA. Moreover, functional study showed that knockdown of HNRNPA2B1 inhibited the proliferation, migration, and invasion of ESCA. Mechanistically, we found that knockdown of HNRNPA2B1 inhibited the expression of de novo fatty acid synthetic enzymes, ACLY and ACC1, and subsequently suppressed cellular lipid accumulation. In conclusion, our study provides critical clues to understand the role of m⁶A and its regulators in ESCA. Moreover, HNRNPA2B1 functions as an oncogenic factor in promoting ESCA progression via up-regulation of fatty acid synthesis enzymes ACLY and ACC1, and it may be a promising prognostic biomarker and therapeutic target for human ESCA.

N6-methyladenosine (m6A) RNA modification in cancer stem cells

Cancer stem cells (CSCs), a unique subset of undifferentiated cells with stem cell-like properties, have emerged as driving forces in mediating tumor growth, metastasis, and therapeutic resistance. Recent advances have highlighted that N6-methyladenosine (m6A) RNA modification plays an important role in cancer biology and CSCs. Dynamic m6A decoration has been demonstrated to be involved in CSC generation and maintenance, governing cancer progression and therapeutic resistance. In this review, we provide the first overview of the current knowledge of m6A modification implicated in CSCs and their impact on CSC properties, tumor progression, and responses to treatment. Finally, we also highlight the potential of m6A machinery as novel targets for cancer therapeutics. The involvement of m6A modification in CSCs provides a new direction for exploring cancer pathogenesis and inspires the development of effective strategies to fully eliminate both cancer cells and CSCs.

WTAP Gene Variants Confer Hepatoblastoma Susceptibility: A Seven-Center Case-Control Study

Hepatoblastoma is a rare disease, and its etiology remains to be revealed. Wilms tumor suppressor-1-associated protein (WTAP) plays a critical role in tumorigenesis. However, whether single nucleotide polymorphisms (SNPs) of the WTAP gene predispose to hepatoblastoma risk awaits to be investigated. With the use of the TaqMan assay, we evaluated the genotype frequencies of three WTAP SNPs (rs7766006 G > T, rs9457712 G > A, and rs1853259 A > G) in Chinese children with 313 hepatoblastoma patients and 1,446 controls. Among these three SNPs, only the rs7766006 T allele exhibited a significant association with hepatoblastoma risk (GT versus GG: adjusted odds ratio [OR] = 0.70, 95% confidence interval [CI] = 0.53-0.92, p = 0.009; GT/TT versus GG: adjusted OR = 0.73, 95% CI = 0.57-0.95, p = 0.017). Combined analysis indicated that subjects with two risk genotypes showed significantly higher hepatoblastoma risk, compared to individuals without a risk genotype (adjusted OR = 1.38, 95% CI = 1.02-1.88, p = 0.037). The stratified analysis revealed that the rs1853259 GG genotype, the rs7766006 GT/TT genotype, and two risk genotypes modified hepatoblastoma risk in certain subgroups. The significant results were validated by haplotype analyses and false-positive report probability analyses. Furthermore, the expression quantitative trait locus analysis indicated that rs7766006 T was associated with decreased expression of WTAP mRNA. Collectively, our results suggest that WTAP SNPs may be genetic modifiers for the development of hepatoblastoma.

Downregulation of m6 A reader YTHDC2 promotes tumor progression and predicts poor prognosis in non-small cell lung cancer

Background

m⁶A modification affects the pathological progress of many diseases by affecting RNA stability and translocation. YTHDC2, a m⁶A reader, is associated with multiple cancers; however, little is known of its role in non‐small cell lung cancer (NSCLC).

Methods

The GEPIA, Oncomine and GEO databases were analyzed to assess expression of YTHDC2 in NSCLC patients. Quantitative polymerase chain reaction, western blot and immunohistochemistry were used to detect YTHDC2 expression in different NSCLC cell lines (H1299, H460, H292 and A549) and patients. The effects of YTHDC2 on NSCLC cell lines (A549 and H1299) proliferation and migration were employed using CCK8 and transwell assays. The potential target RNAs of YTHDC2 were obtained from the POSTAR database. Functional enrichment analysis of YTHDC2 targeted RNAs was performed using the Metascape database.

Results

GEPIA, Oncomine and GEO databases showed low expression of YTHDC2 in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) patients. YTHDC2 expression was significantly decreased in different NSCLC cell lines and our clinical samples. Moreover, low expression of YTHDC2 was significantly associated with poor differentiation, lymph node metastasis, tumor size and stage. In addition, YTHDC2 could suppress the proliferation and migration ability of A549 and H1299 cell lines. Kaplan‐Meier Plotter database analysis revealed that patients with low level of YTHDC2 had a significantly poor prognosis. Finally, functional enrichment analysis of YTHDC2 targeted RNAs indicated several enriched pathways related to cancer.

Conclusions

These findings elucidate that YTHDC2 suppresses tumorigenesis in NSCLC, indicating that YTHDC2 may be a promising therapeutic target for NSCLC.

Key points

Significant findings of the study

This study demonstrated that the downregulation of YTHDC2 promotes tumor progression and predicts poor prognosis in non‐small cell lung cancer (NSCLC).

What this study adds

YTHDC2 might be a promising therapeutic target for non‐small cell lung cancer (NSCLC).

Comparative epigenetics in animal physiology: An emerging frontier

The unprecedented access to annotated genomes now facilitates the investigation of the molecular basis of epigenetic phenomena in phenotypically diverse animals. In this critical review, we describe the roles of molecular epigenetic mechanisms in regulating mitotically and meiotically stable spatiotemporal gene expression, phenomena that provide the molecular foundation for the intra-, inter-, and trans-generational emergence of physiological phenotypes. By focusing principally on emerging comparative epigenetic roles of DNA-level and transcriptome-level epigenetic mark dynamics in the emergence of phenotypes, we highlight the relationship between evolutionary conservation and innovation of specific epigenetic pathways, and their interplay as a priority for future study. This comparative approach is expected to significantly advance our understanding of epigenetic phenomena, as animals show a diverse array of strategies to epigenetically modify physiological responses. Additionally, we review recent technological advances in the field of molecular epigenetics (single-cell epigenomics and transcriptomics and editing of epigenetic marks) in order to (1) investigate environmental and endogenous factor dependent epigenetic mark dynamics in an integrative manner; (2) functionally test the contribution of specific epigenetic marks for animal phenotypes via genome and transcript-editing tools. Finally, we describe advantages and limitations of emerging animal models, which under the Krogh principle, may be particularly useful in the advancement of comparative epigenomics and its potential translational applications in animal science, ecotoxicology, ecophysiology, climate change science and wild-life conservation, as well as organismal health.

m6A modification in RNA: biogenesis, functions and roles in gliomas

The chemical modification of RNA is a newly discovered epigenetic regulation mechanism in cells and plays a crucial role in a variety of biological processes. N6-methyladenine (m6A) mRNA modification is the most abundant form of posttranscriptional RNA modification in eukaryotes. Through the development of m6A RNA sequencing, the relevant molecular mechanism of m6A modification has gradually been revealed. It has been found that the effect of m6A modification on RNA metabolism involves processing, nuclear export, translation and even decay. As the most common malignant tumour of the central nervous system, gliomas (especially glioblastoma) have a very poor prognosis, and treatment efficacy is not ideal even with the application of high-intensity treatment measures of surgery combined with chemoradiotherapy. Exploring the origin and development mechanisms of tumour cells from the perspective of tumour biogenesis has always been a hotspot in the field of glioma research. Emerging evidence suggests that m6A modification can play a key role in gliomas through a variety of mechanisms, providing more possibilities for early diagnosis and targeted therapy of gliomas. The aim of the present review is to focus on the research progress regarding the association between m6A modification and gliomas. And to provide a theoretical basis according to the currently available literature for further exploring this association. This review may provide new insights for the molecular mechanism, early diagnosis, histologic grading, targeted therapy and prognostic evaluation of gliomas.

m6A modification in RNA: biogenesis, functions and roles in gliomas

The chemical modification of RNA is a newly discovered epigenetic regulation mechanism in cells and plays a crucial role in a variety of biological processes. N6-methyladenine (m6A) mRNA modification is the most abundant form of posttranscriptional RNA modification in eukaryotes. Through the development of m6A RNA sequencing, the relevant molecular mechanism of m6A modification has gradually been revealed. It has been found that the effect of m6A modification on RNA metabolism involves processing, nuclear export, translation and even decay. As the most common malignant tumour of the central nervous system, gliomas (especially glioblastoma) have a very poor prognosis, and treatment efficacy is not ideal even with the application of high-intensity treatment measures of surgery combined with chemoradiotherapy. Exploring the origin and development mechanisms of tumour cells from the perspective of tumour biogenesis has always been a hotspot in the field of glioma research. Emerging evidence suggests that m6A modification can play a key role in gliomas through a variety of mechanisms, providing more possibilities for early diagnosis and targeted therapy of gliomas. The aim of the present review is to focus on the research progress regarding the association between m6A modification and gliomas. And to provide a theoretical basis according to the currently available literature for further exploring this association. This review may provide new insights for the molecular mechanism, early diagnosis, histologic grading, targeted therapy and prognostic evaluation of gliomas.

On the Way to Understanding the Interplay between the RNA Structure and Functions in Cells: A Genome-Wide Perspective

RNAs adopt specific structures in order to perform their biological activities. The structure of RNA is an important layer of gene expression regulation, and can impact a plethora of cellular processes, starting with transcription, RNA processing, and translation, and ending with RNA turnover. The development of high-throughput technologies has enabled a deeper insight into the sophisticated interplay between the structure of the cellular transcriptome and the living cells environment. In this review, we present the current view on the RNA structure in vivo resulting from the most recent transcriptome-wide studies in different organisms, including mammalians, yeast, plants, and bacteria. We focus on the relationship between the mRNA structure and translation, mRNA stability and degradation, protein binding, and RNA posttranscriptional modifications.

Multiomics profile and prognostic gene signature of m6A regulators in uterine corpus endometrial carcinoma

Uterine corpus endometrial carcinoma (UCEC) is the most common type of gynecologic malignancy worldwide. Despite advances in the treatments of UCEC, its incidence and mortality rates are still increasing. N6-methyladenosine (m6A) is the most common form of RNA modification and has attracted increasing interest in cancer pathogenesis and progression. Thus, we aimed to identify the landscape of m6A regulators and build a prognostic gene signature in UCEC. In this study, we first analyzed copy number variations (CNVs), single nucleotide variations (SNVs) and gene expression profiles as well as matched clinical information of UCEC patients from The Cancer Genome Atlas (TCGA) database. Next, we determined that CNVs in m6A regulatory genes had a significant negative impact on patient survival. The mRNA expression levels of a total of 16 m6A regulators were significantly correlated with different CNV patterns. Using univariate Cox regression analysis, IGF2BP1, KIAA1429, IGF2BP3, YTHDF3, and IGF2BP2 were found to be closely associated with UCEC patient survival outcomes. Based on the least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression models, we built a 3-gene (IGF2BP3, KIAA1429 and IGF2BP1) signature of m6A regulators with prognostic value in UCEC that could effectively predict patient prognosis (log-rank test p-value < 0.0001). In addition, risk scores were significantly different between patients stratified by tumor stage, SNV, and CNV. Multivariate Cox regression analysis suggested that risk score might be an independent prognostic indicator for the overall survival of patients with UCEC (p-value < 0.05). Gene enrichment analysis indicated that high IGF2BP1 gene expression is associated with cytoplasmic stress granules. KIAA1429 gene expression is associated with cellular nucleic acid metabolism. The expression of the IGF2BP3 gene is associated with RNA binding processes. In conclusion, we determined that genetic alterations in m6A regulatory genes could be effective and reliable biomarkers for UCEC prognosis prediction.

The Potential Roles of RNA N6-Methyladenosine in Urological Tumors

N6-methyladenosine (m⁶A) is regarded as the most abundant, prevalent and conserved internal mRNA modification in mammalian cells. M⁶A can be catalyzed by m⁶A methyltransferases METTL3, METTL14 and WTAP (writers), reverted by demethylases ALKBH5 and FTO (erasers), and recognized by m⁶A -binding proteins such as YTHDF1/2/3, IGF2BP1/2/3 and HNRNPA2B1 (readers). Emerging evidence suggests that m⁶A modification is significant for regulating many biological and cellular processes and participates in the pathological development of various diseases, including tumors. This article reviews recent studies on the biological function of m⁶A modification and the methylation modification of m⁶A in urological tumors.

The contribution of WTAP gene variants to Wilms tumor susceptibility

Wilms tumor is the most frequently occurring pediatric renal malignancy. Wilms tumor suppressor-1-associated protein (WTAP) is a vital component of N6-methyltransferase complex involved in tumorigenesis. However, the roles of WTAP gene single nucleotide polymorphisms (SNPs) in Wilms tumor risk have not been clarified to date. We successfully genotyped three WTAP gene SNPs using TaqMan assay in 405 Wilms tumor patients and 1197 cancer-free controls of Chinese children. Odds ratios (ORs) and 95% confidence intervals (CIs) were applied to determine the effects of WTAP gene SNPs on Wilms tumor risk. Carriers of the rs1853259 G variant are less susceptible to developing Wilms tumor, with an adjusted OR of 0.78 (AG vs. AA: 95% CI = 0.61-0.995, P = 0.046). Single locus analysis of rs9457712 G > A and rs7766006 G > T, as well as the combined analysis of risk genotypes, failed to unveil an association with Wilms tumor risk, respectively. Stratified analysis of the three SNPs and their combined risk effects showed more significant relationships with Wilms tumor risk under certain subgroups. In all, we found weak evidence of the association between WTAP gene SNPs and the risk of Wilms tumor. Further replication studies with greater sample size and different ethnicities are necessary to verify our findings.

Multifaceted Functions and Novel Insight Into the Regulatory Role of RNA N6-Methyladenosine Modification in Musculoskeletal Disorders

RNA modifications have emerged as key regulators of transcript expression in diverse physiological and pathological processes. As one of the most prevalent types of RNA modifications, N⁶-methyladenosine (m⁶A) has become the highlight in modulation of various diseases through interfering RNA splicing, translation, nuclear export, and decay. In many cases, the detailed functions of m⁶A in cellular processes and diseases remain unclear. Notably, recent studies have determined the relationship between m⁶A modification and musculoskeletal disorders containing osteosarcoma, osteoarthritis, rheumatoid arthritis, osteoporosis, etc. Herein, this review comprehensively summarizes the recent advances of m⁶A modification in pathogenesis and progression of musculoskeletal diseases. Specifically, the underlying molecular mechanisms, detection technologies, regulatory functions, clinical implications, and future perspectives of m⁶A in musculoskeletal disorders are discussed, with the aim to provide a novel insight into their association.