YTHDC1 is downregulated by the YY1/HDAC2 complex and controls the sensitivity of ccRCC to sunitinib by targeting the ANXA1-MAPK pathway
Journal of Experimental & Clinical Cancer Research volume 41, Article number: 250 (2022)
Tyrosine kinase inhibitors (TKIs) such as sunitinib are multitarget antiangiogenic agents in clear cell renal cell carcinoma (ccRCC). They are widely used in the treatment of advanced/metastatic renal cancer. However, resistance to TKIs is common in the clinic, particularly after long-term treatment. YTHDC1 is the main nuclear reader protein that binds with m6A to regulate the splicing, export and stability of mRNA. However, the specific role and corresponding mechanism of YTHDC1 in renal cancer cells are still unclear.
The Cancer Genome Atlas (TCGA) dataset was used to study the expression of YTHDC1 in ccRCC. Cell counting kit-8 (CCK-8), wound healing, Transwell and xenograft assays were applied to explore the biological function of YTHDC1 in ccRCC. Western blot, quantitative real time PCR (RT‒qPCR), RNA immunoprecipitation PCR (RIP-qPCR), methylated RIP-qPCR (MeRIP-qPCR) and RNA sequencing (RNA-seq) analyses were applied to study the YY1/HDAC2/YTHDC1/ANXA1 axis in renal cancer cells. The CCK-8 assay and xenograft assay were used to study the role of YTHDC1 in determining the sensitivity of ccRCC to sunitinib.
Our results demonstrated that YTHDC1 is downregulated in ccRCC tissues compared with normal tissues. Low expression of YTHDC1 is associated with a poor prognosis in patients with ccRCC. Subsequently, we showed that YTHDC1 inhibits the progression of renal cancer cells via downregulation of the ANXA1/MAPK pathways. Moreover, we also showed that the YTHDC1/ANXA1 axis modulates the sensitivity of tyrosine kinase inhibitors. We then revealed that HDAC2 inhibitors resensitize ccRCC to tyrosine kinase inhibitors through the YY1/HDAC2 complex. We have identified a novel YY1/HDAC2/YTHDC1/ANXA1 axis modulating the progression and chemosensitivity of ccRCC.
We identified a novel YY1/HDAC2/YTHDC1/ANXA1 axis modulating the progression and chemosensitivity of ccRCC.
Clear cell renal cell carcinoma (ccRCC) is the main type of renal carcinoma . Although surgical treatment is the mainstay of treatment for renal cancer, patients with unresectable advanced renal cancer or recurrent renal cancer rely on systemic treatment, typically antiangiogenic or immunotherapeutic treatment, to prolong their lives [2, 3]. Tyrosine kinase inhibitors (TKIs) such as sunitinib are multitarget antiangiogenic agents in renal carcinoma cells . TKIs are widely used in the treatment of advanced/metastatic renal cancer . However, resistance to TKIs is common in the clinic, particularly after long-term treatment . The mechanism of drug resistance may be related to proangiogenic signaling pathways, changes in the tumor microenvironment, increased tumor invasion and metastasis, microRNA-mediated drug resistance or the activation of other signaling pathways . However, the factors associated with TKI resistance in ccRCC have not been fully clarified. Thus, further study of the mechanism of TKI resistance would be helpful for the development of new antitumor strategies for ccRCC.
N6-Methyladenosine (m6A) is the most abundant modification of internal eukaryotic mRNA . m6A is associated with mRNA metabolism, folding, processing, export, translation and stability . It has been indicated that m6A modification is closely related to various types of tumors such as lung cancer, acute myeloid leukemia (AML), and hepatocellular carcinoma (HCC) . Renal cancer, particularly clear cell renal cell carcinoma (ccRCC), is correlated with m6A modification . Zhou et al. found that changes in m6A-regulating factors were correlated with pathological stage in ccRCC . There are three groups of molecules that regulate m6A modification: m6A “writers” (METTL14, METTL3, and WTAP), m6A “erasers” (ALKBH5 and FTO) and m6A “readers” (YTHDC1, HNRNPC, HNRNPA2B, YTHDF1, YTHDF2 and YTHDF3) . m6A readers are signal transducers that regulate mRNA modification and are dysregulated in renal cell carcinoma [11,12,13]. Felix et al. analyzed the expression of six m6A readers in ccRCC and found that only the mRNA and protein levels of YTHDC1 were consistently downregulated in ccRCC . Of note, YTHDC1 is the main nuclear reader protein that binds with m6A to regulate the splicing, export and stability of mRNA [13,14,15,16]. However, the specific role and corresponding mechanism of YTHDC1 in renal cancer cells are still unclear.
Consistent with previous findings, we demonstrated that YTHDC1 is downregulated in ccRCC and that low expression of YTHDC1 is correlated with a poor prognosis. We then showed that YTHDC1 participates in inhibiting the progression of ccRCC. Interestingly, we found that YTHDC1 represses the expression of ANXA1, which was previously reported to modulate TKI resistance in ccRCC . Moreover, we revealed that the HDAC2/YY1 complex regulates the expression of YTHDC1 in ccRCC and that HDAC2 inhibitors could sensitize ccRCC cells to TKI treatment through YTHDC1. Thus, we identified a novel HDAC2/YY1/YTHDC1/ANXA1 axis in ccRCC.
Materials and methods
Cell lines and cell culture
The renal cell carcinoma cell lines 786-O (#CL-0010, Procell Life Science & Technology, Wuhan, China) and A498 (#CL-0254, Procell Life Science & Technology, Wuhan, China) were identified by short tandem repeat (STR) profiling using Procell Life Science & Technology. The 786-O cells were cultured in RPMI-1640 medium (Gibco, USA), and the A498 cells were cultured in MEM (PM150410, Procell Life Science & Technology) supplemented with 10% fetal bovine serum (FBS) (AC03L055, Shanghai Life-iLab Biotech, China) and incubated at 37 °C in 5% CO2. Sunitinib-resistant 786-O cells (786-O R) were generated as previously reported .
The chemicals used were as follows: The HDAC2 inhibitors included santacruzamate A (CAY10683) (# S7595, Selleck, China), LY3214996 (#S8534, Selleck, China), GSK1120212 (#S2673, Selleck, China), and actinomycin D (#S8964, Selleck, China). The siRNA was obtained from RiboBio (Guangzhou, China). The shRNAs were purchased from GeneCopoeia (USA). The sequences of the siRNA and shRNA are provided in Table S1.
Western blot analysis and antibodies
The western blot analysis procedure was reported previously. In brief, the cells were lysed and the proteins extracted with RIPA buffer (#P0013, Beyotime, China). The protein content was quantified by the BCA method. Proteins were boiled and subjected to sodium dodecyl sulfate‒polyacrylamide gel electrophoresis (SDS‒PAGE) separation. The proteins were transferred onto 0.45 μm polyvinylidene fluoride membranes (Millipore, USA) and incubated with the corresponding primary antibodies and secondary antibodies. Protein signals were visualized using ECL detection reagent (Thermo Fisher Scientific, USA) and ChemiDoc XRS (Bio-Rad Laboratories, USA). The primary antibodies used were as follows: YTHDC1 (#29441-1-AP, Proteintech, 1:500 dilution), ANXA1 (#55018-1-AP, Proteintech, 1:2000 dilution), YY1 (#66281-1-Ig, Proteintech, 1:1000 dilution), HDAC2 (#12922-3-AP, Proteintech, 1:1000 dilution), and beta actin (#66009-1-Ig, Proteintech, 1:5000 dilution).
Immunohistochemistry was performed with primary antibodies against YTHDC1 (#29441-1-AP, Proteintech, 1:1000, ANXA1 (21990-1-AP, Proteintech; 1:10000 dilution) and TMA slides (#KD1921, Avilabio, China). The extent of binding was evaluated based on the proportion of positively stained tumor cells and the staining intensity score. The staining intensity score was as follows: 0 (no staining), 1 (weak staining = light yellow), 2 (moderate staining = yellow brown), and 3 (strong staining = brown). A total score was obtained by multiplying the intensity score by the proportion of cells stained.
Quantitative real-time PCR (RT‒qPCR) and chromatin immunoprecipitation (ChIP)-qPCR analysis
The details of RT‒qPCR were described previously . In brief, TRIzol reagent (Thermo Fisher Scientific, USA) was used to extract total RNA from cells. A reverse transcription kit and PCR kit (#RR037A PrimeScript™ RT reagent Kit and #RR430A TB Green™ Fast qPCR Mix, respectively, Takara Bio Inc. Shigo, Japan) were used to perform the RT‒qPCR assay. The primer sequences for RT‒PCR are provided in Table S2. A chromatin extraction kit (Abcam, ab117152, USA) and the ChIP Kit Magnetic-One Step (Abcam, ab156907, USA) were used to perform ChIP‒qPCR. The primer sequences for ChIP‒PCR are provided in Table S3.
Cell counting Kit-8 (CCK-8) assay
The 786° or A498 cells were seeded in 96-well plates; approximately 104 cells were seeded per well. After culturing for 24 h at 37 °C in 5% CO2, the cells were divided into several groups with different treatments; each group had at least 3 repetitions. Ten microliters of CCK-8 reagent (#C0037, Beyotime, China) was added to each well and incubated for 1 h under the above conditions. The absorbance at 450 nm was measured by a microplate reader.
Nude mouse xenograft assay
Because the gender of the mice had no effect on the results of the study, we chose half males and half females for the experiment. BALB/C-nu/nu mice (Hunan SJA Laboratory Animal Company), 6 weeks old and 22–24 g, were individually housed in the animal center of the Second Xiangya Hospital under 60 ± 3% humidity at 22 ± 0.5 °C; the light/dark cycle was controlled automatically. Animals had free access to food and water. The renal carcinoma cells were injected subcutaneously on the left side of the back of the mouse (1 × 107 cells per mouse). The length and width of the tumor were measured with a Vernier caliper every 2 d, and the tumor volume was calculated using the formula (L × W2)/2. Mice were sacrificed at the appropriate time, and the tumors were collected for further study. All animal experiments were reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of the Second Xiangya Hospital, Central South University (animal license number 2,021,095). The animal experiment complied with the National Institutes of Health Guide for the Care and Use of Laboratory Animals (NIH Publications No. 8023, revised 1978). Because a minimum of four mice per group is necessary to compare differences between two groups, in the nude mouse xenograft assays, six mice per group were originally established for both the 786-O and A498 cell lines. However, in the A498 group, two mice died unnaturally during the course of the experiment. Thus, in the in vivo experimental design, the 786-O cell line comparison was 6 vs. 6, and the A498 cell line comparison was 4 vs. 4.
RNA immunoprecipitation (RIP) and Methylated (Me) RIP-qPCR
Whole cells were lysed with 1 × RIPA buffer (Beyotime Biotech, China), and 10% of the lysates were collected as input. The primary antibody or IgG and protein A + G beads (#P2029, Beyotime, China) were then added to the rest of the cell lysate and incubated on a shaker at 4 ℃ for 12 h. The beads were washed with 1 × RIPA buffer, resurrected with protein kinase K buffer and incubated with shaken RNAiso Plus (TaKaRa, 9109) at 55 ℃ for 30 min. RT‒qPCR was performed using the reverse transcription kit and the polymerase chain reaction kit according to the manufacturer’s instructions. The primer sequences are provided in Table S4. The Magna M6A MeRIP kit (# A-17-10499, A&D Technology, Beijing, China) was used to perform MeRIP-qPCR according to the manufacturer’s protocols. The primer sequences are provided in Table S5.
The experimental data came from three independent experiments and are presented as the mean ± standard error of the mean (mean ± SEM). GraphPad Prism 5 software was used to calculate the P value using an unpaired two-sided Student’s t test to compare values between two groups or a one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparisons post hoc test to compare values between more than two groups. The sample size (n) for each statistical analysis is provided in the figure legends. Differences were considered statistically significant when the P values were less than 0.05. In all cases, the significance of differences was indicated as follows: *P < 0.05; **P < 0.01; ***P < 0.001; and not significant (ns), P > 0.05.
Downregulation of YTHDC1 promotes the progression of ccRCC in vivo and in vitro
We first studied the clinical characteristics of YTHDC1 in malignant tumors. Interestingly, we showed that YTHDC1 was profoundly downregulated in ccRCC with a hazard ratio (HR) of less than 1.0 and P values less than 0.05 by analyzing the ENCORI dataset (https://starbase.sysu.edu.cn/panCancer.php) and the Timer 2.0 dataset (http://timer.cistrome.org/) (Fig. 1A and B). Analysis of the TCGA dataset also demonstrated that low YTHDC1 expression was associated with a poor prognosis in patients with ccRCC (Fig. 1C). These data suggest that YTHDC1 is a potential prognostic predictor in ccRCC. The CancerSEA dataset (http://biocc.hrbmu.edu.cn/CancerSEA/goSearch) was then assessed to predict the cancer-associated role of YTHDC1 in renal cancer cells. We found that YTHDC1 expression was negatively correlated with the proliferation, invasion and metastasis of renal cell carcinoma cells (Fig. 1D). To verify this, YTHDC1 was knocked down by transfection with two independent shRNAs in A498 and 786-O cells (Fig. 1E, F). We showed that YTHDC1 inhibition promoted renal cancer cell proliferation, migration, and invasion in vitro (Fig. 1G and I). Conversely, ectopic overexpression of YTHDC1 impeded the proliferation, migration, and invasion of both A498 and 786-O cells (Fig. 1J-N). In addition, the nude mouse xenograft assay also demonstrated that knockdown of YTHDC1 increased the tumor growth of renal cancer cells (Fig. 2A-J). Thus, our data suggest that YTHDC1 is abnormally downregulated in renal cancer tissues, downregulated YTHDC1 is associated with unfavorable prognosis, and YTHDC1 inhibits the progression of renal cancer in cells and in mice.
YTHDC1 contributes to the inactivation of the ERK/MAPK signaling pathway in ccRCC cells
Next, we explored the underlying mechanism by which YTHDC1 inhibits the progression of renal cancer cells. Transcriptome sequencing was performed in A498 cells after YTHDC1 knockdown (Supplementary Fig. 1A). Reactome enrichment analysis demonstrated that YTHDC1 was negatively correlated with EGFR and FGFR signaling pathways in renal cancer cells, which are reported to be the major driver pathways for the tumorigenesis of renal cancer [18,19,20] (Fig. 3A). Of note, a Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that YTHDC1 inhibition activated the EGFR tyrosine kinase inhibitor resistance pathway in cancer (Fig. 3B). A Gene Set Enrichment Analysis (GSEA) of the TCGA-KIRC dataset demonstrated that low YTHDC1 expression was positively correlated with renal cell carcinoma and the mitogen-activated protein kinase (MAPK) signaling pathway (Fig. 3C). KEGG enrichment analysis indicated that YTHDC1 inhibition resulted in the activation of the MAPK signaling pathway and the tumor necrosis factor (TNF) signaling pathway in cells (Fig. 3D). It is well known that the extracellular signal-regulated kinases (ERK)/MAPK signaling pathway is responsible for the progression of clear cell renal cell carcinoma [21, 22]. Thus, we sought to determine whether YTHDC1 affects the activation of the MAPK signaling pathway in renal cancer cells. Notably, knockdown of YTHDC1 increased the phosphorylation of ERK, which is a marker of activation of the MAPK signaling pathway, in 786-O and A498 cells (Fig. 3E). In contrast, overexpression of YTHDC1 decreased the phosphorylation of ERK in renal cancer cells (Fig. 3 F). We then demonstrated that treatment with ERK or MEK inhibitors (LY3214996 or GSK1120212, respectively) attenuated the effect of YTHDC1 on suppressing renal cancer cell proliferation (Fig. 3G J). Together, these data indicate that YTHDC1 inactivates the MAPK signaling pathway to inhibit the proliferation of renal cancer cells.
YTHDC1 represses ANXA1 expression in ccRCC cells
As an m6A reader, YTHDC1 has been found to regulate mRNA splicing , nuclear export , and stability . We rechecked the RNA-seq data of YTHDC1 and found that Annexin-A1 (ANXA1) was upregulated after knockdown of YTHDC1 in cancer cells (Supplementary Fig. 1 A and B). Although ANXA1 was not associated with prognosis in patients with renal clear cell carcinoma (Supplementary Fig. 1 C and D), ANXA1 is abnormally upregulated in renal cancer specimens, and inhibition of ANXA1 decreases the proliferation and invasion of renal cancer cells . Interestingly, ANXA1 has been suggested to regulate EGFR activity in malignancy . Moreover, ANXA1 contributes to the ERK-NFκB activation loop in breast cancer cells . Because YTHDC1 has been reported to regulate the expression level of target genes by influencing their mRNA stability , we speculated whether ANXA1 was the downstream gene of YTHDC1 in renal cancer cells. First, knockdown of YTHDC1 increased the protein and mRNA levels of ANXA1 in 786-O and A498 cells (Fig. 4A and B). Meanwhile, we also showed that YTHDC1 inhibition led to the up-regulation of ANXA1 in xenografts derived from the 786-O cells (Supplementary Fig. 2 A). In contrast, overexpression of YTHDC1 decreased ANXA1 expression in renal cancer cells (Fig. 4C and D). We then applied the Encyclopedia of RNA Interactomes (ENCORI) dataset to explore the stacked peak regions of YTHDC1 on ANXA1. There were two stacked regions (chr9, 75,773,656–75,773,682; chr9, 75,778,393–75,778,412) on ANXA1 (Supplementary Fig. 1E). We assessed the m6A modification sites of ANXA1 using the RMBase v2.0 (https://rna.sysu.edu.cn/rmbase/index.php) web tool (Supplementary Fig. 1F), and we found that the chr9 75,778,393–75,778,412 region is simultaneously bound by YTHDC1 and tagged with an m6A modification (Supplementary Fig. 1E and F). A RIP assay using YTHDC1 antibodies was performed around the chr9 75,778,393–75,778,412 region. We showed that YTHDC1 binds to ANXA1 in 786-O and A498 cells (Fig. 4E). Moreover, m6A modification was also observed in this region (Fig. 4F). Given that YTHDC1 is also responsible for RNA translocation , we knocked down or overexpressed YTHDC1 but found no changes in the subcellular localization of ANXA1 in either 786-O or A498 cells (Fig. 4G H). Interestingly, we showed that YTHDC1 inhibition increased the stability of ANXA1 mRNA, but overexpression of YTHDC1 decreased the stability of ANXA1 mRNA in 786-O and A498 cells (Fig. 4I J). Moreover, IHC staining of the tissue microarray indicated that YTHDC1 was negatively correlated with ANXA1 in renal cancer (Fig. 4K and L). Furthermore, we stained the nude mouse xenograft sample with ANXA1 and YTHDC1 antibodies. We showed that ANXA1 was upregulated after YTHDC1 knockdown (Supplementary Fig. 2). Therefore, these results suggest that YTHDC1 decreases the mRNA stability of ANXA1 in renal cancer cells.
YTHDC1 inhibits activation of the MAPK signaling pathway by targeting ANXA1 in ccRCC
Because ANXA1 controls MAPK activation in cells [29, 30], we investigated whether YTHDC1 inhibits the MAPK pathway through ANXA1. We constructed renal cancer cells with stable ANXA1 knockdown induced by transfection with ANXA1-specific shRNAs following puromycin selection (Fig. 5A). YTHDC1 was then overexpressed in shControl or shANXA1 cells by transfection with Flag-tagged YTHDC1 plasmids (Fig. 5A). We found that overexpression of YTHDC1 decreased the phosphorylation of ERK in the shControl group but had no obvious effect on the phosphorylation of ERK in the shANXA1 group in renal cancer cells (Fig. 5A). Similarly, we demonstrated that the inhibitory effect on cell proliferation or invasion after overexpression of YTHDC1 was diminished in shANXA1 renal cancer cells (Fig. 5B and D). In contrast, we constructed ANXA1, YTHDC1 knockdown alone, or ANXA1 + YTHDC1 coknockdown renal cancer cells (Fig. 5E). We found that ANXA1 repression could also attenuate the effect of YTHDC1 knockdown on the change in ERK phosphorylation. Furthermore, the CCK-8 assay, colony formation assay, nude mouse xenograft assay (the tumor image is indicated in Fig. 6K) and Transwell assay showed that the change in cell proliferation ability and invasion capability after knockdown of YTHDC1 was diminished after coknockdown of ANXA1 in 786-O and A498 cells (Fig. 5F-K). Thus, our data indicate that the YTHDC1-ANXA1 axis contributes to the activation of the MAPK signaling pathway in renal cancer cells.
The YTHDC1-ANXA1 axis regulates sunitinib sensitivity in ccRCC
Tyrosine kinase inhibitors such as sunitinib are the first-line treatment for advanced renal cell carcinoma [31, 32]. Emerging resistance to sunitinib has become the major obstacle for prolonging the survival of patients with renal cell carcinoma [33, 34]. The KEGG enrichment analysis of the RNA-seq of YTHDC1 showed that YTHDC1 was involved in the EGFR tyrosine kinase inhibitor resistance pathways in A498 cells (Fig. 3B). We previously reported that ANXA1 participates in modulating the sensitivity to sunitinib in ccRCC . The literature has also documented that the aberrant activation of MEK-ERK signaling induces sunitinib resistance in renal cancer [35, 36]. It is not surprising that YTHDC1 knockdown increased the IC50 values of sunitinib in A498 and 786-O cells (Fig. 6A). Ectopic overexpression of YTHDC1 reduced the IC50 values of sunitinib in A498, 786-O and sunitinib-resistant 786-O (786-O R) cells, as reported previously  (Fig. 6B). The CCK-8 assay and Annexin V-FITC/propidium iodide staining assay demonstrated that knockdown of YTHDC1 decreased the sensitivity of renal cancer cells to sunitinib (Fig. 6C and D, Supplementary Fig. 3 A), but overexpression of YTHDC1 led to sunitinib resistance in renal cancer cells (Fig. 6E F, Supplementary Fig. 3B). Moreover, we showed that ANXA1 depletion could attenuate the change in the IC50 values of sunitinib induced by knockdown or overexpression of YTHDC1 in both A498 and 786-O cells (Fig. 6G H). Subsequent in vitro and in vivo cell proliferation studies indicated that knockdown of ANXA1 blocked the effect of YTHDC1 on modulating the sensitivity of renal cancer cells to sunitinib (Fig. 6I M). Together, our results demonstrate that YTHDC1 regulates the sensitivity of renal cancer cells to sunitinib through ANXA1.
The YY1/HDAC2 complex downregulates the expression of YTHDC1 in ccRCC
The above results showed that YTHDC1 played an important role in regulating the progression and drug resistance of renal cancer cells, which suggests that YTHDC1 might be an ideal candidate for the treatment of ccRCC. The following experiment was conducted to explore the regulatory mechanism of YTHDC1 in renal cancer cells. The PROMO dataset (http://alggen.lsi.upc.es/cgi-bin/promo_v3/promo/promoinit.cgi?dirDB=TF_8.3) and the GeneCards dataset (https://www.genecards.org/) were applied to predict the potential transcription factors of YTHDC1 (Fig. 7A and B). Among them, YY1 was found in both datasets (Fig. 7A and B). Next, we examined whether YY1 transcriptionally regulated the expression of YTHDC1 in renal cancer cells. We found that knockdown of YY1 increased the protein and mRNA levels of YTHDC1 in 786-O and A498 cells (Fig. 7C and D). We also showed that YY1 inhibition reduced the phosphorylation of ERK in renal cancer cells (Supplementary Fig. 4 A). In addition, the ChIP–qPCR analysis indicated that YY1 was enriched in the promoter region of YTHDC1, and YY1 inhibition decreased this enrichment (Fig. 7E). Moreover, we cloned a DNA sequence containing the potential binding sequence of YY1 and two mutants to generate three GV592-YTHDC1 promoter plasmids (Fig. 7F). We showed that the luciferase activity of the GV592-YTHDC1 WT promoter plasmid was higher than that of the MUT1 or MUT2 plasmid in both 786-O and A498 cells (Fig. 7G). Moreover, we also demonstrated that the luciferase activity of MUT1 was greater than that of MUT2 (Fig. 7G). Furthermore, we showed that knockdown of YY1 decreased the luciferase activity of GV592-YTHDC1 WT in 786-O cells (Fig. 7H). These data suggest that YY1 decreased the expression of YTHDC1 at the transcriptional level. Analysis of the ChIP-atlas dataset also found that histone deacetylase 2 (HDAC2) colocalized with YY1 in the promoter region of YTHDC1 (Fig. 7B). It has been reported that HDAC2 forms a complex with YY1 to regulate downstream target gene expression . Moreover, HDAC2 typically suppresses target gene expression through the deacetylation of histones [38, 39]. Thus, we speculated whether YTHDC1 was repressed by the YY1/HDAC2 complex. Consistent with previous findings, we showed that HDAC2 repression promoted YTHDC1 expression in 786-O and A498 cells (Fig. 7J and K). In contrast, overexpression of HDAC2 decreased the expression of YTHDC1 in renal cancer cells (Fig. 7 L and M). Moreover, ChIP–qPCR and ChIP-re-ChIP analysis demonstrated that the HDAC2/YY1 complex bound to the promoter region of YTHDC1 in 786-O and A498 cells (Fig. 7N and O). We also showed that coknockdown of HDAC2 and YY1 did not further increase the expression of YTHDC1 in renal cancer cells (Fig. 7P).
HDAC2 inhibitors enhance the sensitivity of ccRCC to sunitinib
To further confirm the regulation of YTHDC1 by the HDAC2/YY1 complex, 786-O and A498 cells were treated with HDAC2-specific inhibitors (CAY10683). We found that CAY10683 treatment increased the protein and mRNA levels of YTHDC1 in renal cancer cells (Fig. 8A and B). We then showed that CAY10683 abrogated the downregulation of YTHDC1 induced by overexpression of HDAC2 (Fig. 8C and D). Moreover, we demonstrated that CAY10683 treatment did not further increase the expression of YTHDC1 in the YY1 knockdown group (Fig. 7E F). In addition, we demonstrated that HDAC2 inhibitor treatment decreased the phosphorylation of ERK in renal cancer cells (Supplementary Fig. 4B). These data suggest that HDAC2 inhibitors also regulate YTHDC1 through the HDAC2/YY1 complex.
Because the downregulation of YTHDC1 contributes to sunitinib resistance in 786-O and A498 cells, we evaluated whether HDAC2 inhibitors overcome sunitinib resistance through YTHDC1. First, we demonstrated that CAY10683 treatment decreased the IC50 values in the shControl group but not in renal cancer cells after YTHDC1 knockdown (Fig. 8G). Next, the CCK-8 assay and nude mouse xenograft assay showed that CAY10683 enhanced the antitumor effect of sunitinib in wild-type renal cancer cells (Fig. 8H-K, Supplementary Fig. 5 A). However, YTHDC1 depletion attenuated the effect of CAY10683 on modulating the sensitivity of sunitinib in vitro and in vivo (Fig. 8H-K, Supplementary Fig. 5 A). Taken together, our data suggest that the HDAC2 inhibitor enhances the antitumor effect of sunitinib through the HDAC2/YY1-YTHDC1 axis in ccRCC (Fig. 8L).
To date, many functions of YTHDC1 have been reported in the literature. For example, YTHDC1 enhances exon inclusion and regulates mRNA splicing by binding with the pre-mRNA splicing factor SRSF3 (SRp20) . Additionally, YTHDC1 is reported to facilitate the nuclear export of mRNAs with m6A modification. Mechanistically, YTHDC1 can cooperate with the NXF1 and TREX mRNA export complexes to promote the export of mRNA out of the nucleus . Recently, YTHDC1 has been found to regulate the stability of nuclear mRNAs, including SQSTM1 and PTEN [24, 40]. Similarly, our results demonstrated that YTHDC1 interacts with ANXA1 mRNA via m6A modification and decreases its stability. We also explored the cancer-related role of YTHDC1 in ccRCC. Consistent with previous findings, YTHDC1 is downregulated in ccRCC tissues and functions as a tumor-inhibiting protein. However, as an m6A reader, the specific function of YTHDC1 in various types of cancer is controversial. YTHDC1 modulates MCM complex-mediated DNA replication and promotes leukemogenesis . However, YTHDC1 was reported to suppress the progression of glioma by inhibiting VPS25 . Moreover, YTHDC1 was also found to inhibit glycolysis in pancreatic cancer through the miR-30d/RUNX1 axis . The reason YTHDC1 performs different functions in tumors may be related to its differential expression in specific tumors and distinct downstream target genes.
In this study, ANXA1 was proven to be the downstream target gene of YTHDC1 in renal cancer cells. Many studies have mentioned that there is a close relationship between ANXA1 and the chemosensitivity of tyrosine kinase inhibitors [17, 43]. Consequently, we have shown that YTHDC1 regulated the sensitivity of sunitinib through ANXA1/MAPK signaling pathways. It has been well documented that the hyperactivation of MAPK/ERK signaling contributes to TKI resistance in various types of cancer [44, 45]. Inhibition of the MAPK/ERK signaling pathways by specific chemicals overcomes the resistance of TKIs . Here, YTHDC1 acted as a negative regulator of the ERK signaling pathway to sensitize ccRCC to TKI inhibitors. Thus, investigating the regulatory mechanism of YTHDC1 would help identify novel chemicals to enhance the antitumor effect of TKIs in renal cancer cells. Interestingly, we found that the HDAC2/YY1 complex suppressed YTHDC1 expression in renal cancer cells. Furthermore, HDAC2 inhibitors could promote the antitumor effect of sunitinib partially through YTHDC1 in renal cancer cells. Moreover, it has been reported that the ERK signaling pathway modulates the activation of HDAC2 . Our results suggest that HDAC2 may also regulate the activation of ERK, but the relationship between the ERK signaling pathway and HDAC2 needs to be further studied. HDAC2 is an epigenetic regulator that controls global gene expression in cells . Although we showed that HDAC2 inhibitors obtained good results in inhibiting tumor growth in cell and animal experiments, HDAC2 is highly homologous to HDAC8, which might result in unexpected side effects due to the paninhibition of HDAC . Similar to the results of DNMT inhibitors (e.g., 5-aza-dC) in clinical trials , HDAC paninhibitors approved by the Food and Drug Administration (FDA) might lead to side effects, including cardiac disturbance (arrhythmias), gastric issues (vomiting), or hematologic imbalance (anemia) . Thus, much attention should be given to the side effects caused by HDAC2 inhibitors in future therapies.
Our results demonstrated that YTHDC1 is downregulated in ccRCC tissues compared with normal tissues. Low expression of YTHDC1 is associated with a poor prognosis in patients with ccRCC. Subsequently, we showed that YTHDC1 inhibits the progression of renal cancer cells via downregulation of the ANXA1/MAPK pathways. Moreover, we also showed that the YTHDC1/ANXA1 axis modulates the sensitivity of tyrosine kinase inhibitors. We also revealed that HDAC2 inhibitors resensitize ccRCC to tyrosine kinase inhibitors through the YY1/HDAC2 complex. Taken together, we identified a novel YY1/HDAC2/YTHDC1/ANXA1 axis modulating the progression and chemosensitivity of ccRCC.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding authors on reasonable request.
Clear cell renal cell carcinoma
Tyrosine kinase inhibitors
The Cancer Genome Atlas
Cell counting kit-8
quantitative real time PCR
RNA immunoprecipitation PCR
acute myeloid leukemia
Encyclopedia of RNA Interactomes
mitogen-activated protein kinase
Institutional Animal Care and Use Committee
sodium dodecyl sulfate‒polyacrylamide gel electrophoresis
tumor necrosis factor
extracellular signal-regulated kinases
YTH domain containing 1
epidermal growth factor receptor
Wettersten HI, Aboud OA, Lara PN Jr, Weiss RH. Metabolic reprogramming in clear cell renal cell carcinoma. Nat Rev Nephrol. 2017;13:410–9.
Atkins MB, Tannir NM. Current and emerging therapies for first-line treatment of metastatic clear cell renal cell carcinoma. Cancer Treat Rev. 2018;70:127–37.
Bedke J, Gauler T, Grunwald V, Hegele A, Herrmann E, Hinz S, Janssen J, Schmitz S, Schostak M, Tesch H, et al. Systemic therapy in metastatic renal cell carcinoma. World J Urol. 2017;35:179–88.
Iacovelli R, Albiges L, Escudier B. Emerging tyrosine kinase inhibitors for the treatment of renal cancer. Expert Opin Emerg Drugs. 2015;20:379–92.
Makhov P, Joshi S, Ghatalia P, Kutikov A, Uzzo RG, Kolenko VM. Resistance to Systemic Therapies in Clear Cell Renal Cell Carcinoma: Mechanisms and Management Strategies. Mol Cancer Ther. 2018;17:1355–64.
Liu N, Pan T. N6-methyladenosine-encoded epitranscriptomics. Nat Struct Mol Biol. 2016;23:98–102.
Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. 2017;18:31–42.
Sun T, Wu R, Ming L. The role of m6A RNA methylation in cancer. Biomed Pharmacother. 2019;112:108613.
Zhou J, Wang J, Hong B, Ma K, Xie H, Li L, Zhang K, Zhou B, Cai L, Gong K. Gene signatures and prognostic values of m6A regulators in clear cell renal cell carcinoma - a retrospective study using TCGA database. Aging. 2019;11:1633–47.
Li A, Chen YS, Ping XL, Yang X, Xiao W, Yang Y, Sun HY, Zhu Q, Baidya P, Wang X, et al. Cytoplasmic m(6)A reader YTHDF3 promotes mRNA translation. Cell Res. 2017;27:444–7.
Jiang X, Liu B, Nie Z, Duan L, Xiong Q, Jin Z, Yang C, Chen Y. The role of m6A modification in the biological functions and diseases. Signal Transduct Target Ther. 2021;6:74.
He L, Li H, Wu A, Peng Y, Shu G, Yin G. Functions of N6-methyladenosine and its role in cancer. Mol Cancer. 2019;18:176.
von Hagen F, Gundert L, Strick A, Klumper N, Schmidt D, Kristiansen G, Tolkach Y, Toma M, Ritter M, Ellinger J. N(6) -Methyladenosine (m(6) A) readers are dysregulated in renal cell carcinoma. Mol Carcinog. 2021;60:354–62.
Zaccara S, Ries RJ, Jaffrey SR. Reading, writing and erasing mRNA methylation. Nat Rev Mol Cell Biol. 2019;20:608–24.
Roundtree IA, Luo GZ, Zhang Z, Wang X, Zhou T, Cui Y, Sha J, Huang X, Guerrero L, Xie P, et al: YTHDC1 mediates nuclear export of N(6)-methyladenosine methylated mRNAs. Elife. 2017;6,:e31311.
Sheng Y, Wei J, Yu F, Xu H, Yu C, Wu Q, Liu Y, Li L, Cui XL, Gu X, et al. A critical role of nuclear m6A reader YTHDC1 in leukemogenesis by regulating MCM complex-mediated DNA replication. Blood. 2021;138:2838–52.
Xiong W, Zhang B, Yu H, Zhu L, Yi L, Jin X. RRM2 Regulates Sensitivity to Sunitinib and PD-1 Blockade in Renal Cancer by Stabilizing ANXA1 and Activating the AKT Pathway. Adv Sci (Weinh). 2021;8:e2100881.
Ma G, Liang Y, Chen Y, Wang L, Li D, Liang Z, Wang X, Tian D, Yang X, Niu H. Glutamine Deprivation Induces PD-L1 Expression via Activation of EGFR/ERK/c-Jun Signaling in Renal Cancer. Mol Cancer Res. 2020;18:324–39.
Wang A, Bao Y, Wu Z, Zhao T, Wang D, Shi J, Liu B, Sun S, Yang F, Wang L, Qu L. Long noncoding RNA EGFR-AS1 promotes cell growth and metastasis via affecting HuR mediated mRNA stability of EGFR in renal cancer. Cell Death Dis. 2019;10:154.
Stoehr CG, Stoehr R, Hartmann A, Hofstaedter F, Junker K, Blaszyk H, Wieland WF, Otto W, Denzinger S, Walter B. Mutational activation of FGFR3: no involvement in the development of renal cell carcinoma. J Cancer Res Clin Oncol. 2012;138:359–61.
Wu H, He D, Biswas S, Shafiquzzaman M, Zhou X, Charron J, Wang Y, Nayak BK, Habib SL, Liu H, Li B. mTOR Activation Initiates Renal Cell Carcinoma Development by Coordinating ERK and p38MAPK. Cancer Res. 2021;81:3174–86.
Wang S, Hao X, He S, Liu C, Wang Q. Suppressive effects of RASAL2 on renal cell carcinoma via SOX2/ERK/p38 MAPK pathway. Exp Ther Med. 2020;20:151.
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.
Liang D, Lin WJ, Ren M, Qiu J, Yang C, Wang X, Li N, Zeng T, Sun K, You L, et al: m(6)A reader YTHDC1 modulates autophagy by targeting SQSTM1 in diabetic skin. Autophagy. 2021;18:1318–37.
Yamanoi M, Yamanoi K, Fujii C, Fukuda MN, Nakayama J. Annexin A1 expression is correlated with malignant potential of renal cell carcinoma. Int J Urol. 2019;26:284–90.
Raulf N, Lucarelli P, Thavaraj S, Brown S, Vicencio JM, Sauter T, Tavassoli M. Annexin A1 regulates EGFR activity and alters EGFR-containing tumour-derived exosomes in head and neck cancers. Eur J Cancer. 2018;102:52–68.
Bist P, Phua QH, Shu S, Yi Y, Anbalagan D, Lee LH, Sethi G, Low BC, Lim LH. Annexin-A1 controls an ERK-RhoA-NFkappaB activation loop in breast cancer cells. Biochem Biophys Res Commun. 2015;461:47–53.
Wang T, Kong S, Tao M, Ju S. The potential role of RNA N6-methyladenosine in Cancer progression. Mol Cancer. 2020;19:88.
Yang YH, Aeberli D, Dacumos A, Xue JR, Morand EF. Annexin-1 regulates macrophage IL-6 and TNF via glucocorticoid-induced leucine zipper. J Immunol. 2009;183:1435–45.
Zhang X, Li X, Li X, Zheng L, Lei L. ANXA1 silencing increases the sensitivity of cancer cells to low-concentration arsenic trioxide treatment by inhibiting ERK MAPK activation. Tumori. 2015;101:360–7.
Aeppli S, Schmaus M, Eisen T, Escudier B, Grunwald V, Larkin J, McDermott D, Oldenburg J, Porta C, Rini BI, et al. First-line treatment of metastatic clear cell renal cell carcinoma: a decision-making analysis among experts. ESMO Open. 2021;6:100030.
Motzer RJ, Hutson TE, Cella D, Reeves J, Hawkins R, Guo J, Nathan P, Staehler M, de Souza P, Merchan JR, et al. Pazopanib versus sunitinib in metastatic renal-cell carcinoma. N Engl J Med. 2013;369:722–31.
Joosten SC, Hamming L, Soetekouw PM, Aarts MJ, Veeck J, van Engeland M, Tjan-Heijnen VC. Resistance to sunitinib in renal cell carcinoma: From molecular mechanisms to predictive markers and future perspectives. Biochim Biophys Acta. 2015;1855:1–16.
Porta C, Paglino C, Grunwald V. Sunitinib re-challenge in advanced renal-cell carcinoma. Br J Cancer. 2014;111:1047–53.
Zhao T, Bao Y, Gan X, Wang J, Chen Q, Dai Z, Liu B, Wang A, Sun S, Yang F, Wang L. DNA methylation-regulated QPCT promotes sunitinib resistance by increasing HRAS stability in renal cell carcinoma. Theranostics. 2019;9:6175–90.
Li XL, Chen XQ, Zhang MN, Chen N, Nie L, Xu M, Gong J, Shen PF, Su ZZ, Weng X, et al. SOX9 was involved in TKIs resistance in renal cell carcinoma via Raf/MEK/ERK signaling pathway. Int J Clin Exp Pathol. 2015;8:3871–81.
Tang W, Zhou W, Xiang L, Wu X, Zhang P, Wang J, Liu G, Zhang W, Peng Y, Huang X, et al. The p300/YY1/miR-500a-5p/HDAC2 signalling axis regulates cell proliferation in human colorectal cancer. Nat Commun. 2019;10:663.
Jung KH, Noh JH, Kim JK, Eun JW, Bae HJ, Xie HJ, Chang YG, Kim MG, Park H, Lee JY, Nam SW. HDAC2 overexpression confers oncogenic potential to human lung cancer cells by deregulating expression of apoptosis and cell cycle proteins. J Cell Biochem. 2012;113:2167–77.
Hu XT, Xing W, Zhao RS, Tan Y, Wu XF, Ao LQ, Li Z, Yao MW, Yuan M, Guo W, et al. HDAC2 inhibits EMT-mediated cancer metastasis by downregulating the long noncoding RNA H19 in colorectal cancer. J Exp Clin Cancer Res. 2020;39:270.
Zhang Z, Wang Q, Zhao X, Shao L, Liu G, Zheng X, Xie L, Zhang Y, Sun C, Xu R. YTHDC1 mitigates ischemic stroke by promoting Akt phosphorylation through destabilizing PTEN mRNA. Cell Death Dis. 2020;11:977.
Zhu X, Yang H, Zhang M, Wu X, Jiang L, Liu X, Lv K. YTHDC1-mediated VPS25 regulates cell cycle by targeting JAK-STAT signaling in human glioma cells. Cancer Cell Int. 2021;21:645.
Hou Y, Zhang Q, Pang W, Hou L, Liang Y, Han X, Luo X, Wang P, Zhang X, Li L, Meng X. YTHDC1-mediated augmentation of miR-30d in repressing pancreatic tumorigenesis via attenuation of RUNX1-induced transcriptional activation of Warburg effect. Cell Death Differ. 2021;28:3105–24.
Chuang MC, Lung JH, Chen YC, Lin YC, Li YC, Hung MS. The Association of Annexin A1 and Chemosensitivity to Osimertinib in Lung Cancer Cells. Cancers (Basel) 2021;13:4106.
Lu Y, Liu Y, Oeck S, Zhang GJ, Schramm A, Glazer PM. Hypoxia Induces Resistance to EGFR Inhibitors in Lung Cancer Cells via Upregulation of FGFR1 and the MAPK Pathway. Cancer Res. 2020;80:4655–67.
Cao M, Nawalaniec K, Ajay AK, Luo Y, Moench R, Jin Y, Xiao S, Hsiao LL, Waaga-Gasser AM. PDE4D targeting enhances anti-tumor effects of sorafenib in clear cell renal cell carcinoma and attenuates MAPK/ERK signaling in a CRAF-dependent manner. Transl Oncol. 2022;19:101377.
Wang W, Xia X, Chen K, Chen M, Meng Y, Lv D, Yang H. Reduced PHLPP Expression Leads to EGFR-TKI Resistance in Lung Cancer by Activating PI3K-AKT and MAPK-ERK Dual Signaling. Front Oncol. 2021;11:665045.
Liu XY, Yao B, Hao JR, Jin L, Gao Y, Yang X, Liu L, Sun XY, Sun N, Gao C. IQGAP1/ERK regulates fear memory formation via histone posttranslational modifications induced by HDAC2. Neurobiol Learn Mem. 2020;171:107210.
Shetty MG, Pai P, Deaver RE, Satyamoorthy K, Babitha KS. Histone deacetylase 2 selective inhibitors: A versatile therapeutic strategy as next generation drug target in cancer therapy. Pharmacol Res. 2021;170:105695.
Yang Y, Hu B, Yang Y, Gong K, Wang H, Guo Q, Tang X, Li Y, Wang J. Rational design of selective HDAC2 inhibitors for liver cancer treatment: computational insights into the selectivity mechanism through molecular dynamics simulations and QM/MM calculations. Phys Chem Chem Phys. 2021;23:17576–90.
Fandy TE. Development of DNA methyltransferase inhibitors for the treatment of neoplastic diseases. Curr Med Chem. 2009;16:2075–85.
This work was supported by grants from the Chinese National Natural Science Foundation Grant No. 82103341 (W.X.).
Ethics approval and consent to participate
The study was conducted in accordance with the principles of the Declaration of Helsinki principles. It was approved by the Animal Use and Care Committees at the Second Xiangya hospital, Central South University.
Consent for publication
All subjects have written informed consent.
The authors have declared that no competing interest exists.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, W., Ye, K., Li, X. et al. YTHDC1 is downregulated by the YY1/HDAC2 complex and controls the sensitivity of ccRCC to sunitinib by targeting the ANXA1-MAPK pathway. J Exp Clin Cancer Res 41, 250 (2022). https://doi.org/10.1186/s13046-022-02460-9