- Research article
- Open Access
- Open Peer Review
Diagnostic role of circulating extracellular matrix-related proteins in non-small cell lung cancer
BMC Cancer volume 18, Article number: 899 (2018)
Interactions between cancer cells and the surrounding microenvironment are crucial determinants of cancer progression. During this process, bi-directional communication among tumor cells and cancer associated fibroblasts (CAF) regulate extracellular matrix (ECM) deposition and remodeling. As a result of this dynamic process, soluble ECM proteins can be released into the bloodstream and may represent novel circulating biomarkers useful for cancer diagnosis. The aim of the present study was to measure the levels of three circulating ECM related proteins (COL11A1, COL10A1 and SPARC) in plasma samples of lung cancer patients and in healthy heavy-smokers controls and test whether such measurements have diagnostic or prognostic value.
Gene expression profiling of lung fibroblasts isolated from paired normal and cancer tissue of NSCLC patients was performed by gene expression microarrays. The prioritization of the candidates for the study of circulating proteins in plasma was based on the most differentially expressed genes in cancer associated fibroblasts. Soluble ECM proteins were assessed by western blot in the conditioned medium of lung fibroblasts and by ELISA assays in plasma samples.
Plasma samples from lung cancer patients and healthy heavy-smokers controls were tested for levels of COL11A1 and COL10A1 (n = 57 each) and SPARC (n = 90 each). Higher plasma levels of COL10A1 were detected in patients (p ≤ 0.001), a difference that was driven specifically by females (p < 0.001). No difference in COL11A1 levels between patients and controls was found. SPARC levels were also higher in plasma patients than controls (p < 0.001) with good performance in discriminating the two groups (AUC = 0.744). No significant association was observed between plasma proteins levels and clinicopathological features or survival.
Soluble factors related to proficient tumor-stroma cross-talk are detectable in plasma of primary lung cancer patients and may represent a valuable complementary diagnostic tool to discriminate lung cancer patients from healthy heavy-smokers individuals as shown for the SPARC protein.
Lung cancer represents a major health issue worldwide due to its high incidence and mortality rates . The disease is often diagnosed at advanced stages when current clinical options are largely ineffective resulting five years survival rates of less than 10% . Identification of circulating molecular biomarkers is therefore critical to improve early detection of lung cancer and represents an important approach with a large clinical potential . Despite the recent improvements in the characterization of different circulating biomarkers such as cell free DNA [4, 5], circulating tumor cells [6,7,8], extracellular vesicles or circulating miRNA [9, 10], the search for optimal biomarkers still remains a challenge. In the past, most studies were centered on identifying potential biomarkers using molecules differentially expressed by tumor cells. However, in the last few years, the concept that each tumor is a complex system composed by both, cancer cells and the surrounding stroma, represented by a variety of cell types such as fibroblasts, immune and endothelial has been well established . This notion has important implications for biomarkers research as novel candidates with diagnostic or prognostic value can potentially be obtained by analyzing molecules produced by stromal cells during their interactions with cancer cells [12, 13].
In particular, activated fibroblasts play a prominent role in lung carcinogenesis due to their abilities to trigger several signaling pathways implicated in tumor formation and metastasis [14,15,16]. It has also been demonstrated that cancer-associated fibroblasts (CAFs) within the reactive stroma are responsible for deposition of elevated amounts of extracellular matrix (ECM) . In physiological condition, ECM provides mechanical and biochemical support to the surrounding cells and is actively involved in cell proliferation and migration . On the other hand, under pathological conditions the interactions between activated fibroblasts and epithelial cells result in production of different growth factors, cytokines and proteases which modify the surrounding ECM, by changing its composition and facilitating pathological alterations such as chronic obstructive pulmonary disease (COPD) or idiopathic pulmonary fibrosis (IPF)  and potentially leading to cancer development [20, 21].
As a consequence of this remodeling, proteins related to ECM are released into blood and could be considered as potential novel circulating biomarkers . Therefore, new biomarkers for early detection of lung cancer derived from remodeling of ECM could be developed by identifying candidate genes and pathways from gene expression profiling of cancer associated fibroblasts (CAF). A recent review started to unravel the key pathways involved in their functional effects highlighting the existence of common mechanisms as well as specificities in different cancer types (breast, prostate and lung cancer) . The results showed that most of the commonly enriched gene sets characterizing tumor-promoting fibroblasts were related to structural ECM molecules and ECM organization (e.g. collagens) and in particular, to the ECM3 signature, an ECM-based signature already found associated with bad prognosis in aggressive breast carcinomas (i.e. grade III) . Mechanistically the deposition of collagens in the surrounding tumor influences cancer cells behavior promoting cancer progression and invasion in several cancer types [25, 26]. However, the utility of any specific collagen fragment as plasma circulating biomarker in lung cancer remains unproven. Several studies have also pointed out the importance of another lung microenvironment-related protein, secreted protein acidic and rich in cysteine (SPARC), a collagen-binding matricellular protein considered as a key player in the tumor progression most likely by supporting crosstalk at the tumor–stroma interface [27, 28]. SPARC has been predominantly detected in the tumor-associated stroma, specifically in ECM produced by activated fibroblasts. Interestingly, the localization of SPARC in NSCLC tissues is linked to disease prognosis. In fact, high levels of SPARC expression within NSCLC tumor tissues, are associated to longer survival, while its absence represents a negative prognostic factor. On the other hand, high expression of SPARC in the stroma is associated with poor overall survival in lung cancer patients . Since changes in the ECM can occur early in cancer progression, in this study we aimed to identify plasma circulating proteins originated by the ECM compartment and to investigate their potential utility as biomarkers for early diagnosis of lung cancer. In addition, we aimed to explore the association of ECM proteins with different clinical parameters and their potential prognostic value.
Patient characteristics and tissue sampling
The current study was approved by the Fondazione IRCCS, Istituto Nazionale dei Tumori Ethics Review Board and included all consecutive patients from whom plasma samples were available and who underwent a complete anatomical resection for primary lung cancer at the Thoracic Surgery Division of the National Cancer Institute of Milan, from January 2012 to July 2014. Healthy heavy smoker controls were enrolled in a lung cancer screening program (clinicaltrial.gov NCT 02247453, www.biomild.org) from January 2013 to January 2016. Written informed consent was obtained from all patients and healthy heavy smokers controls for blood collection. All cases used in this study were confirmed to be primary lung cancer by pathology review. Study participants were mainly heavy smokers (12 non-smokers for the analysis of COL11A1 and COL10A1) and were matched 1:1 to the patient cohorts according to sex and age classes (< 50, 50–54, 55–59, 60–64, 65–69, 70–75, > 75) for the analysis of COL11A1 and COL10A1 and according to sex, smoking history and age classes for the analysis of SPARC. Overall survival was the study outcome of interest, thus patients contributed with their time interval from surgery until the date of death or until 16th January 2017 for survivors. Blood collection was performed shortly before surgery to avoid the impact of surgery in the markers quantification. Plasma extraction was described elsewhere . Briefly, whole blood samples (5–10 ml) were collected as first blood with spray-coated K2EDTA tubes (BD-Becton, Dickinson and Company, Plymouth, UK). Within 2 h, plasma was separated by a first centrifugation step at 2500 RPM at 4 °C for 10 min. The supernatant containing plasma was carefully collected avoiding the fraction closest to the lymphocytic ring. Plasma was then centrifuged a second time at 2500 RPM at 4 °C for 10 min. and the supernatant collected and stored at − 80 C until further.
Establishment of cell cultures and conditioned medium
Cultures of primary cancer-associated (CAF) and fibroblasts derived from normal counterpart of (NF) lung cancer patients were isolated from surgical specimens and cultured as already described . All cell lines were routinely tested to exclude presence of mycoplasma contamination, grown as adherent monolayer and harvested at controlled density. To obtain conditioned medium (CM), cells were grown in controlled conditions in serum free medium at the same density (cells number = 1X106). After 24 h, the CM was collected, centrifuged to eliminate cell debris and stored at − 80 C until further.
RNA purification, microarray and data analysis
Total RNA was extracted from fibroblasts cell cultures using RNA easy kit (Qiagen), followed by a clean-up treatment to remove genomic DNA. RNA purity was assessed with bioanalyzer (Agilent technologies) and concentration of RNA was evaluated by nanodrop 2000c (Thermo Scientific). Each microarray experiment was performed using 300 ng of total RNA. Procedures included first strand synthesis, second strand synthesis, double-strand cDNA clean up, in vitro transcription, cRNA purification and fragmentation. One microgram of biotinylated cRNA were finally applied to each hybridization array, Illumina Human HT-12v4 Expression BeadChip (Illumina, Inc., San Diego, CA, USA) at 58 °C for 18 h. Illumina BeadStudio software version 3.8 was used to obtain the raw data. Class comparison analysis was performed using the limma Bioconductor package . Cancer associated fibroblast were compared with normal fibroblasts and all genes were ranked according to the modified t-statistics values obtained. These ranked gene lists were subjected to a Gene Set Enrichment Analysis (GSEA, v.4.0) to identify Gene Ontology terms or Canonical Pathways (BIOCARTA, KEGG, REACTOME) significantly enriched. Enrichment was considered significant at p-value < 0.05.
Cell lysates, extracellular matrix and conditioned medium preparation
For cell lysates preparation, cell lines were solubilized for 1 h on ice with TNTG lysis buffer containing 50 mM Tris-HCl pH 7.5, 150 mM NaCl, 100 mM NaF, 10 mM sodium pyrophosphate, 10% glycerol, 1% Triton X-100 and protease inhibitor cocktail (Complete Mini, Roche, Basel, Switzerland). For extracellular matrix isolation, cultured cells were treated with the hypotonic buffer NH4OH 20 mM for 20 min. After two washes with Phosphate-buffered saline (PBS) 1X the extracellular matrix present on the plastic plate was recovered with heated loading buffer (Laemmli solution) with the help of a scraper. Protein levels in cell-derived extracellular matrix were normalized with respect to the number of cells seeded and grown in the same conditions of cells used for extracellular matrix recovering. Conditioned media (CMs) were processed with Amicon Ultra-15 Centrifugal Filter Unit with Ultracel-3 membrane (Merk Millipore, Billerica, MA, USA) for concentration of proteins with a molecular mass greater than 5 kDa, according to the manufacturer’s instructions. Concentration factors ranged from 30 to 40X.
Proteins analysis techniques
Western Blot. Protein lysates (20 μg), a fixed volume of solubilized extracellular matrix (20 μl) or a fixed volume of concentrated CM (26 μl) were mixed with loading buffer under reducing conditions, heated for 5 min at 95 °C, loaded on 4–12% precast NuPage SDS-Bis-Tris gels (Life Technologies, Carlsbad, CA, USA). The proteins were then transferred to PVDF membranes (Merk Millipore), stained with Red Ponceau to check loading and membranes saturated for 1 h at room temperature in blocking solution (5% low-fat milk in TBS + 0.1% Tween-20) before probing with the appropriate antibodies. Blots were washed with TBS-0.5% Tween-20 and further incubated with horseradish peroxidase-conjugated secondary antibodies (GE Healthcare, Little Chalfont, UK) for 1 h at room temperature. Western blots were developed using the enhanced chemiluminescence method (GE Healthcare) according to the manufacturer’s instructions. Data were acquired and analyzed using Quantity One 4.6.6 software (Bio-Rad, Hercules, CA, USA). The following primary antibodies were used: COL11A1 (1:500 rabbit polyclonal, NBP1–55803, Novus Biologicals, Littleton, CO USA), COL10A1 (1:1000 rabbit polyclonal, LS-C157654, LSBio LifeSpan Biosciences, Seattle, WA USA), SPARC (5 μg/ml mouse IgG1, 33–5500 Invitrogen), and vinculin (1:1000 mouse monoclonal, hVIN-1 clone, Sigma Aldrich). Vinculin quantification, concentration factor and cell counts were used to normalize ECM proteins in cell lysates, CMs and extracellular matrix preparation, respectively.
Circulating proteins were measured by using commercially available ELISA kits (COL10A1 and COL11A1 LifeSpan Biosciences, Inc., SPARC R&D), according to manufacturer’s instruction. Duplicate measures were performed for each sample. Protein levels were expressed as OD value as measured by Microplate Reader Tecan Infinite® M1000.
After matching, the analysis of the levels of ECM molecules was performed on a set of 57 lung cancer patients and 57 healthy controls for COL11A1 and COL10A1 and on a set of 90 lung cancer patients and 90 healthy controls for SPARC. Raw absorbance values were corrected by exploiting the values of the ELISA standards. The distributions of the absorbance values of ECM molecules in plasma samples of lung cancer patients and healthy controls were compared by using the Wilcoxon test. The association between each molecule levels and clinicopathological variables was investigated using the Wilcoxon test for categorical variables and the Spearman correlation coefficient for continuous variables. A univariable logistic regression model including the molecule was fitted for the comparison between patients and healthy controls, and the area under the ROC curve (AUC) was estimated as a measure of discriminative ability Sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) corresponding to the cutoff of the ROC maximizing the Youden Index, were extrapolated. AUC, sensitivity/specificity, PPV/NPV 95% confidence intervals were estimated using bootstrap procedure (Efron B. An Introduction to the Bootstrap. Chapman and Hall, New York, NY; 1993). Multivariable quantile regression models with median as the reference quantile were implemented to study the association between the molecules and the disease status (tumor vs control), adjusting for possible confounders such as age, packyears, COPD and sex. Moreover, we studied the prognostic effect of ECM molecule levels on overall survival (OS). At univariable analyses each molecule was categorized according to its tertiles and the Kaplan-Meier curves were estimated and statistically compared by means of the log-rank test. Multivariable Cox models were also implemented, where the effect of the molecule was adjusted for age, packyears, COPD (present vs absent) and sex (M vs F); the corresponding Hazard Ratios (HR) were estimated. In all the models the molecules were included as continuous variables using 3-knots restricted cubic splines  The analyses were carried out using R software, version 3.2.0 (http://www.r-project.org/). The test results were considered statistically significant whenever a two-sided p-value below 0.05 was obtained.
ECM-related gene expression profiles are enriched in lung cancer associated fibroblasts
We reasoned that novel potential circulating biomarkers derived from the tumor microenvironment could be developed starting from the discovery of biological pathways defining activated stroma. To identify factors responsible for a proficient cross-talk between fibroblasts and cancer cells in lung tumors, we have previously characterized cancer-associated fibroblasts (CAF) and normal fibroblasts (NF) isolated from patients with lung cancer (n = 60) by gene expression profiles. Gene Set Enrichment Analysis (GSEA) was used to identify groups of genes whose expression was specifically altered in cancer associated fibroblasts . Class comparison analysis identified ECM deposition and remodeling among the most enriched pathways in CAF as demonstrated by the up-regulation of two important collagen isoforms, COL11A1 and COL10A1. Interestingly these two proteins also belong to a gene signature with prognostic value so-called ECM3, originally derived by studying expression changes of extracellular matrix genes in a cohort of breast cancer clinical specimens . The ECM3 signature consists of 58 genes encoding 43 structural ECM proteins we found that 11 out the 58 genes were significantly differentially expressed between CAF and NF, each contributing to an enrichment of the ECM3 signature in the activated fibroblasts phenotype (Table 1). To understand the potential of ECM3-related genes as circulating biomarkers, we have focused on three individual proteins, COL11A1 and COL10A1 that resulted the most differentially expressed genes in CAF, and SPARC which is also present in the ECM3 signature but not significantly overexpressed in CAF. The prioritization of these candidates was based on fibroblasts gene expression profiling data for COL11A1 and COL10A1 and on the existing literature on SPARC in lung tumorigenesis .
ECM-related proteins are enriched in CAF and selectively released in conditioned medium
To validate the ECM-related genes expression patterns at protein level we performed western blot analyses on cell lysates of CAF (n = 3) and NF (n = 3) cell lines established from lung cancer patients and found that COL10A1 and SPARC were slightly enriched in CAF compared to NF (p = 0.171 and p = 0.025 respectively), while no differences in COL11A1 expression levels were detected (p = 0.374) (Fig. 1a).
To evaluate whether fibroblasts could represent the source of circulating ECM-related proteins and whether these factors are preferentially released in the medium conditioned by CAF, we analyzed by western blot conditioned medium obtained in serum-free conditions from CAF and NF (n = 3, respectively). The analysis revealed that CM by CAF was enriched in COL10A1 (p = 0.052) and SPARC (p = 0.004) as compared to NF, while only a trend in the enrichment of COL11A1 protein levels was found (p = 0.085) (Fig. 1b). We concluded that fibroblasts represent one possible source of the secreted proteins and that the differences in the release reflect their expression levels in cell cultures total lysates. To exclude protein deposition in a bound form in the extracellular matrix we also analyzed protein expression in the decellularized matrix without finding any of these proteins, thus indicating their preferential release.
ECM-related proteins as lung cancer biomarkers
To demonstrate the potential clinical utility of ECM-related proteins as lung cancer circulating biomarkers, we analyzed by ELISA assay matched plasma samples from 57 patients (TU) and 57 controls (CTR) for COL11A1 and COL10A1 and from 90 patients (TU) and 90 controls (CTR) for SPARC. The matching was performed according to clinical parameters (sex, age, smoke).
Main characteristics of the subjects included in the present study are shown in Table 2. Overall, the results showed that the levels of COL10A1 and SPARC were significantly higher in lung cancer patients compared to healthy controls (COL10A1: p = 0.001; SPARC: p < 0.001), while no difference in the levels of COL11A1 could be detected (p = 0.270) (Fig. 2). We further explored the levels of the proteins in patients and controls stratifying by sex (COL10A1 and COL11A1: M = 20, F = 37, SPARC: M = 62, F = 28 for both groups). Interestingly we found different levels between patients and controls for COL10A1 only in the female subgroup (p < 0.001) while no differences were found in males (p = 0.490) (Fig. 3a). The difference in SPARC levels between patients and controls remained significant across gender subgroups (p < 0.001) (Fig. 3b) while there were no differences for COL11A1 in gender subgroups (p = 0.232) in female group and p = 0.134 in male group) (Fig. 3c). In multivariable quantile regression analysis, the COL11A1 protein levels resulted as associated with disease status (p < 0.001), while COL10A1 was no longer associated to the disease status (p = 0.814) (Additional file 1: Table S1). No significant association between OS and COL11A1 levels (p = 0.683) or COL10A1 levels (p = 0.960) was found (Additional file 2: Figure S1). Accordingly, the two variables resulted not significant at multivariable Cox analysis (HR = 0.79 [0.25–2.52] p = 0.923 for COL11A1 and HR = 1.06 [0.26–4.35] for COL10A1 p = 0.300).
Reflecting the clinical scenario of resected NSCLC a fraction of patients had received pre-operative chemotherapy (14–16% of patients depending on cohorts, see Table 2). To verify whether chemotherapy (CT) treatment affected circulating ECM molecule levels we considered patients that were treated before surgery and we found that the CT treatment was not associated with COL11A1 and COL10A1 levels (Fig. 4).
Circulating SPARC is an independent diagnostic factor
Based on the inconsistent results on COL11A1 and COL10A1, we focused on SPARC protein. We performed univariable and multivariable analyses to assess the association between plasma SPARC levels and other established clinicopathological parameters such as age, smoking (in terms of packyears), tumor stage and COPD.
As already assessed SPARC was significantly strongly elevated in lung cancer patients compared to healthy controls (p < 0.001) and remained significant across gender subgroups (p < 0.001). In univariable analysis, no association between SPARC protein levels and any of the clinicopathological characteristic was detected. No significant difference in SPARC levels was detected among tumor stages (p = 0.916), suggesting that circulating SPARC levels could represent a potential valuable biomarker in early stage lung cancer (Fig. 5a). Also, the correlation between SPARC levels and packyears resulted as not statistically significant (p = 0.154, Fig. 5b). In multivariable analysis, the protein levels were associated with disease status, even after adjustment for age, packyears, COPD and sex (p ≤ 0.001) (Table 3). Therefore, in both univariable and multivariable analyses, SPARC levels resulted higher in the presence of lung cancer. This is also mirrored by the ROC curve where circulating SPARC levels showed a good capability to discriminate patients from controls (AUC = 0.744, Fig. 6a). The optimal cutoff maximizing the Youden Index of the ROC for SPARC was 0.587 (corrected - OD values) corresponding to a sensitivity of 64.4% [54.4–74.4%] and specificity of 78.9% [70.0–87.8%] (Additional file 3: Table S3), and to a PPV of 73.7% [64.8–83.3%] and a NPV of 72.5% [65.1–83.1%]. Therefore, if SPARC is over the optimal cutoff, the probability of a true positive finding is 73.7%, while for values below the optimal cutoff the probability of a true negative is 72.5%.
Concerning SPARC prognostic value, we did not find any association between OS and SPARC levels in tertiles both in the entire case series (p = 0.990) or considering early stage patients only (p = 0.461) (Fig. 6b, c). Also using the optimal cutoff for SPARC there was not any significant association (p = 0.971). Moreover, the molecule levels resulted not significant at multivariable Cox analysis, i.e. HR = 1.06 [0.57–1.99] (p = 0.593). Interestingly, SPARC levels were significantly lower for pre-treated patients (p = 0.037), indicating that SPARC levels could be an indicator of treatment response (Fig. 4).
In this study, we aimed to determine whether proteins related to the extracellular matrix are released in plasma and can potentially be useful as lung cancer biomarkers. We started our investigation from gene expression profiling of lung fibroblasts which are known to play an important role in cancer progression and ECM remodeling . This approach revealed that ECM related proteins are particularly enriched in primary cultures of CAF from lung cancer patients  and also released in their conditioned culture medium. We then focused on molecules belonging to the ECM3 signature, previously identified in breast cancer and found to be prognostic when overexpressed in the most aggressive tumors . We therefore analyzed two isoforms of collagens (COL11A1 and COL10A1 as the most expressed in lung fibroblasts gene profiling) and SPARC. Collagens represents the most abundant proteins of extracellular matrix and recent studies have highlighted an important role for COL11A1 and COL10A1 in many aspects of neoplastic progression . Previous observations reported that COL11A1 is more expressed in CAF than in normal breast or pancreatic fibroblasts [35,36,37]. Overexpression of COL11A1 was also found in non-small cell lung cancer tissue samples where it correlates with pathological stage, presence of lymph node metastasis, and poor prognosis [36, 38]. Our study showed higher COL11A1 gene expression in CAF compared to NF, but no difference in protein expression in cells or its release in culture medium. In addition, we did not find any difference in protein levels in plasma of patient’s vs heavy-smokers controls highlighting the inadequacy of this protein as potential indicator of pathological features. This observation indicates that COL11A1 is albeit differentially expressed in CAF is not sufficiently discriminatory at the protein level.
Instead, COL10A expression showed remarkable significant difference between controls and lung cancer patients thus constituting a potential diagnostic candidate. However, subgroup analyses showed that this finding was restricted to the female group. In recent years, the importance of gender related biomarkers has gained more attention especially for lung cancer . Several biological processes showed substantial differences between males and females in various hormonal states highlighting their impact on biomarker studies . In this work even if we do not have additional data or insights to potentially explain such differences we underline that gender effects should be considered before starting any biomarkers development study. Our data suggest however that COL10A1 could represent a potentially promising biomarker for lung cancer in females and that its relevance could be potentially explored in gender-related cancers such as breast or ovarian or in specific subgroups in other cancers. The most promising results were obtained from the study on SPARC protein. Although SPARC has recently emerged as a prognostic biomarker in different tumors [41,42,43], its role in lung cancer remains controversial . In non-small cell lung cancer, the localization (stromal or tumoral) of SPARC expression is associated with different disease prognosis. Absence of SPARC expression within the tumors is a negative prognostic factor  while high levels of SPARC protein expression, albeit rare, are associated with longer survival and could be protective against tumor aggressiveness . On the other hand, patients bearing SPARC-positive stroma have significantly poorer overall survival . Other studies found no prognostic impact of stromal SPARC expression . To our knowledge, this is the first study showing the release of high levels of SPARC in plasma of lung cancer patients to indicate that this protein could represent a useful diagnostic biomarker for lung cancer. Based on the previous literature and given the detectable levels of SPARC in conditioned medium of fibroblasts, we can hypothesize that the source of circulating SPARC in our patients is from stromal fibroblasts, rather than from cancer cells, and that protein levels in plasma could reflect changes in a microenvironment that becomes activated and permissive to tumor growth and progression.
SPARC protein levels in plasma were high in all stages of the disease indicating that this could also serve as marker of early lung cancer. Consistently with a previous study  we did not find any association between levels of circulating SPARC and the prognosis of our patients, reinforcing the hypothesis that its presence reflects early changes in the microenvironment more related with initial tumor growth than relapse and metastasis. Therefore, we can speculate that SPARC expression could influence stroma responsiveness during tumor formation, representing an indicator of the disease at its earlier phases. The effect of SPARC on malignant progression may instead depend on other events such as EMT, growth factor or immune modulation [47, 48].
Most importantly, the multivariate analysis also confirmed the significant association between disease status and levels of circulating SPARC. Interestingly, the observation that SPARC could be a promising diagnostic biomarker was confirmed by the ROC analysis, showing an AUC of 0.744, with an optimal cutoff corresponding to 64.4% of sensitivity and 78.9% of specificity. However further studies involving larger numbers of subjects are required to confirm these results.
Interestingly, it has been suggested that ECM modification provides protection against chemotherapy-induced apoptosis and may play a role in the failure of cancer therapy . Since the decrease (in terms of degradation) of ECM related proteins could be a marker of drug response we analyzed in this study patients that were treated with chemotherapy before surgery. Although our results are preliminary due to the low number of cases analyzed we report slightly decreased levels of SPARC protein in these patients indicating that beside a diagnostic marker SPARC may represent a marker of treatment response. Future longitudinal studies are however needed before firm conclusions can be drawn.
In conclusion, this is the first study to test circulating biomarkers related to ECM remodeling as possible diagnostic tools for lung cancer patients. SPARC emerged as the most promising biomarker, but it is possible that other genes identified in the comparative expression analysis of lung fibroblasts could be also useful for diagnostic purpose, either alone or in combination. It is now apparent, in fact, that panels of protein-based and nucleic acids-based cancer biomarkers, as opposed to single biomarkers, will probably be necessary for reliable cancer detection, especially to improve selection of high-risk individuals for CT screening and to distinguish malignant from benign nodules or identify patients with particularly aggressive cancers. Since the classical ELISA method that we used in our study suffer limitations in analysis time, sample size, equipment cost, and is not easily scalable to measure panels of proteins, new bioanalytical technologies should be developed to realize the full potential of protein biomarkers in the clinical setting. Our work represents an explorative study to verify whether proteins derived from ECM could be measured in plasma of lung cancer patients and their utility as circulating biomarkers and provides proof-of-concept on the feasibility and potential of this approach. Additionally, our study supports the concept that stroma-related plasma biomarkers may better fit as early diagnostic biomarkers than those strictly tumor-related. However, large prospective clinical studies are clearly warranted to confirm our preliminary results and further explore the existing potential of ECM-related circulating proteins.
In conclusion, we show here that circulating extracellular matrix related protein could be used as potential biomarkers for early diagnosis of lung cancer. In particular, SPARC emerged as the most promising biomarker representing an indicator of the disease at its earlier phases. Interestingly also COL10A expression showed remarkable significant difference between controls and lung cancer patients thus constituting a potential diagnostic candidate. However, subgroup analyses showed that this finding was restricted to the female group. Our study also highlighted the inadequacy of COL11A1 protein as potential indicator of pathological features. Taken together, our findings sustain the hypothesis that stroma-related plasma biomarkers could represent promising early diagnostic biomarkers related to interactions between incipient tumors and the surrounding microenvironment.
Cancer associated fibroblasts
Gene Set Enrichment Analysis
Non-small cell lung cancer
Secreted protein acidic and rich in cysteine
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, Parkin DM, Forman D, Bray F. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136:E359–86.
Bilfinger T, Keresztes R, Albano D, Nemesure B. Five-year survival among stage IIIA lung Cancer patients receiving two different treatment modalities. Med Sci Monit. 2016;22:2589–94.
Xu-Welliver M, Carbone DP. Blood-based biomarkers in lung cancer: prognosis and treatment decisions. Transl Lung Cancer Res. 2017;6:708–12.
Yi Z, Liu B, Guan X, Ma F. Plasma cell-free DNA and survival in non-small-cell lung cancer: a meta-analysis. Mol Clin Oncol. 2017;7:167–72.
A.A.Chaudhuri, J.J.Chabon, A.F.Lovejoy, A.M.Newman, H.Stehr, T.D.Azad, M.S.Khodadoust, M.S.Esfahani, C.L.Liu, L.Zhou, F.Scherer, D.M.Kurtz, C.Say, J.N.Carter, D.J.Merriott, J.C.Dudley, M.S.Binkley, L.Modlin, S.K.Padda, M.F.Gensheimer, R.B.West, J.B.Shrager, J.W.Neal, H.A.Wakelee, B.W.Loo, A.A.Alizadeh, M.Diehn, Early detection of molecular residual disease in localized lung cancer by circulating tumor DNA profiling, Cancer Discov.2017) CD-0716.
Xu T, Shen G, Cheng M, Xu W, Shen G, Hu S. Clinicopathological and prognostic significance of circulating tumor cells in patients with lung cancer: a meta-analysis. Oncotarget. 2017;8:62524–36.
Gallo M, De LA, Maiello MR, D’Alessio A, Esposito C, Chicchinelli N, Forgione L, Piccirillo MC, Rocco G, Morabito A, Botti G, Normanno N. Clinical utility of circulating tumor cells in patients with non-small-cell lung cancer. Transl Lung Cancer Res. 2017;6:486–98.
Gasiorowski L, Dyszkiewicz W, Zielinski P. In-vivo isolation of circulating tumor cells in non-small cell lung cancer patients by CellCollector. Neoplasma. 2017;64(6):938–44.
Jin X, Chen Y, Chen H, Fei S, Chen D, Cai X, Liu L, Lin B, Su H, Zhao L, Su M, Pan H, Shen L, Xie D, Xie C. Evaluation of tumor-derived Exosomal miRNA as potential diagnostic biomarkers for early-stage non-small cell lung Cancer using next-generation sequencing. Clin Cancer Res. 2017;23:5311–9.
Sozzi G, Boeri M, Rossi M, Verri C, Suatoni P, Bravi F, Roz L, Conte D, Grassi M, Sverzellati N, Marchiano A, Negri E, La VC, Pastorino U. Clinical utility of a plasma-based miRNA signature classifier within computed tomography lung Cancer screening: a correlative MILD trial study. J Clin Oncol. 2014;32:768–73.
Bissell MJ, Radisky D. Putting tumours in context. Nat Rev Cancer. 2001;1:46–54.
Sund M, Kalluri R. Tumor stroma derived biomarkers in cancer. Cancer Metastasis Rev. 2009;28:177–83.
Giussani M, Landoni E, Merlino G, Turdo F, Veneroni S, Paolini B, Cappelletti V, Miceli R, Orlandi R, Triulzi T, Tagliabue E. Extracellular matrix proteins as diagnostic markers of breast carcinoma. J Cell Physiol. 2018;233:6280–90.
Bremnes RM, Donnem T, Al-Saad S, Al-Shibli K, Andersen S, Sirera R, Camps C, Marinez I, Busund LT. The role of tumor stroma in cancer progression and prognosis: emphasis on carcinoma-associated fibroblasts and non-small cell lung cancer. J Thorac Oncol. 2011;6:209–17.
Kuchnio A, Moens S, Bruning U, Kuchnio K, Cruys B, Thienpont B, Broux M, Ungureanu AA, Leite de OR, Bruyere F, Cuervo H, Manderveld A, Carton A, Hernandez-Fernaud JR, Zanivan S, Bartic C, Foidart JM, Noel A, Vinckier S, Lambrechts D, Dewerchin M, Mazzone M, Carmeliet P. The Cancer cell oxygen sensor PHD2 promotes metastasis via activation of Cancer-associated fibroblasts. Cell Rep. 2015;12:992–1005.
Ohlund D, Elyada E, Tuveson D. Fibroblast heterogeneity in the cancer wound. J Exp Med. 2014;211:1503–23.
Hynes RO. The extracellular matrix: not just pretty fibrils. Science. 2009;326:1216–9.
White ES. Lung extracellular matrix and fibroblast function. Ann. Am. Thorac. Soc. 2015;12(Suppl 1):S30–3. https://doi.org/10.1513/AnnalsATS.201406-240MG. S30-S33
Burgess JK, Mauad T, Tjin G, Karlsson JC, Westergren-Thorsson G. The extracellular matrix - the under-recognized element in lung disease? J Pathol. 2016;240:397–409.
Pickup MW, Mouw JK, Weaver VM. The extracellular matrix modulates the hallmarks of cancer. EMBO Rep. 2014;15:1243–53.
Levental KR, Yu H, Kass L, Lakins JN, Egeblad M, Erler JT, Fong SF, Csiszar K, Giaccia A, Weninger W, Yamauchi M, Gasser DL, Weaver VM. Matrix crosslinking forces tumor progression by enhancing integrin signaling. Cell. 2009;139:891–906.
Willumsen N, Bager CL, Leeming DJ, Smith V, Christiansen C, Karsdal MA, Dornan D, Bay-Jensen AC. Serum biomarkers reflecting specific tumor tissue remodeling processes are valuable diagnostic tools for lung cancer. Cancer Med. 2014;3:1136–45.
Gandellini P, Andriani F, Merlino G, D’Aiuto F, Roz L, Callari M. Complexity in the tumour microenvironment: cancer associated fibroblast gene expression patterns identify both common and unique features of tumour-stroma crosstalk across cancer types. Semin. Cancer biol. Revised version submitted. 2015;35:96–106.
Triulzi T, Casalini P, Sandri M, Ratti M, Carcangiu ML, Colombo MP, Balsari A, Menard S, Orlandi R, Tagliabue E. Neoplastic and stromal cells contribute to an extracellular matrix gene expression profile defining a breast cancer subtype likely to progress. PLoS One. 2013;8:e56761.
Chen Y, Terajima M, Yang Y, Sun L, Ahn YH, Pankova D, Puperi DS, Watanabe T, Kim MP, Blackmon SH, Rodriguez J, Liu H, Behrens C, Wistuba II, Minelli R, Scott KL, Sanchez-Adams J, Guilak F, Pati D, Thilaganathan N, Burns AR, Creighton CJ, Martinez ED, Zal T, Grande-Allen KJ, Yamauchi M, Kurie JM. Lysyl hydroxylase 2 induces a collagen cross-link switch in tumor stroma. J Clin Invest. 2015;125:1147–62.
M. Giussani, G. Merlino, V. Cappelletti, E. Tagliabue, M.G. Daidone, Tumor-extracellular matrix interactions: Identification of tools associated with breast cancer progression, Semin Cancer Biol 35:3–10. doi: https://doi.org/10.1016/j.semcancer.2015.09.012. Epub;%2015 Sep 28. (2015) 3–10.
Hung JY, Yen MC, Jian SF, Wu CY, Chang WA, Liu KT, Hsu YL, Chong IW, Kuo PL. secreted protein acidic and rich in cysteine (SPARC) induces cell migration and epithelial mesenchymal transition through WNK1/snail in non-small cell lung cancer. Oncotarget. 2017;8:63691–702.
Tichet M, Prod’Homme V, Fenouille N, Ambrosetti D, Mallavialle A, Cerezo M, Ohanna M, Audebert S, Rocchi S, Giacchero D, Boukari F, Allegra M, Chambard JC, Lacour JP, Michiels JF, Borg JP, Deckert M, Tartare-Deckert S. Tumour-derived SPARC drives vascular permeability and extravasation through endothelial VCAM1 signalling to promote metastasis. Nat Commun. 2015;6:6993. https://doi.org/10.1038/ncomms7993. 6993
Koukourakis MI, Giatromanolaki A, Brekken RA, Sivridis E, Gatter KC, Harris AL, Sage EH. Enhanced expression of SPARC/osteonectin in the tumor-associated stroma of non-small cell lung cancer is correlated with markers of hypoxia/acidity and with poor prognosis of patients. Cancer Res. 2003;63:5376–80.
Boeri M, Verri C, Conte D, Roz L, Modena P, Facchinetti F, Calabro E, Croce CM, Pastorino U, Sozzi G. MicroRNA signatures in tissues and plasma predict development and prognosis of computed tomography detected lung cancer. Proc Natl Acad Sci U S A. 2011;108:3713–8.
Andriani F, Bertolini G, Facchinetti F, Baldoli E, Moro M, Casalini P, Caserini R, Milione M, Leone G, Pelosi G, Pastorino U, Sozzi G, Roz L. Conversion to stem-cell state in response to microenvironmental cues is regulated by balance between epithelial and mesenchymal features in lung cancer cells. Mol Oncol. 2016;10:253–71.
M.E. Ritchie, B. Phipson, D. Wu, Y. Hu, C.W. Law, W. Shi, G.K. Smyth, limma powers differential expression analyses for RNA-sequencing and microarray studies, nucleic acids res. %20;43 (2015) e47.
Durrleman S, Simon R. Flexible regression models with cubic splines. Stat Med. 1989;8:551–61.
Jia D, Liu Z, Deng N, Tan TZ, Huang RY, Taylor-Harding B, Cheon DJ, Lawrenson K, Wiedemeyer WR, Walts AE, Karlan BY, Orsulic S. A COL11A1-correlated pan-cancer gene signature of activated fibroblasts for the prioritization of therapeutic targets. Cancer Lett. 2016;382:203–14.
Garcia-Pravia C, Galvan JA, Gutierrez-Corral N, Solar-Garcia L, Garcia-Perez E, Garcia-Ocana M, Del Amo-Iribarren J, Menendez-Rodriguez P, Garcia-Garcia J, de Jr LT, Simon-Buela L, Barneo L. Overexpression of COL11A1 by cancer-associated fibroblasts: clinical relevance of a stromal marker in pancreatic cancer. PLoS One. 2013;8:e78327.
Shen L, Yang M, Lin Q, Zhang Z, Zhu B, Miao C. COL11A1 is overexpressed in recurrent non-small cell lung cancer and promotes cell proliferation, migration, invasion and drug resistance. Oncol Rep. 2016;36:877–85.
M.Karaglani, I.Toumpoulis, N.Goutas, N.Poumpouridou, D.Vlachodimitropoulos, S.Vasilaros, I.Rizos, C.Kroupis, Development of novel real-time PCR methodology for quantification of COL11A1 mRNA variants and evaluation in breast cancer tissue specimens, BMC Cancer 15:694. doi: https://doi.org/10.1186/s12885-015-1725-8. (2015) 694–1725.
Chong IW, Chang MY, Chang HC, Yu YP, Sheu CC, Tsai JR, Hung JY, Chou SH, Tsai MS, Hwang JJ, Lin SR. great potential of a panel of multiple hMTH1, SPD, ITGA11 and COL11A1 markers for diagnosis of patients with non-small cell lung cancer. Oncol Rep. 2006;16:981–8.
de Torres JP, Casanova C, Pinto-Plata V, Varo N, Restituto P, Cordoba-Lanus E, Baz-Davila R, guirre-Jaime A, Celli BR. Gender differences in plasma biomarker levels in a cohort of COPD patients: a pilot study, PLoS. ONE. 2011;6:e16021.
Dougherty SM, Mazhawidza W, Bohn AR, Robinson KA, Mattingly KA, Blankenship KA, Huff MO, McGregor WG, Klinge CM. Gender difference in the activity but not expression of estrogen receptors alpha and beta in human lung adenocarcinoma cells, Endocr. Relat Cancer. 2006;13:113–34.
Li Z, Li AD, Xu L, Bai DW, Hou KZ, Zheng HC, Qu XJ, Liu YP. SPARC expression in gastric cancer predicts poor prognosis: results from a clinical cohort, pooled analysis and GSEA assay. Oncotarget. 2016;7:70211–22.
Chen ZY, Zhang JL, Yao HX, Wang PY, Zhu J, Wang W, Wang X, Wan YL, Chen SW, Chen GW, Liu YC. Aberrant methylation of the SPARC gene promoter and its clinical implication in Gastric Cancer. Sci Rep. 2014;4:7035. https://doi.org/10.1038/srep07035. 7035
Wang HY, Li YY, Shao Q, Hou JH, Wang F, Cai MB, Zeng YX, Shao JY. Secreted protein acidic and rich in cysteine (SPARC) is associated with nasopharyngeal carcinoma metastasis and poor prognosis. J Transl. Med. 2012;10:27. https://doi.org/10.1186/1479-5876-10-27. 27–10
Wong SL, Sukkar MB. The SPARC protein: an overview of its role in lung cancer and pulmonary fibrosis and its potential role in chronic airways disease. Br J Pharmacol. 2017;174:3–14.
Y.Huang, J.Zhang, Y.Y.Zhao, W.Jiang, C.Xue, F. Xu, H.Y. Zhao, Y. Zhang, L.P. Zhao, Z.H. Hu, z.W. Yao, Q.Y. Liu, L. Zhang, SPARC expression and prognostic value in non-small cell lung cancer, Chin J Cancer 31 (2012) 541–548.
Edlund K, Lindskog C, Saito A, Berglund A, Ponten F, Goransson-Kultima H, Isaksson A, Jirstrom K, Planck M, Johansson L, Lambe M, Holmberg L, Nyberg F, Ekman S, Bergqvist M, Landelius P, Lamberg K, Botling J, Ostman A, Micke P. CD99 is a novel prognostic stromal marker in non-small cell lung cancer. Int J Cancer. 2012;131(10):2264–73.
Tripodo C, Burocchi A, Piccaluga PP, Chiodoni C, Portararo P, Cappetti B, Botti L, Gulino A, Isidori A, Liso A, Visani G, Martelli MP, Falini B, Pandolfi PP, Colombo MP, Sangaletti S. Persistent immune stimulation exacerbates genetically driven myeloproliferative disorders via stromal remodeling. Cancer Res. 2017;77:3685–99.
Sangaletti S, Tripodo C, Santangelo A, Castioni N, Portararo P, Gulino A, Botti L, Parenza M, Cappetti B, Orlandi R, Tagliabue E, Chiodoni C, Colombo MP. Mesenchymal transition of high-grade breast carcinomas depends on extracellular matrix control of myeloid suppressor cell activity. Cell Rep. 2016;17:233–48.
Junttila MR, de Sauvage FJ. Influence of tumour micro-environment heterogeneity on therapeutic response. Nature. 2013;501:346–54.
The authors thank Claudio Citterio for data management and Paola Suatoni for handling of samples. The work was supported by the Italian Association for Cancer Research (Special Program “Innovative Tools for Cancer Risk Assessment and early Diagnosis”, 5 × 1000).
The work was supported by the Italian Association for Cancer Research [Investigator Grants No. 14318 and 18812 to GS, 16847 to LR and 12162 to UP, GS and LR (Special Program “Innovative Tools for Cancer Risk Assessment and Early Diagnosis”, 5 × 1000)]. The funding agencies were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Experiments were approved by the Independent Ethics Committee of the INT (18/12/2012–0016337). Ethics committees approved written informed consent, which was obtained from all the participants prior participating to the study.
Consent for publication
All the authors of this paper have approved the final version of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Quantile regression multivariable analysis of association between COL11A1 and COL10A1 levels and the characteristics of 57 lung cancer patients and 57 healthy controls. (DOCX 15 kb)
Figure S1. COL11A and COL10A levels are not associated with overall survival. Kaplan Meier overall survival (OS) curves for COL11A1 and COL10A1 levels divided into three categories (≤1st tertile, 1st-2nd tertile, >2nd tertile) on the entire case series p = log-rank test p-value. (PPTX 51 kb)
Table S2. Sensitivity and specificity data for SPARC cutoffs. (DOCX 18 kb)
About this article
Cite this article
Andriani, F., Landoni, E., Mensah, M. et al. Diagnostic role of circulating extracellular matrix-related proteins in non-small cell lung cancer. BMC Cancer 18, 899 (2018). https://doi.org/10.1186/s12885-018-4772-0
- Lung cancer
- Extracellular matrix (ECM)
- Circulating biomarkers
- Secreted protein acidic and rich in cysteine
- Cancer associated fibroblasts