Post-transcriptional RNA Modification Research Articles

YTHDF1 Aggravates the Progression of Cervical Cancer Through m6A-Mediated Up-Regulation of RANBP2.

N6-methyladenosine (m6A) is the most common post-transcriptional modification of RNA in eukaryotes, which has been demonstrated to play important roles in various cancers. YTHDF1 acts as a crucial m6A “reader” and regulates the fate of m6A modified mRNA. However, its role in cervical cancer remains unknown. In this study, we showed that YTHDF1 was highly expressed in cervical cancer, and was closely associated with the poor prognosis of cervical cancer patients. YTHDF1 knockdown suppressed the growth, migration and invasion, and induced apoptosis of cervical cancer cells. Moreover, YTHDF1 knockdown inhibited tumorigenesis of cervical cancer cells in vivo. Through combined on-line data analysis of RIP-seq, meRIP-seq and Ribo-seq upon YTHDF1 knockdown, RANBP2 was identified as the key target of YTHDF1 in cervical cancer cells. YTHDF1 regulated RANBP2 translation in an m6A-dependent manner without effect on its mRNA expression. RANBP2 potentiated the growth, migration and invasion of cervical cancer cells. Our study demonstrated the oncogenic role of YTHDF1 in cervical cancer by regulating RANBP2 expression and YTHDF1 represents a potential target for cervical cancer therapy.

DNN-m6A: A Cross-Species Method for Identifying RNA N6-Methyladenosine Sites Based on Deep Neural Network with Multi-Information Fusion.

As a prevalent existing post-transcriptional modification of RNA, N6-methyladenosine (m6A) plays a crucial role in various biological processes. To better radically reveal its regulatory mechanism and provide new insights for drug design, the accurate identification of m6A sites in genome-wide is vital. As the traditional experimental methods are time-consuming and cost-prohibitive, it is necessary to design a more efficient computational method to detect the m6A sites. In this study, we propose a novel cross-species computational method DNN-m6A based on the deep neural network (DNN) to identify m6A sites in multiple tissues of human, mouse and rat. Firstly, binary encoding (BE), tri-nucleotide composition (TNC), enhanced nucleic acid composition (ENAC), K-spaced nucleotide pair frequencies (KSNPFs), nucleotide chemical property (NCP), pseudo dinucleotide composition (PseDNC), position-specific nucleotide propensity (PSNP) and position-specific dinucleotide propensity (PSDP) are employed to extract RNA sequence features which are subsequently fused to construct the initial feature vector set. Secondly, we use elastic net to eliminate redundant features while building the optimal feature subset. Finally, the hyper-parameters of DNN are tuned with Bayesian hyper-parameter optimization based on the selected feature subset. The five-fold cross-validation test on training datasets show that the proposed DNN-m6A method outperformed the state-of-the-art method for predicting m6A sites, with an accuracy (ACC) of 73.58–83.38% and an area under the curve (AUC) of 81.39–91.04%. Furthermore, the independent datasets achieved an ACC of 72.95–83.04% and an AUC of 80.79–91.09%, which shows an excellent generalization ability of our proposed method.

Characterizing RNA Pseudouridylation by Convolutional Neural Networks

Pseudouridine (Ψ) is the most prevalent post-transcriptional RNA modification and is widespread in small cellular RNAs and mRNAs. However, the functions, mechanisms, and precise distribution of Ψs (especially in mRNAs) still remain largely unclear. The landscape of Ψs across the transcriptome has not yet been fully delineated. Here, we present a highly effective model based on a convolutional neural network (CNN), called PseudoUridyLation Site Estimator (PULSE), to analyze large-scale profiling data of Ψ sites and characterize the contextual sequence features of pseudouridylation. PULSE, consisting of two alternatively-stacked convolution and pooling layers followed by a fully-connected neural network, can automatically learn the hidden patterns of pseudouridylation from the local sequence information. Extensive validation tests demonstrated that PULSE can outperform other state-of-the-art prediction methods and achieve high prediction accuracy, thus enabling us to further characterize the transcriptome-wide landscape of Ψ sites. We further showed that the prediction results derived from PULSE can provide novel insights into understanding the functional roles of pseudouridylation, such as the regulations of RNA secondary structure, codon usage, translation, and RNA stability, and the connection to single nucleotide variants. The source code and final model for PULSE are available at https://github.com/mlcb-thu/PULSE.

