Dataset Information

Genome-wide identification and characterization of circular RNA m6A modification in pancreatic cancer

ABSTRACT:

Background

N6-methyladenosine (m6A) is the most abundant modification of RNA in eukaryotic cells and play critical roles in cancer. While most related studies focus on m6A modifications in linear RNA, transcriptome-wide profiling and exploration of m6A modification in circular RNAs in cancer is still lacking.

Methods

For the detection of m6A modification in circRNAs, we developed a new bioinformatics tools called Circm6A and applied it to the m6A-seq data of 77 tissue samples from 58 individuals with pancreatic ductal adenocarcinoma (PDAC).

Results

Circm6A performs better than the existing circRNA identification tools, which achieved highest F1 score among these tools in the detection of circRNAs with m6A modifications. By using Circm6A, we identified a total of 8807 m6A-circRNAs from our m6A-seq data. The m6A-circRNAs tend to be hypermethylated in PDAC tumor tissues compared with normal tissues. The hypermethylated m6A-circRNAs were associated with a significant gain of circRNA-mRNA coexpression network, leading to the dysregulation of many important cancer-related pathways. Moreover, we found the cues that hypermethylated m6A-circRNAs may promote the circularization and translation of circRNAs.

Conclusions

These comprehensive findings further bridged the knowledge gaps between m6A modification and circRNAs fields by depicting the m6A-circRNAs genomic landscape of PDAC patients and revealed the emerging roles played by m6A-circRNAs in pancreatic cancer. Circm6A is available at https://github.com/canceromics/circm6a.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-01002-w.

SUBMITTER: Ye Y

PROVIDER: S-EPMC8605608 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background:N6-methyladenosine (m6A) modification holds an important position in tumorigenesis and metastasis because it can change gene expression and even function in multiple levels including RNA splicing, stability, translocation and translation. In present study, we aim to conducted comprehensive investigation on m6A RNA methylation regulators and m6A-related genes in pancreatic cancer and their association with survival time. Methods:Based on Univariate Cox regression analysis, protein-protein interaction analysis, LASSO Cox regression, a risk prognostic model, STRING, Spearman and consensus clustering analysis, data from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) database was used to analyze 15 m6A RNA methylation regulators that were widely reported and 1,393 m6A-related genes in m6Avar. Results:We found that 283 candidate m6A RNA methylation-related genes and 4 m6A RNA methylation regulatory factors, including RNA binding motif protein 15 (RBM15), methyltransferase like 14 (METTL14), fat mass and obesity-associated protein (FTO), and ?-ketoglutarate-dependent dioxygenase AlkB homolog 5 (ALKBH5), differed significantly among different stages of the American Joint Committee on Cancer (AJCC) staging system. Protein-protein interaction analysis indicated epidermal growth factor receptor (EGFR), plectin-1 (PLEC), BLM RecQ like helicase (BLM), and polo like kinase 1 (PLK1) were closely related to other genes and could be considered as hub genes in the network. The results of LASSO Cox regression and the risk prognostic model indicated that AJCC stage, stage T and N, KRAS mutation status and x8q23.3 CNV fragment mutation differed significantly between the high-risk and the low-risk subgroups. The AUCs of 1 to 5 years after surgery were all more than 0.7 and increased year by year. Finally, we found KRAS mutation status and AJCC stage differed significantly among these groups after TCGA samples divided into subgroups with k=7. Moreover, we identified four m6A RNA methylation related genes expressed significantly differently among these seven subgroups, including collagen type VII alpha 1 chain (COL7A1), branched chain amino acid transaminase 1 (BCAT1), zinc finger protein 596 (ZNF596), and PLK1. Conclusions:Our study systematically analyzed the m6A RNA methylation related genes, including expression, protein-protein interaction, potential function, and prognostic value and provides important clues to further research on the function of RNA m6A methylation and its related genes in pancreatic cancer.

Project description:Eukaryotic messenger RNAs (mRNAs) are often modified with methyl groups at the N6 position of adenosine (m6A), and these changes are interpreted by YTH domain-containing proteins to regulate the metabolism of m6A-modified mRNAs. Although alfalfa (Medicago sativa) is an established model organism for forage development, the understanding of YTH proteins in alfalfa is still limited. In the present investigation, 53 putative YTH genes, each encoding a YT521 domain-containing protein, were identified within the alfalfa genome. These genes were categorized into two subfamilies: YTHDF (49 members) and YTHDC (four members). Each subfamily demonstrates analogous motif distributions and domain architectures. Specifically, proteins encoded by MsYTHDF genes incorporate a single domain structure, while those corresponding to MsYTH5, 8, 12, 16 who are identified as members of the MsYTHDC subfamily, exhibit CCCH-type zinc finger repeats at their N-termini. It is also observed that the predicted aromatic cage pocket that binds the m6A residue of MsYTHDC consists of a sequence of two tryptophan residues and one tyrosine residue (WWY). Conversely, in MsYTHDF, the binding pocket comprises two highly conserved tryptophan residues and either one tryptophan residue (WWW) or tyrosine residue (WWY) in MsYTHDF.Through comparative analysis of qRT-PCR data, we observed distinct expression patterns in specific genes under abiotic stress, indicating their potential regulatory roles. Notably, five genes (MsYTH2, 14, 26, 27, 48) consistently exhibit upregulation, and two genes (MsYTH33, 35) are downregulated in response to both cold and salt stress. This suggests a common mechanism among these YTH proteins in response to various abiotic stressors in alfalfa. Further, integrating qRT-PCR with RNA-seq data revealed that MsYTH2, MsYTH14, and MsYTH16 are highly expressed in leaves at various development stages, underscoring their potential roles in regulating the growth of these plant parts. The obtained findings shed further light on the biological functions of MsYTH genes and may aid in the selection of suitable candidate genes for future genetic enhancement endeavors aimed at improving salt and cold tolerance in alfalfa.