TAK1-mediated phosphorylation of PLCE1 represses PIP2 hydrolysis to impede esophageal squamous cancer metastasis

TAK1, a serine/threonine protein kinase, has been identified as a key regulator in a wide variety of cellular processes. However, its function and involved mechanism in cancer metastasis are still not well understood. Here, we found that knockdown of TAK1 promoted esophageal squamous cancer cell (ESCC) migration and invasion, whereas overexpression of TAK1 resulted in an opposite outcome. Moreover, these in vitro findings could be recapitulated in vivo in a xenograft metastasis mouse model. Mechanistically, co-immunoprecipitation combined with mass spectrometry demonstrated that TAK1 interacted with phospholipase C epsilon 1 (PLCE1), and phosphorylated PLCE1 at serine 1060 (S1060). Functional studies revealed that phosphorylation at S1060 in PLCE1 resulted in decreased enzyme activity, leading to a repression on PIP2 hydrolysis. As a result, the degradation products of PIP2 including DAG and inositol IP3 were reduced, which thereby suppressed signal transduction in the axis of PKC/GSK-3β/β-Catenin. Consequently, cancer metastasis related genes were impeded by TAK1. Overall, our data indicate that TAK1 plays a negative role in ESCC metastasis, which depends on TAK1 induced phosphorylation of PLCE1 at S1060.


Introduction
Esophageal cancer (EC) is the seventh most common malignancy worldwide with more than 600,000 new cases diagnosed annually (Sung, Ferlay et al., 2021).Clinical survey indicate that EC causes more than 500,000 deaths in 2020 globally (Sung et al., 2021), which is mainly due to its high aggressiveness and poor survival rate.The primary histological subtypes of EC include esophageal squamous cell carcinoma and esophageal adenocarcinoma (He, Xu et al., 2021, Lagergren, Smyth et al., 2017).Esophageal squamous cell carcinoma is a major subtype in China that account for more than 90% of EC (He et al., 2021).Due to the lack of specific biomarkers, the early diagnosis rate of ESCC is very low.Currently, esophageal squamous cell carcinoma diagnosis mainly relies on gastroscopy.By this way, most patients are found to be in locally advanced or metastatic stages when been diagnosed.The 5-year survival rate is approximately 15-25% (Pennathur, Gibson et al., 2013).Surgery has been increasingly used to treat esophageal squamous cell carcinoma at early stages.Neoadjuvant chemotherapy or chemoradiotherapy has become the standard treatment for locally advanced esophageal squamous cell carcinoma (He et al., 2021, Lagergren et al., 2017).In recent decades, therapeutic strategies have made great progress towards esophageal squamous cell carcinoma, and yet recurrence and metastasis are still significant challenges (Mao, Zeng et al., 2021).Therefore, it is urgent to deeply understand the mechanisms involved in cell metastasis, and thus try to explore new druggable targets for treating esophageal squamous cell carcinoma.
Transforming growth factor β-activated kinase 1 (TAK1), a member of the mitogen-activated protein kinase (MAPK) kinase kinase (MAP3K) family, is encoded by the gene Map3k7.As a serine/threonine protein kinase, TAK1 has been shown to play an integral role in inflammatory signal transduction across multiple pathways (Sakurai, 2012).By this mean, TAK1 was found to be involved in a wide variety of cellular processes including cell survival, cell migratory and invasive capacity, inflammation, immunity regulation, and tumorigenesis (Cho, Shim et al., 2021, Mukhopadhyay & Lee, 2020).Although TAK1 has been widely studied in cancer progression, its precise role remains controversial.For example, by acting as a tumor suppressor, TAK1 has been shown to negatively associated with tumor progression in several human cancers including prostate cancer (Huang, Tang et al., 2021), hepatocellular carcinoma (HCC) (Tan, Zhao et al., 2020, Wang, Zhang et al., 2021), cervical cancer (Guan, Lu et al., 2017) and certain blood cancers (Guo, Zhang et al., 2019).On the contrary, TAK1 has been shown to promote tumor progression across a range of human cancers, including colon cancer, ovarian, lung cancer and breast cancer (Augeri, Langenfeld et al., 2016, Cai, Shi et al., 2014, Xu, Niu et al., 2022).The above findings indicate that TAK1 may play a dual role in tumor initiation, progression, and metastasis.In our previous study, TAK1 has been shown to phosphorylate RASSF9 at serine 284 to inhibit cell proliferation by targeting the RAS/MEK/ERK axis in ESCC (Shi, Ju et al., 2021).Up to date, there is no report regarding the precise role of TAK1 on ESCC metastasis.
In this study, we examined the potential role of TAK1 in cell metastasis in ESCC.We found that TAK1 negatively regulates cell migration and invasion in ESCC, in which PLCE1 is a downstream target of TAK1.In detail, TAK1 phosphorylates PLCE1 at serine 1060 (S1060) to inhibit its enzyme activity, leading to decreased IP3 and DAG, both are the products of PLCE1 catalyzed reactions.As a result, IP3/DAG triggered signal transduction in the axis of PKC/GSK-3β/β-Catenin was blunted.All these outcomes induced by TAK1 caused a repression on cell migration and invasion in ESCC.