N6-methyladenosine modification of MALAT1 promotes metastasis via reshaping nuclear speckles.

N6-methyladenosine (m6A), one of the most prevalent RNA post-transcriptional modifications, is involved in numerous biological processes. In previous studies, the functions of m6A were typically identified by perturbing the activity of the methyltransferase complex. Here, we dissect the contribution of m6A to an individual-long noncoding RNA-metastasis-associated lung adenocarcinoma transcript 1 (MALAT1). The mutant MALAT1 lacking m6A-motifs significantly suppressed the metastatic potential of cancer cells both invitro and invivo in mouse. Super-resolution imaging showed that the concatenated m6A residues on MALAT1 acted as a scaffold for recruiting YTH-domain-containing protein 1 (YTHDC1) to nuclear speckles. We further reveal that the recognition of MALAT1-m6A by YTHDC1 played a critical role in maintaining the composition and genomic binding sites of nuclear speckles, which regulate the expression of several key oncogenes. Furthermore, artificially tethering YTHDC1 onto m6A-deficient MALAT1 largely rescues the metastatic potential of cancer cells.

Alteration of mRNA 5-Methylcytosine Modification in Neurons After OGD/R and Potential Roles in Cell Stress Response and Apoptosis.

Epigenetic modifications play an important role in central nervous system disorders. As a widespread posttranscriptional RNA modification, the role of the m5C modification in cerebral ischemia-reperfusion injury (IRI) remains poorly defined. Here, we successfully constructed a neuronal oxygen-glucose deprivation/reoxygenation (OGD/R) model and obtained an overview of the transcriptome-wide m5C profiles using RNA-BS-seq. We discovered that the distribution of neuronal m5C modifications was highly conserved, significantly enriched in CG-rich regions and concentrated in the mRNA translation initiation regions. After OGD/R, modification level of m5C increased, whereas the number of methylated mRNA genes decreased. The amount of overlap of m5C sites with the binding sites of most RNA-binding proteins increased significantly, except for that of the RBM3-binding protein. Moreover, hypermethylated genes in neurons were significantly enriched in pathological processes, and the hub hypermethylated genes RPL8 and RPS9 identified by the protein-protein interaction network were significantly related to cerebral injury. Furthermore, the upregulated transcripts with hypermethylated modification were enriched in the processes involved in response to stress and regulation of apoptosis, and these processes were not identified in hypomethylated transcripts. In final, we verified that OGD/R induced neuronal apoptosis in vitro using TUNEL and western blot assays. Our study identified novel m5C mRNAs associated with ischemia-reperfusion in neurons, providing valuable perspectives for future studies on the role of the RNA methylation in cerebral IRI.

Loss of Cnot6l Impairs Inosine RNA Modifications in Mouse Oocytes.

Mammalian oocytes must degrade maternal transcripts through a process called translational mRNA decay, in which maternal mRNA undergoes translational activation, followed by deadenylation and mRNA decay. Once a transcript is translationally activated, it becomes deadenylated by the CCR4-NOT complex. Knockout of CCR4-NOT Transcription Complex Subunit 6 Like (Cnot6l), a deadenylase within the CCR4-NOT complex, results in mRNA decay defects during metaphase I (MI) entry. Knockout of B-cell translocation gene-4 (Btg4), an adaptor protein of the CCR4-NOT complex, results in mRNA decay defects following fertilization. Therefore, mechanisms controlling mRNA turnover have significant impacts on oocyte competence and early embryonic development. Post-transcriptional inosine RNA modifications can impact mRNA stability, possibly through a translation mechanism. Here, we assessed inosine RNA modifications in oocytes, eggs, and embryos from Cnot6l-/- and Btg4-/- mice, which display stabilization of mRNA and over-translation of the stabilized transcripts. If inosine modifications have a role in modulating RNA stability, we hypothesize that in these mutant backgrounds, we would observe changes or a disruption in inosine mRNA modifications. To test this, we used a computational approach to identify inosine RNA modifications in total and polysomal RNA-seq data during meiotic maturation (GV, MI, and MII stages). We observed pronounced depletion of inosine mRNA modifications in samples from Cnot6l-/-, but not in Btg4-/- mice. Additionally, analysis of ribosome-associated RNA revealed clearance of inosine modified mRNA. These observations suggest a novel mechanism of mRNA clearance during oocyte maturation, in which inosine-containing transcripts decay in an independent, but parallel mechanism to CCR4-NOT deadenylation.

