Skip to main content
Download PDF

Proteomic profiling identifies markers for inflammation-related tumor–fibroblast interaction

Clinical Proteomics volume 14, Article number: 33 (2017) Cite this article

Abstract

Background

Cancer associated fibroblasts are activated in the tumor microenvironment and contribute to tumor progression, angiogenesis, extracellular matrix remodeling, and inflammation.

Methods

To identify proteins characteristic for fibroblasts in colorectal cancer we used liquid chromatography-tandem mass spectrometry to derive protein abundance from whole-tissue homogenates of human colorectal cancer/normal mucosa pairs. Alterations of protein levels were determined by two-sided t test with greater than threefold difference and an FDR of < 0.05. Public available datasets were used to predict proteins of stromal origin and link protein with mRNA regulation. Immunohistochemistry confirmed the localization of selected proteins.

Results

We identified a set of 24 proteins associated with inflammation, matrix organization, TGFβ receptor signaling and angiogenesis mainly originating from the stroma. Most prominent were increased abundance of SerpinB5 in the parenchyme and latent transforming growth factor β-binding protein, thrombospondin-B2, and secreted protein acidic-and-cysteine-rich in the stroma. Extracellular matrix remodeling involved collagens type VIII, XII, XIV, and VI as well as lysyl-oxidase-2. In silico analysis of mRNA levels demonstrated altered expression in the tumor and the adjacent normal tissue as compared to mucosa of healthy individuals indicating that inflammatory activation affected the surrounding tissue. Immunohistochemistry of 26 tumor specimen confirmed upregulation of SerpinB5, thrombospondin B2 and secreted protein acidic-and-cysteine-rich.

Conclusions

This study demonstrates the feasibility of detecting tumor- and compartment-specific protein-signatures that are functionally meaningful by proteomic profiling of whole-tissue extracts together with mining of RNA expression datasets. The results provide the basis for further exploration of inflammation-related stromal markers in larger patient cohorts and experimental models.

Background

Accumulating evidence powerfully demonstrates that reactive fibroblasts are central players in physiological and pathological processes [21, 42, 50, 52]. Activation patterns consist of alterations in the cytoskeleton, production of growth factors as well as cytokines, modulation of the extracellular matrix (ECM) and a migratory phenotype. Cancer associated fibroblasts (CAF) are activated by growth factors derived from the cancer cells (e.g. transforming growth factor β—TGFβ) and develop a wound healing phenotype [58]. They are also stimulated by the inflammatory environment which is a hallmark of cancer [38]. Especially cells of the innate immune system may be pro-tumorigenic, pro-angiogenic and pro-metastatic [23, 32, 60]. This resulted in the description of cancer as a wound that does not heal [25]. In this chronic inflammatory environment, endothelial cells and fibroblasts have been shown to express cell type-specific pro-inflammatory signatures [33, 76] that can be modeled by exposure of cells to interleukin 1β (IL1β) in vitro [75]. Comparison of fibroblasts obtained from the skin, the lung and the bone marrow demonstrated that these signatures are tissue-specific and that tissue-specific characteristics are retained in cancer-associated fibroblasts [76].

In colorectal cancer (CRC) such activated fibroblasts have been identified as a driving force of tumor development [14, 41, 52] and both TGFβ and inflammation are activating factors. High expression of TGFβ is a well-established characteristic in CRC [18, 64, 83] and has been identified as a marker of poor prognosis [64]. Inflammation is already apparent in premalignant lesions derived either from chronic inflammatory bowel disease giving rise to typically inflammation-driven tumors [65] or from upregulation of cyclooxygenase-2 in colonic polyps [26, 36]. The inflammatory environment has been shown to support expansion of tumor-initiating cells by providing eicosanoids and interleukins [7, 40, 88]. These small subpopulations of tumor-initiating cells have been identified in CRC [28, 63] and were reported to be involved in therapy resistance and metastasis [24].

This paper aims to determine whether a deeper insight into the connective tissue processes in CRC and the contribution of cancer-associated fibroblasts can be gained by using tissue-proteomics. For this purpose, we have undertaken a liquid chromatography-tandem mass spectrometry (LC–MS/MS) based proteome analysis of human CRC tissue specimens and paired normal intestine. We focused on stromal, inflammation related proteins, verified selected markers using immunohistochemistry (IHC) and linked those candidates to the Consensus Molecular Subtype signatures by screening publicly available RNA expression datasets.

Methods

Tissue acquisition

The study was approved by the Ethics Commission of the Medical University of Vienna (EK 1659/2012). Patients suffering from colorectal cancer who underwent surgery at the General Hospital of Vienna were asked for their informed consent to use part of the resected tissue for analysis of tumor protein markers. Tumor tissue and normal mucosa 15 cm or as far as possible from the tumor (control tissue) were excised by the pathologist and stored at − 80 °C. For protein extraction and LC–MS/MS analysis 6 tissue pairs were used, 3 of them from stage II tumors without any sign of invasion (low-stage) and 3 from tumors that had already spread to the lymph-nodes or the peritoneum (high-stage). Each tissue was analysed individually and the abundance data pooled for statistical analysis. Formalin-fixed paraffin-embedded tissue sections were analyzed by IHC.

For immunohistochemistry formalin-fixed paraffin-embedded tissue sections were obtained from the original 6 patients as well as additional 20 tumors (stages II, III and IV) and 10 normal mucosa resection margins.

Extraction and digestion of proteins

Frozen tissue samples were incubated in sample buffer (7.5 M urea, 1.5 M thiourea, 4% CHAPS, 0.05% SDS, 100 mM dithiothreitol) for 10 min on ice. Subsequently, proteins were extracted by means of an ultrasonic stick. Protein concentrations were determined using a Bradford assay (Bio-Rad-Laboratories, Germany). Thereafter, in-solution digestion of proteins was performed with trypsin (Roche Diagnostics, Germany), as described previously [9, 74]. Briefly, 20 µg of protein were concentrated on a pre-washed 10 kDa molecular weight cut-off filter (Pall Austria Filter GmbH, Vienna, Austria). Upon reduction with dithiothreitol (5 mg/ml dissolved in 8 M guanidinium hydrochloride in 50 mM ammonium bicarbonate buffer, pH 8) and alkylation with iodoacetamide (10 mg/ml in 8 M guanidinium hydrochloride in 50 mM ammonium bicarbonate buffer), proteins were digested enzymatically overnight at 37 °C using trypsin (Roche Diagnostics, Germany). After digestion of proteins, peptide samples were cleaned up using C-18 spin columns (Pierce, Thermo Fisher Scientific), dried and stored until further LC–MS/MS analyses.

LC–MS/MS analysis

Samples were reconstituted in 5 µl 30% formic acid containing 10 fmol each of 4 synthetic standard peptides and then immediately diluted with 40 µl mobile phase A (98% H2O, 2% acetonitrile, and 0.1% formic acid). The synthetic peptides [Glu1-Fribrinopeptide B, EGVNDNEEGFFSAR; M28, TTPAVLDSDGSYFLYSK; HK0, VLETKSLYVR; HK1, VLETK(ε-AC)SLYVR] were spiked into each sample as an internal quality control for monitoring LC–MS instrument stability. Five microliters of the solution were injected into the nano HPLC-system (Dionex Ultimate 3000) loading peptides on a 2 cm × 75 µm C18 Pepmap100 pre-column (Thermo Fisher Scientific) at a flow rate of 10 µl/min using mobile phase A. Afterwards, peptides were eluted to a 50 cm × 75 µm Pepmap100 analytical column (Thermo Fisher Scientific) at a flow rate of 300 nl/min, using a gradient from 8 to 40% mobile phase B (80% acetonnitrile, 20% H2O, 0.1% formic acid) over 235 min. The nano-HPLC system was coupled to a QExactive orbitrap with a nanospray ion source (Thermo Fisher Scientific). MS scans were performed in the range from m/z 400 to 1400 at a resolution of 70,000 (at m/z = 200), MS/MS scans at a resolution of 17,500 (at m/z = 200), using a top 12 method and applying HCD fragmentation at 30% normalized collision energy.

Protein data interpretation

