do not necessarily reflect the views of UKDiss.com.
Stromal TRIM28-associated signalling pathway modulation within the colorectal cancer microenvironment
Background: Stromal gene expression patterns predict patient outcomes in colorectal cancer. TRIM28 is a transcriptional co-repressor that regulates an abundance of genes through the KRAB domain family of transcription factors. We have previously shown that stromal expression of TRIM28 is a marker of disease relapse and poor survival in colorectal cancer. Here, we perform a differential epithelium-stroma proteomic network analysis to characterise signalling pathways associated with TRIM28 within the tumour microenvironment.
Methods: Reverse phase protein arrays were generated from laser capture micro-dissected carcinoma and stromal cells from fresh frozen colorectal cancer tissues. Phosphorylation and total protein levels were measured for 30 cancer-related signalling pathway endpoints. Strength and direction of associations between signalling endpoints were identified using Spearman’s rank-order correlation analysis and compared to TRIM28 levels. Expression status of TRIM28 in tumour epithelium and stromal fibroblasts was assessed using IHC in formalin fixed tissue and the epithelium to stroma Protein Expression Ratio method.
Results: We found distinct proteomic networks in the epithelial and stromal compartments which were linked to expression levels of TRIM28. Low levels of TRIM28 in tumour stroma (High Epithelium:Stroma Ratio) were found in 10 out of 19 cases. Upon proteomic network analyses, these stromal High Ratio cases revealed moderate signalling pathway similarity exemplified by 59 significant Spearman correlations (ρ > 0.75, P ≤ 0.01). Furthermore, low levels of stromal TRIM28 correlated with elevated MDM2 levels in tumour epithelium (p=0.01) and COX-2 levels in tumour stroma (p=0.002). Low TRIM28 epithelium to stroma ratios were associated with elevated levels of caspases 3 and 7 in stroma (p=0.041 and p=0.036) and an increased signalling pathway similarity in stromal cells with 82 significant Spearman correlations (ρ > 0.75, P ≤ 0.01).
Conclusions: By dissecting TRIM28-associated pathways in stromal fibroblasts and epithelial tumour cells, we performed comprehensive proteomic analyses of molecular networks within the tumour microenvironment.We found modulation of several signalling pathways associated with TRIM28, which may be attributed to the pleiotropic properties of TRIM28 through its translational suppression of the family of KRAB-domain transcription factors in tumour stromal compartments.
Keywords: TRIM28; colorectal cancer; epithelium; stroma; prognosis, tumour antigen
The cross-talk between epithelial cells and the non-epithelial cellular component of tumour stroma exerts substantial influence on the severity and aggressiveness of cancer (Quail and Joyce 2013). Stromal fibroblasts, in particular, contribute to the cross-talk through secretion of cytokines and growth factors thereby triggering signaling pathways in tumour cells (Bhowmick, Neilson et al. 2004, Calon, Espinet et al. 2012). In colorectal cancer, tumour stroma is known to impact on patient outcomes (Mesker, Junggeburt et al. 2007), with stromal genetic instability being a key contributing factor (Ishiguro, Yoshida et al. 2006). Recent transcriptomic profiling studies confirmed these earlier findings and have defined extensive stromal gene expression profiles linked to poor prognosis and to the induction of epithelial–mesenchymal transition (EMT) through TGF-ß signalling in colorectal cancer (Calon, Lonardo et al. 2015). Since EMT is associated with distinct gene profiles in colorectal cancer patients with poor outcomes (De Sousa, Wang et al. 2013) and is thought to be regulated by a wide range of transcription factors (Lamouille, Xu et al. 2014), the impact of the stromal compartment on the tumour microenvironment might be determined by transcriptional control of a wide range of genes in stromal fibroblasts.
Tripartite motif-containing 28 (TRIM28) is a nuclear corepressor involved in transcriptional regulation of a large number of transcription factors belonging to the Krüppel-associated box (KRAB) repressor domain family of zinc finger proteins (Peng, Feldman et al. 2002, Huntley, Baggott et al. 2006). TRIM28 has been recently shown to exhibit transcriptional control over a plethora of regulatory networks and cell re-programming pathways (Brattas, Jonsson et al. 2017, Czerwinska, Shah et al. 2017, Miles, de Vries et al. 2017). TRIM28 is an essential component of the fibroblast transcription site-1 activator transcription complex required for induction of EMT (Venkov, Link et al. 2007) and has been shown in cancer to further contribute to EMT via regulation of E- and N-cadherins (Chen, Munoz-Antonia et al. 2014). In addition, TRIM28 can mediate communication between cells through regulation of senescence in fibroblasts causing the induction of secretory phenotype which includes the secretion of interleukins and other pro-inflammatory molecules (Santos and Gil 2014).
We have previously shown that TRIM28 is overexpressed in colorectal cancer (Kijanka, Hector et al. 2010) and that its stromal expression is an independent marker of disease relapse and patient survival (Fitzgerald, Sheehan et al. 2013). More recently, tumour promoting effects of TRIM28 were observed in glioma and breast cancer (Qi, Cai et al. 2016, Wei, Cheng et al. 2016). In this paper, we sought to further investigate how differences in expression patterns of TRIM28 in epithelial and stromal tissue compartments in colorectal cancer can influence patient survival. Since regulatory mechanisms of the corepressor TRIM28 are likely to be imposed on pan-transcriptional level through co-repression or co-activation of over 700 distinct KRAB-domain transcription factors (Huntley, Baggott et al. 2006, Imbeault, Helleboid et al. 2017), in this study we take a proteomic systems approach which allows us to investigate signalling pathway activation differentially in tumour epithelium and stromal compartments. We here employ laser capture microdissection to isolate enriched populations of epithelial and stromal cells and perform a reverse phase protein array-based proteomic analysis of selected cancer signalling endpoints to investigate the underlying molecular interactions at play in relation to TRIM28 within the tumour microenvironment.
Materials and methods
Patient cohort and tissue specimens
The study was approved by Ethics (Medical) Research Committee at Beaumont Hospital, Dublin, Ireland. Informed consent was obtained from all patients. Patients with a history of cancer and neo-adjuvant treatment were excluded. A total of 19 cases with late-stage (Stage III and IV) colorectal cancer (CRC) were investigated. All patients were diagnosed with CRC between 2012 and 2013. The median age of the patients at the time of diagnosis was 67 (range 47-88 years) with 10 male and 9 female patients. In total, 15 patients had colonic carcinoma, whilst 4 had rectal carcinomas. Clinical and pathological parameters of the patient cohort are shown in Table 1.
Both, fresh-frozen and formalin-fixed and paraffin embedded (FFPE) tissue samples were obtained for each case. FFPE tissues were used for immunohistochemical analysis, while fresh-frozen tissues were used for laser capture microdissection (LCM) and reverse phase protein array (RPPA) analysis. A pathologist identified and collected areas of invasive carcinoma from the tumour mass for fresh-frozen preservation and FFPE. Local resection and standard fixation protocols were carried out in all cases. Each FFPE block was sectioned and stained with haematoxylin and eosin (H&E) and graded by a consultant pathologist to confirm pathological stage and grade of the tumours. All fresh tissue samples for RPPA analysis were snap-frozen in liquid nitrogen and processed uniformly and rapidly to ensure preservation of molecular endpoints. The time from removal of a colectomy specimen to snap-freezing of samples was < 20 mins. Fresh-frozen tissue samples were stored at −80°C.
Immunohistochemistry staining and assessment
Immunohistochemistry (IHC) and scoring were carried out using an anti-TRIM28 rabbit monoclonal antibody (mAb) (C42G12, Cell Signalling Technology Inc., Danvers, MA, USA), as previously described (Fitzgerald, Sheehan et al. 2013). The degree of nuclear TRIM28 staining was evaluated for epithelial and stromal tissue separately and scored as follows: absence of staining (score = 0), weak staining (score = 1+), moderate staining (score = 2+) and strong staining (score = 3+). TRIM28 epithelium to stroma ratios were computed based on a Protein Expression Ratio method using IHC scores as previously described (Fitzgerald, Sheehan et al. 2013). Briefly, high TRIM28 Expression Ratio was defined as at least 2 units of difference in staining intensity (e.g. epithelium strong [3+] and stroma weak [1+], or epithelium moderate [2+] and stroma negative ). Low TRIM28 Expression Ratio was defined as 1 or 0 units of difference in staining intensity (e.g. epithelium moderate [2+] and stroma weak [1+], or epithelium weak [1+] and stroma weak [1+]). A previous study in our lab has shown that the inter-observer variability of IHC scoring is as low as 7% (Kay, Barry Walsh et al. 1996). In cases where there were discrepancies between the scorers, a consensus was reached after a joint review using a multi-headed microscope.
Laser capture microdissection and reverse phase protein arrays
Laser capture microdissection (LCM) was performed to isolate separate populations of epithelial and stromal cells for cell signalling analysis as described previously (Emmert-Buck, Bonner et al. 1996, Espina, Heiby et al. 2007, Fitzgerald, Sheehan et al. 2015). Briefly, consecutive 8µm thick frozen tissue sections were cut using a cryostat for each sample and mounted on glass slides. Then, using an infrared-based laser capture system (Arcturus XT, Life Technologies, San Francisco, CA, USA), approximately 20,000 epithelial and stromal cells were removed for each frozen tissue sample. To account for possible heterogeneity within the tissue sample, multiple separate areas within each patient sample were microdissected and no attempt was made to target specific regions within the tumour. Tissue processing and preparation of tissue lysates have been described previously (Wulfkuhle, Aquino et al. 2003, Sheehan, Calvert et al. 2005).
Reverse phase protein arrays (RPPA) were generated as previously described (Espina, Wulfkuhle et al. 2007, Mueller, Liotta et al. 2010, Fitzgerald, Sheehan et al. 2015). Each array contained epithelial and stromal cell lysates for all 19 cases and each lysate was printed in a 2-fold dilution curve representing undiluted lysate, 1:2, 1:4 and 1:8 dilutions. Control lysates were printed in a 2-fold dilution curve. All RPPAs were baked for 2 hours at 80°C to allow fixation and then stored with desiccant at −20°C. Quality control samples were printed on the RPPA to ensure protein deposition and immunostaining reactivity (Gulmann, Sheehan et al. 2006). These included A431 cell lines (±EGF stimulation; BD Pharmingen, San Diego, CA, USA) and a pooled sample of CRC cases.
RPPA immunostaining, image acquisition and data analysis
RPPA slides were blocked (I-Block, Applied Biosystems) for 2 hours prior to immunostaining. Immunostaining was conducted on a Dako Autostainer (Catalysed Signal Amplification (CSA) kit, Dako). In total, 30 primary antibodies specific to known signalling endpoints were used to measure phosphorylation and protein levels using RPPAs (see Supplementary Table 1). Each slide was incubated with a single primary antibody at room temperature for 30 minutes. The negative control slide was incubated with antibody diluent (Dako) instead of a primary antibody. For normalization purposes, a slide was incubated with anti-ssDNA antibody (1:15,000; IBL International GmbH). The secondary antibodies were goat anti-rabbit (1:10,000; Vector Laboratories), or rabbit anti-mouse IgG (1:10; Dako) depending on the host species of the primary antibody. Amplification was achieved using horseradish peroxidase-mediated biotinyl tyramide with chromogenic detection (diaminobenzidine; Dako). In order to determine the total protein concentration in each sample, two RPPA slides were stained with Sypro Ruby Protein Blot Stain (Molecular Probes, Eugene, OR, USA) and visualized with NovaRay Image Acquisition Software (Alpha Innotech, San Leandro, CA, USA).
RPPA slides were scanned on a flatbed scanner (UMAX 2100XL; white balance 255, black 0, middle tone 1.37, 600 dpi, 14 bit) and spot intensity was analysed using the Image Quant v5.2 software package (Molecular Dynamics). To normalize data, the relative intensity for each protein spot was divided by the ssDNA relative intensity for the corresponding spot (Chiechi, Mueller et al. 2012) and data reduction was carried out using a VBA Excel macro, RPPA Analysis Suite (Mueller 2013). One case did not meet the quality criteria as the spot intensities were not significantly different from the negative control and therefore was excluded from proteomic network analysis. All staining intensities are shown in Supplementary Table 2.
Proteomic network analysis
Proteomic network analysis was performed as previously shown (Chiechi, Novello et al. 2013, Fitzgerald, Sheehan et al. 2015). Spearman ρ correlation analysis with ρ ≥ 0.75 and p ≤ 0.01 was used to build proteomic network graphs (Gephi 0.9, The Gephi Consortium, Paris, France, www.gephi.org). Proteomic networks were drafted based on 18 CRC cases (one excluded) and classified into two groups based on their epithelial to stromal TRIM28 Expression Ratio score; TRIM28 Low Ratio group (n=9) representing good patient prognosis and TRIM28 High Ratio group (n=9) representing poor patient prognosis (Fitzgerald, Sheehan et al. 2013). In order to identify the differences in the proteomic signature of the cohorts, significant correlations that were common in both TRIM28 Ratio groups in the same tissue compartment (epithelial or stromal) were excluded and only exclusive significant correlations were used to construct the proteomic network graphs. For example, correlations in the epithelial compartment that were common in both TRIM28 High and Low ratio cohorts, were excluded and only correlations that were unique to the TRIM28 High or Low Ratio groups were used to construct the proteomic network graphs for each group.
All tests were analysed using SPSS 21.0 software (SPSS, Chicago, IL, USA). The Spearman rank correlation coefficient, ρ, was calculated for each protein pair in the RPPA cohort; ρ ≥ 0.75 with p ≤ 0.01 was considered significant. An independent-samples t-test was used to compare the expression levels of each endpoint between groups, the findings were considered statistically significant at p < 0.05.
We have previously shown that the assessment of TRIM28 expression levels, when it encompasses epithelial and stromal tissue compartments, can serve an independent marker of relapse and patient survival in colorectal cancer (Fitzgerald, Sheehan et al. 2013). In order to assess stromal and epithelial TRIM28 levels in this cohort, we separately scored TRIM28 in epithelial carcinoma cells and stromal fibroblasts in 19 FFPE tissues using IHC. Strong epithelial staining (3+) for TRIM28 was found in 8 cases, moderate staining (2+) in 9 and weak staining (1+) in 2. Cells in the stromal compartment showed moderate TRIM28 staining (2+) in 2 cases, weak staining (1+) in 12 cases and absence of staining (0) in 5 cases (Figure 1).
The tissue-compartment Protein Expression Ratio method developed for reciprocal tumour microenvironment assessment (Fitzgerald, Sheehan et al. 2013), allowed us to divide the 19-case CRC cohort in 2 groups; the TRIM28 High Epithelium:Stroma Ratio group, which is mainly characterised by low levels of TRIM28 in tumour stroma (Figure 1A) and the Low Epithelium:Stroma Ratio group, which shows more similar TRIM28 expression levels in both, epithelial and stromal compartments (Figure 1B). Using the IHC scoring data, the Ratio method resulted in 10 cases being assigned to the TRIM28 High Ratio group and 9 cases to the Low Ratio group (Figure 1C). In the High Ratio group, one case had high epithelial TRIM28 levels and absence of TRIM28 in stroma expression, 6 cases had high levels of epithelial TRIM28 and weak staining in the stroma and 3 cases had moderate levels of TRIM28 in the tumour epithelium with an absence of staining for TRIM28 in stroma (Figure 1C). Most cases in the TRIM28 Low Ratio group had weak to moderate expression in both, epithelium and stroma. One case had strong epithelial expression and another case presented with absence of stromal TRIM28 staining. The distribution of IHC scores allotted to both TRIM28 Expression Ratio groups is shown in Figure 1 C.
Figure 1.Epithelial to stromal TRIM28 Expression Ratios in colorectal cancer. Panel (A) shows moderate (2+) expression of TRIM28 in epithelial cells and absence (0) of staining in stromal compartments resulting in a high epithelial to stromal TRIM28 Expression Ratio score (x200). Panel (B) shows moderate (2+) expression of TRIM28 in both epithelial and stromal cells resulting in a low epithelial to stromal TRIM28 Expression Ratio score (x200). Panel (C) shows a graph of distributions of all TRIM28 Expression Ratios cases. Absence of staining (0), weak (1+), moderate (2+), or strong (3+) TRIM28 staining intensities were found in both epithelial and stromal tissue compartments. The TRIM28 Expression Ratios are labelled by (○) epithelial and (■) stromal symbols, which are linked with straight lines for matched cases. The numbers above the connecting lines depict the total number of cases for each specific distribution of epithelial to stromal TRIM28 Expression Ratios. All tissues were stained using IHC with an anti-TRIM28 C42G12 antibody.
High Ratio TRIM28 stromal networks are linked to tumour progression
We next sought to understand the underlying molecular interplay associated with TRIM28 Expression Ratios and signalling pathways linked to cancer. RPPAs were immunostained against 30 signalling endpoints and the resulting quantitative data were used for proteomic network analyses in 9 TRIM28 High Ratio and 9 TRIM28 Low Ratio cases (Supplementary Table 2).The Spearman’s rank-order correlation analysis revealed that 55 of 435 possible protein pair correlations were shared (ρ > 0.75, P ≤ 0.01) between stroma of TRIM28 High and Low Ratio cases (Supplementary Table 4).
In addition, the High Ratio TRIM28 stromal proteomic network revealed a moderate number of 59 exclusive signalling endpoints associations, which were not present in TRIM28 Low Ratio stoma (Supplementary Table 4). As shown in Figure 2A, the stromal High Ratio network is divided into three sub-networks. Each sub-network is dominated by nodes linked to tumour progression, including Acetyl-CoA, LC3B and COX-2 as prominent endpoints.
Low Ratio TRIM28 stromal networks are linked to apoptosis
The Low Ratio TRIM28 stromal network showed an overall higher number of exclusive signalling endpoints with 82 significant associations when compared to High Ratio TRIM28 stoma (Supplementary Table 4). As shown in Figure 2B, the network is also divided in three sub-groups, these however are considerably different to the High Ratio networks and are dominated by endpoints linked to apoptosis, such as the activated caspases 3 & 7 and PP2A. Such association with apoptotic proteins is not seen in the TRIM28 High Ratio stromal network and shows that the proteomic architecture of the stromal tissue compartment is significantly different in TRIM28 High and Low Ratio patients and hence may account for differences in patient outcomes seen in both colorectal cancer groups.
Figure 2: Stromal proteomic networks for TRIM28 High (A) and Low (B) Ratio cases. Both, TRIM28 High (A) and TRIM28 Low Ratio (B) stromal tissue proteomic networks consist of three main sub-networks (green, pink and light green). Dominant nodes in the High Ratio sub-networks are proteins linked to tumour progression and nodes prominent in the Low ratio sub-networks are proteins linked to apoptosis. Nodes represent proteins quantified by RPPA; the bigger the node, the more significant correlations relative to that protein. Each line connecting 2 nodes represents a significant correlation between the nodes; the thicker the line, the higher the Spearman ρ correlation. Proteins are grouped on the basis of Spearman ρ values and the number of connections among a group of nodes; strongly correlated nodes are represented close to each other and with the same colour.
High Ratio TRIM28 epithelial networks are associated with cellular pro-survival pathways
In order to obtain a complete picture of potential molecular tumour microenvironment characteristics, we performed an equivalent Spearman rank correlation mapping for the epithelial compartments of the patient cohort. The 30-endpoint Spearman rank analysis revealed that 87 of 435 possible protein pair correlations were shared between epithelia of TRIM28 High and Low Ratio cases (Supplementary Table 3), indicating more similar signalling patterns between both groups when compared to the stroma counterpart.
Fifty-five protein pairs correlated (ρ > 0.75, P ≤ 0.01) exclusively in the epithelium of TRIM28 High Ratio cases and were entirely absent in TRIM28 Low Ratio cases (Supplementary Table 3). Detailed graphical representation of those High Ratio-exclusive associations of protein expression and phosphorylation in epithelial compartments is shown in Figure 3A. The TRIM28 High Ratio epithelial proteomic network is associated with several pro-survival endpoints including cell-cell adhesion (E-Cadherin) and autophagy (AKT, Survivin).
Low Ratio TRIM28 epithelial networks are associated with pro-apoptotic signalling
The proteomic network analysis of tumour epithelium in the TRIM28 Low Ratio group revealed statistically significant protein pair correlations in 69 signalling endpoints (ρ > 0.75, P ≤ 0.01) (Supplementary Table 3). Those pairs were exclusive to the epithelium of TRIM28 Low Ratio cases and were absent in TRIM28 Low Ratio epithelium (Fig. 3B). The proteomic network is dominated by pro-apoptotic proteins represented, among others, by activated Caspase-7 and the anti-apoptotic regulator BCL-2. Those nodes and protein correlations are noticeably diminished in the High Ratio epithelial proteomic network. In addition, the proto-oncogene Lck, another prominent node in the TRIM28 Low Ratio epithelial network, was also identified as a dominant node in TRIM28 Low Ratio stromal network (Figure 2B & 3B).
Figure 3: Epithelial proteomic networks for TRIM28 High (A) and Low (B) Ratio cases. Both,TRIM28 High (A) and TRIM28 Low Ratio (B) epithelial tumour tissue proteomic networks consist of three main sub-networks (green, pink and light green). Dominant nodes in the High Ratio networks represent pro-survival pathways. Prominent nodes in the Low Ratio proteomic networks represent pro-apoptotic pathways. Nodes represent proteins quantified by RPPA; the bigger the node, the more significant correlations relative to that molecule. Each line connecting 2 nodes represents a significant correlation between the nodes; the thicker the line, the higher the Spearman ρ correlation. Proteins are grouped on the basis of Spearman ρ values and the number of connections among a group of nodes; strongly correlated nodes are represented close to each other and with the same colour.
Stromal MDM2 and COX-2 levels in High Ratio TRIM28 colorectal cancer
We next compared expression levels of 30 signalling proteins between tumour epithelial and stromal compartments in either TRIM28 High or Low Ratio cases. Our analysis revealed that expression levels of MDM2 are significantly lower in stroma of TRIM28 High Ratio cases when compared to corresponding epithelial compartments (p = 0.010, Fig. 4A). No statistically significant differences in MDM2 levels were found in Low Ratio cases (p = 0.923, Fig. 4A). In contrast, stromal expression levels of COX-2 were significantly elevated in High Ratio cases when compared to tumour epithelium (p = 0.002, Fig. 4B). There is no significant differences in COX-2 intensities between epithelium and stroma of TRIM28 Low Ratio cases (p = 0.465, Fig. 4B). These findings further substantiate a molecular disparity between TRIM28 High and Low Ratio cases. Our earlier findings, for instance, show that TRIM28 may interact with MDM2 to influence the levels of the tumour suppressor p53 and this may account for their significantly poorer 5-year overall survival and 5-year ‘recurrence-free’ survival seen in TRIM28 High Ratio colorectal cancer patients (Fitzgerald, Sheehan et al. 2013).
Figure 4: High Ratio TRIM28 stromal tissue shows low levels of MDM2 and high levels of COX-2. Box-plot diagrams of (A) MDM2 and (B) COX-2 represent RPPA intensity levels in the epithelium and stroma of TRIM28 High and Low Ratio cases. MDM2 is downregulated in stroma of TRIM28 High Ratio cases (p = 0.010, A). COX2 is elevated in stroma in TRIM28 High Ratio cases were found (p = 0.002, B). The median (line within the box), mean (centre of the diamond), 25th and 75th percentiles and maximum and minimum values are all displayed on each box-plot.
Stromal Caspases 3 and 7 levels in High Ratio TRIM28 colorectal cancer
The expression levels of all investigated signalling endpoints were then analysed across epithelial and stromal data regardless of their TRIM28 Ratio scores. We found significantly less activated caspase-3 and caspase-7 expressed in the stroma of TRIM28 High Ratio patients compared to stroma of Low Ratio patients (p = 0.041, Fig. 5A; p = 0.036, Fig. 5B). Activated caspases are the main executioners of apoptosis in the cell. Since the reduction of stromal apoptosis was previously shown to be an independent prognostic factor for poorer overall survival and disease-recurrence in CRC (Koelink, Sier et al. 2009), the lower levels of activated caspases seen in the stroma may be another contributing factor to the poorer outcomes in colorectal cancer patients with TRIM28 High Ratios.
Figure 5: Caspase-3 and Caspase-7 protein levels are lower in the stroma of TRIM28 High Ratio patients. (A)Caspase-3 and (B) Caspase-7 RPPA intensity levels in the epithelium and stroma of TRIM28 High and Low Ratio patients are represented by box-plot diagrams. There is significantly less active caspase-3 and -7present in the stroma of TRIM28 High Ratio patients than in the stroma of TRIM28 low Ratio patients (p = 0.04 and p = 0.036, respectively). The median (line within the box), mean (centre of the diamond), 25th and 75th percentiles and maximum and minimum values are all displayed on each box-plot.
TRIM28 is a marker of disease relapse and poor survival in colorectal cancer that engages different molecular mechanism in context of the tumour microenvironment. Since TRIM28 can regulate several hundred different genes in its role a transcriptional co-repressor, we investigated its role within the tumour microenvironment using a differential proteomics approach analysing LCM-enriched patient tissue. The analysis of signalling pathways in stromal fibroblasts and epithelial tumour cells revealed distinctive molecular networks associated with TRIM28. The modulation of signalling pathways found in this study may be attributed the pleiotropic properties of TRIM28, specifically impacting tumour-associated stromal compartments.
Altered protein expression in stromal cells has be shown in many malignancies including lung, prostate, breast cancer and colorectal cancer to be prognostic of patient survival (Ogawa, Takenaka et al. 2004, Koelink, Sier et al. 2009, Sloan, Ciocca et al. 2009, Hagglof, Hammarsten et al. 2010). We have previously shown that high TRIM28 expression ratios between stromal and epithelial compartments of colorectal carcinomas are associated with poor patient survival and disease recurrence (Fitzgerald, Sheehan et al. 2013). The findings of this current study suggest that the aggressiveness of TRIM28 High Ratio cancers in mediated through activation of tumour progression pathways in stromal compartments and cellular pro-survival pathways in tumour epithelium.
The pro-survival mechanisms associated with TRIM28 High Ratio epithelial compartments were highlighted central involvement of cell-cell adhesion (E-cadherin) and autophagy (AKT, survivin) pathways. Activated AKT mediates downstream responses, including cell survival, growth, proliferation, cell migration and angiogenesis by phosphorylating a range of intracellular proteins. Survivin is a member of the apoptosis (IAP) inhibitor family that inhibits caspase activation, thereby leading to negative regulation of apoptosis or programmed cell death. Activated caspase 3 & 7 are noticeably diminished in this protein interaction network. In contrast, the TRIM28 Low Ratio epithelial proteomic network is associated with pro-apoptotic proteins including activated Caspase-7 and the apoptotic regulator BCL-2.
Our earlier studies propose that during the epithelial to mesenchymal transformation of tumour cells, some carcinoma cells can take on characteristics of stromal fibroblasts (Sheehan, Gulmann et al. 2007). EMT is thought to be required physiologically during embryogenesis, but its persistence in tumour cells is suggested to play a role in the promotion of an invasive phenotype. Fibroblasts produced by EMT express a gene encoding fibroblast-specific protein 1 (FSP1), which is regulated by a promoter element called fibroblast transcription site–1 (FTS-1) (Strutz, Okada et al. 1995). TRIM28 interacts with FTS-1 resulting in the transcriptional activation of genes encoding the EMT proteome suggesting that TRIM28 plays a key role in fibroblast formation (Venkov, Link et al. 2007). During EMT epithelial cells may lose many epithelial phenotypes including the loss of specialized cell-cell contacts, which is related to the downregulation of E-cadherin (Xiong, Hong et al. 2012, Yu, Sun et al. 2015) and E-cadherin is the dominant node in the High Ratio epithelial network. Taken together, these results suggest that some cancer cells in these patients have undergone EMT.
We demonstrate that there is significantly more epithelial than stromal MDM2 expression in TRIM28 High Ratio patients and this is in line with our previous finding that TRIM28 is overexpressed in epithelial cancer tissue. MDM2 plays an important role in the regulation of p53 by: (i) by inhibiting its transcriptional activity, (ii) by controlling its subcellular localization, and (iii) by modulating its protein stability (Elabd, Meroni et al. 2016). A number of TRIM proteins modulate the abundance and the activity of p53 and/or MDM2 and, therefore, can influence p53’s tumour-suppressive activity (Yue, Liu et al. 2014). TRIM28 interacts with MDM2 resulting in enhanced deacetylation of p53 and reduced transcriptional activity (Wang, Ivanov et al. 2005). antagonism of the MDM2-p53 interaction activates p53 signalling leading to a regression in human tumours in preclinical cancer models (Tovar, Graves et al. 2013). Our results suggest that high epithelial expression of both TRIM28 and MDM2, as seen in TRIM28 High Ratio patients, results in the decrease of p53’s tumour supressing activity. This finding is relevant as it supports our previous finding that there is less TRIM28 in the stroma of TRIM28 High Ratio patients. MDM2 is a RING domain ubiquitin E3 ligase and a major regulator of the tumour suppressor p53. Importantly, TRIM28 was previously identified as an MDM2-binding protein and shown to form a complex with MDM2 and p53 in vivo (Wang, Ivanov et al. 2005). The binding is mediated by the N-terminal coiled-coil domain of TRIM28 and the central acidic domain of MDM2.
Lck in one of the dominant nodes in both of the TRIM28 Low Ratio networks, and is absent from both of the TRIM28 High Ratio networks. Lck overexpression has previously been linked to improved survival (Elsberger, Fullerton et al. 2010). In recent years, cancer immunotherapy has been gaining greater significance in the treatment of cancer, in particular in adoptive cell transfer (ACT). The principle of ACT is that T-cells are genetically engineered to recognise antigens on the surface of cancer cells and destroy these cells (Palucka and Banchereau 2012). Lck plays an essential role in the selection and maturation of developing T-cells in the thymus and in the function of mature T-cells (Nika, Soldani et al. 2010). The role of Lck in T-cell immunity may affect cancer cells in a similar manner to ACT, thereby contributing to low cancer recurrence rates seen in TRIM28 Low Ratio patients (Fitzgerald, Sheehan et al. 2013).
COX-2 plays an important role in the inflammatory response and is generally unexpressed under normal conditions in most cells, but elevated levels are found during inflammation. The overexpression of COX-2 is associated with various types of cancer and the expression levels are generally proportional to their aggressiveness (Ristimäki, Sivula et al. 2002, Peng, Zhou et al. 2013). We have previously shown that the expression levels of COX-2 is related to lymph node metastasis, advanced Dukes staging, and poorer long-term outcome for patients with colorectal cancer (Sheehan, Sheahan et al. 1999). In this study we show that COX-2 is overexpressed in the stroma of TRIM28 High Ratio patients. Low COX-2 expression in epithelial cells and upregulation in stroma has also been shown to be indicative of tumour progression in laryngeal squamous cell carcinoma (Kourelis, Vandoros et al. 2009). Taken together these results highlight the importance of evaluating the expression levels of biomarkers in the tumour microenvironment, as well as the epithelial cells, as the tumour microenvironment can also have substantial prognostic significance.
In conclusion, our findings suggest that the proteomic architecture of the stromal tissue compartment is significantly different in TRIM28 High and Low Ratio patients, favouring apoptotic processes in the Low Ratio patients, which may be ultimately linked to better outcome seen in this patient cohort. By dissecting the effects of TRIM28 in stromal fibroblasts and epithelial tumour cells, we were able to elucidate the complex relationship between stromal and epithelial compartments in CRC and the proteomic data discovered offers insights into the tumour suppressive and tumour promoting effects of the TRIM28 protein in the microenvironment of CRC patients.
Brattas, P. L., M. E. Jonsson, L. Fasching, J. Nelander Wahlestedt, M. Shahsavani, R. Falk, A. Falk, P. Jern, M. Parmar and J. Jakobsson (2017). “TRIM28 Controls a Gene Regulatory Network Based on Endogenous Retroviruses in Human Neural Progenitor Cells.” Cell Rep 18(1): 1-11.
Calon, A., E. Espinet, S. Palomo-Ponce, D. V. Tauriello, M. Iglesias, M. V. Cespedes, M. Sevillano, C. Nadal, P. Jung, X. H. Zhang, D. Byrom, A. Riera, D. Rossell, R. Mangues, J. Massague, E. Sancho and E. Batlle (2012). “Dependency of colorectal cancer on a TGF-beta-driven program in stromal cells for metastasis initiation.” Cancer Cell 22(5): 571-584.
Calon, A., E. Lonardo, A. Berenguer-Llergo, E. Espinet, X. Hernando-Momblona, M. Iglesias, M. Sevillano, S. Palomo-Ponce, D. V. Tauriello, D. Byrom, C. Cortina, C. Morral, C. Barcelo, S. Tosi, A. Riera, C. S. Attolini, D. Rossell, E. Sancho and E. Batlle (2015). “Stromal gene expression defines poor-prognosis subtypes in colorectal cancer.” Nat Genet 47(4): 320-329.
Chiechi, A., C. Mueller, K. M. Boehm, A. Romano, M. S. Benassi, P. Picci, L. A. Liotta and V. Espina (2012). “Improved data normalization methods for reverse phase protein microarray analysis of complex biological samples.” Biotechniques 0(0): 1-7.
Chiechi, A., C. Novello, G. Magagnoli, E. F. Petricoin, 3rd, J. Deng, M. S. Benassi, P. Picci, I. Vaisman, V. Espina and L. A. Liotta (2013). “Elevated TNFR1 and serotonin in bone metastasis are correlated with poor survival following bone metastasis diagnosis for both carcinoma and sarcoma primary tumors.” Clin Cancer Res 19(9): 2473-2485.
Czerwinska, P., P. K. Shah, K. Tomczak, M. Klimczak, S. Mazurek, B. Sozanska, P. Biecek, K. Korski, V. Filas, A. Mackiewicz, J. N. Andersen and M. Wiznerowicz (2017). “TRIM28 multi-domain protein regulates cancer stem cell population in breast tumor development.” Oncotarget 8(1): 863-882.
De Sousa, E. M. F., X. Wang, M. Jansen, E. Fessler, A. Trinh, L. P. de Rooij, J. H. de Jong, O. J. de Boer, R. van Leersum, M. F. Bijlsma, H. Rodermond, M. van der Heijden, C. J. van Noesel, J. B. Tuynman, E. Dekker, F. Markowetz, J. P. Medema and L. Vermeulen (2013). “Poor-prognosis colon cancer is defined by a molecularly distinct subtype and develops from serrated precursor lesions.” Nat Med 19(5): 614-618.
Elsberger, B., R. Fullerton, S. Zino, F. Jordan, T. J. Mitchell, V. G. Brunton, E. A. Mallon, P. G. Shiels and J. Edwards (2010). “Breast cancer patients’ clinical outcome measures are associated with Src kinase family member expression.” Br J Cancer 103(6): 899-909.
Espina, V., J. D. Wulfkuhle, V. S. Calvert, E. F. Petricoin and L. A. Liotta (2007). Reverse Phase Protein Microarrays for Monitoring Biological Responses. Cancer Genomics and Proteomics: Methods and Protocols. P. B. Fisher. Totowa, NJ, Humana Press: 321-336.
Fitzgerald, S., K. M. Sheehan, V. Espina, A. O’Grady, R. Cummins, D. Kenny, L. Liotta, R. O’Kennedy, E. W. Kay and G. S. Kijanka (2015). “High CerS5 expression levels associate with reduced patient survival and transition from apoptotic to autophagy signalling pathways in colorectal cancer.” The Journal of Pathology: Clinical Research 1(1): 54-65.
Fitzgerald, S., K. M. Sheehan, A. O’Grady, D. Kenny, R. O’Kennedy, E. W. Kay and G. S. Kijanka (2013). “Relationship between epithelial and stromal TRIM28 expression predicts survival in colorectal cancer patients.” J Gastroenterol Hepatol 28(6): 967-974.
Fitzgerald, S., K. M. Sheehan, A. O’Grady, D. Kenny, R. O’Kennedy, E. W. Kay and G. S. Kijanka (2013). “Relationship between epithelial and stromal TRIM28 expression predicts survival in colorectal cancer patients.” Journal of Gastroenterology and Hepatology 28(6): 967-974.
Gulmann, C., K. M. Sheehan, E. W. Kay, L. A. Liotta and E. F. Petricoin, 3rd (2006). “Array-based proteomics: mapping of protein circuitries for diagnostics, prognostics, and therapy guidance in cancer.” J Pathol 208(5): 595-606.
Hagglof, C., P. Hammarsten, A. Josefsson, P. Stattin, J. Paulsson, A. Bergh and A. Ostman (2010). “Stromal PDGFRbeta expression in prostate tumors and non-malignant prostate tissue predicts prostate cancer survival.” PLoS One 5(5): e10747.
Huntley, S., D. M. Baggott, A. T. Hamilton, M. Tran-Gyamfi, S. Yang, J. Kim, L. Gordon, E. Branscomb and L. Stubbs (2006). “A comprehensive catalog of human KRAB-associated zinc finger genes: insights into the evolutionary history of a large family of transcriptional repressors.” Genome Res 16(5): 669-677.
Kay, E. W., C. J. Barry Walsh, D. Whelan, A. O’Grady and M. B. Leader (1996). “Inter-observer variation of p53 immunohistochemistry–an assessment of a practical problem and comparison with other studies.” Br J Biomed Sci 53(2): 101-107.
Kijanka, G., S. Hector, E. W. Kay, F. Murray, R. Cummins, D. Murphy, B. D. MacCraith, J. H. Prehn and D. Kenny (2010). “Human IgG antibody profiles differentiate between symptomatic patients with and without colorectal cancer.” Gut 59(1): 69-78.
Koelink, P. J., C. F. M. Sier, D. W. Hommes, C. B. H. W. Lamers and H. W. Verspaget (2009). “Clinical significance of stromal apoptosis in colorectal cancer.” British Journal of Cancer 101(5): 765-773.
Kourelis, K., G. Vandoros, T. Kourelis, T. Papadas, P. Goumas and G. Sotiropoulou-Bonikou (2009). “Low COX2 in tumor and upregulation in stroma mark laryngeal squamous cell carcinoma progression.” Laryngoscope 119(9): 1723-1729.
Mesker, W. E., J. M. Junggeburt, K. Szuhai, P. de Heer, H. Morreau, H. J. Tanke and R. A. Tollenaar (2007). “The carcinoma-stromal ratio of colon carcinoma is an independent factor for survival compared to lymph node status and tumor stage.” Cell Oncol 29(5): 387-398.
Miles, D. C., N. A. de Vries, S. Gisler, C. Lieftink, W. Akhtar, E. Gogola, I. Pawlitzky, D. Hulsman, E. Tanger, M. Koppens, R. L. Beijersbergen and M. van Lohuizen (2017). “TRIM28 is an Epigenetic Barrier to Induced Pluripotent Stem Cell Reprogramming.” Stem Cells 35(1): 147-157.
Mueller, D. C. (2013). Reverse Phase Protein Microarray Analysis Suite, http://capmm.gmu.edu/rpma-analysis-suite.
Nika, K., C. Soldani, M. Salek, W. Paster, A. Gray, R. Etzensperger, L. Fugger, P. Polzella, V. Cerundolo, O. Dushek, T. Höfer, A. Viola and O. Acuto (2010). “Constitutively Active Lck Kinase in T Cells Drives Antigen Receptor Signal Transduction.” Immunity 32(6): 766-777.
Ogawa, E., K. Takenaka, K. Yanagihara, M. Kurozumi, T. Manabe, H. Wada and F. Tanaka (2004). “Clinical significance of VEGF-C status in tumour cells and stromal macrophages in non-small cell lung cancer patients.” Br J Cancer 91(3): 498-503.
Peng, H., I. Feldman and F. J. Rauscher Iii (2002). “Hetero-oligomerization Among the TIF Family of RBCC/TRIM Domain-containing Nuclear Cofactors: A Potential Mechanism for Regulating the Switch Between Coactivation and Corepression.” Journal of Molecular Biology 320(3): 629-644.
Peng, L., Y. Zhou, Y. Wang, H. Mou and Q. Zhao (2013). “Prognostic Significance of COX-2 Immunohistochemical Expression in Colorectal Cancer: A Meta-Analysis of the Literature.” PLoS ONE 8(3): e58891.
Ristimäki, A., A. Sivula, J. Lundin, M. Lundin, T. Salminen, C. Haglund, H. Joensuu and J. Isola (2002). “Prognostic significance of elevated cyclooxygenase-2 expression in breast cancer.” Cancer research 62(3): 632-635.
Sheehan, K. M., V. S. Calvert, E. W. Kay, Y. Lu, D. Fishman, V. Espina, J. Aquino, R. Speer, R. Araujo, G. B. Mills, L. A. Liotta, E. F. Petricoin, 3rd and J. D. Wulfkuhle (2005). “Use of reverse phase protein microarrays and reference standard development for molecular network analysis of metastatic ovarian carcinoma.” Mol Cell Proteomics 4(4): 346-355.
Sheehan, K. M., C. Gulmann, G. S. Eichler, J. N. Weinstein, H. L. Barrett, E. W. Kay, R. M. Conroy, L. A. Liotta and E. F. Petricoin, III (2007). “Signal pathway profiling of epithelial and stromal compartments of colonic carcinoma reveals epithelial-mesenchymal transition.” Oncogene 27(3): 323-331.
Sloan, E. K., D. R. Ciocca, N. Pouliot, A. Natoli, C. Restall, M. A. Henderson, M. A. Fanelli, F. D. Cuello-Carrion, F. E. Gago and R. L. Anderson (2009). “Stromal cell expression of caveolin-1 predicts outcome in breast cancer.” Am J Pathol 174(6): 2035-2043.
Tovar, C., B. Graves, K. Packman, Z. Filipovic, B. H. M. Xia, C. Tardell, R. Garrido, E. Lee, K. Kolinsky, K.-H. To, M. Linn, F. Podlaski, P. Wovkulich, B. Vu and L. T. Vassilev (2013). “MDM2 Small-Molecule Antagonist RG7112 Activates p53 Signaling and Regresses Human Tumors in Preclinical Cancer Models.” Cancer Research 73(8): 2587-2597.
Venkov, C. D., A. J. Link, J. L. Jennings, D. Plieth, T. Inoue, K. Nagai, C. Xu, Y. N. Dimitrova, F. J. Rauscher and E. G. Neilson (2007). “A proximal activator of transcription in epithelial-mesenchymal transition.” J Clin Invest 117(2): 482-491.
Wei, C., J. Cheng, B. Zhou, L. Zhu, M. A. Khan, T. He, S. Zhou, J. He, X. Lu, H. Chen, D. Zhang, Y. Zhao and J. Fu (2016). “Tripartite motif containing 28 (TRIM28) promotes breast cancer metastasis by stabilizing TWIST1 protein.” Sci Rep 6: 29822.
Wulfkuhle, J. D., J. A. Aquino, V. S. Calvert, D. A. Fishman, G. Coukos, L. A. Liotta and E. F. Petricoin, 3rd (2003). “Signal pathway profiling of ovarian cancer from human tissue specimens using reverse-phase protein microarrays.” Proteomics 3(11): 2085-2090.
Xiong, H., J. Hong, W. Du, Y.-w. Lin, L.-l. Ren, Y.-c. Wang, W.-y. Su, J.-l. Wang, Y. Cui, Z.-h. Wang and J.-Y. Fang (2012). “Roles of STAT3 and ZEB1 Proteins in E-cadherin Down-regulation and Human Colorectal Cancer Epithelial-Mesenchymal Transition.” Journal of Biological Chemistry 287(8): 5819-5832.
Yu, J.-M., W. Sun, F. Hua, J. Xie, H. Lin, D.-D. Zhou and Z.-W. Hu (2015). “BCL6 induces EMT by promoting the ZEB1-mediated transcription repression of E-cadherin in breast cancer cells.” Cancer Letters 365(2): 190-200.