Fumonisin B1 alters global m6A RNA methylation and epigenetically regulates Keap1-Nrf2 signaling in human hepatoma (HepG2) cells.

FB1 is a common contaminant of cereal grains that affects human and animal health. It has become increasingly evident that epigenetic changes are implicated in FB1 toxicity. N6-methyladenosine (m6A), the most abundant post-transcriptional RNA modification, is influenced by fluctuations in redox status. Since oxidative stress is a characteristic of FB1 exposure, we determined if there is cross-talk between oxidative stress and m6A in FB1-exposed HepG2 cells. Briefly, HepG2 cells were treated with FB1 (0, 5, 50, 100, 200µM; 24h) and ROS, LDH and m6A levels were quantified. qPCR was used to determine the expression of m6A modulators, Nrf2, Keap1 and miR-27b, while western blotting was used to quantify Keap1 and Nrf2 protein expression. Methylation status of Keap1 and Nrf2 promoters was assessed and RNA immunoprecipitation quantified m6A-Keap1 and m6A-Nrf2 levels. FB1 induced accumulation of intracellular ROS (p ≤ 0.001) and LDH leakage (p ≤ 0.001). Elevated m6A levels (p ≤ 0.05) were accompanied by an increase in m6A "writers" [METLL3 (p ≤ 0.01) and METLL14 (p ≤ 0.01)], and "readers" [YTHDF1 (p ≤ 0.01), YTHDF2 (p ≤ 0.01), YTHDF3 (p ≤ 0.001) and YTHDC2 (p ≤ 0.01)] and a decrease in m6A "erasers" [ALKBH5 (p ≤ 0.001) and FTO (p ≤ 0.001)]. Hypermethylation and hypomethylation occurred at Keap1 (p ≤ 0.001) and Nrf2 (p ≤ 0.001) promoters, respectively. MiR-27b was reduced (p ≤ 0.001); however, m6A-Keap1 (p ≤ 0.05) and m6A-Nrf2 (p ≤ 0.01) levels were upregulated. This resulted in the ultimate decrease in Keap1 (p ≤ 0.001) and increase in Nrf2 (p ≤ 0.001) expression. Our findings reveal that m6A RNA methylation can be modified by exposure to FB1, and a cross-talk between m6A and redox regulators does occur.

N6-methyladenosine (m6A) is an endogenous A3 adenosine receptor ligand

About 150 post-transcriptional RNA modifications have been identified in all kingdoms of life. During RNA catabolism, most modified nucleosides are resistant to degradation and are released into the extracellular space. In this study, we explored the physiological role of these extracellular modified nucleosides and found that N6-methyladenosine (m6A), widely recognized as an epigenetic mark in RNA, acts as a ligand for the human adenosine A3 receptor, for which it has greater affinity than unmodified adenosine. We used structural modeling to define the amino acids required for specific binding of m6A to the human A3 receptor. We alsodemonstrated that m6A was dynamically released in response to cytotoxic stimuli and facilitated type I allergy invivo. Our findings implicate m6A as a signaling molecule capable of activating G protein-coupled receptors (GPCRs) and triggering pathophysiological responses, a previously unreported property of RNA modifications.

RNA Post-Transcriptional Modification Mapping Data Analysis Using RNA Framework.

RNA post-transcriptional modifications (PTMs) are progressively gaining relevance in the study of coding-independent functions of RNA. RNA PTMs act as dynamic regulators of several aspects of the RNA physiology, from translation to half-life. Rising interest is supported by the advance of high-throughput techniques enabling the detection of these modifications on a transcriptome-wide scale. To this end, here we illustrate the usefulness of RNA Framework, a comprehensive toolkit for the analysis of RNA PTM mapping experiments, by reanalyzing two published transcriptome-scale datasets of N1-methyladenosine (m1A) and pseudouridine (Ψ) mapping, based on two different experimental strategies.

Assaying the Molecular Determinants and Kinetics of RNA Pseudouridylation by H/ACA snoRNPs and Stand-Alone Pseudouridine Synthases.

