Identification of m6A-Related Biomarkers Associated with Prognosis of Colorectal Cancer

BACKGROUND: Colorectal cancer (CRC) is the second most deadly cancer in the world according to GLOBOCAN 2020 data. Accumulating evidence suggests that RNA methylation modification is also misregulated in human cancers and may be a potential ideal target for cancer treatment. MATERIAL/METHODS: m6A-related differentially expressed genes (DEGs) were identified from colon adenocarcinoma and rectum adenocarcinoma esophageal carcinoma patients with different pathological stages. The protein-protein interaction (PPI) network construction, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs were conducted. Cox regression analysis was applied to the screening of m6A-related DEGs significantly associated with the overall survival (OS), and those selected genes were used for LASSO regression analysis to construct prognostic signature and calculate patients’ risk scores. RESULTS: We identified 673 m6A-related DEGs from CRC patients in different pathologic stages, and 146 of them were associated with OS. CTNNB1, TRIM37, RAB7A, CASC5/KNL1, CENPE, CCNB1, UBE2H, HSPA8, KIF1A, and FBXW4 were hub genes of the PPI network. Nine m6A-related genes were screened out to build the prognostic risk model. TNM stage, vascular invasion, and the risk score were independently related to the OS of CRC patients. CONCLUSIONS: Nine candidate m6A-related mRNA biomarkers (LRRC17, NFKB1, NOS2, PCDHB2, RAB7A, RPS6KA1, RRNAD1, TLE6, and UBE2H) were found to be closely related to the clinicopathology and prognosis of colorectal cancer, indicating that they could be potential prognostic biomarkers for patients with colorectal cancer.