Identification of proteins and label-free quantification (LFQ) were performed using the MaxQuant 1.5.2.8 software including the Andromeda search engine and the Perseus statistical analysis package version 1.5.2.3 [19, 20], searching against the UniProt database for human proteins (version 102,014 with 20,195 entries, restricted to reviewed entries only). Search criteria included a peptide mass tolerance of 25 ppm, an MS/MS match tolerance of 20 ppm, a maximum of two missed cleavages and at least one unique peptide per protein. Carbamidomethylation of cysteines was set as fixed modification, whereas methionine oxidation as well as N-terminal protein acetylation as variable modifications. Furthermore, match between runs was performed using a 5 min match time window. For peptides and proteins, a false discovery rate (FDR) of less than 0.01 was applied. Protein regulation was determined by comparing the LFQ values for each individual protein in the different samples using Perseus, normalizing to the same initial protein amount of 20 µg. Regarding the protein inference problem, indistinguishable proteins sharing the same peptides were summed up into protein groups. For quantitative analysis, a protein had to be found in at least 4 of 6 replicates in a group and imputation of missing values were performed based on normal distribution. Changes in protein abundance values were determined by a two-sided t test, considering proteins as significantly altered when the abundance difference was at least threefold with p < 0.05. Additionally, to highlight the most robust changes in protein abundance, we performed a permutation-based FDR correction applying a global FDR < 0.05. Proteins meeting this additional criterion are marked with a “+” for multiparameter (MP) significance in the tables und used for further study.

Additionally, raw files were analyzed using Proteome Discoverer 1.4 (Thermo Fisher Scientific, Austria) utilizing Mascot 2.5 (Matrix Science, UK) in order to enable upload of mass spectrometric data to a publicly available repository. Therefore, protein identification was performed by searching against the SwissProt Database (version 11/2015 with 20.193 entries) with mass tolerance at the MS1 level of 50 ppm and 100 mmu at the MS2 level allowing for up to two missed cleavages per peptide. Peptide modifications set were carbamido-methylation for cysteines as fixed modification as wells as methionine oxidation and protein N-terminal acetylation. Data was submitted to the ProteomeXchange Consortium via the PRIDE partner repository and can be accessed via http://www.proteomeexchange.org with the identifier PXD006776 [85].

Biological context of protein subset

To put the 24 protein subset into biological context, we analysed GO-term enrichment of biological processes using the Cytoscape plugin ClueGO in combination with CluePedia [10, 11]. Detected proteins were used as universe, p values were corrected according to Benjamini–Hochberg.

Immunohistochemical analysis

Tumors were fixed, dehydrated and embedded in paraffin. After sectioning, slides were dewaxed by incubation at 65 °C for 10 min and subsequent xylol treatment. Rehydration was performed in a gradient of ethanol and ddH2O, then endogenous peroxidases were inactivated with 0.3% (v/v) H2O2 in 1× PBS and antigens retrieved in a steamer with 10 mM citrate-buffer pH 6. Slides were washed in 1× PBS containing 0.1% Tween 20 (PBST), blocked for 5 min at room temperature (RT) with Ultra V Block (Thermo Fisher Scientific) and incubated with primary antibodies for 30 min at RT. If not otherwise stated, primary antibodies were diluted in PBST containing 1% goat serum (Dako) as follows: THBS2 (TA590658, Origene) 1:150, αSMA (1A4, Dako) 1:100, SerpinB5 (sc-271,694, Santa Cruz) 1:100, SPARC (D10F10, Cell Signaling) 1:2400 in SignalStain Antibody Diluent (Cell Signaling). Bound antibodies were detected by incubation for 10 min with Primary Antibody Enhancer (Thermo Fisher Scientific), followed by 15 min with HRP Polymer (Thermo Fisher Scientific), visualized by incubation for 2 min with DAB substrate (Dako), and counterstained with hematoxylin solution.

For analysis, slides were digitalized using a microscopic slidescanner (Pannoramic Midi, 3DHistech) with a 40× objective and 5 randomly chosen areas analyzed using Tissue Studio® (Definiens®) histomorphometric software (Additional file 1: Figure S1). Marker intensity thresholds were set to 0.25/0.5/0.75 for low/mid/high DAB detection. The percentage of stained area compared to total area for both epithelial and stromal compartments of tumor/normal tissue sections were statistically analyzed with GraphPad Prism 6 by Mann–Whitney U test and considered significant when p < 0.5.

Computational analysis of differential expression

For comparing gene expression between the subtypes and healthy colon mucosa, published data sets were obtained from ArrayExpress as used in the consensus subtype classification of CRC [34]. Additionally, the data set GSE44076 [70], also obtained from ArrayExpress, was analyzed independently.

All analysis were performed with the statistics program R. The datasets were pre-processed and normalized as follows: Data from the consensus subtype classification was normalized using fRMA. Comparison was made between the consensus molecular subtypes (CMS2-CMS1, CMS3-CMS1… etc.). For GSE44076, normalization was performed using RMA (as no fRMA vectors were available for the hgu219 platform) followed by comparison of data between mucosa, normal and tumor.

For all datasets, limma package (Bioconductor.org) was used in R to perform linear modelling and create contrasts for pair-wise comparisons. After model fitting, empirical Bayesian statistics was used for analysis of differential expression [77]. FDR detection was performed with the Benjamini–Hochberg method.

Results

Protein abundance in tumor tissue

Surgical specimens of both tumor and normal mucosa were obtained from 6 colon cancer patients. With regard to tumor stage, there were 3 localised tumors (stage II, low-stage) and 3 invasive tumors (stage III and IV, high-stage). To avoid artefacts caused by cell isolation procedures whole tissue specimens were homogenized and extracted using sample buffer to yield total protein samples for proteome analysis by LC–MS/MS. We identified a total of 4864 proteins in the analyzed samples (Fig. 1a; Additional file 2: Table S1).

Fig. 1
figure 1

Proteom analysis of CRC. a Protein abundances in tumor versus normal tissue samples. Volcano plot representing differences in LFQ values (log2FC, fold-changes in logarithmic scale to the base of two) of proteins including corresponding p values (logarithmic scale). Red labeled proteins are higher abundant in the tumor tissue whereas green labeled proteins are higher abundant in the normal tissue samples. b Stromal contribution. Proteomics derived data of significantly regulated proteins were plotted against stromal contribution of their respective mRNA [41]. Red indicates upregulated proteins of stromal origin. Blue indicates the 2 stromal proteins THBS2 and SPARC plus the parenchymal SerpinB5 that were chosen for IHC staining

Full size image

Statistical analysis of the pooled results revealed 122 proteins that were increased in abundance and 81 proteins that were significantly decreased in the tumor tissues compared to normal mucosa. When the comparison was made separately for low-stage versus normal and high-stage versus normal, alterations were similar for the majority of proteins with very few exceptions. Relevant for our analysis were fibroblast activation protein (FAP) and the macrophage marker MRC1 that only achieved multi-parameter significance in low-stage tumors versus normal tissue. Both proteins were included in spite of this restriction, because they are classical cell type markers. In addition, the mast cell markers TPSAB1, CPA3, and CMA1 were significantly different in high-stage tumors compared to low-stage (Additional file 2: Table S1).

To determine whether proteins with altered abundance originated from the tumor epithelium or from the connective tissue, we turned to the analysis of gene expression data published by Isella et al. [41], who identified stromal contribution (SC) to CRC xenografts by their host origin. Alignment with the SC calculated by Isella et al. was possible for 187 out of 203 proteins identified in our analysis (Additional file 3: Table S2). Of those proteins 115 had elevated levels in the tumor and 73 had higher abundance in the normal tissue. The alignment revealed that the majority of proteins that were elevated in the tumor (104/115) were mostly of parenchymal origin (SC < 50%) (Fig. 1b). This included SerpinB5, a protein induced by TGFβ signaling and/or cell stress [12] as the strongest increased protein. Only 11 of the high-abundance proteins were of stromal origin (highlighted in Fig. 1b; Table 1). Prominent among those proteins was latent-transforming growth factor β-binding protein 2 (LTBP2)—an indicator of TGFβ signaling that has been previously identified as a product of CRC-associated fibroblasts obtained from a spontaneous mouse colon cancer model [82]. Another protein whose level was elevated in the tumor was insulin-like growth factor-binding protein 7 (IGFBP7) that has been recently described as a tumor stroma and epithelial-to-mesenchymal transition (EMT) marker in various epithelial cancers [67]. The fibroblast activation protein (FAP) was increased more than fourfold. The elevated stromal proteins thrombospondin-2 (THBS2) and secreted protein acidic and cysteine-rich (SPARC) have been previously identified as parts of a fibroblast-specific inflammation signature [75]. In addition, collagens and lysyloxidase homolog 2 (COL8A1, COL12A1, LOXL2) were found elevated as well as adipocyte enhancer-binding protein 1 (AEBP1) that is involved in wound healing and the endothelial activation marker endoglin (ENG).

