Abstract
N6-methyladenosine (m6A) modification on RNA plays an important role in tumorigenesis and metastasis, which could change gene expression and even function at multiple levels such as RNA splicing, stability, translocation, and translation. In this study, we aim to conduct a comprehensive analysis on m6A RNA methylation-related genes, including m6A RNA methylation regulators and m6A RNA methylation-modified genes, in liver hepatocellular carcinoma, and their relationship with survival and clinical features. Data, which consist of the expression of widely reported m6A RNA methylation-related genes in liver hepatocellular carcinoma from The Cancer Genome Atlas (TCGA), were analyzed by one-way ANOVA, Univariate Cox regression, a protein–protein interaction network, gene enrichment analysis, feature screening, a risk prognostic model, correlation analysis, and consensus clustering analysis. In total, 405 of the m6A RNA methylation-related genes were found based on one-way ANOVA. Among them, DNA topoisomerase 2-alpha (TOP2A), exodeoxyribonuclease 1 (EXO1), ser-ine/threonine-protein kinase Nek2 (NEK2), baculoviral IAP repeat-containing protein 5 (BIRC5), hyaluronan mediated motility receptor (HMMR), structural maintenance of chromosomes protein 4 (SMC4), bloom syndrome protein (BLM), ca-sein kinase I isoform epsilon (CSNK1E), cytoskeleton-associated protein 5 (CKAP5), and inner centromere protein (INCENP), which were m6A RNA methylation-modified genes, were recognized as the hub genes based on the protein–protein interaction analysis. The risk prognostic model showed that gender, AJCC stage, grade, T, and N were significantly different between the subgroup with the high and low risk groups. The AUC, the evaluation parameter of the prediction model which was built by RandomForest, was 0.7. Furthermore, two subgroups were divided by consensus clustering analysis, in which stage, grade, and T differed. We identified the important genes expressed significantly among two clusters, including uridine-cytidine kinase 2 (UCK2), filensin (BFSP1), tubulin-specific chaperone D (TBCD), histone-lysine N-methyltransferase PRDM16 (PRDM16), phosphorylase b ki-nase regulatory subunit alpha (PHKA2), serine/threonine-protein kinase BRSK2 (BRSK2), Arf-GAP with coiled-coil (ACAP3), general transcription factor 3C polypep-tide 2 (GTF3C2), and guanine nucleotide exchange factor MSS4 (RABIF). In our study, the m6A RNA methylation-related genes in liver hepatocellular carcinoma were analyzed systematically, including the expression, interaction, function, and prognostic values, which provided an important theoretical basis for m6A RNA methylation in liver cancer. The nine important m6A-related genes could be prognostic markers in the survival time of patients.
Highlights
Liver hepatocellular carcinoma (LIHC) is one of the common malignant tumors of the digestive system, mainly caused by hepatitis virus infection
Isoform epsilon (CSNK1E), cytoskeleton-associated protein 5 (CKAP5), and inner centromere protein (INCENP), which were m6A RNA methylation-modified genes, were recognized as the hub genes based on the protein–protein interaction analysis
All 1190 m6A RNA methylation-related genes, which were related to LIHC, were selected with AJCC stage and RNA-seq expression data
Summary
Liver hepatocellular carcinoma (LIHC) is one of the common malignant tumors of the digestive system, mainly caused by hepatitis virus infection. It is a leading cause of cancerrelated deaths worldwide, with few effective therapeutic options [1]. Liver hepatocellular carcinoma is insidious and develops rapidly, and most patients have already progressed to the middle and late stages when they are diagnosed, and have lost the opportunity for surgery [2]. Finding new treatment and prognostic targets for LIHC will be beneficial to patients. With the deepening of research, epigenetic modifications such as DNA methylation, histone acetylation, and RNA modification have all been proved to participate in the development and progression of tumors, which have been recognized as new treatment and prognostic targets. The dynamic balance of m6A methylation modification is maintained in the cell by m6A RNA methylation writers, readers, and erasers [4]
Sign-in/Register to view highlights & summary Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