TAK1 negatively regulates ESCC migration and invasion
In our previous study, it has been shown that TAK1 expression was reduced in esophageal squamous tumor tissues, and TAK1 inhibits ESCC proliferation (Shi et al., 2021).To examine whether TAK1 affects the epithelial-mesenchymal transition (EMT) process in ESCC, we increased TAK1 expression in ECA-109 cells by transfecting a plasmid expressing Map3k7 (TAK1 gene name) and confirmed the overexpression of TAK1 (Fig. 1A).Due to the promising roles of epidermal growth factor (EGF) on the EMT (Lu, Ghosh et al., 2003), it was used to trigger the EMT in ESCC.As shown in Fig. 1B, EGF treatment induced spindle-shaped cell morphology in ECA-109 cells, which was markedly prevented by TAK1.These data imply that TAK1 is a negative regulator of the EMT in ESCC.To address this prediction, we performed transwell assay and found that TAK1 repressed cell migration and invasion in ECA-109 cells (Fig. 1C; SI Appendix, Fig. S1A).In addition, wound healing assay further confirmed the negative impact of TAK1 on cancer cell migration (Fig. 1D; SI Appendix, Fig. S1B).Next, we examined whether TAK1 affects the EMT related gene expression.Our data showed that TAK1 increased E-cadherin and ZO-1, two epithelial molecules, while mesenchymal molecules such as N-cadherin, Vimentin, and ZEB1 were reduced by TAK1 (Fig. 1E; SI Appendix, Fig. S1C).The qRT-PCR data further confirmed these changes (SI Appendix, Fig. S1D).
Next, we performed loss-of-function studies to examine the effects of TAK1 on cell migration and invasion in ESCC.The knockdown efficiency of Map3k7 siRNA on TAK1 expression was reported previously (Shi et al., 2021).Since the EMT process involves dynamic and spatial regulation of the cytoskeleton (Fife, McCarroll et al., 2014, Li & Wang, 2020), we therefore analyzed the expression of F-Actin by immunofluorescence.As shown in Fig. 1F, TAK1 knockdown induced F-Actin expression.Meanwhile, TAK1 knockdown promoted spindle-shaped mesenchymal morphology in ECA-109 cells (Fig. 1G).In accordance to these changes, the cell invasion and migration were stimulated as evidenced by the data from wound healing and transwell assays (Fig. 1H; SI Appendix, Fig. S2A-B).The epithelial markers including E-cadherin and Claudin-1 were decreased by TAK1 knockdown; while the mesenchymal markers such as N-cadherin, Vimentin, Snail, and Slug were increased (Fig. 1I; SI Appendix, Fig. S2C).At mRNA level, Cdh1 and Cldn1 were inhibited by TAK1 knockdown, while Cdh2 and Snail1 were promoted (SI Appendix, Fig. S2D).To verify TAK1 knockdown on cell migration and invasion, we also downregulated TAK1 expression by using lentivirus carrying Map3k7 shRNA (LV-Map3k7 shRNA) (SI Appendix, Fig. S3A, B).In LV-Map3k7 shRNA transduced cells, cell migration and invasion were obviously increased (SI Appendix, Fig. S3C-E).The epithelial markers including E-cadherin, ZO-1, and Claudin-1 were repressed by TAK1 knockdown; while the mesenchymal markers such as N-cadherin, Slug, and ZEB-1 were activated (SI Appendix, Fig. S3F, G).The qRT-PCR data showed similar changes in these EMT genes (SI Appendix, Fig. S3H).Additionally, TAK1 knockdown was accomplished by CRISPR-Cas9 using Map3k7 guide RNA (gRNA) (SI Appendix, Fig. S4A).Again, we observed cell migration and invasion were enhanced by Map3k7 gRNA (SI Appendix, Fig. S4B-D).The expressions in E-cadherin, ZO-1, and Claudin-1 were inhibited by Map3k7 gRNA, while the expressions in N-cadherin, Vimentin, Snail, and Slug were activated (SI Appendix, Fig. S4E, F).Similar changes in these EMT related gene expression were also observed in qRT-PCR analysis (SI Appendix, Fig. S4G).Furthermore, we inhibited TAK1 activity by Oxo, NG25, or Takinib and observed that all these treatments induced morphological changes in ECA-109 cells from round shape to spindle-like shape (SI Appendix, Fig. S5A).Cell migration and invasion were also stimulated by TAK1 inhibition (SI Appendix, Fig. S5B-E).Taken together, these findings indicate that TAK1 is a negative regulator of cell migration and invasion in ESCC.

