Disruption of ETV6 leads to TWIST1-dependent progression and resistance to epidermal growth factor receptor tyrosine kinase inhibitors in prostate cancer
Molecular Cancer volume 17, Article number: 42 (2018)
ETS variant gene 6 (ETV6) is a putative tumor suppressor and repressed by epidermal growth factor receptor (EGFR) signaling in prostate cancer. Since EGFR antagonists seem ineffective in castration-resistant prostate cancer (CRPC), we aim to study the role of ETV6 in the development of drug resistance.
Etv6 target gene was validated by ChIP and promoter reporter assays. Correlation of ETV6 and TWIST1 was analyzed in human clinical datasets and tissue samples. Migration, invasion, and metastasis assays were used to measure the cellular responses after perturbation of ETV6 -TWIST1 axis. Proliferation and tumor growth in xenograft model were performed to evaluate the drug sensitivities of EGFR-tyrosine kinase inhibitors (TKIs).
ETV6 inhibits TWIST1 expression and disruption of ETV6 promotes TWIST1-dependent malignant phenotypes. Importantly, ETV6 is required to the anti-proliferation effects of EGFR-TKIs, partly due to the inhibitory function of ETV6 on TWIST1. We also found that EGFR-RAS signaling is tightly controlled by ETV6, supporting its role in TKI sensitivity.
Our study demonstrates that disruption of ETV6 contributes to EGFR-TKI resistance, which is likely due to derepression of TWIST1 and activation of EGFR-RAS signaling. Our results implicate ETV6 as a potential marker for predicting efficacy of an EGFR-targeted anticancer approach. Combination treatment of TWIST1 inhibitors could sensitize the anti-proliferation effects of EGFR-TKIs.
ETS variant gene 6 (ETV6), which belongs to the E26 transformation-specific (ETS) transcription factor family, was suggested to play a role as a tumor suppressor due to observed deletions in last stage, advanced prostate cancer [1,2,3,4]; however, its biological functions remain unclear. We demonstrated that ETV6 exhibits antitumor effects suppressing proliferation and metastatic progression, and found that epidermal growth factor receptor (EGFR) signaling inhibits ETV6 through a microRNA-mediated mechanism . EGFR signaling is involved in prostate cancer progression [6,7,8,9]; however, single-agent therapy using an EGFR tyrosine kinase inhibitor (TKI) was ineffective in castration-resistant prostate cancer (CRPC) . Since resistance to EGFR antagonists eventually develops and remains a challenging phenomenon , it is possible that loss of ETV6 function promotes the development of drug resistance in prostate cancer.
EGFR signaling follows three general steps: a ligand-monomeric EGFR interaction, dimerization (either homodimer or heterodimer) coupled with autophosphorylation through tyrosine kinase activity, and signal transduction for tumor-related properties . Therefore, mutations altering either the EGFR structure or activities of downstream components (e.g., KRAS) are recognized as key mechanisms responsible for resistance. However, several studies showed that the epithelial-to-mesenchymal transition (EMT) also plays a critical role in drug resistance consistent with the cancer-stem-cell (CSC) hypothesis [12, 13]. Indeed, the EMT-based gene profile is a powerful predictor of resistance to EGFR inhibitors . Also, one important transcription factor of the EMT, TWIST1 , was found to be associated with EGFR activation [16,17,18,19] and attributable to drug resistance [20,21,22]. In prostate cancer, Twist family BHLH transcription factor 1 (TWIST1) was found to be highly expressed in 90% of prostate cancers compared to 6.7% of benign hyperplasia  and is involved in the development of CRPC . Consistent with the CSC hypothesis, TWIST1 promotes tumor sphere formation, a functional indication of the self-renewal ability and CSC population .
Both PTEN and TP53 are frequently disrupted in prostate cancer; in addition, TP53 mutations occur in half of all CRPC [24, 25]. Following our earlier studies of ETV6 , we continued to investigate the molecular mechanism underlying its antitumor effects by utilizing prostate cancer cells derived from a prostate-specific Pten/Trp53 double-knockout mouse [24, 26]. We demonstrated that Etv6 associates at the promoter region of Twist1 and suppresses its transcription in a sequence-dependent manner. In human prostate cancer cells, ETV6 also inhibits TWIST1 expression and ETV6-knockdown can promote TWIST1-dependent malignant phenotypes. Importantly, perturbation of ETV6-TWIST1 axis can contribute to development of drug resistance. Prostate cancer cells with ETV6-knockdown are insensitive to TKIs while exogenous expression of ETV6 restores the anti-proliferative effects in the TKI-resistant RasB1 cell line, which expresses a mutated RAS oncogene [27, 28]. We also found an inhibitory circuit between ETV6 and EGFR-RAS signaling; therefore, there could be multiple mechanisms accounting for the drug-sensitizing effect of ETV6. Our results provide a molecular mechanism by which ETV6 suppresses tumor progression through transcriptional regulation of TWIST1 and disruption of EGFR-RAS signaling.
Cells, constructs, and reagents
The mouse AC1, AC3, C1, and C2 cell lines were isolated from PbCre4+;Ptenfl/fl;TP53fl/fl Luc + mouse prostate tumors and were established as previously described [24, 26]. AC1 and AC3 cells were cultured in PrEGM medium (Lonza, Walkersville, MD, USA); C1 cells were cultured in PrEGM/DHT with 5% serum and 5% 3 T3-conditioned medium; C2 cells were cultured in PrEGM/DHT with 5% 3 T3-conditioned medium. The mouse wild-type (WT) prostatic basal cell line was provided by Dr. Lei Fang (NCI/NIH, Bethesda, MD, USA) and was cultured in WIT-P medium (Stemgent, San Diego, CA, USA) as previously described. DU145, PC3, LNCaP, and 22RV1 human prostate cancer cell lines were obtained from ATCC (Rockville, MD, USA). The metastatic RasB1 cell line was previously characterized and used to study molecular mechanisms of prostate cancer metastasis in multiple peer-reviewed articles [27,28,29,30,31,32,33]. All human prostate cancer cell lines were cultured in RPMI 1640 medium supplemented with 10% fetal bovine serum (FBS). RasB1 and PC3 cells with stable expression of ETV6 were established by transfection with an ETV6 complementary (c)DNA-encoding or empty pCDH-CMV-MCS-EF1-Puro vector (System Biosciences, Palo Alto, CA, USA); 2 × 105 cells were seeded and transfected with 5 μg DNA and selected with puromycin for 1 month. Mouse and human ON-TARGETplus SMARTpool siRNAs (scrambled and ETV6) and a human shRNA vector (LacZ and ETV6) were from Dharmacon (Thermo Scientific, Waltham, MA, USA) and the RNAi Core Lab (Academia Sinica, Taipei, Taiwan), respectively. Transient transfections of plasmids and siRNAs were carried out using the X-tremeGENE HP DNA transfection reagent (Roche, CA, USA) or Lipofectamine RNAiMAX (Invitrogen, Carlsbad, CA, USA). Cells were treated with EGFR inhibitors, CI1033 (10 ng/ml) and AG1478 (10 μM) for 24 h in medium containing 10% serum. For EGF treatment, cells were subjected to serum-starvation for 24 h, followed by the addition of 100 ng/ml EGF for 24 h also in serum-free medium. The EGF was from R&D Systems (Minneapolis, MN, USA), and the EGFR inhibitors (CI1033 and AG1478) were from Selleck (Houston, TX, USA). The mouse Etv6-binding site was located upstream of mouse Twist1 on chromosome 12: 33957354 at GRCm38. The Twist1-red fluorescent protein (RFP) reporter containing the mouse Twist1 promoter with the Etv6 response element was constructed using a Clone-it Enzyme free Lentivectors Kit (System Biosciences). ETV6 response element mutations were made using a Site-Directed Mutagenesis System kit (Invitrogen). All primers used for these constructs are listed in Additional file 1; Table S1. All constructs were verified by a DNA sequence analysis.
Quantitative real-time reverse-transcription (qRT)-polymerase chain reaction (PCR)
An qRT-PCR was used to measure Etv6, Cdh1, Twist1, and Vim in mouse cell lines or ETV6 and TWIST1 expressions in human prostate cancer cell lines. Total RNA was isolated using the mirVana PARIS RNA isolation system (Thermo Scientific, Waltham, MA, USA). For RT, 3 μg of total RNA was used with the SuperScript III kit (Invitrogen). Samples containing primer pairs were mixed in SYBR green PCR master mix (Applied Biosystems, Waltham, MA, USA), and the amplification program was as follows: initial 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. All reactions were normalized to mouse Gapdh or human GAPDH expression and run in triplicate. All primers used for the PCR are listed in Additional file 1; Table S2.
Chromatin Immunoprecipitation (ChIP) assay
ChIP assays were performed using the EZ magna ChIP A kit (Millipore, Billerica, MD, USA) with a modified protocol. For small interfering (si)RNA treatment, 107 AC1 cells in 10-cm dishes were transfected with mouse scrambled or ETV6 siRNAs for 48 h. Cells were cross-linked with 1% formaldehyde in culture medium at room temperature for 15 min and then quenched with the addition of 1 ml of 10× glycine. Cells were washed twice with cold phosphate-buffered saline (PBS) containing a protease inhibitor (Roche) and centrifuged at 105 rpm. Cell pellets were resuspended in 0.5 ml of cell Lysis Buffer (BioRad, Hercules, CA, USA) and incubated on ice for 15 min. Nuclei were collected by centrifugation at 105 rpm and 4 °C for 10 min and resuspended in nuclear lysis buffer. Genomic DNA was sheared by a microtip during sonication (Branson Sonifier 250, Germany) following 15 cycles of a 20-s burst then 1 min of cooling on ice. This procedure resulted in DNA fragments sized approximately 100~ 300 bp. Sheared chromatin was aliquoted to perform immunoprecipitation with a control immunoglobulin G (IgG) antibody (Santa Cruz Biotechnology, Santa Cruz, CA, USA) or antibodies against ETV6 or Gapdh at 4 °C overnight. A qPCR was performed in triplicate with 2 μl of eluted chromatin. ChIP antibodies and PCR primers are listed in Additional file 1; Table S3.
Promoter reporter assay
For promoter reporter assays, AC1 or AC3 cells in 12-well plates (5 × 104 cells/well) were transiently transfected with 1 μg of the mouse Twist1-RFP reporter containing the Etv6 response element (RE). Cells were also pretreated with 100 nM siRNA (scrambled, mouse Etv6) or 1 μg DNA (empty vector or Etv6-expressing vector, OriGene, Rockville, MD, USA) by transfection. The promoter function was analyzed using fluorescence-activated cell sorting (FACS, BD Biosciences, San Jose, CA, USA), and relative median fluorescent intensity (MFI) values were measured as previously described . The MFI value for the RFP was measured by FACS using FACSDiva software (BD Biosciences) and was normalized to the value of the vehicle. Three independent experiments were run with triplicate samples.
Correlation analyses using human gene expression datasets
To compare ETV6 expression levels with prostate cancer progression and with TWIST1 expression levels, we used mRNA expression data from human prostate cancer databases of the Taylor dataset  and the Cancer Genome Atlas (TCGA). The study using the Taylor dataset was conducted under Memorial-Sloan Kettering Cancer Center (MSKCC) Institutional Review Board approval on 98 primary and 13 metastatic prostate cancer samples in addition to 28 normal prostate samples. Analysis of TCGA dataset was performed with 48 normal solid tissues.
Western blot analysis
Cells grown on 6-well plates (106 cells/well) were lysed in 150 μl RIPA buffer containing complete protease inhibitors (Roche) and phosphatase inhibitors (Roche), 25 mM β-glycerophosphate, 10 mM sodium fluoride, and 1 mM sodium vanadate. Twenty micrograms of protein was separated per lane by sodium dodecylsulfate (SDS)-gel electrophoresis. After being transferred to polyvinylidene difluoride membrane, blots were blocked with 5% BSA in PBST. Primary antibodies were incubated overnight at 4 °C, and secondary antibodies were incubated at room temperature for 1 h as indicated in Additional file 1; Table S4.
Twenty-two cases of prostatic adenocarcinoma were collected from the Taipei Medical University Joint human biological database (Taipei, Taiwan), and approved by the Taipei Medical University-Joint Institutional Review Board (approval no.: 201311034). RNA was extracted from dissected tissue containing greater than 70% tumor cell content. The method for separating the specimens into two groups of ‘low’ (TWIST1_L) and ‘high’ TWIST1 (TWIST1_H) expressions was pre-decided by half of the number of patients according to TWIST1 levels by an RT-qPCR.
Migration and invasion assay
For the migration and invasion assay, metastatic RasB1 and PC3 cells were stably transfected with the ETV6 expression vector or an empty vector. DU145 cells were stably transfected with a LacZ or ETV6 shRNA vector. LNCaP and 22RV1 cells were transiently transfected with SMARTpool ETV6 siRNA or control scrambled siRNA. Cells were resuspended at a concentration of 2.5 × 105 cells/ml in serum-free medium. Matrigel™ for the invasion assay was purchased from BD Biosciences (San Jose, CA, USA). Matrigel-coated transwell dishes were prepared by adding 200 μl of 10-fold serum-free medium-diluted Matrigel. In total, 2.5 × 105 cells/well in serum-free medium was plated above the Matrigel. The lower chamber was filled with 600 μl of serum-containing medium or serum-free medium with the addition of 200 ng/μL EGF. Cells that had invaded the Matrigel-coated transwells in response to EGF after 12 h were fixed and stained with a 0.5% crystal violet fixative solution for 15 min. Invaded cells on the underside of the membrane were counted and quantified by an enzyme-linked immunosorbent assay (ELISA) reader at OD 550 nm for each replicate in triplicate. The migration assay used transwells without Matrigel, and cells were fixed and stained as described for the invasion assay.
Animal work was performed in accordance with a protocol approved by the Taipei Medical University Animal Care and Use Committee (Taipei, Taiwan). For the metastasis analysis, DU145 cells harboring a LacZ or ETV6 shRNA vector were subjected to intracardiac injections into 5-week-old male nude mice (National Laboratory Animal Center, Taipei, Taiwan; six mice/group) at 105 cells per mouse. For survival studies, mice were euthanized when one of the following situations applied: 10% loss of body weight, paralysis, or head tilting. Hematoxylin and eosin (H&E) staining of brain tissues of mice was performed on day 80 after the injection as previously described . To analyze tumorigenesis, 5-week-old male nude mice (National Laboratory Animal Center; five mice/group) were subcutaneously injected with 106 RasB1 cells harboring an empty or ETV6 expression vector in 50% Matrigel™ (BD Biosciences). Subcutaneous tumors were harvested and measured from mice after treatment with 20 mg/kg CI1033 or DMSO as the control for 1 month as previously described .
DU145 cells were stably transfected with a LacZ or ETV6 or TWIST1 short hairpin (sh)RNA vector, or RasB1 cells were stably transfected with an empty or ETV6 or TWIST1 expression vector, and seeded at a density of 2 × 103 cells/well in 96-well plates. Cells were treated with 0, 0.1, 0.5, 1, 2, 5, and 10 nM CI1033 or 0, 0.1, 0.5, 1, 2, 5, and 10 μM AG1478 for 24 h, and analyzed using a Cell Proliferation Assay Kit (Promega, Madison, WI, USA) according to the manufacturer’s protocol.
All data are presented as the mean ± standard error of the mean (SEM). Statistical calculations were performed with GraphPad Prism analytical tools. Differences between individual groups were determined by Student’s t-test or a one-way analysis of variance (ANOVA) followed by Bonferroni’s post-test for comparisons among three or more groups. The association between ETV6 and TWIST1 expressions was compared using a Chi-squared test. p values of < 0.05 were considered statistically significant.
Etv6 recognizes the promoter of Twist1 and suppresses its transcription
When investigating the role of ETV6 in tumor progression, we hypothesized that ETV6 was inhibitory to the EMT. Our earlier study concluded that loss of Pten and TP53, two common genetic lesions in prostate cancer, promotes the EMT and cell plasticity in a mouse prostate-specific Pten/Trp53 double-knockout mouse model [24, 26]. When comparing the expression between one clonally derived cell line from this model and that of normal mouse prostate tissues, many EMT transcription factors were increased, while Etv6 was comparatively reduced (Additional file 1; Figure S1A). We further investigated a panel of mouse prostate cancer cell lines derived from the mouse model (Fig. 1A), and specifically analyzed the expression pattern of Twist1 due to its association with prostate cancer [22, 23]. Consistent with our hypothesis, Etv6 was positively correlated with an epithelial marker (Cdh1) (Pearson correlation, r = 0.966) but negatively with Twist1 (Pearson correlation, r = − 0.9002) (Fig. 1B). We further focused on a comparison between AC3 and AC1, since AC3, but not AC1, exhibits the transforming growth factor (TGF)-β-induced EMT . Again, both Etv6 and Cdh1 were reduced, while the mesenchymal markers (Twist1 and Vim) were increased in AC3 (Fig. 1C). In summary, our results showed that Etv6 is negatively associated with the EMT.
Because Etv6 is a transcription repressor , we searched for a consensus sequence on the promoter of Twist1 and identified one Etv6 response element (RE) at − 1286 (Fig. 1D). We performed chromatin immunoprecipitation (ChIP) assays and found an enriched signal at the Etv6 RE only when using an antibody against Etv6, but not the antibody against Gapdh, suggesting that the signal was specific to Etv6 (Fig. 1E). When using the same Etv6-specific antibody, we observed no enrichment at a site containing no Etv6-consensus sequence (Non-Etv6 RE), suggesting that the binding was sequence-specific (Fig. 1E). In addition, the ChIP signal derived from the Etv6 RE was significantly reduced after Etv6-knockdown using Etv6-specific siRNA (siEtv6), further supporting that binding was Etv6-dependent (Fig. 1F). We thus analyzed the cis-effect of the Etv6 RE on transcription by performing a reporter assay. The reporter activity of the construct containing the WT Etv6 RE was increased after Etv6 knockdown (scr. vs. siEtv6, Fig. 1G), while it was reduced after expression of exogenous Etv6 (EV vs. Etv6, Fig. 1H). The reporter construct containing a mutated (Mut) Etv6 RE was not responsive to different ETV6 statuses (Fig. 1G, H). These results suggest that Etv6 inhibits the EMT, partly through suppression of Twist1 by a physical interaction at the promoter of Twist1.
ETV6 is inversely correlated with TWIST1 in human prostate cancer
To further confirm the negative relationship between ETV6 and TWIST1 in human prostate cancer, we analyzed two public prostate cancer datasets. Many EMT drivers either showed positive or no significant correlation, instead of a negative one (Fig. 2A); however, only TWIST1 showed a negative correlation with ETV6 in both datasets, supporting ETV6-specific regulation (Fig. 2B, C). In addition, in the Taylor prostate dataset , mean expression of the ETV6 gene was significantly lower in primary tumors and further reduced in metastatic tumors compared those in normal tissues (Fig. 2D). On the contrary, the distribution profile was completely reversed with TWIST1 (Fig. 2E). Consistent with the profile analyzed from the Taylor dataset, both the messenger (m)RNA and protein levels of ETV6 were lower in the metastatic RasB1 cell line than those in non-metastatic cells (22RV1, LNCaP, and DU145) (Fig. 2F, G). Again, TWIST1 was negatively associated with ETV6 and increased only in metastatic RasB1 cells (Fig. 2F, G). We further investigated the relationship between ETV6 and TWIST1 in human prostate cancer tissues collected from the Taipei Medical University Joint human biological database (approval no.: 201,311,034, Taipei, Taiwan). We divided samples into two groups based on relative TWIST1 levels (TWIST1_H vs. TWIST1_L) and found a negative association between ETV6 and TWIST1 (Fig. 2H, left panel). An inverse correlation was also demonstrated in the same set of tissue samples (Fig. 2H, right panel). In summary, we concluded that TWIST1 is negatively associated with ETV6 and is involved in tumor progression in human prostate cancer.
ETV6 tightly controls TWIST1 expression and serves as a mediator of EGFR-TWIST1 signaling in human prostate cancer cells
ETV6 exhibits antitumor effects and can be negatively regulated by EGFR signaling . Since TWIST1 was shown to be a downstream effector following EGFR activation [16,17,18,19], ETV6 could be a crucial component in the pathway. To test this possibility, we confirmed that EGFR activity negatively regulated ETV6 mRNA (Fig. 3A), and that using either an EGFR inhibitor (CI1033) or overexpression of ETV6 could reduce TWIST1 mRNA in metastatic RasB1 cells (Fig. 3B). We observed the same effects by monitoring their protein levels (Fig. 3C). In non-metastatic prostate cancer cells, which express more ETV6 compared to RasB1 (Fig. 2G), ETV6 knockdown efficiently increased TWIST1 at both the mRNA and protein levels (Fig. 3D, E), suggesting that TWIST1 is tightly controlled by ETV6. To determine whether ETV6 is involved in EGFR-TWIST1 signaling, we found that treatments modulating EGFR activities (i.e., EGF and CI1033) no longer affected TWIST1 following ETV6 knockdown (siETV6, Fig. 3F). Our results support a novel EGFR-ETV6-TWIST1 pathway in that ETV6 serves as a gatekeeper to maintain TWIST1 at low levels in prostate cancer.
ETV6-knockdown leads to TWIST1-dependent malignant progression
We demonstrated that ETV6 efficiently suppressed metastasis of prostate cancer ; however, the underlying mechanism remained unclear. Based on current findings, decreasing TWIST1 expression could account for the antitumor effects of ETV6. To test this idea, we first confirmed the metastasis-related functions of TWIST1 by monitoring a metastatic PC3 cell line stably expressing exogenous ETV6, followed by transient expression of TWIST1 (Fig. 4A). Although ETV6 did suppress the malignant abilities (both migration and invasion), overexpression of TWIST1 clearly enhanced those in the presence of ETV6 (ETV6 vs. ETV6 + TWIST1, Fig. 4B, C). We observed the same malignant effects of TWIST1 when performing these experiments in RasB1 cells (Additional file 1; Figure S1B).
Consistent with the antitumor function of ETV6, stable ETV6-knockdown in DU145 cells efficiently promoted the malignant abilities (shLacZ vs. shETV6, Additional file 1; Figure S1C). Again, we confirmed that TWIST1 was induced after ETV6-knockdown (Fig. 4D). When simulating the metastasis process by delivering this pair of cells into mice by an intracardiac injection, ETV6-knockdown indeed reduced the life expectancy (Fig. 4E) and promoted metastasis according to tumor masses detected in the brain (arrow, Fig. 4F). Since DU145 was derived from a human prostate adenocarcinoma metastasizing to the brain , ETV6 knockdown could reactivate the metastatic properties of DU145 to the brain. To test whether the malignant effects following ETV6-knockdown were due to increased TWIST1, we successfully reduced both ETV6 and TWIST1 by specific siRNAs in a non-metastatic LNCaP cell line (Fig. 4G). Compared to ETV6-knockdown alone, which increased migration and invasion, additional TWIST1-knockdown reduced both functions to background levels (siETV6 vs. siETV6 + siTWIST1, Fig. 4H, I). We observed the same results when performing experiments using 22RV1 cells (Additional file 1; Figure S1D-F). In summary, we concluded that ETV6-knockdown leads to derepression of TWIST1 which contributes to tumor progression.
ETV6-TWIST1 signaling is involved in the development of resistance to EGFR antagonists
The EGFR is a validated target for cancer therapy; however, resistance to EGFR inhibitors eventually evolves . Since EGFR-based therapeutics showed no beneficial effects in prostate cancer [10, 36], it is important to determine whether the ETV6-TWIST1 axis plays a role in the development of drug resistance. We showed the anti-proliferative effects using one EGFR kinase inhibitor (AG1478) in DU145 cells; however, following ETV6-knockdown, cells became resistant (Fig. 5A). Interestingly, simultaneous knockdown of both ETV6 and TWIST1 recovered the anti-proliferative effect of the inhibitor (shETV6 + siTWIST1, Fig. 5A), consistent with our hypothesis that inducing TWIST1 after disruption of ETV6 contributes to malignant progression. The same experiment was performed using another EGFR inhibitor (CI1033, Additional file 1; Figure S1G).
The RasB1 cell line is derived from DU145 by introducing a mutant RAS and thus is resistant to EGFR TKIs [27, 28]. We demonstrated again that an EGFR antagonist (CI1033) did not suppress proliferation of RasB1, but overexpression of ETV6 clearly recovered the anti-proliferative function (EV vs. ETV6, Fig. 5B). However, RasB1 was resistant to CI1033 again after co-expression of both ETV6 and TWIST1, supporting the role of TWIST1 in the development of drug resistance (Fig. 5B). The same conclusion was derived using another TKI (AG1478, Additional file 1; Figure S1H). We also used a mouse xenograft model with the subcutaneous injection of RasB1 cells and tested the effect of ETV6 on CI1033. After the mice developed tumors, the mice were treated with CI1033 for 4 weeks. We found that CI1033 did not significantly reduce tumors (EV/DMSO vs. EV/CI1033); however, cells expressing ETV6 showed drastic antitumor effects (EV/DMSO vs. ETV6/DMSO, Fig. 5C), consistent with our previous findings . Further CI1033 treatment in mice injected with RasB1 cells expressing exogenous ETV6 completely suppressed tumor growth in three of five mice tested (ETV6/CI1033, Fig. 5C, D).
A mutual inhibition circuit exists between EGFR-RAS signaling and ETV6
In addition to the inhibitory role of ETV6 on TWIST1, we investigated the effect of ETV6 on EGFR-RAS signaling. We found that exogenous ETV6 negatively regulated the phosphorylation status of the EGFR (p-EGFR, Fig. 5E, left) in metastatic RasB1 cells. Interestingly, exogenous ETV6 also efficiently suppressed the phosphorylation status of extracellular signal-regulated kinase 1 and 2 (ERK1/2), a downstream signal transducer of EGFR-RAS signaling , while ETV6-knockdown in non-metastatic 22RV1 and DU145 cells increased the phosphorylation signal (p-ERK1/2, Fig. 5E, right). Thus, we hypothesized a mutual inhibition between EGFR-RAS signaling and ETV6. The inhibitory effect of ETV6 on EGFR-RAS signaling could also explain an earlier study showing that ETV6 can suppress RAS-induced transformation in an NIH3T3 cell model . Based on our aggregated results, we propose a working model that disruption of ETV6 contributes to tumor progression and TKI-resistance through derepression of TWIST1 and activation of EGFR-RAS signaling (Fig. 5F).
An earlier study showed that the genomic locus containing ETV6 is the most common translocation site in leukemia . While one allele can fuse to over 30 different genomic loci, the other usually undergoes deletion, which accounts for the frequently observed loss of heterozygosity . Therefore, ETV6 was considered a tumor suppressor, even in certain types of solid tumors . In prostate cancer, many genes from the ETS family participate in fusion transcripts with transmembrane serine protease isoform 2 (TMPRSS2); in fact, three ETS members (ERG, ETV1, and ETV4) contribute to about 80% of TMPRSS2 fusion [41, 42]. However, ETV6 has not yet been reported to be involved in chromosome translocation or fusion with TMPRSS2; instead, frequent deletions were observed in late stages, metastatic prostate cancer [1,2,3,4]. It was shown that following androgen signaling, both the androgen receptor (AR) and topoisomerase II beta (TOP2B) were localized to TMPRSS2-ERG genomic breakpoints, followed by TOP2B-mediated DNA breakage and recombination . Therefore, genomic rearrangements in prostate cancer, although frequently occurring, are likely dependent on the nature of the AR and TOP2B and are restricted to certain ETS loci except for ETV6. Since genomic lesions of ETV6 are common in leukemia, which does not rely on androgen, deletion of ETV6 might not be dependent on AR signaling.
Androgen deprivation therapy (ADT) is a standard procedure in prostate cancer; however, patients eventually develop metastatic castration-resistant prostate cancer. Therefore, patients could benefit from combining therapeutic approaches with different mechanisms. One would consider EGFR-targeted therapy since that the majority of prostate cancers are derived from an epithelial origin and are associated with elevated activities of the EGFR family [6,7,8,9]; however, using either an EGFR antagonist or anti-HER2 antibody did not achieve therapeutic effectiveness in CRPC [10, 36]. This might reflect the challenging issue for EGFR-targeted therapy in general since resistance inevitably occurs even though the antagonists are used in many types of cancer . Our proposed model that disruption of ETV6 leads to TKI resistance via derepression of both TWIST1 and EGFR-RAS signaling (Fig. 5F), provides one explanation to this issue. Since ETV6 is frequently deleted in late stage, malignant prostate cancer [1,2,3,4], the majority of CRPC could be lack of ETV6 function, consequently, irresponsive to TKIs. Prostate cancer still at earlier or hormone-sensitive stages containing intact ETV6 activity could be responsive to EGFR-targeted antagonists as monotherapy. Patients might benefit from prescreening their ETV6 statuses (genetics or expression levels) before treated with EGFR-TKIs.
Our results support a model that the EGFR facilitates tumor malignancy by reducing ETV6, which enhances TWIST1 activities. In addition, we previously reported that activation of EGFR signaling facilitates bone metastasis of prostate cancer through EGFR-mediated transcriptional suppression of microRNA-1 (miR-1) . Reduced miR-1 can increase TWIST1 function since miR-1 targets the 3′ untranslated region of TWIST1 and destabilizes its mRNA . Furthermore, an earlier study also showed that EGFR activation can trigger signal transducer and activator of transcription 3 (STAT3)-dependent transcription of TWIST1 and EMT . Based on conclusions from our studies and others, TWIST1 appears to be a key mediator that promotes malignant phenotypes; therefore, EGFR signaling utilizes multiple mechanisms to increase TWIST1 expression, including suppression of its negative regulators (miR-1 and ETV6) and activation of its positive regulator (e.g. STAT3) following EGFR activation. Therefore, designing novel TWIST1 inhibitors could sensitize the anti-proliferation effects and reduce the risk of resistance when choosing EGFR-TKIs.
Our studies provide a novel and testable hypothesis that connects ETV6-TWIST1 signaling to EGFR-TKI resistance during prostate cancer progression. Information about the genetic or expressional statuses of ETV6 could be valuable for decision making in personalized medicine when considering EGFR-based therapeutics. EGFR antagonists could achieve better and more-sustainable antitumor responses in combination with TWIST1-targeted agents.
androgen deprivation therapy
castration-resistant prostate cancer
epidermal growth factor receptor
ETS variant gene 6
signal transducer and activator of transcription 3
transforming growth factor β
tyrosine kinase inhibitor
Twist family BHLH transcription factor 1
Kibel AS, Schutte M, Kern SE, Isaacs WB, Bova GS. Identification of 12p as a region of frequent deletion in advanced prostate cancer. Cancer Res. 1998;58:5652–5.
Kibel AS, Faith DA, Bova GS, Isaacs WB. Mutational analysis of ETV6 in prostate carcinoma. Prostate. 2002;52:305–10.
Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS, Arora VK, Kaushik P, Cerami E, Reva B, et al. Integrative genomic profiling of human prostate cancer. Cancer Cell. 2010;18:11–22.
Kibel AS, Freije D, Isaacs WB, Bova GS. Deletion mapping at 12p12-13 in metastatic prostate cancer. Genes Chromosomes Cancer. 1999;25:270–6.
Tsai YC, Chen WY, Siu MK, Tsai HY, Yin JJ, Huang J, Liu YN. Epidermal growth factor receptor signaling promotes metastatic prostate cancer through microRNA-96-mediated downregulation of the tumor suppressor ETV6. Cancer Lett. 2017;384:1–8.
Traish AM, Morgentaler A. Epidermal growth factor receptor expression escapes androgen regulation in prostate cancer: a potential molecular switch for tumour growth. Br J Cancer. 2009;101:1949–56.
Cai CQ, Peng Y, Buckley MT, Wei J, Chen F, Liebes L, Gerald WL, Pincus MR, Osman I, Lee P. Epidermal growth factor receptor activation in prostate cancer by three novel missense mutations. Oncogene. 2008;27:3201–10.
Schlomm T, Kirstein P, Iwers L, Daniel B, Steuber T, Walz J, Chun FH, Haese A, Kollermann J, Graefen M, et al. Clinical significance of epidermal growth factor receptor protein overexpression and gene copy number gains in prostate cancer. Clin Cancer Res. 2007;13:6579–84.
Di Lorenzo G, Tortora G, D'Armiento FP, De Rosa G, Staibano S, Autorino R, D'Armiento M, De Laurentiis M, De Placido S, Catalano G, et al. Expression of epidermal growth factor receptor correlates with disease relapse and progression to androgen-independence in human prostate cancer. Clin Cancer Res. 2002;8:3438–44.
Canil CM, Moore MJ, Winquist E, Baetz T, Pollak M, Chi KN, Berry S, Ernst DS, Douglas L, Brundage M, et al. Randomized phase II study of two doses of gefitinib in hormone-refractory prostate cancer: a trial of the National Cancer Institute of Canada-clinical trials group. J Clin Oncol. 2005;23:455–60.
Ciardiello F, Tortora G. EGFR antagonists in cancer treatment. N Engl J Med. 2008;358:1160–74.
Shibue T, Weinberg RA. EMT, CSCs, and drug resistance: the mechanistic link and clinical implications. Nat Rev Clin Oncol. 2017;14(10):611–29.
Singh A, Settleman J. EMT, cancer stem cells and drug resistance: an emerging axis of evil in the war on cancer. Oncogene. 2010;29:4741–51.
Byers LA, Diao L, Wang J, Saintigny P, Girard L, Peyton M, Shen L, Fan Y, Giri U, Tumula PK, et al. An epithelial-mesenchymal transition gene signature predicts resistance to EGFR and PI3K inhibitors and identifies Axl as a therapeutic target for overcoming EGFR inhibitor resistance. Clin Cancer Res. 2013;19:279–90.
Puisieux A, Brabletz T, Caramel J. Oncogenic roles of EMT-inducing transcription factors. Nat Cell Biol. 2014;16:488–94.
Chang YS, Chen WY, Yin JJ, Sheppard-Tillman H, Huang J, Liu YN. EGF receptor promotes prostate cancer bone metastasis by Downregulating miR-1 and activating TWIST1. Cancer Res. 2015;75:3077–86.
Lo HW, Hsu SC, Xia W, Cao X, Shih JY, Wei Y, Abbruzzese JL, Hortobagyi GN, Hung MC. Epidermal growth factor receptor cooperates with signal transducer and activator of transcription 3 to induce epithelial-mesenchymal transition in cancer cells via up-regulation of TWIST gene expression. Cancer Res. 2007;67:9066–76.
Pallier K, Cessot A, Cote JF, Just PA, Cazes A, Fabre E, Danel C, Riquet M, Devouassoux-Shisheboran M, Ansieau S, et al. TWIST1 a new determinant of epithelial to mesenchymal transition in EGFR mutated lung adenocarcinoma. PLoS One. 2012;7:e29954.
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.
Cosset E, Hamdan G, Jeanpierre S, Voeltzel T, Sagorny K, Hayette S, Mahon FX, Dumontet C, Puisieux A, Nicolini FE, Maguer-Satta V. Deregulation of TWIST-1 in the CD34+ compartment represents a novel prognostic factor in chronic myeloid leukemia. Blood. 2011;117:1673–6.
Ruan D, He J, Li CF, Lee HJ, Liu J, Lin HK, Chan CH. Skp2 deficiency restricts the progression and stem cell features of castration-resistant prostate cancer by destabilizing Twist. Oncogene. 2017;
Kwok WK, Ling MT, Lee TW, Lau TC, Zhou C, Zhang X, Chua CW, Chan KW, Chan FL, Glackin C, et al. Up-regulation of TWIST in prostate cancer and its implication as a therapeutic target. Cancer Res. 2005;65:5153–62.
Shiota M, Yokomizo A, Tada Y, Inokuchi J, Kashiwagi E, Masubuchi D, Eto M, Uchiumi T, Naito S. Castration resistance of prostate cancer cells caused by castration-induced oxidative stress through Twist1 and androgen receptor overexpression. Oncogene. 2010;29:237–50.
Martin P, Liu YN, Pierce R, Abou-Kheir W, Casey O, Seng V, Camacho D, Simpson RM, Kelly K. Prostate epithelial Pten/TP53 loss leads to transformation of multipotential progenitors and epithelial to mesenchymal transition. Am J Pathol. 2011;179:422–35.
Sircar K, Yoshimoto M, Monzon FA, Koumakpayi IH, Katz RL, Khanna A, Alvarez K, Chen G, Darnel AD, Aprikian AG, et al. PTEN genomic deletion is associated with p-Akt and AR signalling in poorer outcome, hormone refractory prostate cancer. J Pathol. 2009;218:505–13.
Liu YN, Abou-Kheir W, Yin JJ, Fang L, Hynes P, Casey O, Hu D, Wan Y, Seng V, Sheppard-Tillman H, et al. Critical and reciprocal regulation of KLF4 and SLUG in transforming growth factor beta-initiated prostate cancer epithelial-mesenchymal transition. Mol Cell Biol. 2012;32:941–53.
Yin J, Pollock C, Tracy K, Chock M, Martin P, Oberst M, Kelly K. Activation of the RalGEF/Ral pathway promotes prostate cancer metastasis to bone. Mol Cell Biol. 2007;27:7538–50.
Siu MK, Abou-Kheir W, Yin JJ, Chang YS, Barrett B, Suau F, Casey O, Chen WY, Fang L, Hynes P, et al. Loss of EGFR signaling regulated miR-203 promotes prostate cancer bone metastasis and tyrosine kinase inhibitors resistance. Oncotarget. 2014;5:3770–84.
Yin JJ, Zhang L, Munasinghe J, Linnoila RI, Kelly K. Cediranib/AZD2171 inhibits bone and brain metastasis in a preclinical model of advanced prostate cancer. Cancer Res. 2010;70:8662–73.
Yin J, Liu YN, Tillman H, Barrett B, Hewitt S, Ylaya K, Fang L, Lake R, Corey E, Morrissey C, et al. AR-regulated TWEAK-FN14 pathway promotes prostate cancer bone metastasis. Cancer Res. 2014;74:4306–17.
Siu MK, Tsai YC, Chang YS, Yin JJ, Suau F, Chen WY, Liu YN. Transforming growth factor-beta promotes prostate bone metastasis through induction of microRNA-96 and activation of the mTOR pathway. Oncogene. 2015;34:4767–76.
Liu YN, Yin J, Barrett B, Sheppard-Tillman H, Li D, Casey OM, Fang L, Hynes PG, Ameri AH, Kelly K. Loss of androgen-regulated MicroRNA 1 activates SRC and promotes prostate cancer bone metastasis. Mol Cell Biol. 2015;35:1940–51.
Chen WY, Liu SY, Chang YS, Yin JJ, Yeh HL, Mouhieddine TH, Hadadeh O, Abou-Kheir W, Liu YN. MicroRNA-34a regulates WNT/TCF7 signaling and inhibits bone metastasis in Ras-activated prostate cancer. Oncotarget. 2015;6:441–57.
Lopez RG, Carron C, Oury C, Gardellin P, Bernard O, Ghysdael J. TEL is a sequence-specific transcriptional repressor. J Biol Chem. 1999;274:30132–8.
Stone KR, Mickey DD, Wunderli H, Mickey GH, Paulson DF. Isolation of a human prostate carcinoma cell line (DU 145). Int J Cancer. 1978;21:274–81.
Solit DB, Rosen N. Targeting HER2 in prostate cancer: where to next? J Clin Oncol. 2007;25:241–3.
Roberts PJ, Der CJ. Targeting the Raf-MEK-ERK mitogen-activated protein kinase cascade for the treatment of cancer. Oncogene. 2007;26:3291–310.
Rompaey LV, Potter M, Adams C, Grosveld G. Tel induces a G1 arrest and suppresses Ras-induced transformation. Oncogene. 2000;19:5244–50.
Shurtleff SA, Buijs A, Behm FG, Rubnitz JE, Raimondi SC, Hancock ML, Chan GC, Pui CH, Grosveld G, Downing JR. TEL/AML1 fusion resulting from a cryptic t(12;21) is the most common genetic lesion in pediatric ALL and defines a subgroup of patients with an excellent prognosis. Leukemia. 1995;9:1985–9.
Hatta Y, Takeuchi S, Yokota J, Koeffler HP. Ovarian cancer has frequent loss of heterozygosity at chromosome 12p12.3-13.1 (region of TEL and Kip1 loci) and chromosome 12q23-ter: evidence for two new tumour-suppressor genes. Br J Cancer. 1997;75:1256–62.
Gutierrez-Hartmann A, Duval DL, Bradford AP. ETS transcription factors in endocrine systems. Trends Endocrinol Metab. 2007;18:150–8.
Tomlins SA, Laxman B, Dhanasekaran SM, Helgeson BE, Cao X, Morris DS, Menon A, Jing X, Cao Q, Han B, et al. Distinct classes of chromosomal rearrangements create oncogenic ETS gene fusions in prostate cancer. Nature. 2007;448:595–9.
Haffner MC, Aryee MJ, Toubaji A, Esopi DM, Albadine R, Gurel B, Isaacs WB, Bova GS, Liu W, Xu J, et al. Androgen-induced TOP2B-mediated double-strand breaks and prostate cancer gene rearrangements. Nat Genet. 2010;42:668–75.
We thank Ai-Ning Shao, Kuo-Ching Jiang, Wei-Hau Chen, and Hong-Yuan Tsai (Taipei Medical University) for technical assistance.
This work was supported by the Ministry of Science and Technology of Taiwan to YCT (MOST104–2320-B-038-055-MY3), WYC (MOST106–2320-B-038-057), and YNL (MOST104–2314-B-038-045-MY3 and MOST105–2628-B-038 -006 -MY3), by Taipei Medical University-Wan Fang Hospital to WYC (105TMU-WFH-04), by the National Health Research Institutes of Taiwan to YNL (NHRI-EX107-10702BI), and by the Health and Welfare Surcharge of Tobacco Products to YNL (MOHW106-TDU-B-212-144001).
Availability of data and materials
The human clinical datasets utilized in the study are included in the Methods sections.
Ethics approval and consent to participate
We followed the guideline of the Taipei Medical University-Joint Institutional Review Board (approval no.: 201,311,034).
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.
About this article
Cite this article
Tsai, YC., Zeng, T., Abou-Kheir, W. et al. Disruption of ETV6 leads to TWIST1-dependent progression and resistance to epidermal growth factor receptor tyrosine kinase inhibitors in prostate cancer. Mol Cancer 17, 42 (2018). https://doi.org/10.1186/s12943-018-0785-1