Deciphering the Divergent Gene Expression Landscapes of m6A/m5C/m1A Methylation Regulators in Hepatocellular Carcinoma Through Single-Cell and Bulk RNA Transcriptomic Analysis
DOI: https://doi.org/10.2147/jhc.s448047
2023-12-28
Journal of Hepatocellular Carcinoma
Abstract:Hang-Tsung Liu, 1, &ast Cheng-Shyuan Rau, 2, &ast Yueh-Wei Liu, 3 Ting-Min Hsieh, 1 Chun-Ying Huang, 1 Peng-Chen Chien, 4 Hui-Ping Lin, 4 Chia-Jung Wu, 4 Pei-Chin Chuang, 5 Ching-Hua Hsieh 4 1 Department of Trauma Surgery, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung, 83301, Taiwan; 2 Department of Neurosurgery, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung, 83301, Taiwan; 3 Department of General Surgery, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung, 83301, Taiwan; 4 Department of Plastic Surgery, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung, 83301, Taiwan; 5 Department of Medical Research, Kaohsiung Chang Gung Memorial Hospital, Kaohsiung, 83301, Taiwan &astThese authors contributed equally to this work Correspondence: Pei-Chin Chuang; Ching-Hua Hsieh, Tel +886-7-3454746, Email ; Introduction: RNA modifications mediated by the m6A, m1A, and m5C regulatory genes are crucial for the progression of malignancy. This study aimed to explore the expression of regulator genes for m6A/m5C/m1A methylation at the single-cell level and to validate their expression in cancerous and adjacent para-cancerous liver tissues of adult patients with HCC who underwent tumor resection. Methods: The bulk sequencing from The Cancer Genome Atlas (TCGA) database and the single-cell RNA sequencing (scRNA-seq) data obtained from the Gene Expression Omnibus (GEO) database were used to identify the dysregulated m6A/m5C/m1A genes for hepatocellular carcinoma (HCC). A real-time polymerase chain reaction (real-time PCR) was used to measure the expression of dysregulated m6A/m5C/m1A genes in collected human HCC tissues and compared with adjacent para-cancerous liver tissues. Immune cell infiltration with these significantly expressed methylation-related genes was evaluated using Timer2.0. Results: A discrepancy in m6A/m5C/m1A gene expression was observed between bulk sequencing and scRNA-seq. The clustered heatmap of the scRNA-seq-identified dysregulated m6A/m5C/m1A genes in TCGA cohort revealed heterogeneous expression of these methylation regulators within the cancer, whereas their expression in the adjacent liver tissues was more homogeneous. The real-time PCR validated the significant overexpression of DNMT1, NSUN5, TRMT6, IGF2BP1, and IGFBP3, which were identified using scRNA-seq, and IGFBP2, which was identified using bulk sequencing. These dysregulated methylation genes are mainly correlated with the infiltration of natural killer cells. Discussion: This study suggests that cellular diversity inside tumors contributes to the discrepancy in the expression of methylation regulator genes between traditional bulk sequencing and scRNA-seq. This study identified five regulatory genes that will be the focus of further studies regarding the function of m6A/m5C/m1A in HCC. Keywords: hepatocellular carcinoma, HCC, methylation, N6-methyladenosine, m6A, N1-methyladenosine, m1A, 5-methylcytidine, m5C, bulk sequencing, single-cell RNA sequencing, scRNA-seq Hepatocellular carcinoma (HCC) is the most common primary liver cancer worldwide, and is distinguished by its rapid development, high invasiveness, and high mortality rate. 1 Apart from well-explored epigenetic modifications such as DNA methylation, RNA modifications have emerged as significant contributors to the regulation of gene expression and cellular functions. These post-transcriptional modifications dynamically affect various aspects of RNA metabolism, including RNA processing, export, translation, and stability, and play crucial roles in cancer behaviors. 2–8 Research has demonstrated that the development of HCC is governed by intricate genetic and epigenetic mechanisms, in which abnormal RNA modifications play a pivotal role in causing variations in gene expression, thereby influencing the occurrence and progression of HCC. 9–15 Methylation involves a complex and dynamically reversible series of structural alterations induced by various enzymes. The discovery of N6-methyladenosine (m6A) has sparked interest in understanding its biological significance. 16,17 Research has revealed that m6A plays essential roles in various cellular processes, including RNA splicing and alternative splicing, RNA translation and stability, and regulation of gene expression. 18,19 Additionally, the presence of 5-methylcytidine (m5C) is markedly lower than that of m6A, with levels appro -Abstract Truncated-
oncology