Posttranscriptional modifications of RNA play an important role in promoting the maturation and functional diversity of many RNA species. Accordingly, understanding the enzymes and mechanisms that underlie RNA modifications is an important aspect in advancing our knowledge of the continually expanding RNA modification field. However, of the more than 160 currently identified RNA modifications, a large portion remains without quantitative detection assays for their biochemical characterization. Here, we describe the tritium release assay as a convenient tool allowing for the quantitative assessment of in vitro RNA pseudouridylation by stand-alone or box H/ACA RNA-guided pseudouridine synthases. This assay enables quantification of RNA pseudouridylation over a time course to effectively compare pseudouridylation activity between different substrates and/or different recombinant enzymes as well as to determine kinetic parameters. With the help of a quench-flow apparatus, the tritium release assay can be adapted for rapid kinetic measurements under single-turnover conditions to dissect reaction mechanisms. As examples, we show the formation of pseudouridines by a reconstituted Saccharomyces cerevisiae H/ACA small ribonucleoprotein (snoRNP) and an Escherichia coli stand-alone pseudouridine synthase.

Metabolomic Profiles and Heart Failure Risk in Black Adults: Insights From the Jackson Heart Study.

Heart failure (HF) is a heterogeneous disease characterized by significant metabolic disturbances; however, the breadth of metabolic dysfunction before the onset of overt disease is not well understood. The purpose of this study was to determine the association of circulating metabolites with incident HF to uncover novel metabolic pathways to disease. We performed targeted plasma metabolomic profiling in a deeply phenotyped group of Black adults from the JHS (Jackson Heart Study; n=2199). We related metabolites associated with incident HF to established etiological mechanisms, including increased left ventricular mass index and incident coronary heart disease. Furthermore, we evaluated differential associations of metabolites with HF with preserved ejection fraction versus HF with reduced ejection fraction. Metabolites associated with incident HF included products of posttranscriptional modifications of RNA, as well as polyamine and nitric oxide metabolism. A subset of metabolite-HF associations was independent of well-established HF pathways such as increased left ventricular mass index and incident coronary heart disease and included homoarginine (per 1 SD increase in metabolite level, hazard ratio, 0.77; P=1.2×10-3), diacetylspermine (hazard ratio, 1.34; P=3.4×10-3), and uridine (hazard ratio, 0.79; P, 3×10-4). Furthermore, metabolites involved in pyrimidine metabolism (orotic acid) and collagen turnover (N-methylproline) among others were part of a distinct metabolic signature that differentiated individuals with HF with preserved ejection fraction versus HF with reduced ejection fraction. The integration of clinical phenotyping with plasma metabolomic profiling uncovered novel metabolic processes in nontraditional disease pathways underlying the heterogeneity of HF development in Black adults.

Translational adaptation to heat stress is mediated by RNA 5-methylcytosine in Caenorhabditis elegans.

Methylation of carbon‐5 of cytosines (m5C) is a post‐transcriptional nucleotide modification of RNA found in all kingdoms of life. While individual m5C‐methyltransferases have been studied, the impact of the global cytosine‐5 methylome on development, homeostasis and stress remains unknown. Here, using Caenorhabditis elegans, we generated the first organism devoid of m5C in RNA, demonstrating that this modification is non‐essential. Using this genetic tool, we determine the localisation and enzymatic specificity of m5C sites in the RNome in vivo. We find that NSUN‐4 acts as a dual rRNA and tRNA methyltransferase in C. elegans mitochondria. In agreement with leucine and proline being the most frequently methylated tRNA isoacceptors, loss of m5C impacts the decoding of some triplets of these two amino acids, leading to reduced translation efficiency. Upon heat stress, m5C loss leads to ribosome stalling at UUG triplets, the only codon translated by an m5C34‐modified tRNA. This leads to reduced translation efficiency of UUG‐rich transcripts and impaired fertility, suggesting a role of m5C tRNA wobble methylation in the adaptation to higher temperatures.

Detection of 5-formylcytosine in Mitochondrial Transcriptome.

Posttranscriptional RNA modifications have recently emerged as essential posttranscriptional regulators of gene expression. Here we present two methods for single nucleotide resolution detection of 5-formylcytosine (f5C) in RNA. The first relies on chemical protection of f5C against bisulfite treatment, the second method is based on chemical reduction of f5C to hm5C. In combination with regular bisulfite treatment of RNA, the methods allow for precise mapping of f5C. The protocol is used for f5C detection in mtDNA-encoded RNA, however, it can be straightforwardly applied for transcriptome-wide analyses.

