Skip to main content

Correlation of the tumor escape phenotype with loss of PRELP expression in melanoma

Abstract

Background

Despite immunotherapies having revolutionized the treatment of advanced cutaneous melanoma, effective and durable responses were only reported in a few patients. A better understanding of the interaction of melanoma cells with the microenvironment, including extracellular matrix (ECM) components, might provide novel therapeutic options. Although the ECM has been linked to several hallmarks of cancer, little information is available regarding the expression and function of the ECM protein purine-arginine-rich and leucine-rich protein (PRELP) in cancer, including melanoma.

Methods

The structural integrity, expression and function of PRELP, its correlation with the expression of immune modulatory molecules, immune cell infiltration and clinical parameters were determined using standard methods and/or bioinformatics.

Results

Bioinformatics analysis revealed a heterogeneous, but statistically significant reduced PRELP expression in available datasets of skin cutaneous melanoma when compared to adjacent normal tissues, which was associated with reduced patients’ survival, low expression levels of components of the MHC class I antigen processing machinery (APM) and interferon (IFN)-γ signal transduction pathway, but increased expression of the transforming growth factor (TGF)-β isoform 1 (TFGB1) and TGF-β receptor 1 (TGFBR1). In addition, a high frequency of intra-tumoral T cells directly correlated with the expression of MHC class I and PRELP as well as the T cell attractant CCL5 in melanoma lesions. Marginal to low PRELP expression levels were found in the 47/49 human melanoma cell lines analysis. Transfection of PRELP into melanoma cell lines restored MHC class I surface expression due to transcriptional upregulation of major MHC class I APM and IFN-γ pathway components. In addition, PRELP overexpression is accompanied by high CCL5 secretion levels in cell supernatant, an impaired TGF-β signaling as well as a reduced cell proliferation, migration and invasion of melanoma cells.

Conclusions

Our findings suggest that PRELP induces the expression of MHC class I and CCL5 in melanoma, which might be involved in an enhanced T cell recruitment and immunogenicity associated with an improved patients’ outcome. Therefore, PRELP might serve as a marker for predicting disease progression and its recovery could revert the tumorigenic phenotype, which represents a novel therapeutic option for melanoma.

Background

Melanoma, a malignant skin tumor, is estimated as the fifth most common cancer worldwide, with a steadily increasing incidence over the past decades. Although early-stage melanoma has a good prognosis after surgery, relatively small melanoma often exerts a high metastatic potential with a low 5 years survival rate of 15.7% before 2011 [1]. Due to its reported immunogenicity, different T cell-based immunotherapies have been developed for the treatment of melanoma. In particular, immune checkpoint inhibitors (ICPi) targeting the programmed cell death protein 1 (PD-1), PD-1 ligand (PD-L1) and/or the cytotoxic T lymphocyte antigen 4 (CTLA4) have significantly improved the outcome of patients with advanced melanoma [2, 3]. Since only a limited number of patients have a long-term response to ICPi, biomarkers for predicting response to immunotherapy and resistance development are urgently required. Recently, surface expression of HLA class I has been shown to serve as a predictive biomarker for ICPi response [4] or acquired resistance [5,6,7], suggesting that loss or downregulation of HLA class I expression might drive immune escape of melanoma. Indeed, Seliger and co-authors [8] provided one of the first evidences about structural alterations in components of the MHC class I antigen processing machinery (APM) and the interferon (IFN) signaling pathway in melanoma cells resulting in a loss of MHC class I surface expression. This could be reverted by gene transfer of APM and IFN components, which was associated with increased CD8+ T cell responses. Although genetic abnormalities are rare, the frequent downregulation of MHC class I expression in melanoma is mainly due to a deregulation of MHC class I APM molecules and/or molecules involved in IFN signaling at distinct levels [9,10,11]. It is noteworthy that a link between a low MHC class I surface antigens and low numbers of CD8+ T cells exists, which is associated with progression of melanoma [12, 13].

Since one primary goal in the field of tumor immunology is to increase the immunogenicity of tumor cells by re-establishing and maintaining MHC class I surface expression, the identification of key molecules or substances overriding these tumor intrinsic escape routes will help to improve durable elimination of cancer cells by CD8+ T cells.

In this context, proteoglycans (PGs) as crucial constituents of the extracellular matrix (ECM), which differ between tissues, developmental stages and (patho) physiological conditions, have been suggested to be involved in the immunogenicity of tumor cells by orchestrating cellular processes mediated by cell–cell and cell–matrix interactions [14] exhibiting tumor-suppressing and tumor-initiating properties [15,16,17]. Recently, a member of the small leucine-rich repeat proteoglycan (SLRP) family, named biglycan (BGN), has been shown as a regulator of MHC class I expression in HER‐2/neu- [18, 19] and K-RAS- [20] transformed cancer cells. However, little information exists about the role of other SLRPs on HLA class I expression, including the purine-arginine-rich and leucine-rich repeat protein (PRELP), also known as MST161, MSTp161, prolargin or SLRR2A [21]. PRELP consists of a core protein and interacts with collagen fibrils, heparin and heparan sulfate [22,23,24,25,26,27], proteins of the complement system, various membrane proteins [28, 29] as well as signaling pathways like the wnt and TGF-β [30]. The latter resulted in the suppression of Smad2 phosphorylation, leading to the inhibition of the TGF-β pathway [31]. PRELP is constitutively expressed in healthy cartilage, lung, kidney, liver and skin, secreted into the ECM [32] and involved in maintaining normal cellular structures and epithelial cell integrity. Next to its physiological function, PRELP has been reported as a biomarker in different human cancers, but its role is controversially discussed dependent on the tumor type analyzed [33,34,35] regarding its correlation with the patients’ survival [34, 36] and prognostic potential [30, 36,37,38,39,40,41]. While overexpression of PRELP has been reported to have tumor suppressive activity in solid tumors characterized by a reduced invasion, anchorage-independent growth, metastasis formation and tumorgenicity [42], PRELP expression is involved in the pathobiology of chronic lymphocytic leukemia [43]. Since no information exists for skin cancer, we investigated the expression, function and clinical relevance of PRELP in melanoma. Based on data from The Cancer Genome Atlas (TCGA), melanoma patients were divided into a PRELPhigh and PRELPlow group and associated with clinical parameters, expression of HLA class I APM and IFN-γ signal pathway components, CCL5 as well as the tumor immune infiltration. These results were further confirmed in melanoma cell lines and the in vitro activity effect of PRELP on cell growth and immunogenicity was assessed.

Methods

Cell lines, treatment and transfection

Murine PRELPlow B16 F10 melanoma cells were purchased from the American Tissue Culture Collection. The different human melanoma cell lines were provided by S. Ferrone (Harvard Medical School, Boston, USA). Murine B16 F10 cells and the human melanoma cell lines were maintained in Eagles modified essential medium (EMEM, Lonza, Basel, Switzerland) supplemented with 10% fetal calf serum (FCS), 2 mM glutamine, 100 U/ml penicillin and 100 µg/ml streptomycin (PAA; Pasching, Austria) and cultured at 37 °C in 5% CO2 humified air.

PRELP-overexpressing (PRELPhigh) murine B16 F10 cells and human overexpressing cells (Buf 1088) as well as mock controls (PRELPlow/mock) were generated as recently described for BGN [20]. Briefly, different melanoma cell lines were stably transfected with a PRELP expression vector according to Recktenwald and co-authors [44]. All stable PRELP transfectants and vector controls were maintained in complete EMEM supplemented with 1 mg/ml G418 (PAA).

Determination of cell proliferation and migration

The growth properties of PRELPlow and PRELPhigh melanoma cell lines were determined as recently described [45, 46]. Cell proliferation was analyzed after 48 h in triplicates using the cell proliferation kit II (Roche Applied Science, Penzberg, Germany) according to the manufacturers' instructions.

For determining cell migration, 5 × 104 cells were plated into a trans-well chamber using a gradient of 0.5 to 10% FCS as an attractant in the lower chamber. After 24 h, the number of migrated cells into the lower chamber was determined using the Cell Titer Glo Luminescence cell viability assay (Promega) according to the manufacturer’s instructions. The luminescence was normalized against the luminescence cells that were directly seeded into the bottom of the trans-well plate. The results are expressed as % of migrated cells from at least three independent experiments using triplicates.

RNA isolation, reverse transcription and quantitative real-time RT-PCR

Total cellular RNA was isolated from cell lines using the NucleoSpin RNA II kit (Macherey–Nagel, Dueren, Germany) followed by reverse transcription of 2 µg total RNA into cDNA and PCR as previously described [47] using target-specific primers listed in Additional file 1: Table S1. Experiments were independently repeated three times.