TAK1 phosphorylates PLCE1 at serine 1060
The above data clearly showed that TAK1 negatively regulates cell migration and invasion in ESCC.To reveal the involved mechanism, we performed co-immunoprecipitation combined with mass spectrometry to identify the potential downstream target of TAK1.As previously reported, 24 proteins were found to be phosphorylated in the immunocomplex (Shi et al., 2021).Of them, phospholipase C epsilon 1 (PLCE1) caused our attention due to its essential roles in cell growth, migration, and metastasis in various human cancers (Abnet et al., 2010, Chen, Wang et al., 2019, Chen, Xin et al., 2020, Gu, Zheng et al., 2018, Kadamur & Ross, 2013, WangZhou et al., 2010).According to the mass spectrometry data, the serine residue at 1060 (S1060) in PLCE1 was phosphorylated (Fig. 2A).Currently, there is no antibody against phosphorylated PLCE1 at S1060 (p-PLCE1 S1060), we therefore generated an antibody for detecting p-PLCE1 S1060.To verify the antibody specificity, S1060 in PLCE1 was mutated into alanine (S1060A).As shown in Fig. 2B, wild type (WT) PLCE1 could be phosphorylated by TAK1 at S1060, while PCLE1 S1060A was resistant to TAK1 induced phosphorylation at S1060.Moreover, p-PLCE1 S1060 induced by TAK1 could be weakened by (5Z)-7-Oxozeaenol (Oxo), a potent inhibitor of TAK1 (Fig. 2C).Of note, phosphorylated TAK1 (p-TAK1) was markedly decreased in the presence of Oxo, indicating TAK1 was inhibited (Fig. 2C).To further confirm these findings, other TAK1 inhibitors such as Takinib and NG25 were employed, and similar changes in p-PLCE1 and p-TAK1 were observed (Fig. 2D, E).
Collectively, these results indicate that PLCE1 is a downstream target of TAK1, and PLCE1 S1060 was phosphorylated by TAK1.To further confirm this notion, we analyzed TAK1 and p-PCLE1 expression in clinical samples.As shown Fig. 2F, TAK1 expression was reduced in tumor tissues as compared to their respective adjacent normal tissues.In accordance to these changes in TAK1, p-PLCE1 was also decreased in tumor tissues (Fig. 2G).The Pearson correlation tests showing TAK1 expression was positively correlated with p-PLCE1 (Fig. 2H).These changes in TAK1 and p-PLCE1 were recapitulated by western blot data (Fig.

TAK1 phosphorylates PLCE1 to inhibit cell migration and invasion in ESCC
It has been well documented that PLCE1 plays a key role in cancer progression (Abnet et al., 2010, Chen et al., 2019, Chen et al., 2020, Gu et al., 2018, Wang et al., 2010).Therefore, we predicted that the negative impact of TAK1 on cell migration and invasion in ESCC may rely on PLCE1.To verify this prediction, we first examined whether PLCE1 affects cell migration and invasion.In ECA-109 cells, PLCE1 overexpression significantly increased cell migration and invasion (Fig. 3A-C; SI Appendix, Fig. S6A, B).In accordance to these changes, the epithelial marker E-cadherin was repressed by PLCE1, while the mesenchymal markers including Vimentin, Snail, Slug, and ZEB-1 were activated (Fig. 3D; SI Appendix, Fig. S6C).The qRT-PCR data also revealed that Cdh1 was downregulated by PLCE1, while Vim, Snail1, and Snail2 were upregulated (SI Appendix, Fig. S6D).To further confirm these findings, we performed PLCE1 knockdown by using siRNA.As shown in Fig. 3E, Plce1 expression was markedly reduced by all three tested siRNAs.Similar changes were also observed in the western blot data (Fig. 3F).Of these siRNAs, siRNA-2# exhibited the best knockdown efficiency and it was chosen in the following experiments.On the contrary, as compared to PLCE1 overexpression, PLCE1 knockdown impeded cell migration and invasion (Fig. 3G, H

S9A-D).
Due to PLCE1 is a lipid hydrolase, we next asked whether TAK1-induced phosphorylation in PLCE1 affects its enzyme activity.To this aim, we transfected ECA-109 cells with a plasmid bearing Myc-PLCE1 or TAK1, and PLCE1 was captured by pulldown using anti-Myc beads, which was then subjected to enzyme activity assay.As shown in Fig. 4A, the PLCE1 activity was reduced in the presence of TAK1; however, the PLCE1 S1060A activity was not affected (Fig. 4B), indicating TAK1 induced phosphorylation at S1060 repressed PLCE1 activity.As a lipid hydrolase, PLCE1 catalyzes PIP2 hydrolysis to produce IP3 and DAG, both are the second messengers in diverse cellular events (Kadamur & Ross, 2013).Therefore, we examined these two products to verify the inhibitory effect of TAK1 on PLCE1.Our data showed that the productions of IP3 and DAG were increased by PLCE1, which could be largely counteracted by TAK1 (Fig. 4C-E).As a messenger, IP3 can induce endoplasmic reticulum to release Ca 2+ into the cytoplasm (Kadamur & Ross, 2013).Hence, we analyzed the cytoplasmic Ca 2+ by Fluo-4 AM staining.As expected, the signals for the cytoplasmic Ca were obviously increased by PLCE1 in ECA-109 cells, and this trend was weakened by TAK1 (Fig. 4F, G).Meanwhile, we also directly detected the cytoplasmic Ca 2+ with a fluorospectrophotometer and similar changes were observed (Fig. 4H).Flow cytometry analysis also evidenced the increase in the cytoplasmic Ca 2+ induced by PLCE1 was prevented by TAK1 (Fig. 4I).All these evidences indicate that TAK1 inhibits PLCE1 activity, which is likely due to TAK1 induced phosphorylation in PLCE1 at S1060.To further verify this notion, a series of functional experiments were performed.For instance, we observed that TAK1 could reverse PLCE1-induced cell migration and invasion; however, TAK1 failed to affect PLCE1 S1060A-induced cell migration and invasion (SI Appendix, Fig. S10A-D).Overall, these data clearly indicate that PLCE1 has a positive impact on cell migration and invasion in ESCC.By inducing phosphorylation at S1060, TAK1 inhibits PLCE1 enzyme activity and thereby counteracts PLCE1-induced cell migration and invasion in ESCC.At this regard, PLCE1 is a downstream target of TAK1 for transducing its inhibitory effects on cell migration and invasion in ESCC.