Table 1 Proteins of mainly stromal origin that are increased in abundance
Full size table

Forty-five of the 73 proteins that were decreased in abundance in the tumor came mostly from the stroma (SC > 50% according to [41]) including ECM constituents and mast cell markers. In addition, ECM proteins of epithelial origin were also altered in a tumor-specific manner (Table 2). Among the collagen I associated FACIT collagens [27], the type XII α1-chain was increased. The type XIV α1-chain was decreased, but narrowly missed the threefold change threshold (not shown). Among the widespread network collagens [27], the type VI α5-chain was decreased. At the same time type VIII α1-chain was increased and the basement membrane component laminin B shifted from subunit B3 to subunit B2. In addition, LOXL2 was increased. Matrix metalloproteinases (MMP) were not strongly changed: for MMP 1, 8, and 9 differences missed significance (p > 0.05) and MMP 2 and 14 were increased but missed the threshold (not shown). The small leucine-repeat proteoglycans decorin (DCN) and asporin (ASPN) were strongly decreased.

Table 2 Alterations in ECM constituents
Full size table

The mast cell markers TPSAB1, CPA3 and CMA1 were decreased in abundance especially in the low-stage tumors together with reduced levels of the macrophage marker MRC1 (Table 3). In high-stage tumors the abundance of these proteins increased and was significantly higher than in low-stage, but not higher than in normal tissue.

Table 3 Macrophage and mast cell markers
Full size table

In summary, this defines a set of 24 proteins, 20 of those originating from the stroma, playing important roles in ECM-organization, angiogenesis, TGFβ signaling and inflammation (Fig. 2). Gene ontology (GO)-term enrichment analysis based on biological processes grouped most of the 24 proteins into ECM-organization (GO ID: GO:0030198, 14 proteins, p < 0.0001,), angiogenesis (GO ID: GO:0001525, 8 proteins, p < 0.0001) and TGFβ receptor signaling pathways (GO ID: GO:0007179, 3 proteins, p = 0.0017) (Additional file 4: Table S3). Due to the mostly indirect role of inflammation, the specific GO-term “inflammatory response” (GO ID: GO:0006954) and child terms were not covered, instead processes connected to inflammation like “response to corticosteroid” (GO ID: GO:0031960, 3 proteins, p < 0.0001), “blood vessel development” (GO ID: GO:0001568, 8 proteins, p < 0.0001) or “endothelial cell migration” (GO ID: GO:0043542, 4 proteins, p < 0.0001) were included.

Fig. 2
figure 2

Proteomic alterations associated with inflammation, ECM organization, TGFβ receptor signaling pathway and angiogenesis. Schematic presentation of identified regulated proteins and their functional annotation according to literature. Most proteins were confirmed and grouped by GO-term enrichment analysis of biological processes using the Cytoscape plugin Cluego. Significant associations were found with ECM-organization (14 proteins, p > 0.0001), angiogenesis (8 proteins, p < 0.0001) and TGFβ receptor signaling pathways (3 proteins, p = 0.0017). Inflammation was indirectly covered by terms like “response to corticosteroid” (3 proteins, p < 0.0001), “blood vessel development” (8 proteins, p < 0.0001) and “endothelial cell migration” (4 proteins, p < 0.0001). Blue represent proteins with lower abundance in the tumor, while red indicates upregulation compared to adjacent normal mucosa

Full size image

Transcriptome profile in published data sets

To assess whether the alterations we observed were similar at the transcriptional level, we turned to transcriptome datasets available from Sanz-Pamplona et al. [70] who analyzed gene expression of 98 stage II CRC tumors and their corresponding normal mucosa as well as of mucosa from 50 healthy individuals. We compared the published RNA levels of parenchymal SerpinB5, as well as 6 upregulated stromal proteins (THBS2, SPARC, FAP, COL12A1, COL8A1, LTBP2, Fig. 3a–f). Significant induction of RNA levels was found for all 7 markers in the tumor samples as compared with the normal mucosa indicating that the upregulated proteins found by LC–MS/MS originated from cells of the tumor and its microenvironment (Fig. 3a–g). Interestingly, increased mRNA expression was already observed in normal tissue of tumor patients compared to healthy donors for THBS2, SPARC, LTBP2 (Fig. 3a, b, f) as well as the collagens COL8A1 and COL12A1 (Fig. 3d, e). For the mast cell marker CMA1 and the macrophage marker MRC1, mRNA-levels in the normal mucosa of CRC patients were elevated compared to both healthy and tumor tissue. In the tumor tissue, CMA1 mRNA was low throughout, while MRC1 mRNA levels showed high inter-individual variations (Fig. 3h, i). Overall, mRNA level changes correlated well with protein abundance alterations for 19 of the 24 proteins in our marker set (Additional file 1: Figures S2, S3). Exceptions were the laminin chains, collagen XXVIII, ASPN and AEBP1.

Fig. 3
figure 3

mRNA expression of tissues derived from healthy volunteers and CRC patients [70]. For analyzing gene expression of stage II CRC tumors, we obtained the dataset GSE44076 consisting of tissue expression data of 50 healthy volunteers (healthy) as well as 98 CRC samples (tumor) and paired normal adjacent mucosa (normal). a–f upregulated stromal markers; g: SerpinB5. Overall, expression increased from healthy to normal and again from normal to tumor tissue; h, i for mast cell and macrophage markers, expression was highest in the normal tissue compared to both healthy and tumor tissue, with the tumor showing the lowest expression. Statistical analysis for differential gene expression was perfomed by using empirical Bayesian statistics with FDR detection according to the Benjamini–Hochberg method. *p < 0.05; **p < 0.01; ***p < 0.001; ns not significant

Full size image

RNA expression data from the CRC consensus classification [34] were used as a second source for gene expression profiles. The connective tissue markers we found increased at protein level (Table 1) were also increased on the RNA level in CSM4 tumors that are characterized by a high connective tissue content when compared to any other subtype. Figure 3 shows the regulation pattern for FAP, THBS2, SPARC, COL8A1, and COL12A1 (Fig. 4a–e). SerpinB5, originating from the parenchyma, was not enriched in CMS4 tumors (Fig. 4f).

Fig. 4
figure 4

Consensus Molecular Subtype Classification of upregulated stromal markers and SerpinB5 [34]. For comparing gene expression between the subtypes and selected markers, raw data of published data sets were obtained as used in the consensus subtype classification of CRC [34], processed, analyzed and differences plotted as log2-foldchange (log2-FC). a–e expression of 5 selected stromal markers that were upregulated on protein level in our data set. f results for SerpinB5, a parenchymal derived marker. Statistical analysis of log2-foldchanges were performed by using empirical Bayesian statistics. Adjusted p values were < 0.001 except for those marked non-significant (ns)

Full size image

Tissue localization of proteins

From the list of proteins increased in the tumor, 3 were chosen for verification of tissue localization. Those were SerpinB5, because it was the strongest elevated marker and should localize to the parenchyma; SPARC and THBS2, because they have previously been described as prominent parts of a fibroblast-specific inflammation signature [75] and their SC was > 70%. Tissue sections of 26 CRC tissues (tumor stage II–IV) and 16 normal mucosa specimen were stained by IHC and quantified separately for the epithelial and stromal compartment (Additional file 1: Figure S1).

