- Open Access
m6A mRNA methylation regulates CTNNB1 to promote the proliferation of hepatoblastoma
Molecular Cancer volume 18, Article number: 188 (2019)
- The Correction to this article has been published in Molecular Cancer 2020 19:24
N6-Methyladenosine (m6A) modification has been implicated in many biological processes. It is important for the regulation of messenger RNA (mRNA) stability, splicing, and translation. However, its role in cancer has not been studied in detail. Here we investigated the biological role and underlying mechanism of m6A modification in hepatoblastoma (HB).
We used Reverse transcription quantitative real-time PCR (RT-qPCR) and Western blotting to determine the expression of m6A related factors. And we clarified the effects of these factors on HB cells using cell proliferation assay, colony formation, apoptotic assay. Then we investigated of methyltransferase-like 13 (METTL3) and its correlation with clinicopathological features and used xenograft experiment to check METTL3 effect in vivo. m6A-Seq was used to profiled m6A transcriptome-wide in hepatoblastoma tumor tissue and normal tissue. Finally, methylated RNA immunoprecipitation (MeRIP) assay, RNA remaining assay to perform the regulator mechanism of MEETL3 on the target CTNNB1 in HB.
In this research, we discovered that m6A modifications are increased in hepatoblastoma, and METTL3 is the main factor involved with aberrant m6A modification. We also profiled m6A across the whole transcriptome in hepatoblastoma tumor tissues and normal tissues. Our findings suggest that m6A is highly expressed in hepatoblastoma tumors. Also, m6A is enriched not only around the stop codon, but also around the coding sequence (CDS) region. Gene ontology analysis indicates that m6A mRNA methylation contributes significantly to regulate the Wnt/β-catenin pathway. Reduced m6A methylation can lead to a decrease in expression and stability of the CTNNB1.
Overall our findings suggest enhanced m6A mRNA methylation as an oncogenic mechanism in hepatoblastoma, METTL3 is significantly up-regulated in HB and promotes HB development. And identify CTNNB1 as a regulator of METTL3 guided m6A modification in HB.
Hepatoblastoma (HB) is the most common pediatric liver cancer. It is an embryonal neoplasms and mostly can be diagnosed during the first three years of life. It originates from undifferentiated hepatic progenitor cells, and undergo a malignant transformation during embryogenesis [1, 2]. Typical therapeutic strategies such as combined surgery and chemotherapy have demonstrated improved outcomes for children with HB. However, the prognosis for patients with advanced or chemotherapy-refractory disease is still very poor [3,4,5]. Components of the Wnt/β-catenin pathway are frequently mutated and overactive in solid malignancies and promote tumor development . In case of HB, CTNNB1 encoding β-catenin, is the most recurrently mutated driving proto-oncogene gene with 50–90% frequency. The point mutations and in-frame deletion of the exon3 in CTNNB1, has been reported as the primary cause of HB. In-frame deletions or missense mutations within exon 3 are gain-of-function mutations can lead to a degradation-resistant β-catenin protein that accumulates in the nucleus, binds to the TCF4/LEF − 1 transcription factor, and drives the activation of target genes such as Jun, c-Myc and Cyclin D1 [7,8,9]. Moreover, products of the loss-of-function somatic mutations of tumor suppressor genes AXIN1 and AXIN2 affecting the β-catenin degradation have also been reported in HB .
m6A is considered as the most common internal modification in messenger RNAs (mRNAs). These modifications are generally present on adenosine embedded in the consensus sequence G [G > A]m6 AC[U > A > C] across the mRNA transcripts, especially near the stop codons . Recent studies have demonstrated that m6A modifications in mRNAs or non-coding RNAs can influence RNA fate and functions such as mRNA stability, splicing, transport, localization and translation, primary microRNA processing, and RNA-protein interactions [12,13,14]. Also, m6A modifications are critical for most bioprocesses including tissue development, Deoxyribonucleic acid (DNA) damage response, sex determination, and tumorigenesis [11, 15, 16]. The reversible process of m6A formation occurs through a multicomponent methyltransferase complex consisting of the methyltransferase-like 13 (METTL3)/methyltransferase like 14 (METTL14) heterodimeric catalytic core and a regulatory subunit, Wilms’ tumor 1-associating protein (WTAP), KIAA1429, RNA Binding Motif Protein 15 (RBM15) and Zinc finger CCCH domain-containing protein 13 (ZC3H13), [17, 18]. METTL3, METTL14, and WTAP constitute the core of the methyltransferase complex. METTL3 possesses an active methyltransferase domain that catalyzes the conversion of adenosine (A) to m6A, METTL14 acts as the RNA-binding platform, and therefore, the METTL3-METTL14 forms a stable heterodimer and highly conserved in mammals. The heterodimer is required for the methylation activity, WTAP interacts with METTL3 and METTL14 to modulate levels of RNA transcripts [19,20,21,22]. The reversible process is conducted by m6A eraser fat mass and obesity-associated protein (FTO) and alkB homolog 5 (ALKBH5). The FTO and ALKBH5 serve as two mammalian RNA demethylases to reverse the m6A modifications [23,24,25]. m6A modification needs to be recognized by reader proteins in order to exert their biological functions. YT521-B homology (YTH) domain family proteins are reported to be m6A readers. Among these, YTHDF1 as the reader, located at the 5′ untranslated region of m6A has been suggested to promote mRNA translation efficiency in a cap-independent manner . Another m6A reader YTHDF2 binds to m6A at the 3′ untranslated region through its C-terminal YTD domain and localizes the targeted mRNA to processing, bodies for accelerated degradation through its N-terminal domain . Furthermore, heterogeneous nuclear ribonucleoprotein (hnRNP) could also serve as a potential nuclear m6A reader that can influence mRNA localization and alternative splicing . Given the crucial role of m6A modification in regulating gene expression and various biological processes, it is reasonable to speculate that aberrant m6A modification might also be involved in carcinogenesis. Previous research had shown that m6A eraser protein ALKBH5 has the function of regulating cancer stem cell properties . FTO has been reported to serve as an oncogene in acute myeloid leukemia . METTL3 depletion can leads to apoptosis through the p53 signaling, also it can promote liver cancer (HCC) progression through YTHDF2 dependent post-transcriptional silencing of suppressor of cytokine signaling 2 (SOCS2), [29, 30]. However, research on the role of METTL3 and the m6A machinery in cancers including hepatoblastoma still in its early stage, indicating an urgent need to head light on future clinical cancer therapy and drug developments.
In the present study we assessed levels of m6A, its key catalytic genes and profiled m6A transcriptome-wide in HB tissues. The level of m6A was increased in HB, and METTL3 was found to be the main factor involved in aberrant m6A modification. Furthermore, we also explored the target gene CTNNB1 and addressed the mechanism of methyltransferase METTL3 participates modulate CTNNB1 in HB.
Materials and methods
HB specimens and cell lines
Primary HB tumors and their nontumor counterparts from 25 patients were used. The study was approved by the Shanghai Children’s Medical Center affiliated to Shanghai Jiaotong University School of Medicine Institutional Review Board. Informed consent was obtained from each participant. Human HB cell line Huh6 was purchased from Cell Bank of Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China). Human HB cell line HepG2, human normal liver cell lines (QSG-7701, HL-7702) and HEK293T were purchased from American Type Culture Collection (ATCC) (Maryland, USA). HepG2 cells were cultured in Minimum Essential Medium (MEM) (Gibco, USA) supplemented with 10% (v/v) Fetal Bovine Serum (FBS) (Gibco, USA) and 1% (v/v) antibiotics (penicillin streptomycin) (Gibco, USA), Huh6, HEK293T and HL-7702 cells were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM) (Gibco, USA) with 10% Fetal Bovine Serum (FBS) (Gibco, USA) and 1% (v/v) antibiotics (penicillin streptomycin) (Gibco, USA), QSG-7701 cells were cultured in RPMI-1640 Medium (Gibco, USA) supplemented with 10% (v/v) Fetal Bovin Serum (FBS) (Gibco, USA) and 1% (v/v) antibiotics (penicillin streptomycin) (Gibco, USA). Cells were cultured at 37 °C in an atmosphere with 5% CO2.
Measurement of m6A level
m6A level in total RNA were detected using a commercial m6A RNA methylation quantification kit (Abcam, UK). Briefly, 200 ng total RNA was added into each well and the capture antibody solution and detection antibody solution were added according to the manufacturer’s protocol. The absorbance of each well at a wavelength of 450 nm was colorimetrically measured for m6A level.
High throughput m6A sequencing
Total RNA was extracted from HB tumor and normal tissues. RNA was tested for quality using nanodrop and gel electrophoresis. RNA was randomly fragmented to ~ 200 nucleotides by RNA fragmentation reagents. Fragmented RNA was incubated for 2 h at 4 C with m6A antibody (Synaptic System, Cat. No. 202003, diluted to 0.5 μg/ul) for immunoprecipitation following the standard protocol of the Magna methylated RNA immune precipitation (MeRIP) m6A kit (Merk Millipore, MA). The mixture was then incubated with beads and eluted with elution buffer (1 × IP buffer and 6.7 mM m6A). Eluted RNA was precipitated by 75% ethanol. The eluted RNA was treated with RNasin (Ambion, Cat No. AM2694) according to the manufacturer’s instructions. TruSeq Stranded mRNA Sample Pre Kit (Illumina) was used to construct the library from immunoprecipitated RNA and input RNA following to a published protocol. Sequencing was performed on an Illumina HiSeq machine with 2 × 100 cycles Solexa paired-end sequencing.
RNA extraction and quantitative PCR (qPCR) analysis
Total RNA was isolated using TRIzol reagent (Invitrogen, USA) following the manufacturer’s protocol and quantified by nanodrop 2000. Briefly, cells were lysed with TRIzol reagent at room temperature for more than 5 min and chloroform was added. After centrifugation at 12,000 × rpm for 15 min the aqueous phase was collected and mixed with isopropanol before centrifuging at 12,000 × rpm for 10 min at 4 °C. RNA was dissolved in RNase-free water and recovered by centrifuge. For analysis of mRNA expression, 200–500 ng of RNA was converted to cDNA using PrimeScript™ RT Regent Kit (TAKARA, Dalian, China). Quantitative real-time PCR using KAPA SYBR® FAST qPCR Kit Master Mix (2×) Universal (Applied Biosystems, USA) was performed on a 7500 Fast Real-time PCR system (Applied Biosystems). Total RNA was isolated for qPCR analysis. Quantitative PCR primers sequence are listed in Additional file 3: Table S1. GAPDH was used as a normalizing gene in all experiments.
Western blot assay
Equal amounts (around 0.2 g) of tissues and counted cells were homogenized in RIPA buffer (0.1 g/mL) (Beyotime, China) by using a Speed-Mill PLUS homogenizer (Analytik Jena, Jena, Germany). The homogenate was then centrifuged at 12,000 × rpm for 20 mins at 4 °C, and the supernatant proteins were collected. After the measurement of protein concentration, the proteins were boiled in 1× sodium dodecyl sulfate (SDS) sample buffer (100 μl for 1× 106 cells) for 5 mins, and then equal amounts of protein (50 μg) were loaded and separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel and transferred to polyvinylidene fluoride (PVDF) membrane (GE Healthcare, Buckinghamshire, UK). The membrane was then incubated with specific primary antibodies (Additional file 4: Table S2) and fluorescently conjugated secondary antibodies (Licor, USA), followed by the detection with the Odyssey system (TAITEC Co., Saitama, Japan).
Cell proliferation assay
The cell proliferation was assessed by CCK-8 (Beyotime, China) following the manufacturer’s instructions. Briefly, cells were seeded on a 96-well plate in triplicates at a density of 1000 cells/200 μl. The medium in each well was replaced with 100 μl of fresh medium supplemented with 10 μl of CCK-8 reagent. At indicated time points and incubated at 37 °C for 3–4 h. The optical density (absorbance) at a wavelength of 450 nm was determined by a multiplate reader (Bio Tek, Vermont, USA).
Cells were plated into 12-well plates at a density of 1000 cells. After incubation for approximately 7 days, the cells in 12-well plates were fixed with 4% paraformaldehyde (PFA) and stained with 1% crystal violet solution to visualize their colony-forming ability.
For apoptosis assays, FITC Annexin V apoptosis Detection Kit (BD Biosciences, San Diego, CA) was used following the manufacturer’s specifications. Briefly, the hepatoblastoma cells (2 × 105 cells/plate) in 6-well plates cells were incubated for 48 h. Cells were collected by mild trypsinization washed twice with cold PBS, stained with FITC-Annexin V and propidium iodide (PI) on ice for 5 mins, and subjected to flow cytometric analysis using a BD LSRFortessa analyzer (BD Biosciences).
Quantification of mRNA methylation with m6A-IP and RT-qPCR
m6A-IP enrichment followed by RT-qPCR to quantify the changes in m6A methylation of the target gene was performed using Magna MeRIP m6A Kit (Millipore, MA) following the manufacturer’s instructions. 5 μg of fragmented mRNA extracted from HepG2 or Huh6 stable cells was incubated with 5 μg m6A antibody (202,003, Synaptic Systems) or mouse IgG (CS200621, Millipore)-conjugated beads in 500 μl 1 × IP buffer for 4 h at 4 °C. Methylated RNA was eluted by free m6A from the beads and purified with RNeasy Mini kit (217,004, Qiagen). One tenth of the fragmented RNA was saved as an input control for standardization. The relevant enrichment of m6A of CTNNB1 in each sample was analyzed by RT-qPCR.
LC-MS/MS quantification of m6A in total RNA
Total RNA was isolated using TRIzol reagent (Invitrogen, USA) following to the manufacturer’s instruction. 1 μg of total RNA was digested by 4 μl nuclease P1 (Sigma, USA) in 40 μl buffer solution (10 mM Tris-HCl pH 7.0, 100 mM NaCl, 2.5 mM ZnCl2) at 37 °C for 12 h, followed by incubating with 1 μl alkaline phosphatase (Sigma, USA) at 37 °C for 2 h. RNA solution was diluted to 100 μl and injected into LC-MS/MS. The nucleosides were separated by reverse phase high-performance liquid chromatography on an Agilent C18 column, coupled with mass spectrometry detection using AB SCIEX QTRAP 5500. The m6A levels were calculated as the ratio of m6A to A based on the calibrated concentrations according to the standard curve obtained from pure nucleoside standards running with the same batch of samples.
Oligonucleotide transfection and lentivirus transduction
Synthetic siRNA oligonucleotides specific for regions in the METTL3, WTAP, FTO, YTHDF2, CTNNB1 mRNA were designed and synthesized by GenePharma (Shanghai, China). The sequences that gave successful knockdown were generalized in Additional file 5: Table S3. Cells were transfected with the oligonucleotides using Lipofectamine 20,000 (Invitrogen) as per the manufacturer’s instructions. To construct METTL3 stable knockdown cell lines the shRNA against METTL3 and the control shRNA obtained from General Biosystems (Anhui, China) were used. Lentivirus were packaged in HEK293T cells through co-transfecting every shRNA construct with packing vectors (PsPAX2, pMD2.G) into HEK293T cells, the lentivirus particles were harvested at 24 h, 48 h and 72 h and directly infected HepG2 and Huh6 cells under polybrene (Santa Cruz Biotechnology) for 12–24 h. Then, positively transfected HepG2 and Huh6 cells were selected for 7–10 days using 2 μg/ml puromycin (Invitrogen).
Animal studies were approved by the Institutional Animal Care and Use Committee of Shanghai Children’s Medical Center, Shanghai, China. Four-week-old male nude mice (n = 7) purchased from Shanghai Super-B&K Laboratory Animal Crop (Shanghai, China) were used. 5 × 106 cells were subcutaneously injected into the left or right flank of each mouse. Tumor volume was measured every 2 days and calculated using the following formula: Volume (mm3) = Length (mm) × Width2 (mm2) /2. Mice were sacrificed after 26 days and tumor weighed was noted.
RNA stability assays
To detect the lifetime of CTNNB1, cells were incubated with actinomycin (5 mg/ml) to terminate the transcription. Samples were collected at 0, 30, 60, 120 min post-termination. The total RNA was extracted and the remaining CTNNB1 was determined by RT-qPCR.
All the data were presented as mean ± standard deviation. Differences between two groups or multiple groups were analyzed by Student’s t-test and ANOVA, respectively. Pearson’s X2 test was used to examine the relationships between METTL3 expression levels and clinicopathological characteristics. Kaplan-Meier analysis and log-rank test were used to evaluate the differences in patient survival. The receiver operating characteristic (ROC) curves were generated to evaluate the diagnostic value. Statistical analyses were performed using the SPSS software package ver20.0 (SPSS, Inc., Chicago, IL, USA) and GraphPad Prism 7.0 (GraphPad Software, La Jolla, CA, USA). Data were considered statistically significant as follows: *p-value< 0.05, **p-value< 0.01, ***p-value< 0.001 and ****p-value< 0.0001.
The abnormal m6A modification in HB and the functional roles of METTL3, WTAP, FTO, YTHDF2 in HB cells
To explore the potential role of m6A modification in HB, we first examined m6A levels in the tumor and normal tissues. The elevation in the m6A mRNA level was identified in tumor tissues (Fig. 1a) and HB cells (Additional file 1: Figure S1A). Additionally, we also evaluated the expression of m6A writers, erasers and readers in tumor and adjacent normal hepatic tissues by quantitative PCR coupled with reverse transcription (RT-qPCR) and Western blotting. We found that most of the tumors exhibit significantly up-regulated METTL3, WTAP, FTO and YTHDF2 levels when compared against adjacent normal tissues (Fig. 1b and c). In contrast, no significant difference was observed in the expression of METTL14, KIAA1429 and ALKBH5 between the tumor and normal tissues (Additional file 1: Figure S1B). Overall, these results confirmed that m6A modification was indeed up-regulated in tumor tissue.
To investigate the functional roles of METTL3, WTAP, FTO, YTHDF2 in HB. We established METTL3, WTAP, FTO, YTHDF2 knockdown HepG2 cells with two independent siRNA sequences. RT-qPCR and Western blotting were used to investigate the successful knockdown of METTL3, WTAP, FTO, YTHDF2 in HepG2 cells (Fig. 1d). The knockdown of METTL3, WTAP, FTO, YTHDF2 genes remarkably suppressed HB cell proliferation (Fig. 1e) and inhibited their colony-forming ability (Fig. 1f). Furthermore, annexin V-PI double staining assay revealed that METTL3, WTAP, FTO, YTHDF2 knockdown can also induce apoptosis (Fig. 1g). Considering the catalytic function of these genes and the up-regulated level of m6A modification. We finally choose METTL3 as the candidate molecule as a marker for the aberrant m6A modification in HB. It may act as an oncogene promoting HB proliferation and can also serve as a crucial active component of the m6A methyltransferase having a positive correlation with m6A levels.
METTL3 as a potential diagnostic and prognostic biomarker for HB patients
As m6A modification exhibits an increasing tend with HB (Fig. 1a), a possible link between METTL3 and HB development was evaluated. Immunohistochemistry on 25 pairs of HB patient samples was performed. Based upon the grading criteria (Fig. 2a) and the Pearson’s X2 tests (Fig. 2b), findings revealed an increase in METTL3 levels in HB (Fig. 2c). There were consistent with our previous observations at the mRNA and protein levels (Fig. 1b, Fig. 1c). In the next step, patients’ clinical data was preprocessed to explore the relationship between METTL3 expression and their clinicopathological characteristics. RT-qPCR method was used to detected the METTL3 expression in tissue samples. Based upon the median ΔCT value it was divided into high and low group. Kaplan Meier analysis revealed that patients with increased METTL3 expression had frequent recurrence and had poor survival (Fig. 2d). ROC curve further confirmed the diagnostic value of METTL3 (AUC = 0.8928; P-value< 0.0001) (Fig. 2e). Taken together, these results suggest that METTL3 and could serve as a promising biomarker for the diagnosis and prognosis of HB patients.
METTL3 promotes HB tumor growth in vivo
Further to verify the oncogenic potential of METTL3 in relation with HB, we performed a subcutaneous implantation experiment in nude mice to test the effect of METTL3 knockdown in HB tumorigenicity. We used the Huh6 cells infected with lentivirus expressing shRNA against METTL3 (shMETTL3–1) and its negative control (shNC). The RT-qPCR analysis confirmed that METTL3 expression was significantly downregulated in Huh6 cells infected with shMETTL3–1 compared to the negative control (Fig. 3a). From the in vivo experiment, we observed that stable knockdown of METTL3 can cause significant reduction in tumor size and weight with respect to the nontarget shRNA control (Fig. 3b and c). Further more, we confirmed decreased expression of CTNNB1 downstream transcription factor TCF4 and target gene Cyclin D1 in METTL3 knockout HB cell Huh6 and tumor-bearing tissues of mice (Fig. 3d and e). Collectively, these experiments support the oncogenic role of METTL3 in HB tumorigenesis and development.
Widespread m6A modification of HB and identifies important biological pathways with altered methylation in HB tumors
In order to obtain a transcriptome-wide m6A map in HB, we performed m6A-seq analysis of five tumor and normal tissues. All five tumors exhibited high amount of total m6A levels. m6A-seq identified 4071 peaks in 1089 genes in tumor tissues. Among the 1089 genes, 525genes were common between the normal and tumor tissues (Fig. 4a). We also analyzed a possible relationship between m6A peaks and gene number in the normal and tumor tissues (Fig. 4b). To determine if the m6A peaks had consensus sequence of RRACH (where RR Represents purine, A is m6A and H is a non-guanine base), the top 1000 most significant peaks were analyzed,we found at least two motifs among the significant peaks representing the most common consensus in normal tissue and tumor tissue (Fig. 4c). We further investigated the distribution pattern of m6A across the whole transcriptome for both normal and tumor tissues. A similar pattern of m6A distribution was observed when the RNA species were divided into stop codon, coding sequence (CDS), start codon regions of mRNAs and non-coding RNAs. Our findings suggest that the reads from m6A-IP are highly enriched around the stop codon and CDS in both normal and tumor tissues (Fig. 4d). Further to confirm the preferential locations of m6A on transcripts, we investigated the metagene profiles of m6A peaks. Consistent with the distribution of reads, m6A peaks were abundant near the CDS (38%) and stop codon (26.2%) followed by the 3’UTR (19.4%) and then start codon (11.2%) in normal tissues. In tumor, the m6A peaks were abundant near the CDS (36.7%) and stop codon (29.5%) followed by the 3’UTR (19.4%) and start codon (11.2%) (Fig. 4e).
Further to explore the relationship between m6A peaks and mRNA expression differentially expressed genes in normal and tumor tissues from the m6A-Seq data were identified. There were 1440 strongly expressed genes in normal as compared to 2868 genes in tumor. Specific genes containing high level m6A peaks in normal were identified and labeled in blue. Similarly, in tumor these were marked in red color (Fig. 4f). Furthermore, we also observed a linear relationship between expression and percentage methylation of all genes in the bin (Fig. 4g). We extended our analysis to the whole transcriptome. The findings suggest that m6A modification leads to more genes upregulation in the tumor as compared to normal (Fig. 4h). The positions of m6A peaks involved in upregulated genes were also evaluated and it was observed that a large amount of these peaks were located at the 5′ end of the corresponding genes. We separated m6A-containing genes into PeakStart and PeakStop groups and checked how the genes expression is related to the locations of m6A peaks. Findings indicate that genes in the PeakStart and PeakStop category have overall higher expression levels and correlate well with m6A peaks in tumor (Fig. 4i).
In order to determine the effect of m6A in HB, we selected strongly expressed genes with up-regulated m6A methylation and identified the KEGG pathway. The KEGG pathway analyses were performed to organize and identify differentially activated biological processes between normal and tumor tissues. The significant genes in the tumor were related with viral carcinogenesis, proteoglycans in cancer, protein processing in endoplasmic reticulum and Wnt signaling pathway and so on. As HB is a Wnt/β-catenin-driven malignancy, we hypothesized that up-regulated m6A methylation might promote tumor growth through activation of the Wnt/β-catenin signaling pathway. Therefore, Wnt/β-catenin signaling pathway was chosen for further research (Fig. 4j). Indeed, several genes such as CTNNB1, CCND1 and NKD1 involved in the Wnt/β-catenin signaling pathway showed increased m6A methylation in tumor tissues with respect to normal (Fig. 4k). Figure 4l showed the m6A abundances of CTNNB1 were markedly increased HB tumor tissues. Overall these results are supportive of the fact that m6A has tendency of demonstrating a positive correlation with gene expression in large fraction of HB transcripts and activate m6A methylation activates the Wnt/β-catenin signaling pathway.
m6A methylation regulates activation of CTNNB1
Among CTNNB1, CCND1 and NKD1 we chose CTNNB1 as our candidate gene, as CTNNB1 is a key component of Wnt/β-catenin signaling pathway, and almost 50–90% of HB has the CTNNB1 mutation. We first confirmed CTNNB1 up regulation in HB tumor tissues (Fig. 5a and b). Knockdown CTNNB1 reduced the HB cells viability (Fig. 5c, Additional file 2: Figure S2A) and induced apoptosis (Fig. 5d, Additional file 2: Figure S2B, S2C). Furthermore, we also noticed a significantly positive correlation between CTNNB1 and METTL3 expression levels (Fig. 5e). To evaluate the relationship between METTL3 and CTNNB1, we first confirmed METTL3 was high-regulated in Huh6 and HepG2 cells (Addition file 2: Figure S2D), then we used Huh6 and HepG2 cells infected with lentivirus expressing shRNA against METTL3 (shMETTL3), or its negative control (shNC). Data from RT-qPCR and Western blotting analysis confirmed that METTL3 expression was significantly down-regulated in Huh6 and HepG2 cells infected with shMETTL3 with respect to the negative control cells (Fig. 5f, Additional file 2: Figure S2E). The results also demonstrated that reduced METTL3 expression can decrease the CTNNB1 expression (Fig. 5g, Additional file 2: Figure S2F). We then evaluated if knockdown METTL3 has a negative effect on m6A level using a colorimetric quantification assay (Fig. 5h, Additional file 2: Figure S2G). Next, we evaluated CTNNB1 in the METTL3 knockdown cells by m6A immunoprecipitation (m6A-IP) followed by RT-qPCR, which confirmed the reduction of CTNNB1 m6A methylation (Fig. 5i, Additional file 2: Figure S2H). The stability of CTNNB1 expression were evaluated by METTL3 knockdown cells. The results demonstrated that reduced m6A methylation can decrease the CTNNB1 stability (Fig. 5j, Addition file 2: Figure S2I). Taken together, the findings indicate that CTNNB1 is a downstream target of METTL3 in HB.
The discovery and mapping of m6A demethylases in mammalian systems suggests that m6A methylation of mRNA is a reversible and dynamic process with regulatory functions. Through a functional interplay between m6A methyltransferase and demethylases, the dynamic m6A modification is implicated in a variety of cellular processes, such as RNA splicing, protein translation, and stem cell renewal. Recently developed, antibody-based, high-throughput sequencing technology has allowed researchers to precisely map m6A sites to obtain further insights. Although, m6A is considered as the most prevalent alteration in human mRNA. Very few studies have focused upon the role of m6A modification in HB. In the present study we observed an increased tendency of m6A modification in HB, suggesting its potentially involvement in the malignancy. Genes METTL3, WTAP, FTO and YTHDF2 were up-regulated, and they all seem to play an essential role in promoting HB growth. Among the m6A related proteins, METTL3 and WTAP were chosen for further research as they were the m6A writers and have increased tendency of m6A modification. However, WTAP lacks methyltransferase domains and shows no catalytic activity of m6A modification. Collectively, it was found that METTL3 was the main factor for the aberrant m6A modification in HB.
The m6A levels vary in different tissues and cell lines. The loss of m6A can promote self-renewal and hinder differentiation of mouse embryonic stem cells . For instance, knockdown of mettl3 had led to reduced self-renewal abilities in mouse embryonic stem cells . However, a study had demonstrated that m6A modification was decreased in hepatocellular carcinoma, especially in metastatic hepatocellular carcinoma, and METTL14 was the main factor involved in aberrant m6A modification . These studies indicate that m6A modification is indeed a complicated regulation network in different tissues.
Previous studies have shown that m6A plays a critical role in HB. Our findings also revealed transcriptome-wide m6A distributions in HB normal and tumor tissues. We found that m6A is highly conserved across normal and tumor tissues. The m6A-Seq data suggest that m6A distribution in normal and tumor mRNA is enriched around the CDS which is distinct from mammals. Also, both in normal and tumor tissues the enrichment of m6A correlates positively with the overall upregulation of mRNA expression level. Nevertheless, there were different methylated mRNA between normal and tumor tissues. Methylation levels of 243 and 200 m6A mRNAs were upregulated in tumor tissue and normal tissues. Which provide a starting roadmap for uncovering m6A functions in HB proliferation.
The Wnt/β-catenin pathway is frequently mutated and overactive in solid malignancies promoting tumor development. In HB, driving and most recurrently mutated proto-oncogene with 50–90% frequency is CTNNB1. In this study, we discovered that m6A mRNA methylation regulates the Wnt/β-catenin pathway significantly. In HB, m6A methylation normally induces Wnt/β-catenin pathway activity by promoting the expression of CTNNB1. Using cultured HB cells, we revealed that reduced m6A methylation leads to a decrease in CTNNB1 expression and stability. We also used m6A experiment to confirm that decreased m6A methylation further decreased CTNNB1 m6A modification enrichment. Our m6A-Seq results from HB tumors and normal tissues along with our mechanistic studies in HB cells indicate that regulation of CTNNB1 m6A modification as an important mediator of these changes for cell proliferation. In our future studies we will be exploring the relationship between the m6A modification and CTNNB1 mutation. Also the question if m6A modification can induce mutated CTNNB1 to accumulate in the nucleus and promote the binding the TCF/LEF transcription factor will be elucidated further.
In summary, our study delineates a METTL3-m6A module that regulates the process of CTNNB1 with respect to its functional implications in HB proliferation. Our results suggest that m6A modification plays a critical role in tumor biology. Also, the m6A modification of CTNNB1 might be a general growth mechanism affecting a range of other biological processes, that warrants further exploration.
Availability of data and materials
All data generated or analyzed during the current study are included in this published article (and its supplementary information files).
AlkB homoogue 5
Dulbecco’s Modified Eagle’s Medium
Fetal Bovin Serum
Fat mass and obesity-associatd protein
Minimum Essential Medium
Methylated RNA immunoprecipitation
Receiver operating characteristic
Reverse transcription quantitative real-time PCR
Sodium dodecyl sulfate
Sodium dodecyl sulfate-polyacrylamide gel electrophoresis
Suppressor Of Cytokine Signaling 2
Wilms tumor 1-assocated protein
Finger CCCH domain-containing protein 13
Trobaugh-Lotrario AD, Katzenstein HM. Chemotherapeutic approaches for newly diagnosed hepatoblastoma: past, present, and future strategies. Pediatr Blood Cancer. 2012;59:809–12.
Sharma D, Subbarao G, Saxena R. Hepatoblastoma. Semin Diagn Pathol. 2017;34:192–200.
Meyers RL, Tiao G, de Ville de Goyet J, Superina R, Aronson DC. Hepatoblastoma state of the art: pre-treatment extent of disease, surgical resection guidelines and the role of liver transplantation. Curr Opin Pediatr. 2014;26:29–36.
Venkatramani R, Stein JE, Sapra A, Genyk Y, Jhaveri V, Malogolowkin M, Mascarenhas L. Effect of neoadjuvant chemotherapy on resectability of stage III and IV hepatoblastoma. Br J Surg. 2015;102:108–13.
Yuan XJ, Wang HM, Jiang H, Tang MJ, Li ZL, Zou X, Fang YJ, Pan C, Tou JF, Zhang KR, et al. Multidisciplinary effort in treating children with hepatoblastoma in China. Cancer Lett. 2016;375:39–46.
Russell JO, Monga SP. Wnt/beta-catenin signaling in liver development, homeostasis, and pathobiology. Annu Rev Pathol. 2018;13:351–78.
Crippa S, Ancey PB, Vazquez J, Angelino P, Rougemont AL, Guettier C, Zoete V, Delorenzi M, Michielin O, Meylan E. Mutant CTNNB1 and histological heterogeneity define metabolic subtypes of hepatoblastoma. EMBO Mol Med. 2017;9:1589–604.
Lopez-Terrada D, Gunaratne PH, Adesina AM, Pulliam J, Hoang DM, Nguyen Y, Mistretta TA, Margolin J, Finegold MJ. Histologic subtypes of hepatoblastoma are characterized by differential canonical Wnt and notch pathway activation in DLK+ precursors. Hum Pathol. 2009;40:783–94.
Ashihara E, Takada T, Maekawa T. Targeting the canonical Wnt/beta-catenin pathway in hematological malignancies. Cancer Sci. 2015;106:665–71.
Taniguchi K, Roberts LR, Aderca IN, Dong X, Qian C, Murphy LM, Nagorney DM, Burgart LJ, Roche PC, Smith DI, et al. Mutational spectrum of beta-catenin, AXIN1, and AXIN2 in hepatocellular carcinomas and hepatoblastomas. Oncogene. 2002;21:4863–71.
Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, Cesarkas K, Jacob-Hirsch J, Amariglio N, Kupiec M, et al. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 2012;485:201–6.
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. 2014;505:117–20.
Xiao W, Adhikari S, Dahal U, Chen YS, Hao YJ, Sun BF, Sun HY, Li A, Ping XL, Lai WY, et al. Nuclear m (6) a reader YTHDC1 regulates mRNA splicing. Mol Cell. 2016;61:507–19.
Ma JZ, Yang F, Zhou CC, Liu F, Yuan JH, Wang F, Wang TT, Xu QG, Zhou WP, Sun SH. METTL14 suppresses the metastatic potential of hepatocellular carcinoma by modulating N (6) -methyladenosine-dependent primary MicroRNA processing. Hepatology. 2017;65:529–43.
Pan Y, Ma P, Liu Y, Li W, Shu Y. Multiple functions of m (6) a RNA methylation in cancer. J Hematol Oncol. 2018;11:48.
Tong J, Flavell RA, Li HB. RNA m (6) a modification and its function in diseases. Front Med. 2018;12:481–9
Yang Y, Hsu PJ, Chen YS, Yang YG. Dynamic transcriptomic m (6) a decoration: writers, erasers, readers and functions in RNA metabolism. Cell Res. 2018;28:616–24.
Wen J, Lv R, Ma H, Shen H, He C, Wang J, Jiao F, Liu H, Yang P, Tan L, et al. Zc3h13 regulates nuclear RNA m (6) a methylation and mouse embryonic stem cell self-renewal. Mol Cell. 2018;69:1028–38 e1026.
Liu J, Yue Y, Han D, Wang X, Fu Y, Zhang L, Jia G, Yu M, Lu Z, Deng X, et al. A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat Chem Biol. 2014;10:93–5.
Sorci M, Ianniello Z, Cruciani S, Larivera S, Ginistrelli LC, Capuano E, Marchioni M, Fazi F, Fatica A. METTL3 regulates WTAP protein homeostasis. Cell Death Dis. 2018;9:796.
Wang P, Doxtader KA, Nam Y. Structural basis for cooperative function of Mettl3 and Mettl14 Methyltransferases. Mol Cell. 2016;63:306–17.
Wu Y, Xie L, Wang M, Xiong Q, Guo Y, Liang Y, Li J, Sheng R, Deng P, Wang Y, et al. Mettl3-mediated m (6) a RNA methylation regulates the fate of bone marrow mesenchymal stem cells and osteoporosis. Nat Commun. 2018;9:4772.
Ben-Haim MS, Moshitch-Moshkovitz S, Rechavi G. FTO: linking m6A demethylation to adipogenesis. Cell Res. 2015;25:3–4.
Xu C, Liu K, Tempel W, Demetriades M, Aik W, Schofield CJ, Min J. Structures of human ALKBH5 demethylase reveal a unique binding mode for specific single-stranded N6-methyladenosine RNA demethylation. J Biol Chem. 2014;289:17299–311.
Zheng G, Dahl JA, Niu Y, Fedorcsak P, Huang CM, Li CJ, Vagbo CB, Shi Y, Wang WL, Song SH, et al. ALKBH5 is a mammalian RNA demethylase that impacts RNA metabolism and mouse fertility. Mol Cell. 2013;49:18–29.
Meyer KD, Patil DP, Zhou J, Zinoviev A, Skabkin MA, Elemento O, Pestova TV, Qian SB, Jaffrey SR. 5′ UTR m (6) a promotes cap-independent translation. Cell. 2015;163:999–1010.
Zhang S, Zhao BS, Zhou A, Lin K, Zheng S, Lu Z, Chen Y, Sulman EP, Xie K, Bogler O, et al. m (6) a Demethylase ALKBH5 maintains Tumorigenicity of Glioblastoma stem-like cells by sustaining FOXM1 expression and cell proliferation program. Cancer Cell. 2017;31:591–606 e596.
Li Z, Weng H, Su R, Weng X, Zuo Z, Li C, Huang H, Nachtergaele S, Dong L, Hu C, et al. FTO plays an oncogenic role in acute myeloid leukemia as a N (6)-Methyladenosine RNA Demethylase. Cancer Cell. 2017;31:127–41.
Zhou J, Wan J, Gao X, Zhang X, Jaffrey SR, Qian SB. Dynamic m (6) a mRNA methylation directs translational control of heat shock response. Nature. 2015;526:591–4.
Chen M, Wei L, Law CT, Tsang FH, Shen J, Cheng CL, Tsang LH, Ho DW, Chiu DK, Lee JM, et al. RNA N6-methyladenosine methyltransferase-like 3 promotes liver cancer progression through YTHDF2-dependent posttranscriptional silencing of SOCS2. Hepatology. 2018;67:2254–70.
Batista PJ, Molinie B, Wang J, Qu K, Zhang J, Li L, Bouley DM, Lujan E, Haddad B, Daneshvar K, et al. m (6) a RNA modification controls cell fate transition in mammalian embryonic stem cells. Cell Stem Cell. 2014;15:707–19.
Wang Y, Li Y, Toth JI, Petroski MD, Zhang Z, Zhao JC. N6-methyladenosine modification destabilizes developmental regulators in embryonic stem cells. Nat Cell Biol. 2014;16:191–8.
This work was supported by Program of Shanghai Academic/Technology Research Leader (18XD1402600),Municipal Human Resources Development Program for Outstanding Leaders in Medical Disciplines in Shanghai (2017BR036),Shanghai Municipal Education Commission -Gaofeng Clinical Medicine Grant Support (20171926), Science and Technology Development.
Fund of Pudong New Area of Shanghai (PKJ2017-Y03),National Natural Science Foundation of China (81572330),National Natural Science Foundation of China (81871727),National Natural Science Foundation of China (81472624), National Natural Science Foundation of China (81772941). Key Program of the National Natural Science Foundation (81930066). Shanghai anti-cancer association eagle plan (SACA-CY1B05). Shanghai municipal health and family planning commission (20174Y0228). National Natural Science Foundation of China (81702741), Natural Science Foundation of Science and Technology Commission of Shanghai (17ZR1421700).
Ethics approval and consent to participate
All the collection of specimens and animal handling in this work was reviewed and approved by the Medical Ethics Committee of the Shanghai Children’s Medical Center.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Li Liu and Jing Wang contributed equally to this work.
Additional file 1: Figure S1. The m6A mRNA level in HB cells and the expression of METTL14, KIAA1429 and ALKBH5 in HB tumor tissues. (A) LC-MS/MS quantification of the m6A/A ratio in total RNA isolated from hepatocytes and HB cells (one-way analysis of variance, Dunnett’s test). (B) The expression of METTL14, KIAA1429 and ALKBH5 were quantified by RT-qPCR in 10 paired HB tumor and normal tissues (paired samples t-test). ns, no significant difference. **p-value< 0.01,***p-value< 0.001.
Additional file 2: Figure S2. Reduced m6A methylation decreases CTNNB1 expression and stability in HepG2 cells. (A) CCK8 assay was used to evaluate the viability of CTNNB1 knockdown in HepG2 cells measured (one-way analysis of variance, Dunnett’s test). (B) The apoptosis rate in CTNNB1 knockdown in HepG2 cells were measured by FACS assays. (C) Histogram of flow cytometry assays from 3 independent experiments (one-way analysis of variance, Dunnett’s test). (D) METTL3 expression in hepatocytes and HB cells was evaluated by RT-qPCR (one-way analysis of variance, Dunnett’s test). (E) METTL3 stable knockdown in HepG2 cells were generated by lentiviral-based shRNA expression. METTL3 knockdown efficiency was confirmed at both the mRNA and protein levels by RT-qPCR and western blotting (one-way analysis of variance, Dunnett’s test). (F) RT-qPCR was used to test the expression of CTNNB1 upon METTL3 knockdown in HepG2 cells (one-way analysis of variance, Dunnett’s test). (G) Relative m6A level in knockdown METTL3 in HepG2 cells (independent-samples t-test). (H) m6A-IP combined with RT-qPCR was used to quantify the relative m6A modified level of CTNNB1 upon METTL3 depletion in HepG2 cells (independent-samples t test). (I) Lifespans of CTNNB1 expression in cells transfected with the shMETTL3 in HepG2 cells. Relative mRNA levels were quantified by RT-qPCR. *p-value< 0.05,**p-value< 0.01,***p-value< 0.001, ****p-value<0.0001.
About this article
Cite this article
Liu, L., Wang, J., Sun, G. et al. m6A mRNA methylation regulates CTNNB1 to promote the proliferation of hepatoblastoma. Mol Cancer 18, 188 (2019). https://doi.org/10.1186/s12943-019-1119-7
- RNA m6A methylation
- Wnt/β-catenin pathway