Zhenyuan Han1,2, Biao Yang3, Qin Wang4, Yuhua Hu1, Yuqiong Wu5,6, Zhen Tian7. 1. Department of Oral Pathology, Shanghai Ninth People's Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200011, China. 2. National Clinical Research Center for Oral Diseases, Shanghai, 200011, China. 3. Department of Neurosurgery, Huashan Hospital of Fudan University, Shanghai, 200040, China. 4. Clinical Translational Research Center, Shanghai Pulmonary Hospital, School of Life Sciences and Technology, Tongji University, Shanghai, 200092, China. 5. National Clinical Research Center for Oral Diseases, Shanghai, 200011, China. wyq192131@126.com. 6. Department of Prosthodontics, Shanghai Ninth People's Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200011, China. wyq192131@126.com. 7. Department of Oral Pathology, Shanghai Ninth People's Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200011, China. tianzhen_256@126.com.
Abstract
BACKGROUND: Invasive malignant pleomorphic adenoma (IMPA) is a highly invasive parotid gland tumor and lacks effective therapy. N6-Methyladenosine (m6A) is the most prevalent post-transcriptional modification of mRNAs in eukaryotes and plays an important role in the pathogenesis of multiple tumors. However, the significance of m6A-modified mRNAs in IMPA has not been elucidated to date. Hence, in this study, we attempted to profile the effect of IMPA in terms of m6A methylation in mRNA. METHODS: Methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were utilized to acquire the first transcriptome-wide profiling of the m6A methylome map in IMPA followed by bioinformatics analysis. RESULTS: In this study, we obtained m6A methylation maps of IMPA samples and normal adjacent tissues through MeRIP-seq. In total, 25,490 m6A peaks associated with 13,735 genes were detected in the IMPA group, whereas 33,930 m6A peaks associated with 18,063 genes were detected in the control group. Peaks were primarily enriched within coding regions and near stop codons with AAACC and GGAC motifs. Moreover, functional enrichment analysis demonstrated that m6A-containing genes were significantly enriched in cancer and metabolism relevant pathways. Furthermore, we identified a relationship between the m6A methylome and the RNA transcriptome, indicating a mechanism by which m6A modulates gene expression. CONCLUSIONS: Our study is the first to provide comprehensive and transcriptome-wide profiles to determine the potential roles played by m6A methylation in IMPA. These results may open new avenues for in-depth research elucidating the m6A topology of IMPA and the molecular mechanisms governing the formation and progression of IMPA.
BACKGROUND:Invasive malignant pleomorphic adenoma (IMPA) is a highly invasive parotid gland tumor and lacks effective therapy. N6-Methyladenosine (m6A) is the most prevalent post-transcriptional modification of mRNAs in eukaryotes and plays an important role in the pathogenesis of multiple tumors. However, the significance of m6A-modified mRNAs in IMPA has not been elucidated to date. Hence, in this study, we attempted to profile the effect of IMPA in terms of m6A methylation in mRNA. METHODS: Methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were utilized to acquire the first transcriptome-wide profiling of the m6A methylome map in IMPA followed by bioinformatics analysis. RESULTS: In this study, we obtained m6A methylation maps of IMPA samples and normal adjacent tissues through MeRIP-seq. In total, 25,490 m6A peaks associated with 13,735 genes were detected in the IMPA group, whereas 33,930 m6A peaks associated with 18,063 genes were detected in the control group. Peaks were primarily enriched within coding regions and near stop codons with AAACC and GGAC motifs. Moreover, functional enrichment analysis demonstrated that m6A-containing genes were significantly enriched in cancer and metabolism relevant pathways. Furthermore, we identified a relationship between the m6A methylome and the RNA transcriptome, indicating a mechanism by which m6A modulates gene expression. CONCLUSIONS: Our study is the first to provide comprehensive and transcriptome-wide profiles to determine the potential roles played by m6A methylation in IMPA. These results may open new avenues for in-depth research elucidating the m6A topology of IMPA and the molecular mechanisms governing the formation and progression of IMPA.