TAK1 inhibits PLCE1-induced signal transduction in the PKC/GSK-3β/β-Catenin axis
PLCE1 hydrolyzes PIP2 to produce two important second messengers DAG and IP3 that trigger Ca 2+ release from the endoplasmic reticulum to the cytoplasm to activate PKC (Harden, Hicks et al., 2009, Kadamur & Ross, 2013), suggesting PKC activation is responsible for PLCE1-induced cell migration and invasion in ESCC.To address this prediction, we treated cells with 2-APB, an IP3 receptor (IP3R) inhibitor, and found that 2-APB treatment almost completely reversed PLCE1-induced the intracellular Ca 2+ (SI Appendix, Fig. S11A-E).
We therefore asked whether PLCE1-induced cell growth behaviors are depended on PKC activation.To this end, Midostaurin was used to inhibit PKC and our data showed that Midostaurin almost completely abolished cell migration and invasion induced by PLCE1 (SI Appendix, Fig. S16A-D).
Therefore, we predicted that PLCE1 stimulates cell migration and invasion via the axis of PKC/GSK-3β/β-Catenin. To test this prediction, ECA-109 cells were treated with IP3 receptor

TAK1 inhibition promotes ESCC metastasis in vivo
To examine TAK1-regulated ESCC metastasis in vivo, a xenograft model with nude mice was employed.Mice was injected with ECA109 cells via the tail vein, and then mice were daily treated with Takinib (50 mg/kg) or corn oil (Vehicle) for consecutive 15 days.Eight weeks later, mice were sacrificed for analyzing cancer cell metastasis.We found 4 of 6 mice in the Takinib group developed lung metastasis, while only 1 of 6 mice in the Vehicle group had lung metastasis.Moreover, the number of metastatic nodules in lung tissues was higher in Takinib treated mice (Fig. 6A-C).In addition, the effects of TAK1 inhibition on the axis of PKC/GSK-3β/β-Catenin was examined.As shown in Fig. 6D, TAK1 inhibition by Takinib activated PKC, leading to increased p-GSK-3.As a result, p--Catenin was reduced and MMP2 expression was upregulated (Fig. 6D, E).Of note, p-TAK1 was decreased by Takinib, indicating TAK1 activity was successfully inhibited (Fig. 6D, E).In accordance to this change, p-PLCE1 was decreased in Takinib-treated mice (Fig. 6D, E).Hence, these data indicate that TAK1 represses ESCC metastasis in vivo, and this benefit is likely due to TAK1-mediated phosphorylation in PLCE1 at S1060.

