- Open Access
PLCε regulates prostate cancer mitochondrial oxidative metabolism and migration via upregulation of Twist1
Journal of Experimental & Clinical Cancer Researchvolume 38, Article number: 337 (2019)
Metabolic rewiring is a common feature of many cancer types, including prostate cancer (PCa). Alterations in master genes lead to mitochondrial metabolic rewiring and provide an appealing target to inhibit cancer progression and improve survival. Phospholipase C (PLC)ε is a regulator of tumor generation and progression. However, its role in mitochondrial metabolism remains unclear.
The GEO, The Cancer Genome Atlas, and the GTEx databases were used to determine Twist1 mRNA levels in tumors and their non-tumor counterparts. Fifty-five PCa and 48 benign prostatic hypertrophy tissue samples were tested for the presence of PLCε and Twist1 immunohistochemically. An association between PLCε and Twist1 was determined by Pearson’s correlation analysis. PLCε was knocked down with a lentiviral short hairpin RNA. Mitochondrial activity was assessed by measuring the oxygen consumption rate. Western blotting analyses were used to measure levels of PPARβ, Twist1, phosphorylated (p)-Twist1, p-MEK, p-ERK, p-P38, and p-c-Jun N-terminal kinase (JNK). Cells were treated with inhibitors of MEK, JNK, and P38 MAPK, and an agonist and inhibitor of peroxisome proliferator activated receptor (PPAR) β, to evaluate which signaling pathways were involved in PLCε-mediated Twist1 expression. The stability of Twist1 was determined after blocking protein synthesis with cycloheximide. Reporter assays utilized E-cadherin or N-cadherin luciferase reporters under depletion of PLCε or Twist1. Transwell assays assessed cell migration. Finally, a nude mouse tumor xenograft assay was conducted to verify the role of PLCε in tumor formation.
Our findings revealed that the expression of PLCε was positively associated with Twist1 in clinical PCa samples. PLCε knockdown promoted mitochondrial oxidative metabolism in PCa cells. Mechanistically, PLCε increased phosphorylation of Twist1 and stabilized the Twist1 protein through MAPK signaling. The transcriptional activity of Twist1, and the Twist1-mediated epithelial-to-mesenchymal transition, cell migration, and transcription regulation, were suppressed by PLCε knockdown and by blocking PPARβ nuclear translocation. The tumor xenograft assay demonstrated that PLCε depletion diminished PCa cell tumorigenesis.
These findings reveal an undiscovered physiological role for PLCε in the suppression of mitochondrial oxidative metabolism that has significant implications for understanding PCa occurrence and migration.
Prostate cancer (PCa), a leading cause of cancer mortality in males worldwide,is associated with rapid treatment resistance, metastasis, and death . PCa metastasis is the main cause of death . Therefore, defining the molecular basis of metastasisis urgently required within a changing tumor microenvironment.
Cells continually balance their metabolism according to the availability of extracellular nutrients and the evolving microenvironment to meet the needs of their current state and function [3, 4]. Mitochondria, the metabolic cores of cells, play an essential role in this dynamic . Moreover, metabolic rewiring, a common feature of many cancer types, is an appealing target to inhibit cancer progression and improve survival .
Numerous reports indicate that the main trigger for metabolic reprogramming is an alteration in cancer genes that remodels the signaling pathways involving one or a few enzymes . Peroxisome proliferator-activated receptor (PPAR) gamma co-activator-1 alpha (PGC-1α) was first identified as a PPARγ-interacting protein in brown fat . Expression of PGC-1α powerfully induces mitochondrial oxidative metabolism and uncoupling proteins, and it can interact with a diverse range of transcription factors to control metabolic functions [8, 9]. Although, mitochondrial functions modulated by PGC-1α positively impact the in-vivo growth of breast cancer cell lines , human tumor transcriptional data, cell culture, and in-vivo mouse experiments revealed that a PGC1α-modulated transcriptional program promotes oxidative metabolism and a universal catabolic state to restrain PCa growth and metastasis . Another report also demonstrated low oxygen consumption of prostatic tissue in general . Overall, investigations of mitochondrial metabolism in PCa are increasing.
Phospholipase C (PLC)ε, belongings to the phosphoinositide-specific PLC family, is a multifunctional signaling protein and combines both PLC and guanine nucleotide exchange factor activities [13, 14]. PLC isozymes catalyze the cleavage of polyphosphoinositides, such as phosphatidylinositol 4,5-bisphosphate, into two intracellular second messengers, diacylglycerol and inositol 1,4,5-trisphosphate, that regulate protein kinase C activity and calcium mobilization, respectively . Importantly, specific to the Ras-associating and CDC25A domains, PLCε leads to Rap1phosphorylation and B-raf activation, causing an enzyme-induced cascade reaction of mitogen-activated protein kinases (MAPKs) .
Even several studies that analyzed human colorectal tumors showed markedly reduced levels of PLCε [17,18,19]. An increasing number of studies have demonstrated that PLCε is involved in the generation and progression of tumors [20,21,22], and a single nucleotide polymorphism of esophageal squamous cell carcinoma in the PLCε gene indicates that it potentially serves as an oncogene leading to cellular transformation . However, the functions of PLCε in mitochondrial metabolismin cancer are not yet fully established.
Twist1, a highly conserved transcription factor that contains a basic helix-loop-helixdomain was initially recognized as a mesoderm-inducing factor in Drosophila , and it is well-known as a major inducer of epithelial-to-mesenchymal transition (EMT) . Aberrant Twist1 overexpression facilitates EMT, cell motility, and invasive activity, and enhances some features of cancer stem cells via control of downstream gene expression [26, 27]. Recent studies have highlighted that the increase of Twist1 expression in patients is associated with the progression of several human malignancies and poor survival [25, 28, 29].
Accumulating evidence demonstrates that Twist1 exerts multiple activities during carcinogenesis . Notably, Pan et al.found that the nuclear receptor PPARβ mediates the actions of PGC-1α and suppresses mitochondrial metabolism and uncoupling via Twist1. This revealed an unexpected physiological role of Twist1in the maintenance of energy homeostasis and control of metabolism . Furthermore, to our knowledge, the exact mechanisms of Twist1 in modulating tumor metabolism remain elusive.
In this study, we investigated the role of PLCε in PCa metabolic control and attempted to find potential connections between PLCε and Twist1. Our hypothesis was that PLCε knockdown would accelerate PGC-1α-mediated mitochondrial metabolism by modulating Twist1.
Patient and specimens
A total of 48 cases of benign prostatic hypertrophy (BPH) samples and 55 cases of PCa samples were collected from patients at Department of Urinary Surgery, the First Affiliated Hospital of Chongqing Medical University, between 2015 and 2017. All of the tissue samples were diagnosed to be PCa or BPH by histological identification and stored in liquid nitrogen before the experiments. Meanwhile another 97 patient serum samples were collected between 2017 and 2018. All patients shared a similar socioeconomic status and provided informed consent.
The tissue specimens were embedded in paraffin overnight and were then sliced into 5-mm-thick sections. Immunohistochemical (IHC) analysis of prostate tissues were performed using antibodies against PLCε (santa cruz) and Twist1 (Absin). The protocols are the same as we published before .
Human prostate cancer cells (PC3, DU145, LNCap) are purchased by two agent companies from the Cell Resource Center of Shanghai Academy of Life Sciences, Chinese Academy of Sciences. The PC3 and LNCap cells were purchased from Shanghai Zhong Qiao Xin Zhou Biotechnology Company. The DU145 cells were purchased from Shanghai Biowing Applied Biotechnology Company. Bicalutamide-resistant (Bica-R) cells was constructed early in our group . The LNCap cells, one of the androgen-dependent prostate cancer cell strains, were treated with bicalutamide (10 μM) (Selleck Chemicals, Houston, TX, USA), respectively for at least 6 months to generate bicalutamide-resistant cells. Cells were cultured in RPMI-1640 (Gibco, USA) medium supplemented with 10% fetal bovine serum (Gibco, USA), 1% penicilin/streptomycin ((Life Technologies, Carlsbad CA) in a humidified incubator of 5% CO2 at 37 °C.
Plasmid and lentiviral vector
Twist1 (Myc-DDK-tagged) vector was buyed from OriGene Technologies (RC202920). 7TFP CDH1 (E-cadherin) reporter was purchased from Addgene (http://www.addgene.org). Twist1 reporter (S717559) was purchased from Active Motif. the N-cadherin reporter was generated as described previously . The first intron (+ 746 to + 3156 bp) from the human N-cadherin gene was obtained by PCR using primers 5′- AAATTTGAGCTCGGCTCTAGGGGCTGGATT-3′ (forward) and 5′-GGTTGGAGATCTTGTTGTTCGGGCGTGTAA-3′(reverse). The sequences of LV-sh-NC, LV-sh-PLCε#1, LV-sh-PLCε#2, LV-sh-PLCε#3, vetor-sh-NC, vetor-sh-Twist1#1, vetor-sh-Twist2#, vetor-sh-Twist1#3 and vetor-sh-PPARβ#1 (all transcripts from Gene Pharma Company, Shanghai, China) are showed in the following: LV-sh-NC-F, TTCTCCGAACGTGTCACGT; LV-sh-NC-R, AAGAGGCTTGCACAGTGCA. LV-sh-PLCε#1-F, GCAGATATCTGATGCCATTGC; LV-sh-PLCε#1-R, CGTCTATAGACTACGGTCCCG; LV-sh-PLCε#2-F; GCTTCTTAACACGGGACTTGG.LV-sh-PLCε#2-R, CGAAGAATTGTGACCTGAACC. LV-sh-PLCε#3-F, GGTTCTCTCCTAGAAGCAACC; LV-sh-PLCε#3-R, CCAAGAGAGGATCTTCGTTGG. vetor-sh-NC-F, CACCGTTCTCCGAACGTGTCACGTTTCAAGAGAACGTGA CACGTTCGGAGAATTTTTTG; vetor-sh-NC-R, GATCCAAAAAATTCTCCGAACGTGTCACGTTCTCTTGAAACGTGACACGTTCGGAGAAC. vetor-sh-Twist1#1-F, CACCGAAGCTGAGCAAGATTCAGACCTTCAAGAGAGGTCTGAATCTTGCTCAGCTTTTTTTTG; vetor-sh-Twist1#1-R, GATCCAAAAAAAAGCTGAGCAAGATTCAGACCTCTCTTGAAGGTCTGAATCTTGATCAGCTTC.vetor-sh-Twist1#2-F, CCGGTGCTGGACTCCAAGATGGCAAGCTCGAGCTTGCCATCTTGGAGTCCAGCTTTTTG; vetor-sh-Twist1#2-R, AATTCAAAAAGCTGGACTCCAAGATGGCAAGCTCGAGCTTGCCATCTTGGAGTCCAGCA. vetor-sh-Twist1#3-F, CCGGTCCTGAGCAACAGCGAGGAAGACTCGAGTCTTCCTCGCTGTTGCTCAGGTTTTTG; vetor-sh-Twist1#3-R, AATTCAAAAACCTGAGCAACAGCGAGGAAGACTCGAGTCTTCCTCGCTGTTGCTCAGGA. vetor-sh-PPARβ#1, TTTCCAAAAAAGAAGGCCCGCAGCATCCTTCTCTTGAAAGGATGCTGCGGGCCTTCTGGG; vetor-sh-PPARβ#2(sc-36305-V) from Santa Cruz Biotechnology.
Antibodies, cellular protein extracts and Western blotting
Antibodies used were anti-PLCε, anti-PPARβ (santa cruz); anti-PPARα, anti-PPARγ, anti-Twist1, anti- phosphorylated (p)-Twist1 (Ser68), anti-PGC-1α, anti- CPT1B, anti-ERRα, anti-UCP-1, anti-ACADM, anti-MEKK3, anti-MEK4, anti-MEKK2, anti-ERK1/2, anti-p-MEK (S218/S222), anti-MEK, anti-p-JNK (Thr183/Tyr185), anti-JNK (Abcam); anti-p-ERK1/2 (Thr202/Tyr204), anti-p-P38, anti-P38, anti-B-Raf, N-Cadherin, E-Cadherin, Vimentin (Cell Signaling Technology); anti-Ras, anti-H3(wanlei); β-actin(Immunoway biotechnology); Anti-DDK (FLAG) (ORIGENE). The total protein, cell cytoplasmic protein and nuclear protein extracts and western blot assays were performed with the above antibodies by standard protocols we described before .
Total RNA was isolated from cells using TRIzol regents (invitrogen). Reverse transcription was performed with 1 μg RNA with a reverse transcriptase core kit (Eurogentec). For the real-time PCR analysis, aliquots of double-stranded cDNA were amplified using a SYBR Green PCR Kit (TaKaRa). The sequences of the primers were as follows: PLCε sense, 5′-GGAGAATCCTCGGTAG and anti-sense, 5′-GGTTGTCAGCGTATGTCC; Twist1 sense, 5′-TCCGCGTCCCACTAGCA and anti-sense, 5′-AGTTATCCAGCTCCAGAGTCTCTAGAC; PPARβ sense, 5′-GCCGCCCTACAACGAGATCA and anti-sense, 5′-CCACCAGCAGTCCGTCTTTGT; β-actin sense, 5′-GGGACCTGACTGACTACCTC and anti-sense, 5′-ACGAGACCACCTTCAACTCCAC. All gene expression experiments were repeated at least 3 times.
Cells grown on a sterile cover slip were washed with PBS, fixed in 4% paraformaldehyde, permeabilized with 1% Triton X-100 and subsequently subjected to stain with anti-Twist (Abcam) at 4 °C overnight.
Measurement of oxygen consumption rates
The oxygen consumption rates (OCR) of cells were detected by an oxygen electrode assay with 10 mM glucose. First, baseline cellular oxygen consumption (basal respiration) was measured. After oligomycin (1 μM), OCR depicts ATP-linked respiration (mitochondrial ATP production) by subtracting the oligomycin rate from baseline cellular OCR. fluoro-carbonyl cyanide phenylhydrazone (FCCP, 0.2 μM) collapses the inner membrane gradient, driving the electron transport chain (ETC) to function at maximal rate, and the maximal respiration is derived by subtracting non-mitochondrial respiration from the FCCP OCR. Finally, antimycin A/rotenone (Rot/AA, 0.5 μM) were added to inhibit ETC function, indicating the non-mitochondrial respiration. The mitochondrial respiratory reserve capacity was calculated by subtracting basal from maximal respiration. Data were calculated and graphs were plotted using Agilent Seahorse Wave Desktop software and report generator, MS Excel and GraphPad Prism. All experiments were performed in triplicates.
Transwell assay was were performed as previously described .
Dual luciferase reporter assay
Cells were transfected with 1 μg human E-cadherin and N-cadherin reporter and co-transfected with 1 μg sh-Twist1 vetor. Cells were collected 48 h after transfection and analyze during the dual-luciferase reporter assay system (Promega, Madison, WI, USA). The pRL-SV40 vector that shows constitutive expression of Renilla luciferase was co-transfected as an internal control to correct for differences in transfection. Three independent were performed in triplicate.
PC3 cells (4 × 106 cells) infected with sh-NC lentiviral, sh-PLCε lentiviral and vetor-sh-Twist1 were suspended in a 1:1 dilution of Matrigel and implanted subcutaneously into 6- to 8-week-old male nude mice (Chongqing Medical University Laboratory Animal Center). The growth of the primary tumour was monitored by external caliper measurement every 2–3 days. Six weeks later, the mice were sacrificed and the tumour tissues were surgically removed, measured, and fixed for histology studies.
Data are presented as mean ± SD. Statistical analysis was performed using a one-way analysis of variance and a Student’s t test (two-tailed) was used to compare two groups with the SPSS (Version 17.0) and GraphPad (Prism 5.0) software programs. Throughout this study: *p < 0.05; **p < 0.01; and *** p < 0.001.
Prognostic significance of Twist1 and PLCε expression in PCa
Analysis of the GEO (GSE21034), The Cancer Genome Atlas (TCGA), and the GTEx databases revealed high levels of Twist1 expression in PCa patients (Fig. 1a and Additional file 1: Figure S1a, P < 0.05). The GEO database also showed that Twist1 mRNA levels were significantly higher in metastatic PCa specimens than primary tumors (Fig. 1b, P < 0.05). Examining 496 human prostate samples from TCGA showed that the expression of PLCε was significantly associated with the Gleason grade (Fig. 1c, P < 0.05). Furthermore, high Twist1 and PLCε mRNA levels were associated with a shorter survival time in PCa patients (Fig. 1d, P = 0.021; 1e, P = 0.078). Data from GTEx and cBioPortal (http://www.cbioportal.org) revealed a positive correlation between the expression of Twist1 and PLCε (Fig. 1f and Additional file 1: Figure S1b, P < 0.001).
Fifty-five PCa tissue samples from patients treated with radical prostatectomy, and 48 BPH samples as controls, were tested for the presence of PLCε and Twist1. IHC analysis showed that approximately 70.9% (39/55) of the PCa samples were positive for PLCε staining, vs. 29.1% (16/55) of the BPH samples. Twist1 staining was positive in 85.5% (47/55) of the PCa samples vs. 14.5% (8/55) of the BPH samples (Fig. 1g). Semi-quantitative staining scores revealed significantly increased PLCε and Twist1 in PCa tissues, compared with that in BPH tissues (Fig. 1h, i, P < 0.001). This finding was confirmed by western blot analyses (Fig. 2a-c).
A positive association between PLCε and Twist1 IHC staining scores by Pearson’s correlation analysis was observed (Fig. 1j, n = 103, r = 0.8037, P < 0.001). Semi-quantitative analysis showed that PLCε and Twist1 protein expression in PCa tissue increased with increases in the degree of malignancy (Fig. 2d, e). Furthermore, a positive statistical association was found between the protein levels of PLCε and Twist1 (Fig. 2f, P < 0.001). Similar findings were observed with PLCε and Twist1 mRNA expression in PCa tissue (Fig. 2g-i). These results also demonstrated that the expression of PLCε and Twist1 was significantly associated with the Gleason grade in various clinicopathological characteristics of PCa patients (Table 1).
We next investigated whether high expression of Twist1 could be a prognostic marker of PCa due to its association with the metastasis of various cancers. Using a distinct group of clinical patient serum samples, we observed that serum Twist1 mRNA expression was higher in PCa than in BPH, and was also higher in metastatic than in non-metastatic PCa (Fig. 3a, P < 0.001). Furthermore, a high Twist1 mRNA serum level significantly correlated with the Gleason grade in various clinicopathological characteristics (Table 2).
Knockdown of PLCε suppresses Twist1 mRNA and protein expression in PCa cell lines
Four human PCa cell lines (PC3, DU145, LNCap, and Bica-R) were examined by western blotting. The results showed that these celllines exhibited varying levels of PLCε and Twist1 (Fig. 3b), and in response to bicalutamide-resistance, Twist1 expression increased in LNCap. PLCε expression was also upregulated in Bica-R cells, further indicating a possible correlation between PLCε and Twist1 levels. Therefore, subsequent experiments generally used the three cells lines PC3, DU145, and Bica-R.
To silence PLCε, PCa cells were infected with three lentiviral short hairpin (sh) RNAs specific for PLCε (LV-sh-PLCε#1, LV-sh-PLCε#2, and LV-sh-PLCε#3). LV-sh-PLCε#3 markedly reduced the mRNA and protein expression of both PLCε and Twist1 (Fig. 3c, e, f).We also found that vector-sh-Twist1#1 effectively downregulated Twist1 mRNA and protein expression without affecting PLCε expression (Fig. 3d, g, h).
PLCε expression negatively correlates with PGC-1α-mediated mitochondrial oxidative metabolism and uncoupling through Twist1
We next investigated whether downregulating PLCε in PCa cells had any effects on mitochondrial oxidative metabolism and uncoupling through inhibition of Twist1 expression. As shown in Fig. 4a and b, and Additional file 1: Figure S2a and b, stable knockdown of PLCε increased the expression of PGC-1α, CPT1B, ERRα, UCP-1, and ACADM, and decreased the expression of Twist1.
Based on these findings, we examined the mitochondrial activity of PCa cell lines by measuring the OCR. The results showed that knockdown of PLCε increased the OCR and ATP production in PCa cells (Fig. 4c, d and Additional file 1: Figure S2c, d). We hypothesized that PLCε regulated mitochondrial oxidative metabolism and uncoupling specifically through Twist1. Western blotting and OCR analyses showed that overexpressing Twist1 in sh-PLCε-transfected PCa cells abolished the sh-PLCε-mediated increases in metabolism, uncoupling, OCR, and ATP production (Fig. 4e-h and Additional file 1: Figure S2e-h). Interestingly, the combination of sh-PLCε and sh-Twist1 promoted metabolism and uncoupling significantly more than either shRNA alone (Fig. 4i-l and Additional file 1: Figure S2i-l). Together, these data support our hypothesis that the PLCε/Twist1 axis suppresses PCa cell mitochondrial oxidative metabolism and uncoupling.
PLCε depletion decreases phosphorylation of Twist1 on serine 68 and reduces stability of Twist1 protein via MAPKs
At the gene level, the expression of Twist1 is modulated by several upstream regulators through multiple pathways, depending on the cancer type and tissue context . For example, signal transducer and activator of transcription 3 binds directly to the human Twist1 promoter and activates its transcriptional activity in human breast cancer, hepatocellular carcinoma, and gastric cancer cell lines [36,37,38]. Moreover, NF-κB, Ras, and transforming growth factor-β signaling pathways enhance Twist1 expression [39,40,41]. At the post-translational level, the MAPK pathway decreases E3-mediated ubiquitination and stabilizes Twist1 without altering mRNA expression through a considerable increase in Twist1 serine68 phosphorylation in breast cancer cells . PLCε can also regulate Ras expression .
To explore the mechanisms through which PLCε modulate Twist1 expression, we analyzed the MAPK signaling pathway, the main pathway upstream of Twist. Western blotting showed decreased levels of Twist 1 and phosphorylated (p)Twist1 (Fig. 5a, b), as well as decreased p-MEK, p-ERK, p-P38, and p-c-Jun N-terminal kinase (JNK), and some key molecules in the MAPK pathway (Fig. 5c, d and Additional file 1: Figure S3a, b) in PCa cells treated with PLCε-shRNA. We subsequently utilized 1 nM trametinib (a specific inhibitor of MEK), 5 nM JNK-In-8 (a specific inhibitor of JNK), or 0.5 μM SB203580 (a specific inhibitor of P38 MAPK) to evaluate whether these signaling pathways were involved in PLCε-mediated Twist1 expression. A concentration achieving 50% inhibition efficiency was selected for each inhibitors (Additional file 1: Figure S3c-e). The results showed that all three inhibitors effectively antagonized their corresponding targets and decreased Twist1 protein levels in the absence or presence of PLCε but had little effect on Twist1 mRNA levels (Fig. 5e-j, Additional file 1: Figures S3f-j and S4a-c).
Finally, we assessed the stability of the Twist1 protein in PC3 cells after blocking protein synthesis with 100 μM cycloheximide. The half-life of Twist1 in DDK-Twist1-transfected PC3 cells was 5.78 h (Fig. 6a). This suggests that Twist1 protein has a short half-life and its stability is affected by ubiquitin modification. Interestingly, sh-PLCε-mediated inhibition of Twist1 protein levels could be reversed with 25 μM MG132 (a proteasome inhibitor) (Fig. 6b). Collectively, these finding indicated that the regulation of Twist1 by PLCε is probably proteasome-dependent, and that PLCε regulates Twist1 proteasomal degradation via MAPKs.
PLCε regulates Twist1 through nuclear translocation of PPARβ
In PCa cells, knockdown of PLCε decreased the PPARβ protein level but had little effect on PPARα and PPARγ (Fig. 6c, d and Additional file 1: Figure S4d, e). PCa cells were transfected with two shRNAs specific for PPARβ (vector-sh-PPARβ#1 and vector-sh-PPARβ#2). Vector-sh-PPARβ#2 effectively downregulated PPARβ mRNA and protein expression (Fig. 6e, f and Additional file 1: Figure S4f). Interestingly, the combination of sh-PLCε and sh-PPARβ inhabited Twist1 mRNA and protein expression significantly more than either shRNA alone (Fig. 6g, h and Additional file 1: Figure S5c-e). PPARβ is a transcription factor that mediates lipid and fatty acid metabolism. Immunofluorescence staining showed that PLCε depletion markedly decreased the PPARβ signal in the nucleus (Fig. 6j). To confirm the distribution of PPARβ within the nucleus, we performed a nuclear-cytoplasmic separation. Western blotting of each fraction showed that the nuclear distribution of PPARβ was notably decreased in PCa cells (Fig. 6k). We surmised that PLCε depletion inhibits PPARβ by crippling its nuclear translocation.
The Twist1 luciferase reporter assay showed that PLCε and PPARβ knockdown suppressed Twist1 activities in PC3 cells. This indicated that PPARβ is a direct transcriptional target of Twist1 (Fig. 6i). Treatment of sh-PLCε-transfected PCa cells with1 nM GW501516, a highly selective agonist of PPARβ, prevented the PLCε depletion-mediated decreases in Twist1 mRNA and protein, while GSK3787 (1 μM), an inhibitor of PPARβ, enhanced PLCε depletion-mediated Twist1 repression (Fig. 7a, b and Additional file 1: Figure S5a-c).
To further support the hypothesis that PLCε promotes Twist1 transcription though PPARβ, assays were performed in PC3 cells with E-cadherin or N-cadherin promoter-driven luciferase reporters after depletion of PLCε or Twist1. The results demonstrated that PLCε and Twist1 knockdown enhanced E-cadherin and decreased N-cadherin reporter activities. Furthermore, GW501516 prevented and GSK3787 enhanced these effects (Fig. 7e). Collectively, these results indicate that PLCε induces Twist1 mRNA, protein, and transcriptional activity through PPARβ.
PLCε depletion inhibits EMT via PPARβ-mediated Twist1 downregulation
To investigate whether PLCε regulates Twist1-mediated EMT and cell migration, western blotting analyses were performed. The results revealed that PLCε depletion inhibited the expression of Twist1, PPARβ, N-cadherin, and vimentin, and increased the expression of E-cadherin. GW501516 antagonized and GSK3787 promoted the effects of PLCε depletion in PC3 and DU145 cells (Fig. 7c, d and Additional file 1: Figure S5f, g). Finally, Transwell migration assays revealed that PLCε depletion inhibited the migration of PCa cells, and GW501516 antagonized and GSK3787 promoted this effect (Fig. 7f, g).
PLCε depletion enhances the inhibitory effect of Twist1 knockdown on PCa tumorigenesis in vivo
An in vivo nude mouse tumor xenograft assay was conducted to verify the role of PLCε in tumor formation. Consistent with the in vitro findings, the PLCε- and Twist1-depleted groups both displayed less prominent tumorigenesis than the sh-NC group; the mean tumor volumes and weights were both lower in the PLCε- and Twist1-depleted groups than in the sh-NC group. Combined knockdown inhibited the growth and weight of tumors more than either knockdown alone (Fig. 8a-c). In addition, IHC analyses of the mouse tumors revealed that PLCε and Twist1 were expressed at low levels and PGC-1α was expressed at a high level in the sh-PLCε group compared to that in the sh-NC group. Downregulation of Twist1 alone decreased Twist1 and increased PGC-1α levels, while the combined treatment group expressed less PLCε and Twist1, and more PGC-1α than other three groups (Fig. 8d). These results support the tumor suppressor function of PLCε depletion in PCa.
The results presented here reveal that PLCε exerts tumor promotional activity in PCa cell lines through effects on mitochondrial oxidative metabolism. We demonstrated that PLCε depletion could regulate Twist1 mRNA levels and protein stability through two possible mechanisms, and thus promote PGC-1α-mediated mitochondrial oxidative metabolism. Moreover, the oxidative metabolic program elicited by PGC-1α can prevent tumor growth and migration.
PLCε and Twist1are associated with multiple types of cancers, including PCa [43,44,45]. We used patient prostate and serum samples to provide a comprehensive analysis of PLCε and Twist1 expression in normal (i.e., BPH) and cancerous (i.e., PCa) tissues. The results indicated that high expression of PLCε in PCa tissues correlated with low survival and a high Gleason score. Significantly elevated Twist1 expression in prostate tissues and serum from PCa patients also correlated significantly with low survival and a high Gleason score. Moreover, a statistical association was found between PLCε and Twist1 expression in PCa. Thus, measurements of PLCε and Twist1 might have diagnostic and/or prognostic value in PCa.
PLCε is required for the calcium release mechanism  that is responsible for the activity of the mTOR and MAPK pathways. Importantly, it is widely believed that mitochondrial oxidative metabolism is controlled by both the mTOR and MAPK pathways [47, 48]. Thus, we hypothesized that PLCε regulated mitochondrial oxidative metabolism in PCa. In our studies, we knocked down PLCε in PCa cell lines using lentiviral sh-PLCε. PLCε depletion decreased Twist1 mRNA and protein levels. This implied that PLCε may regulate Twist1 mRNA levels. PLCε depletion also increased the expression of PGC-1α, CPT1B, ERRα, UCP-1, and ACADM, and promoted the OCR. By transfecting cells with a Twist1 (Myc-DDK-tagged) vector and vector-sh-Twist1, we clearly found that PLCε regulated PCa mitochondrial oxidative metabolism by targeting Twist1.
MAPKs can promote breast tumor cell EMT and metastasis via phosphorylation and stabilization of Twist1, and PLCε can regulate key molecules involved in MAPK expression . We observed that knocking down PLCε decreased p-Twist1expression significantly, and an inhibitor of the MAPK pathway suppressed both Twist1 and p-Twist1 levels. Interestingly, sh-PLCε-mediated inhibition of Twist1 protein levels could be reversed with MG132. These results suggest that PLCε depletion suppresses the phosphorylation of Twist1 and accelerates its ubiquitination, thus decreasing Twist1 protein levels via all three MAPK subfamily members.
Interestingly, PLCε depletion decreased Twist1 mRNA and protein levels, which implies that PLCε regulates Twist1 mRNA levels in another way. PPARβ, a member of the PPAR nuclear hormone receptor super family, is involved in brain lipid metabolism, proliferation of anterior fat cells, fat formation, embryonic inhibition, macrophage cholesterol homeostasis, and tumor formation after being activated by ligands . Some studies suggest that PPARβ is highly expressed in PCa cells and is associated with cell cycle regulation and proliferation . PLC influences the expression of PPARβ through Ca2+/cytosolic phospholipase A2, thus affecting glucose metabolism [51, 52]. Following PLCε knockdown, we observed a decreased level of PPARβ alone. Using GW501516 and GSK3787, an agonist and inhibitor of PPARβ, respectively, we showed that PLCε might regulate the transcription of Twist1, as well as Twist1-mediated EMT, cell migration, and transcriptional regulation; these are well-known features of Twist1 . Importantly, immunofluorescence and western blot analyses indicated that PLCε depletion probably blocked PPARβ nuclear translocation.
The results of our study have elucidated a probable role for PLCε in antagonizing mitochondrial oxidative metabolism in PCa. We also identified that these changes are mediated by Twist1. Moreover, PLCε may regulate Twist1 expression through several pathways, including MAPK-mediated ubiquitination and PPARβ-mediated transcription.
Availability of data and materials
Literature collection was performed using PubMed. Statistical analysis were executed by using SPSS 17.0 software (IBM, Chicago, IL, USA) and GraphPad (Prism 5.0). Raw and processed data are stored in corresponding author and are available upon request. The RNAseq datasets and clinical information of Twist1 and PLCε in PCa were retrieved from TCGA database (http://www.cancergenome.nih.gov), cBioPortal for Cancer Genomics (http://www.cbioportal.org/), GEO Datasets (https://www.ncbi.nlm.nih.gov/gds) and GTEx portal (https://www.gtexportal.org).
Benign prostatic hypertrophy
Electron transport chain
Luoro-carbonyl cyanide phenylhydrazone
Jun N-terminal kinase
Mitogen-activated protein kinases
Oxygen consumption rates
Peroxisome proliferator-activated receptor gamma co-activator 1 alpha
Peroxisome proliferator activated receptor
Beltran H, Tomlins S, Aparicio A, Arora V, Rickman D, Ayala G, Huang J, True L, Gleave ME, Soule H, et al. Aggressive variants of castration-resistant prostate cancer. Clin Cancer Res. 2014;20:2846–50.
Bostwick DG. Staging prostate cancer--1997: current methods and limitations. Eur Urol. 1997;32(Suppl 3):2–14.
Herst PM, Grasso C, Berridge MV. Metabolic reprogramming of mitochondrial respiration in metastatic cancer. Cancer Metastasis Rev. 2018;37:643-53.
Wallace M, Metallo CM. PGC1alpha drives a metabolic block on prostate cancer progression. Nat Cell Biol. 2016;18:589–90.
Liberti MV, Locasale JW. The Warburg effect: how does it benefit Cancer cells? Trends Biochem Sci. 2016;41:211–8.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324:1029–33.
Puigserver P, Wu Z, Park CW, Graves R, Wright M, Spiegelman BM. A cold-inducible coactivator of nuclear receptors linked to adaptive thermogenesis. Cell. 1998;92:829–39.
Michael LF, Wu Z, Cheatham RB, Puigserver P, Adelmant G, Lehman JJ, Kelly DP, Spiegelman BM. Restoration of insulin-sensitive glucose transporter (GLUT4) gene expression in muscle cells by the transcriptional coactivator PGC-1. Proc Natl Acad Sci U S A. 2001;98:3820–5.
Lin J, Handschin C, Spiegelman BM. Metabolic control through the PGC-1 family of transcription coactivators. Cell Metab. 2005;1:361–70.
Audet-Walsh E, Papadopoli DJ, Gravel SP, Yee T, Bridon G, Caron M, Bourque G, Giguere V, St-Pierre J. The PGC-1alpha/ERRalpha Axis represses one-carbon metabolism and promotes sensitivity to anti-folate therapy in breast Cancer. Cell Rep. 2016;14:920–31.
Torrano V, Valcarcel-Jimenez L, Cortazar AR, Liu X, Urosevic J, Castillo-Martin M, Fernandez-Ruiz S, Morciano G, Caro-Maldonado A, Guiu M, et al. The metabolic co-regulator PGC1alpha suppresses prostate cancer metastasis. Nat Cell Biol. 2016;18:645–56.
Schopf B, Schafer G, Weber A, Talasz H, Eder IE, Klocker H, Gnaiger E. Oxidative phosphorylation and mitochondrial function differ between human prostate tissue and cultured cells. FEBS J. 2016;283:2181–96.
Bunney TD, Katan M. Phospholipase C epsilon: linking second messengers and small GTPases. Trends Cell Biol. 2006;16:640–8.
Bunney TD, Katan M. PLC regulation: emerging pictures for molecular mechanisms. Trends Biochem Sci. 2011;36:88–96.
Singer WD, Brown HA, Sternweis PC. Regulation of eukaryotic phosphatidylinositol-specific phospholipase C and phospholipase D. Annu Rev Biochem. 1997;66:475–509.
Jin TG, Satoh T, Liao Y, Song C, Gao X, Kariya K, Hu CD, Kataoka T. Role of the CDC25 homology domain of phospholipase Cepsilon in amplification of Rap1-dependent signaling. J Biol Chem. 2001;276:30301–7.
Sorli SC, Bunney TD, Sugden PH, Paterson HF, Katan M. Signaling properties and expression in normal and tumor tissues of two phospholipase C epsilon splice variants. Oncogene. 2005;24:90–100.
Wang X, Zhou C, Qiu G, Yang Y, Yan D, Xing T, Fan J, Tang H, Peng Z. Phospholipase C epsilon plays a suppressive role in incidence of colorectal cancer. Med Oncol. 2012;29:1051–8.
Danielsen SA, Cekaite L, Agesen TH, Sveen A, Nesbakken A, Thiis-Evensen E, Skotheim RI, Lind GE, Lothe RA. Phospholipase C isozymes are deregulated in colorectal cancer--insights gained from gene set enrichment analysis of the transcriptome. PLoS One. 2011;6:e24419.
Bunney TD, Baxendale RW, Katan M. Regulatory links between PLC enzymes and Ras superfamily GTPases: signalling via PLCepsilon. Adv Enzym Regul. 2009;49:54–8.
Zhou RM, Li Y, Wang N, Liu BC, Chen ZF, Zuo LF. PLC-epsilon1 gene polymorphisms significantly enhance the risk of esophageal squamous cell carcinoma in individuals with a family history of upper gastrointestinal cancers. Arch Med Res. 2012;43:578–84.
Wang Y, Wu X, Ou L, Yang X, Wang X, Tang M, Chen E, Luo C. PLCepsilon knockdown inhibits prostate cancer cell proliferation via suppression of notch signalling and nuclear translocation of the androgen receptor. Cancer Lett. 2015;362:61–9.
Wang LD, Zhou FY, Li XM, Sun LD, Song X, Jin Y, Li JM, Kong GQ, Qi H, Cui J, et al. Genome-wide association study of esophageal squamous cell carcinoma in Chinese subjects identifies susceptibility loci at PLCE1 and C20orf54. Nat Genet. 2010;42:759–63.
Thisse B, el Messal M, Perrin-Schmitt F. The twist gene: isolation of a Drosophila zygotic gene necessary for the establishment of dorsoventral pattern. Nucleic Acids Res. 1987;15:3439–53.
Yang J, Mani SA, Donaher JL, Ramaswamy S, Itzykson RA, Come C, Savagner P, Gitelman I, Richardson A, Weinberg RA. Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell. 2004;117:927–39.
Martin A, Cano A. Tumorigenesis: Twist1 links EMT to self-renewal. Nat Cell Biol. 2010;12:924–5.
Mani SA, Guo W, Liao MJ, Eaton EN, Ayyanan A, Zhou AY, Brooks M, Reinhard F, Zhang CC, Shipitsin M, et al. The epithelial-mesenchymal transition generates cells with properties of stem cells. Cell. 2008;133:704–15.
Qin Q, Xu Y, He T, Qin C, Xu J. Normal and disease-related biological functions of Twist1 and underlying molecular mechanisms. Cell Res. 2012;22:90–106.
Lee KE, Bar-Sagi D. Oncogenic KRas suppresses inflammation-associated senescence of pancreatic ductal cells. Cancer Cell. 2010;18:448–58.
Parajuli P, Kumar S, Loumaye A, Singh P, Eragamreddy S, Nguyen TL, Ozkan S, Razzaque MS, Prunier C, Thissen JP, et al. Twist1 Activation in Muscle Progenitor Cells Causes Muscle Loss Akin to Cancer Cachexia. Dev Cell. 2018;45:712–725.e6.
Pan D, Fujimoto M, Lopes A, Wang YX. Twist-1 is a PPARdelta-inducible, negative-feedback regulator of PGC-1alpha in brown fat metabolism. Cell. 2009;137:73–86.
Quan Z, He Y, Luo C, Xia Y, Zhao Y, Liu N, Wu X. Interleukin 6 induces cell proliferation of clear cell renal cell carcinoma by suppressing hepaCAM via the STAT3-dependent up-regulation of DNMT1 or DNMT3b. Cell Signal. 2017;32:48–58.
Du Z, Li L, Sun W, Wang X, Zhang Y, Chen Z, Yuan M, Quan Z, Liu N, Hao Y, et al. HepaCAM inhibits the malignant behavior of castration-resistant prostate cancer cells by downregulating notch signaling and PF-3084014 (a gamma-secretase inhibitor) partly reverses the resistance of refractory prostate cancer to docetaxel and enzalutamide in vitro. Int J Oncol. 2018;53:99–112.
Alexander NR, Tran NL, Rekapally H, Summers CE, Glackin C, Heimark RL. N-cadherin gene expression in prostate carcinoma is modulated by integrin-dependent nuclear translocation of Twist1. Cancer Res. 2006;66:3365–9.
Zhao Z, Rahman MA, Chen ZG, Shin DM. Multiple biological functions of Twist1 in various cancers. Oncotarget. 2017;8:20380–93.
Cheng GZ, Zhang WZ, Sun M, Wang Q, Coppola D, Mansour M, Xu LM, Costanzo C, Cheng JQ, Wang LH. Twist is transcriptionally induced by activation of STAT3 and mediates STAT3 oncogenic function. J Biol Chem. 2008;283:14665–73.
Zhang C, Guo F, Xu G, Ma J, Shao F. STAT3 cooperates with twist to mediate epithelial-mesenchymal transition in human hepatocellular carcinoma cells. Oncol Rep. 2015;33:1872–82.
Hsu KW, Hsieh RH, Huang KH, Fen-Yau Li A, Chi CW, Wang TY, Tseng MJ, Wu KJ, Yeh TS. Activation of the Notch1/STAT3/twist signaling axis promotes gastric cancer progression. Carcinogenesis. 2012;33:1459–67.
Jiang J, Kosman D, Ip YT, Levine M. The dorsal morphogen gradient regulates the mesoderm determinant twist in early Drosophila embryos. Genes Dev. 1991;5:1881–91.
Satoh K, Hamada S, Kimura K, Kanno A, Hirota M, Umino J, Fujibuchi W, Masamune A, Tanaka N, Miura K, et al. Up-regulation of MSX2 enhances the malignant phenotype and is associated with twist 1 expression in human pancreatic cancer cells. Am J Pathol. 2008;172:926–39.
Thuault S, Valcourt U, Petersen M, Manfioletti G, Heldin CH, Moustakas A. Transforming growth factor-beta employs HMGA2 to elicit epithelial-mesenchymal transition. J Cell Biol. 2006;174:175–83.
Hong J, Zhou J, Fu J, He T, Qin J, Wang L, Liao L, Xu J. Phosphorylation of serine 68 of Twist1 by MAPKs stabilizes Twist1 protein and promotes breast cancer cell invasiveness. Cancer Res. 2011;71:3980–90.
Han KY, Chen PN, Hong MC, Hseu YC, Chen KM, Hsu LS, Chen WJ. Naringenin attenuated prostate Cancer invasion via reversal of epithelial-to-mesenchymal transition and inhibited uPA activity. Anticancer Res. 2018;38:6753–8.
Yang Y, Wu XH. Study on the influence of metformin on castration-resistant prostate cancer PC-3 cell line biological behavior by its inhibition on PLCepsilon gene-mediated Notch1/Hes and androgen receptor signaling pathway. Eur Rev Med Pharmacol Sci. 2017;21:1918–23.
Wang X, Fan Y, Du Z, Fan J, Hao Y, Wang J, Wu X, Luo C. Knockdown of phospholipase Cepsilon (PLCepsilon) inhibits cell proliferation via phosphatase and Tensin homolog deleted on chromosome 10 (PTEN)/AKT signaling pathway in human prostate Cancer. Med Sci Monit. 2018;24:254–63.
Kondratskyi A, Yassine M, Kondratska K, Skryma R, Slomianny C, Prevarskaya N. Calcium-permeable ion channels in control of autophagy and cancer. Front Physiol. 2013;4:272.
Corazao-Rozas P, Guerreschi P, Andre F, Gabert PE, Lancel S, Dekiouk S, Fontaine D, Tardivel M, Savina A, Quesnel B, et al. Mitochondrial oxidative phosphorylation controls cancer cell’s life and death decisions upon exposure to MAPK inhibitors. Oncotarget. 2016;7:39473–85.
Cunningham JT, Rodgers JT, Arlow DH, Vazquez F, Mootha VK, Puigserver P. mTOR controls mitochondrial oxidative function through a YY1-PGC-1alpha transcriptional complex. Nature. 2007;450:736–40.
Duval C, Chinetti G, Trottein F, Fruchart JC, Staels B. The role of PPARs in atherosclerosis. Trends Mol Med. 2002;8:422–30.
Jarvis MC, Gray TJ, Palmer CN. Both PPARgamma and PPARdelta influence sulindac sulfide-mediated p21WAF1/CIP1 upregulation in a human prostate epithelial cell line. Oncogene. 2005;24:8211–5.
Lo SH, Lee KS, Chen LJ, Cheng JT, Chen CH. Increase of PPARdelta by dopamine mediated via DA-1 receptor-linked phospholipase C pathway in neonatal rat cardiomyocytes. Auton Neurosci. 2013;177:211–6.
Suh HN, Huong HT, Song CH, Lee JH, Han HJ. Linoleic acid stimulates gluconeogenesis via Ca2+/PLC, cPLA2, and PPAR pathways through GPR40 in primary cultured chicken hepatocytes. Am J Physiol Cell Physiol. 2008;295:C1518–27.
7TFP CDH1 reporter was a gift from Bob Weinberg (Addgene plasmid # 91704; http://n2t.net/addgene:91704; RRID:Addgene_91704). The authors thank Dr. Bob Weinberg.
This study is one of the projects of the National Natural Science Foundation of China (NSFC) (Grant No.81802543).
In our study, All samples were collected with informed consent and the experiments were approved by the institutional ethics committee of Chongqing medical university. All animal experiments were carried out with approval from the institutional ethics committee of the Chongqing Medical University, strictly following the guidelines set for the usage of animals by this committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. High expression levels of PLCε and Twist1 in PCa tissues. (a) RNA-seq mRNA expression data from the TCGA and GTEx database was used to compare Twist1 expression between in PCa tumors (T) (n = 492) and their non-tumor counterparts (N) (n = 152). *p < 0.05. (b) Correlation between the expression of Twist1 and PLCε in PCa patients from cBioPortal database. Figure S2. PLCε expression negatively correlates with PGC-1α-mediated mitochondrial oxidative metabolism and uncoupling through Twist1 in DU145 and Bica-R cells. Figure S3. PLCε depletion decreases phosphorylation of Twist1 on serine 68 and reduces stability of Twist1 protein via MAPKs in DU145 and Bica-R cells. Figure S4 PLCε depletion decreases Twist1 protein via MAPKs and alleviates the PPARβ expression. Figure S5 PLCε depletion decreases Twist1 by PPARβ. (DOCX 3776 kb)