![Open Access](/resources/images/iconopenaccess.png)
Novel prognostic implications of YTH domain family 2 in resected hepatocellular carcinoma
- Authors:
- Published online on: May 19, 2021 https://doi.org/10.3892/ol.2021.12799
- Article Number: 538
-
Copyright: © Nakagawa et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
Hepatocellular carcinoma (HCC) is the main type of primary liver cancer and is one of the most common malignancies with poor survival (1). Hepatectomy is a potentially curative treatment, but the recurrence rate of HCC after surgery is remarkably high at approximately 70% (2). Therefore, further understanding of the molecular mechanisms of HCC development and recurrence is required.
More than 100 types of chemical modifications have been identified in RNA (3). Recently, internal modifications of mRNA have received attention for their roles in mRNA metabolism. N6-methyladenosine (m6A) is the most abundant modification of mRNA in eukaryotes and was first reported in the 1970s (4). Recent evidence suggests that m6A has various functions in RNA metabolism, such as pre-mRNA splicing, 3′-end processing, nuclear export, translation regulation, regulation of mRNA decay and noncoding RNA processing (5–7). Furthermore, m6A methylation has been revealed to have crucial roles in the initiation and progression of cancer (8).
m6A readers are the proteins that recognize and bind to m6A sites and thereby elicit multiple effects (9). YT521-B homology domain family 2 (YTHDF2) was the first identified m6A binding protein (10). YTHDF2 weakens mRNA stability by recognizing m6A, while YT521-B homology domain family 1 (YTHDF1) promotes mRNA translation efficiency (11). In solid cancer, YTHDF2 and YTHDF1 have been reported to have roles as both tumour promoters and suppressors (12–14). However, the significance of YTHDF1 and YTHDF2 in oncogenesis remains unclear.
In the current study, we assessed the expression of YTHDF1 and YTHDF2 in both resected HCC tissues and paired normal liver tissues collected from patients who underwent surgery with curative intent. We also sought to discover novel prognostic implications of m6A readers that could be used to predict prognosis in patients with resected HCC.
Patients and methods
Patients and samples
A total of 177 frozen tumour specimens and paired paratumor noncancerous tissues were collected from patients with HCC who underwent surgery at Nagoya University Hospital (Nagoya City, Japan) between January 1998 and April 2014. All fresh tissues were immediately frozen in liquid nitrogen and stored at −80°C until use. Patient characteristics are summarized in Table I. After surgery, all patients were monitored via blood examinations, ultrasonography, and computed tomography once every six months. Angiography was performed for further information whenever recurrence was suspected. The median follow-up duration of all patients was 48.8 months (range, 0.3 to 191 months). This study and all procedures were approved by the Institutional Review Board at Nagoya University (Nagoya City, Japan), and all patients provided written informed consent. All clinical investigations were conducted in accordance with the principles of the Declaration of Helsinki.
RNA isolation and RT-qPCR
Total RNA was extracted from tissue samples using a Qiagen miRNeasy mini-kit (Qiagen, Hilden, Germany). We used DNAse to avoid contamination, and RNA quality was analysed by a NanoDrop (Thermo Scientific Fisher, Waltham, MA, USA). Total RNA was converted to complementary DNA by reverse transcription with M-MLV Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA). This total cDNA was used as a template for the next step of quantitative PCR (qPCR). qPCR was performed using SYBR Premix Ex Taq II (Takara Clontech, Kyoto, Japan) under the following conditions: Denaturation at 95°C for 10 sec and 40 cycles of denaturation at 95°C for 5 sec and annealing/extension at 60°C for 30 sec. The SYBR Green signal was detected in real-time using a StepOne Plus Real-Time PCR System (Life Technologies, Carlsbad, CA, USA). The relative quantification method was used, and the expression level of each gene was normalized to the expression level of the control gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) for each sample. The relative gene expression levels were determined using the comparative threshold cycle (2−ΔCT) method.
The PCR primers used in the current study were specific for the 78-base-pair fragment of YTHDF1 (sense, 5′-TCCATCTTCGACGACTTTGCT-3′; antisense, 5′-TCGACTCTGCCGTTCCTTG-3′) and for the 50-base-pair fragment of YTHDF2 (sense, 5′-GAGGATCTGAGAGCCATGTCG-3′; antisense, 5′-ATTTTGTACTGCTCCAAGAGGC-3′). GAPDH primers (sense, 5′-GAGTCCACTGGCGTCTTCAC-3′; antisense, 5′-GTTCACACCCATGACGAACA-3′) were used to quantify the expression in each sample as an internal control. The primers were designed as intron spanning. All qPCR experiments were performed in duplicate, including the template-omitted negative controls.
Acquisition of publicly available data
Normalized TCGA RNA-sequencing data of HCC were downloaded from the Broad GDAC Firehose (http://gdac.broadinstitute.org/, accessed on January 1st, 2020). This dataset consists of 50 noncancerous cases and 360 HCC cases, including seven HCC cases mixed with hepatocholangiocarcinoma and two cases with fibrolamellar carcinoma. Of the 360 cases, there were 266 cases with recurrence-free survival (RFS) information and 336 cases with overall survival (OS) information.
Statistical analysis
Continuous variables are expressed as the median (range), and the expression of each target gene was compared by a Wilcoxon signed-rank test and paired t-test. Categorical variables were compared using χ2 or Fisher's exact tests, as appropriate. The OS and RFS rates at each point of the follow-up time were estimated using the Kaplan-Meier method and compared using a log-rank test. A Cox proportional hazard regression model was used to perform univariate analysis and multivariate analysis for OS and RFS. In the multivariate analysis, variables that showed statistical significance in the univariate Cox proportional hazard regression were included. All statistical analyses were performed using R version 3.5.3 (http//www.r-project.org/), and P<0.05 obtained using two-tailed tests was considered to indicate statistical significance.
Results
YTHDF1 and YTHDF2 in resected specimens from HCC patients
First, expression analyses of m6A readers were conducted with our surgically resected specimens. The expression levels of YTHDF1 and YTHDF2 were measured by qPCR. The expression of YTHDF1 and YTHDF2 was not significantly different between tumour tissues and noncancerous tissues (P=0.93 and P=0.7, respectively, Fig. 1A). Fig. 1B shows individual changes in YTHDF1 and YTHDF2 expression in paired analysis. Based on the results obtained by qPCR, 177 HCC cases were divided into two groups according to YTHDF1 and YTHDF2 expression in tumour tissues. We selected the cut-off values that showed the best statistical difference. Clinical features of the groups stratified by YTHDF1 and YTHDF2 expression are shown in Table II. In HCC tissues, low YTHDF1 expression was associated with a lack of capsule (65% vs. 81%, P=0.02), whereas low YTHDF2 expression was associated with septal formation (73% vs. 53%, P=0.02).
![]() | Table II.Clinical features of 177 patients with hepatocellular carcinoma according to YTHDF1 and YTHDF2 expression. |
Prognostic significance of YTHDF1 and YTHDF2 in resected HCC cases
Next, the effects of the expression levels on RFS and OS were evaluated. In HCC tissues, both high YTHDF1 expression and high YTHDF2 expression were significantly correlated with shorter RFS (YTHDF1: MST=34.0 vs. 19.0 months, P=0.014; YTHDF2: MST=30.1 vs. 12.9 months, P=0.0032, Fig. 2A), whereas YTHDF1 expression and YTHDF2 expression were not correlated with OS (YTHDF1: MST=99.4 vs. 70.2 months, P=0.74; YTHDF2: MST=98.4 vs. 64.1 months, P=0.28, Fig. 2B).
YTHDF1 and YTHDF2 expression levels and their correlation with HCC prognosis in a publicly available dataset
We analysed the expression levels of YTHDF1 and YTHDF2 in HCC and noncancerous tissues using a TCGA RNA-sequence dataset. This analysis revealed that the expression of YTHDF1 was significantly higher in HCC tumour tissues and that YTHDF2 expression was significantly lower in HCC tumour tissues than in noncancerous tissues (YTHDF1, P=0.0011; YTHDF2, P=0.025, Fig. 3A and B). We then confirmed the prognostic impact of YTHDF1 and YTHDF2 expression in resected HCC patients using the same TCGA dataset. HCC cases were divided into two groups according to the YTHDF1 and YTHDF2 expression in HCC tissues in the normalized RNA-sequencing data. We also selected the cut-off values that showed the best statistical difference. This analysis revealed that the patients with high YTHDF1 expression had significantly worse RFS (MST=754 vs. 489 days, P=0.006, Fig. 3C), and the patients with high YTHDF2 expression tended to have worse RFS (MST=636 vs. 315 days, P=0.06, Fig. 3D). In addition, the patients with high YTHDF1 expression had significantly worse OS (MST=2456 vs. 1,372 days, P=0.015, Fig. 3E), and the patients with high YTHDF2 expression had significantly worse OS (MST=2,456 vs. 724 days, P=0.0002, Fig. 3F).
Cox regression analysis of HCC survival
Since the survival curves showed that YTHDF1 and YTHDF2 expression levels in HCC tissues were correlated with RFS, we performed Cox proportional hazards analyses to further investigate the prognostic value of YTHDF1 and YTHDF2 expression. The multivariate analysis identified serosal invasion (hazard ratio (HR): 2.39, 95% confidence interval (95% CI): 1.30–4.42, P=0.005), portal vein or hepatic vein invasion (HR, 2.82, 95% CI: 1.26–6.28, P=0.01) and YTHDF2 expression in HCC tissues (HR, 1.85, 95% CI: 1.09–3.15, P=0.02) as significant independent factors for RFS (Table III) and AFP (HR, 1.79, 95% CI: 1.10–2.92, P=0.02), serosal invasion (HR, 1.99, 95% CI: 1.17–3.34, P=0.01), and portal vein or hepatic vein invasion (HR, 3.02, 95% CI: 1.38–6.61, P=0.006) as significant independent factors for OS (Table IV). Consequently, high expression of YTHDF2 in HCC tissues was significantly associated with recurrence after HCC surgery.
![]() | Table III.Univariate and multivariate Cox proportional-hazard regression analysis of recurrence free survival in patients with hepatocellular carcinoma. |
![]() | Table IV.Univariate and multivariate cox proportional-hazard regression analysis of overall survival in patients with hepatocellular carcinoma. |
Discussion
In this study, we primarily evaluated the clinical effects of m6A readers in resected HCC patients. The expression of YTHDF1 and YTHDF2 in HCC tissues was correlated with tumour recurrence. Furthermore, YTHDF2 was an independent prognostic factor in resected HCC patients. Members of the YT521-B homology (YTH) domain family, including YTHDF1, YTHDF2, YTHDF3, YTHDC1 and YTHDC2, all have a conserved m6A-binding domain and preferentially bind to m6A-modified RNA at the RRm6ACH consensus sequence (15). YTHDF2, the first characterized m6A reader, accelerates the decay of m6A-modified transcripts by facilitating the recruitment of the CCR4-NOT complex directly (10). In contrast, YTHDF1 was initially demonstrated to bind to m6A sites around the stop codon and then cooperate with the translation initiation machinery to improve the translation efficacy of target RNAs in mammals (11).
We first compared the expression levels of YTHDF1 and YTHDF2 in HCC and noncancerous liver tissues from patients who underwent hepatectomy with curative intent at our institution. Neither YTHDF1 expression nor YTHDF2 expression was significantly different between HCC tissues and noncancerous tissues. In the paired analysis, patients with high expression in tumour tissues tended to have low expression in noncancerous tissues. In another cohort in the TCGA RNA-sequence dataset, however, the expression of YTHDF1 was found to be significantly higher in HCC tissues, and YTHDF2 was significantly lower in HCC tissues. Li et al reported that upregulation of YTHDF2 was observed in TCGA prostate cancer tissues compared with normal controls (16). Bai et al also reported that YTHDF1 is significantly upregulated in tumour compared with adjacent normal tissues in colorectal cancer (17). In the TCGA dataset, there were relatively small numbers of noncancerous tissues available, but we studied the expression of YTHDF1 and YTHDF2 in both tumour and noncancerous tissues from 177 HCC patients. This might have caused the discrepancy in results between the TCGA cases and our cases. In addition, our data showed that the reader expression in tumour tissues was different from that in noncancerous tissues. This change could be more important than the absolute value. In addition, the prognostic analysis stratified by YTHDF1 and YTHDF2 expression in our study revealed that the expression of these two m6A readers in HCC tissues is not associated with OS. On the other hand, high expression of both YTHDF1 and YTHDF2 in HCC tissues was associated with significantly worse RFS. In the public dataset, high expression of YTHDF1 or YTHDF2 was associated with worse prognosis than low expression. In particular, our Cox regression analysis showed that YTHDF2 was an independent risk factor for recurrence in resected HCC. Thus, YTHDF1 and YTHDF2 might have inherent effects in HCC carcinogenesis and influence the long-term outcome after HCC resection, for example, by causing sporadic recurrence.
Evidence of RNA modifications in cancer development and progression has been increasing. The RNA methyltransferase METTL3 is the first characterized component of the m6A methyltransferase complex. METTL3 promotes tumour proliferation and invasion in several cancers (18–22). The m6A demethylases FTO and ALKBH5 were identified in the 2010s. FTO and ALKBH5 also play an important role in human cancer (23–25). However, the functions of YTHDF1 and YTHDF2 in HCC have not been uncovered. Zhao et al reported that YTHDF1 played a vital role in the regulation of HCC metabolism (26). Qu et al reported that m6A RNA methylation modulators, including YTHDF1, affected OS in HCC patients (27). YTHDF2 was able to degrade both tumour promoter and suppressor gene mRNAs. Zhang et al reported that YTHDF2 promotes the cancer stem cell liver phenotype and cancer metastasis by modulating the m6A methylation of OCT4 mRNA (28). In contrast, YTHDF2 may act as a tumour suppressor to repress cell proliferation and growth by destabilizing EGFR mRNA in HCC (29). Further investigations are required to reveal the role of YTHDF1 and YTHDF2 in HCC.
Although we showed important aspects of YTHDF1 and YTHDF2, there are some inherent limitations to the present study. First, more data are necessary because we used specimens from a single institute in this study. Second, more detailed molecular mechanisms through which specific m6A methylation enhances HCC development need to be discovered. Further investigation is necessary before the clinical utility of our findings can be determined.
In conclusion, our study revealed that high YTHDF2 expression in HCC tissues is related to cancer recurrence. Our results may pave the way for discovering the clinical utility of m6A methylation and associated genes in HCC therapy in the future.
Acknowledgements
The authors would like to thank Dr Raju Kandimalla (Bristol Myers Squibb, San Diego, CA, USA) for useful discussions and Ms. Yoko Nishikawa (Nagoya University Graduate School of Medicine, Nagoya, Japan) for technical assistance with the experiments.
Funding
The present study was partly supported by the Japan Society for the Promotion of Science KAKENHI Grant-in-Aid for Scientific Research (grant nos. 18H06176 and 19K18110).
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Authors' contributions
NN, FS and YK conceived and designed the study. FS and YK provided financial support. SY and YK provided administrative support. SY, FS, YS, YI, HT, MH, MKa, CT, GN and MKo provided study materials and patients. SY, FS, YS, YI, HT, MH, MKa, CT, GN and MKo assisted with analyses and manuscript preparation. NN, FS, KT and YS collected and assembled the data. NN and KT performed quantitative PCR data analysis and interpretation. NN, FS and YK wrote the manuscript and confirmed the authenticity of all the raw data. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The study and all procedures were approved by the Institutional Review Board at Nagoya University (approval no. 2013-0295; Nagoya, Japan), and all patients provided written informed consent. All clinical investigations were conducted in accordance with the principles of the Declaration of Helsinki.
Patient consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
References
Siegel RL, Miller KD and Jemal A: Cancer statistics, 2020. CA Cancer J Clin. 70:7–30. 2020. View Article : Google Scholar : PubMed/NCBI | |
Tabrizian P, Jibara G, Shrager B, Schwartz M and Roayaie S: Recurrence of hepatocellular cancer after resection: Patterns, treatments, and prognosis. Ann Surg. 261:947–955. 2015. View Article : Google Scholar : PubMed/NCBI | |
Roundtree IA, Evans ME, Pan T and He C: Dynamic RNA modifications in gene expression regulation. Cell. 169:1187–1200. 2017. View Article : Google Scholar : PubMed/NCBI | |
Adams JM and Cory S: Modified nucleosides and bizarre 5′-termini in mouse myeloma mRNA. Nature. 255:28–33. 1975. View Article : Google Scholar : PubMed/NCBI | |
Camper SA, Albers RJ, Coward JK and Rottman FM: Effect of undermethylation on mRNA cytoplasmic appearance and half-life. Mol Cell Biol. 4:538–543. 1984. View Article : Google Scholar : PubMed/NCBI | |
Finkel D and Groner Y: Methylations of adenosine residues (m6A) in pre-mRNA are important for formation of late simian virus 40 mRNAs. Virology. 131:409–425. 1983. View Article : Google Scholar : PubMed/NCBI | |
Gilbert WV, Bell TA and Schaening C: Messenger RNA modifications: Form, distribution, and function. Science. 352:1408–1412. 2016. View Article : Google Scholar : PubMed/NCBI | |
Ma S, Chen C, Ji X, Liu J, Zhou Q, Wang G, Yuan W, Kan Q and Sun Z: The interplay between m6A RNA methylation and noncoding RNA in cancer. J Hematol Oncol. 12:1212019. View Article : Google Scholar : PubMed/NCBI | |
Allis CD and Jenuwein T: The molecular hallmarks of epigenetic control. Nat Rev Genet. 17:487–500. 2016. View Article : Google Scholar : PubMed/NCBI | |
Du H, Zhao Y, He J, Zhang Y, Xi H, Liu M, Ma J and Wu L: YTHDF2 destabilizes m(6)A-containing RNA through direct recruitment of the CCR4-NOT deadenylase complex. Nat Commun. 7:126262016. View Article : Google Scholar : PubMed/NCBI | |
Wang X, Zhao BS, Roundtree IA, Lu Z, Han D, Ma H, Weng X, Chen K, Shi H and He C: N(6)-methyladenosine modulates messenger RNA translation efficiency. Cell. 161:1388–1399. 2015. View Article : Google Scholar : PubMed/NCBI | |
Chen J, Sun Y, Xu X, Wang D, He J, Zhou H, Lu Y, Zeng J, Du F, Gong A and Xu M: YTH domain family 2 orchestrates epithelial-mesenchymal transition/proliferation dichotomy in pancreatic cancer cells. Cell Cycle. 16:2259–2271. 2017. View Article : Google Scholar : PubMed/NCBI | |
Nishizawa Y, Konno M, Asai A, Koseki J, Kawamoto K, Miyoshi N, Takahashi H, Nishida N, Haraguchi N, Sakai D, et al: Oncogene c-Myc promotes epitranscriptome m6A reader YTHDF1 expression in colorectal cancer. Oncotarget. 9:7476–7486. 2018. View Article : Google Scholar : PubMed/NCBI | |
Shi Y, Fan S, Wu M, Zuo Z, Li X, Jiang L, Shen Q, Xu P, Zeng L, Zhou Y, et al: YTHDF1 links hypoxia adaptation and non-small cell lung cancer progression. Nat Commun. 10:48922019. View Article : Google Scholar : PubMed/NCBI | |
Wang X, Lu Z, Gomez A, Hon GC, Yue Y, Han D, Fu Y, Parisien M, Dai Q, Jia G, et al: N6-methyladenosine-dependent regulation of messenger RNA stability. Nature. 505:117–120. 2014. View Article : Google Scholar : PubMed/NCBI | |
Li J, Xie H, Ying Y, Chen H, Yan H, He L, Xu M, Xu X, Liang Z, Liu B, et al: YTHDF2 mediates the mRNA degradation of the tumor suppressors to induce AKT phosphorylation in N6-methyladenosine-dependent way in prostate cancer. Mol Cancer. 19:1522020. View Article : Google Scholar : PubMed/NCBI | |
Bai Y, Yang C, Wu R, Huang L, Song S, Li W, Yan P, Lin C, Li D and Zhang Y: YTHDF1 regulates tumorigenicity and cancer stem cell-like activity in human colorectal carcinoma. Front Oncol. 9:3322019. View Article : Google Scholar : PubMed/NCBI | |
Han J, Wang JZ, Yang X, Yu H, Zhou R, Lu HC, Yuan WB, Lu JC, Zhou ZJ, Lu Q, et al: METTL3 promote tumor proliferation of bladder cancer by accelerating pri-miR221/222 maturation in m6A-dependent manner. Mol Cancer. 18:1102019. View Article : Google Scholar : PubMed/NCBI | |
Lin S, Choe J, Du P, Triboulet R and Gregory RI: The m(6)A Methyltransferase METTL3 promotes translation in human cancer cells. Mol Cell. 62:335–345. 2016. View Article : Google Scholar : PubMed/NCBI | |
Miao W, Chen J, Jia L, Ma J and Song D: The m6A methyltransferase METTL3 promotes osteosarcoma progression by regulating the m6A level of LEF1. Biochem Biophys Res Commun. 516:719–725. 2019. View Article : Google Scholar : PubMed/NCBI | |
Wang Q, Chen C, Ding Q, Zhao Y, Wang Z, Chen J, Jiang Z, Zhang Y, Xu G, Zhang J, et al: METTL3-mediated m6A modification of HDGF mRNA promotes gastric cancer progression and has prognostic significance. Gut. 69:1193–1205. 2020. View Article : Google Scholar : PubMed/NCBI | |
Xia T, Wu X, Cao M, Zhang P, Shi G, Zhang J, Lu Z, Wu P, Cai B, Miao Y and Jiang K: The RNA m6A methyltransferase METTL3 promotes pancreatic cancer cell proliferation and invasion. Pathol Res Pract. 215:1526662019. View Article : Google Scholar : PubMed/NCBI | |
Li J, Han Y, Zhang H, Qian Z, Jia W, Gao Y, Zheng H and Li B: The m6A demethylase FTO promotes the growth of lung cancer cells by regulating the m6A level of USP7 mRNA. Biochem Biophys Res Commun. 512:479–485. 2019. View Article : Google Scholar : PubMed/NCBI | |
Niu Y, Lin Z, Wan A, Chen H, Liang H, Sun L, Wang Y, Li X, Xiong XF, Wei B, et al: RNA N6-methyladenosine demethylase FTO promotes breast tumor progression through inhibiting BNIP3. Mol Cancer. 18:462019. View Article : Google Scholar : PubMed/NCBI | |
Zhang S, Zhao BS, Zhou A, Lin K, Zheng S, Lu Z, Chen Y, Sulman EP, Xie K, Bögler O, et al: m6A Demethylase ALKBH5 maintains tumorigenicity of glioblastoma stem-like cells by sustaining FOXM1 expression and cell proliferation program. Cancer Cell. 31:591–606.e6. 2017. View Article : Google Scholar : PubMed/NCBI | |
Zhao X, Chen Y, Mao Q, Jiang X, Jiang W, Chen J, Xu W, Zhong L and Sun X: Overexpression of YTHDF1 is associated with poor prognosis in patients with hepatocellular carcinoma. Cancer Biomark. 21:859–868. 2018. View Article : Google Scholar : PubMed/NCBI | |
Qu N, Qin S, Zhang X, Bo X, Liu Z, Tan C, Wen G and Jiang H: Multiple m6A RNA methylation modulators promote the malignant progression of hepatocellular carcinoma and affect its clinical prognosis. BMC Cancer. 20:1652020. View Article : Google Scholar : PubMed/NCBI | |
Zhang C, Huang S, Zhuang H, Ruan S, Zhou Z, Huang K, Ji F, Ma Z, Hou B and He X: YTHDF2 promotes the liver cancer stem cell phenotype and cancer metastasis by regulating OCT4 expression via m6A RNA methylation. Oncogene. 39:4507–4518. 2020. View Article : Google Scholar : PubMed/NCBI | |
Zhong L, Liao D, Zhang M, Zeng C, Li X, Zhang R, Ma H and Kang T: YTHDF2 suppresses cell proliferation and growth via destabilizing the EGFR mRNA in hepatocellular carcinoma. Cancer Lett. 442:252–261. 2019. View Article : Google Scholar : PubMed/NCBI |