SerpinB5 staining was very weak in normal mucosa (Fig. 5a). In the tumors, staining was strong and almost exclusively localized in epithelial cells (Fig. 5b). Quantification of staining intensity showed a significant increase in the tumor both when total tissue was assessed and when only the epithelial cells were scored (Fig. 5c, d). Specifically, all tumors harbored SerpinB5-positiv cells representing 61.3 ± 20.6% of the parenchyme. This represents a highly significant increase as compared to normal tissue (10.5 ± 9.2%).

Fig. 5
figure 5

Tissue localization of selected proteins. Serial sections of normal mucosa and tumor tissue were stained using antibodies directed against SerpinB5 (a, b), THBS2 (e, f), SPARC (i, j), and αSMA (m, n). a, e, i, m representative area of normal intestinal mucosa. b, f, j, n representative area of tumor tissue. Scale bars correspond to 500 µm. A selected area was magnified ×20 and is shown in the inserts. Staining intensity for all images was quantified using Definiens software. The diagrams depict the pooled quantification obtained from 26 tumor samples and corresponding normal tissue with regard to overall tumor and normal tissue (c) or tumor and normal epithelium for SerpinB5 (d). For THBS2, quantification results are depicted for total tumor and normal areas (g) or tumor and normal stromal compartments (h). For SPARC, quantification is presented for tumor and normal stromal compartments (k) and specifically for low/high staining intensities (l). *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns: not significant, according to Mann–Whitney U test

Full size image

In contrast, staining for THBS2 was seen in both the tumor and the stroma. In normal colon tissue, it was highest in the blood vessel walls and the muscularis mucosae. Otherwise staining was evenly distributed between the stromal and epithelial compartment (Fig. 5e). In the 6 tumors used in the proteome analysis, IHC showed increased staining in the stroma (Fig. 5f; Additional file 1: Figure S4). In the larger patient cohort used for IHC validation we found higher inter-specimen variability than for SerpinB5 or SPARC. In some of the specimen, staining intensity increased in the epithelium as well as in the stroma. This is also reflected by the quantification results that produced significant increase in staining intensity for both total tumor and stroma, but also large standard deviations (Fig. 5g, h).

Overall SPARC staining intensities were weak. In the normal colon, signals came mostly from the submucosal blood vessel endothelial cells. Both the epithelial cells and the stromal cells of the mucosa were almost all negative (Fig. 5i). In the tumor stroma, staining increased compared to the normal stroma (Fig. 5j) but the difference did not reach significance (Fig. 5k). Higher magnification revealed that staining was strongest in individual cells of the tumor microenvironment (Fig. 5j insert). This population was observed in all 26 tumor specimen and accounted for 10.66 ± 5.15% of the stained area in the tumors, while it was almost absent (0.46 ± 0.64%) in normal tissue (Fig. 5l).

αSMA was stained as a well-established fibroblast activation marker and was positive in proximity to the epithelium and around blood vessels in normal tissue (Fig. 5m). In the tumor it was clearly elevated throughout the stromal compartment (Fig. 5n) similar to the staining pattern of THSB2 (Fig. 5f).

Discussion

The microenvironment of CRC is characterized by TGFβ as well as inflammation that both contribute to CAF activation [42]. Our results now show in CRC a set of 24 proteins (Fig. 2) of mainly stromal origin that are associated with TGFβ, inflammation, matrix remodeling and wound healing as confirmed by GO-term enrichment analysis. Of these proteins, 20 came from the tumor stroma and 4 originated from the epithelial compartment. The epithelial proteins were included as markers of cell stress (SerpinB5) or indicators of ECM alterations (laminins, collagen type XXVIII). Of the stromal proteins, 11 were elevated in tumor tissue and were part of a fibroblast signature correlating with poor prognosis identified by Calon et al. [14]. Nine stromal proteins decreased in abundance were ECM constituents involved in matrix remodeling and markers of macrophages and mast cells.

Methodically, our analysis used whole tissue extracts that reflected all cell populations in the tissue. The origin of altered proteins was determined by linking our data with known expression profiles obtained from pure cell populations [75, 76, 82] and the stromal contribution map published by Isella et al. [41]. As the initial purpose of the study was to determine the feasibility of this strategy, we only used very few patient samples and no stage-specific analysis was intended. Proteins of interest were validated by comparing them to published mRNA expression datasets [34, 70] and by confirming the predicted tissue localization of SerpinB5, THBS2 and SPARC in a larger set of 26 tumors and 16 normal tissue samples.

For SerpinB5 and the elevated stromal proteins, protein/mRNA levels correlated extremely well. This was surprising, because protein abundance is determined not only by transcription but also by post-transcriptional regulation, translational control mechanisms and differences in protein stability. In a comparative proteo-genomic study of the colon and rectum the average Spearman’s correlation coefficient (rs) between protein and RNA regulation was only 0.23 [93]. Looking at specific proteins, correlation was above average for THBS2 and SPARC (rs = 0.357 and 0.4169 respectively, p < 0.001). The highest correlation was found for SerpinB5 (rs = 0.8142, p < 0.001).

In our proteome analysis the protein with the highest level in colorectal tumors as compared to normal intestinal tissue was SerpinB5. Staining of tissue sections showed the protein was nearly absent in the normal mucosa and localized almost exclusively in the tumor epithelium. In individual tumors, SerpinB5-positive cell populations ranged from 14 to 96% of all epithelial cells. While SerpinB5 was characterized as a tumor suppressor with anti-invasive and anti-angiogenic functions in several tumor types [8, 12], analysis of colorectal tumors indicates increased protein levels in the tumor and an association with more aggressive disease and worse prognosis [46, 84]. This feature is also correlated with nuclear localization of the protein at sites of tumor budding [46]. On the other hand SerpinB5 is already observed in premalignant lesions—specifically serrated polyps [66] and in inflammatory bowel disease [15]. The gene may be upregulated by TGFβ [12] which is a characteristic feature of most CRCs [18] or it may be induced by the cellular stress observed in both CRC and active inflammation of IBD lesions [46, 84]. In our dataset evidence of inflammation comes from the increased levels of THBS2 and SPARC in the tumor stroma that are both induced by IL1β in fibroblasts in vitro [75].

The largest increase of a stromal protein was observed for LTBP2, a large ECM protein associated with elastin that binds latent TGFβ and therefore modulates TGFβ-signaling [29, 56]. Its expression is induced by TGFβ [3], similar to SerpinB5, which can be regarded as an indicator of a TGFβ-activated stroma. More recently it has been identified as a marker of poor prognosis in cervical [62], head-and-neck [37] and pancreatic cancer [87]. It was identified as part of a CAF-specific signature in fibroblasts isolated from a spontaneous mouse colon cancer model [82].

THBS2 is a multifunctional extracellular glycoprotein that is produced by most cell types. It has functions in inflammation, inhibits angiogenesis and mediates ECM assembly [1, 13]. Up-regulated THSB2 has been described to inhibit tissue repair due to aberrant fibroblast migration and adhesion [5]. On one hand, THBS2 is mostly described as an anti-angiogenic, anti-metastatic factor in cancer [13, 79]. On the other hand, it is involved in epithelial-to-mesenchymal transition in breast cancer by enhancing AXL-dependent activation of niche fibroblasts by the fibroblastoid tumor cells [22]. In CRC THSB2 mRNA levels were reported to be increased as compared to normal tissue: analysis of published gene expression data from the Cancer Genome Atlas indicates that expression increases with tumor stage and node involvement [90]. Our LC–MS/MS analysis found increased abundance of the protein in tumors. The IHC analysis demonstrated that the protein is localized in both the tumor parenchyma and the stroma as predicted from the stromal contribution map [41]. IHC revealed higher intensity staining in the tumors, but also high variability between individual tumors. Whether these differences are pathophysiologically and clinically relevant, needs to be investigated in a follow-up study with a larger patient cohort. In addition, investigation of microarray datasets derived from 98 colorectal patients and 50 healthy volunteers revealed that THBS2 mRNA is increased not only in the stage II tumors but also in normal adjacent mucosa of cancer patients compared to healthy control individuals.