Determination of the APM promoter activity by luciferase (luc) assays

TAP1, TAP2, TAPBP and PSMB9 promoter sequences were amplified from genomic DNA and then cloned into the pGl3 luciferase (luc) vector (Promega, Fitchburg, USA) as recently described [20, 48]. For transient transfections, 1 × 105 cells were cultured in 100 µl OptiMEM (Invitrogen), followed by transfection with 0.3 µg promoter constructs and 0.016 µg β-galactosidase (β-gal) vector using Lipofectamine 2000 (Invitrogen, Waltham, USA) as transfection reagent according to the manufacturer’s instructions. 48 h after transfection, the luc activity was determined by adding the luc substrate (Promega) using a luminometer and normalized to the transfection efficiency determined by ß-gal enzyme activity. Experiments were independently done three times using triplicates.

Western blot analysis

For Western blot analysis, 5 × 106 cells were harvested, and proteins were extracted. 30 µg protein/lane was loaded in 10% SDS-PAGE gels, transferred onto nitrocellulose membranes (Schleicher & Schuell, Dassel, Germany) and stained with Ponceau S as previously described [47]. Membranes were incubated overnight at 4 °C with primary monoclonal antibodies (mAbs) directed against PRELP (Invitrogen, Waltham, USA), β-actin (Sigma-Aldrich, St. Louis, USA) and/or glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (Cell Signaling Technology, Danvers, USA) as a loading control, respectively, followed by incubation for one hour with a horseradish peroxidase linked secondary antibody (Ab; Cell Signaling Technology) and developed using the ECL method. Chemiluminescence signals were visualized by the Lumi-Light Western Blotting Substrate (Roche Diagnostics) and recorded with a LAS3000 system (Fuji, Tokyo, Japan).

Flow cytometry

For measuring MHC class I surface expression in murine and human melanoma cell lines, the following mAbs were used: the anti-H-2D mAb (Cedarlane Laboratories LTD, Burlington, Canada) and the anti-HLA class I-specific mAb (Beckman Coulter). Flow cytometric analysis was performed as previously described [48]. Briefly, 5 × 105 cells were incubated with the appropriate amount of the respective antibodies at 4 °C for 30 min before the stained cells were measured on a NAVIOS (Becton Dickinson, Franklin Lakes, USA) and subsequently analyzed with the Kaluza Analysis Software (Beckman Coulter). The data are presented as mean specific fluorescence intensity (MFI) from three independent experiments.

For cell cycle analysis, 1 × 106 cells were cultured in 0.1% FCS for 48 h followed by their cultivation in 10% FCS before cell cycle analysis [49]. Then, cells were fixed with ice-cold 70% ethanol and stained with propidium iodide (Sigma). The cells were subjected to flow cytometry using the Kaluza Analysis Software (Beckman Coulter). Three independent experiments were performed.

Quantification of CCL5 by ELISA

For quantification of CCL5 in the supernatants of PRELPlow and PRELPhigh melanoma cells, 2 × 104 cells were seeded in 24 well plates. After an overnight incubation to allow cell adhesion, the medium was changed and cells were incubated for an additional 72 h. The supernatants were then harvested and the amount of CCL5 released by the melanoma cells was assessed by an ELISA (ELISA MAX™ Deluxe set, BioLegend) according to the manufacturer's instructions.

CD107a degranulation assay

Peripheral blood mononuclear cells (PBMC) were purified from healthy human donors. Buffy coat was prepared by a gradient density centrifugation, then stimulated for 18 h with 1 ng/ml IL-12, 5 ng/ml IL-15 (Immunotools, Friesoythe, Germany) and 50 ng/ml IL-18 (Biovision, Milpitas, CA, USA) in X-vivo15 (Lonza) medium, before they were incubated together with PRELPlow vs. PRELPhigh target cells. For the CD107a degranulation assay [50], the anti-CD107a Ab was added after 1 h of co-culture, while staining the cells with Abs directed against CD3, CD16 and CD56 to identify NK cells was performed after 4 h [19].

Bioinformatics analyses of public datasets