PLCE1 facilitates ESCC metastasis in vivo
To further verify the function of PLCE1 in ESCC metastasis, we generated a stable cell line with PLCE1 knockdown by using lentivirus bearing PLCE1 shRNA (LV-shPLCE1) (SI Appendix, Fig. S18A-B).These cells with low PLCE1 expression were injected into nude mice via the tail vein to produce a mouse xenograft tumor model.Eight weeks later, mice were sacrificed, the nodule number and incidence rate were reduced in mice injected with LV-shPLCE1 cells (Fig. 7A-C).In addition, western blot analysis showed that PLCE1 silencing inhibited signal transduction in the PKC/GSK-3β/β-Catenin axis, as evidenced by reduced p-PKC and p-GSK-3β, and increased p-β-Catenin (Fig. 7D, E).In accordance to these changes, the expression of MMP2 was decreased (Fig. 7D, E).These results further confirmed the notion that PLCE1 plays a key role in cancer cell metastasis in ESCC.
Furthermore, TAK1 represses the transcription of human telomerase; and activate tumor suppressor protein LKB1, all these evidences indicate that TAK1 is a tumor suppressor (Adhikari, Xu et al., 2007, Fujiki, Miura et al., 2007, Xie, Zhang et al., 2006).Consistent with these findings, in our previous study, we found that TAK1 expression was reduced in esophageal squamous tumors compared to adjacent normal tissues, and TAK1 negatively correlates with esophageal squamous tumor patient survival (Shi et al., 2021).In this study, we extended our research on TAK1 and observed that TAK1 inhibits cell migration and invasion in ESCC.In contrast to these findings, TAK1 has been shown to play a positive role in tumor cell proliferation, migration, invasion, colony formation, and metastasis, especially in breast cancer, pancreatic cancer, and non-small cell lung cancer (Kim, Kim et al., 2023, Santoro, Zanotto et al., 2020, Tripathi, Shin et al., 2019).Therefore, TAK1 may play pleiotropic roles in different types of cancer cells.
As a kinase, TAK1 is considered as a central regulator of tumor cell proliferation, migration and invasion and has been demonstrated to play tumor suppression or activation role via the NF-κB and MAPK activation (Mukhopadhyay & Lee, 2020, Wang et al., 2021, Zhang, Cheng et al., 2023).However, the precise molecular and cellular mechanisms by which TAK1 regulates ESCC metastasis remain unclear.As previously reported, by co-immunoprecipitation coupled with mass spectrometry (MS/MS), we identified RASSF9 as a downstream target of TAK1 for transducing its repression on cell proliferation in ESCC (Shi et al., 2021).By phosphorylating RASSF9 at S284, TAK1 negatively regulates cell proliferation in ESCC (Shi et al., 2021).
Meanwhile, we also found that PLCE1 was presented in the immunocomplex, which was phosphorylated at S1060 in the presence of TAK1 in ECA-109 cells.As a member of the human phosphoinositide-specific PLC family, PLCE1 can be activated by a variety of intracellular and extracellular signal molecules including hormones, cytokines, neurotransmitters, growth factors, etc. (Kadamur & Ross, 2013).Numerous studies have shown that PLCE1 plays a key role in cancer development and progression through various pathways (Abnet et al., 2010, Fan et al., 2019, Ghosh, Nataraj et al., 2021, He, Wang et al., 2016, Yue et al., 2019).Moreover, PLCE1 has been characterized as a susceptibility gene in ESCC (Abnet et al., 2010, Wang et al., 2010).All these evidences suggest that PLCE1 is a potential downstream target for transducing the negative impact of TAK1 on cell migration and invasion in ESCC.
Once activated, PLCE1 catalyzes the hydrolysis of PIP2  Smrcka, Brown et al., 2012).Consistent with these notions, in this study, we observed that PLCE1-overexpressing cells exhibit higher production of intracellular IP3, DAG, and intracellular Ca 2+ .As a consequence, we found that PLCE1 overexpression accelerates cell migration and invasion, while PLCE1 knockdown results in an opposite phenotype.Moreover, in the presence of TAK1, PLCE1-induced cell migration and invasion were largely counteracted.By considering TAK1 is a protein kinase, therefore, we proposed that TAK1 phosphorylates PLCE1 at S1060 to inhibit its enzyme activity and thus impedes cell migration and invasion in ESCC.Indeed, PLCE1 activity was blunted in the presence of TAK1, whereas mutated PLCE1 (PLCE1 S1060A) was not affected.It has been shown that PKC positively regulates the expression and stability of β-Catenin in a GSK-3β-dependent manner (Duong et al., 2017, Gwak et al., 2006, Ryu & Han, 2015).Therefore, we investigated the impact of PLCE1 on the PKC/GSK-3β/β-Catenin axis.As expected, our data revealed that PLCE1 promotes cell migration and invasion in ESCC by activating the PKC/GSK-3β/β-Catenin pathway, results in a series of phosphorylation modification in PKC and GSK-3β.As a result, phosphorylated β-Catenin decreased, leading to higher stability of β-Catenin, which eventually promotes the EMT process by affecting epithelial and mesenchymal gene expression.
In summary, our findings indicate TAK1 is a negative regulator in cell migration and invasion in ESCC.Mechanistically, TAK1 phosphorylates PLCE1 at S1060 to inhibit its phospholipase activity, leading to reductions in IP3 and DAG.Consequently, PKC activity was blunted, which results in decreased phosphorylated GSK-3β and increased phosphorylated β-Catenin.As a consequence, the stability of β-Catenin was decreased and its transcriptional activity was blunted.By this mean, the epithelial marker gene expression was upregulated by TAK1, while the mesenchymal marker gene expression was downregulated.All these outcomes eventually cause a breakdown in cell migration and invasion in ESCC.Hence, our data revealed a new facet of TAK1 in the EMT process in ESCC by inhibiting PLCE1 activity and its downstream signal transduction in the axis of PKC/GSK-3/-Catenin.Moreover, TAK1, together with its downstream target PLCE1, are potential druggable targets for developing agents for dealing with ESCC.

