TAZ Expression as a Prognostic Indicator in Colorectal Cancer

Abstract

The Hippo pathway restricts the activity of transcriptional coactivators TAZ (WWTR1) and YAP. TAZ and YAP are reported to be overexpressed in various cancers, however, their prognostic significance in colorectal cancers remains unstudied. The expression levels of TAZ and YAP, and their downstream transcriptional targets, AXL and CTGF, were extracted from two independent colon cancer patient datasets available in the Gene Expression Omnibus database, totaling 522 patients. We found that mRNA expressions of both TAZ and YAP were positively correlated with those of AXL and CTGF (p<0.05). High level mRNA expression of TAZ, AXL or CTGF significantly correlated with shorter survival. Importantly, patients co-overexpressing all 3 genes had a significantly shorter survival time, and combinatorial expression of these 3 genes was an independent predictor for survival. The downstream target genes for TAZ-AXL-CTGF overexpression were identified by Java application MyStats. Interestingly, genes that are associated with colon cancer progression (ANTXR1, EFEMP2, SULF1, TAGLN, VCAN, ZEB1 and ZEB2) were upregulated in patients co-overexpressing TAZ-AXL-CTGF. This TAZ-AXL-CTGF gene expression signature (GES) was then applied to Connectivity Map to identify small molecules that could potentially be utilized to reverse this GES. Of the top 20 small molecules identified by connectivity map, amiloride (a potassium sparing diuretic,) and tretinoin (all-trans retinoic acid) have shown therapeutic promise in inhibition of colon cancer cell growth. Using MyStats, we found that low level expression of either ANO1 or SQLE were associated with a better prognosis in patients who co-overexpressed TAZ-AXL-CTGF, and that ANO1 was an independent predictor of survival together with TAZ-AXL-CTGF. Finally, we confirmed that TAZ regulates Axl, and plays an important role in clonogenicity and non-adherent growth in vitro and tumor formation in vivo. These data suggest that TAZ could be a therapeutic target for the treatment of colon cancer.

About the Authors:

Hiu-Fung Yuen

Contributed equally to this work with: Hiu-Fung Yuen, Cian M. McCrudden, Yu-Han Huang

Affiliation: Institute of Molecular and Cell Biology, Agency for Science, Technology and Research (A*STAR), Singapore, Singapore

Cian M. McCrudden

Contributed equally to this work with: Hiu-Fung Yuen, Cian M. McCrudden, Yu-Han Huang

Affiliation: Centre for Cancer Research and Cell Biology (CCRCB), Queen's University Belfast, Belfast, United Kingdom

Yu-Han Huang

Contributed equally to this work with: Hiu-Fung Yuen, Cian M. McCrudden, Yu-Han Huang

Affiliation: Institute of Molecular and Cell Biology, Agency for Science, Technology and Research (A*STAR), Singapore, Singapore

Jill M. Tham

Affiliation: Institute of Molecular and Cell Biology, Agency for Science, Technology and Research (A*STAR), Singapore, Singapore

Xiaoqian Zhang

Affiliation: Institute of Molecular and Cell Biology, Agency for Science, Technology and Research (A*STAR), Singapore, Singapore

Qi Zeng

Affiliations Institute of Molecular and Cell Biology, Agency for Science, Technology and Research (A*STAR), Singapore, Singapore, Department of Biochemistry, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore

Shu-Dong Zhang

* E-mail: [email protected] (SZ); [email protected] (WH)

Affiliation: Centre for Cancer Research and Cell Biology (CCRCB), Queen's University Belfast, Belfast, United Kingdom

WanJin Hong

* E-mail: [email protected] (SZ); [email protected] (WH)

Affiliations Institute of Molecular and Cell Biology, Agency for Science, Technology and Research (A*STAR), Singapore, Singapore, School of Pharmaceutical Sciences, Xiamen University, Xiamen, Fujian, People's Republic of China

Introduction

The Hippo pathway plays an important role in cell proliferation, organ size control, and cancer development and progression [1], [2], [3], [4]. YAP and TAZ are both transcriptional co-activators that are inhibited by the Hippo pathway [1], [2], [3], [4]. Aberrant inactivation of the Hippo pathway and/or overexpression of TAZ and YAP results in transcriptional activation of their downstream targets [1], [2], [3], [4].

YAP overexpression induces cell proliferation and epithelial mesenchymal transition (EMT), and inhibits apoptosis and contact inhibition [5], [6], [7]. Transcriptional activation of epidermal growth factor receptor ligand amphiregulin may account for YAP-mediated induction of cell proliferation, especially under serum-depletion [7], while YAP also cooperates with Myc to promote cell proliferation [8]. Recently, YAP has been shown to play a critical role in stem cell biology. It is induced during pluripotent stem cell reprogramming, whilst silencing of YAP reduces the pluripotency of embryonic stem cells [9]. YAP promotes ovarian cancer progression, and high levels of nuclear expression are inversely associated with patient survival [10]. In particular, YAP is associated with clear cell ovarian tumors, an ovarian malignancy subtype with poor prognosis [11]. YAP has also been shown to play an oncogenic role in esophageal squamous cell carcinoma [12]. In liver cancer, microRNA-mediated inhibition of YAP inhibits tumor characteristics including cell proliferation and invasion [13]. Conversely, there are reports showing an opposite, tumor suppressive, role of YAP in promoting p73-mediated apoptosis [14], [15]. In breast and head and neck cancers, YAP has been shown to act as a tumor suppressor in certain circumstances [14], [16].

TAZ is structurally homologous to YAP, is likewise inhibited by the Hippo pathway, and also promotes EMT-mediated cancer progression [17], [18], [19]. TAZ regulates mesenchymal stem cell differentiation by modulating Runx2- and PPARgamma-dependent gene expression [20], as well as stem cell self-renewal through controlling localization of Smad [21]. TAZ plays an important role in the progression of breast [22], [23] and non-small cell lung cancer [24], [25]. Importantly, TAZ confers cancer stem cell-related traits on breast cancer cells, further highlighting its importance in tumor initiation and progression [26]. TAZ is also overexpressed in papillary thyroid carcinoma [27].

TAZ and YAP have been shown to interact with several transcriptional factors [2], [3], [4], with the TEAD family of transcriptional factors (TEAD1-4) being the most relevant in cell proliferation and cancer progression [19], [28], [29]. The X-ray crystal structures of YAP-TEAD complexes have been resolved and the proposed interaction is supported by and consistent with functional analysis [30], [31], [32], showing that YAP-TEAD complexes activates gene transcription.

YAP expression was observed in colon adenocarcinoma [33], [34], [35]. It is overexpressed in human colon cancer specimens and overexpression of YAP promotes cell proliferation and survival in colon cancer cells [36]. Recent findings show that knockdown of TAZ results in a decrease in cell proliferation in culture and tumor growth in vivo [26].

Despite evidence suggesting the potential implication of YAP and TAZ in colon cancer progression, their prognostic significance in colorectal cancer is unknown. In this study, we analyzed the mRNA expression of YAP and TAZ, and two of its downstream target genes, AXL and CTGF, in two independent colon cancer patient cohorts comprising 522 patients. We found that TAZ, but not YAP, is a prognostic marker in colon cancer progression. Furthermore, TAZ-AXL-CTGF co-overexpression, which defines both the expression of TAZ and its transcriptional activity on target gene expression, is a novel prognostic indicator, that is independent of tumor grade and stage, for colon cancer patients. The role of TAZ in colon cancer cell proliferation and oncogenesis was validated by functional study.

Materials and Methods

Extraction of clinical and microarray gene expression data from colon cancer patient datasets