Histone mRNA is subject to 3' uridylation and re-adenylation in Aspergillus nidulans.

The role of post-transcriptional RNA modification is of growing interest. One example is the addition of non-templated uridine residues to the 3' end of transcripts. In mammalian systems, uridylation is integral to cell cycle control of histone mRNA levels. This regulatory mechanism is dependent on the nonsense-mediated decay (NMD) component, Upf1, which promotes histone mRNA uridylation and degradation in response to the arrest of DNA synthesis. We have identified a similar system in Aspergillus nidulans, where Upf1 is required for the regulation of histone mRNA levels. However, other NMD components are also implicated, distinguishing it from the mammalian system. As in human cells, 3' uridylation of histone mRNA is induced upon replication arrest. Disruption of this 3' tagging has a significant but limited effect on histone transcript regulation, consistent with multiple mechanisms acting to regulate mRNA levels. Interestingly, 3' end degraded transcripts are also subject to re-adenylation. Both mRNA pyrimidine tagging and re-adenylation are dependent on the same terminal-nucleotidyltransferases, CutA, and CutB, and we show this is consistent with the in vitro activities of both enzymes. Based on these data we argue that mRNA 3' tagging has diverse and distinct roles associated with transcript degradation, functionality and regulation.

Emerging role of N4-acetylcytidine modification of RNA in gene regulation and cellular functions.

Post-transcriptional chemical modification of RNA is rapidly emerging as a key player in regulating gene expression and has propelled the development of 'epitranscriptomics' or 'RNA epigenetics' as a frontier area of research. Several RNA modifications are known to decorate RNAs and impact its structure and function. One such recently discovered modification is acetylation of RNA i.e. N4-acetylcytidine (ac4C) chemical modification. N4-acetylcytidine is an ancient and evolutionarily conserved modification, which maps to a wide spectrum of RNAs from archaea bacteria to humans. This modification results in a variety of functional outcomes which impact normal development and disease. In this review, we summarize the recent progress, emerging methods, biological implications and the future challenges for ac4C modification.

Detecting N6-methyladenosine sites from RNA transcriptomes using random forest

N6-methyladenosine (m6A) modifications are one the most frequently occurring RNA post transcriptional modifications. These modifications perform vital roles in different biological processes, including, localization and translation of proteins, X chromosome inactivation, cell stability, microRNA regulation, and reprogramming etc. Any abnormal change in m6A sites may lead to several abnormalities, including, cancer, brain-related disorders and many other life threatening diseases. Precise detection of m6A modifications is crucial for the diagnosis and treatment of these diseases. Existing methods suffer from the problem of inefficient detection of m6A sites, especially in yeast transcriptomes (due to varied structure) and inability of the computational techniques to capture the encoded information surrounding the m6A sites. In this work, we propose a novel method (called m6A-pred predictor) that utilizes a fusion of characteristics including, statistical, and chemical properties of the nucleotides, to precisely predict the presence of m6A sites in RNA sequences. The fusion of multiple types of features results in a high dimensional vector which is further optimized using an evolutionary algorithm. Finally, the random forest classifier is used to detect m6A sites by using the most discriminative features. The results, benchmarked on yeast transcriptomes, indicate that m6A-pred predictor outperforms all the previously reported predictors, notably, with an accuracy value of 78.58%, specificity value of 79.65% and Matthews correlation coefficient of 0.5717.

M5CPred-SVM: a novel method for predicting m5C sites of RNA

BackgroundAs one of the most common post-transcriptional modifications (PTCM) in RNA, 5-cytosine-methylation plays important roles in many biological functions such as RNA metabolism and cell fate decision. Through accurate identification of 5-methylcytosine (m5C) sites on RNA, researchers can better understand the exact role of 5-cytosine-methylation in these biological functions. In recent years, computational methods of predicting m5C sites have attracted lots of interests because of its efficiency and low-cost. However, both the accuracy and efficiency of these methods are not satisfactory yet and need further improvement.ResultsIn this work, we have developed a new computational method, m5CPred-SVM, to identify m5C sites in three species, H. sapiens, M. musculus and A. thaliana. To build this model, we first collected benchmark datasets following three recently published methods. Then, six types of sequence-based features were generated based on RNA segments and the sequential forward feature selection strategy was used to obtain the optimal feature subset. After that, the performance of models based on different learning algorithms were compared, and the model based on the support vector machine provided the highest prediction accuracy. Finally, our proposed method, m5CPred-SVM was compared with several existing methods, and the result showed that m5CPred-SVM offered substantially higher prediction accuracy than previously published methods. It is expected that our method, m5CPred-SVM, can become a useful tool for accurate identification of m5C sites.ConclusionIn this study, by introducing position-specific propensity related features, we built a new model, m5CPred-SVM, to predict RNA m5C sites of three different species. The result shows that our model outperformed the existing state-of-art models. Our model is available for users through a web server at https://zhulab.ahu.edu.cn/m5CPred-SVM.

