CYLD were significantly downregulated in ovarian cancer samples and lower CYLD expression indicated poor prognosis.
To identify the role of CYLD in cancer development, the expression levels of CYLD in different types of cancer were detected to be generally lower in pan-cancer used UCSC Xena database, with especially in TCGA database of ovarian cancer, the expression levels of CYLD were significantly decreased (Supplemental Fig.1). Then the expression levels of CYLD were determined in tissue samples of 10 normal ovarian (Normal) and 33 epithelial ovarian cancer (Tumor) using immunohistochemistry, and the tumor samples had much lower CYLD expression levels and lower IHC sores compared with Normal samples (Fig.1A). In addition, the protein and mRNA expression levels of CYLD were significantly downregulated in epithelial ovarian cancer (Tumor) samples compared with the normal ovarian (Normal) tissues (Fig.1B). We further studied the relationship between CYLD expression and the clinicopathological characteristics in ovarian cancer tissues, and conveyed that there were no significant correlation in the clinical characteristics, such as age, ascites volume and metastasis site. Interestingly, the CYLD expression level decreased with the increase of the clinical stage, as well as higher CYLD expression indicated more sensitive to cisplatin treatment (Supplemental Tab.4). And the patients with lower CYLD expression levels had unfavorable outcomes compared to those with higher CYLD expression group (Fig.1C). Furthermore, CYLD expression was explored with the GEPIA database, and the results confirmed that CYLD was a tumor suppressor in ovarian cancer and the expression level decreased in higher clinical stage(Fig.1D). The prognostic roles of CYLD in ovarian cancer were further evaluated using Kaplan-Meier plotter database, and higher CYLD expression had better survival rate(Fig.1E). Thus, we firstly suggested that CYLD was a tumor suppressor in ovarian cancer, and the patients with lower CYLD expression had poor prognosis in the patient cohort.
CYLD could be an essential DDP sensitizer biomarker.
Since resistance to DDP treatment is a major clinical challenge in ovarian cancer therapy, therefore, we considered that CYLD may be associated with cisplatin resistance in ovarian cancer. Then CYLD expression levels were detected to be reduced in 24 tumor tissues lesions from cisplatin sensitivity group, compared with resistance group (Fig.2A). And the mRNA expression levels of CYLD were significantly downregulated in resistance samples compared to sensitive samples (Fig.2B).
We further selected 2 cases (Case 1 from cisplatin resistant group and Case 2 from cisplatin sensitive group ) to illustrate the crucial role of CYLD in cisplatin resistant ovarian cancer. Here, using computed tomography (CT) analysis with tumor tissues of 2 clinical patients, Case 1 showed pelvic shadow in CT image at 4 months after treatment (Fig.2C, right panel), Case 2 showed no recurrence as suggested with the CT images at six months after treatment (Fig.2D, right panel). Then regular review was checked by serum biomarker Ca125 levels and imaging examination every 2 months. Two months later, Ca125 levels were found to be raised in Case 1, but not changed in Case 2, which also suggested platinum resistance recurrence in Case 1 (Fig.2E). The mRNA expression levels of CYLD in Case 1 and Case 2 also confirmed that CYLD expression levels were significantly lower in platinum resistant ovarian cancer tissues (Fig. 2F). In addition, immuno-histochemistry staining for CYLD of these tumor tissues were conducted using the surgical specimens, the CYLD expression levels showed negative in Case 1 patient with platinum resistance recurrence, but positive in platinum sensitive Case 2 patient (Fig. 2G). Thus, we firstly identified that CYLD expression was attenuated in ovarian cancer patient tissues when the cancer became DDP resistance, and CYLD could be an essential DDP sensitizer biomarker.
CYLD knockdown rendered cells resistant to DDP, while force expression of CYLD promoted drug sensitivity.
To identify the expression of CYLD in DDP resistance ovarian cancer, we used the parental and drug resistant cell lines (OVCAR3/OVCAR3-DDP, A2780/A2780-DDP) to conduct analysis, and found that the protein expression levels of CYLD were significantly downregulated in DDP resistant cell lines compared to parental cells (Fig. 3A). Similarly, much lower levels of CYLD were identified in other types of OC cell lines (Caov3, SKOV3, A2780 and OVCAR3) compared to human immortalized ovarian epithelial cell line (IOSE386) (Fig.3B). These results showed that CYLD expression was suppressed in progression of OC cisplatin resistance, but whether it is an initiating factor of drug resistance still needed to be investigated. To test whether CYLD suppression is sufficient in affecting chemoresistance, we showed that CYLD knockdown in OVCAR3 and A2780 cells were sufficient to promote cells more resistance to DDP treatment compared to control cells by CCK-8 assay (Fig.3C-3D). Additionally, the CYLD-overexpressed OVCAR3-DDP and A2780-DDP cells were established. Overexpression of CYLD was sufficient to increase cell sensitivity to DDP treatment compared with control group, making the DDP resistant cell lines sensitive to DDP (Fig.3E-3F). These findings suggested that CYLD expression is essential in regulating DDP resistance.
CYLD knockdown inhibited cell apoptosis activities and promoted drug efflux.
Since previous studied reported that loss of CYLD expression may cause cisplatin resistance through NF-κB hyperactivation [22, 23], we wondered whether CYLD regulated DDP resistance through regulating cellular apoptosis in ovarian cancer. The apoptosis rates were analyzed in OVCAR3 and A2780 cells treated with 5μM DDP or PBS, and the percentages of apoptotic cells in CYLD knockdown group significantly decreased compared to that in the control group with DDP treatment (Fig.4A, Supplemental Fig.2A-2B). Then we found that the protein expression levels of pro-apoptosis factor Bax were downregulated in CYLD knockdown group, while the anti-apoptotic proteins Bcl-XL and Bcl-2 were upregulated (Fig.4B). In addition, the percentage of apoptosis in A2780-DDP with overexpressed CYLD group was remarkably higher than that of A2780-DDP, and the result was similar to that of A2780 (Fig.4C, Supplemental Fig.2C). Further studies suggested that compared to OVCAR3-DDP/A2780-DDP, Bax levels were upregulated in CYLD-overexpressed DDP resistant cells and OVCAR3/A2780 cells, while Bcl-XL levels were downregulated, respectively (Fig.4D). These results indicated that CYLD mediated DDP resistance through apoptosis via Bax/Bcl-XL pathway.
Base on the effect of CYLD in sensitizing cells to apoptosis might be cell type specific[20, 37]. Therefore, we considered whether there were other molecular regulatory pathways that caused drug resistance in ovarian cancer. Thus, we adopted the intracellular levels of Rhodamine 123, a fluorescent drug that can be excreted through drug transporters on the surface of cell membrane, and measured transporter activities between CYLD knockdown OVCAR3 cells and control. After incubated with different concentrations of Rhodamine 123, the fluorescence intensity of CYLD knockdown group was higher than control group, and the difference was more remarkable at the concentration of 2 μM (Fig.4E). Then the fluorescence intensities of cells at different time points after incubated Rhodamine 123 were also measured (Supplemental Fig. 3A), which showed that the fluorescence intensities were significantly different in 12h group, while there was no difference after 16h treatment (Supplemental Fig. 3B). These results suggested that CYLD knockdown promoted drug to pump intracellularly. To further examine the potential targets of CYLD in ABC transporter superfamily, we found that ABCB1 were significantly induced in CYLD-knockdown cells (Fig.4F), while ABCC1, ABCG2 and ABCG9 were not much difference(Supplemental Fig.3C-3D). In addition, ABCB1 in OVCAR3-DDP/A2780-DDP cells were showed to be greatly higher than those in other cells (OVCAR3-DDP/A2780-DDP with CYLD overexpression and OVCAR3/A2780 cells) (Fig.4G). Since CYLD is a deubiquitinase, we examined whether CYLD regulates ABCB1 expression through ubiquitination, but there is no binding between CYLD and ABCB1 by IP experiment, which suggested that CYLD may indirectly regulate ABCB1(Supplemental Fig.3E). To further gain insight into the molecular mechanism, we found that key intermediate signaling molecule p-65(ser536) in CYLD knockdown cells were induced compared with control cells (Fig.4H). Furthermore, two binding sites on the promoter region of ABCB1 were predicted in JASPAR according to the binding motif of p65, and the site 1(~1839-1848) was the direct binding site of p65 in ABCB1 promoter (Fig.4I). Hence, ABCB1 was a pivotal downstream mediator of CYLD-regulating DDP resistance.
HER3 negatively regulated CYLD expression via phosphorylation of STAT3
To further investigate the molecular mechanisms in regulation of CYLD, the mRNA expression levels of the ERBB family (EGFR, HER2, HER3, HER4) were detected in patient samples. There was a significantly negative relationship between the expression levels of HER3 and CYLD (R=-0.4333, P=0.0389, Fig.5A), while other ERBB family members had no significant correlations (Supplemental Fig.4A). Therefore, HER3 may be potential upstream regulator of CYLD in ovarian cancer. To further confirm the role of HER3 in CYLD regulation, we found that the protein expression levels of CYLD were decreased in HER3 overexpressed OVCAR3/A2780 cells, compared to control cells (Fig.5B). In addition, the expression levels of CYLD were higher in HER3 knockdown OVCAR3/A2780 cells compared to control cells (Fig.5C). These findings further suggested that HER3 may negatively regulate CYLD expression. To examine the mechanism that how HER3 attenuated CYLD’s expression, HER3-overexpressed cells showed much higher expression levels of key intermediate signaling molecule p-STAT3(T705) compared with control cells (Fig.5D), while other candidate molecules, p-c-Jun, p-p38 and p-Erk1/2 showed no significantly changed (Supplemental Fig.4B). Moreover, the expression levels of p-STAT3 in HER3-silenced OVCAR3/A2780 cells were significantly decreased compared to the control (Fig.5D). Furthermore, the immunoprecipitation assay showed that HER3 protein directly bound and interacted with CYLD protein (Fig.5E). Thus, these results showed that HER3 inhibited CYLD expression via phosphorylation of STAT3.
HER3 functioned as crucial upstream regulator in DDP resistance and HER3 inhibitor facilitated DDP sensitivity.
To further confirm the role of HER3 in ovarian cancer drug resistance, the protein levels of HER3 in OVCAR-DDP and A2780-DDP cells were showed higher compared to parental cells (Fig. 6A). The HER3-silenced cell lines showed more sensitive to drug treatment, while HER3-overexpressed cell lines found to be more resistant to DDP treatment (Fig.6B-6C). These results illustrated that HER3 played an important role in drug resistance of ovarian cancer cells.
In order to test whether HER3 inhibitor could improve the DDP sensitivity, the drug combination effect of HER3 inhibitor (TX1-85-1) and DDP in A2780-DDP cell line were analyzed. To identify combination effect, Chou-Talalay method were used to detect the synergistic effect of DDP (0, 0.1, 0.3, 0.5, 0.7, 1, 3, 5, 7, 10μM) in combination with TX1-85-1(0, 0.3, 0.7, 1, 3, 7μM). The Combination Index (CI) data of these two drugs was showed in Supplemental Fig. 4C. The Fraction Affected-combination index (Fa) analyzed the effects of the combination with different doses of TX1-85-1 and DDP (Fig.6D). The synergistic effect of these two drugs was described in a Scatter plot, and the fractions affected-combination index (Fa-CI) plot were almost below 1, that indicated that these two drugs had synergistic effect (Fig.6E). The synergistic effect suggested that inhibition of HER3 improved drug sensitivity in DDP resistance OC cells. Then according to the CI results, we selected two drug concentrations (DDP 7μM, TX1-85-1 3 and 7μM), which was the best combination values in all CI data, for further apoptosis analysis. The percentages of apoptosis cells treated with DDP and TX1-85-1 group were significantly higher than those in the cells treated with each single drug (Fig.6F). The same trends were also revealed in OVCAR3-DDP cells (Supplemental Figs. 4D-4F), and dual-treatment of TX1-85-1 and DDP alleviated drug resistance. Hence, these results strongly demonstrated that HER3 promoted DDP resistance in OC cells, and HER3 inhibitor enhanced the drug sensitivity in DDP resistance OC cells.
CYLD knockdown increased DDP resistance by inhibiting apoptosis and inducing ABCB1 expression in vivo.
To investigate whether CYLD knockdown promotes tumor growth and DDP resistance, we established a xenograft model by BALB/C nude mice with CYLD-knockdown cells and control cells: OVCAR3-shNC and OVCAR3-shCYLD (Fig.7A). The results suggested that tumors developed from OVCAR3-shCYLD group grew much faster than those of control group after DDP treatment (Fig.7B-7D). Furthermore, the survival times of the mice with OVCAR3-shCYLD cells were remarkably shorter than those of the mice in the control group (Fig.7E). These findings demonstrated that CYLD knockdown promoted DDP resistance of ovarian tumors with shorter survival time. The tumor samples developed from OVCAR3-shCYLD group showed increasing expression levels of Bcl-XL and ABCB1, but decreasing expression levels of Bax (Fig.7F-7G). Base on Bcl-XL and Bax are vital proteins in apoptosis pathway, these results showed that CYLD knockdown affected the apoptotic pathway and elevated ABCB1 expression which could promote drug resistance through drug efflux in vivo. These findings showed that inhibition of CYLD in ovarian cancer cells induced DDP resistance by inhibiting apoptosis pathway and inducing ABCB1 expression in vivo.
Overexpression of CYLD partially reversed DDP resistance, and dual-treatment with DDP and Verapamil inhibited tumorigenesis.
Based on the results above, we further established a xenograft model with CYLD-overexpressed A2780-DDP cells and its control cells: A2780 and A2780-DDP (Fig.8A). The tumors of mice in CYLD-overexpressed group and A2780 group grew slower than tumors in the A2780-DDP group after DDP treatment (Fig.8B-8C). These findings suggested that overexpressed CYLD partially reversed DDP resistance in OC.
Then we wondered whether ABCB1 inhibitor Verapamil could reverse DDP resistance caused by CYLD knockdown in OC. The synergistic effects of DDP (0, 0.1, 0.3, 0.5, 0.7, 1, 3, 5, 7, 10μM) in combination with Verapamil (0, 1, 3, 7, 10, 30μM) were detected in CYLD knockdown OVCAR3 cells. The results of the CI data of these two drugs were showed in Supplemental Fig. 5A. The synergistic effects of these two drugs in OVCAR3 cells were also described in a Scatter plot (Fig.8D), and the Fa-CI plot which showed these two drugs had synergistic effects (the value was almost below 1) (Fig.8E). According to the CI results, these two drugs’ concentrations (DDP 3μM, Verapamil 3μM) may have the best combination values among all the CI data, and the two concentrations were selected for further apoptosis experiment (Supplemental Fig. 5B). The apoptosis rates of this combinatorial group in CYLD knockdown OVCAR3 cells were significantly higher than those of each single treatment (Fig. 8F). We also obtained similar results in A2780-shCYLD cells (Supplemental Fig. 5C-5F). To further confirm the synergistic effect of DDP and Verapamil, we established a xenograft model with OVCAR3-shCYLD cells. As shown in Fig.8G-8I, treatment with DDP or verapamil alone did not effectively inhibit tumor growth, while dual-inhibitor treatment greatly inhibited tumor growth. In addition, the mice body weights were similar in all 4 groups, suggesting that the combination regimen at the indicated dose did not cause any toxicity in mice (Fig.8J). In conclusion, these findings demonstrated that CYLD is an essential DDP sensitizer biomarker in ovarian cancer, and the HER3 inhibitor and ABCB1 inhibitor combination greatly increased cisplatin sensitivity in resistant OC cells (Fig.8K).