Microarray and bioinformatic analysis reveal the parental genes of m6A modified circRNAs as novel prognostic signatures in colorectal cancer

BACKGROUND: Accumulating evidences have revealed that the abnormal N6-methyladenosine (m6A) modification is closely associated with the occurrence, development, progression and prognosis of cancer. It is noteworthy that m6A modification is widely existed in circRNAs and found its key biological functions in regulating circRNAs metabolism. However, the role of m6A modified circRNAs in colorectal cancer (CRC) remains unknown. To better understand the role of circRNAs in the pathogenesis of CRC, we focus on the relationship between m6A-modified circRNAs and their parental genes. METHODS: Arraystar m6A-circRNA epitranscriptomic microarray was used to identify differentially m6A modified circRNAs between CRC and the control group. In addition, TCGA-COAD and GSE106582 cohort were used to identify differentially expressed mRNAs. In this study, we screened the parental genes for which both circRNAs and mRNAs were down-regulated further to analyze, including gene expression, survival prognosis, enrichment analysis. Additionally, Western Blotting was used to further validate the role of the parental gene in CRC. RESULTS: We found that 1405 significantly downregulated circRNAs in CRC by our microarray data. Moreover, we obtained 113 parental genes for which both circRNAs and mRNAs were down-regulated to analyze the relationship with the prognosis of CRC based on TCGA-COAD cohort. And we identified nine potential prognostic genes, including ABCD3, ABHD6, GAB1, MIER1, MYOCD, PDE8A, RPS6KA5, TPM1 and WDR78. And low expression of these genes was associated with poor survival prognosis of the patients with CRC. In addition, we found that TPM1 is downregulated in CRC by western blotting experiment. And the calcium-signaling pathway may involve the process of the CRC progression. CONCLUSIONS: We identified nine potential prognostic genes, after analyzed the relationship between the parental genes of m6A modified circRNAs and the progression of CRC. Above all, our study further validated TPM1 can serve as a potentail signature for CRC patients.