Cancer Biother Radiopharm. 2020 Oct 14.
Purpose: N6-methyladenosine (m6A) methylation was the most abundant internal modification on messenger RNAs in eukaryotes. This study intended to explore the role of m6A methylation in endometrial cancer (EC). Materials and Methods: The m6A-sequencing data "GSE93911" of human EC were downloaded from Gene Expression Omnibus database. Hisat2 software and MACS2 were used to perform the alignment of reads and m6A methylation peak calling, and the peaks were annotated using Chipseeker. Then, differential m6A methylation peaks between normal and tumor samples were analyzed, followed by the functional enrichment analysis of the differentially methylated genes in promoter and 3' untranslated region (UTR) using Clusterprofiler. Based on the 450K methylated chip data, gene expression and clinical data in The Cancer Genome Atlas, the differentially methylated genes were verified, followed by Cox univariate/multivariate regression analysis and survival analysis. Finally, a risk prognosis model was constructed. Results: The m6A peak number was decreased in EC. The distribution of m6A peaks was highly enriched near transcriptional start site, in promoter, UTR, intron and exon, followed by distal intergenic. A total of 581 differentially methylated genes (361 hyper- and 220 hypomethylated genes) were identified in promoter and UTR regions that were enriched in insulin resistance (IR) and extracellular matrix (ECM). A total of 181 genes with significant differential expressions and differential methylation site in EC were selected. Of which, 31 genes were correlated with survival, and an 11-gene risk prognosis model was identified, including GDF7, BNC2, SLC8A1, B4GALNT3, DHCR24, ESRP1, HOXB9, IGSF9, KIAA1324, MSnX1, and PHGDH. Conclusion: The m6A methylation regulated EC progression by targeting the genes related to IR and ECM. A 11-gene risk prognosis model was identified to predict survival of patients with EC.
Keywords: Msh homeobox 1; N6-methyladenosine methylation; endometrial cancer