The transcriptome profiles and clinical data from healthy individuals and cancer patients were obtained from the TCGA dataset (http://cancergenome.nih.gov/) and microarray data of the NCBI Gene Expression Omnibus (NCBI GEO) [51]. The datasets were available on GENT2 [52], HCMDB [53], UALCAN [54], R2 Genomics (http://r2.amc.nl), UCSC Xena (UCSC Xena (http://xena.ucsc.edu) platform and cBioPortal [55]. The number of human samples, cancers vs. adjuvant non-tumorigenic tissues and melanoma cell lines analyzed [56,57,58,59,60] are presented in Additional file 1: Table S2.

Using the Omic Horizon Expression database [61], PRELP expression was downloaded as transcripts per kilobase million (TPM) from multiple tissues and different sub-tissues according to the anatomical site. PRELP gene expression data were downloaded from the NCBI GEO database generated by the two microarray platforms Affymetrix U133A or U133Plus2. Then, the data were stratified into cancer and normal samples across 72 tissues by GENT2 system and relative PRELP expression was analyzed and correlated to cancer stage, age and clinicopathologic features. The comparison between cancer and normal tissues are presented as Log2FC (fold change calculates for logarithm to base 2) and p values are given.

PrognoScan (http://www.prognoscan.org/) was used to evaluate the relationship between gene expression and the survival of patients. Based on the gene expression levels, samples were divided into PRELPhigh and PRELPlow expression groups and the differences in the risk and the overall survival (OS) between both groups (cut off 0.42) were analyzed and cumulative survival plots and Kaplan–Meier curves were constructed with the log-rank test. By using UCSC Xena, melanoma patients were categorized based on the global PRELP expression. The global PRELP expression ranges from log2 5.7 to 13. In the PRELP global range, PRELPhigh refers to log2 < 10.6; the patients with log2 of < 8.1 were grouped as PRELPlow. The Kaplan–Meier results are presented as log-rank test values and p values are given.

The OncoPrint algorithm in cBioPortal database [55] was used to analyze genetic alterations, such as amplification, deep deletions and mutations. A concise and compact graphical summary was generated for genomic alterations in PRELP across cutaneous melanoma samples. PRELP correlation analyses were performed with the R2: Genomics Analysis and Visualization Platform and correlation statistics were shown as Pearson correlation coefficients. The correlation of tumor immune cell infiltration and PRELP expression was analyzed as a co-variable to generate corrected multivariate Cox proportional hazard models. Results from the Cox regression models are presented as Hazard Ratios (HR) of z-scores and p values. In addition, the Tumor Immune Estimation Resource (TIMER) [62] was used to explore the associations of PRELP expression with the levels of immune cell infiltrates by employing the xCell algorithm [63]. The infiltration of CD8+ T cell subpopulations in tumor samples and PRELP were used for Cox univariate and multivariable regression analysis (p < 0.05), generating an independent Cox model.

Statistical analysis

Microsoft Excel-Office 365 and R (programming language) were used for student’s t-test and one-way ANOVA. A p-value of < 0.05 was considered as statistically significant (*p < 0.05; **p < 0.01; ***p < 0.001). The survival curve was derived from the Kaplan–Meier method. The log-rank test was used to compare the survival rate.

Results

Reduced PRELP expression in melanoma

Members of the SLRP family exhibit pro-tumorigenic and anti-tumorigenic activities, dependent on the tissue, context, localization and tumor entity [14, 64,65,66], but limited information exists about their effects on the immunogenicity of tumor cells [26, 27]. A tissue-wide gene expression profile analysis of PRELP demonstrated a heterogeneous expression among 123 healthy human tissues and sub-tissues using the Omic Horizon Expression database. PRELP was highly expressed in cartilage (2963.97 TPM), aorta (1774.93 TPM), brain-superficial zone (1059.59 TPM), kidney-inner medulla (891.05 TPM) and Achilles tendon (840.84 TPM), whereas skin showed a moderate PRELP expression (110.45 TPM; Additional file 1: Fig. S1A). Since the expression of PRELP has only been analyzed in some tumor types, like bladder, pancreatic, colorectal, hepatocellular carcinoma and retinoblastoma [30, 36, 38, 41, 67], PRELP expression was determined in different cancers using the two independent databases GENT2–U133Plus2 and GENT2–U133A consisting of 5.487 normal tissues and 35.806 cancer samples. As shown in Fig. 1A, PRELP expression was generally downregulated in all solid cancers (log2FC = − 0.758 and p < 0.001), but there exist tumor-dependent differences in the extent of downregulation, which was more pronounced in melanoma (log2FC = − 1.583 and p < 0.001) when compared to corresponding normal tissues. However, PRELP is heterogeneously expressed in melanoma lesions leading to the assignment of PRELPlow and PRELPhigh melanoma. These data are in line with cell line database analyses demonstrating downregulation of PRELP expression in skin cancer cells (log2FC = − 0.77 and p < 0.031) when compared to their normal counterparts, such as melanocytes and keratinocytes (GENT2—skin cancer cell lines) (Fig. 1B). In addition, 48 melanoma cell lines were tested for PRELP mRNA expression using qRT-PCR. Only 2/48 melanoma cell lines express high PRELP mRNA levels (Additional file 1: Fig. S1B), but the underlying mechanism of the lack of PRELP expression has not yet been determined.

Fig. 1
figure 1

Correlation of PRELP expression in tumors of distinct origin and healthy controls and its clinical relevance. A Reduced PRELP expression in solid tumors compared to normal healthy tissues. PRELP expression was determined in 5.487 normal tissue samples and 35.806 cancer tissues. B Decreased PRELP expression in skin cancer cell lines compared to its normal counterparts. Box plots comparing PRELP expression in skin cancer (247 samples) vs. normal skin cell lines (32 samples) are shown. The data are expressed as log2FC-fold change in a logarithmic scale (base 2). C Distinct PRELP expression patterns in primary and metastatic skin tissues. A box plot comparing PRELP expression between metastases vs. primary tumors is shown. D Increased overall survival of melanoma patients with high PRELP expression levels in tumors. The Kaplan–Meier survival curve of PRELPhigh and PRELPlow melanoma patients from the GEO (ID GSE98394; [68]) melanoma dataset was analyzed and presented as a Kaplan–Meier plot as described in “Methods”. The statistical significance is presented as *p < 0.05; **p < 0.01; ***p < 0.001

Clinical relevance of PRELP expression in melanoma

To get in-depth insights into the role of PRELP expression in cutaneous skin cancers, PRELP transcription was analyzed in an extensive series of primary and metastatic melanoma lesions and normal controls and compared to clinical parameters of melanoma patients (Tumor Skin Cutaneous Melanoma—TCGA, Additional file 1: Table S2). The PRELP mRNA expression levels were lower in primary melanoma lesions than in normal skin (log2FC = − 1.245 and p = 2.485e−3) and further reduced in metastasis (log2FC = − 0.558 and p = 3.047e−2) (Fig. 1C), but independent of nodal metastasis status, disease stage, patient’s age and gender (Additional file 1: Fig. S1C; GEO ID: GSE7553) [69]. Kaplan–Meier analysis revealed that PRELP expression levels correlated with the overall survival (OS) of melanoma patients with a reduced patients’ survival rate (HR = − 0.2635; p = 0.0121) in PRELPlow melanoma lesions (Fig. 1D).

Link of PRELP expression with the immune escape phenotype

In order to determine whether PRELP expression is inversely associated with an immune escape phenotype, the HLA class I expression was analyzed in PRELPhigh vs. PRELPlow melanoma lesions by in silico analysis using three different human melanoma data sets: (i) mixed melanoma metastasis (83 samples; GEO ID: GSE8401) and (ii) tumor skin cutaneous melanoma (480 samples; TCGA) (Additional file 1: Table S3). A statistically significant positive correlation to PRELP was found for HLA class I antigens (Fig. 2A) and TAPBP (Fig. 2B) in two datasets for at least one proteasome subunit (Additional file 1: Table S3). By comparison of the HLA-A expression in 133 skin cancer cell lines (GENT2—skin cancer cell lines) with PRELP, a significant positive correlation, Pearson correlation coefficient, r = 0.37 and p = 1.98E−09 (Additional file 1: Fig. S2A). In addition, a positive correlation was also detected between PRELP and IRF5 expression (Fig. 2C). The data from melanoma lesions were in line with a strong downregulation of MHC class I surface expression in murine and human melanoma cells, which exhibit low to marginal PRELP expression (PRELPlow/neg) in 46/48 melanoma cell lines analyzed. Thus, a link between high PRELP expression levels and increased immunogenicity of melanoma cells is shown in vitro and in vivo.

Fig. 2
figure 2

Correlation of PRELP expression with HLA class I expression. AC Correlation of PRELP expression with HLA-A, TAPBP and IRF5. Using the GEO (ID GSE8401) melanoma dataset and R2 Genomics for PRELP expression was correlated to HLA-A (A), TAPBP (B) and IRF5 (C) mRNA levels. DF Correlation of PRELPhigh vs. PRELPlow melanoma with the expression of HLA-A and overall survival. Data of Tumor Skin Cutaneous Melanoma—TCGA were analyzed and results are depicted as Kaplan–Meier plots. Correlation of the OS for HLA-A expression independent of PRELP expression (D), PRELPhigh (E) and PRELPlow (F) samples are shown in Kaplan–Meier plots

Further analysis of the Tumor Skin Cutaneous Melanoma—TCGA dataset revealed that a high HLA-A expression level favored an increased OS of melanoma patients (p = 0.0006; log-rank test = 11.76) regardless of PRELP expression levels (log2 5.7 to 13) (Fig. 2D), but PRELPhigh expression (log2 > 10.6) was accompanied by high levels of HLA class I expression and increased patients’ OS (p = 0.01168; log-rank test = 6.359) (Fig. 2E). In contrast, in PRELPlow melanoma (log2 < 8.1), HLA-A failed to predict the survival of melanoma patients (p = 0.8903; log-rank test = 0.019; Fig. 2F).

Reversion of the tumorigenic phenotype by PRELP overexpression in melanoma cell lines

Transfection of PRELPlow/neg B16F10 and Buf1088 cells with a PRELP expression vector restored PRELP mRNA (Additional file 1: Fig. S3) and protein expression (Fig. 3A). When compared to PRELPlow/neg B16F10 and Buf1088 cells PRELP overexpressing (PRELPhigh) B16F10 and Buf1088 transfectants exhibited altered growth properties, such as reduced cell proliferation (Fig. 3B) and reduced migration determined by trans-well invasion assay (Fig. 3C) accompanied by a cell cycle arrest in the G0/G1 phase (Fig. 3D).

Fig. 3
figure 3

Restoration of PRELP expression in melanoma cells and its effect on growth properties. Murine and human PRELPlow melanoma cells by transfection. Restoration of PRELP expression. A PRELP expression in PRELPhigh vs. PRELPlow expressing melanoma cells was determined by Western blot analysis as described in “Methods”. A representative WB using an anti-PRELP Ab is shown and the PRELP protein expression (37 kD) is marked. B Reduced proliferation of PRELPhigh vs. PRELPlow expressing melanoma cells. The growth properties of PRELPhigh vs. PRELPlow B16F10 and Buf1088 cells were determined as described in “Methods”. The proliferation rates of PRELPlow and PRELPhigh melanoma cells were correlated to that of PRELPlow vector controls, which were set to “100”. The data are represented in bar charts as the mean of three independent experiments. C Altered migration capacity of PRELPlow vs. PRELPhigh melanoma cells. The migration rate of PRELPlow vs. PRELPhigh B16F10 and Buf1088 cells was determined by ATP-based fluorescence as described in “Methods” and normalized to the seeding control. The graph represents the % of migrated PRELPlow and PRELPhigh melanoma cells as a mean of three independent experiments. D Altered cell cycle in PRELPlow vs. PRELPhigh melanoma cells. Cell cycle distribution of PRELPlow vs. PRELPhigh B16F10 and Buf1088 was assessed by flow cytometry as described in “Methods”. The data were presented as the % of cells in the different cell cycle phases as the mean of three independent experiments demonstrating a cell cycle arrest in the G0/G1 phase of PRELPhigh melanoma cells. The statistical significance is presented as *p < 0.05; **p < 0.01; ***p < 0.001

Furthermore, PRELPhigh B16F10 and human Buf1088 cells showed a strong upregulation of MHC class I surface expression (Fig. 4A), which was accompanied by an increased expression of major components of the MHC class I APM, like the transporter-associated with antigen processing (TAP)1 and TAP2, tapasin (TAPBP), β2-microglobulin (β2-m) and IFN-γ-induced proteasome subunits (Fig. 4B).

Fig. 4
figure 4

Upregulation of MHC class I APM and IFN signaling components and altered NK cell response by restoration of PRELP in PRELPlow/neg B16F10 and Buf1088 cells. A PRELP-mediated upregulation of MHC class I surface antigen expression. MHC class I surface expression was assessed by flow cytometry in PRELPlow/neg vs. PRELPhigh B16F10 and Buf1088 cells as described in “Methods”. The data are represented as histograms of MHC class I surface expression of PRELPlow/neg vs. PRELPhigh B16F10 and Buf1088 cells. B PRELP-mediated upregulation of MHC class I APM components expression in melanoma cells. The mRNA expression levels of MHC class I APM components in PRELPlow/neg vs. PRELPhigh melanoma cells were analyzed by qPCR as described in “Methods”. The results are expressed in bar charts representing the mean of three independent experiments as relative mRNA expression levels of selected APM components in PRELPlow/neg and PRELPhigh B16F10 and Buf1088 cells. C. Transcriptional downregulation of selected MHC class I APM components in PRELPlow/neg melanoma cells. APM promoter activity in PRELPlow/high, mock transfectants and two independent PRELP transfectants of B16F10 cells was determined by luciferase assays as described in “Methods”. The data were normalized to β-gal activity and presented in a bar chart as the mean of the relative luc activity of at least three independent experiments. Error bars indicate the standard error. The statistical significance is presented as *p < 0.05; **p < 0.01

Since the components of the IFN-γ signal transduction and NLRC5 have been linked to MHC class I APM component expression [9, 70], these molecules were also analyzed in PRELPlow vs. PRELPhigh murine and human melanoma cell lines cells. PRELPhigh melanoma cells express increased levels of NLRC5 (Fig. 5A) as well as selected IFN-γ signal transduction molecules, e.g., IRF1 (Fig. 5B), IRF5, STAT1 and STAT2 (Fig. 5C).

Fig. 5
figure 5

Altered expression of IFN-γ signal components, NLRC5 and TGF-β pathway molecules in PRELPhigh vs. PRELPlow/neg melanoma cells. A Increased mRNA levels of IRF1 and NLRC5 in PRELPhigh vs. PRELPlow melanoma cells were assessed by qPCR as described in “Methods” for IRF1 and NLRC5 expression. The qPCR data are expressed in bar charts relative to parental cells (set 1) and represent the mean ± SE format from three independent experiments. B A representative WB using an anti-IRF1 Ab is shown and the PRELP protein expression is marked. C The mRNA levels of IRF5, STAT1 and STAT2 were determined by qPCR as described in “Methods”. The qPCR data are shown in bar charts relative to parental cells (set 1) and represent the mean ± SE format from three independent experiments. D Altered TGF-β signaling components in PRELPhigh B16F10 cells. The expression levels of TGFB1, TRFBR1 and SMAD2 were determined by qPCR. The data are expressed in bar charts relative to parental cells (set 1) and represent the mean ± SE format from three independent experiments. E Reduced NK cell activity in PRELPhigh vs. PRELPlow Buf1088 cells. CD107a degranulation assay was performed as described in “Methods” by co-culture with NK cells from three different donors with PRELPlow vs. PRELPhigh Buf1088 cells. The mean ± SE of the CD107a degranulation of PRELPlow vs. PRELPhigh Buf1088 cells using NK cells representing total NK cell activity are shown. The statistical significance is presented as *p < 0.05; **p < 0.01; ***p < 0.001

Based on the known effect of PRELP on the TGF-β signaling, a possible link between PRELP and MHC class I expression with this pathway was also investigated. As expected, the TGF-β receptor (TGFBR1), its ligand (TGF-β1) and SMAD2 were downregulated in PRELPhigh Buf1088 cells (Fig. 5D). The functional impact of the PRELP-mediated MHC class I upregulation on NK cell responses was analyzed in Buf1088 cells using a CD107 degranulation assay. As shown in Fig. 5E, a reduced NK cell recognition of PRELPhigh compared to PRELPlow Buf1088 cells was found.

Underlying molecular mechanisms of low PRELP expression and PRELP-mediated downregulation of MHC class I in melanoma

To understand the molecular mechanisms of low PRELP expression in melanoma, the frequency of structural alterations in the PRELP gene was determined in 287 melanoma samples of the TCGA dataset (Tumor Skin Cutaneous Melanoma—TCGA) with available mutation and copy number alteration data. As shown in Additional file 1: Fig. S4, mutations only occurred in 20/287 melanoma samples (7%) with 7 missense mutations with unknown significance distributed over the whole gene. Based on these bioinformatics data, it was hypothesized that the low PRELP expression levels in melanoma might be mainly due to deregulation rather than genomic abnormalities.

To determine the PRELP-mediated downregulation of MHC class I APM component expression, the promoter activity of selected APM components was analyzed by luciferase (luc) reporter assays in PRELPlow vs. PRELPhigh melanoma cells. Transient transfection of the APM promoter luc constructs and vector controls in parental PRELPlow/neg B16F10 cells and PRELPhigh B16F10 transfectants demonstrated a strong induction of luc activity in PRELPhigh compared to PRELPlow/neg B16F10 cells, but with differences in the promoter activity levels depending on the APM component analyzed (Fig. 4C).

Link of PRELP expression with immune cell infiltration and CCL5 in melanoma cell lines

In order to determine whether there is a link between high PRELP expressions associated with high HLA class I levels and immune cell infiltration, immune cell infiltration was determined in correlation to PRELP expression. As shown in Fig. 6A, high expression of PRELP in tumors correlated with a high number of tumor-infiltrating CD8+ T cells in skin cutaneous (HR = − 4.079, p = 1.25e−04) and metastatic melanoma (HR = − 3.758; p = 6.19e−04) (Fig. 6B).

Fig. 6
figure 6

Link of CD8+ T cell infiltration levels with PRELP and CCL5 expression in melanoma. A, B Data sets of tumor skin cutaneous melanoma (A) and metastatic melanoma (B) were analyzed by bioinformatics and linked to immune cell infiltration. The prognostic signature of the infiltration of CD8+ T cell subpopulations and the expression of PRELP profiled by xCell algorithm is represented as the Cox model. C A correlation plot was displayed between CCL5 and CD8+ T infiltration in a cohort of 480 melanoma samples of the TCGA dataset. D Secretion of CCL5 in PRELPhigh vs. PRELPlow melanoma cells. Supernatants of PRELPhigh and PRELPlow/negative Buf1088 cells were analyzed for the production of the chemokine CCL5 cells using ELISA as described in “Methods”. Representative data from three independent experiments are shown. The statistical significance is presented as **p < 0.01

In-depth analysis of T cell subpopulations revealed that high PRELP expression levels in melanoma were significantly correlated to high infiltrations of naive (HR = − 2.810; p = 1.36e−02 and HR = − 2.805; p = 1.72e−02), central memory (HR = − 4.326; p = 4.54e−05 and HR = − 3.973; p = 2.71e−04) and effector memory (HR = − 2.807; p = 1.25e−02 and HR = − 2.676; p = 2.3e−02) CD8+ T cells in both tumor skin cutaneous (Additional file 1: Fig. S5A) and metastatic melanoma (Additional file 1: Fig. S5B) datasets.

Comparable to tumor-infiltrating CD8+ T cells, high PRELP expression in SKCM was correlated with an increased number of tumor-infiltrating monocytes (HR = − 3.077; p = 5.87E-03 and HR = − 2.832; p = 1.59E−02), macrophages (HR = − 3.934; p = 3.01E−04 and HR = -3.46; p = 2.23E−03), myeloid dendritic cells (HR = − 2.619; p = 2.30E−02 and HR = − 2.375; p = 5.47E−02) in both skin cutaneous and metastatic melanoma datasets (Additional file 1: Table S4). The number of tumor-infiltrating B cells and tumor PRELP expression were associated with a decreased risk of skin cutaneous melanoma (HR = − 2.533; p = 2.81E−02), which was not significant for metastatic melanoma (HR = − 1.995; p = 1.25E−01). In contrast to the associations described, a high frequency of common myeloid progenitors was associated with an increased risk (HR = 2.284; p = 4.84E−02) in PRELPhigh skin cutaneous melanoma patients. The frequency of other infiltrating cell subpopulations listed in Additional file 1: Table S4, including CD4+ T cell subpopulations, innate immune cells, cells of the hematopoietic lineage and stroma cells, lack significant clinical relevance. Since the chemokine CCL5 is known to recruit CD8+ T lymphocytes to the site of inflammation [71], the frequency of CD8+ T infiltration in tumors was correlated with the CCL5 expression in a cohort of 480 melanoma samples of the TCGA dataset. A strong positive correlation (Rho = 0.733; p = 3.76e−78) between CCL5 and CD8+ T infiltration was found (Fig. 6C). Based on these data, the effect of PRELP overexpression on CCL5 secretion in Buf1088 melanoma cells was determined using the ELISA. When compared to PRELPlow parental cells, PRELPhigh Buf1088 showed higher secretion levels of CCL5 into the cell supernatant (Fig. 6D). This is consistent with the TCGA data of global PRELP (log2 5.7 to 13) (Additional file 1: Fig. S6A) and PRELPhigh (log2 > 10.6) (Additional file 1: Fig. S6B) melanoma with a higher OS (global PRELP; p = 0.00006; log-rank test = 17.72; PRELPhigh; p = 0.002; log-rank test = 9.472) upon high CCL5 expression. In contrast, CCL5 has no prognostic relevance in PRELPlow patients (p = 0.163; log-rank test = 1.950; log2 < 8.1) (Additional file 1: Fig. S6C).

Discussion

PRELP has been shown to affect the immunogenicity of tumors by upregulating MHC class I surface expression and activating IFN signaling, which could modulate tumor development as well as the TME leading to an enhanced CD8+ T cell infiltration. Based on the activity of PRELP on tumor cell function described in this study, this SLRP is suggested as a candidate for anti-cancer therapy. However, so far there exist only a few experimental studies that consider the therapeutic role of this molecule. Since upregulation of PRELP has been demonstrated to be associated with anti-tumoral activity, PRELP should be reconstituted in melanoma, which might be a novel form of anti-cancer strategy. In addition, compounds should be identified, which upregulate PRELP expression leading to anti-cancer activity. Next to its therapeutic potential, our bioinformatics pan-cancer analysis as well as our functional results demonstrated the potential use of PRELP as prognostic and therapeutic markers. Due to a dual role of some SLRPs in cancer, further characterisation of key ECM molecules in tumor progression or tumor suppression is required [72,73,74].

During the last decade, increased information exists about the role and function of ECM components in tumors, which are altered at the biomechanical, biochemical, architectural and topographic level [14]. The ECM has been further shown to be an essential and dynamic part of the TME [75] and undergoes remodeling mediated by several matrix-degrading enzymes during normal and pathologic conditions. A deregulation of the ECM composition and structure has been associated with the development of tumor progression by a detachment of tumor cells from each other, from adjacent immune and stromal cells and contributes to increased cell proliferation, migration and invasion as well as reduced patients’ survival [76]. In addition, ECM components affect the EMT and dissemination of cells, the composition of the TME and tumor immunogenicity [76, 77].

In general, SLRPs were first correlated with the regulation of innate immune responses, thereby generating a pro-inflammatory TME, which could trigger tumorigenesis [78]. However, some SLRPs, an oncogenic and/or tumor-suppressive role has been discussed. For example, a link between the BGN-mediated downregulation, oncogenic transformation and a reduced immunogenicity due to low MHC class I surface expression levels was reported. This was caused by transcriptional suppression of major APM components and an increased TGF-β signaling [18]. In addition, the SLRP decorin also exhibited anti-tumorigenic properties by downregulating the release of IL-10 thereby inhibiting tumor growth [79]. It can also alter the tumor stroma and modulate the tumor-associated inflammation [80].

Next to BGN, the expression of PRELP was also downregulated in K-RAS- and/or HER-2/neu-transformed fibroblasts [44]. Based on these results, a tumor-suppressive activity of PRELP in tumor cells has been suggested. Indeed, our results demonstrated a PRELP-mediated induction of the expression of MHC class I and secretion of CCL5 in melanoma cells (Fig. 7), which might contribute (i) to the recruitment of effector T cells into the tumor tissues and (ii) to an increased recognition by CD8+ T cells.

Fig. 7
figure 7

Schematic diagram of PRELP activity in hallmarks of cancer. The reduced PRELP expression in melanoma cells is associated with low expression levels of MHC class I APM components and members of the IFN-γ signal transduction as well as CCL5, which could be reverted by reconstitution of PRELP. PRELP induces the expression of MHC class I and secretion of CCL5 in melanoma cells, which contributes to the recruitment of effector T cells into the tumor tissues and increased recognition by CD8+ T cells (created with BioRender.com)

This is in line with recent reports demonstrating a tumor inhibitory capacity of PRELP [38, 40], which might be at least partially due to its ability to directly bind TGFB1 [30, 31, 34]. Interestingly, our study demonstrated a reduced expression of TGFB1, TFGBR1 and SMAD2 in PRELPhigh vs PRELPlow melanoma cells.

Furthermore, bioinformatics analyses of high throughput TCGA and GEO data showed a reduced PRELP expression occurring at a high frequency in solid cancers compared to adjacent normal tissues, including primary melanoma lesions, which was even more pronounced in melanoma metastases. Furthermore, low levels of PRELP expression were associated with a low T cell infiltration and low HLA class I expression suggesting that an impaired PRELP expression of tumor cells represents an immune escape phenotype.

This is further in line with a reduced expression of CCL5 in PRELPlow melanoma, which is a potent chemoattractant for T cells. CCL5 plays an important role in recruitment of effector T cells to the tumor site. Furthermore, its expression shows a positive correlation with both increased T cell infiltration and the survival of patients with cancer [81]. Indeed, a positive correlation exists between increased numbers of tumor-infiltrating monocytes and CD8+ T cells and the expression of CCL5 in melanoma, but also in other types of cancer, which further helped the recruitment of CD8+ effector T cells [71, 82, 83]. Thus, the expression of PRELP appears to be linked with the composition of the TME and the cellular immune responses.

Concerning the underlying molecular mechanisms leading to PRELP downregulation, genomic abnormalities in PRELP, such as mutations and deletions, were only found at a low frequency (7%) suggesting deregulation as a major mechanism of impaired PRELP expression in melanoma (Additional file 1: Fig. S3). This could occur at the epigenetic, transcriptional as well as post-transcriptional levels. However, the underlying mechanisms of PRELP expression in 2/48 melanoma cell lines have not been identified. Restoring PRELP expression might have anti-tumoral activity and could be combined with immunotherapies to enhance treatment efficacy.

There exists increasing evidence about the clinical relevance of PRELP as a tumor suppressor in solid tumors and its use as a prognostic marker [30, 38, 84]. In this study, the PRELP expression status was significantly associated with the patients’ prognosis. Melanoma patients with low PRELP expression had a shorter OS, which was not associated with the metastatic status, age and gender. These data were confirmed by in vitro analysis of human and murine melanoma cell lines and were in line with downregulation of PRELP in CRC specimens compared to adjacent normal mucosa [41].

Since proliferation, invasion and metastasis formation are significant features of melanoma, the effect of PRELP overexpression on growth properties was analyzed, which might also lead to identifying novel therapeutic targets controlling disease progression. In accordance with recently published results for hepatocellular carcinoma [38], PRELP overexpression inhibited cell proliferation, migration and invasion of murine and human melanoma cell lines. However, the inhibitory role of PRELP on melanoma growth has to be investigated in more detail and might be associated with the modulation of multiple signal transduction pathways, such as β-catenin or NF-KB signaling [85, 86]. PRELP treatment has been shown to suppress cancer progression by inhibiting the TGF-β and EGF pathways leading to the control of the EMT [30]. Interestingly, an altered TGF-β signaling was also shown in BGN-overexpressing cells, which was associated with an increased SMAD2 expression due to the inhibition of miR-21. This results in an upregulation of MHC class I APM components and an increased immunogenicity of BGNhigh vs. BGNlow cells.

In sum, this study extends the knowledge about the role of SLRPs regarding their tumor suppressive activity [19]. Not only BGN, but also PRELP has a positive immune modulatory potential by increasing the expression of HLA class I APM and IFN-γ signaling components suggesting a tumor suppressive activity of PRELP in melanoma. Vice versa, reduced PRELP expression in melanoma cells is associated with low expression levels of MHC class I APM components and members of the IFN-γ signal transduction as well as CCL5, which could be reconstituted by PRELP overexpression (Fig. 7). Since genetic alterations were not identified, the impaired PRELP expression in melanoma might be mainly due to its deregulation rather than structural abnormalities. In addition, high PRELP expression is significantly correlated with increased patients’ survival and CD8+ T cell infiltration, which might be associated with enhanced T cell responses in melanoma. Therefore, PRELP is a prognostic biomarker and might be a potential (immune) therapeutics for enhancing CD8+ T cell responses, thereby leading to novel approaches for the treatment of melanoma patients.

Availability of data and materials

All data generated or analyzed during this study are included either in this article or in the supplementary information files. The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

Ab:

Antibody

APM:

Antigen processing machinery

BM:

Basement membrane

β2-m:

β2-Microglobulin

CCL5:

Chemokine ligand 5

CPTAC:

Clinical Proteomic Tumor Analysis Consortium

CRC:

Colorectal carcinoma

CTL:

Cytotoxic T lymphocyte

CTLA4:

Cytotoxic T lymphocyte antigen 4

ECM:

Extracellular matrix

EMT:

Epithelial mesenchymal transition

FCS:

Fetal calf serum

FPKM:

Fragments per kilo base per million mapped reads

GADPH:

Glyceride aldehyd-3-phosphate dehydrogenase

GAG:

Glycosaminoglycan

HC:

Heavy chain

HCC:

Hepatocellular carcinoma

HLA:

Human leukocyte antigen

HR:

Hazard ratio

HRP:

Horse reddish peroxidase

ICP:

Immune checkpoint

ICPi:

Immune checkpoint inhibitor

IFN:

Interferon

IRF:

Interferon regulated factor

JAK:

Janus kinase

LILR:

Leukocyte immunoglobulin receptor

LMP:

Low molecular weight protein

LRR:

Leucine-rich repeat

luc:

Luciferase

mAb:

Monoclonal antibody

MFI:

Mean fluorescence intensity

MHC:

Major histocompatibility complex

MITF:

Microphthalmia-associated transcription factor

NF-kB:

Nuclear factor-kappa B

NK:

Natural killer

OS:

Overall survival

PD1:

Programed death receptor 1

PD-L1:

Programed death ligand 1

PG:

Proteoglycan

PRELP:

Purine-arginine-rich and leucine-rich protein

SKCM:

Skin cutaneous melanoma

SLRP:

Small leucine rich proteoglycan

STAT:

Signal transducer and activator of transcription

TAP:

Transporter associated with antigen processing

TAPBP:

Tapasin

TCGA:

The Cancer Genome Atlas

TGF-β:

Transforming growth factor β

TIMER:

Tumor Immune Estimation Resource

TME:

Tumor microenvironment

TPM:

Transcripts per kilo base million

References

  1. Pollack LA, Li J, Berkowitz Z, Weir HK, Wu XC, Ajani UA, et al. Melanoma survival in the United States, 1992 to 2005. J Am Acad Dermatol. 2011;65(5 Suppl 1):S78-86.

    PubMed  PubMed Central  Google Scholar 

  2. Eggermont AMM, Crittenden M, Wargo J. Combination immunotherapy development in melanoma. Am Soc Clin Oncol Educ Book. 2018;38:197–207.

    PubMed  Google Scholar 

  3. Larkin J, Chiarion-Sileni V, Gonzalez R, Grob JJ, Rutkowski P, Lao CD, et al. Five-year survival with combined nivolumab and ipilimumab in advanced melanoma. N Engl J Med. 2019;381(16):1535–46.

    CAS  PubMed  Google Scholar 

  4. Ladanyi A, Papp E, Mohos A, Balatoni T, Liszkay G, Olah J, et al. Role of the anatomic site in the association of HLA class I antigen expression level in metastases with clinical response to ipilimumab therapy in patients with melanoma. J Immunother Cancer. 2020;8(1): e000209.

    PubMed  PubMed Central  Google Scholar 

  5. Sucker A, Zhao F, Pieper N, Heeke C, Maltaner R, Stadtler N, et al. Acquired IFNgamma resistance impairs anti-tumor immunity and gives rise to T-cell-resistant melanoma lesions. Nat Commun. 2017;8:15440.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Donia M, Harbst K, van Buuren M, Kvistborg P, Lindberg MF, Andersen R, et al. Acquired immune resistance follows complete tumor regression without loss of target antigens or IFNgamma signaling. Cancer Res. 2017;77(17):4562–6.

    CAS  PubMed  Google Scholar 

  7. Andersen R, Westergaard MCW, Kjeldsen JW, Muller A, Pedersen NW, Hadrup SR, et al. T-cell responses in the microenvironment of primary renal cell carcinoma-implications for adoptive cell therapy. Cancer Immunol Res. 2018;6(2):222–35.

    CAS  PubMed  Google Scholar 

  8. Seliger B, Ritz U, Abele R, Bock M, Tampe R, Sutter G, et al. Immune escape of melanoma: first evidence of structural alterations in two distinct components of the MHC class I antigen processing pathway. Cancer Res. 2001;61(24):8647–50.

    CAS  PubMed  Google Scholar 

  9. Respa A, Bukur J, Ferrone S, Pawelec G, Zhao Y, Wang E, et al. Association of IFN-gamma signal transduction defects with impaired HLA class I antigen processing in melanoma cell lines. Clin Cancer Res. 2011;17(9):2668–78.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Rodriguez T, Mendez R, Del Campo A, Jimenez P, Aptsiauri N, Garrido F, et al. Distinct mechanisms of loss of IFN-gamma mediated HLA class I inducibility in two melanoma cell lines. BMC Cancer. 2007;7:34.

    PubMed  PubMed Central  Google Scholar 

  11. Sucker A, Zhao F, Real B, Heeke C, Bielefeld N, Mabetaen S, et al. Genetic evolution of T-cell resistance in the course of melanoma progression. Clin Cancer Res. 2014;20(24):6593–604.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Kageshita T, Hirai S, Ono T, Hicklin DJ, Ferrone S. Down-regulation of HLA class I antigen-processing molecules in malignant melanoma: association with disease progression. Am J Pathol. 1999;154(3):745–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Lee JH, Shklovskaya E, Lim SY, Carlino MS, Menzies AM, Stewart A, et al. Transcriptional downregulation of MHC class I and melanoma de-differentiation in resistance to PD-1 inhibition. Nat Commun. 2020;11(1):1897.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Cox TR. The matrix in cancer. Nat Rev Cancer. 2021;21(4):217–38.

    CAS  PubMed  Google Scholar 

  15. Vitale D, Kumar Katakam S, Greve B, Jang B, Oh ES, Alaniz L, et al. Proteoglycans and glycosaminoglycans as regulators of cancer stem cell function and therapeutic resistance. FEBS J. 2019;286(15):2870–82.

    CAS  PubMed  Google Scholar 

  16. Iozzo RV, Schaefer L. Proteoglycan form and function: a comprehensive nomenclature of proteoglycans. Matrix Biol. 2015;42:11–55.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Neill T, Schaefer L, Iozzo RV. Oncosuppressive functions of decorin. Mol Cell Oncol. 2015;2(3): e975645.

    PubMed  PubMed Central  Google Scholar 

  18. Subbarayan K, Leisz S, Wickenhauser C, Bethmann D, Massa C, Steven A, et al. Biglycan-mediated upregulation of MHC class I expression in HER-2/neu-transformed cells. Oncoimmunology. 2018;7(4): e1373233.

    PubMed  PubMed Central  Google Scholar 

  19. Subbarayan K, Massa C, Lazaridou MF, Ulagappan K, Seliger B. Identification of a novel miR-21-3p/TGF-beta signaling-driven immune escape via the MHC class I/biglycan axis in tumor cells. Clin Transl Med. 2021;11(3): e306.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Subbarayan K, Massa C, Leisz S, Steven A, Bethmann D, Biehl K, et al. Biglycan as a potential regulator of tumorgenicity and immunogenicity in K-RAS-transformed cells. Oncoimmunology. 2022;11(1):2069214.

    PubMed  PubMed Central  Google Scholar 

  21. Matsushima N, Ohyanagi T, Tanaka T, Kretsinger RH. Super-motifs and evolution of tandem leucine-rich repeats within the small proteoglycans–biglycan, decorin, lumican, fibromodulin, PRELP, keratocan, osteoadherin, epiphycan, and osteoglycin. Proteins. 2000;38(2):210–25.

    CAS  PubMed  Google Scholar 

  22. Bengtsson E, Morgelin M, Sasaki T, Timpl R, Heinegard D, Aspberg A. The leucine-rich repeat protein PRELP binds perlecan and collagens and may function as a basement membrane anchor. J Biol Chem. 2002;277(17):15061–8.

    CAS  PubMed  Google Scholar 

  23. Bengtsson E, Neame PJ, Heinegard D, Sommarin Y. The primary structure of a basic leucine-rich repeat protein, PRELP, found in connective tissues. J Biol Chem. 1995;270(43):25639–44.

    CAS  PubMed  Google Scholar 

  24. Pietraszek-Gremplewicz K, Karamanou K, Niang A, Dauchez M, Belloy N, Maquart FX, et al. Small leucine-rich proteoglycans and matrix metalloproteinase-14: key partners? Matrix Biol. 2019;75–76:271–85.

    PubMed  Google Scholar 

  25. McEwan PA, Scott PG, Bishop PN, Bella J. Structural correlations in the family of small leucine-rich repeat proteins and proteoglycans. J Struct Biol. 2006;155(2):294–305.

    CAS  PubMed  Google Scholar 

  26. Matsushima N, Takatsuka S, Miyashita H, Kretsinger RH. Leucine rich repeat proteins: sequences, mutations, structures and diseases. Protein Pept Lett. 2019;26(2):108–31.

    CAS  PubMed  Google Scholar 

  27. Matsushima N, Miyashita H, Kretsinger RH. Sequence features, structure, ligand interaction, and diseases in small leucine rich repeat proteoglycans. J Cell Commun Signal. 2021;15:519–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Happonen KE, Furst CM, Saxne T, Heinegard D, Blom AM. PRELP protein inhibits the formation of the complement membrane attack complex. J Biol Chem. 2012;287(11):8092–100.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Kosuge H, Nakakido M, Nagatoishi S, Fukuda T, Bando Y, Ohnuma SI, et al. Proteomic identification and validation of novel interactions of the putative tumor suppressor PRELP with membrane proteins including IGFI-R and p75NTR. J Biol Chem. 2021;296: 100278.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Papadaki V, Asada K, Watson JK, Tamura T, Leung A, Hopkins J, et al. Two secreted proteoglycans, activators of urothelial cell–cell adhesion, negatively contribute to bladder cancer initiation and progression. Cancers. 2020;12(11):3362.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Chacon-Solano E, Leon C, Carretero M, Garcia M, Sanchez-Dominguez R, Quero F, et al. Mechanistic interrogation of mutation-independent disease modulators of RDEB identifies the small leucine-rich proteoglycan PRELP as a TGF-beta antagonist and inhibitor of fibrosis. Matrix Biol. 2022;111:189–206.

    CAS  PubMed  Google Scholar 

  32. Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.

    PubMed  Google Scholar 

  33. Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, et al. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia. 2004;6(1):1–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Chiavarina B, Ronca R, Otaka Y, Sutton RB, Rezzola S, Yokobori T, et al. Fibroblast-derived prolargin is a tumor suppressor in hepatocellular carcinoma. Oncogene. 2022;41(10):1410–20.

    CAS  PubMed  Google Scholar 

  35. Sethi MK, Downs M, Shao C, Hackett WE, Phillips JJ, Zaia J. In-depth matrisome and glycoproteomic analysis of human brain glioblastoma versus control tissue. Mol Cell Proteom. 2022;21: 100216.

    CAS  Google Scholar 

  36. Chen R, Dawson DW, Pan S, Ottenhof NA, de Wilde RF, Wolfgang CL, et al. Proteins associated with pancreatic cancer survival in patients with resectable pancreatic ductal adenocarcinoma. Lab Invest. 2015;95(1):43–55.

    CAS  PubMed  Google Scholar 

  37. Castells X, Acebes JJ, Boluda S, Moreno-Torres A, Pujol J, Julia-Sape M, et al. Development of a predictor for human brain tumors based on gene expression values obtained from two types of microarray technologies. OMICS. 2010;14(2):157–64.

    CAS  PubMed  Google Scholar 

  38. Hong R, Gu J, Niu G, Hu Z, Zhang X, Song T, et al. PRELP has prognostic value and regulates cell proliferation and migration in hepatocellular carcinoma. J Cancer. 2020;11(21):6376–89.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Ning X, Deng Y. Identification of key pathways and genes influencing prognosis in bladder urothelial carcinoma. Onco Targets Ther. 2017;10:1673–86.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Iuga C, Seicean A, Iancu C, Buiga R, Sappa PK, Volker U, et al. Proteomic identification of potential prognostic biomarkers in resectable pancreatic ductal adenocarcinoma. Proteomics. 2014;14(7–8):945–55.

    CAS  PubMed  Google Scholar 

  41. Liu X, Xu Y, Meng Q, Zheng Q, Wu J, Wang C, et al. Proteomic analysis of minute amount of colonic biopsies by enteroscopy sampling. Biochem Biophys Res Commun. 2016;476(4):286–92.

    CAS  PubMed  Google Scholar 

  42. Rucci N, Capulli M, Ventura L, Angelucci A, Peruzzi B, Tillgren V, et al. Proline/arginine-rich end leucine-rich repeat protein N-terminus is a novel osteoclast antagonist that counteracts bone loss. J Bone Miner Res. 2013;28(9):1912–24.

    CAS  PubMed  Google Scholar 

  43. Mikaelsson E, Osterborg A, Jeddi-Tehrani M, Kokhaei P, Ostadkarampour M, Hadavi R, et al. A proline/arginine-rich end leucine-rich repeat protein (PRELP) variant is uniquely expressed in chronic lymphocytic leukemia cells. PLoS ONE. 2013;8(6): e67601.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Recktenwald CV, Mendler S, Lichtenfels R, Kellner R, Seliger B. Influence of Ki-ras-driven oncogenic transformation on the protein network of murine fibroblasts. Proteomics. 2007;7(3):385–98.

    CAS  PubMed  Google Scholar 

  45. Wulfanger J, Biehl K, Tetzner A, Wild P, Ikenberg K, Meyer S, et al. Heterogeneous expression and functional relevance of the ubiquitin carboxyl-terminal hydrolase L1 in melanoma. Int J Cancer. 2013;133(11):2522–32.

    PubMed  Google Scholar 

  46. Rodriguez LG, Wu X, Guan JL. Wound-healing assay. Methods Mol Biol. 2005;294:23–9.

    PubMed  Google Scholar 

  47. Leisz S, Schulz K, Erb S, Oefner P, Dettmer K, Mougiakakos D, et al. Distinct von Hippel-Lindau gene and hypoxia-regulated alterations in gene and protein expression patterns of renal cell carcinoma and their effects on metabolism. Oncotarget. 2015;6(13):11395–406.

    PubMed  PubMed Central  Google Scholar 

  48. Lennicke C, Rahn J, Bukur J, Hochgrafe F, Wessjohann LA, Lichtenfels R, et al. Modulation of MHC class I surface expression in B16F10 melanoma cells by methylseleninic acid. Oncoimmunology. 2017;6(6): e1259049.

    PubMed  Google Scholar 

  49. Stehle F, Schulz K, Fahldieck C, Kalich J, Lichtenfels R, Riemann D, et al. Reduced immunosuppressive properties of axitinib in comparison with other tyrosine kinase inhibitors. J Biol Chem. 2013;288(23):16334–47.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Massa C, Seliger B. Fast dendritic cells stimulated with alternative maturation mixtures induce polyfunctional and long-lasting activation of innate and adaptive effector cells with tumor-killing capabilities. J Immunol. 2013;190(7):3328–37.

    CAS  PubMed  Google Scholar 

  51. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 2013;41(Database issue):D991–5.

    CAS  PubMed  Google Scholar 

  52. Park SJ, Yoon BH, Kim SK, Kim SY. GENT2: an updated gene expression database for normal and tumor tissues. BMC Med Genom. 2019;12(Suppl 5):101.

    Google Scholar 

  53. Zheng G, Ma Y, Zou Y, Yin A, Li W, Dong D. HCMDB: the human cancer metastasis database. Nucleic Acids Res. 2018;46(D1):D950–5.

    CAS  PubMed  Google Scholar 

  54. Chen F, Chandrashekar DS, Varambally S, Creighton CJ. Pan-cancer molecular subtypes revealed by mass-spectrometry-based proteomic characterization of more than 500 human cancers. Nat Commun. 2019;10(1):5679.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.

    PubMed  PubMed Central  Google Scholar 

  56. Raskin L, Fullen DR, Giordano TJ, Thomas DG, Frohm ML, Cha KB, et al. Transcriptome profiling identifies HMGA2 as a biomarker of melanoma progression and prognosis. J Invest Dermatol. 2013;133(11):2585–92.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Xu L, Shen SS, Hoshida Y, Subramanian A, Ross K, Brunet JP, et al. Gene expression changes in an animal melanoma model correlate with aggressiveness of human melanoma metastases. Mol Cancer Res. 2008;6(5):760–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. Jerby-Arnon L, Shah P, Cuoco MS, Rodman C, Su MJ, Melms JC, et al. A cancer cell program promotes T cell exclusion and resistance to checkpoint blockade. Cell. 2018;175(4):984-997.e24.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Cancer Genome Atlas N. Genomic classification of cutaneous melanoma. Cell. 2015;161(7):1681–96.

    Google Scholar 

  60. Toma G, Lemnian IM, Karapetian E, Grosse I, Seliger B. Transcriptional analysis of total CD8(+) T cells and CD8(+)CD45RA(−) memory t cells from young and old healthy blood donors. Front Immunol. 2022;13: 806906.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. Qi H, Song S, Wang P. ImmuMethy, a database of DNA methylation plasticity at a single cytosine resolution in human blood and immune cells. Database. 2022. https://doi.org/10.1093/database/baac020.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–14.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18(1):220.

    PubMed  PubMed Central  Google Scholar 

  64. Diehl V, Huber LS, Trebicka J, Wygrecka M, Iozzo RV, Schaefer L. The role of decorin and biglycan signaling in tumorigenesis. Front Oncol. 2021;11: 801801.

    CAS  PubMed  PubMed Central  Google Scholar 

  65. Hu X, Villodre ES, Larson R, Rahal OM, Wang X, Gong Y, et al. Decorin-mediated suppression of tumorigenesis, invasion, and metastasis in inflammatory breast cancer. Commun Biol. 2021;4(1):72.

    CAS  PubMed  PubMed Central  Google Scholar 

  66. Subbarayan K, Seliger B. Tumor-dependent effects of proteoglycans and various glycosaminoglycan synthesizing enzymes and sulfotransferases on patients’ outcome. Curr Cancer Drug Targets. 2019;19(3):210–21.

    CAS  PubMed  Google Scholar 

  67. Hopkins J, Asada K, Leung A, Papadaki V, Davaapil H, Morrison M, et al. PRELP regulates cell–cell adhesion and EMT and inhibits retinoblastoma progression. Cancers. 2022;14(19):4926.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. Badal B, Solovyov A, Di Cecilia S, Chan JM, Chang LW, Iqbal R, et al. Transcriptional dissection of melanoma identifies a high-risk subtype underlying TP53 family genes and epigenome deregulation. JCI Insight. 2017;2(9): e92102.

    PubMed  PubMed Central  Google Scholar 

  69. Riker AI, Enkemann SA, Fodstad O, Liu S, Ren S, Morris C, et al. The gene expression profiles of primary and metastatic melanoma yields a transition point of tumor progression and metastasis. BMC Med Genom. 2008;1:13.

    Google Scholar 

  70. Meissner TB, Li A, Biswas A, Lee KH, Liu YJ, Bayir E, et al. NLR family member NLRC5 is a transcriptional regulator of MHC class I genes. Proc Natl Acad Sci USA. 2010;107(31):13794–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  71. Liu J, Li F, Ping Y, Wang L, Chen X, Wang D, et al. Local production of the chemokines CCL5 and CXCL10 attracts CD8+ T lymphocytes into esophageal squamous cell carcinoma. Oncotarget. 2015;6(28):24978–89.

    PubMed  PubMed Central  Google Scholar 

  72. Appunni S, Rubens M, Ramamoorthy V, Anand V, Khandelwal M, Sharma A. Biglycan: an emerging small leucine-rich proteoglycan (SLRP) marker and its clinicopathological significance. Mol Cell Biochem. 2021;476(11):3935–50.

    CAS  PubMed  Google Scholar 

  73. Berdiaki A, Giatagana EM, Tzanakakis G, Nikitovic D. The landscape of small leucine-rich proteoglycan impact on cancer pathogenesis with a focus on biglycan and lumican. Cancers. 2023;15(14):3549.

    CAS  PubMed  PubMed Central  Google Scholar 

  74. Neill T, Schaefer L, Iozzo RV. Decoding the matrix: instructive roles of proteoglycan receptors. Biochemistry. 2015;54(30):4583–98.

    CAS  PubMed  Google Scholar 

  75. Theocharis AD, Skandalis SS, Gialeli C, Karamanos NK. Extracellular matrix structure. Adv Drug Deliv Rev. 2016;97:4–27.

    CAS  PubMed  Google Scholar 

  76. Gerarduzzi C, Hartmann U, Leask A, Drobetsky E. The matrix revolution: matricellular proteins and restructuring of the cancer microenvironment. Cancer Res. 2020;80(13):2705–17.

    CAS  PubMed  Google Scholar 

  77. Walker C, Mojares E, Del Rio Hernandez A. Role of extracellular matrix in development and cancer progression. Int J Mol Sci. 2018;19(10):3028.

    PubMed  PubMed Central  Google Scholar 

  78. Tzanakakis G, Neagu M, Tsatsakis A, Nikitovic D. Proteoglycans and immunobiology of cancer-therapeutic implications. Front Immunol. 2019;10:875.

    CAS  PubMed  PubMed Central  Google Scholar 

  79. Merline R, Moreth K, Beckmann J, Nastase MV, Zeng-Brouwers J, Tralhao JG, et al. Signaling by the matrix proteoglycan decorin controls inflammation and cancer through PDCD4 and microRNA-21. Sci Signal. 2011;4(199):ra75.

    PubMed  PubMed Central  Google Scholar 

  80. Gubbiotti MA, Vallet SD, Ricard-Blum S, Iozzo RV. Decorin interacting network: a comprehensive analysis of decorin-binding partners and their versatile functions. Matrix Biol. 2016;55:7–21.

    CAS  PubMed  PubMed Central  Google Scholar 

  81. Hong M, Puaux AL, Huang C, Loumagne L, Tow C, Mackay C, et al. Chemotherapy induces intratumoral expression of chemokines in cutaneous melanoma, favoring T-cell infiltration and tumor control. Cancer Res. 2011;71(22):6997–7009.

    CAS  PubMed  Google Scholar 

  82. Harlin H, Meng Y, Peterson AC, Zha Y, Tretiakova M, Slingluff C, et al. Chemokine expression in melanoma metastases associated with CD8+ T-cell recruitment. Cancer Res. 2009;69(7):3077–85.

    CAS  PubMed  Google Scholar 

  83. Uehara J, Ohkuri T, Kosaka A, Ishibashi K, Hirata Y, Ohara K, et al. Intratumoral injection of IFN-beta induces chemokine production in melanoma and augments the therapeutic efficacy of anti-PD-L1 mAb. Biochem Biophys Res Commun. 2017;490(2):521–7.

    CAS  PubMed  Google Scholar 

  84. Dozen A, Shozu K, Shinkai N, Ikawa N, Aoyama R, Machino H, et al. Tumor suppressive role of the PRELP gene in ovarian clear cell carcinoma. J Pers Med. 2022;12(12):1999.

    PubMed  PubMed Central  Google Scholar 

  85. Li H, Cui Y, Luan J, Zhang X, Li C, Zhou X, et al. PRELP (proline/arginine-rich end leucine-rich repeat protein) promotes osteoblastic differentiation of preosteoblastic MC3T3-E1 cells by regulating the beta-catenin pathway. Biochem Biophys Res Commun. 2016;470(3):558–62.

    CAS  PubMed  Google Scholar 

  86. Rucci N, Rufo A, Alamanou M, Capulli M, Del Fattore A, Ahrman E, et al. The glycosaminoglycan-binding domain of PRELP acts as a cell type-specific NF-kappaB inhibitor that impairs osteoclastogenesis. J Cell Biol. 2009;187(5):669–83.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank Maria Heise for excellent secretarial help.

Funding

Open Access funding enabled and organized by Projekt DEAL. This study was supported by grants from the Wilhelm Sander-Stiftung (BS, 2019.076.1) and an intramural grant of the Martin Luther University Halle-Wittenberg, Medical Faculty (HS).

Author information

Authors and Affiliations

Authors

Contributions

HS was performing most of the experiments together with KS. CM performed the NK cell assay, AM and CV made the Western blot analysis. All four authors prepared the figures for the manuscript. BS designed the project and the experiments, was responsible for the coordination of the experimental work and wrote the manuscript. All authors discussed and corrected the manuscript.

Corresponding author

Correspondence to Barbara Seliger.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors read the final version of the manuscript and gave consent for publication.

Competing interests

There are no competing interests. The authors have no relevant financial or non-financial interests to disclose.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Figure S1.

PRELP expression in neoplastic and non-neoplastic human tissues and cells. Figure S2. Correlation plot of HLA-A and PRELP in 133 skin cancer cell lines (GENT2—skin cancer cell lines). Figure S3. Reconstitution of PRELP expression in melanoma cells. Overexpression of PRELP in melanoma cell lines was obtained after transfection of a PRELP expression vector in murine and human PRELPlow melanoma cells. Transfection with a mock vector served as control. Figure S4. Distribution of gene mutations of PRELP in melanoma. The Skin Cutaneous Melanoma (TCGA) dataset was analysed for genetic alterations in PRELP. Figure S5. Association between PRELP expression, CD8+ T cell infiltration and activation as well as patients’ survival. Figure S6. A–C. Correlation of PRELPhigh vs. PRELPlow samples with the expression of CCL5 and overall survival. Table S1. Primers used for qPCR analyses. Table S2. Datasets used and number of samples analyzed. Table S3. Correlation of PRELP expression with HLA class I and APM components using two different melanoma datasets. Table S4. Exploring the association between PRELP and immune infiltrates with the clinical outcome in SKCM datasets.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schäfer, H., Subbarayan, K., Massa, C. et al. Correlation of the tumor escape phenotype with loss of PRELP expression in melanoma. J Transl Med 21, 643 (2023). https://doi.org/10.1186/s12967-023-04476-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-023-04476-x

Keywords