Two colon cancer patient datasets, GSE14333 [37] and GSE17538 [38], available in the Gene Expression Omnibus (GEO) Database (http://www.ncbi.nlm.nih.gov/gds) were included in this study. The GEO website has standardized URLs for its individual datasets, e.g. the overall summary information about the microarray dataset GSE14333 can be accessed at http://www.ncbi.nlm.gov/geo/query/acc.cgi?acc=GSE14333. For each GEO data series, links are provided at the bottom of the page to the Series Matrix File(s), which contain the expression values for each gene (probeset) and each microarray. The URLs to the Series Matrix File(s) are also standardized. For GSE14333, the URL was ftp://ftp.ncbi.nlm.nih.gov/pub/geo/DATA/SeriesMatrix/GSE14333. The files in gzip format were then unzipped to the tab-delimited text format, which contain detailed information for statistical analysis. The GSE14333 and GSE17538 datasets are the two largest colon cancer patient datasets on the database, and comprise 522 patients, 458 of whose survival data is available in the database. The GSE17538 data series consists of four SubSeries: GSE17536, GSE17537, GSE19072 and GSE19073. GSE19072 and GSE19073 were excluded from this study as they lack clinical data. Microarray gene expression data were retrieved from the data matrices deposited to the GEO database by the original authors. The gene expression levels in GSE14333 and GSE17538 are represented by base-2 logarithm of the MAS5 value and the RMA values, respectively, as adopted by the original authors. R scripting was used to extract the expression values of a small number of genes (probesets) of interest and the clinical data from the data matrixes downloaded from GEO.

Demographic and clinical data of the two patient cohorts

Both age and gender of the patients were available demographic data in the two datasets analyzed in the present study. Patients had a median age of 67 years (Range: 26–92 year-old) and 65.5 years (Range: 23–94 year-old) in GSE14333 and GSE17538, respectively. There were 43% and 46% female, respectively, in cohort GSE14333 and GSE17538. The GSE14333 cohort consists of 290 patients, for whom 226 survival data were available. Fifteen, 32, 31 and 21% of patients had stage A, B, C and D tumors, respectively. The survival data were not available for all the patients of stage D. Survival data were available for all 232 patients in the GSE17538 cohort. Twelve, 31, 33 and 24% patients had tumors of AJCC stage I, II, III and IV, respectively, and 8, 78 and 14% of patients had Grade A, B and C tumors respectively.

Correlations of gene expression levels and clinical data

All statistical analyses were performed using SPSS19.0. The associations between expression levels of genes were analyzed by Spearman's rank test. Expression levels were further divided into high and low levels using median expression level as the cut-off point for Kaplan-Meier survival analysis. Results were compared by log-rank test. Univeriate Cox regression analysis was used to correlate the gene expression levels and patient survival and multivariate Cox regression analysis was used to identify independent predictors for patient survival using a backward stepwise approach with an entry limit of p<0.1 and a removal limit of p>0.05. Patients were divided further into 4 groups based on the expression levels of TAZ, AXL and CTGF; the TAZ-AXL-CTGF-low group consisted of patients who expressed all these 3 genes at low levels; the TAZ-AXL-CTGF-intermediate-low group consisted of patients who expressed one of these three genes at a high level; the TAZ-AXL-CTGF-intermediate-high group consisted of patients who expressed two of these three genes at high levels; the TAZ-AXL-CTGF-high group consisted of patients who expressed all three genes at high levels. The survival time of patients stratified by this grouping method were analyzed by Kaplan-Meier analysis and Cox regression as described above.

Identification of TAZ-AXL-CTGF co-expressing genes

Patients were stratified into four groups based on the expression levels of TAZ, AXL and CTGF as described above. The gene expression patterns of patients in TAZ-AXL-CTGF low subgroup and those in the TAZ-AXL-CTGF high subgroup (whose survival was significantly poorer) were compared. Probesets that were differentially expressed between these two subgroups were identified by 2-sample Welch's T-test. This test was used to avoid the type I error due to unequal variances of the values of probesets between subgroups. Briefly, a Welch's t test was applied to each probeset corresponding to a certain gene in the data matrix using our own Java application MyStats. P values and the differential expression in fold changes for all the probesets were generated as tab-delimited worksheets of Excel for further analysis. The genes were prioritized by ascending p-values. The top 100 probesets were prioritized in both patient datasets, and the genes common to both datasets were analyzed further.

Identification of potential inhibitory compounds targeting TAZ-AXL-CTGF overexpressing colon cancer (Connectivity Map)

Gene expression connectivity mapping was performed using Statistically Significant Connection's Map (sscMap) to identify candidate small molecule compounds that may inhibit the expression of genes that are co-regulated in TAZ-AXL-CTGF co-expressing aggressive colon cancer [39], [40], [41]. Of the probesets identified to be co-expressed with TAZ-AXL-CTGF, 33 were present on the Affymetrix HG-U133A microarray platform, which was used to generate the microarray database for the Connectivity Map [39]. The compiled gene signature was then fed to the Java application sscMap [41] as a query signature, and its association with the 6000 gene expression profiles generated by treating cancer cells with over 1000 small molecules were compared. The gene signature perturbation procedure, which increases the specificity of the output results, was applied as previously described [42]. All the small molecular compounds, that were negatively associated with the TAZ-AXL-CTGF-GES, were sorted and ranked by their p-value, perturbation stability and standardized connection score. The p-value that was considered significant was set at a stringent threshold (p = 1/1309), ensuring that the results generated by sccMap yield only maximally one expected false positive small molecule over the 1309 small molecules tested in the sccMap [42]. The top 20 small molecules were then entered into the Pubmed (www.pubmed.com) search engine together with colon cancer to identify research articles that have described their effects of the particular molecules on treatment of colon cancer.

Identification of therapeutic targets for colon cancer patients overexpressing TAZ-AXL-CTGF

Patients who co-overexpressed TAZ-AXL-CTGF were stratified into two groups based on their survival statuses. Differential expressions of different probesets between patients in the TAZ-AXL-CTGF-alive subgroup and those in the TAZ-AXL-CTGF-deceased subgroup were identified as described above.

Cell culture and retroviral transduction

HCT116 and SW620 cells were obtained from American Type Culture Collection and maintained in F12K/DME medium supplemented with 10% Fetal Bovine Serum (FBS), 10 ug/ml Penicillin/Streptomycin (P/S) (Life technologies, Carlsbad, CA). The amphotropic Phoenix packaging cell line was obtained from the Nolan Laboratory (Stanford University) and maintained in DMEM medium supplemented with 10% FBS and 10 ug/ml P/S. Retroviral infection was performed as previously described [18]. The short hairpin (shRNA) against human TAZ construct (5′-AGGTACTTCCTCAATCACA-3′) carried by the pSuperRetro-puro vector (shTAZ) was used for TAZ knockdown (Oligoengine, Seattle, WA) while the scramble shRNA construct (5′-CCTAAGGTTAAGTCGCCCTCG-3′; shScr) was used for control. Stable cell lines were established by selecting the traduced cells in 2 ug/ml puromycin (Sigma Aldrich, St. Louis, MO).

Western blot analysis

Immunoblotting was performed as previously described [18] using SuperSignal West Pico (Thermo Scientific, Rockford, IL). The commercial TAZ antibody was obtained from Imgenex (San Diego, CA).

Anchorage-independent soft agar and clonogenic assays

For the soft agar assay, 5000 cells resuspended in 0.35% (w/v) agarose in culture medium were overlaid on a solidified 0.5% (w/v) agarose in culture medium. The upper layer was allowed to solidify. Medium with puromycin was added the following day and cells were then incubated at 37°-C with 5% CO2. Fresh medium with puromycin was supplemented and the colonies formed were stained with 1 mg/ml Thiazolyl Blue Tetrazolium Bromide (Sigma Aldrich) for 4 hours in the incubator at 37°C with 5% CO2. The excess dye was removed by destaining multiple times with water and the number of colonies was determined.

For clonogenic assay, 500 cells were seeded per well in triplicate in 6-well plate. The culture medium was refreshed every week and the cells were fixed with 4% paraformaldehyde in PBS after 2 weeks incubation at 37°C with 5% CO2. The fixed cells were stained with 0.5% crystal violet (Sigma Aldrich) in 20% ethanol overnight. The cells were rinsed with water, dried and colony number was analyzed.

Tumorigenesis in nude mice

Hundred ul of a cell suspension of 1.5×107/ml were inoculated subcutaneously in the left and right hind flanks of four-to-six week-old female nude mice. Tumor development was monitored sfter 2 weeks. Mice were then euthanized and the tumors were removed for analysis.

Results

TAZ and YAP mRNA expressions positively correlate with mRNA expression of their downstream target genes, AXL and CTGF

Previously, we and others have shown that AXL and CTGF are two important downstream target genes of TAZ and YAP [23], [43], [44]. In the present study, we investigated whether the mRNA expression levels of the two transcriptional co-activators in the Hippo pathway, TAZ and YAP, correlate with the mRNA expression of AXL and CTGF. In the 290 colon cancer patients from the GSE14333 dataset, TAZ expression was significantly correlated with both AXL (Spearman's rank test, r = 0.547, p<0.001; Figure 1A) and CTGF (r = 0.543, p<0.001; Figure 1B) expressions. YAP mRNA expression was also positively correlated with AXL (r = 0.154, p = 0.009; Figure 1C) and CTGF (r = 0.141, p = 0.016; Figure 1D) mRNA expression in the same dataset, but to a lesser extent. In 232 colon cancer patients from GSE17538, TAZ mRNA expression was significantly positively correlated with both AXL (r = 0.752, p<0.001; Figure 1E) and CTGF (r = 0.686, p<0.001; Figure 1F) mRNA expressions, while YAP mRNA was also significantly positively correlated with mRNA expression of both genes, again to a lesser extent (AXL: r = 0.343, p<0.001; Figure 1G and CTGF: r = 0.387, p<0.001; Figure 1H).

[Figure omitted. See PDF.]

Figure 1. The correlations among mRNA expression of TAZ (WWTR1), YAP (YAP1), AXL and CTGF in colon cancer specimens.

Scatter plots for (A) TAZ mRNA expression versus AXL mRNA expression, (B) TAZ mRNA expression versus CTGF mRNA expression, (C) YAP mRNA expression versus AXL mRNA expression, and (D) YAP mRNA expression versus CTGF mRNA expression in the GSE14333 colon cancer datasets, and (E) TAZ mRNA expression versus AXL mRNA expression, (F) TAZ mRNA expression versus CTGF mRNA expression, (G) YAP mRNA expression versus AXL mRNA expression, and (H) YAP mRNA expression versus CTGF mRNA expression in the GSE17538 colon cancer datasets.

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

TAZ, but not YAP, mRNA expression is a predictor for patient survival

In the 226 patients whose survival data were available from the GSE14333 colon cancer patient cohort, a high level of TAZ mRNA expression was significantly correlated with a shorter survival (high level: mean survival = 72.3 months, 95% Confidence Interval (CI) = 63–81 months; low level: mean survival = 129 months, 95% CI = 121–136 months, p<0.001; Figure 2A). By Cox-regression analysis, TAZ mRNA expression was significantly correlated with survival (Hazards Ratio (HR) = 2.251, 95% CI = 1.626–3.116, p<0.001; Figure 2C) in the GSE14333 colon cancer patient cohort. By multivariate analysis (Figure 2C), TAZ mRNA expression (HR = 2.062, 95% CI = 1.472–3.116, p<0.001) and tumor staging (p<0.001) are both independent predictors of survival in the same cohort. Similarly, a high level of TAZ mRNA expression was significantly correlated with shorter patient survival in the GSE17538 colon cancer patient cohort (high level: mean survival = 84 months, 95% CI = 72–96 months; low level: mean survival = 109 months, 95% CI = 97–120 months, p = 0.011; Figure 2B). By Cox-regression, TAZ mRNA expression is a predictor of survival in this cohort (HR = 1.743, 95% CI = 1.177–2.582, p = 0.006; Figure 2D). Mulivariate Cox-regression analysis also showed that TAZ mRNA expression (HR = 1.998, 95% CI = 1.245–3.205, p = 0.004; Figure 2D) was an independent predictor for survival together with stage (p<0.001) and grade (p = 0.03) of the cancers in this colon cancer patient cohort. On the other hand, YAP mRNA expression did not significantly correlate with patient survival by Kaplan-Meier analysis (GSE14333: p = 0.519; Figure 2E and GSE17538: p = 0.634; Figure 2F) or by Cox-regression analysis (GSE14333: p = 0.673; Figure 2G and GSE17538: p = 0.979; Figure 3H) in either dataset. These results suggest that TAZ mRNA expression is a novel prognostic marker for colon cancer patients, but YAP is not.

[Figure omitted. See PDF.]

Figure 2. The associations between TAZ or YAP, and survival in colon cancer patients.

Kaplan-Meier analyses for TAZ mRNA expression in (A) GSE14333 and (B) GSE17538 colon cancer patient datasets. Univariate and Multivariate Cox regression analyses for TAZ mRNA expression, age, tumor stage and tumor grade in (C) GSE14333 and (D) GSE17538 colon cancer patient datasets. Kaplan-Meier analyses for YAP mRNA expression in (E) GSE14333 and (F) GSE17538 colon cancer patient datasets. Univariate and Multivariate Cox regression analyses for TAZ mRNA expression, age, tumor stage and tumor grade in (G) GSE14333 and (H) GSE17538 colon cancer patient datasets.

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

[Figure omitted. See PDF.]

Figure 3. The associations between AXL or CTGF, and survival in colon cancer patients.

Kaplan-Meier analyses for AXL mRNA expression in (A) GSE14333 and (B) GSE17538 colon cancer patient datasets. Univariate and Multivariate Cox regression analyses for AXL mRNA expression, age, tumor stage and tumor grade in (C) GSE14333 and (D) GSE17538 colon cancer patient datasets. Kaplan-Meier analyses for CTGF mRNA expression in (E) GSE14333 and (F) GSE17538 colon cancer patient datasets. Univariate and Multivariate Cox regression analyses for CTGF mRNA expression, age, tumor stage and tumor grade in (G) GSE14333 and (H) GSE17538 colon cancer patient datasets.

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

Both AXL and CTGF, downstream target genes of TAZ and YAP, are predictors of patient survival

To investigate whether the functional outcome of the transcriptional program mediated by TAZ and YAP confers prognostic value in colon cancer patients, we analyzed whether the mRNA expression of AXL and CTGF, two of the most well defined target genes of TAZ/YAP-TEAD complexes, correlate with survival of patients in the two colon cancer patient cohorts.

High level mRNA expression of AXL correlated with a shorter colon cancer patient survival time in the GSE14333 cohort, although this correlation was not significant (high level: mean survival = 80 months, 95% CI = 71–88 months; low level: mean survival = 114 months, 95% CI = 101–129 months, p = 0.064; Figure 3A). In the GSE17538 cohort, high level mRNA expression of AXL was significantly correlated with shorter survival (high level: mean survival = 84 months, 95% CI = 72–96 months; low level: mean survival = 104 months, 95% CI = 94–114 months, p = 0.004; Figure 3B). By Cox-regression analysis, AXL mRNA expression was a predictor of patient survival in both the GSE14333 (HR = 1.839, 95% CI = 1.230–2.748, p = 0.003; Figure 3C) and the GSE17538 (HR = 2.158, 95% CI = 1.141–4.081, p = 0.018; Figure 3D) cohorts. In the GSE14333 cohort, AXL mRNA expression (HR = 1.631, 95% CI = 1.076–2.471, p = 0.021; Figure 3C) was an independent predictor of patient survival with tumor staging (p<0.001), while in the GSE17538 cohort, AXL mRNA expression (HR = 3.700, 95% CI = 1.665–8.220, p = 0.001; Figure 3D) was an independent predictor of patient survival together with stage (p<0.001) and grade (p = 0.055) of the cancers.

Similar results were obtained for CTGF. A high level of CTGF mRNA expression was significantly correlated with shorter patient survival in both GSE14333 (high level: mean survival = 87 months, 95% CI = 74–102 months; low level: mean survival = 98 months, 95% CI = 88–108 months, p = 0.012; Figure 3E) and GSE17538 (high level: mean survival = 85 months, 95% CI = 73–97 months; low level: mean survival = 105 months, 95% CI = 95–114 months, p = 0.004; Figure 3F) cohorts. Univariate Cox-regression analysis also showed that CTGF mRNA expression was a predictor for patient survival in both cohorts (GSE14333: HR = 1.667, 95% CI = 1.260–2.232, p<0.001; Figure 3G and GSE17538: HR = 1.433, 95% CI = 1.126–1.825, p = 0.004; Figure 3H). Moreover, CTGF mRNA expression is an independent predictor for survival in both cohorts (GSE14333: HR = 1.543, 95% CI = 1.155–2.061, p<0.001; Figure 3G and GSE17538: HR = 1.902, 95% CI = 1.385–2.612, p<0.001; Figure 3H) together with stage (GSE14333: p<0.001 and GSE17538: p<0.001) and grade (GSE17538: p = 0.053) of the cancers.

TAZ, AXL and CTGF can be used in combination to predict colon cancer patient survival

As described above, expression of TAZ, AXL and CTGF all correlated with colon cancer patient survival. We further investigated whether their expression could be used in combination as a prognostic marker for colon cancer patients. In the GSE14333 cohort, patients whose tumors had low level expression of the three genes had a mean survival time of 117 months (95% CI = 110–123 months), while those whose tumors overexpressed only one of the three genes also had a mean survival of 117 months (95% CI = 101–132 months). Patients whose tumors had a high level expression of two of the three genes had a mean survival time of 65 months (95% CI = 57–74 months), while those whose tumors had a high level expression of all three genes had a mean survival time of 72 months (95% CI = 60–84 months). The survival time of the patients between these four subgroups were significantly different (p = 0.001; Figure 4A). In Cox-regression analysis, using the subgroup containing patients whose tumors had a low expression of all three genes as a reference group, we found that patients whose tumors overexpressed one (HR = 3.953, 95% CI = 1.070–14.607, p = 0.039), two (HR = 6.503, 95% CI = 1.894–22.330, p = 0.003) or all the three genes (HR = 7.656, 95% CI = 2.287–25.628, p = 0.001) had an increasing risk for disease progression (Figure 4C).

[Figure omitted. See PDF.]

Figure 4. The associations between TAZ-AXL-CTGF expression and survival in colon cancer patients.

Patients were divided into 4 groups according to the number of genes that they overexpressed (expressed at above the Median level) among TAZ, AXL and CTGF. Kaplan-Meier analyses for TAZ-AXL-CTGF mRNA expression in (A) GSE14333 and (B) GSE17538 colon cancer patient datasets. Univariate and Multivariate Cox regression analyses for TAZ-AXL-CTGF mRNA expression, age, tumor stage and tumor grade in (C) GSE14333 and (D) GSE17538 colon cancer patient datasets.

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

Similar results were obtained in the GSE17538 patient cohort. Patients whose tumors expressed the three genes at low level had a mean survival time of 108 months (95% CI = 97–119 months), those whose tumors had a high level expression of one of the three genes had a mean survival time of 88 months (95% CI = 72–104 months), those whose tumors overexpressed two of the three genes had a mean survival time of 99 months (95% CI = 78–121 months), while those whose tumors overexpressed all three genes had a mean survival time of 77 months (95% CI = 63–91 months). Increasing incidence of overexpression of these three genes resulted in significantly shorter survival in these patients (Kaplan-Meier analysis, p = 0.01; Figure 4B). Cox-regression analysis using tumors that overexpressed none of the three genes as reference revealed that patients whose tumors overexpressed one (HR = 1.675, 95% CI = 0.75–3.741, p = 0.203), two (HR = 1.770, 95% CI = 0.790 = 3.966, p = 0.166) and all of the three genes (HR = 2.807, 95% CI = 1.482–5.317, p = 0.002) had an increased risk for disease progression (Figure 4D).

Indeed, patients whose tumors expressed TAZ, AXL and CTGF mRNA at low levels (GSE14333: mean survival = 117 months, 95% CI = 110–123; GSE17538: mean survival = 108 months, 95% CI = 95–121 months) had superior survival compared to those whose tumors overexpressed all the three genes (GSE14333: mean survival = 96, 95% CI = 84–108, p<0.001; GSE17538: mean survival = 91 months, 95% CI = 81–101 months, p = 0.037; Figure S1).

Neither the age nor sex of the patients was factor that determined the number of TAZ, AXL and CTGF that overexpressed in either colon cancer patient cohort. TAZ-AXL-CTGF expression was only significantly lower in tumors stage A compared to tumors of other stages (p = 0.022), but were not significantly different between other stages in GSE14333. In GSE17538, TAZ-AXL-CTGF expression was not significantly different between AJCC stages or tumor grades. When the patients were stratified by their stage, we found that TAZ-AXL-CTGF expression was still significantly correlated with poor survival in patients with higher stage tumors. In GSE14333, increasing number of genes overexpressed among TAZ, AXL and CTGF was significantly correlated with a shorter survival (HR = 1.414, 95% CI = 1.031–1.938, p = 0.032) in patients with stage C tumors. Similar results were obtained in GSE17538, in which, TAZ-AXL-CTGF expression was associated with poorer survival in patients with tumors of higher grade and higher stage. TAZ-AXL-CTGF expression was correlated with poorer survival in both grade B (HR = 1.315, 95% CI = 1.026 = 1.685, p = 0.031) and grade C (HR = 1.671, 95% CI = 1.001–2.790, p = 0.05) tumors, as well as AJCC stage 3 (HR = 1.694, 95% CI = 1.075–2.670, p = 0.023) and stage 4 (HR = 1.415, 95% CI = 1.074–1.866, p = 0.014) tumors.

Identification of genes that are differentially expressed between TAZ-AXL-CTGF-high and low patients

We found that 39 genes, including TAZ, AXL and CTGF, were significantly differentially expressed between TAZ-AXL-CTGF-high and low patients in both datasets. Some of these 39 genes have been shown to be induced by overexpressing YAP/TAZ-TEAD complexes in transfected cells [19], [29], including Akt3, CCDC80, FBN1, FRMD6, MSRB3, MYL9, PTRF, TAGLN, TNS1, VCAN and ZEB1 (Tables 1 and 2). Also co-regulated with TAZ-AXL-CTGF were genes involved in EMT, migration and invasion, colon cancer progression, calcium signaling, angiogenesis, cytoskeleton association, membrane trafficking, focal adhesion, Hippo pathway regulation, and MMPs inhibition (Tables 1 and 2), suggesting that these biological processes may play an important role in TAZ-AXL-CTGF-mediated colon cancer progression.

[Figure omitted. See PDF.]

Table 1. Genes that are co-regulated with TAZ-AXL-CTGF expression in GSE14333 colon cancer patient cohort.

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

[Figure omitted. See PDF.]

Table 2. Genes that are co-regulated with TAZ-AXL-CTGF expression in GSE17538 colon cancer patient cohort.

https://doi.org/10.1371/journal.pone.0054211.t002

Individually, as expected because of their strong association with the TAZ-AXL-CTGF expressions, all of the genes in the list (Tables 1 and 2) were significantly correlated with patient survival (p<0.05; data not shown). More importantly, we found that although none of the genes on the list were associated with patient survival in the TAZ-AXL-CTGF-high group of patients (data not shown), MGP, PDLIM3, TAGLN and ZEB2 were predictors of survival in the TAZ-AXL-CTGF-low group of patients, in the combined colon cancer patient cohort (Figure 5). When the two colon cancer patient datasets were combined, as expected, TAZ-AXL-CTGF mRNA expression levels were associated with patient survival (p<0.001; Figure 5A). A high level of MGP (p = 0.01; Figure 5B), PDLIM3 (p = 0.037; Figure 5C), TAGLN (p = 0.044; Figure 5D) and ZEB2 (p = 0.038; Figure 5E) mRNA expression were correlated with a shorter survival time in the TAZ-AXL-CTGF-low group of patients. These results suggest that these four genes may be used as prognostic markers for patients who express low levels of TAZ, AXL and CTGF mRNA.

[Figure omitted. See PDF.]

Figure 5. Kaplan-Meier analyses for the mRNA expression of TAZ-AXL-CTGF and its co-regulated genes in the combined colon cancer dataset.

(A) Kaplan-Meier analysis for TAZ-AXL-CTGF mRNA expression in the combined colon cancer datasets. Kaplan-Meier analyses for mRNA expression of (B) MGP, (C), PDLIM3, (D) TAGLN and (E) ZEB2 in colon cancer patients expressing low levels of TAZ-AXL-CTGF (blue line in (A)) in the combined colon cancer datasets.

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

Identification of potential small molecules that could target the 39 gene signature of TAZ-AXL-CTGF co-expression

Analysis to determine alterations in gene expression following small molecule treatment of cancer cells revealed a total of 257 small molecules that were associated with a gene expression signature significantly correlated with the TAZ-AXL-CTGF gene expression signature (the 25 genes that are presented in bold in Tables 1 and 2 together with TAZ, AXL and CTGF). Of the 257 small molecules, 164 had a 100% perturbation stability and 138 of them are inversely correlated with the TAZ-AXL-CTGF gene expression signature. The results are listed in the supplementary information. The top 20 small molecules were further analyzed through a Pubmed search regarding their effects on treatment of colon cancer. We found that amiloride and tretinoin have yielded 55 and 123 publications, respectively, when coupled with colon cancer in the search engine. Several publications have shown their inhibitory effect on colon cancer growth. Amiloride treatment has been shown to inhibit the growth of colon cancer cells in vitro [45] and in vivo [46]. Importantly, it can sensitize doxorubicin resistant colon cancer cells to treatment with doxorubicin [47], suggesting that amiloride and doxorubicin can be combined to treat doxorubicin resistant colon cancer. Tretinoin, also known as all-trans retinoic acid, has been shown to inhibit proliferation and anchorage-independent growth of colon cancer cells in vitro [48], [49] and in vivo [50], probably through regulating the differentiation state of cancer cells [51].

Identification of potential therapeutic targets for patients overexpressing TAZ-AXL-CTGF

We further investigated which genes could be used to further predict survival in TAZ-AXL-CTGF-high patients; these genes may be potential therapeutic targets specific for this group of patients whose colon cancers are more aggressive. Therefore, we compared gene expression in this group of patients between those patients who were still living and those who were deceased. The two cohorts' top 100 differentially expressed genes (based on survival status) were compared, revealing that ANO1 and SQLE are the two genes commonly differentially expressed between those TAZ-AXL-CTGF-high colon cancer patients who are still alive and those who are deceased. By Kaplan-Meier analysis, we found that these two genes themselves were associated with patient survival in both colon cancer datasets (Figure S2). We then investigated whether the associations between ANO1 or SQLE and patient survival were specific to patients who overexpressed TAZ, AXL and CTGF (i.e. the TAZ-AXL-CTGF-high subgroup). A low level of ANO1 expression was significantly (GSE14333, p = 0.002; Figure 6D and GSE17538, p = 0.007; Figure 6L) associated with better survival in TAZ-AXL-CTGF-high patients in both colon cancer datasets (GSE1433: high level, mean survival = 59 months, 95% CI = 44–73 months, vs. low level, mean survival = 100 months, 95% CI = 86–114 months; GSE17538: high level, mean survival = 66 months, 95% CI = 51–82 months, vs. low level, mean survival = 107 months, 95% CI = 87–127 months), but was not associated with survival in patients with other patterns of expression for TAZ-AXL-CTGF (Figure 6A, B, C, I, J, K). Similarly, a low level of SQLE expression was significantly (GSE14333, p = 0.02; Figure 6H and GSE17538, p = 0.01; Figure 6P) associated with better survival for TAZ-AXL-CTGF-high patients in both colon cancer datasets (GSE1433: high level, mean survival = 55 months, 95% CI = 40–70 months, vs. low level, mean survival = 88 months, 95% CI = 73–103 months; GSE17538: high level, mean survival = 62 months, 95% CI = 45–79 months, vs. low level, mean survival = 86 months, 95% CI = 71–101 months), but was not or less significant in patients with other patterns of expression for TAZ-AXL-CTGF (Figure 6E, F, G, M, N, O). When both ANO1 and SQLE genes were included in multivariate Cox-regression analysis, we found that mRNA expression of ANO1, but not SQLE, was an independent predictor of patient survival together with both tumor stage and TAZ-AXL-CTGF expression in both colon cancer patient datasets (Tables 3 and 4). These results suggest that ANO1 and SQLE mRNA expression may determine the aggressiveness of TAZ-AXL-CTGF-high tumors and that ANO1 mRNA expression could be used in combination with TAZ-AXL-CTGF and tumor stage for better prognostification. Our results also imply that blockade of ANO1 or SQLE mRNA expression or inhibition of these two proteins in this group of patients may prolong survival.

[Figure omitted. See PDF.]

Figure 6. The prognostic significance of ANO1 and SQLE mRNA expression in colon cancer patients stratified by their TAZ-AXL-CTGF mRNA expression.

Kaplan-Meier analyses for ANO1 mRNA expression in patients overexpressing (A) none, (B) one, (C) two and (D) three of the three genes (TAZ, AXL and CTGF) in the GSE14333 colon cancer patient datasets. Kaplan-Meier analyses for SQLE mRNA expression in patients overexpressing (E) none, (F) one, (G) two and (H) three of the three genes in the GSE14333 colon cancer patient datasets. Kaplan-Meier analyses for ANO1 mRNA expression in patients overexpressing (I) none, (J) one, (K) two and (L) three of the three genes in the GSE17538 colon cancer patient datasets. Kaplan-Meier analyses for SQLE mRNA expression in patients overexpressed in the GSE17538 colon cancer patient datasets.

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

[Figure omitted. See PDF.]

Table 3. Cox-regression analysis of GSE14333 datasets.

https://doi.org/10.1371/journal.pone.0054211.t003

[Figure omitted. See PDF.]

Table 4. Cox-regression analysis of GSE17538 datasets.

https://doi.org/10.1371/journal.pone.0054211.t004

The role of TAZ in colon cancer progression in vitro and in vivo

To confirm the bioinformatics analysis, we investigated the effect of TAZ knockdown in two colon cancer cell lines, HCT116 and SW620. TAZ knockdown in these two cell lines abolished the expression of TAZ and down-regulated AXL expression (Figure 7A). This result is in line with our finding in human specimens showing that AXL is a downstream target of TAZ. Knockdown of TAZ also resulted in a significant reduction in the number of colonies formed in both clonogenic and non-adherent soft-agar assays in these two cell lines (Figure 7B and C). Importantly, these in vitro results were replicated in in vivo tumorigenesis assay. Both HCT116-shTAZ and SW620-shTAZ cells formed significantly larger tumors in nude mice compared to HCT116-shScr and SW620-shScr cells, respectively (Figure 7D and 7E, respectively). Our results suggest that TAZ expression is required for a higher cell proliferation, non-adherent growth and tumorigenesis in colon cancer cells, traits that are associated with colon cancer progression.

[Figure omitted. See PDF.]

Figure 7. In vitro and in vivo assays for colon cancer cells expressing scramble shRNA or TAZ shRNA.

(A) Western blot showing that shTAZ specifically knockdown TAZ, but not YAP, and that AXL was down-regulated in shTAZ cells compared to shScr cells. (B) The clonogenicity and non-adherent growth of HCT116 cells expressing shScr or shTAZ were assessed and the number of colonies formed from three repeats was recorded (C) The clonogenicity and non-adherent growth of SW620 cells expressing shScr or shTAZ were assessed and the number of colonies formed from three repeats was recorded. (D) The in vivo tumorigenicity of HCT116 cells expressing shScr or shTAZ was assessed in nude nice, and the tumor formed was excised and weighted (n = 3 in each group). (E) The in vivo tumorigenicity of SW620 cells expressing shScr or shTAZ was assessed in nude nice, and the tumor formed was excised and weighted (n = 3 in each group).

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

Discussion

In the present study, we have shown that TAZ mRNA expression is positively correlated with two of its downstream targets, AXL and CTGF, and that TAZ is significantly associated with poor survival of colon cancer patients in two independent colon cancer datasets, comprising 522 patients. Interestingly, the upregulation of AXL and CTGF, which reflects the increased transcriptional activity of TAZ-TEAD complexes, can be used in combination with TAZ mRNA expression, for better prognostification in these two independent colon cancer patient datasets. Genes that are co-regulated with TAZ-AXL-CTGF overexpression are involved in several important cellular processes, including cell migration, angiogenesis and calcium signaling, as well as others that have already been described as prognostic markers for colon cancer progression. These genes may be upstream factors or downstream effectors of TAZ and the dysregulated Hippo pathway in colon cancers. Importantly, we also identified two potential therapeutic targets, ANO1 and SQLE, for patients with upregulated TAZ-AXL-CTGF expression; downregulation of either of these two genes may greatly improve the survival of TAZ-AXL-CTGF-high colon cancer patients.

An increase in TAZ mRNA may not necessarily correlate with an increase of its transcriptional activity, due to the fact that TAZ could be post-translationally regulated via cytoplasmic sequestration by 14-3-3 [1]. Possible disparity between TAZ mRNA expression and its transcriptional activity may impair the lone utilization of TAZ mRNA expression as a prognostic marker. Previously, we have shown that AXL [43] and CTGF [44] are both downstream targets of the Hippo pathway. In the present study, we have analyzed the transcriptional outcome of TAZ-TEAD complexes in colon cancer using AXL or CTGF alone or in combination with TAZ and found that patients who co-overexpressed all three genes had significantly poorer survival compared to those who had other expression patterns. Our results, therefore, strongly suggest that TAZ, AXL and CTGF can be used in combination for prognostification in colon cancer patients.

Several genes that are related to EMT are overexpressed in high TAZ-AXL-CTGF expressing patients (Tables 1 and 2). These include ACTA2 [52], ZEB1 [53], [54] and ZEB2 [55], [56]. Interestingly, ZEB1 has been shown to be a downstream target of TAZ in retinal pigment epithelial cells [57], suggesting that ZEB1 may act as a downstream effector of TAZ to promote cancer metastasis, while ZEB1 and ZEB2 have also been shown to be prognostic markers in colon cancer [53], [56]. Genes that govern migration and invasion were also differentially expressed (Tables 1 and 2). AKT3, which has been shown to contribute oncogenic functions similar to other AKT isoforms [58], was up-regulated in the TAZ-AXL-CTGF-high group of patients. DDR2 [59], [60], FERMT2 [61], [62], [63] and VCAN [64] also play significant roles in cell adhesion and migration and are upregulated in TAZ-AXL-CTGF-high group of patients.

Colon cancer biomarkers were also differentially expressed in these two groups of patients (Tables 1 and 2), namely EFEMP2, a serum biomarker for early detection of colon cancer [65] and SULF1, a protein important in colon cancer diagnosis [66] and whose serum level is elevated in patients with colon adenomas [67]. Genes that are implicated in angiogenesis were also identified (Tables 1 and 2); both VCAN [68], [69] and ANTXR1 [70], [71], [72] promote angiogenesis for cancer progression. Four genes that are involved in calcium binding or signaling were co-regulated with TAZ-AXL-CTGF in the colon cancer specimens, including FBN1 [73], CALD1 [74], MGP [75] and MYL9 [76].

Interestingly, genes that play a role in regulating the Hippo pathway were also co-regulated with TAZ-AXL-CTGF. The angiomotin family members act as tumor suppressors by inhibiting the oncogenic functions of YAP and TAZ [44], [77], while FRMD6 also acts as an antagonist of YAP by activating Hippo pathway kinases [78]. In this study, we found that AMOLT1 and FRMD6 are co-overexpressed in TAZ-AXL-CTGF positive tumors, suggesting that AMOLT1 and FRMD6 may form a negative regulatory loop with TAZ activation, which requires further investigation in vitro in colon cancer cell line models. In addition, VGLL3, which has been shown to act as a co-activator of TEAD transcription factors is also upregulated in TAZ-AXL-CTGF-high tumors [79].

Treatments targeting TAZ-AXL-CTGF-high cancers are required due to the aggressive nature of this type of colon cancer. We employed two different analyses to identify potential therapeutic agents and gene targets for this type of cancer. In sccMap, we found that amiloride and tretinoin may inhibit the gene expression signature associated with this aggressive type of colon cancer; these molecules have been shown in the literature to provide strong inhibitory effects on colon cancer proliferation in vitro and in vivo. The other 18 small molecules of the top twenty identified by sccMap as listed in the supplementary information have yet to be studied in term of their effect on colon cancer progression driven by the overexpression of TAZ-AXL-CTGF co-overexpression.

ANO1, also named DOG1, has been shown to be ubiquitously expressed in gastrointestinal stromal tumors [80] and its overexpression is correlated with development of distant metastases in head and neck squamous cell carcinoma [81]. Recently, ANO1 was shown to promote tumorigenesis and cancer progression via activating MAPK [82]. In the same study, the authors demonstrated that pharmacological inhibition of ANO1 resulted in cancer cell death, however the role of ANO1 in colon cancer progression has not been examined. In the present study, we found that a low level expression of ANO1 in the aggressive TAZ-AXL-CTGF-high subgroup of colon cancer patients is associated with prolonged survival (high level to low level of ANO1 expression; mean survival from 59 to 100 months and from 66 to 107 month in GSE14333 and GSE17538 patient datasets, respectively), suggesting that pharmacological inhibition of ANO1 may represent a novel therapeutic approach for this group of patients with aggressive colon cancer. ANO1 is a calcium ion-activated chloride channel. Interestingly, multiple genes involved in calcium signaling (FBN1, CALD1, MGP and MYL9) were co-overexpressed with TAZ-AXL-CTGF in colon cancer. Due to the fact that calcium signaling plays an important role in cancer progression [83], [84], [85], further investigation into the relationship between ANO1 and calcium signaling in TAZ-AXL-CTGF-mediated cancer progression is warranted. Nonetheless, our results provide a clue to the involvement of calcium signaling and Ca2+-activated CI− channels in colon cancer progression mediated by the Hippo pathway.

Little is known about how SQLE promotes cancer progression. It has been shown to be overexpressed in lung squamous cell carcinoma by suppression subtractive hydridization [86], in pancreatic cancer by genome-wide analysis using microarray based techniques [87], and in prostate cancer progression by bioinformatics analysis [88]. Its overexpression in breast cancer is correlated with decreased distant metastasis-free survival [89]. These results show that SQLE promotes cancer progression in multiple types of cancer. Again, its role in colon cancer progression was undefined. In the present study, we also found that a high level expression of SQLE is associated with poorer survival in TAZ-AXL-CTGF-high patients to be associated with poorer survival (high level to low level of SQLE expression; mean survival from 55 to 88 months and from 62 to 102 month in GSE14333 and GSE17538 patient datasets, respectively), suggesting that SQLE may be a novel therapeutic target for this group of patients.

In conclusion, this study has shown that TAZ-AXL-CTGF in combination may be a novel prognostic indicator for colon cancer progression, and that their overexpression is associated with increased expression of genes that are associated with colon cancer progression. Furthermore, ANO1 and SQLE overexpression may further define a poorer prognosis for colon cancer patients overexpressing TAZ-AXL-CTGF.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Colon cancer patients expressing low levels of TAZ, AXL and CTGF had superior survival. Patients were stratified into two groups; those whose tumors expressed TAZ, AXL and CTGF mRNA at low level (solid line) and those whose tumors expressed at least one of TAZ, AXL and CTGF at high level (dotted line). Kaplan-Meier analyses for these two subgroups of patients in (A) GSE14333 and (B) GSE17538 colon cancer datasets.

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

(TIF)

Figure S2.

The associations between ANO1 or SQLE, and survival in colon cancer patients. Kaplan-Meier analyses for (A) ANO1 and (B) SQLE mRNA expression in the GSE14333 colon cancer patient dataset. Kaplan-Meier analyses for (C) ANO1 and (D) SQLE mRNA expression in the GSE17538 colon cancer patient dataset.

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

(TIF)

Author Contributions

Conceived and designed the experiments: HFY SDZ WH. Performed the experiments: HFY CMM YHH JMT XZ SDZ. Analyzed the data: HFY SDZ WH. Contributed reagents/materials/analysis tools: QZ SDZ WH. Wrote the paper: HFY CMM SDZ WH.

Citation: Yuen H-F, McCrudden CM, Huang Y-H, Tham JM, Zhang X, Zeng Q, et al. (2013) TAZ Expression as a Prognostic Indicator in Colorectal Cancer. PLoS ONE 8(1): e54211. https://doi.org/10.1371/journal.pone.0054211

References

1. Chan SW, Lim CJ, Chen L, Chong YF, Huang C, et al. (2011) The Hippo pathway in biological control and cancer development. J Cell Physiol 226: 928–939.

2. Zhao B, Li L, Lei Q, Guan KL (2010) The Hippo-YAP pathway in organ size control and tumorigenesis: an updated version. Genes Dev 24: 862–874.

3. Pan D (2010) The hippo signaling pathway in development and cancer. Dev Cell 19: 491–505.

4. Harvey K, Tapon N (2007) The Salvador-Warts-Hippo pathway - an emerging tumour-suppressor network. Nat Rev Cancer 7: 182–191.

5. Overholtzer M, Zhang J, Smolen GA, Muir B, Li W, et al. (2006) Transforming properties of YAP, a candidate oncogene on the chromosome 11q22 amplicon. Proc Natl Acad Sci U S A 103: 12405–12410.

6. Zhao B, Wei X, Li W, Udan RS, Yang Q, et al. (2007) Inactivation of YAP oncoprotein by the Hippo pathway is involved in cell contact inhibition and tissue growth control. Genes Dev 21: 2747–2761.

7. Zhang J, Ji JY, Yu M, Overholtzer M, Smolen GA, et al. (2009) YAP-dependent induction of amphiregulin identifies a non-cell-autonomous component of the Hippo pathway. Nat Cell Biol 11: 1444–1450.

8. Neto-Silva RM, de Beco S, Johnston LA (2010) Evidence for a growth-stabilizing regulatory feedback mechanism between Myc and Yorkie, the Drosophila homolog of Yap. Dev Cell 19: 507–520.

9. Lian I, Kim J, Okazawa H, Zhao J, Zhao B, et al. (2010) The role of YAP transcription coactivator in regulating stem cell self-renewal and differentiation. Genes Dev 24: 1106–1118.

10. Hall CA, Wang R, Miao J, Oliva E, Shen X, et al. (2010) Hippo pathway effector Yap is an ovarian cancer oncogene. Cancer Res 70: 8517–8525.

11. Zhang X, George J, Deb S, Degoutin JL, Takano EA, et al. (2011) The Hippo pathway transcriptional co-activator, YAP, is an ovarian cancer oncogene. Oncogene 30: 2810–2822.

12. Muramatsu T, Imoto I, Matsui T, Kozaki K, Haruki S, et al. (2011) YAP is a candidate oncogene for esophageal squamous cell carcinoma. Carcinogenesis 32: 389–398.

13. Liu AM, Poon RT, Luk JM (2010) MicroRNA-375 targets Hippo-signaling effector YAP in liver cancer and inhibits tumor properties. Biochem Biophys Res Commun 394: 623–627.

14. Yuan M, Tomlinson V, Lara R, Holliday D, Chelala C, et al. (2008) Yes-associated protein (YAP) functions as a tumor suppressor in breast. Cell Death Differ 15: 1752–1759.

15. Basu S, Totty NF, Irwin MS, Sudol M, Downward J (2003) Akt phosphorylates the Yes-associated protein, YAP, to induce interaction with 14-3-3 and attenuation of p73-mediated apoptosis. Mol Cell 11: 11–23.

16. Ehsanian R, Brown M, Lu H, Yang XP, Pattatheyil A, et al. (2010) YAP dysregulation by phosphorylation or DeltaNp63-mediated gene repression promotes proliferation, survival and migration in head and neck cancer subsets. Oncogene 29: 6160–6171.

17. Lei QY, Zhang H, Zhao B, Zha ZY, Bai F, et al. (2008) TAZ promotes cell proliferation and epithelial-mesenchymal transition and is inhibited by the hippo pathway. Mol Cell Biol 28: 2426–2436.

18. Chan SW, Lim CJ, Guo K, Ng CP, Lee I, et al. (2008) A role for TAZ in migration, invasion, and tumorigenesis of breast cancer cells. Cancer Res 68: 2592–2598.

19. Zhang H, Liu CY, Zha ZY, Zhao B, Yao J, et al. (2009) TEAD transcription factors mediate the function of TAZ in cell growth and epithelial-mesenchymal transition. J Biol Chem 284: 13355–13362.

20. Hong JH, Hwang ES, McManus MT, Amsterdam A, Tian Y, et al. (2005) TAZ, a transcriptional modulator of mesenchymal stem cell differentiation. Science 309: 1074–1078.

21. Varelas X, Sakuma R, Samavarchi-Tehrani P, Peerani R, Rao BM, et al. (2008) TAZ controls Smad nucleocytoplasmic shuttling and regulates human embryonic stem-cell self-renewal. Nat Cell Biol 10: 837–848.

22. Zhao D, Zhi X, Zhou Z, Chen C (2012) TAZ antagonizes the WWP1-mediated KLF5 degradation and promotes breast cell proliferation and tumorigenesis. Carcinogenesis 33: 59–67.

23. Lai D, Ho KC, Hao Y, Yang X (2011) Taxol resistance in breast cancer cells is mediated by the hippo pathway component TAZ and its downstream transcriptional targets Cyr61 and CTGF. Cancer Res 71: 2728–2738.

24. Zhou Z, Hao Y, Liu N, Raptis L, Tsao MS, et al. (2011) TAZ is a novel oncogene in non-small cell lung cancer. Oncogene 30: 2181–2186.

25. Xie M, Zhang L, He CS, Hou JH, Lin SX, et al. (2012) Prognostic Significance of TAZ Expression in Resected Non-Small Cell Lung Cancer. J Thorac Oncol

26. Cordenonsi M, Zanconato F, Azzolin L, Forcato M, Rosato A, et al. (2011) The Hippo transducer TAZ confers cancer stem cell-related traits on breast cancer cells. Cell 147: 759–772.

27. de Cristofaro T, Di Palma T, Ferraro A, Corrado A, Lucci V, et al. (2011) TAZ/WWTR1 is overexpressed in papillary thyroid carcinoma. Eur J Cancer 47: 926–933.

28. Chan SW, Lim CJ, Loo LS, Chong YF, Huang C, et al. (2009) TEADs mediate nuclear retention of TAZ to promote oncogenic transformation. J Biol Chem 284: 14347–14358.

29. Zhao B, Ye X, Yu J, Li L, Li W, et al. (2008) TEAD mediates YAP-dependent gene induction and growth control. Genes Dev 22: 1962–1971.

30. Chen L, Chan SW, Zhang X, Walsh M, Lim CJ, et al. (2010) Structural basis of YAP recognition by TEAD4 in the hippo pathway. Genes Dev 24: 290–300.

31. Li Z, Zhao B, Wang P, Chen F, Dong Z, et al. (2010) Structural insights into the YAP and TEAD complex. Genes Dev 24: 235–240.

32. Tian W, Yu J, Tomchick DR, Pan D, Luo X (2010) Structural and functional analysis of the YAP-binding domain of human TEAD2. Proc Natl Acad Sci U S A 107: 7293–7298.

33. Steinhardt AA, Gayyed MF, Klein AP, Dong J, Maitra A, et al. (2008) Expression of Yes-associated protein in common solid tumors. Hum Pathol 39: 1582–1589.

34. Konsavage WM Jr, Kyler SL, Rennoll SA, Jin G, Yochum GS (2012) Wnt/beta-catenin signaling regulates Yes-associated protein (YAP) gene expression in colorectal carcinoma cells. J Biol Chem 287: 11730–11739.

35. Avruch J, Zhou D, Bardeesy N (2012) YAP oncogene overexpression supercharges colon cancer proliferation. Cell Cycle 11: 1090–1096.

36. Zhou D, Zhang Y, Wu H, Barry E, Yin Y, et al. (2011) Mst1 and Mst2 protein kinases restrain intestinal stem cell proliferation and colonic tumorigenesis by inhibition of Yes-associated protein (Yap) overabundance. Proc Natl Acad Sci U S A 108: E1312–1320.

37. Jorissen RN, Gibbs P, Christie M, Prakash S, Lipton L, et al. (2009) Metastasis-Associated Gene Expression Changes Predict Poor Outcomes in Patients with Dukes Stage B and C Colorectal Cancer. Clin Cancer Res 15: 7642–7651.

38. Smith JJ, Deane NG, Wu F, Merchant NB, Zhang B, et al. (2010) Experimentally derived metastasis gene expression profile predicts recurrence and death in patients with colon cancer. Gastroenterology 138: 958–968.

39. Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, et al. (2006) The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science 313: 1929–1935.

40. Lamb J (2007) The Connectivity Map: a new tool for biomedical research. Nat Rev Cancer 7: 54–60.

41. Zhang SD, Gant TW (2009) sscMap: an extensible Java application for connecting small-molecule drugs using gene-expression signatures. BMC Bioinformatics 10: 236.

42. McArt DG, Zhang SD (2011) Identification of candidate small-molecule therapeutics to cancer by gene-signature perturbation in connectivity mapping. PLoS One 6: e16382.

43. Xu MZ, Chan SW, Liu AM, Wong KF, Fan ST, et al. (2011) AXL receptor kinase is a mediator of YAP-dependent oncogenic functions in hepatocellular carcinoma. Oncogene 30: 1229–1240.

44. Chan SW, Lim CJ, Chong YF, Pobbati AV, Huang C, et al. (2011) Hippo pathway-independent restriction of TAZ and YAP by angiomotin. J Biol Chem 286: 7018–7026.

45. Szolgay-Daniel E, Carlsson J, Zierold K, Holtermann G, Dufau E, et al. (1991) Effects of amiloride treatment on U-118 MG and U-251 MG human glioma and HT-29 human colon carcinoma cells. Cancer Res 51: 1039–1044.

46. Tatsuta M, Iishi H, Baba M, Uehara H, Nakaizumi A (1995) Chemoprevention by amiloride of experimental carcinogenesis in rat colon induced by azoxymethane. Carcinogenesis 16: 941–942.

47. Pannocchia A, Revelli S, Tamponi G, Giorgianni A, Todde R, et al. (1996) Reversal of doxorubicin resistance by the amiloride analogue EIPA in multidrug resistant human colon carcinoma cells. Cell Biochem Funct 14: 11–18.

48. Hoosein NM, Brattain DE, McKnight MK, Childress KE, Chakrabarty S, et al. (1988) Comparison of the antiproliferative effects of transforming growth factor-beta, N,N-dimethylformamide and retinoic acid on a human colon carcinoma cell line. Cancer Lett 40: 219–232.

49. Kim EJ, Kang YH, Schaffer BS, Bach LA, MacDonald RG, et al. (2002) Inhibition of Caco-2 cell proliferation by all-trans retinoic acid: role of insulin-like growth factor binding protein-6. J Cell Physiol 190: 92–100.

50. Narahara H, Tatsuta M, Iishi H, Baba M, Uedo N, et al. (2000) K-ras point mutation is associated with enhancement by deoxycholic acid of colon carcinogenesis induced by azoxymethane, but not with its attenuation by all-trans-retinoic acid. Int J Cancer 88: 157–161.

51. Rai K, Sarkar S, Broadbent TJ, Voas M, Grossmann KF, et al. (2010) DNA demethylase activity maintains intestinal cells in an undifferentiated state following loss of APC. Cell 142: 930–942.

52. Kalluri R, Weinberg RA (2009) The basics of epithelial-mesenchymal transition. J Clin Invest 119: 1420–1428.

53. Spaderna S, Schmalhofer O, Hlubek F, Berx G, Eger A, et al. (2006) A transient, EMT-linked loss of basement membranes indicates metastasis and poor survival in colorectal cancer. Gastroenterology 131: 830–840.

54. Wellner U, Schubert J, Burk UC, Schmalhofer O, Zhu F, et al. (2009) The EMT-activator ZEB1 promotes tumorigenicity by repressing stemness-inhibiting microRNAs. Nat Cell Biol 11: 1487–1495.

55. Beltran M, Puig I, Pena C, Garcia JM, Alvarez AB, et al. (2008) A natural antisense transcript regulates Zeb2/Sip1 gene expression during Snail1-induced epithelial-mesenchymal transition. Genes Dev 22: 756–769.

56. Kahlert C, Lahes S, Radhakrishnan P, Dutta S, Mogler C, et al. (2011) Overexpression of ZEB2 at the invasion front of colorectal cancer is an independent prognostic marker and regulates tumor invasion in vitro. Clin Cancer Res 17: 7654–7663.

57. Liu Y, Xin Y, Ye F, Wang W, Lu Q, et al. (2010) Taz-tead1 links cell-cell contact to zeb1 expression, proliferation, and dedifferentiation in retinal pigment epithelial cells. Invest Ophthalmol Vis Sci 51: 3372–3378.

58. Mende I, Malstrom S, Tsichlis PN, Vogt PK, Aoki M (2001) Oncogenic transformation induced by membrane-targeted Akt2 and Akt3. Oncogene 20: 4419–4423.

59. Labrador JP, Azcoitia V, Tuckermann J, Lin C, Olaso E, et al. (2001) The collagen receptor DDR2 regulates proliferation and its elimination leads to dwarfism. EMBO Rep 2: 446–452.

60. Olaso E, Labrador JP, Wang L, Ikeda K, Eng FJ, et al. (2002) Discoidin domain receptor 2 regulates fibroblast proliferation and migration through the extracellular matrix in association with transcriptional activation of matrix metalloproteinase-2. J Biol Chem 277: 3606–3613.

61. Montanez E, Ussar S, Schifferer M, Bosl M, Zent R, et al. (2008) Kindlin-2 controls bidirectional signaling of integrins. Genes Dev 22: 1325–1330.

62. Tu Y, Wu S, Shi X, Chen K, Wu C (2003) Migfilin and Mig-2 link focal adhesions to filamin and the actin cytoskeleton and function in cell shape modulation. Cell 113: 37–47.

63. Shen Z, Ye Y, Dong L, Vainionpaa S, Mustonen H, et al. (2012) Kindlin-2: a novel adhesion protein related to tumor invasion, lymph node metastasis, and patient outcome in gastric cancer. Am J Surg 203: 222–229.

64. Said N, Sanchez-Carbayo M, Smith SC, Theodorescu D (2012) RhoGDI2 suppresses lung metastasis in mice by reducing tumor versican expression and macrophage infiltration. J Clin Invest 122: 1503–1518.

65. Yao L, Lao W, Zhang Y, Tang X, Hu X, et al. (2012) Identification of EFEMP2 as a Serum Biomarker for the Early Detection of Colorectal Cancer with Lectin Affinity Capture Assisted Secretome Analysis of Cultured Fresh Tissues. J Proteome Res

66. Babel I, Barderas R, Diaz-Uriarte R, Moreno V, Suarez A, et al. (2011) Identification of MST1/STK4 and SULF1 proteins as autoantibody targets for the diagnosis of colorectal cancer by using phage microarrays. Mol Cell Proteomics 10: M110 001784.

67. Matusiewicz M, Krzystek-Korpacka M, Diakowska D, Grabowski K, Augoff K, et al. (2008) Serum sulfatase activity is more elevated in colonic adenomas than cancers. Int J Colorectal Dis 23: 383–387.

68. Koyama H, Hibi T, Isogai Z, Yoneda M, Fujimori M, et al. (2007) Hyperproduction of hyaluronan in neu-induced mammary tumor accelerates angiogenesis through stromal cell recruitment: possible involvement of versican/PG-M. Am J Pathol 170: 1086–1099.

69. Zheng PS, Wen J, Ang LC, Sheng W, Viloria-Petit A, et al. (2004) Versican/PG-M G3 domain promotes tumor growth and angiogenesis. FASEB J 18: 754–756.

70. Chaudhary A, Hilton MB, Seaman S, Haines DC, Stevenson S, et al. (2012) TEM8/ANTXR1 blockade inhibits pathological angiogenesis and potentiates tumoricidal responses against multiple cancer types. Cancer Cell 21: 212–226.

71. Fernando S, Fletcher BS (2009) Targeting tumor endothelial marker 8 in the tumor vasculature of colorectal carcinomas in mice. Cancer Res 69: 5126–5132.

72. Duan HF, Hu XW, Chen JL, Gao LH, Xi YY, et al. (2007) Antitumor activities of TEM8-Fc: an engineered antibody-like molecule targeting tumor endothelial marker 8. J Natl Cancer Inst 99: 1551–1555.

73. Kettle S, Yuan X, Grundy G, Knott V, Downing AK, et al. (1999) Defective calcium binding to fibrillin-1: consequence of an N2144S change for fibrillin-1 structure and function. J Mol Biol 285: 1277–1287.

74. Mayanagi T, Sobue K (2011) Diversification of caldesmon-linked actin cytoskeleton in cell motility. Cell Adh Migr 5: 150–159.

75. Farzaneh-Far A, Proudfoot D, Weissberg PL, Shanahan CM (2000) Matrix gla protein is regulated by a mechanism functionally related to the calcium-sensing receptor. Biochem Biophys Res Commun 277: 736–740.

76. Szczesna-Cordary D, Guzman G, Zhao J, Hernandez O, Wei J, et al. (2005) The E22K mutation of myosin RLC that causes familial hypertrophic cardiomyopathy increases calcium sensitivity of force and ATPase in transgenic mice. J Cell Sci 118: 3675–3683.

77. Zhao B, Li L, Lu Q, Wang LH, Liu CY, et al. (2011) Angiomotin is a novel Hippo pathway component that inhibits YAP oncoprotein. Genes Dev 25: 51–63.

78. Angus L, Moleirinho S, Herron L, Sinha A, Zhang X, et al. (2012) Willin/FRMD6 expression activates the Hippo signaling pathway kinases in mammals and antagonizes oncogenic YAP. Oncogene 31: 238–250.

79. Helias-Rodzewicz Z, Perot G, Chibon F, Ferreira C, Lagarde P, et al. (2010) YAP1 and VGLL3, encoding two cofactors of TEAD transcription factors, are amplified and overexpressed in a subset of soft tissue sarcomas. Genes Chromosomes Cancer 49: 1161–1171.

80. West RB, Corless CL, Chen X, Rubin BP, Subramanian S, et al. (2004) The novel marker, DOG1, is expressed ubiquitously in gastrointestinal stromal tumors irrespective of KIT or PDGFRA mutation status. Am J Pathol 165: 107–113.

81. Ayoub C, Wasylyk C, Li Y, Thomas E, Marisa L, et al. (2010) ANO1 amplification and expression in HNSCC with a high propensity for future distant metastasis and its functions in HNSCC cell lines. Br J Cancer 103: 715–726.

82. Duvvuri U, Shiwarski DJ, Xiao D, Bertrand C, Huang X, et al. (2012) TMEM16A, induces MAPK and contributes directly to tumorigenesis and cancer progression. Cancer Res

83. Clapham DE (2007) Calcium signaling. Cell 131: 1047–1058.

84. Monteith GR, McAndrew D, Faddy HM, Roberts-Thomson SJ (2007) Calcium and cancer: targeting Ca2+ transport. Nat Rev Cancer 7: 519–530.

85. Prevarskaya N, Skryma R, Shuba Y (2011) Calcium in tumour metastasis: new roles for known actors. Nat Rev Cancer 11: 609–618.

86. Liu Y, Sun W, Zhang K, Zheng H, Ma Y, et al. (2007) Identification of genes differentially expressed in human primary lung squamous cell carcinoma. Lung Cancer 56: 307–317.

87. Harada T, Chelala C, Crnogorac-Jurcevic T, Lemoine NR (2009) Genome-wide analysis of pancreatic cancer using microarray-based techniques. Pancreatology 9: 13–24.

88. Zhao H, Logothetis CJ, Gorlov IP (2010) Usefulness of the top-scoring pairs of genes for prediction of prostate cancer progression. Prostate Cancer Prostatic Dis 13: 252–259.

89. Helms MW, Kemming D, Pospisil H, Vogt U, Buerger H, et al. (2008) Squalene epoxidase, located on chromosome 8q24.1, is upregulated in 8q+ breast cancer and indicates poor clinical outcome in stage I and II disease. Br J Cancer 99: 774–780.