SPARC is a secreted factor that is produced in organs undergoing rapid proliferation or remodeling [47]. It is essential for wound healing [6] and up-regulated in endothelial cells by exposure to VEGF [44]. On the molecular level, its main effect is modulation of the ECM and cell adhesion [57]. In tumors, its role seems to be context-dependent with tumor suppressor characteristics in urothelial [68] and pancreatic cancer [71], anti-metastatic impact in prostate cancer [72] and pro-tumorigenic effects in breast cancer and glioma cell models [30, 53]. In melanoma, SPARC has been shown to enhance tumorigenesis [48, 49] by inhibiting cytotoxic anti-tumor response [4]. In murine breast cancer cell models, the protein was reported to induce an immunosuppressive environment and to interact with myeloid suppressor cells to induce EMT [69].

The cellular origin of SPARC protein also seems important for its function, as specifically stromal SPARC is associated with tumor progression in pancreatic cancer [35], while tumor-derived SPARC increased vascular permeability in melanoma [81]. In CRC, stromal SPARC was found to be increased as compared to normal colon and increasing from stage I to stage IV. It also was a predictive marker for good overall survival [17]. In our specimens SPARC staining was strongest in individual cells in the tumor-microenvironment that represented about 10% of the cells in the tumor stroma. Co-localization of αSMA with SPARC and the spindeloid shape of positive cells points at activated fibroblasts as protein source. However, only a subset of all αSMA-positive cells were also expressing SPARC, indicating that only a fraction of CAFs produced SPARC in our specimens.

Our analysis of public available datasets revealed that SPARC, THBS2, FAP, and LTBP2 were closely related to CMS4, a subtype which is marked by increased CAF abundance and an inflammatory as well as immunotolerant microenvironment [34]. For 3 of the 4 genes (SPARC, THBS2, LTBP2) mRNA was increased not only in the tumors but also in normal adjacent mucosa compared to healthy mucosa. In addition, CMA1 and MRC1 were upregulated in patient normal mucosa compared to both healthy and tumor tissue. Taken together with upregulation of THBS2 and SPARC by IL1β, this suggests that inflammatory activation of stromal fibroblasts is not only an integral part of CRC, but even spreads into the surrounding normal mucosa.

The decreased levels of MRC1 and the mast cell markers in the low-stage tumors were a surprising result, because tumor-associated mast cells are reported to be increased in CRC [80, 92]. They are also considered indicative of advanced, aggressive tumors and poor prognosis [2, 16, 31, 54]. Here, analysis of mRNA data sets suggests that macrophages and mast cells accumulate in the normal mucosa of tumor patients as compared to healthy individuals, thus increasing abundance in the normal adjacent mucosa. We also found MRC1 and CMA1 significantly increased in high-stage versus low-stage tumors, which is in agreement with reports on macrophage and mast cell invasion of CRC mentioned above. Due to small sample size and high inter-individual differences, further investigations with a larger cohort are needed for a more precise and accurate analysis of these markers.

With regard to the ECM, we did not observe any alterations in the highly abundant collagens type I and IV. In our analysis, protein abundance was increased for collagen I-associated FACIT collagen type XII that has been described as a product of CRC-associated CAFs at the invasion front [43] and the network collagen Type VIII, which is inductive of angiogenesis [39]. By contrast, collagen type XIV was reduced in abundance as well as collagen type VI that belongs to the widespread network collagens [27]. Alteration of collagen type VI affected the α1-, α2-, and α3-chains (abundance < 0.5-fold as compared to control; Additional file 2: Table S1) as well as the α5-chain that we found decreased more than threefold. This observation was surprising, because collagen type VI α3 had been previously described as a marker upregulated in CRC and associated with advanced tumors and poor prognosis [61]. In addition, we observed high LOXL2 that are markers of a stiff cancer-associated matrix that affects cell adhesion, cell migration, cancer stem cell characteristics and EMT [59, 89, 91]. The upregulated proteolytic enzyme FAP and SPARC are also signs of a tumor-specific matrix [45]. MMPs that have previously been found upregulated in CRC and are considered to be prognostic factors [55, 86] were only moderately altered in our proteome analysis. The reason may be that these secreted proteins do not remain localized in the tumor, but may spread to the adjacent normal tissue and can even be found increased in patient serum [78, 94]. In many tumors the cancer-specific matrix also includes increased ASPN and DCN, due to their induction by e.g. TGFβ or stiff matrix. This was not observed in our tumor tissues, which can be regarded as an additional marker of inflammation, as production of ASPN and DCN is suppressed by IL1β [73] and low abundance of these markers was also observed in triple-negative inflammatory breast cancer CAFs [51].

Taken together the alterations we describe are characteristic for ECM-remodeling CAFs as defined by Kalluri [42] and correlates well with the CAF signature obtained from explant cultures by Torres et al. [82]. However, we did not find any growth factors and chemokines that would drive either angiogenesis or tumor cell survival and metastasis. This may be due to limitations of tissue proteomics with regard to low abundant proteins. In conclusion, we describe a protein signature reflecting inflammatory activation of CRC-associated CAFs and activation by TGFβ. The marker set we found in the colon share THBS2 and SPARC with dermal fibroblasts [75] and FAP with bone marrow fibroblasts [76]. It is distinctly different from signatures found in melanoma and lung cancer associated fibroblasts [75, 76] and from the wound healing signature observed in breast cancer [33].

Conclusions

This study demonstrates the feasibility of detecting tumor- and compartment-specific protein-signatures that are functionally meaningful by proteomic profiling of whole-tissue extracts together with mining of RNA expression datasets. In the context of our study of CRC-stroma interactions, the results provide the basis for further exploration of inflammation-related stromal markers in larger patient cohorts and experimental models.

At a wider scope, the same methods can be used to define other marker sets with different functional annotations for CRC, but also for other malignancies. This can be developed into an efficient method of tumor classification on the protein level.

Abbreviations

AEBP1:

adipocyte enhancer-binding protein 1

ASPN:

asporin

CAF:

cancer associated fibroblast

CMA:

chymase 1

COL:

collagen

CPA3:

carboxypeptidase 3

CRC:

colorectal cancer

DCN:

decorin

ECM:

extracellular matrix

EMT:

epithelial-to-mesenchymal transition

ENG:

endoglin

FAP:

fibroblast activation protein

FDR:

false discovery rate

HPLC:

high-performance liquid chromatography

IGFBP7:

insulin-like growth factor binding protein 7

IHC:

immunohistochemistry

IL1β:

interleukin 1 beta

LC–MS/MS:

liquid chromatography tandem mass spectrometry

LOXL2:

lysyloxidase 2 homologue

LTBP2:

latent transforming growth factor beta binding protein 2

MMP:

matrix metalloprotease

MP:

multiparameter

MRC1:

mannose receptor C-type 1

rs :

correlation coefficient rho

SC:

stromal contribution

SPARC:

secreted protein acidic and cysteine-rich

TGFβ:

transforming growth factor beta

THBS2:

thrombospondin 2

TPSAB1:

tryptase alpha/beta 1

