- Open Access
DYRK1A reinforces epithelial-mesenchymal transition and metastasis of hepatocellular carcinoma via cooperatively activating STAT3 and SMAD
Journal of Biomedical Science volume 29, Article number: 34 (2022)
Hepatocellular carcinoma (HCC) accounts for the majority of liver cancer cases, while metastasis is considered the leading cause of HCC-related death. However, the currently available treatment strategies for efficient suppression of metastasis are limited. Therefore, novel therapeutic targets to inhibit metastasis and effectively treat HCC are urgently required.
Wound healing and Transwell assays were used to determine the migration and invasion abilities of HCC cells in vitro. Quantitative real-time PCR (qRT-PCR), protein array, immunofluorescence, and immunoprecipitation experiments were used to study the mechanism of DYRK1A-mediated metastasis. A tail vein metastasis model and H&E staining were utilized to assess metastatic potential in vivo.
The results of the current study demonstrated that dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) was upregulated in HCC tissues compared with normal liver tissues. Additionally, the level of DYRK1A was increased in primary HCC tissues of patients with metastasis compared with those of patients without metastasis, and DYRK1A overexpression correlated with worse outcomes in liver cancer patients. Gain- and loss-of-function studies suggested that DYRK1A enhanced the invasion and migration abilities of HCC cells by promoting epithelial-mesenchymal transition (EMT). Regarding the promoting effect of DYRK1A on cell invasion, the results showed that DYRK1A was coexpressed with TGF-β/SMAD and STAT3 signalling components in clinical tumour samples obtained from patients with HCC. DYRK1A also activated TGF-β/SMAD signalling by interacting with tuberous sclerosis 1 (TSC1) and enhanced metastasis of HCC cells by activating STAT3. Furthermore, DYRK1A promoted EMT by cooperatively activating STAT3/SMAD signalling.
Overall, the present study not only uncovered the promoting effect of DYRK1A on HCC metastasis and revealed the mechanism but also provided a new approach to predict and treat metastatic HCC.
The high mortality rate of liver cancer is primarily due to metastasis and postsurgical recurrence . The majority of patients with hepatocellular carcinoma (HCC), the most common type of liver cancer, are diagnosed at an advanced stage of the disease with distant metastases, when the available treatment options are limited [2, 3]. Therefore, identifying key molecules involved in HCC progression is of great importance for the development of novel therapeutic strategies. Epithelial-mesenchymal transition (EMT) is a key step in the early stages of cancer metastasis . The Janus kinase (JAK)/signal transducer and activator of transcription 3 (STAT3) and TGF-β1-Smad2/3 signalling pathways, important cascades for signal transduction from plasma membrane receptors to the nucleus, synergistically augment EMT and promote metastasis . These two signalling pathways participate in crosstalk during EMT in HCC: STAT3 regulates Smad3-mediated signalling induced by TGF-β, while Smad3 regulates STAT3-mediated signalling induced by various cytokines . Thus, it is critical to develop novel therapeutic molecules to target these two signalling pathways linked to EMT and metastasis in HCC.
Dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A), located in the Down’s syndrome critical region, plays pivotal roles in neurodegenerative diseases, such as Alzheimer’s disease and Down’s syndrome . Furthermore, it has been suggested that targeting DYRK1A is a potential strategy for treating cancer . For instance, a study revealed that DYRK1A can modulate the expression of c-mesenchymal-epithelial transition receptor (c-MET) and drive tumour growth, and it was therefore considered a target for treating pancreatic ductal adenocarcinoma . In addition to mediating tumour growth, DYRK1A is also involved in the regulation of several cellular processes associated with cancer progression, such as cancer stemness maintenance and cancer cell invasion .
Tumour metastasis, the major cause of mortality in patients with cancer, is a multistep process that involves tumour cell invasion through the basal membrane into blood vessels followed by extravasation and colonization of host organs [11, 12]. Angiogenesis represents the process by which new blood capillaries are formed from preexisting vasculature, thus promoting primary tumour growth and metastasis . A study demonstrated that DYRK1A positively regulated VEGF-dependent nuclear factor of activated T cells (NFAT) transcriptional responses in endothelial cells and was involved in angiogenesis . In addition, DYRK1A enhanced the migration ability of glioblastoma cells by activating NFATC1 in vitro . Although current evidence has highlighted the importance of DYRK1A in tumour metastasis, its promoting effects on tumour metastasis, as well as the underlying mechanism, need further investigation. Furthermore, EMT is involved in metastatic dissemination of carcinomas . The present study first demonstrated that DYRK1A can promote HCC cell EMT and metastasis both in vitro and in vivo.
Materials and methods
Transforming growth factor-β (TGF-β; catalogue number: 100–21) was obtained from PeproTech (Rocky Hill, NJ, USA). We obtained harmine (catalogue number: HY-N0737A) from MedChemExpress (Monmouth Junction, NJ, USA). Anti-N-cadherin (14215S), anti-E-cadherin (3195P), anti-DYRK1A (2771), anti-p-Smad3 (9520), anti-Slug (9585) and anti-Snail (3879S) antibodies were obtained from Cell Signaling Technology (Danvers, MA, USA). Anti-β-actin (47778), anti-Smad2/3 (133098), and anti-STAT3 (482) antibodies were purchased from Santa Cruz Biotechnology (Dallas, Texas, USA). The anti-p-STAT3 (Tyr705) (ab76315) and anti-tuberous sclerosis 1 (TSC1; ab270967) antibody was obtained from Abcam (Cambridge, Cambridgeshire, UK). The anti-α-tubulin (AF0001) antibody was purchased form Beyotime Institute of Biotechnology (Shanghai, China).
HepG2 cells and Hep3B cells were purchased from Shanghai Institute of Biochemistry and Cell Biology and cultured in Dulbecco’s modified Eagle medium (DMEM) containing 10% foetal bovine serum (FBS). The human foetal hepatocyte line L-02 was cultured in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% FBS.
Wound healing assay
Treated HCC cells (8 × 105 cells/well) were incubated in 24-well plates. Following 24 h of incubation, the cells were pretreated with 5 μM mitomycin C for 3 h to suppress proliferation. Subsequently, the HCC cell monolayer was scratched using a sterile 10-µL pipette tip and was then washed with PBS. After 24 h, images of the migrated cells were acquired under a microscope. The migrating cells in five distinct scratched areas per well were counted, and the migration rate was calculated.
Transwell migration and invasion assays
Briefly, following serum starvation overnight, HCC cells (1 × 105 cells/mL) in serum-free medium were seeded in the Transwell chambers, with or without a Matrigel coating on the insert membrane, for the invasion or migration assay, respectively. Then, 600 μL of DMEM supplemented with 20% FBS was added to the lower compartments of the Transwell chambers. Following incubation at indicated times, cells on the bottom surface of the membrane were fixed and stained with crystal violet solution containing 20% methanol for 30 min. Quantification of migrated and invaded cells was performed with ImageJ software.
Transfection and DYRK1A short hairpin RNA (shRNA) construction
HCC cells were incubated in six-well plates (2 × 105 cells/well) overnight, and transfection of siRNA or plasmids was then performed with jetPRIME (Polyplus, NY, USA) following the manufacturers’ protocols. The sense sequences of the DYRK1A siRNAs were 5′-AUGGAGCUAUGGACGUUAATT-3′ (DYRK1A siRNA-1) and 5′-AAACUCGAAUUCAACCUUATT-3′ (DYRK1A siRNA-2), that of the STAT3 siRNA was 5′- CCACUUUGGUGUUUCAUAATT-3′, that of the TSC-1 siRNA was 5′-CGGCUGAUGUUGUUAAAUATT-3′, and that of the negative control siRNA was 5′-UUCUCCGAACGUGUCACGUTT-3′.
The following oligonucleotide sequence (forward primer) was used to synthesize the DYRK1A shRNA constructs (the targeted sequence is bolded). 5′-CCGGGTTCGGCTTGCACCGTCATTTCTCGAGAAATGACGGTGCAAGCCGAACTTTTTG-3′. The reverse primer sequence was 5′-AATTCAAAAAGTTCGGCTTGCACCGTCATTTCTCGAGAAATGACGGTGCAAGCCGAAC-3′. The designed shRNA contained 58 bp with 21 bp of perfectly matched nucleotide base pairs connected by a 6 bp loop [CTCGAG]. The shRNA targeting DYRK1A was synthesized and inserted into the pLKO.1-TRC Cloning vector empty vector according to the manufacturer’s instructions. The DYRK1A shRNA and NC shRNA were transfected into HepG2 cells using Lipofectamine 3000 reagent (Invitrogen, Shanghai, China) according to the manufacturer’s instructions. After transfection, HepG2 cells were selected using puromycin (10 μg/ml) for 5 days. Puromycin-resistant colonies were then single-cell cloned and expanded. The expression of DYRK1A was determined by western blot analysis.
HCC cells were first lysed in RIPA buffer, and samples containing ~ 20 μg of protein were then separated by SDS–PAGE and transferred to a PVDF membrane. Subsequently, the membrane was blocked using 5% (w/v) skim milk, incubated with primary antibodies, and then incubated with the corresponding secondary antibodies. Finally, the blots were visualized using an imaging system.
For immunofluorescence analysis, treated HCC cells were incubated in a culture dish overnight at 37 °C. Following incubation, HCC cells were incubated with 4% paraformaldehyde for 15 min and permeabilized using 0.3% Triton X-100 for 30 min. After blocking using 1% BSA for 30 min, the cells were incubated with specific antibodies at 4 °C for 24 h and were then stained with the corresponding isotype-specific secondary antibodies for 1 h in the dark. Nuclei were stained with 0.1% DAPI solution, and staining was analysed using a fluorescence microscope.
Total RNA was isolated with TRIzol, sedimented with isopropyl alcohol, washed with 70% ethanol two times, and reverse transcribed into cDNA. SYBR Green-based qRT–PCR was used to examine mRNA expression levels. The primers were as follows: DYRK1A, 5′-TCTGGGTATTCCACCTGCTC-3′ (forward) and 5′-GTCCTCCTGTTTCCACTCCA-3′ (reverse); GAPDH, 5′-GAGTCAACGGATTTGGTCGT-3′ (forward) and 5′-TTGATTTTGGAGGGATCTCG-3′ (reverse); N-cadherin, 5′-CAACTTGCCAGAAAACTCCAGG-3′ (forward) and 5′-ATGAAACCGGGCTATCTGCTC-3′ (reverse); E-cadherin, 5′-GAACGCATTGCCACATACAC-3′ (forward) and 5′-GAATTCGGGCTTGTTGTCAT-3′ (reverse).
An immunoprecipitation assay was performed to evaluate protein–protein interactions. Briefly, HCC cells were lysed using an NP-40 solution supplemented with protease inhibitor cocktail. After centrifugation at 12,000×g for 15 min, the lysates were incubated with 20 μL of protein A/G beads for 1 h. Subsequently, following centrifugation, the supernatants were incubated first with 5 μL of a specific primary antibody for 10 h and then with protein A/G magnetic beads for 2 h. Finally, the beads were analysed using western blotting.
Tail vein injection model of metastasis in mice
Four- to six-week-old female BALB/c nude mice were purchased from Shanghai SLAC laboratory animal company and randomly divided into two groups (n = 6 mice per group). Subsequently, 1 × 106 HCC cells suspended in 100 μL of PBS were injected into the tail veins of nude mice. Sixty days after injection, the nude mice were sacrificed to evaluate cell metastasis in vivo. Briefly, tissues were embedded in paraffin and were then stained with haematoxylin and eosin (H&E) to evaluate cell metastasis. Lung metastasis was detected by staining with Bouin’s solution and H&E. Neutral balsam was used to preserve the samples. The experiments were performed in compliance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals.
Data are presented as the mean ± standard deviation values. Two‐tailed Student’s t‐test was used to determine significant differences. *p < 0.05; **p < 0.01; ***p < 0.001.
DYRK1A plays an essential role in the malignant progression and metastasis of HCC
To investigate the key genes involved in HCC metastasis, we established a series of HepG2 cells with different invasion abilities by serial selection in a Transwell assay (Additional file 1: Fig. S1a). RNA sequencing was performed to compare the gene transcription levels between highly invasive HepG2 cells (P5) and weakly invasive HepG2 cells (P0), and DYRK1A was found to be overexpressed in P5 compared with P0 HepG2 cells (Additional file 6: Table S1). Furthermore, bioinformatic analysis using the Human Protein Atlas database revealed that DYRK1A was upregulated in HCC tissues compared with normal tissues (Fig. 1a and b, Additional file 2: Fig. S2a and b; p < 0.001) [17, 18]. Moreover, DYRK1A was overexpressed in HCC cells compared with normal foetal hepatocytes (Fig. 1c). Additionally, further bioinformatic analysis with the UALCAN database showed that the mRNA level of DYRK1A was increased in HCC tissues compared with normal liver tissues (Fig. 1d). The level of DYRK1A was correlated with the lymph node metastasis status in HCC patients (Fig. 1e) . DYRK1A was also upregulated in primary metastatic HCC tissues compared with nonmetastatic HCC tissues (Fig. 1f; p < 0.05) . Additionally, alteration of DYRK1A expression was correlated with poor outcomes in patients with liver cancer (n = 365: n = 9 in the altered group and n = 356 in the unaltered group; p < 0.05; Fig. 1g) [21, 22]. Moreover, high DYRK1A levels were associated with shorter progression-free survival and progression/metastasis-free survival times in liver cancer patients (Additional file 1: Fig. S1b) . Therefore, we hypothesized that DYRK1A upregulation could be associated with malignant progression and especially metastasis in HCC. Gene–gene Pearson correlation analysis with DYRK1A in HCC clinical samples and functional enrichment analysis revealed that DYRK1A expression was strongly associated with the expression of genes enriched in the cell migration, cell motility and cell–cell junction ontologies, suggesting that DYRK1A could regulate the invasion of HCC cells (Additional file 1: Fig. S1c and Additional file 7: Table S2) .
DYRK1A promotes the migration and invasion abilities of HCC cells
DYRK1A has been reported to promote cell proliferation and tumour growth in several cancer types [9, 25]. However, the effect of DYRK1A on the proliferation of HCC cells remains unclear. Proliferation was evaluated in HepG2 and Hep3B cells transfected with a DYRK1A small interfering RNA (siRNA) to knock down DYRK1A expression (Fig. 2a). DYRK1A knockdown did not efficiently suppress cell proliferation compared with that in the control siRNA group, and DYRK1A overexpression also could not promote the proliferation of HCC cells (Additional file 3: Fig. S3a and b). However, DYRK1A knockdown inhibited the migration and invasion of HCC cells (Fig. 2b and c). Additionally, treatment of HepG2 cells with harmine, a DYRK1A inhibitor, notably attenuated their migration and invasion in a dose-dependent manner but had little effect on their proliferation (Fig. 2d and Additional file 3: Fig. S3c). Consistent with these findings, the results of wound healing assays indicated the potent antimetastatic effect of DYRK1A depletion in HCC cells (Fig. 2e and f). To further evaluate the role of DYRK1A in metastasis in vivo, we investigated whether DYRKA knockdown can suppress HCC metastasis in vivo. To this end, DYRKA-knockdown HepG2 cells and parental HepG2 cells were intravenously injected into nude mice. As expected, the number of pulmonary metastatic nodules was decreased in the DYRK1A-knockdown cell group compared with the control group (Fig. 2g and h). In contrast, ectopic DYRK1A overexpression clearly enhanced the migration and invasion of HCC cells (Fig. 3a and b). In addition, wound healing assays showed that DYRK1A overexpression significantly increased the migration ability of HCC cells (Fig. 3c). These findings indicated that DYRK1A could promote the metastasis of HCC cells.
DYRK1A promotes EMT in HCC cells
Epithelial cells go through EMT to acquire mesenchymal characteristics following the dissociation of cell–cell junctions which is crucial in cancer invasion and metastasis, and the hallmark of EMT is the upregulation of N-cadherin followed by the downregulation of E-cadherin . Our data indicated that siRNA- and harmine-induced DYRK1A suppression increased the level of the epithelial marker E-cadherin and reduced the level of the mesenchymal marker N-cadherin, suggesting that DYRK1A silencing can inhibit the EMT process in HCC cells (Fig. 3d and e). In contrast, ectopic DYRK1A overexpression significantly upregulated N-cadherin and downregulated E-cadherin (Fig. 3f). Moreover, in HepG2 and Hep3B cells, DYRK1A silencing significantly increased E-cadherin and reduced N-cadherin mRNA expression, and DYRK1A overexpression increased N-cadherin and decreased E-cadherin mRNA expression (Fig. 3g and h). These findings indicated that DYRK1A could enhance the metastatic potential of HCC cells by promoting EMT.
The TGF-β/SMAD and JAK/STAT3 pathways are highly enriched in HCC clinical samples with increased DYRK1A expression
To better understand the mechanism underlying the effect of DYRK1A on HCC progression, gene set enrichment analysis (GSEA) was performed in HCC samples with low and high DYRK1A expression. GSEA of pathways in the WikiPathways and PANTHER pathway databases was performed using the LinkedOmics platform (Fig. 4a and b; Additional file 8: Table S3 and Additional file 9: Table S4) . As shown in Fig. 4c, bioinformatic analysis with the WikiPathway and PANTHER pathway databases revealed that the TGF-β/SMAD, IL-6/JAK/STAT3 and EGF/EGFR pathways were significantly enriched in the samples with high DYRK1A expression. A previous study by our laboratory demonstrated that DYRK1A could regulate STAT3/EGFR/MET signalling and sensitize EGFR wild-type non-small cell lung cancer (NSCLC) cells to AZD9291, thus indicating that DYRK1A could regulate EGFR via STAT3 . The results of the present study showed that the expression of DYRK1A in clinical samples obtained from patients with HCC could be associated with the activity of DNA-binding transcription factors. Therefore, the SMAD and STAT3 transcription factors could play a critical role in regulating tumorigenesis (Additional file 1: Fig. S1b). These results suggested that DYRK1A could regulate metastasis and EMT via the TGF-β/SMAD and IL-6/JAK/STAT3 signalling pathways (Fig. 4d and e).
DYRK1A promotes TGF-β/SMAD signalling by interacting with TSC1
To gain mechanistic insights into the effects of DYRK1A on activating TGF-β/SMAD signalling and promoting metastasis, the proteins interacting with DYRK1A were analysed using mass spectrometry; in addition, genes coexpressed with DYRK1A, SMAD2 and SMAD3 in liver cancer were identified with the UALCAN platform as those with Pearson correlation coefficient (r) values greater than 0.5 (Additional file 10: Table S5) [19, 28]. To identify the key genes associated with DYRK1A-mediated promotion of TGF-β/SMAD signalling, we generated a Venn diagram, which showed that there were 25 overlapping genes among those coexpressed with DYRK1A, SMAD2 and SMAD3 in liver cancer and DYRK1A-interacting proteins (Fig. 5a) . Furthermore, among the 25 overlapping genes, three genes (DCAF7, EPS15, and TSC1) were overexpressed both in liver cancer tissue compared with normal tissue and in primary HCC tissues of patients with metastasis compared with primary HCC tissues of patients without metastasis, as determined using the online analysis tool HCMDB, indicating that these genes might be involved in the development of HCC (Additional file 4: Fig. S4a). In addition, TSC1 overexpression resulted in poor outcomes in liver cancer patients (p < 0.05), but DCAF7 and EPS15 might not be associated with poor outcomes in liver cancer patients (both p > 0.05; Additional file 4: Fig. S4b) [23, 30]. More importantly, a previous study reported that TSC1 activates TGF-β-Smad2/3 signalling and promotes the EMT process . Thus, we sought to investigate the role of the interaction between DYRK1A and TSC1 in promoting the activation of TGF-β/SMAD signalling and metastasis.
Immunoprecipitation assays verified that DYRK1A can interact with TSC1 in HCC cells (Fig. 5b). Furthermore, DYRK1A and TSC1 were coexpressed in clinical tumour samples from patients with HCC (Spearman correlation analysis; r = 0.61, p < 9.4e–39; Pearson correlation analysis, r = 0.62, p < 0.0001; Fig. 5c) . Therefore, it was hypothesized that DYRK1A could activate TGF-β/SMAD signalling by interacting with TSC1. Thus, the associations of the expression of TGF-β/SMAD signalling-related genes and the expression of DYRK1A with tumour metastasis were first investigated. As shown in Fig. 5d, DYRK1A was coexpressed with TGF-β/SMAD signalling-related genes, such as TGF-β receptor 1, TGF-β receptor 2, SMAD2 and SMAD3, in clinical samples from patients with HCC . Additionally, immunoprecipitation assays revealed that DYRK1A can interact with TSC1 and SMAD2/3, indicating that the DYRK1A-TSC1-SMAD2/3 complex might play a vital role in HCC tumorigenesis (Fig. 5e). Furthermore, DYRK1A knockdown in HCC cells by siRNA or inhibition with harmine suppressed the expression of TSC1/TGF-β/SMAD signalling-related genes, such as TSC1, snail, SMAD2 and SMAD3 (Fig. 5f and g). Meanwhile, TSC1 knockdown downregulated DYRK1A and SMAD2/3 (Additional file 5: Fig. S5). In contrast, DYRK1A overexpression increased the expression levels of TSC1/TGF-β/SMAD signalling-related genes (Fig. 5h). Finally, DYRK1A silencing or cell treatment with harmine attenuated the ability of TGF-β to activate the SMAD signalling pathway (Fig. 5i and j). Overall, these results suggested that DYRK1A could activate the TGF-β/SMAD pathway by interacting with TSC1 and SMAD2/3.
DYRK1A promotes EMT and enhances the metastatic potential of HCC cells through the TSC1/TGF-β/SMAD pathway
HCC cell treatment with harmine significantly attenuated the TGF-β-triggered migration of HepG2 and Hep3B cells (Fig. 6a). Subsequently, our results demonstrated that DYRK1A silencing or cell treatment with harmine inhibited TGF-β-induced EMT (Fig. 6b and c). In addition, TSC1 silencing restored the suppressive effect of DYRK1A knockdown on the EMT process in HCC cells (Fig. 6d). Furthermore, TSC1 knockdown enhanced the inhibitory effect of DYRK1A depletion on the migration ability of HCC cells (Fig. 6e). Moreover, TSC1 silencing reversed the enhancement of the EMT process and SMAD3 activation induced by DYRK1A overexpression in HCC cells (Fig. 6f). In addition, TSC1 and DYRK1A overexpression could enhance the activation of EMT process and SMAD3 compared with either TSC1 or DYRK1A overexpression (Fig. 6g). Therefore, these results indicated that DYRK1A could enhance EMT and the metastatic potential of HCC cells via the TSC1/TGF-β/SMAD signalling pathway.
DYRK1A promotes EMT and the metastatic potential of HCC cells by activating STAT3
A previous study by our laboratory showed that DYRK1A could affect the antitumor activity of AZD9291 by regulating STAT3 signalling in NSCLC . Therefore, it was hypothesized that DYRK1A could enhance HCC cell metastasis by regulating STAT3. DYRK1A was found to be coexpressed with several STAT3 target genes, including zinc finger E-box-binding homeobox 1, VEGFA and heat shock protein 90, in clinical tumour samples from patients with HCC (Fig. 7a). Moreover, the interaction between DYRK1A and STAT3 in HCC cells was verified by immunoprecipitation, as shown in Fig. 7b. Furthermore, DYRK1A deletion by siRNA or harmine significantly inhibited the activation of STAT3 in HCC cells (Fig. 7c and d). In contrast, DYRK1A overexpression notably upregulated STAT3 (Tyr705) phosphorylation (Fig. 7e), while STAT3 depletion attenuated DYRK1A-induced EMT in HCC cells (Fig. 7f). Furthermore, STAT3 overexpression abrogated the inhibitory effect of DYRK1A depletion on the migration of HCC cells (Fig. 7g and h). These findings indicated that DYRK1A could promote EMT and the metastatic potential of HCC cells by activating STAT3.
DYRK1A promotes EMT by cooperatively activating the STAT3/SMAD axis
Here, we hypothesized that DYRK1A-induced EMT was mediated not only by STAT3 but also potentially by TSC1/TGF-β/SMAD signalling. To assess whether DYRK1A can promote the interaction between STAT3 and SMAD2/3, an immunofluorescence assay was performed. DYRK1A depletion not only decreased the nuclear translocation of STAT3 and SMAD2/3 but also inhibited the colocalization of STAT3 and SMAD3 in the nucleus in HCC cells (Fig. 8a). However, as shown in Fig. 8b, the direct interaction of STAT3 and SMAD2/3 was not observed both on empty vector and DYRK1A overexpression groups in HCC cells. These data indicated that DYRK1A promoted the activation of STAT3 and SMAD2/3 independently, and the interaction of STAT3 and SMAD2/3 could not be observed in DYRK1A overexpression group. However, STAT3 overexpression restored DYRK1A depletion-induced SMAD3 downregulation in HCC cells (Fig. 8c). These results indicated that DYRK1A could promote EMT by activating a cooperative STAT3/SMAD gene transcription program.
DYRK1A has been reported to play a multifaceted role in tumorigenesis. Several studies have indicated that DYRK1A acts as a tumour suppressor gene, while others have highlighted its pro-oncogenic activity . For instance, DYRK1A could activate NFATC1 to increase the migration ability of glioblastoma cells, while it could play an important role in maintaining cancer stemness in oral squamous cell carcinoma [15, 33]. In contrast, another study showed that DYRK1A could act as a tumour suppressor by downregulating c-Myc in acute myeloid leukaemia cells . Therefore, it was hypothesized that the effect of DYRK1A on tumorigenesis was highly dependent on the cellular context . Previous studies from our laboratory demonstrated that DYRK1A could act as an oncogene in lung tumorigenesis [27, 36, 37]. Therefore, exploring the role of DYRK1A in HCC tumorigenesis is of great importance. Selective inhibitors of DYRK1A have been extensively explored and are considered potential candidates for treating multiple diseases, including Down’s syndrome, Alzheimer’s disease, Parkinson’s disease and cancer . Therefore, uncovering the effect of DYRK1A on HCC tumorigenesis could help identify cancer patients who are most likely to benefit from treatment with DYRK1A inhibitors. The current study first demonstrated that DYRK1A was upregulated in HCC tissues compared with normal liver tissues. Furthermore, DYRK1A overexpression enhanced the metastatic ability of HCC cells both in vitro and in vivo. These findings indicated that DYRK1A could function as an oncogene in HCC tumorigenesis. A recent study by Li et al.  showed that DYRK1A could directly bind to trophinin-associated protein to promote cell cycle progression in HCC cells by activating the Akt/GSK-3β signalling pathway. The results of the current study highlighting the role of DYRK1A in promoting HCC tumorigenesis could also be supported by the above study.
The JAK/STAT3 and TGF-β/SMAD signalling pathways play critical roles in regulating cancer cell metastasis and can synergistically enhance EMT and the metastatic ability of cancer cells . Therefore, it is important to develop molecular drugs that can regulate both the JAK/STAT3 and TGF-β/SMAD signalling pathways to treat patients with metastatic cancer. Herein, DYRK1A was found to be coexpressed with members of the TGF-β/SMAD and STAT3 signalling pathways in clinical tumour samples from patients with HCC. Furthermore, DYRK1A activated TGF-β/SMAD signalling by interacting with TSC1 and enhanced the metastatic ability of HCC cells by interacting and activating STAT3. Recently, Bhansali et al.  demonstrated that DYRK1A could regulate the development of B-cell acute lymphoblastic leukaemia by binding to and phosphorylating STAT3. This work was consistent with and supported our findings. In addition, the results of the current study suggested that DYRK1A could induce EMT by promoting a cooperative STAT3/SMAD gene transcription program. Although STAT3 and SMAD are attractive therapeutic targets for cancer, they have been shown to be notoriously difficult to target with small molecule inhibitors in clinical trials [41, 42]. In addition, the results of the present work suggested that DYRK1A depletion could simultaneously inhibit the activation of JAK/STAT3 and TGF-β/SMAD signalling. Therefore, DYRK1A could be considered a novel therapeutic target for metastatic HCC. When STAT3 is aberrantly activated in tumours, STAT3 selectively and directly interacts with Smad3, sequesters Smad3 from the Smad nucleoprotein complex and thus antagonizes TGF-β signalling . Thus, the physically interaction of STAT3 and SMAD3 can attenuate the process of EMT, our data indicated that the interaction of STAT3 and SMAD3 could not be observed in HCC cells, and DYRK1A overexpression did not affect the interaction of STAT3 and SMAD3 in HCC cells. These data suggested that DYRK1A cooperatively activated STAT3 and SMAD signalling in an independent way. However, the crosstalk between STAT3 and SMAD signalling activated by DYRK1A could be observed in HCC cells, which is consistence with literature . For example, STAT3 overexpression upregulated SMAD2/3 in HCC cells, while STAT3 silencing exhibited the opposite effect (Additional file 5: Fig. S5b). However, the mechanism of crosstalk between STAT3 and Smad signalling activated by DYRK1A may need further investigation.
It has been reported that the interaction of TSC1 with the TGF-β receptor complex and SMAD2/3 is required for TGF-β1-mediated SMAD2/3 phosphorylation and EMT . Furthermore, PI3K activation and phosphorylation of Akt at Thr308 can result in phosphorylation and inhibition of the TSC1 complex. Therefore, Akt was identified as the first kinase to directly phosphorylate and negatively regulate the TSC1 complex [44, 45]. Mature DYRK1A, a dual-specificity protein kinase, phosphorylates substrates on serine or threonine (S/T) residues . Several kinases have been reported to phosphorylate TSC1 on S/T residues. For example, polo-like kinase 1 can phosphorylate TSC1 at S467 and S578 to destabilize TSC1 . The current study revealed that DYRK1A can directly bind to TSC1 and positively regulate its expression, thus suggesting that TSC1 could be a novel substrate of DYRK1A in HCC cells. Additionally, DYRK1A can bind to and regulate SMAD2/3 by directly interacting with TSC1 in HCC cells. However, the binding sites of between DYRK1A and TSC1 should be further determined. It has been reported that DYRK1A increases the stability of its substrates, such as NFATC1, by phosphorylating and interfering with their ubiquitin-mediated proteasomal degradation . Our data indicated that DYRK1A enhanced the expression of TSC1 and that TSC1 enhanced the expression of DYRK1A. Thus, we hypothesized that DYRK1A phosphorylated and inhibited the ubiquitin-mediated proteasomal degradation of TSC1, subsequently increasing the stability of TSC1. However, further investigation might be needed to verify this hypothesis.
Collectively, this work suggested that DYRK1A could be a potential biomarker for the diagnosis and prognosis of patients with HCC. In addition, DYRK1A could promote HCC tumorigenesis by enhancing the metastatic ability of HCC cells. Furthermore, DYRK1A promoted HCC cell EMT and metastasis by activating a cooperative STAT3/SMAD gene transcription program (Fig. 9). Therefore, this study not only uncovered the effect of DYRK1A on promoting HCC metastasis and revealed the underlying mechanism but also provided a novel therapeutic approach for improving the prognosis and treatment of patients with metastatic HCC.
Availability of data and materials
The data presented in the study are included in the article and additional material.
Dual-specificity tyrosine phosphorylation-regulated kinase 1A
C-mesenchymal-epithelial transition receptor
Transforming growth factor-β
Tuberous sclerosis 1
Fetal bovine serum
Hematoxylin and eosin
Small interfering RNA
Gene set enrichment analysis
Serine or threonine
Dulbecco’s modified Eagle medium
Signal transducer and activator of transcription 3
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(2):529–43.
Xu M, Wang J, Xu Z, Li R, Wang P, Shang R, Cigliano A, Ribback S, Solinas A, Pes GM, Evert K, Wang H, Song X, Zhang S, Che L, Pascale RM, Calvisi DF, Liu Q, Chen X. SNAI1 promotes the cholangiocellular phenotype, but not epithelial-mesenchymal transition, in a murine hepatocellular carcinoma model. Cancer Res. 2019.
Starley BQ, Calcagno CJ, Harrison SA. Nonalcoholic fatty liver disease and hepatocellular carcinoma: a weighty connection. Hepatology. 2010;51(5):1820–32.
Liu RY, Zeng Y, Lei Z, Wang L, Yang H, Liu Z, Zhao J, Zhang HT. JAK/STAT3 signaling is required for TGF-beta-induced epithelial-mesenchymal transition in lung cancer cells. Int J Oncol. 2014;44(5):1643–51.
Wang B, Liu T, Wu JC, Luo SZ, Chen R, Lu LG, Xu MY. STAT3 aggravates TGF-beta1-induced hepatic epithelial-to-mesenchymal transition and migration. Biomed Pharmacother. 2018;98:214–21.
Itoh Y, Saitoh M, Miyazawa K. Smad3-STAT3 crosstalk in pathophysiological contexts. Acta Biochim Biophys Sin. 2018;50(1):82–90.
Czarna A, Wang J, Zelencova D, Liu Y, Deng X, Choi HG, Zhang T, Zhou W, Chang JW, Kildalsen H, Seternes OM, Gray NS, Engh RA, Rothweiler U. Novel scaffolds for dual specificity tyrosine-phosphorylation-regulated kinase (DYRK1A) inhibitors. J Med Chem. 2018;61(17):7560–72.
Laham AJ, Saber-Ayad M, El-Awady R. DYRK1A: a down syndrome-related dual protein kinase with a versatile role in tumorigenesis. Cell Mol Life Sci CMLS. 2021;78(2):603–19.
Luna J, Boni J, Cuatrecasas M, Bofill-De Ros X, Nunez-Manchon E, Gironella M, Vaquero EC, Arbones ML, de la Luna S, Fillat C. DYRK1A modulates c-MET in pancreatic ductal adenocarcinoma to drive tumour growth. Gut. 2019;68(8):1465–76.
Lee SB, Frattini V, Bansal M, Castano AM, Sherman D, Hutchinson K, Bruce JN, Califano A, Liu G, Cardozo T, Iavarone A, Lasorella A. An ID2-dependent mechanism for VHL inactivation in cancer. Nature. 2016;529(7585):172–7.
Steinbichler TB, Dudas J, Riechelmann H, Skvortsova II. The role of exosomes in cancer metastasis. Semin Cancer Biol. 2017;44:170–81.
van Zijl F, Krupitza G, Mikulits W. Initial steps of metastasis: cell invasion and endothelial transmigration. Mutat Res. 2011;728(1–2):23–34.
Lou W, Liu J, Gao Y, Zhong G, Chen D, Shen J, Bao C, Xu L, Pan J, Cheng J, Ding B, Fan W. MicroRNAs in cancer metastasis and angiogenesis. Oncotarget. 2017;8(70):115787–802.
Rozen EJ, Roewenstrunk J, Barallobre MJ, Di Vona C, Jung C, Figueiredo AF, Luna J, Fillat C, Arbones ML, Graupera M, Valverde MA, de la Luna S. DYRK1A kinase positively regulates angiogenic responses in endothelial cells. Cell Rep. 2018;23(6):1867–78.
Liu H, Sun Q, Chen S, Chen L, Jia W, Zhao J, Sun X. DYRK1A activates NFATC1 to increase glioblastoma migration. Cancer Med. 2021;10:6416.
Jolly MK, Ware KE, Gilja S, Somarelli JA, Levine H. EMT and MET: necessary or permissive for metastasis? Mol Oncol. 2017;11(7):755–69.
Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson A, Kampf C, Sjostedt E, Asplund A, Olsson I, Edlund K, Lundberg E, Navani S, Szigyarto CA, Odeberg J, Djureinovic D, Takanen JO, Hober S, Alm T, Edqvist PH, Berling H, Tegel H, Mulder J, Rockberg J, Nilsson P, Schwenk JM, Hamsten M, von Feilitzen K, Forsberg M, Persson L, Johansson F, Zwahlen M, von Heijne G, Nielsen J, Ponten F. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.
Uhlen M, Zhang C, Lee S, Sjostedt E, Fagerberg L, Bidkhori G, Benfeitas R, Arif M, Liu Z, Edfors F, Sanli K, von Feilitzen K, Oksvold P, Lundberg E, Hober S, Nilsson P, Mattsson J, Schwenk JM, Brunnstrom H, Glimelius B, Sjoblom T, Edqvist PH, Djureinovic D, Micke P, Lindskog C, Mardinoglu A, Ponten F. A pathology atlas of the human cancer transcriptome. Science. 2017;357(6352)
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, Varambally S. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.
Wei L, Lian B, Zhang Y, Li W, Gu J, He X, Xie L. Application of microRNA and mRNA expression profiling on prognostic biomarker discovery for hepatocellular carcinoma. BMC Genomics. 2014;15(Suppl 1):S13.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, Cerami E, Sander C, Schultz N. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, Antipin Y, Reva B, Goldberg AP, Sander C, Schultz N. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.
Menyhart O, Nagy A, Gyorffy B. Determining consistent prognostic biomarkers of overall survival and vascular invasion in hepatocellular carcinoma. Royal Society open science. 2018;5(12): 181006.
Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63.
Guo X, Zhang D, Zhang X, Jiang J, Xue P, Wu C, Zhang J, Jin G, Huang Z, Yang J, Zhu X, Liu W, Xu G, Cui Z, Bao G. Dyrk1A promotes the proliferation, migration and invasion of fibroblast-like synoviocytes in rheumatoid arthritis via down-regulating Spry2 and activating the ERK MAPK pathway. Tissue Cell. 2018;55:63–70.
Loh CY, Chai JY, Tang TF, Wong WF, Sethi G, Shanmugam MK, Chong PP, Looi CY. The E-cadherin and N-cadherin switch in epithelial-to-mesenchymal transition: signaling, therapeutic implications, and challenges. Cells. 2019;8(10):1118.
Li YL, Ding K, Hu X, Wu LW, Zhou DM, Rao MJ, Lin NM, Zhang C. DYRK1A inhibition suppresses STAT3/EGFR/Met signalling and sensitizes EGFR wild-type NSCLC cells to AZD9291. J Cell Mol Med. 2019;23(11):7427–37.
Guard SE, Poss ZC, Ebmeier CC, Pagratis M, Simpson H, Taatjes DJ, Old WM. The nuclear interactome of DYRK1A reveals a functional role in DNA damage repair. Sci Rep. 2019;9(1):6539.
Bardou P, Mariette J, Escudie F, Djemiel C, Klopp C. jvenn: an interactive Venn diagram viewer. BMC Bioinformatics. 2014;15:293.
Lanczky A, Gyorffy B. Web-based survival analysis tool tailored for medical research (KMplot): development and implementation. J Med Internet Res. 2021;23(7): e27633.
Thien A, Prentzell MT, Holzwarth B, Klasener K, Kuper I, Boehlke C, Sonntag AG, Ruf S, Maerz L, Nitschke R, Grellscheid SN, Reth M, Walz G, Baumeister R, Neumann-Haefelin E, Thedieck K. TSC1 activates TGF-beta-Smad2/3 signaling in growth arrest and epithelial-to-mesenchymal transition. Dev Cell. 2015;32(5):617–30.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.
Martin CE, Nguyen A, Kang MK, Kim RH, Park NH, Shin KH. DYRK1A is required for maintenance of cancer stemness, contributing to tumorigenic potential in oral/oropharyngeal squamous cell carcinoma. Exp Cell Res. 2021;405(1): 112656.
Liu Q, Liu N, Zang S, Liu H, Wang P, Ji C, Sun X. Tumor suppressor DYRK1A effects on proliferation and chemoresistance of AML cells by downregulating c-Myc. PLoS ONE. 2014;9(6): e98853.
Fernandez-Martinez P, Zahonero C, Sanchez-Gomez P. DYRK1A: the double-edged kinase as a protagonist in cell growth and tumorigenesis. Mol Cell Oncol. 2015;2(1): e970048.
Li Y, Zhou D, Xu S, Rao M, Zhang Z, Wu L, Zhang C, Lin N. DYRK1A suppression restrains Mcl-1 expression and sensitizes NSCLC cells to Bcl-2 inhibitors. Cancer Biol Med. 2020;17(2):387–400.
Chen JG, Chen WQ, Zhang SW, Zheng RS, Zhu J, Zhang YH. Incidence and mortality of liver cancer in China: an analysis on data from the National Registration System between 2003 and 2007. Zhonghua liu xing bing xue za zhi Zhonghua liuxingbingxue zazhi. 2012;33(6):547–53.
Lee Walmsley D, Murray JB, Dokurno P, Massey AJ, Benwell K, Fiumana A, Foloppe N, Ray S, Smith J, Surgenor AE, Edmonds T, Demarles D, Burbridge M, Cruzalegui F, Kotschy A, Hubbard RE. Fragment-derived selective inhibitors of dual-specificity kinases DYRK1A and DYRK1B. J Med Chem. 2021;64(13):8971–91.
Li L, Wei JR, Song Y, Fang S, Du Y, Li Z, Zeng TT, Zhu YH, Li Y, Guan XY. TROAP switches DYRK1 activity to drive hepatocellular carcinoma progression. Cell Death Dis. 2021;12(1):125.
Bhansali RS, Rammohan M, Lee P, Laurent AP, Wen Q, Suraneni P, Yip BH, Tsai YC, Jenni S, Bornhauser B, Siret A, Fruit C, Pacheco-Benichou A, Harris E, Besson T, Thompson BJ, Goo YA, Hijiya N, Vilenchik M, Izraeli S, Bourquin JP, Malinge S and Crispino JD. DYRK1A regulates B cell acute lymphoblastic leukemia through phosphorylation of FOXO1 and STAT3. J Clin Invest. 2021;131(1)
Li YL, Wu LW, Zeng LH, Zhang ZY, Wang W, Zhang C, Lin NM. ApoC1 promotes the metastasis of clear cell renal cell carcinoma via activation of STAT3. Oncogene. 2020;39(39):6203–17.
Alvarez MA, Freitas JP, Mazher Hussain S, Glazer ES. TGF-beta inhibitors in metastatic pancreatic ductal adenocarcinoma. J Gastrointest Cancer. 2019;50(2):207–13.
Wang G, Yu Y, Sun C, Liu T, Liang T, Zhan L, Lin X, Feng XH. STAT3 selectively interacts with Smad3 to antagonize TGF-beta signalling. Oncogene. 2016;35(33):4388–98.
Altomare DA, Testa JR. Perturbations of the AKT signaling pathway in human cancer. Oncogene. 2005;24(50):7455–64.
Kurebayashi Y, Nagai S, Ikejiri A, Koyasu S. Recent advances in understanding the molecular mechanisms of the development and function of Th17 cells. Genes Cells. 2013;18(4):247–65.
Becker W, Sippl W. Activation, regulation, and inhibition of DYRK1A. FEBS J. 2011;278(2):246–56.
Li Z, Kong Y, Song L, Luo Q, Liu J, Shao C, Hou X, Liu X. Plk1-mediated phosphorylation of TSC1 enhances the efficacy of rapamycin. Can Res. 2018;78(11):2864–75.
The study was funded by National Natural Science Foundation of China [81702887, 2017], Zhejiang Provincial Natural Science Foundation of China [LTY21H160001, LY21H160017, 2021], Huadong Medicine Joint Funds of the Zhejiang Provincial Natural Science Foundation of China [HDMY22H160421, 2022], Public-service Technology Research Plan of Zhejiang Province [LGF21H310002, 2021], Zhejiang Provincial Medical and Health Technology Project [2021433724, 2021; 2020RC026, 2020], Key Laboratory of Clinical Cancer Pharmacology and Toxicology Research of Zhejiang Province [2020E10021, 2020], Zhejiang Provincial Program for the Cultivation of High-level Innovative Health Talents [2020-18, 2020].
Ethics approval and consent to participate
The experiments were performed in compliance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals.
Consent for publication
The authors have declared that no competing interest exists.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
DYRK1A might be involved in the metastasis of HCC cells. (a) Highly invasive HepG2 cells were established by serial selection via a Transwell assay. (b) The effect of DYRK1A on the progression-free survival and progression/metastasis-free survival of liver cancer patients. The data were collected from Kaplan–Meier Plotter (http://kmplot.com/analysis/index.php?p=background). Gene: DYRK1A; Survival: PFS, n = 370 (left panel); Survival: PFS/RFS, Vascular invasion: none, n = 205 (right panel) (c) Data collected from the LinkedOmics platform (www.linkedomics.org/admin.php) are shown. Sample cohort: TCGA_LIHC; Institute: UNC; Data type: RNAseq; Platform: HiSeq RNA; Attribute: DYRK1A; Statistical methods: Pearson correlation test; Patients: 371; Tool: overrepresentation enrichment analysis; Functional database: Gene ontology analysis (biological process, cellular component and molecular function).
Expression of DYRK1A in HCC samples. Compared with that in (a) normal liver tissues, DYRK1A expression was upregulated in (b) HCC tissues.
DYRK1A suppression failed to suppress cell proliferation. (a) HCC cells were incubated with the DYRK1A overexpression plasmid or empty vector for 48 h, and the expression of DYRK1A was detected (upper panel). HCC cells were transfected with DYRK1A plasmid or empty vector for 24 h at 6-well plates, then transferred to 96-well plated for the indicated times, and finally SRB assay was performed (lower panel). (b) HCC cells were incubated with DYRK1A siRNA or control siRNA for 24 h at 6-well plates, then transferred to 96-well plated for the indicated times, and finally SRB assay was performed. (c) HCC cells were incubated with harmine, and an SRB assay was performed.
TSC1 was involved in DYRK1A-mediated promotion of metastasis. (a) The online analysis tool HCMDB was used to investigate the expression of 25 overlapping genes in the indicated tissues (http://hcmdb.i-sanger.com/index). (b) Prognostic value of TSC1, DCAF7 and EPS15 in patients with liver cancer (http://kmplot.com/analysis/index.php?p=background).
TSC1 knockdown downregulated SMAD2/3 in HCC cells. (a) HCC cells were transfected with TSC1 siRNA or control small interfering RNA for 48 h, and the expression levels of the indicated proteins were assessed by western blotting. (b) HCC cells were incubated with the STAT3 overexpression plasmid or STAT3 siRNA for 48 h, and western blot analysis was then performed.
RNA sequencing analysis of gene expression in highly invasive and weakly invasive HepG2 cells. RNA sequencing was performed to compare gene transcription levels between P5 and P0.
Overrepresentation enrichment analysis of DYRK1A in clinical samples obtained from patients with hepatocellular carcinoma.
GSEA of hepatocellular carcinoma samples with low and high DYRK1A expression using the LinkedOmics platform by WikiPathways enrichment.
GSEA of hepatocellular carcinoma samples with low and high DYRK1A expression using the LinkedOmics platform by PANTHER pathways enrichment.
The gene lists of DYRK1A interacting and co-expressed proteins and SMAD2/3 positively co-expressed proteins.
About this article
Cite this article
Li, Yl., Zhang, Mm., Wu, Lw. et al. DYRK1A reinforces epithelial-mesenchymal transition and metastasis of hepatocellular carcinoma via cooperatively activating STAT3 and SMAD. J Biomed Sci 29, 34 (2022). https://doi.org/10.1186/s12929-022-00817-y
- Hepatocellular carcinoma