Human esophageal squamous tumor specimens
Human esophageal squamous tumor specimens were obtained from the Affiliated Hospital of Nantong University.Human sample collection for research was conducted in accordance with the recognized ethical guideline of Declaration of Helsinki and approved by the Ethics Committee of Affiliated Hospital of Nantong University.

Generation of antibodies against phospho-PLCE1 (S1060)
Rabbit polyclonal antibody that recognizes phospho-PLCE1 (S1060) was raised against c-WSARNPS(p)PGTSAK peptide at Absin (Shanghai, China).Briefly, the phosphorylated polypeptide c-WSARNPS(p)PGTSAK was synthesized as the target peptide.Then two rabbits were immunized with the target peptide antigen.Meanwhile, non-phosphorylated polypeptide c-WSARNPSPGTSAK was synthesized as a control.The rabbits were sacrificed and blood was collected for further purification.The prepared antibody was verified and stored at -20°C until use.

CRISPR-Cas9 mediated TAK1 deletion
Map3k7 was deleted by the method of CRISPR-Cas9 with the procedures describe previously (Shi et al., 2021).Briefly, we designed an integrated vector based on OriP/EBNA1 (epiCRISPR), which bears gRNA, Cas9 and puromycin resistance genes.ECA-109 cells were transfected with the plasmids using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions.The cells were cultured in a medium containing 2.5 μg/ml puromycin for 1 week.The surviving stable cells were used for further experiments.The gRNA sequence is listed in Table S1.