References

  1. Adams JC, Lawler J. The thrombospondins. Int J Biochem Cell Biol. 2004;36:961–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Afik R, Zigmond E, Vugman M, Klepfish M, Shimshoni E, Pasmanik-Chor M, et al. Tumor macrophages are pivotal constructors of tumor collagenous matrix. J Exp Med. 2016;213:2315–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Ahmed W, Kucich U, Abrams W, Bashir M, Rosenbloom J, Segade F, et al. Signaling pathway by which TGF-beta1 increases expression of latent TGF-beta binding protein-2 at the transcriptional level. Connect Tissue Res. 1998;37:263–76.

    Article  CAS  PubMed  Google Scholar 

  4. Alvarez MJ, Prada F, Salvatierra E, Bravo AI, Lutzky VP, Carbone C, et al. Secreted protein acidic and rich in cysteine produced by human melanoma cells modulates polymorphonuclear leukocyte recruitment and antitumor cytotoxic capacity. Can Res. 2005;65:5123–32.

    Article  CAS  Google Scholar 

  5. Bancroft T, Bouaouina M, Roberts S, Lee M, Calderwood DA, Schwartz M, et al. Up-regulation of thrombospondin-2 in Akt1-null mice contributes to compromised tissue repair due to abnormalities in fibroblast function. J Biol Chem. 2015;290:409–22.

    Article  CAS  PubMed  Google Scholar 

  6. Basu A, Kligman LH, Samulewicz SJ, Howe CC. Impaired wound healing in mice deficient in a matricellular protein SPARC (osteonectin, BM-40). BMC Cell Biol. 2001;2:15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Bellamkonda K, Chandrashekar NK, Osman J, Selvanesan BC, Savari S, Sjölander A. The eicosanoids leukotriene D4 and prostaglandin E2 promote the tumorigenicity of colon cancer-initiating cells in a xenograft mouse model. BMC Cancer. 2016;16:425.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Berardi R, Morgese F, Onofri A, Mazzanti P, Pistelli M, Ballatore Z, et al. Role of maspin in cancer. Clin Transl Med. 2013;2:8.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Bileck A, Kreutz D, Muqaku B, Slany A, Gerner C. Comprehensive assessment of proteins regulated by dexamethasone reveals novel effects in primary human peripheral blood mononuclear cells. J Proteome Res. 2014;13:5989–6000.

    Article  CAS  PubMed  Google Scholar 

  10. Bindea G, Galon J, Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics. 2013;29:661–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-into decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25:1091–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Bodenstine TM, Seftor RE, Khalkhali-Ellis Z, Seftor EA, Pemberton PA, Hendrix MJ. Maspin: molecular mechanisms and therapeutic implications. Cancer Metastasis Rev. 2012;31:529–51.

    Article  CAS  PubMed  Google Scholar 

  13. Bornstein P. Thrombospondins function as regulators of angiogenesis. J Cell Commun Signal. 2009;3:189–200.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Calon A, Lonardo E, Berenguer-Llergo A, Espinet E, Hernando-Momblona X, Iglesias M, et al. Stromal gene expression defines poor-prognosis subtypes in colorectal cancer. Nat Genet. 2015;47:320–9.

    Article  CAS  PubMed  Google Scholar 

  15. Cao D, Wilentz RE, Abbruzzese JL, Ho L, Maitra A. Aberrant expression of maspin in idiopathic inflammatory bowel disease is associated with disease activity and neoplastic transformation. Int J Gastrointest Cancer. 2005;36:39–46.

    Article  CAS  PubMed  Google Scholar 

  16. Chen X, Churchill MJ, Nagar KK, Tailor YH, Chu T, Rush BS, et al. IL-17 producing mast cells promote the expansion of myeloid-derived suppressor cells in a mouse allergy model of colorectal cancer. Oncotarget. 2015;6:32966.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Chew A, Salama P, Robbshaw A, Klopcic B, Zeps N, Platell C, et al. SPARC, FOXP3, CD8 and CD45 correlation with disease recurrence and long-term disease-free survival in colorectal cancer. PLoS ONE. 2011;6:e22047.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Coffey RJ, Goustin AS, Soderquist AM, Shipley GD, Wolfshohl J, Carpenter G, et al. Transforming growth factor α and β expression in human colon Cancer lines: implications for an autocrine model. Can Res. 1987;47:4590–4.

    CAS  Google Scholar 

  19. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26:1367–72.

    Article  CAS  PubMed  Google Scholar 

  20. Cox J, Mann M. 1D and 2D annotation enrichment: a statistical method integrating quantitative proteomics with complementary high-throughput data. BMC Bioinform. 2012;13(Suppl 16):S12.

    Article  CAS  Google Scholar 

  21. De Wever O, Demetter P, Mareel M, Bracke M. Stromal myofibroblasts are drivers of invasive cancer growth. Int J Cancer. 2008;123:2229–38.

    Article  PubMed  Google Scholar 

  22. del Pozo Martin Y, Park D, Ramachandran A, Ombrato L, Calvo F, Chakravarty P, et al. Mesenchymal cancer cell-stroma crosstalk promotes niche activation, epithelial reversion, and metastatic colonization. Cell Rep. 2015;13:2456–69.

    Article  PubMed  PubMed Central  Google Scholar 

  23. DeNardo DG. Andreu P and Coussens LM Interactions between lymphocytes and myeloid cells regulate pro- versus anti-tumor immunity. Cancer Metastasis Rev. 2010;29:309–16.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Dieter SM, Ball CR, Hoffmann CM, Nowrouzi A, Herbst F, Zavidij O, et al. Distinct types of tumor-initiating cells form human colon cancer tumors and metastases. Cell Stem Cell. 2011;9:357–65.

    Article  CAS  PubMed  Google Scholar 

  25. Dvorak HF. Tumors wounds that do not heal—redux. Cancer Immunol Res. 2015;3:1–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Eberhart CE, Coffey RJ, Radhika A, Giardiello FM, Ferrenbach S, DuBois RN. Up-regulation of cyclooxygenase 2 gene expression in human colorectal adenomas and adenocarcinomas. Gastroenterology. 1994;107:1183–8.

    Article  CAS  PubMed  Google Scholar 

  27. Fang M, Yuan J, Peng C, Li Y. Collagen as a double-edged sword in tumor progression. Tumor Biol. 2014;35:2871–82.

    Article  CAS  Google Scholar 

  28. Gemei M, Mirabelli P, Di Noto R, Corbo C, Iaccarino A, Zamboli A, et al. CD66c is a novel marker for colorectal cancer stem cell isolation, and its silencing halts tumor growth in vivo. Cancer. 2013;119:729–38.

    Article  CAS  PubMed  Google Scholar 

  29. Gibson MA, Hatzinikolas G, Davis EC, Baker E, Sutherland GR, Mecham RP. Bovine latent transforming growth factor beta 1-binding protein 2: molecular cloning, identification of tissue isoforms, and immunolocalization to elastin-associated microfibrils. Mol Cell Biol. 1995;15:6932–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Gilles C, Bassuk JA, Pulyaeva H, Sage EH, Foidart JM, Thompson EW. SPARC/osteonectin induces matrix metalloproteinase 2 activation in human breast cancer cell lines. Cancer Res. 1998;58:5529–36.

    CAS  PubMed  Google Scholar 

  31. Giusca SE, Caruntu ID, Cimpean AM, Avadanei RE, Balica AR, Jitariu AA, et al. Tryptase-positive and CD117 positive mast cells correlate with survival in patients with liver metastasis. Anticancer Res. 2015;35:5325–31.

    CAS  PubMed  Google Scholar 

  32. Grivennikov SI, Greten FR, Karin M. Immunity, inflammation, and cancer. Cell. 2010;140:883–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Groessl M, Slany A, Bileck A, Gloessmann K, Kreutz D, Jaeger W, et al. Proteome profiling of breast cancer biopsies reveals a wound healing signature of cancer-associated fibroblasts. J Proteome Res. 2014;13:4773–82.

    Article  CAS  PubMed  Google Scholar 

  34. Guinney J, Dienstmann R, Wang X, de Reynies A, Schlicker A, Soneson C, et al. The consensus molecular subtypes of colorectal cancer. Nat Med. 2015;21:1350–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Gundewar C, Sasor A, Hilmersson KS, Andersson R, Ansari D. The role of SPARC expression in pancreatic cancer progression and patient survival. Scand J Gastroenterol. 2015;50:1170–4.

    Article  CAS  PubMed  Google Scholar 

  36. Gupta RA, DuBois RN. Cyclooxygenase-2 inhibitor therapy for the prevention of esophageal adenocarcinoma in Barrett’s esophagus. J Natl Cancer Inst. 2002;94:406–7.

    Article  PubMed  Google Scholar 

  37. Han L, Tang MM, Xu X, Jiang B, Huang J, Feng X, et al. LTBP2 is a prognostic marker in head and neck squamous cell carcinoma. Oncotarget. 2016;7:45052–9.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74.

    Article  CAS  PubMed  Google Scholar 

  39. Hansen NUB, Willumsen N, Sand JMB, Larsen L, Karsdal MA, Leeming DJ. Type VIII collagen is elevated in diseases associated with angiogenesis and vascular remodeling. Clin Biochem. 2016;49:903–8.

    Article  CAS  PubMed  Google Scholar 

  40. Huynh PT, Beswick EJ, Coronado YA, Johnson P, O’Connell MR, Watts T, et al. CD90(+) stromal cells are the major source of IL-6 which supports cancer stem-like cells and inflammation in colorectal cancer. Int J Cancer. 2016;138:1971–81.

    Article  CAS  PubMed  Google Scholar 

  41. Isella C, Terrasi A, Bellomo SE, Petti C, Galatola G, Muratore A, et al. Stromal contribution to the colorectal cancer transcriptome. Nat Genet. 2015;47:312–9.

    Article  CAS  PubMed  Google Scholar 

  42. Kalluri R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 2016;16:582–98.

    Article  CAS  PubMed  Google Scholar 

  43. Karagiannis GS, Petraki C, Prassas I, Saraon P, Musrap N, Dimitromanolakis A, et al. Proteomic signatures of the desmoplastic invasion front reveal collagen type XII as a marker of myofibroblastic differentiation during colorectal cancer metastasis. Oncotarget. 2012;3:265–85.

    Article  Google Scholar 

  44. Kato Y, Lewalle J-M, Baba Y, Tsukuda M, Sakai N, Baba M, et al. Induction of SPARC by VEGF in human vascular endothelial cells. Biochem Biophys Res Commun. 2001;287:422–6.

    Article  CAS  PubMed  Google Scholar 

  45. Kessenbrock K, Plaks V, Werb Z. Matrix metalloproteinases: regulators of the tumor microenvironment. Cell. 2010;141:52–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Kim JH, Cho NY, Bae JM, Kim KJ, Rhee YY, Lee HS, et al. Nuclear maspin expression correlates with the CpG island methylator phenotype and tumor aggressiveness in colorectal cancer. Int J Clin Exp Pathol. 2015;8:1920–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Lane TF, Sage EH. The biology of SPARC, a protein that modulates cell-matrix interactions. FASEB J. 1994;8:163–73.

    CAS  PubMed  Google Scholar 

  48. Ledda F, Bravo AI, Adris S, Bover L, Mordoh J, Podhajcer OL. The expression of the secreted protein acidic and rich in cysteine (SPARC) is associated with the neoplastic progression of human melanoma. J Invest Dermatol. 1997;108:210–4.

    Article  CAS  PubMed  Google Scholar 

  49. Ledda MF, Adris S, Bravo AI, Kairiyama C, Bover L, Chernajovsky Y, et al. Suppression of SPARC expression by antisense RNA abrogates the tumorigenicity of human melanoma cells. Nat Med. 1997;3:171–6.

    Article  CAS  PubMed  Google Scholar 

  50. Liao D, Luo Y, Markowitz D, Xiang R, Reisfeld RA. Cancer associated fibroblasts promote tumor growth and metastasis by modulating the tumor immune microenvironment in a 4T1 murine breast cancer model. PLoS ONE. 2009;4:e7965.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Maris P, Blomme A, Palacios AP, Costanza B, Bellahcene A, Bianchi E, et al. Asporin is a fibroblast-derived TGF-beta1 inhibitor and a tumor suppressor associated with good prognosis in breast cancer. PLoS Med. 2015;12:e1001871.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Marsh T, Pietras K, McAllister SS. Fibroblasts as architects of cancer pathogenesis. Biochim Biophys Acta. 2013;1832:1070–8.

    Article  CAS  PubMed  Google Scholar 

  53. McClung HM, Thomas SL, Osenkowski P, Toth M, Menon P, Raz A, et al. SPARC upregulates MT1-MMP expression, MMP-2 activation, and the secretion and cleavage of galectin-3 in U87MG glioma cells. Neurosci Lett. 2007;419:172–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Mehdawi L, Osman J, Topi G, Sjolander A. High tumor mast cell density is associated with longer survival of colon cancer patients. Acta Oncol. 2016;55:1434–42.

    Article  CAS  PubMed  Google Scholar 

  55. Mook OR, Frederiks WM, Van Noorden CJ. The role of gelatinases in colorectal cancer progression and metastasis. Biochim Biophys Acta. 2004;1705:69–89.

    CAS  PubMed  Google Scholar 

  56. Moren A, Olofsson A, Stenman G, Sahlin P, Kanzaki T, Claesson-Welsh L, et al. Identification and characterization of LTBP-2, a novel latent transforming growth factor-beta-binding protein. J Biol Chem. 1994;269:32469–78.

    CAS  PubMed  Google Scholar 

  57. Nagaraju GP, Dontula R, El-Rayes BF, Lakka SS. Molecular mechanisms underlying the divergent roles of SPARC in human carcinogenesis. Carcinogenesis. 2014;35:967–73.

    Article  CAS  PubMed  Google Scholar 

  58. Öhlund D, Elyada E, Tuveson D. Fibroblast heterogeneity in the cancer wound. J Exp Med. 2014;211:1503–23.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Paszek MJ, Zahir N, Johnson KR, Lakins JN, Rozenberg GI, Gefen A, et al. Tensional homeostasis and the malignant phenotype. Cancer Cell. 2005;8:241–54.

    Article  CAS  PubMed  Google Scholar 

  60. Qian BZ, Pollard JW. Macrophage diversity enhances tumor progression and metastasis. Cell. 2010;141:39–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Qiao J, Fang C-Y, Chen S-X, Wang X-Q, Cui S-J, Liu X-H, et al. Stroma derived COL6A3 is a potential prognosis marker of colorectal carcinoma revealed by quantitative proteomics. Oncotarget. 2015;6:29929–46.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Ren Y, Lu H, Zhao D, Ou Y, Yu K, Gu J, et al. LTPB2 acts as a prognostic factor and promotes progression of cervical adenocarcinoma. Am J Transl Res. 2015;7:1095–105.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. Ricci-Vitiani L, Lombardi DG, Pilozzi E, Biffoni M, Todaro M, Peschle C, et al. Identification and expansion of human colon-cancer-initiating cells. Nature. 2007;445:111–5.

    Article  CAS  PubMed  Google Scholar 

  64. Robson H, Anderson E, James RD, Schofield PF. Transforming growth factor β 1 expression in human colorectal tumours: an independent prognostic marker in a subgroup of poor prognosis patients. Br J Cancer. 1996;74:753–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Rogler G. Chronic ulcerative colitis and colorectal cancer. Cancer Lett. 2014;345:235–41.

    Article  CAS  PubMed  Google Scholar 

  66. Rubio CA, Kaufeldt A, Bjork J, Jaramillo E. Maspin, a marker of serrated colorectal polyps. Anticancer Res. 2015;35:4139–44.

    CAS  PubMed  Google Scholar 

  67. Rupp C, Scherzer M, Rudisch A, Unger C, Haslinger C, Schweifer N, et al. IGFBP7, a novel tumor stroma marker, with growth-promoting effects in colon cancer through a paracrine tumor-stroma interaction. Oncogene. 2015;34:815–25.

    Article  CAS  PubMed  Google Scholar 

  68. Said N. Roles of SPARC in urothelial carcinogenesis, progression and metastasis. Oncotarget. 2016;7:67574–85.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Sangaletti S, Tripodo C, Santangelo A, Castioni N, Portararo P, Gulino A, et al. Mesenchymal transition of high-grade breast carcinomas depends on extracellular matrix control of myeloid suppressor cell activity. Cell Rep. 2016;17:233–48.

    Article  CAS  PubMed  Google Scholar 

  70. Sanz-Pamplona R, Berenguer A, Cordero D, Mollevi DG, Crous-Bou M, Sole X, et al. Aberrant gene expression in mucosa adjacent to tumor reveals a molecular crosstalk in colon cancer. Mol Cancer. 2014;13:46.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Sato N, Fukushima N, Maehara N, Matsubayashi H, Koopmann J, Su GH, et al. SPARC/osteonectin is a frequent target for aberrant methylation in pancreatic adenocarcinoma and a mediator of tumor-stromal interactions. Oncogene. 2003;22:5021–30.

    Article  CAS  PubMed  Google Scholar 

  72. Sharma S, Xing F, Liu Y, Wu K, Said N, Pochampally R, et al. Secreted protein acidic and rich in cysteine (SPARC) mediates metastatic dormancy of prostate cancer in bone. J Biol Chem. 2016;291:19351–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Simkova D, Kharaishvili G, Slabakova E, Murray PG, Bouchal J. Glycoprotein asporin as a novel player in tumour microenvironment and cancer progression. Biomed Pap Med Fac Univ Palacky Olomouc Czech Repub. 2016;160:467–73.

    PubMed  Google Scholar 

  74. Slany A, Bileck A, Kreutz D, Mayer RL, Muqaku B, Gerner C. Contribution of Human fibroblasts and endothelial cells to the hallmarks of inflammation as determined by proteome profiling. Mol Cell Proteomics MCP. 2016;15:1982–97.

    Article  CAS  PubMed  Google Scholar 

  75. Slany A, Bileck A, Kreutz D, Mayer RL, Muqaku B, Gerner C. Contribution of human fibroblasts and endothelial cells to the hallmarks of inflammation as determined by proteome profiling. Mol Cell Proteomics. 2016;15:1982–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Slany A, Meshcheryakova A, Beer A, Ankersmit HJ, Paulitschke C, Gerner C. Plasticity of fibroblasts demonstrated by tissue-specific and function-related proteome profiling. Clin Proteomics. 2014;11:41.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:1–25.

    Article  Google Scholar 

  78. Sonnante AM, Correale M, Linsalata M, Di Leo A, Guerra V. Circulating levels of matrix metalloproteinase-9 in patients with colorectal cancer. Scand J Gastroenterol. 2000;35:671–2.

    Article  CAS  PubMed  Google Scholar 

  79. Sun R, Wu J, Chen Y, Lu M, Zhang S, Lu D, et al. Down regulation of Thrombospondin2 predicts poor prognosis in patients with gastric cancer. Mol Cancer. 2014;13:225.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Suzuki S, Ichikawa Y, Nakagawa K, Kumamoto T, Mori R, Matsuyama R, et al. High infiltration of mast cells positive to tryptase predicts worse outcome following resection of colorectal liver metastases. BMC Cancer. 2015;15:840.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Tichet M, Prod’Homme V, Fenouille N, Ambrosetti D, Mallavialle A, Cerezo M, et al. Tumour-derived SPARC drives vascular permeability and extravasation through endothelial VCAM1 signalling to promote metastasis. Nat Commun. 2015;6:6993.

    Article  CAS  PubMed  Google Scholar 

  82. Torres S, Bartolome RA, Mendes M, Barderas R, Fernandez-Acenero MJ, Pelaez-Garcia A, et al. Proteome profiling of cancer-associated fibroblasts identifies novel proinflammatory signatures and prognostic markers for colorectal cancer. Clin Cancer Res. 2013;19:6006–19.

    Article  CAS  PubMed  Google Scholar 

  83. Tsamandas AC, Kardamakis D, Ravazoula P, Zolota V, Salakou S, Tepetes K, et al. The potential role of TGFβ1, TGFβ2and TGFβ3 protein expression in colorectalcarcinomas. Strahlenther Onkol. 2004;180:201–8.

    Article  PubMed  Google Scholar 

  84. Umekita Yoshihisa, Souda Masakazu, Yoshida Souda. Expression of maspin in colorectal cancer. In Vivo. 2006;20:797–800.

    CAS  PubMed  Google Scholar 

  85. Vizcaino JA, Deutsch EW, Wang R, Csordas A, Reisinger F, Rios D, et al. ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat Biotechnol. 2014;32:223–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Wagenaar-Miller RA, Gorden L, Matrisian LM. Matrix metalloproteinases in colorectal cancer: is it worth talking about? Cancer Metastasis Rev. 2004;23:119–35.

    Article  CAS  PubMed  Google Scholar 

  87. Wang C, Wang G, Zhang L, Pan J, Wei Y. Latent transforming growth factor beta binding protein 2 (LTBP2) as a novel biomarker for the diagnosis and prognosis of pancreatic carcinoma. Med Sci Monit. 2017;23:3232–9.

    Article  PubMed  PubMed Central  Google Scholar 

  88. Wang D, Fu L, Sun H, Guo L, DuBois RN. Prostaglandin E2 promotes colorectal cancer stem cell expansion and metastasis in mice. Gastroenterology. 2015;149(1884–1895):e1884.

    Article  Google Scholar 

  89. Wang TH, Hsia SM, Shieh TM. Lysyl oxidase and the tumor microenvironment. Int J Mol Sci. 2016;18:62.

    Article  PubMed Central  Google Scholar 

  90. Wang X, Zhang L, Li H, Sun W, Zhang H, Lai M. THBS2 is a potential prognostic biomarker in colorectal cancer. Sci Rep. 2016;6:33366.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Werb Z, Lu P. The role of stroma in tumor development. Cancer J (Sudbury, Mass.). 2015;21:250–3.

    Article  CAS  Google Scholar 

  92. Wu X, Zou Y, He X, Yuan R, Chen Y, Lan N. Tumor-infiltrating mast cells in colorectal cancer as a poor prognostic factor. Int J Surg Pathol. 2013;21:111–20.

    Article  PubMed  Google Scholar 

  93. Zhang B, Wang J, Wang X, Zhu J, Liu Q, Shi Z, et al. Proteogenomic characterization of human colon and rectal cancer. Nature. 2014;513:382–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Zucker S, Hymowitz M, Conner C, Zarrabi HM, Hurewitz AN, Matrisian L, et al. Measurement of matrix metalloproteinases and tissue inhibitors of metalloproteinases in blood and tissues. Clinical and experimental applications. Ann N Y Acad Sci. 1999;878:212–27.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

