A Histone Deacetylase Inhibitor Suppresses Epithelial-Mesenchymal Transition and Attenuates Chemoresistance in Biliary Tract Cancer

Abstract

Epithelial-mesenchymal transition (EMT) is involved in the characteristics of malignancy, such as invasion, metastasis, and chemoresistance. In biliary tract cancer (BTC), EMT is induced by transforming growth factor-beta 1 (TGF-β1). The EMT is reversible; therefore, it is conceivable that it could be related to some epigenetic changes. We focused on histone deacetylase (HDAC) inhibitors as regulators of TGF-β1 signaling, and investigated their effect on EMT and chemoresistance. We employed four BTC cell lines (MzChA-1, gemcitabine-resistant MzChA-1, TFK-1, and gemcitabine-resistant TFK-1) and used vorinostat as the HDAC inhibitor. The relative mRNA expression of an epithelial marker (CDH1) and mesenchymal markers (CDH2, vimentin, SNAI1) were measured by qRT-PCR to evaluate factors associated with EMT. MTT (3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay was performed to evaluate the chemoresistance of each cell line. In addition, NOD/SCID mice were used to evaluate the effect of vorinostat in vivo. In the parent MzChA-1 and TFK-1 cell lines, TGF-β1 induced EMT and chemoresistance; while vorinostat inhibited the EMT and chemoresistance induced by TGF-β1. In gemcitabine-resistant cell lines that highly expressed TGF-β1, vorinostat inhibited EMT and attenuated chemoresistance. We showed that vorinostat inhibits nuclear translocation of SMAD4 which is a signaling factor of TGF-β1, and this is one of the mechanisms by which vorinostat regulates EMT. We also showed that vorinostat attenuates the binding affinity of SMAD4 to the CDH1-related transcription factors SNAI1, SNAI2, ZEB1, ZEB2, and TWIST. Furthermore, combination therapy with vorinostat and gemcitabine improved survival time in the mice xenografted with gemcitabine resistant MzChA-1 cells. In conclusion, vorinostat regulated TGF-β1-induced EMT and chemoresistance through inhibition of SMAD4 nuclear translocation.

About the Authors:

Takuya Sakamoto

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Shogo Kobayashi

Affiliations Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan, Department of Surgery, Osaka Medical Center for Cancer and Cardiovascular diseases, Nakamichi 1-3-3, Higashinari-ku, Osaka, Osaka 537–8511, Japan

Daisaku Yamada

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Hiroaki Nagano

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Akira Tomokuni

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Yoshito Tomimaru

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Takehiro Noda

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Kunihito Gotoh

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Tadafumi Asaoka

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Hiroshi Wada

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Koichi Kawamoto

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Shigeru Marubashi

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Hidetoshi Eguchi

* E-mail: [email protected]

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Yuichiro Doki

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Masaki Mori

Affiliation: Department of Surgery, Graduate School of Medicine, Osaka University, Yamadaoka 2–2 (E2), Suita, Osaka 565–0871, Japan

Introduction

Biliary tract cancer (BTC), which has an increasing incidence worldwide, has a poor prognosis because it is difficult to diagnose in the early stages. The 5-year survival rate is less than 30% [1–3]. At diagnosis, many patients have unresectable BTC due to local disease spread and/or metastases. Gemcitabine (GEM)-based chemotherapy is often used; however, the antitumor effect is insufficient, and the median survival time is limited to approximately 10 months [4,5]. To improve the prognosis of BTC patients, the molecular biology of the disease has been studied, including the inflammatory cytokines, epithelial-mesenchymal transition (EMT), DNA-repair system, cell signaling, and apoptosis [6–9]. We previously showed that EMT was observed at the invasion front and in regional lymph node metastases with inflammatory cytokine expression in BTC [10]. Several kinds of cytokines reportedly induce EMT, which is related to malignant processes such as invasion, metastasis, and chemoresistance (S1 Table) [11–15]. Transforming growth factor-beta 1 (TGF-β1) induces EMT and chemoresistance, and chemoresistant BTC cells produce IL-6 and TGF-β1 [10]; therefore, we considered that these cytokines could play an important role in promoting EMT in BTC.

Consecutive rounds of EMT and mesenchymal-epithelial transition (MET) occur during tumor development, suggesting that EMT could be reversible and associated with epigenetic changes [16]. We focused on histone deacetylase (HDAC) as a factor related to epigenetic changes. In our microarray data, the mRNA levels of class I HDACs (HDAC1, HDAC2, HDAC3, and HDAC8) were elevated in GEM-resistant BTC cells. Moreover, the HDAC activity was increased by exposure to TGF-β1. These findings suggest that HDACs could be associated with EMT, and inhibition of HDACs may suppress EMT through effects on TGF-β1 signaling. We also examined SMAD4, the key signaling factor of TGF-β1 for EMT, and found it was directly stimulated by TGF-β1 signaling, with translocation into the nucleus where it upregulates EMT-related transcription factors such as SNAI1, SNAI2, ZEB1, ZEB2, and TWIST [17–21]. Additionally, in normal hepatocytes and mesothelial cells, inhibition of HDACs attenuated the nuclear translocation of SMAD4 [22,23]. In this study, we hypothesized that HDAC inhibitors might regulate EMT through the TGF-β1 signaling pathway in BTC cells.

The aim of this study was to investigate the effect of the HDAC inhibitor vorinostat on the TGF-β1 signaling pathway and chemoresistance, which is related to EMT in BTC cells. First, we investigated the effect of vorinostat on TGF-β1-induced EMT. Then, we investigated the effect of vorinostat on chemoresistant BTC cells. Subsequently, we showed the influence of vorinostat on SMAD4 nuclear translocation in BTC cells. Finally, we evaluated the efficacy of vorinostat combined with GEM for BTC with using mice xenografted with MzChA-1 or GEM-resistant MzChA-1 cells. These results showed the effectiveness of a HDAC inhibitor in regulating EMT and chemoresistance in vitro and in vivo. We also identified one of the mechanisms by which a HDAC inhibitor suppresses EMT and attenuates chemoresistance.

Materials and Methods

Cell lines, cultures, and drugs

We used the human BTC cell lines MzChA-1 and TFK-1. The MzChA-1 cells were kindly provided by Prof. Gregory J. Gores of the Mayo Clinic, Rochester, Minnesota, United States of America (USA) [24–27]. The MzChA-1 cells were propagated in Dulbecco’s modified Eagle’s medium.

The TFK-1 cell line was obtained from the Riken BioResource Center in Japan. These cells were propagated in RPMI 1640 medium. Each medium was supplemented with 10% fetal bovine serum and 1% penicillin–streptomycin. All cell lines were incubated at 37°C in a humidified incubator with 5% CO₂.

All cell lines were treated with 5 ng/ml of recombinant human TGF-β1 (Pepro Tech Inc., Rocky Hill, NJ, USA), 100 nM (MzChA-1 and GEM-resistant MzChA-1) or 300 nM (TFK-1 and GEM-resistant TFK-1) of the histone deacetylase inhibitor vorinostat (Selleckchem, Houston, TX, USA), or dimethyl sulfoxide (as a control) for 72 h. In the experiment with TGF-β1, each medium was changed to serum-free medium 24 h after passage to preclude the influence of TGF-β1 in serum. The TGF-β1 and vorinostat were added to each cell line after the change of medium. The antibodies described in S2 Table were used for western blot analysis, immunohistochemistry, and immunocytochemistry.

Establishment of GEM-resistant MzChA-1 (MzChA-1_GR) and TFK-1 (TFK-1_GR) clones

The GEM-resistant MzChA-1 and TFK-1 cell lines were established through exposure to increasing concentrations of GEM (from 0.2 ng/μl to 80 ng/μl for MzChA-1, and from 0.2 ng/μl to 40 ng/μl for TFK-1) over 3 months. After confirming that these cell lines attained more resistance to GEM than the parent cell lines, the MzChA-1_GR and TFK-1_GR cell clones were established by limiting dilution.

Assay of HDAC activity

The MzChA-1 cells were cultured with or without 5 ng/ml TGF-β1 and 300 nM vorinostat for 72 h. Then, nuclear proteins were extracted with the Nuclear Extract Kit (Active Motif, Carlsbad, CA, USA). Subsequently, the activity of HDAC was measured with the HDAC activity assay kit (BioVision, Palo Alto, USA). In this experiment, we used HeLa nuclear extract as a positive control and distilled water as a negative control. All procedures were conducted according to the manufacturer’s recommendations.

Immunocytochemistry

Cells were fixed with 4% paraformaldehyde for 15 min at room temperature, permeabilized with 0.1% Triton X-100, and blocked with 1% BSA for 10 min at room temperature. Subsequently, the cells were stained with the indicated antibodies.

Quantitative reverse transcription-polymerase chain reaction

Total RNA was isolated from cultured cells using Trizol reagent (Invitrogen, Carlsbad, CA, USA). Complementary DNA was synthesized from 3.0 μg of total RNA with the SuperScript first-strand synthesis system (Invitrogen) according to the manufacturer’s protocol. Quantitative real-time polymerase chain reactions (qRT-PCRs) were conducted with the LightCycler-Fast-Start DNA Master SYBR Green I kit (Roche Applied Science, Indianapolis, IN, USA) with gene-specific oligonucleotide primers, as shown in S3 Table. Amplifications were performed in triplicate using the LightCycler System (Roche Applied Science) in accordance with the manufacturer’s protocol. A melting curve analysis was performed to distinguish specific products from non-specific products and primer dimers. The relative expression was calculated as the ratio of specific mRNA to endogenous β-actin (ACTB) mRNA in each sample. All qRT-PCRs were performed in triplicate, and the results were presented as the mean ± standard deviation.

Growth inhibition assays with GEM therapy

Each cell line (5×10³ cells/well) was seeded in a 96-well plate and incubated for 24 h. Then, the medium was changed to serum-free medium. Subsequently, the cells were exposed to GEM (1-120ng/mL) with the indicated drugs (TGF-β1, vorinostat, or dimethyl sulfoxide) for 72h. These assays were repeated at least three times, and similar results were obtained each time. The proportion of MTT (3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide)-positive cells incubated without drugs was defined as 100% viability.

Transfection of small interfering RNA (siRNA)

We used TGF-β siRNA (Invitrogen) and scrambled oligonucleotide siRNA (Sc) as a negative control. We performed siRNA transfection with Lipofectamine RNAiMAX (Invitrogen), according to the manufacturer’s protocols [28]. Transfection efficiency was confirmed by qRT-PCR.

Western blot analysis

Western blot analysis was performed as described previously [28]. Briefly, a whole cell lysate was extracted from BTC cells with RIPA Buffer (Thermo Fisher Scientific, Inc., Rockford, IL, USA) and nuclear proteins were extracted with the EpiQuick Nuclear Extraction Kit I (Epigentek Group Inc., Brooklyn, NY, USA), according to the manufacturer’s protocol. Aliquots (12 μg) of total or nuclear proteins were electrophoresed on sodium dodecyl sulfate-polyacrylamide gels containing 10% Tris–HCl (Bio-Rad Laboratories Inc., Hercules, CA, USA). The separated proteins were transferred to polyvinylidene difluoride membranes (Millipore) and incubated with each primary antibody overnight at 4°C.

Chromatin immunoprecipitation (ChIP)

The ChIP assay was performed with using a ChIP-IT Express Enzymatic kit (Active Motif) as previously described [29]. Briefly, 1 × 10⁷ cells were cross-linked with 1% formaldehyde for 10 min at room temperature and quenched by adding glycine. To harvest chromatin-DNA complexes, the cell lysate was treated with an enzymatic shearing kit (Active Motif), according to the manufacturer's specifications. For immunoprecipitation, 120 μg of chromatin-DNA complexes were incubated with Protein G magnetic beads linked to anti-SMAD4 antibody (10 μg, Cell Signaling Technology, Danvers, Massachusetts, USA), or anti-IgG antibody (10 μg, Active Motif). Eluates were used as templates for PCR. Equal amounts of anti-IgG or pre-immune chromatin-DNA complex were used as controls. The primer sets used for PCR amplification were: SNAI1: 5′-GCTGTCACACCCGGCACCAAG-3′ and 5′-GGCGGCTTGAAATGCCACGG-3′, SNAI2: 5′-ATGCGTGTGAAGTGCTTAGCATAGT-3′ and 5′-CACTCAGTGCCCAACAGTGTGT-3′, ZEB1: 5′-TTTCGGGAAGTTAAAATGTTTG-3′ and 5′-ATCCTGCTTCATCTGCCTGA-3′, ZEB2: 5′-TACGCCTGCGCTGTGACCTA-3′ and 5′-ACTCACTGGACCCGCCTCAG-3′, and TWIST: 5′-AGTCTCCTCCGACCGCTTCCTG-3′ and 5′-CTCCGTGCAGGCGGAAAGTTTGG-3′.

Mouse xenograft tumor model

The MzChA-1 or MzChA-1_GR cells (1.0 × 10⁶ cells per mouse) were injected into the subcutaneous tissue on the back of five- to six-week-old male nude mice (CLEA Japan Inc., Tokyo, Japan), and tumors were allowed to develop. When the tumors grew to the volume of 200 mm³, GEM (125 mg/kg, once a week) with or without vorinostat (60 mg/kg, 5 consecutive days per week) was administered by intraperitoneal injection. To investigate the effect of vorinostat on survival time, treatment was continued until appropriate euthanasia. Resected specimens of tumors were used for immunohistochemistry to investigate the expression of SMAD4.

Immunohistochemistry

Immunohistochemical studies on 10 resected tumor specimens from the mice xenografted with MzChA-1_GR cells. Briefly, formalin-fixed, paraffin-embedded tissues were deparaffinized, boiled for antigen retrieval, incubated with a SMAD4 antibody for 1 hour at room temperature, then visualized with avidin–biotin complex reagents (Vector Laboratory Inc., Burlingame, CA, USA) and diaminobenzidine. All sections were counterstained with hematoxylin. All immunohistochemical studies were evaluated in serial sections with each antibody.

Ethical statements

All animal experiments were carried out in strict accordance with the Institutional Animal Care and Use Committee at the Osaka University (approval number: 20–087), and studies were conducted in compliance with institutional guidelines. In the animal experiments, all invasive treatments such as xenograft transplantation, gavage, and euthanasia were performed under general anesthesia with sevoflurane. The mice were euthanized at the onset of clinical signs such as severe cachexia, significant weight loss, or inactivity.

Statistical analysis

The data are presented as the mean ± standard deviation (SD) from at least three independent experiments. Statistical analysis was performed using Student’s t-test or Fisher’s exact test for categorical data. The unpaired Student’s t-test was used to examine differences in growth inhibitory effects in vitro. P-values < 0.05 were considered statistically significant.

Results

Establishment of Mz-ChA-1_GR and TFK-1_GR cell lines

We established the GEM-resistant cell lines, MzChA-1_GR and TFK-1_GR, to investigate the effect of HDAC inhibitors on chemoresistance. The MzChA-1_GR cells were resistant to GEM (IC₅₀ for GEM > 100 ng/mL, P < 0.01 versus the parent MzChA-1 cells, S1A Fig) and TFK-1_GR cells were resistant to GEM (IC₅₀ for GEM > 100 ng/mL, P < 0.01 versus the parent TFK-1 cells, S1B Fig). These GEM-resistant cells had spindle-like shapes.

Microarray analysis

A microarray analysis was performed using the TORAY 3D-Gene^® to compare the expression profiles of MzChA-1 cells with those of MzChA-1_GR cells and investigate the factors related to epigenetics. The results showed that the expression of class I HDACs (HDAC-1, HDAC-2, HDAC-3, and HDAC-8) were higher in MzChA-1_GR cells than parent MzChA-1 cells (Table 1). Vorinostat inhibits these class I HDACs, so we considered that vorinostat might be effective in chemoresistant BTC cells.

[Figure omitted. See PDF.]

Table 1. The expression of class I HDACs in MzChA-1 (parent) and MzChA-1_GR cells, determined by microarray analysis.

https://doi.org/10.1371/journal.pone.0145985.t001

The mRNA expression of class I HDACs and HDAC activity in MzChA-1 and MzChA-1_GR cells

To validate the microarray analysis results, qRT-PCR for class I HDACs was performed for MzChA-1 and MzChA-1_GR cells. All class I HDACs showed higher expression in MzChA-1_GR cells compared with MzChA-1 cells (P<0.05, Fig 1A). The HDAC activity was significantly higher in MzChA-1 cells exposed to TGF-β1 for 72 h and MzChA-1_GR cells compared with parent MzChA-1 cells (P<0.05, Fig 1B). In addition, vorinostat suppressed the HDAC activity in both MzChA-1 cells exposed to TGF-β1 and MzChA-1_GR cells (P<0.05, Fig 1C).

[Figure omitted. See PDF.]

Fig 1. The mRNA expression of class I HDACs and HDAC activity in MzChA-1 and MzChA-1_GR cells.

Values represent the mean ± S.D. *P < 0.05. All experiments were conducted at least three times. (A) Comparison of the expression of class I HDACs in MzChA-1 and MzChA-1_GR cells by qRT-PCR. (B) The effect of TGF-β1 on HDAC activity in MzChA-1 cells and the HDAC activity in MzChA-1_GR cells, activity was measured with the HDAC activity assay kit. (C) The changes in HDAC activity in MzChA-1 and MzChA-1_GR cells treated with TGF-β1 and vorinostat. In the experiments for panels (B) and (C), the cells were incubated with or without 5 ng/ml TGF-β1 or 100 nM vorinostat for 72 h.

https://doi.org/10.1371/journal.pone.0145985.g001

The relationship between TGF-β1 expression and chemoresistance

To clarify the significance of TGF-β1 in EMT and chemoresistance, we investigated the expression of TGF-β1 in parent MzChA-1 and MzChA-1_GR cells. The results showed that the expression level of TGF-β1 was significantly higher in MzChA-1_GR cells than parent MzChA-1 cells (Fig 2A). Then, we transfected with siTGF-β to investigate whether EMT and chemoresistance could be suppressed by TGF-β knock down in MzChA-1_GR cells. The result showed that MzChA-1_GR cells transfected with siTGF-β had increased the mRNA expression of CDH1 and decreased the mRNA expression of CDH2, vimentin, and SNAI1 (Fig 2B). Moreover, transfection with siTGF-β attenuated chemoresistance in MzChA-1_GR cells (Fig 2C). These results suggested the significance of TGF-β in EMT and chemoresistance in BTC.

[Figure omitted. See PDF.]

Fig 2. The effect of TGF-β1 knockdown using TGF-β small interfering RNA (siRNA) on EMT and chemoresistance in MzChA-1_GR cells.

MzChA-1_GR cells were transfected with scrambled oligonucleotide siRNA (negative control) or TGF-β1 siRNA. All experiments were conducted at least three times. Values represent the mean ± S.D. *P < 0.05. (A) Comparison of the expression of TGF-β1 in MzChA-1 and MzChA-1_GR cells. (B) The changes of EMT-related mRNA expression as a result of TGF-β siRNA transfection in MzChA-1_GR cells. (C) The effect of TGF-β siRNA transfection on chemoresistance in MzChA-1_GR cells. Growth inhibition assays were performed for transfected and non-transfected cells treated with GEM.

https://doi.org/10.1371/journal.pone.0145985.g002

The effect of vorinostat on TGF-β1-induced EMT and chemoresistance

We investigated the effect of vorinostat on TGF-β1-induced EMT and chemoresistance in BTC cells. In MzChA-1 cells, TGF-β1 caused morphological changes of the cells, from valvate-like shapes to spindle-like shapes. Moreover, TGF-β1 down-regulated the expression of CDH1 as determined by immunofluorescence staining (Fig 3A). In contrast, vorinostat prevented these changes caused by TGF-β1 (Fig 3A). Furthermore, TGF-β1 decreased the mRNA expression of CDH1 and increased the mRNA expression of CDH2, vimentin (VIM), and SNAI1. In contrast, vorinostat inhibited the changes in mRNA expression caused by TGF-β1 (P<0.05, Fig 3B). Additionally, vorinostat attenuated GEM resistance induced by TGF-β1 (P<0.05, Fig 3C). In TFK-1 cells, vorinostat prevented morphological changes and the down-regulation of CDH1 caused by TGF-β1 (Fig 3D). Vorinostat also inhibited the changes in mRNA expression related to EMT (P<0.05, Fig 3E), and attenuated chemoresistance to GEM caused by TGF-β1 (P<0.05, Fig 3F). Thus, vorinostat suppressed EMT and attenuated chemoresistance caused by TGF-β1.

[Figure omitted. See PDF.]

Fig 3. Influence of vorinostat on TGF-β1-induced EMT and chemoresistance in BTC cell lines.

In each experiment, the cells were incubated with or without 5 ng/ml of TGF-β1 and 100 nM of vorinostat for MzChA-1 cells and 300 nM of vorinostat for TFK-1 cells for 72 h. Values represent the mean ± S.D. *P < 0.05. Scale bars: 100 μm. All experiments were conducted at least three times. (A) Representative cell morphological changes induced by TGF-β and vorinostat in MzChA-1 cells. Immunofluorescence for CDH1 (red) was performed in MzChA-1 cells. Nuclear staining (blue) was performed with Hoechst. (B) The effect of vorinostat on the EMT-related mRNAs expression in MzChA-1 cells. (C) The effect of vorinostat on chemoresistance induced by TGF-β1 exposure in MzChA-1 cells. The growth inhibition assays were performed for MzChA-1 cells treated with GEM. (D) Representative cell morphological changes induced by TGF-β and vorinostat in TFK-1 cells. Immunofluorescence for CDH1 (red) was performed in TFK-1 cells. Nuclear staining (blue) was performed with Hoechst. (E) The effect of vorinostat on the EMT-related mRNAs expression in TFK-1 cells. (F) The effect of vorinostat on chemoresistance induced by TGF-β1 exposure in TFK-1 cells. Growth inhibition assays were performed for TFK-1 cells treated with GEM.

https://doi.org/10.1371/journal.pone.0145985.g003

The influence of vorinostat on MzChA-1_GR and TFK-1_GR cells

In MzChA-1_GR cells, vorinostat changed the morphology from spindle-like shapes to valvate-like shapes and enhanced the staining of CDH1, as determined by immunofluorescence staining (Fig 4A). The CDH1 mRNA expression tended to increase, and CDH2, VIM, and SNAI1 were significantly decreased by vorinostat (P<0.05, Fig 4B). In addition, vorinostat attenuated GEM chemoresistance in MzChA-1_GR cells (P<0.05, Fig 4C).

[Figure omitted. See PDF.]

Fig 4. Influence of vorinostat on chemoresistant BTC cells.

In each experiment, MzChA-1_GR cells were incubated with or without 100 nM vorinostat and TFK-1_GR cells with or without 300nM vorinostat for 72 h. Values represent the mean ± S.D. *P < 0.05. Scale bars: 100 μm. All experiments were conducted at least three times. (A) Representative morphological change induced by vorinostat in MzChA-1_GR cells. Immunofluorescence for CDH1 (red) was performed in MzChA-1 cells. Nuclear staining (blue) was performed with Hoechst. (B) The effect of vorinostat on the EMT-related mRNAs expression in MzChA-1_GR cells. (C) The effect of vorinostat on chemoresistance in MzChA-1_GRcells. Growth inhibition assays were performed for MzChA-1_GR cells treated with GEM. (D) Representative morphological changes induced by vorinostat in TFK-1_GR cells. Immunofluorescence for CDH1 (red) was performed in TFK-1_GR cells. Nuclear staining (blue) was performed with Hoechst. (E) The effect of vorinostat on the EMT-related mRNAs expression in TFK-1_GR cells. (F) The effect of vorinostat on chemoresistance in TFK-1_GR cells. Growth inhibition assays were performed for TFK-1_GR cells treated with GEM.

https://doi.org/10.1371/journal.pone.0145985.g004

In TFK-1_GR cells, vorinostat changed the morphology and enhanced the staining of CDH1, as measured by immunofluorescence staining similar to MzChA-1_GR cells (Fig 4D). In addition, vorinostat increased the expression of CDH1 mRNA, and decreased the mRNA expression of CDH2, VIM, and SNAI1 (P<0.05, Fig 4E). GEM resistance was also attenuated by vorinostat in TFK-1_GR cells (P<0.05, Fig 4F). These results suggested that vorinostat suppressed EMT and attenuated chemoresistance in GEM resistant BTC cells.

The effect of vorinostat on SMAD4 nuclear translocation in MzChA-1 cells

In our previous study, SMAD4 was the key molecule regulating TGF-β1-induced EMT [10]. Therefore, we investigated the effect of vorinostat on the expression of SMAD4 by western blot analysis. In whole cell lysates, there was no significant difference in the expression of SMAD4 protein between the TGF-β1 group and the TGF-β1 plus vorinostat group (Fig 5A, left). However, TGF-β1 upregulated the expression of the SMAD4 nucleoprotein compared with the control, and vorinostat inhibited this upregulation of the SMAD4 nucleoprotein caused by TGF-β1 (Fig 5A, right). Furthermore, we used the ChIP method to analyze the molecular mechanism of SMAD4. While TGF-β1 enhanced the binding affinity of SMAD4 to the CDH1-related transcriptional factors SNAI1, SNAI2, ZEB1, ZEB2, and TWIST, vorinostat inhibited the TGF-β1-enhanced binding affinity of SMAD4 to these transcriptional factors (Fig 5B).

[Figure omitted. See PDF.]

Fig 5. Inhibition of the nuclear translocation of SMAD4 by vorinostat in MzChA-1 cells.

In each experiment, cells were incubated with or without 5 ng/ml of TGF-β1 and 100 nM of vorinostat for 72 h. Scale bars: 100 μm. All experiments were conducted at least three times. (A) The change in SMAD4 expression in whole cell lysates (left panel) and in the nucleus (right panel) caused by TGF-β and vorinostat. (B) ChIP assay showing expression of SNAI1, SNAI2, ZEB1, ZEB2, and TWIST, which are regulatory elements of CDH1. (C) Immunofluorescence of SMAD4 (green) was performed in MzChA-1 cells. The effect of TGF-β1 and vorinostat on SMAD4 nuclear translocation were investigated. Nuclear staining (blue) was performed with Hoechst.

https://doi.org/10.1371/journal.pone.0145985.g005

Immunocytochemistry demonstrated that SMAD4 stained strongly in the nuclei after treatment with TGF-β1 (Fig 5C). In contrast, SMAD4 was weakly stained in the nuclei after treatment with TGF-β1 plus vorinostat (Fig 5C). Based on these results, vorinostat regulated EMT through inhibition of SMAD4 nuclear translocation.

The effect of vorinostat on the TGF-β signaling factors, SMAD2, SMAD3, and Jun N-terminal kinase (JNK)

We also investigated the effect of vorinostat on SMAD2 and SMAD3 because they are closely related to SMAD4 nuclear translocation in TGF-β signaling. We found TGF-β1 induced SMAD2 and SMAD3 phosphorylation in both whole cell lysates and nuclear proteins (Fig 6). In addition, vorinostat inhibited TGF-β1-induced nuclear translocation of phosphorylated SMAD2 (p-SMAD2) and SMAD3 (p-SMAD3); while, vorinostat had no effect on the nuclear translocation of SMAD2 and SMAD3 (Fig 6). These results suggest that vorinostat can inhibit nuclear translocation of p-SMAD2, p-SMAD3, and SMAD4 at the same time. We also investigated the activation of JNK because recent reports showed that JNK activation plays a central role in TGF-β mediated EMT. The result showed TGF-β1 upregulated the expression of phosphorylated JNK; however, vorinostat had no effect on JNK (Fig 6). Accordingly, the effect of vorinostat on EMT seems to be independent of the JNK pathway in BTC.

[Figure omitted. See PDF.]

Fig 6. The effect of TGF-β1 and vorinostat on SMAD2, SMAD3, and JNK in MzChA-1.

In each experiment, cells were incubated with or without 5 ng/ml of TGF-β1 and 100 nM of vorinostat for 72 h. All experiments were conducted at least three times. The changes in SMAD2, p-SMAD2, SMAD3, and p-SMAD3 expression in whole cell lysates (left panel) and in the nucleus (right panel) are shown. The expression of JNK and p-JNK in whole cell lysates are also shown (left panel).

https://doi.org/10.1371/journal.pone.0145985.g006

The effect of vorinostat on the mouse xenografted tumor model

To confirm the in vitro results, we evaluated the effect of vorinostat on mice xenografted with MzChA-1 or MzChA-1_GR cells. In the mice xenografted with MzChA-1 cells, combination therapy with vorinostat and GEM did not improve the survival time compared with a single administration of GEM (Fig 7A). In contrast, combination therapy with vorinostat and GEM improved the survival time compared with a single administration of GEM in the mice xenografted with MzChA-1_GR cells (log-rank test, P<0.01, Fig 7B).

[Figure omitted. See PDF.]

Fig 7. In vivo analysis of mice xenografted with MzChA-1 and MzChA-1_GR cells.

NOD/SCID mice xenografted with MzChA-1 and MzChA-1_GR cells were employed as in vivo models. When the tumor volume was over 200 mm³, the mice were treated with the indicated drugs (GEM [125 mg/kg, once a week] and/or vorinostat [60 mg/kg, 5 consecutive days per week]). The prognostic value was evaluated by the Kaplan–Meier method and a log-rank test. (A) Survival curve for mice xenografted with MzChA-1 cells. (B) Survival curve for mice xenografted with MzChA-1_GR cells. (C) Immunohistochemistry for SMAD4 on the resected tumor specimens of the mice xenografted with MzChA-1_GR cells.

https://doi.org/10.1371/journal.pone.0145985.g007

Subsequently, we performed immunohistochemistry on the resected tumor specimens of the mice to evaluate the expression of SMAD4 in vivo. The result showed that SMAD4 was strongly stained in the nuclei of tumor cells when the mice were treated with a single administration of GEM. In contrast, SMAD4 was weakly stained in the nuclei of tumor cells when the mice were treated with GEM and vorinostat (Fig 7C). These results were the same in both primary and metastatic sites (pulmonary). Thus, vorinostat improved survival time in the mice xenografted with chemoresistant cells when they were treated with GEM.

Discussion

In many cancers, several genes related to tumorigenesis and tumor progression are upregulated, and genes related to tumor suppression are downregulated. One of the factors that regulate this gene expression is posttranslational modification of the histones. Methylation and acetylation of histones are well known post-translational modifications [30–32]. Recently, much attention has been paid to the methylation and acetylation of histones for the treatment of cancers [33,34]. In this study, we focused on the ability of HDAC inhibitors to regulate EMT. We used vorinostat, which inhibits class I HDACs and is often used as a HDAC inhibitor in cancer clinical trials, as shown in S4 Table [35–39]. In several cancers, such as non-small cell lung cancer and breast cancer, combination chemotherapy with vorinostat was safe and showed potential for improving the efficacy [38,40].

We focused on SMADs, particularly SMAD4 which is a signaling factor of TGF-β1, as one of the mechanisms that vorinostat regulates EMT. It is often reported that TGF-β induces EMT and chemoresistance [10]. Activated type I TGF-β1 receptors phosphorylate SMAD2 and SMAD3, which form a heteromeric complex with SMAD4 [41,42]. This complex translocates to the nucleus and targets a variety of DNA binding proteins to regulate transcriptional responses [41,42]. Thus, transcriptional regulators of CDH1, such as SNAI, ZEB, and TWIST, are upregulated [42]. Then, CDH1 is downregulated and EMT occurs [43,44]. Thus, inhibiting nuclear translocation of SMAD2, SMAD3, and SMAD4 complexes could regulate TGF-β1-related EMT. First, we investigated whether vorinostat changes the mRNA expression of these SMADs by qRT-PCR. No changes in these SMADs were observed after exposure to vorinostat (S2 Fig). However, a few studies have shown aneffect of HDAC inhibitors on SMAD4 [22,23]. Kaimori et al. showed that histone deacetylase inhibition suppresses the TGF-β1-induced EMT in hepatocytes via inhibition of SMAD4 nuclear translocation [22]. Chung et al. showed that the HDAC inhibitor, m-carboxycinnamic acid, inhibited SMAD4 nuclear translocation in human pleural mesothelial cells [23]. However, there are no studies on the effect of HDAC inhibitors on SMAD4 nuclear translocation in any cancer. Our investigation showed that vorinostat inhibited SMAD4 nuclear translocation induced by TGF-β1 in BTC. We also showed that vorinostat inhibited nuclear translocation of p-SMAD2 and p-SMAD3. These results suggest that vorinostat inhibits nuclear translocation of the p-SMAD2, p-SMAD3, and SMAD4 complex.

To clarify whether vorinostat sufficiently inhibited EMT related to TGF-β signaling, we investigated the effect of vorinostat and siTGF-β transfection in MzChA-1_GR cells. We found vorinostat with siTGF-β transfection had no additional effect on EMT and chemoresistance (data not shown). This result suggests that the effect of vorinostat is similar to transfection with siTGF-β in regulating EMT and chemoresistance.

In vivo experiments showed that vorinostat with GEM treatment improved the survival time of the mice xenografted with MzChA-1_GR cells compared with single administration of GEM. However, there was no significant change in tumor development between the vorinostat with GEM treatment group and GEM treatment group (S3 Fig). One of the reasons why vorinostat could improve the survival rate might be that vorinostat suppressed tumor metastasis. In fact, the tumors in metastatic sites tended to be large in the GEM treatment group compared with the vorinostat with GEM treatment group. Immunohistochemistry showed that SMAD4 was more strongly stained in the nuclei of cancer cells in the GEM treatment group compared with the vorinostat with GEM treatment group. These results suggest that vorinostat regulates tumor progression and metastasis by inhibiting nuclear translocation of SMAD4 in vivo.

Some studies reported that HDAC inhibitors have a similar function as TGF-β via activation of CDK inhibitors. One of the representative CDK inhibitors related to EMT is p21, and TGF-β induces p21 through a p53-independent mechanism [45]. The HDAC inhibitors are reported to activate p21 in breast cancer cells [46]. Several studies showed that p21 inhibited EMT [47]. Therefore, we investigated the change in p21 expression by TGF-β1 and vorinostat. The result showed that only TGF-β1 induced p21 expression (S4 Fig). In addition, another CDK inhibitor, p15, which is reported to be induced by TGF-β1 was similarly induced (S4 Fig).These results indicate the effect of vorinostat in regulating the activation of CDK inhibitors might differ depending on the kind of cancer. This issue is controversial and more investigation would be needed in other types of cancer.

We used vorinostat in this study; however, other HDAC inhibitors, such as trichostatin A, panobinostat, entinostat, and valproic acid have been used to regulate EMT, and some of them have been used in clinical trials [48–51]. In the future, these drugs will need to be tested to determine if they regulate EMT and chemoresistance, as suggested by the results of this study. In conclusion, the HDAC inhibitor, vorinostat, regulated EMT and chemoresistance in BTC via inhibition of nuclear translocation of SMAD4 in vitro and in vivo.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. Growth inhibition assays for GEM in parent cell lines and GEM-resistant cell lines.

Comparison of chemoresistance in parent cell lines and established GEM-resistant cell lines by growth inhibition assays. Values represent the mean ± S.D. *P < 0.05. (A) Growth inhibition assay in MzChA-1 and MzChA-1_GR cells. (B) Growth inhibition assay in TFK-1 and TFK-1_GR cells.

https://doi.org/10.1371/journal.pone.0145985.s001

(TIF)

S2 Fig. The effect of TGF-β and vorinostat on the SMAD2, SMAD3, and SMAD4 mRNAs expression in MzChA-1 cells.

Each cell line was cultured with or without 5 ng/ml of TGF-β1 and 100 nM of vorinostat for 72 h. Then, SMAD2, SMAD3, and SMAD4 were assessed by qRT-PCR. Values represent the mean ± S.D. All experiments were conducted at least three times.

https://doi.org/10.1371/journal.pone.0145985.s002

(TIF)

S3 Fig. Tumor progression in the mice xenografted with MzChA-1_GR cells.

Tumor progression was assessed by tumor volume. Values represent the mean ± S.D.

https://doi.org/10.1371/journal.pone.0145985.s003

(TIF)

S4 Fig. The effect of TGF-β and vorinostat on the expression of p15 and p21 in MzChA-1 cells.

Each cell line was cultured with or without 5 ng/ml of TGF-β1 and 100 nM of vorinostat for 72 h. Then, p15 and p21 were assessed by qRT-PCR. Values represent the mean ± S.D. *P < 0.05. All experiments were conducted at least three times.

https://doi.org/10.1371/journal.pone.0145985.s004

(TIF)

S1 Table. Cytokines associated with EMT (Table A). References for S1A Table (Table B).

https://doi.org/10.1371/journal.pone.0145985.s005

(DOCX)

S2 Table. Antibodies used in western blot analysis, immunohistochemistry, and immunocytochemistry.

https://doi.org/10.1371/journal.pone.0145985.s006

(DOCX)

S3 Table. Primers used in qRT-PCR.

https://doi.org/10.1371/journal.pone.0145985.s007

(DOCX)

S4 Table. Clinical trials involving HDAC inhibitors (Table A). References for S4A Table (Table B).

https://doi.org/10.1371/journal.pone.0145985.s008

(DOCX)

Acknowledgments

We thank the members of our laboratory for technical advice and helpful discussion.

Author Contributions

Conceived and designed the experiments: MM YD HE HN SK. Performed the experiments: TS DY AT YT TN KG TA HW KK SM. Analyzed the data: TS SK DY. Contributed reagents/materials/analysis tools: DY HW SM. Wrote the paper: TS SK DY.

Citation: Sakamoto T, Kobayashi S, Yamada D, Nagano H, Tomokuni A, Tomimaru Y, et al. (2016) A Histone Deacetylase Inhibitor Suppresses Epithelial-Mesenchymal Transition and Attenuates Chemoresistance in Biliary Tract Cancer. PLoS ONE 11(1): e0145985. https://doi.org/10.1371/journal.pone.0145985

References

1. Jang JY, Kim SW, Park DJ, Ahn YJ, Yoon YS, Choi MG, et al. Actual long-term outcome of extrahepatic bile duct cancer after surgical resection. Ann Surg. 2005; 241: 77–84. pmid:15621994

2. Murakami S, Ajiki T, Okazaki T, Ueno K, Kido M, Matsumoto I, et al. Factors affecting survival after resection of intrahepatic cholangiocarcinoma. Surg Today. 2014; 44: 1847–54. pmid:24452507

3. Kim HJ, Kim CY, Hur YH, Koh YS, Kim JC, Kim HJ, et al. The prognostic factors for survival after curative resection of distal cholangiocarcinoma: perineural invasion and lymphovascular invasion. Surg Today. 2014; 44: 1879–86. pmid:24535697

4. Kobayashi S, Nagano H, Marubashi S, Wada H, Eguchi H, Takeda Y, et al. Treatment of borderline cases for curative resection of biliary tract cancer. J Surg Oncol. 2011; 104: 499–503. pmid:21538362

5. Valle J, Wasan H, Palmer DH, Cunningham D, Anthoney A, Maraveyas A, et al. Cisplatin plus gemcitabine versus gemcitabine for biliary tract cancer. N Engl J Med. 2010; 362: 1273–81. pmid:20375404

6. Kobayashi S, Werneburg NW, Bronk SF, Kaufmann SH, Gores GJ. Interleukin-6 contributes to Mcl-1 up-regulation and TRAIL resistance via an Akt-signaling pathway in cholangiocarcinoma cells. Gastroenterology. 2005; 128: 2054–65. pmid:15940637

7. Nitta T, Mitsuhashi T, Hatanaka Y, Miyamoto M, Oba K, Tsuchikawa T, et al. Prognostic significance of epithelial-mesenchymal transition-related markers in extrahepatic cholangiocarcinoma: comprehensive immunohistochemical study using a tissue microarray. Br J Cancer. 2014; 111: 1363–72. pmid:25077440

8. Nakashima S, Kobayashi S, Nagano H, Tomokuni A, Tomimaru Y, Asaoka T, et al. BRCA/Fanconi Anemia Pathway Implicates Chemo-resistance to Gemcitabine in Biliary Tract Cancer. Cancer Sci. 2015

9. Mott JL, Kobayashi S, Bronk SF, Gores GJ. mir-29 regulates Mcl-1 protein expression and apoptosis. Oncogene. 2007; 26: 6133–40. pmid:17404574

10. Yamada D, Kobayashi S, Wada H, Kawamoto K, Marubashi S, Eguchi H, et al. Role of crosstalk between interleukin-6 and transforming growth factor-beta 1 in epithelial-mesenchymal transition and chemoresistance in biliary tract cancer. Eur J Cancer. 2013; 49: 1725–40. pmid:23298711

11. Zavadil J, Böttinger EP. TGF-bold italic beta and epithelial-to-mesenchymal transitions. Oncogene. 2005; 24: 5764–74. pmid:16123809

12. Gonzalez DM, Medici D. Signaling mechanisms of the epithelial-mesenchymal transition. Sci Signal. 2014 7: re8. pmid:25249658

13. Farrell J, Kelly C, Rauch J, Kida K, García-Muñoz A, Monsefi N, et al. HGF induces epithelial-to-mesenchymal transition by modulating the mammalian hippo/MST2 and ISG15 pathways. J Proteome Res. 2014; 13: 2874–86. pmid:24766643

14. Rokavec M, Öner MG, Li H, Jackstadt R, Jiang L, Lodygin D, et al. IL-6R/STAT3/miR-34a feedback loop promotes EMT-mediated colorectal cancer invasion and metastasis. J Clin Invest. 2014; 124: 1853–67. pmid:24642471

15. Ogunwobi OO, Puszyk W, Dong HJ, Liu C. Epigenetic upregulation of HGF and c-Met drives metastasis in hepatocellular carcinoma. PLoS One. 2013; 8: e63765. pmid:23723997

16. Lamouille S, Xu J, Derynck R. Molecular mechanisms of epithelial-mesenchymal transition. Nat Rev Mol Cell Biol. 2014; 15: 178–96. pmid:24556840

17. Nieto MA. The snail superfamily of zinc-finger transcription factors. Nat Rev Mol Cell Biol. 2002; 3: 155–66. pmid:11994736

18. Chanmee T, Ontong P, Mochizuki N, Kongtawelert P, Konno K, Itano N. Excessive hyaluronan production promotes acquisition of cancer stem cell signatures through the coordinated regulation of Twist and the transforming growth factor β (TGF-β)-Snail signaling axis. J Biol Chem. 2014; 289: 26038–56. pmid:25077968

19. Hotz B, Arndt M, Dullat S, Bhargava S, Buhr HJ, Hotz HG. Epithelial to mesenchymal transition: expression of the regulators snail, slug, and twist in pancreatic cancer. Clin Cancer Res. 2007; 13: 4769–76. pmid:17699854

20. Kurahara H, Takao S, Maemura K, Mataki Y, Kuwahata T, Maeda K, et al. Epithelial-mesenchymal transition and mesenchymal-epithelial transition via regulation of ZEB-1 and ZEB-2 expression in pancreatic cancer. J Surg Oncol. 2012; 105: 655–61. pmid:22213144

21. von Burstin J, Eser S, Paul MC, Seidler B, Brandl M, Messer M, et al. E-cadherin regulates metastasis of pancreatic cancer in vivo and is suppressed by a SNAIL/HDAC1/HDAC2 repressor complex. Gastroenterology. 2009; 137: 361–71. pmid:19362090

22. Kaimori A, Potter JJ, Choti M, Ding Z, Mezey E, Koteish AA. Histone deacetylase inhibition suppresses the transforming growth factor beta1-induced epithelial-to-mesenchymal transition in hepatocytes. Hepatology. 2010; 52: 1033–45. pmid:20564330

23. Chung CL, Sheu JR, Chen WL, Chou YC, Hsiao CJ, Hsiao SH, et al. Histone deacetylase inhibitor m-carboxycinnamic acid bis-hydroxamide attenuates plasminogen activator inhibitor-1 expression in human pleural mesothelial cells. Am J Respir Cell Mol Biol. 2012; 46: 437–45. pmid:22033265

24. Kobayashi S, Werneburg NW, Bronk SF, Kaufmann SH, Gores GJ, et al. Interleukin-6 contributes to Mcl-1 up-regulation and TRAIL resistance via an Akt-signaling pathway in cholangiocarcinoma cells. Gastroenterology. 2005; 128: 2054–65. pmid:15940637

25. Shimizu Y, Demetris AJ, Gollin SM, Storto PD, Bedford HM, Altarac S, et al. Two new human cholangiocarcinoma cell lines and their cytogenetics and responses to growth factors, hormones, cytokines or immunologic effector cells. Int J Cancer. 1992; 52: 252–60 pmid:1355757

26. Murakami T, Yano H, Maruiwa M, Sugihara S, Kojiro M. Establishment and characterization of a human combined hepatocholangiocarcinoma cell line and its heterologous transplantation in nude mice. Hepatology. 1987; 7: 551–6 pmid:3032760

27. Knuth A, Gabbert H, Dippold W, Klein O, Sachsse W, Bitter-Suermann D, et al. Biliary adenocarcinoma. Characterisation of three new human tumor cell lines. J Hepatol. 1985; 1: 579–96 pmid:4056357

28. Yamada D, Kobayashi S, Yamamoto H, Tomimaru Y, Noda T, Uemura M, et al. Role of the hypoxia-related gene, JMJD1A, in hepatocellular carcinoma: clinical impact on recurrence after hepatic resection. Ann Surg Oncol. 2012; 19 Suppl 3: S355–64 pmid:21607773

29. Miyoshi N, Ishii H, Nagai K, Hoshino H, Mimori K, Tanaka F, et al. Defined factors induce reprogramming of gastrointestinal cancer cells. Proc Natl Acad Sci USA. 2010; 107: 40–5 pmid:20018687

30. Imhof A, Wolffe AP. Transcription: gene control by targeted histone acetylation. Curr Biol. 1998; 8: R422–4. pmid:9637914

31. Jenuwein T, Allis CD. Translating the histone code. Science. 2001; 293: 1074–80. pmid:11498575

32. Irvine RA, Lin IG, Hsieh CL. DNA methylation has a local effect on transcription and histone acetylation. Mol Cell Biol. 2002; 22: 6689–96. pmid:12215526

33. Sui X, Zhu J, Zhou J, Wang X, Li D, Han W, et al. Epigenetic modifications as regulatory elements of autophagy in cancer. Cancer Lett. 2015; 360: 106–13. pmid:25687886

34. Rodrigues MF, Carvalho É, Pezzuto P, Rumjanek FD, Amoêdo ND. Reciprocal modulation of histone deacetylase inhibitors sodium butyrate and trichostatin a on the energy metabolism of breast cancer cells. J Cell Biochem. 2015; 116: 797–808. pmid:25510910

35. Munster PN, Thurn KT, Thomas S, Raha P, Lacevic M, Miller A, et al. A phase II study of the histone deacetylase inhibitor vorinostat combined with tamoxifen for the treatment of patients with hormone therapy-resistant breast cancer. Br J Cancer. 2011; 104: 1828–35 pmid:21559012

36. Ramaswamy B, Fiskus W, Cohen B, Pellegrino C, Hershman DL, Chuang E, et al. Phase I-II study of vorinostat plus paclitaxel and bevacizumab in metastatic breast cancer: evidence for vorinostat-induced tubulin acetylation and Hsp90 inhibition in vivo. Breast Cancer Res Treat. 2012; 132: 1063–72 pmid:22200869

37. Fakih MG, Groman A, McMahon J, Wilding G, Muindi JR. A randomized phase II study of two doses of vorinostat in combination with 5-FU/LV in patients with refractory colorectal cancer. Cancer Chemother Pharmacol. 2012; 69: 743–51 pmid:22020318

38. Tu Y, Hershman DL, Bhalla K, Fiskus W, Pellegrino CM, Andreopoulou E, et al. A phase I-II study of the histone deacetylase inhibitor vorinostat plus sequential weekly paclitaxel and doxorubicin-cyclophosphamide in locally advanced breast cancer. Breast Cancer Res Treat. 2014;146: 145–52. pmid:24903226

39. Krug LM, Kindler HL, Calvert H, Manegold C, Tsao AS, Fennell D, et al. Vorinostat in patients with advanced malignant pleural mesothelioma who have progressed on previous chemotherapy (VANTAGE-014): a phase 3, double-blind, randomised, placebo-controlled trial. Lancet Oncol. 2015; 16: 447–56 pmid:25800891

40. Reguart N, Rosell R, Cardenal F, Cardona AF, Isla D, Palmero R, et al. Phase I/II trial of vorinostat (SAHA) and erlotinib for non-small cell lung cancer (NSCLC) patients with epidermal growth factor receptor (EGFR) mutations after erlotinib progression. Lung Cancer. 2014; 84: 161–7. pmid:24636848

41. Kawabata M, Miyazono K. Signal transduction of the TGF-beta superfamily by Smad proteins. J Biochem. 1999; 125: 9–16. pmid:9880789

42. Shi Y, Massagué J. Mechanisms of TGF-beta signaling from cell membrane to the nucleus. Cell. 2003; 113: 685–700. pmid:12809600

43. Cieply B, Farris J, Denvir J, Ford HL, Frisch SM, et al. Epithelial-mesenchymal transition and tumor suppression are controlled by a reciprocal feedback loop between ZEB1 and Grainyhead-like-2. Cancer Res. 2013; 73: 6299–309. pmid:23943797

44. Talbot LJ, Bhattacharya SD, Kuo PC. Epithelial-mesenchymal transition, the tumor microenvironment, and metastatic behavior of epithelial malignancies. Int J Biochem Mol Biol. 2012; 3: 117–36. pmid:22773954

45. Datto MB, Li Y, Panus JF, Howe DJ, Xiong Y, Wang XF. Transforming growth factor beta induces the cyclin-dependent kinase inhibitor p21 through a p53-independent mechanism. Proc Natl Acad Sci U S A. 1995; 92: 5545–9 pmid:7777546

46. Bellucci L, Dalvai M, Kocanova S, Moutahir F, Bystricky K. Activation of p21 by HDAC inhibitors requires acetylation of H2A.Z. PLoS One. 2013; 8: e54102 () pmid:23349794

47. Li XL, Hara T, Choi Y, Subramanian M, Francis P, Bilke S, et al. A p21-ZEB1 complex inhibits epithelial-mesenchymal transition through the microRNA 183-96-182 cluster. Mol Cell Biol. 2014; 34: 533–50 pmid:24277930

48. Sun S, Han Y, Liu J, Fang Y, Tian Y, Zhou J, et al. Trichostatin A targets the mitochondrial respiratory chain, increasing mitochondrial reactive oxygen species production to trigger apoptosis in human breast cancer cells. PLoS One. 2014; 9: e91610. pmid:24626188

49. Abaza MS, Bahman AM, Al-Attiyah RJ. Valproic acid, an anti-epileptic drug and a histone deacetylase inhibitor, in combination with proteasome inhibitors exerts antiproliferative, pro-apoptotic and chemosensitizing effects in human colorectal cancer cells: underlying molecular mechanisms. Int J Mol Med. 2014; 34: 513–32. pmid:24899129

50. Cassier PA, Floquet A, Penel N, Derbel O, Bui N'guyen B, Guastalla JP, et al. The histone deacetylase inhibitor panobinostat is active in patients with advanced pretreated ovarian sex-cord tumors. Ann Oncol. 2014; 25: 1074–5. pmid:24651409

51. Ngamphaiboon N, Dy GK, Ma WW, Zhao Y, Reungwetwattana T, DePaolo D, et al. A phase I study of the histone deacetylase (HDAC) inhibitor entinostat, in combination with sorafenib in patients with advanced solid tumors. Invest New Drugs. 2015; 33: 225–32. pmid:25371323