Deciphering the map of METTL14-mediated lncRNA m6A modification at the transcriptome-wide level in breast cancer.

Abstract

Emerging studies have demonstrated the critical role of RNA m6A methylation in tumor progression, whereas lncRNA m6A modification profiles in breast cancer remain largely unknown. Our previous study has shown that METTL14 accelerates breast cancer migration and invasion in an m6A-dependent manner, making it critical to analyze METTL14-mediated m6A modification at a transcriptome-wide scale in breast cancer. Here, we performed MeRIP-seq analysis in METTL14 overexpressed and control MDA-MB-231 cells. Conjoint analysis of MeRIP-seq and RNA-seq data was used to select lncRNAs with m6A methylation and differential expression. Finally, the screened lncRNA was verified by MeRIP-PCR and its function was studied via transwell assay. Our results determined that high expression of METLL14 results in 3996 hypermethylation peaks from 3107 transcripts, and 4100 hypomethylation peaks from 2918 transcripts. Furthermore, conjoint analysis of MeRIP-seq and RNA-seq data identified 25 lncRNAs with discrepant methylation and simultaneously discrepant expression, among which the top 10 differentially expressed LncRNAs were AC026401.3, CYTOR, LINC01943, AC084125.2, FLJ20021, LINC00472, and NORAD, MALAT1, AL161431.1, and LINC01764. Moreover, over-expressed METTL14 stimulated the m6A modification of AC084125.2, while decreasing its expression. Compared to adjacent tissues, AC084125.2 was lowly expressed in tumors and could be used as a biomarker in the diagnosis of breast cancer. Meanwhile, AC084125.2 inhibited the migration and invasion of cancer cells. In conclusion, METTL14-mediated m6A modification of lncRNAs, which might provide reference for future intervention in tumor progression.