Study design: CG, BM, DD, AB. Experiments: DD, AB. Tissue Acquisition and pathological evaluation: AB. Morphometric analysis: DD, GT. Analysis of published data sets: ZNE, TM, DD. Writing of the manuscript: CG, BM, DD, AB. Reviewing of the manuscript: All authors read and approved the final manuscript.

Acknowledgements

The authors gratefully acknowledge expert technical assistance by Xenia Hudec.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Protein data was submitted to the ProteomeXchange Consortium via the PRIDE partner repository and can be accessed via http://www.proteomeexchange.org with the identifier PXD00.

Consent for publication

All authors have seen the manuscript and agree to its publication.

Ethics approval and consent to participate

The study was approved by the Ethics Commission of the Medical University of Vienna (EK 1659/2012). Patients suffering from colorectal cancer who underwent surgery at the General Hospital of Vienna were asked for their informed consent to use part of the resected tissue for analysis of tumor protein markers.

Funding

This project was supported by grants from the Austrian National Bank (15655) and the Austrian Science Fund (P26218).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Author notes

  1. Daniel Drev and Andrea Bileck have contributed equally to this work

  2. Christopher Gerner and Brigitte Marian have contributed equally and are co-corresponding authors

Authors and Affiliations

  1. Department of Medicine 1, Institute of Cancer Research and Comprehensive Cancer Center, Medical University of Vienna, Vienna, Austria

    Daniel Drev, Zeynep N. Erdem, Thomas Mohr, Gerald Timelthaler & Brigitte Marian

  2. Institute of Analytical Chemistry, University of Vienna, Vienna, Austria

    Andrea Bileck & Christopher Gerner

  3. Clinical Institute of Pathology, Medical University of Vienna, Vienna, Austria

    Andrea Beer

Authors
  1. Daniel Drev
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Andrea Bileck
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Zeynep N. Erdem
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Thomas Mohr
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Gerald Timelthaler
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Andrea Beer
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Christopher Gerner
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Brigitte Marian
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Christopher Gerner or Brigitte Marian.

Additional files

12014_2017_9168_MOESM1_ESM.docx

Additional file 1: Figures S1: Epithelia/Stromal separation. Figure S2. mRNA levels for proteins increased in the tumor. Figure S3. mRNA levels of proteins decreased in the tumor. Figure S4. Quantification of THBS2 staining in IHC slides from the original 6 patients.

Additional file 2. Table S1. Proteins identified in colorectal tumors and normal mucosa.

Additional file 3. Table S2. Significantly altered proteins in colorectal tumors.

Additional file 4. Table S3. GO-term enrichment analysis based on biological process.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Drev, D., Bileck, A., Erdem, Z.N. et al. Proteomic profiling identifies markers for inflammation-related tumor–fibroblast interaction. Clin Proteom 14, 33 (2017). https://doi.org/10.1186/s12014-017-9168-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12014-017-9168-7

Keywords

Clinical Proteomics

ISSN: 1559-0275