Next-Generation RNA-Sequencing of Serum Small Extracellular Vesicles Discovers Potential Diagnostic Biomarkers for Dementia With Lewy Bodies

ObjectiveThere is an urgent clinical need for identifying blood-based diagnostic biomarkers for Dementia with Lewy Bodies (DLB). Transcriptomic studies have reported unique RNA changes in postmortem DLB brains. Small extracellular vesicles (SEV) that transport RNA between brain and peripheral circulation enable identifying molecular changes in living human brain. Hence, we aimed to identify differentially expressed RNA in serum SEVs from people with DLB. MethodsWe investigated serum SEV total RNA profiles in people with DLB (n = 10) and age and gender matched comparisons (n = 10) using next-generation RNA-sequencing. SEVs were separated by ultracentrifugation with density gradient and were characterized by nanoparticle analysis and western blotting. We verified the differential expression levels of identified differentially expressed genes (DEG) using high-throughput qPCR. Functional implications of identified DEG were evaluated using Ingenuity pathway analyses. ResultsWe identified 846 nominally significant DEG including 30 miRNAs in DLB serum SEVs. We identified significant downregulation of proinflammatory genes, IL1B, CXCL8, and IKBKB. Previously reported postmortem DLB brain DEGs were significantly enriched (χ2=4.99; df=1; p = 0.03) among the identified DEGs, and the differential expression of 40 postmortem DLB brain DEGs could be detected in serum SEVs of people living with DLB. Functional pathway and network analyses highlighted the importance of immunosenescence, ubiquitin proteasome system (UPS) dysfunction, DNA repair, and RNA post-transcriptional modification deficits in DLB pathology. ConclusionIdentified DEGs, especially reduced expression levels of inflammation, and UPS-associated RNA, may aid diagnosing DLB, and their biomarker potential warrants further investigation in larger clinical cohorts. Our findings corroborate the absence of chronic neuroinflammation in DLB.

2'-O-Methylation can increase the abundance and lifetime of alternative RNA conformational states.

2′-O-Methyl (Nm) is a highly abundant post-transcriptional RNA modification that plays important biological roles through mechanisms that are not entirely understood. There is evidence that Nm can alter the biological activities of RNAs by biasing the ribose sugar pucker equilibrium toward the C3′-endo conformation formed in canonical duplexes. However, little is known about how Nm might more broadly alter the dynamic ensembles of flexible RNAs containing bulges and internal loops. Here, using NMR and the HIV-1 transactivation response (TAR) element as a model system, we show that Nm preferentially stabilizes alternative secondary structures in which the Nm-modified nucleotides are paired, increasing both the abundance and lifetime of low-populated short-lived excited states by up to 10-fold. The extent of stabilization increased with number of Nm modifications and was also dependent on Mg2+. Through phi-value analysis, the Nm modification also provided rare insights into the structure of the transition state for conformational exchange. Our results suggest that Nm could alter the biological activities of Nm-modified RNAs by modulating their secondary structural ensembles as well as establish the utility of Nm as a tool for the discovery and characterization of RNA excited state conformations.

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(m6A) 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 m6A 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 m6A regulators and clinicopathological characteristics was examined using chi-square test. Results showed a high frequency of copy number alterations (CNAs) of m6A regulatory genes in PAAD patients, but somatic mutations were rarely happened. CNAs and mutations of m6A regulatory genes was associated with patient’s gender, pathologic stage and resected tumor size. Patients with “gain of function” for m6A “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 m6A “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 m6A genomic alterations and worse clinical outcomes. These innovative findings are expected to guide further research on the mechanism of m6A in PAAD.