Protein extraction and western blot analysis
Western blot assay was performed as previously described (Tan, Krueger et al., 2022).Briefly, Total protein were extracted from cells and tissues using protein lysis buffer supplemented with protease and phosphatase inhibitors (Roche Applied Science, Penzberg, Germany).Protein concentrations were assayed by the Pierce™ BCA Protein Assay Kit (Thermo Scientific, Waltham, MA, USA).Subsequently, equal amounts of protein were separated with 6% or 10% SDS-PAGE, and then protein were transferred onto 0.45 μm PVDF membranes.After blocked with 5% BSA in TBST for 1 h, the membranes were incubated with primary antibodies at 4°C overnight.Thereafter, the membranes were incubated with HRP-linked secondary antibodies for 1 h at room temperature.Protein signals were visualized using enhanced chemiluminescence (ECL) reagent (Thermo Fisher Scientific, Waltham, MA, USA) and quantitatively analyzed using ImageJ (NIH).Actin were used as a loading control.

2+
Cells were transfected with plasmid expressing PLCE1 and/or TAK1 for 6 h and then treated with BAPTA-AM (10 µM), 2-APB (10 µM), or Midostaurin (100 nM) for additional 18 h, respectively.Measurement of intracellular Ca 2+ in treated cells was performed using a Fluo-4 AM Assay Kit (Beyotime, Jiangsu, China).In detail, cells were washed with PBS, and incubated with 1 μM Fluo-4 AM in PBS at 37°C for 30 min.Cell were then washed with PBS for

PLCE1 pulldown and phospholipase activity assay
Protein pulldown was performed using a previously reported method (Shi et al., 2021).Briefly, HEK293 cells were transfected with a plasmid expressing PLCE1-Myc or TAK1.24 h post-transfection, cells were harvested for preparing total cell lysates, which was then subjected to protein pulldown assay using the Myc-Tag beads (Sepharose® Bead Conjugate; #55464; CST, Beverley, MA, USA).After 4 h-incubation on a rotator, unbound proteins were removed using the ice-cold wash buffer.Finally, phospholipase activity bound on the beads was detected with a PLC Activity Assay Kit (Jining Shiye, Shanghai, China) with the provided protocol.In rationale, PLC catalyzes the hydrolysis of O-(4-nitrophenyl) choline (NPPC) to produce p-nitrophenol (PNP), which has an absorption maximum at 405 nm.By detecting the increase rate of PNP at 405 nm, the magnitude of PLC enzyme activity can be calculated.

IP3 and DAG assays
IP3 and DAG levels in cells were assayed by a Human Inositol 1,4,5-triphosphate enzyme-linked Immunosorbent Assay (ELISA) Kit and a Human Diacylglycerol commercial ELISA Kit (Mlbio, Shanghai, China) according to the manufacturer's instructions.The absorbance was measured at 450 nm with a microplate reader.

Co-immunoprecipitation and MS/MS spectrometry
Co-immunoprecipitation and MS/MS assay were performed as previously described (Shi et al., 2021).Briefly, ECA-109 cells were transfected with a plasmid expressing Map3k7.36 h post-transfection, cells were harvested for co-immunoprecipitation by using an antibody against TAK1.The resulting immunocomplex was subjected to LC-MS/MS spectrometry analysis (Shanghai Applied Protein Technology Co., Ltd.).MS/MS spectra were searched using MASCOT engine (Matrix Science, London, UK; version 2.4) against the UniProKB human (161584 total entries, downloaded 20180105).

Mouse xenograft metastasis model
For in vivo metastasis assay, 4-week-old male BALB/c nude mice (18-20 g) were purchased from Shanghai Slake Laboratory Animal Co. Ltd. (Shanghai, China) and randomly divided into 2 groups (n = 6 -10 per group).Mice were housed in air-filtered laminar flow cages (5 mice/cage) with 12 h light cycle and adequate food and water.To examine the role of PLCE1 on metastasis, each mouse was intravenously (i.v.) injected with 1 x 10 6 ECA109 cells diluted in 100 µl PBS, which were transduced with LV-shPLCE1 or negative control virus.Additionally, to examine TAK1 inhibition on metastasis, each mouse was intravenously (i.v.) injected with 1 x 10 6 ECA109 cells.Mice were then intraperitoneally injected with Takinib (50 mg/kg/day) or vehicle (corn oil) for 15 days.Mouse body weight was measured once a week.Eight weeks later, mice were sacrificed, and lungs and livers were collected for further analysis.All animal experiments were carried out in accordance with the institutional ethical guidelines of animal care, and were approved by the Animal Experimentation Ethics Committee of the Nantong University (Approval ID: SYXK (SU) 2017-0046).

Statistical analysis
Statistical significance of differences between groups were tested by using analysis of variance (ANOVA) or unpaired Student's t-test.Data were presented as mean ± SD from at least three independent experiments.All statistical analyses were performed using GraphPad Prism version 8.0 (GraphPad Software, La Jolla, CA, USA).* p < 0.05 was accepted as statistical significance.

3
times and further incubated for 20 min to ensure Fluo-4 AM completely transformed into Fluo-4 in cells.Finally, Fluo-4 fluorescent intensity was detected with a confocal laser scanning fluorescence microscope (Olympus BX51, Tokyo, Japan), fluorescence microplate (BioTek Instruments, Inc.), or flow cytometry (FACSCalibur, BD Bioscience, Franklin Lakes, NJ, USA) at excitation of 488 nm and emission of 520 nm to determine the change of intracellular Ca 2+ concentration.

Figure legends Figure 1 .
Figure legends

Figure 3 .
Figure 3. PLCE1 positively regulates ESCC migration and invasion (A) Increased expression of PLCE1 in ECA-109 cells transfected with a plasmid expressing Plce1.(B-C) Increased expression of PLCE1 enhances cell migration and invasion.ECA-109 cells were transfected with the plasmid carrying Plce1.Twenty-four h-post transfection, cells were subjected to transwell (B) or wound healing (C) assay.n = 5 biologically independent replicates.(D) Increased expression of PLCE1 in ECA-109 cells induces mesenchymal protein marker