Skip to main content

Identification of epigenetic silencing of the SFRP2 gene in colorectal cancer as a clinical biomarker and molecular significance

Abstract

Background

Several studies have suggested secreted frizzled-related protein 2 (SFRP2) gene as a potential clinical biomarker in colorectal cancer (CRC). However, its diagnostic role remains unclear. In this study, we aimed to investigate the significance of SFRP2 methylation levels in a large cohort of biological specimens (including blood, adipose and colonic tissues) from patients with CRC, thereby potentially identifying new biomarker utility.

Methods

We examined the expression (by qPCR) and methylation status (by 450 K DNA array and DNA pyrosequencing) of the SFRP2 gene in healthy participants (N = 110, aged as 53.7 (14.2), 48/62 males/females) and patients with CRC (N = 85, aged 67.7 (10.5), 61/24 males/females), across different biological tissues, and assessing its potential as a biomarker for CRC. Additionally, we investigated the effect of recombinant human SFRP2 (rhSFRP2) as a therapeutic target, on cell proliferation, migration, and the expression of key genes related to carcinogenesis and the Wnt pathway.

Results

Our findings revealed that SFRP2 promoter methylation in whole blood could predict cancer stage (I + II vs. III + IV) (AUC = 0.653), lymph node invasion (AUC = 0.692), and CRC recurrence (AUC = 0.699) in patients with CRC (all with p < 0.05). Furthermore, we observed a global hypomethylation of SFRP2 in tumors compared to the adjacent area (p < 0.001). This observation was validated in the TCGA-COAD and TCGA-READ cohorts, demonstrating overall hypermethylation (both with p < 0.001) and low expression (p < 0.001), as shown in publicly available scRNA-Seq data. Notably, neoadjuvant-treated CRC patients exhibited lower SFRP2 methylation levels compared to untreated patients (p < 0.05) and low promoter SFRP2 methylation in untreated patients was associated with poor overall survival (p < 0.05), when compared to high methylation. Finally, treatment with 5 µg of rhSFRP2 treatment in CRC cells (HCT116 cells) inhibited cell proliferation (p < 0.001) and migration (p < 0.05), and downregulated the expression of AXIN2 (p < 0.01), a gene involved in Wnt signaling pathway.

Conclusions

These findings establish promoter methylation of the SFRP2 gene as a prognostic candidate in CRC when assessed in blood, and as a therapeutic prognostic candidate in tumors, potentially valuable in clinical practice. SFRP2 also emerges as a therapeutic option, providing new clinical and therapeutical avenues.

Introduction

Colorectal cancer (CRC) is the third most common cancer and the fourth leading cause of cancer associated mortality. By 2040, it is estimated that the incidence of CRC will increase up to 63% increase, reflecting a 73% rise in mortality rates [1]. The overall 5-year survival rate for CRC stands at 63%. However, diagnosis at early and localized stage rises the survival rate up to 91%. Nevertheless, advanced cases decrease drastically the overall survival, up to 12–14% [2, 3], clearly positioning CRC as a major public health problem.

Generally, CRC originates from genetic and epigenetic alterations in colonic epithelial cells [4, 5], which play a pivotal role in both carcinogenesis and cancer progression. Epigenetic alterations in cancer are considered the hallmark of human tumors. Specifically, DNA methylation is the most significant epigenetic alterations in CRC, mainly affecting promoter regions of tumor suppressor genes, leading to their inactivation and contributing to uncontrolled cell proliferation [6]. In particular, epigenetic alterations in the secreted frizzled-related protein 2 (SFRP2) gene have been proposed as a potential biomarker in CRC.

SFRP2 is an extracellular Wnt antagonist that binds to Wnt ligands, and inhibit the Wnt signaling. This inhibition prevents β-catenin from translocating from the cytoplasm to the nucleus, affecting cell proliferation, survival, differentiation, and migration [7]. The latest meta-analysis conducted by Yu et al. (2019), which included 83 studies with a total of 21,612 samples, identified a significant association between SFRP2 promoter hypermethylation (in different biological tissues) and an increased risk of CRC, with an overall pooled odds ratio (OR) of 8.41 (95% CI 5.91–11.97) [8]. Another meta-analysis, conducted by Shariatpanahi et al., found that promoter methylation of SFRP2 in stool samples was associated with an increased risk of CRC (OR = 35.3 [18.7–66.8], AUC = 0.94) [9], indicating its potential use as a diagnostic biomarker. However, the sensitivity of SFRP2 methylation as a biomarker for CRC varies widely. Moreover, the level of SFRP2 methylation differs among different tissues. While SFRP2 methylation in stool samples can be useful for cancer diagnosis, SFRP2 methylation in blood samples may serve as a prognostic factor [10]. This variability can be attributed to various factors that may disrupt SFRP2 methylation status. In particular, it has been reported that SFRP2 methylation is significantly higher on the right side of CRC (57.33%) compared to the left side (36.60%) [11]. Therefore, assessing SFRP2 methylation in other tissues may serve as a stable prognostic factor.

Nevertheless, there are limited studies investigating the potential of SFRP2 methylation as a predictor for other cancer outcomes, such as response to therapy or recurrence in other biological and less invasive tissues. This highlights the need for further research to identify its predictive value. Additionally, there are few studies concerning the use of recombinant human SFRP2 (rhSFRP2) in cancer research, emphasizing the urgency of exploring its potential as a therapeutic target. Therefore, the aim of this study is to test whether there is a significant association between SFRP2 DNA methylation and other cancer outcomes, such as survival rates, response to therapy, location, or recurrence in blood and tumor samples. Furthermore, we aimed to evaluate the effect of recombinant human SFRP2 in colorectal cancer cell model (HCT116 cells) to understand its contribution to colorectal carcinogenesis, and to recommend new strategies for its use as a cancer biomarker.

Materials and methods

Participants and study design

This study recruited participants from the “Virgen de la Victoria” University Hospital between 2012 and 2014. A total of 110 healthy participants (without CRC) were included, who underwent hiatus hernia surgery or cholecystectomy, which was the method used to obtain adipose tissue. Additionally, 85 patients diagnosed with CRC were included, diagnosed by a pathology specialist using a biopsy and colonoscopy. Biopsy samples were classified according to histological features by pathologists, according to the “World Health Organization Classification of Tumors of the Digestive System” (2016) [12]. Additional medical records and pathological examinations were reviewed to confirm the diagnosis. All patients with CRC underwent surgery with curative intent, specifically, a hemicolectomy and lower anterior resection with an ileostomy, followed by total meso-colorectal excision. For rectal cases, total mesorectal excision preceded. Patients with CRC were treated with radiation therapy and/or fluoropyrimidine-based chemotherapy. All patients had at least a 5-year follow-up, including clinical visits, checkups, and biochemical variables every three months for the first two years, then every six months starting in the third year. A colonoscopy was performed annually. Exclusion criteria were patients with inflammatory and cardiovascular diseases, hereditary non-polyposis colorectal cancer or familial adenomatous polyposis, type 2 diabetes, insulin resistance, or renal and infectious diseases. Participants who had undergone treatment that altered lipid or glucose metabolism or who consumed > 20 g of ethanol per day were excluded. Written informed consent was obtained from all patients and subjects and reviewed and approved by the Ethics Committees of the “Virgen de la Victoria” University Hospital (Málaga, Spain) (Registration Number: 0311/PI7). The study was conducted in accordance with the guidelines laid down in the Declaration of Helsinki.

Samples included in the study

Blood samples were collected from all participants of the same patients, and serum samples were obtained by centrifugation of the blood samples at 4000 rpm for 15 min at 4 °C. Epiploic adipose tissue samples (surrounding and attached to the colon) were obtained during surgery and immediately frozen in liquid nitrogen for transportation to the laboratory. The adipose tissue samples were then washed with PBS, cut into 100 mg pieces, and stored at −80 °C until processing. Biopsy samples were used to obtain tumor samples containing both the tumor area and the adjacent normal tumor-free area (NAT). These samples were fixed using formalin and embedded in paraffin (formalin-fixed paraffin-embedded or FFPE). Pathologists determined the limits between the areas and calculated the percentage of tumor cells within the tumor area, which was greater than 80% of the tumor cell percentage.

Biochemical variables determination

Fasting glucose, total cholesterol, triglycerides, and high-density lipoprotein cholesterol (HDL-c) were measured using the Dimension Autoanalyzer (Dade Behring Inc., USA). Low-density lipoprotein cholesterol (LDL-c) was calculated using the Friedewald equation [13]. Insulin levels were determined by radioimmunoassay methods using BioSource International Inc. (Camarillo, CA, USA). The homeostasis model assessment of insulin resistance (HOMA-IR) was calculated using the following equation: HOMA-IR = fasting insulin (μIU/mL) x fasting glucose (mmol/L)/22.5 [14].

DNA and RNA extraction

Genomic DNA was extracted from 200 µl of blood samples and HCT116 cells (1 × 104–1 × 105 cells) using the Qiamp DNA mini-Kit (Qiagen GmbH, Hilden, Germany), following the manufacturer's instructions. For adipose tissue, 30 mg of tissue were used, and genomic DNA was extracted using the Qiamp DNA Tissue Kit (Qiagen GmbH, Hilden, Germany) according to the manufacturer's instructions. Genomic DNA from the tumor and NAT area was isolated from 5 to 10 sections of 14 μm of FFPE samples, using the Qiamp DNA FFPE Tissue Kit according to the manufacturer's instructions (Qiagen GmbH, Hilden, Germany). DNA integrity was determined using the NanoDrop ND-1000 (Thermo Fisher Scientific, Indianapolis, USA) at 260/280 and 260/230 ratios and confirmed by electrophoresis in a 1.5% agarose gel in TAE buffer. The Pico Green dsDNA Quantitation Reagent (Invitrogen, Carlsbad, CA) was used to further analyze DNA integrity and quantify DNA concentrations.

For RNA extraction, total RNA was isolated from HCT116 cells (at a density of 1 × 104–1 × 105 cells) using the Total RNA Purification Kit (Norgen Biotek Corp., ON, Canada) according to the manufacturer’s instructions. Total RNA from blood and adipose tissue samples was isolated using the Qiamp RNA Easy Kit (Qiagen GmbH, Hilden, Germany) and the RNeasy Lipid Tissue Mini Kit (Qiagen GmbH, Hilden, Germany), respectively, according to the manufacturer's instructions. Total RNA from FFPE tumor samples was isolated from 5 to 10 sections of 14 μm using the RNeasy FFPE Kit, following the manufacturer's instructions (Qiagen GmbH, Hilden, Germany). RNA quantity and integrity was determined using the NanoDrop ND-1000 (Thermo Fisher Scientific, Indianapolis, USA) at 260/280 and 260/230 ratios. To generate first-strand cDNA synthesis, 1 μg of total extracted RNA and random primers were used with the PrimeScript RT-PCR Kit as indicated by the manufacturer.

Bisulfite reaction and genome-wide DNA methylation analysis

For the genome-wide methylation analysis study, high-quality genomic DNA samples (500 ng) from the tumor area (N = 27) and the NAT area (N = 15) were treated with bisulfite using the EZ-96 DNA Methylation Kit (Zymo Research, Irvine, CA), following the manufacturer's instructions. Subsequently, DNA methylation was analyzed by microarray assays using Infinium Human Methylation 450 K bead chip technology (Illumina, San Diego, CA). Whole-genome amplification and hybridization were then performed using BeadChip, followed by single-base extension and analysis using the HiScan SQ module (Illumina) to assess cytosine methylation states. DNA quality checks, bisulfite modification, hybridization, data normalization, statistical filtering, and value calculation were performed as previously described [15]. For the pyrosequencing reaction, purified DNA (2 μg) from different tissues was used for the bisulfite reaction, using the EpiTect Fast Bisulfite Kit (Qiagen GmbH, Hilden, Germany). The primer sequences for the promoter of the SFRP2 gene were as follows: forward primer: AGAAGTTTTGGGTTTAGTTTATGAT, and reverse primer: CTCACATCTACCCAATATAAAAACTCACCA. A PCR reaction was performed using 0.2 nmol/L of primers. DNA pyrosequencing was carried out using the PyroMark Q96 ID pyrosequencing System (Qiagen GmbH, Hilden, Germany) and the sequencing primer: TGTTGAAYGGTGGTTGG. The methylation average was presented as the percentage of methylated cytosine over the sum of methylated and unmethylated cytosines. Interassay precision (%CV) was 2.5%, and intraassay precision (%CV) was 1.0%. Non-CpG cytosine residues were used as built-in controls to verify bisulfite conversion. We also used unmethylated and methylated DNA as controls in our assay (New England Biolabs, UK).

Gene expression analysis

For the quantification of SFRP2 (Hs.PT.58.20705989), AXIN2 (Hs.PT.58.39305692), CTNNB1 (beta-catenin) (Hs.PT.58.40551289), TIAM1 (Tumor Growth Factor beta 1) (Hs.PT.58.24232612), VEGF (Vascular Endothelial Growth Factor) (Hs.PT.58.21234833) (all related with the Wnt pathway) and PPIA (Peptidylprolyl isomerase A) (Hs.PT.58v.38887593.g) genes, we used commercially available TaqMan primer/probe mix (Integrated DNA Technologies Inc., Madrid, Spain). Gene expression was quantified using Premix Ex Taq (Probe qPCR) (Takara Bio USA, Inc., Mountain View, CA), following the manufacturer's instructions on the QuantStudio 6 Pro (Applied Biosystems, Darmstadt, Germany). Gene expression was normalized by the 2−ΔCt method [16], expressed as the target gene/PPIA ratio.

Cell culture

HCT116 (obtained from ATCC with the reference: CCL-247) human colorectal carcinoma cells were cultured in Dulbecco’s Modified Eagle’s Medium—low glucose (Biowest, France), plus 10% fetal bovine serum (FBS) (Gibco, CA, USA), 1% of l-glutamine and 1% of streptomycin/penicillin at 37 °C and 5% CO2. Cells were periodically checked for mycoplasma using DAPI analysis [17] and PCR using universal primers [18]. Cells are maintained at passages lower than 20. All cell experiments were conducted in triplicates and three independent replicates.

DNA demethylation in vitro

For the 5-Aza-2-deoxycytidine (AZA) (A3656, Sigma Aldrich, Madrid, Spain) treatment as an inhibitor of DNA methyltransferases (DNMTs), cells were treated with 1 µM, 5 µM, and 10 µM of AZA for 72 h (increased concentration to test if the effect is dose-dependent). The media was renewed every 48 h, and RNA extraction was performed after 72 h of treatment.

Cell proliferation assay

Cells were seeded in 96-well plates at a concentration of 5,000 cells/well. After 24 h, the cells were treated with 0.1, 0.5, 1, 2.5, and 5 µg/mL of rhSFRP2 (recombinant human SFRP2) for 48 h (to find the concentration that cause decrease in cell proliferation). After the treatment, the media were removed, and the MTT reagent was added using an MTT Assay Kit (Cell Proliferation, Cat. No. ab211091, Abcam, UK) and incubated for 3 h at 37 ºC. Then, the crystals were dissolved in the solvent reagent, and the absorbance at 590 nm was measured.

Migration assay

A 96-well tissue culture plate was seeded with cells at a density of 50,000 cells/well. After 24 h, a circular wound was created by scratching the monolayer as described in [19]. A p10 pipette tip attached to a vacuum pump was gently pressed perpendicularly onto the cell monolayer, without lateral movement, to detach cells from the substratum and create a circular wound. Cells were then treated with 5 µg/mL of rhSFRP2, and images were captured at 0, 24, and 48 h. Migration was quantified by calculating the gap distance and expressing it as a percentage of control (i.e., 100% wound healing). The gap distance was measured using ImageJ software (http://rsb.info.nih.gov/ij/download.html).

Bioinformatic analysis

For DNA methylation analysis, we used the minfi package and the ChAMP pipeline from Bioconductor with default settings [20, 21]. Probes with p-values above 0.01 were excluded, and normalization and preprocessing were conducted using the minfi package. Probes were stratified by region and bad samples (low intensity) were eliminated using the median of methylated and unmethylated signals for each sample. Probes with single nucleotide polymorphisms at CpG or single-base extension sites and those on sex chromosomes were removed, and M values were obtained for further investigation.

For TCGA data analysis, we used the “TCGAMethylation450k” and “TCGAbiolinks” packages to extract data from The Cancer Genome Atlas (TCGA) database (Cancer Genome Atlas Research) using mainly idat files and information of participants. We obtained methylation data from June 2023, for TCGA-COAD (colon adenocarcinoma) and TCGA-READ (rectal adenocarcinoma) from The Cancer Genome Atlas (TCGA) database (Cancer Genome Atlas Research) (with the primary objective to validate and assess the consistency of our results, as well as increase confidence in the robustness and reproducibility of our findings) [22]. We included IDAT files from 308 cancer cases and 41 NAT from TCGA-COAD, and 94 cancer cases and 10 NAT from TCGA-READ, which were downloaded and processed using the “TCGAbiolinks” R packages [23]. The data were filtered, normalized, and preprocessed following similar protocols as described above (probes with p-values exceeding 0.01 were omitted). Survival analyses were extracted from the OncoDB web tool [24]. We used normalized RNA-sequencing (RNA-seq) expression data from TCGA-COAD and TCGA-READ, obtained using the TCGAbiolinks R package, to identify differentially expressed genes. For the single-cell analysis, we extracted the results from single-cell analysis using 52,609 cells, in June 2023, as published by Lee HO et al. [25], using the EMBL-EBI platform: https://www.ebi.ac.uk/

Statistical analysis

The results are presented as mean ± standard deviation (SD) for continuous variables and as numbers (percentages) for categorical variables. Depending on the normality of the variables, either a Student’s t-test or a Wilcoxon test was used. Pearson correlation coefficients were calculated to evaluate the associations between methylation and anthropometric and biochemical parameters. Overall survival analyses were conducted using Kaplan–Meier curves, and multivariate proportional regression was used to calculate odds ratios (ORs) [95% confidence intervals (CIs)] and area under curve (AUC) for SFRP2 methylation adjusted by age, sex, BMI, HDL, triglycerides and cancer location. Logistic and linear regression analyses were conducted using SFRP2 promoter methylation in tumor, adjusted for age, sex and BMI. All analyses were performed in R version 3.5.1, and statistical significance was set at p < 0.05 [26]. Graphical representations were also created using R.

Results

General clinical and characteristic data of the participants

The baseline biochemical and anthropometric data of healthy participants and patients with CRC are summarized in Table 1. As observed, patients with CRC were older than healthy participants (p < 0.001). Furthermore, glucose and triglyceride levels were higher in patients with CRC compared to healthy participants (p = 0.004 and p = 0.032, respectively). In contrast, HOMA-IR, total cholesterol, LDL-c and HDL-c levels were lower in patients with CRC than in healthy participants (all with p < 0.001). Importantly, there was no significant difference in BMI when comparing both groups.

Table 1 Baseline table of the participants included in the study

DNA methylation of the SFRP2 promoter in whole blood as a biomarker in colorectal cancer

In this study, we measured the expression and methylation levels of the SFRP2 gene in whole blood and adipose tissue. As depicted in Fig. 1A, there were no significant differences in SFRP2 promoter methylation levels analyzed by pyrosequencing between healthy participants and patients with CRC in whole blood. Additionally, we did not observe any significant differences in the prognostic value of promoter SFRP2 methylation levels in whole blood between high and low methylation levels (Fig. 1B). Similar results were observed in SFRP2 gene expression levels between healthy participants and patients with CRC in whole blood (Supplementary Figs. 1A and B).

Fig. 1
figure 1

Promoter SFRP2 methylation in whole blood as a biomarker for colorectal cancer. A We measured the promoter methylation of SFRP2 gene from healthy participants (N = 31) and patients with CRC (N = 62) (p = 0.220). B Kaplan–Meier curve comparing the median of the SFRP2 promoter methylation as low and high methylation. The significance of differences is evaluated with the Log-rank test. We measured the promoter methylation of SFRP2 gene from patients with CRC in C early stage (I + II) (N = 42) and late stage (III + IV) (N = 39), (p = 0.044) D absence (N = 44) and presence (N = 35) of lymph node invasion (p = 0.014) and E absence (N = 58) and presence (N = 24) of recurrence (p = 0.013). Asterisks indicate significant differences between the groups according to the Mann Whitney test (*p < 0.05, **p < 0.01, ***p < 0.001). F, G and H indicate ROC curve of stage (F), lymph node invasion (G) and Recurrence (H) in three models. Model 1 indicate SFRP2 methylation. Model 2 indicates clinical variables, such as age, sex, BMI, HDL, triglycerides, and tumor location. Model 3 includes Model 1 + Model 2. AUC Area under curve, CRC colorectal cancer, ROC receiver operating characteristic curve, Sens sensitivity, Spec specificity

When we examined other cancer outcomes, we found that the promoter SFRP2 methylation levels were significantly lower in patients with CRC in late stage (III + IV) than in those in the early stage (I + II) (p = 0.044) (Fig. 1C), and also when compared with healthy participants (p = 0.035) (Supplementary Fig. 1C), with an AUC of 0.879 (95% CI 0.700–0.942) (Fig. 1F, Supplementary Table 1) when adjusted by clinical variables. Patients with CRC who had lymph node invasion also had lower SFRP2 methylation levels than those without lymph node invasion (p = 0.014) (Fig. 1D), with an AUC of 0.844 (95% CI 0.736–0.951) (Fig. 1G, Supplementary Table 1). Finally, patients with CRC who had recurrence had lower SFRP2 methylation levels than those without recurrence (p = 0.013) (Fig. 1E), with an AUC of 0.821 (95% CI 0.700–0.942) (Fig. 1H, Supplementary Table 1). We also measured the expression and methylation levels of the SFRP2 gene in adipose tissue. As shown in Supplementary Fig. 1D, there were no significant differences in SFRP2 promoter methylation levels between healthy participants and patients with CRC. However, patients with CRC had lower SFRP2 expression levels than healthy participants (p = 0.002) (Supplementary Fig. 1F). Furthermore, we did not observe any significant differences between high and low SFRP2 methylation or expression levels in adipose tissue (Supplementary Fig. 1E and Supplementary Fig. 1G).

DNA Methylation of the SFRP2 gene in tumoral tissue from the same cohort as a biomarker in colorectal cancer

We next evaluated the expression and methylation of the SFPP2 gene in colonic tissue taking advantage of Infinium DNA methylation 450 K data (adjusted and age, sex and BMI). As shown in Fig. 2A, global β methylation of the SFRP2 gene was higher in the tumor area, when compared to the NAT area (p < 0.001). However, we did not find a significant difference between both areas regarding SFRP2 expression (p = 0.171) (Fig. 2B). After that, we next validated these findings in the TCGA-COAD and TCGA-READ. As a result, we discovered that global β methylation of the SFRP2 gene was higher in the tumor area than in the NAT area (p < 0.001) in both TCGA-COAD and -READ (Fig. 2C, F). Furthermore, when we focused on the prognostic value of promoter SFRP2 methylation in tumors and in the TCGA-COAD and TCGA-READ, we found no significant difference between high and low SFRP2 methylation in global β methylation, promoter methylation, and in the TCGA-COAD and TCGA-READ (Supplementary Fig. 2A, Supplementary Fig. 2B, Supplementary Fig. 2C, Supplementary Fig. 2D). In addition, we also found that SFRP2 expression was lower in the tumor area in comparison with the NAT area (p < 0.001) (Fig. 2E, F). Next, we analyzed the single-cell RNA sequencing of CRC tumors. We found that SFRP2 expression was mostly expressed in the NAT, although there was a small fraction of the tumor core and tumor border areas that showed a relative increase in SFRP2 expression. At cell level, we observed that SFRP2 was predominately expressed in stromal cells (located in the NAT area) and myofibroblasts (located in the tumor core and tumor border areas).

Fig. 2
figure 2

Diagnostic and prognostic value of SFRP2 methylation in tumor tissue from colorectal cancer. A Global β methylation was done with significant CpG sites between tumor and NAT. The CpG sites included in the global SFRP2 methylation in colon tissue in our participants were: cg14289246, cg05961809, cg23292160, cg25775322, cg23207990, cg22178613, cg10318528, cg11354906, cg00082664, cg03202804, cg05874561, cg04965141, cg25645268, cg14330641, cg23121156, cg06549216, cg10942078, cg20881942, cg23714408, cg05164933, cg14063488. B SFRP2 expression in colon samples from NAT (N = 11) and tumor area(N = 16). C Global β methylation was done with significant CpG sites between tumor and NAT in the TCGA-COAD. D SFRP2 expression in colon samples from NAT (N = 51) and tumor area (N = 308) in the TCGA-COAD cohort. E Global β methylation was done with significant CpG sites between tumor and NAT in the TCGA-READ. F SFRP2 expression in colon samples from NAT (N = 10) and tumor area(N = 94) in the TCGA-READ cohort. G Single cell RNA-seq analysis of SFRP2 gene in NAT, tumor core and tumor border. SFRP2 was mostly expressed in NAT and a small fraction of tumor core, in stromal cells and myofibroblasts. H We measured the promoter methylation of SFRP2 gene from patients with CRC no treated (N = 57) and treated (N = 24) with neoadjuvant therapy. I Kaplan–Meier curve comparing the median of the SFRP2 promoter methylation as low and high methylation in non-treated patients with CRC. The significance of differences is evaluated with the Log-rank test. Asterisks indicate significant differences between the groups according to the Mann Whitney test (*p < 0.05, **p < 0.01, ***p < 0.001). COAD colon adenocarcinoma, NAT normal adjacent area, READ rectum adenocarcinoma, TCGA The cancer genome atlas

Next, we investigated those variables that can affect the promoter methylation of SFRP2 in patients with CRC. For that, we conducted a linear regression analysis adjusted for age, sex, and BMI. We found that the location of the tumor (colon vs. rectum) was associated with promoter SFRP2 methylation (p < 0.05) (Table 2). For this purpose, we found that the promoter SFRP2 methylation in colon was higher than that in rectum (p = 0.025) (Supplementary Fig. 1E). In addition, we found that neoadjuvant therapy was strongly associated with SFRP2 methylation (p < 0.001) (Table 2). Accordingly, we found that treated CRC patients had lower promoter SFRP2 methylation when compared with non-treated CRC patients (p < 0.001). We discovered in the Kaplan–Meier plot that low SFRP2 gene methylation was associated with poorer survival when compared to high methylation status (p = 0.027) in non-treated CRC patients (Fig. 2I).

Table 2 Association between promoter SFRP2 methylation and clinicopathological variables under linear and logistic regression

Furthermore, in order to find other variables that are related to SFRP2 methylation, we conducted a logistic regression, taking SFRP2 methylation as at low and high values. We found in the Table 2 that higher promoter SFRP2 methylation were more likely to be located in the colon (OR = 0.20, 95% CI 0.06–0.59, p = 0.005), more likely to be higher in those patients that did not have lymph node invasion (OR = 0.36, 95% CI 0.12–0.99, p = 0.05), more likely to be higher in those patients that did not receive neoadjuvant treatment (OR = 0.05, 95% CI 0.02–0.20, p < 0.001) and more likely to be higher in those patients in early stage (I + II stages) (OR = 0.31, 95% CI 0.10–0.86, p = 0.040).

SFRP2 decreases cell proliferation and migration in colorectal cancer cell

We next studied SFRP2 in a cell model, using HCT116 cells. We analyzed the promoter methylation in these cells. We found that the promoter of SFRP2 gene in HCT116 was fully methylated at 7 CpG sites that we analyzed. The percentages of methylation were 93 ± 3.7% in site 1; 89 ± 2.5% in site 2; 97 ± 2.5% in site 3; 99 ± 1.5% in site 4; 92 ± 4.3% in site 5; 98 ± 2.5% in site 6 and 97 ± 3.6% in site 7, being the mean of all CpG sites 95% (Fig. 3A). Then, we studied whether this fully methylated promoter could affect SFRP2 expression. For this, we treated HCT116 cells with AZA, an DNMT inhibitor. Before AZA treatment, we did detect SFRP2 expression in HCT116 cells. However, after AZA treatment, we found that SFRP2 expression was recovered at a dose-dependent effect, increasing at 1, 5 and 10 µM, respectively (p < 0.001) (Fig. 3B).

Fig. 3
figure 3

SFRP2 inhibit cell proliferation and migration, through decreasing AXIN2 expression. A promoter methylation of SFRP2 in HCT116 cells. B Cells were seeded in 6-well plates at the concentration of 50.000 cell/well. Cells were treated with 1, 5 and 10 µM of AZA (5-Aza-2-deoxycytidine), during 72 h. SFRP2 gene expression was measured using qPCR and Probe qPCR kit (Takara, Spain). Gene normalization was used by PPIA. C Cells were treated at crescent concentration of rhSFRP2 at 0.1, 0.5, 1, 2.5 and 5 µg/mL of rhSFRP2 during 48 h. D wound healing Assay was conducted in the presence of 5 µg/mL of rhSFRP2. EH gene expression in cell treated with vehicle or with 5 µg/mL of rhSFRP2 of AXIN2, CTNNB1, TIAM1 and VEGF. AZA azacytidine, CTNNB1 beta-catenin, ND No Detected, PPIA Peptidylprolyl Isomerase A, rhSFRP2 recombinant human secreted frizzled related Protein 2, TIAM1 T-cell lymphoma invasion and metastasis-inducing protein 1, VEGR Vascular endothelial growth factor

As SFRP2 is not expressed in HCT116, we decided to treat the cells with recombinant human SFRP2 (rhSFRP2). After 48 h of treatment, we observed that cell growth was maintained al 100% of cell viability at the concentrations of 0.1, 0.5 and 1 µg/mL of rhSFRP2. At 2.5 µg/mL of rhSFRP2, the cell survival decreased up to 81%, and up to 66% when we treated cells with 5 µg/mL of rhSFRP2, with a significant difference between 5 µg/mL and vehicle (p < 0.001). Then, we decided to use 5 µg/mL of rhSFRP2 concentration for the following experiments (Fig. 3C). Next, we evaluated the role of rhSFRP2 in cell migration. At 24 h after scratching cells, we did not observe a significant difference between cells treated with rhSFRP2 and vehicle. However, after 48 h, we observed that cells migrated less when they were treated with 5 µg/mL of rhSFRP2 (p = 0.017) (Fig. 3D). Furthermore, we found that treatment with rhSFRP2 decreased the expression of AXIN2 (p < 0.01) (Fig. 3F). No effect was observed in the CTNNB1, TIAM1 and VEGF genes (Fig. 3G–I).

Discussion

In our study, we characterized the expression and DNA methylation levels of SFRP2 in whole blood, adipose tissue, and colonic tissues, considering different cancer outcomes. We found that SFRP2 promoter methylation in whole blood is a promising prognostic factor for several cancer outcomes, such as stages, lymph node invasion and CRC recurrence, while SFRP2 promoter methylation in tumors is a potential predictor of treatment response. Additionally, treatment with rhSFRP2 was able to decrease cell proliferation, likely due to repression of AXIN2 in the Wnt signaling pathway, indicating its role as a therapeutic tool. The integration of in vitro experiments, clinical data, and bioinformatics analyses provides a comprehensive approach, regarding molecular mechanisms, potential therapeutic targets, and predictive and therapeutic biomarkers. This integration enhances our understanding of the role of SFRP2 in the CRC context. Consequently, SFRP2 is a potential tumor suppressor gene in CRC, which could be recommended as a therapeutic target for cancer clinical practice.

Several studies have evaluated the use of SFRP2 methylation as a non-invasive biomarker. Using blood as a non-invasive clinical tool, these studies found that SFRP2 methylation was predominantly hypermethylated in blood samples [27], including circulating tumor cells (CTCs) [28], in serum [29], circulating free DNA [30, 31], and plasma [32, 33]. However, contradictory results were observed, showing low sensitivity and specificity, mainly because most of these studies used the bisulfite PCR method to evaluate the methylation status of SFRP2. In our study, utilizing specific SFRP2 promoter DNA pyrosequencing, we did not observe a significant difference in whole blood between healthy participants and patients. This finding contrasts with that of a meta-analysis [34], underscoring the importance of increasing the sample size to draw definitive conclusions regarding its utility as diagnostic biomarkers. However, we found that low SFRP2 methylation was observed in patients with late-stage disease, lymph node invasion, and recurrence. Similarly, Zhang et al. (2015) found that SFRP2 methylation in plasma was significantly correlated with TNM stage and lymph node invasion, as observed in our study [33], indicating that SFRP2 methylation may be a good prognostic biomarker for specific clinical outcomes. In contrast, another study conducted by Bartak et al. (2022) did not find a significant difference in promoter SFRP2 in cell-free DNA between patients with stable disease and those with recurrence or remission [31]. However, further studies are needed to establish a standardized method to measure SFRP2 methylation, involving a large number of participants.

On the other hands, multiple studies have reported that SFRP2 is hypermethylated in tumor tissue from CRC patients [11, 35, 36]. However, the results have shown higher variability. In our study, we determined the global β methylation of the SFRP2 gene in the tumor and NAT area using DNA array. We found that SFRP2 methylation in the tumor area was significantly higher than that in the NAT area, but not for SFRP2 expression, probably due to the sample size. These results were validated using the TCGA-COAD and READ cohorts, and by studying the promoter methylation as reported by our group in a previous study [11]. However, our results did not show that SFRP2 methylation in the tumor had prognostic value. Controversial results were observed; a study performed by Liu et al. (2019) found that SFRP2 hypomethylation in the tumor had significantly shorter survival than those with SFRP2 hypermethylation [36]. In contrast, a study performed by Tang et al. (2011) found that SFRP2 hypermethylation in the tumor had significantly shorter survival than those negative for SFRP2 methylation and was an independent prognostic factor [10]. Therefore, further studies are needed to clarify these observations.

We also found that SFRP2 methylation in tumors is very sensitive to chemotherapy, which could be a potential prognostic factor for therapy monitoring. Treated patients showed lower SFRP2 methylation than those who were not treated. Indeed, untreated patients with lower SFRP2 methylation had poor overall survival compared to those with high SFRP2 methylation. A study conducted by Singh et al. (2018) found that SFRP2 methylation in peripheral blood plasma decreased in those patients after chemotherapy [37]. In addition, another study found that SFRP2 was downregulated in liver metastasis CRC patients [38], indicating a sensitive biomarker for treatment response. However, further studies are needed to understand which variables affect SFRP2 methylation. In our logistic and linear regression analysis, we found that SFRP2 methylation could be affected by therapy, stage, lymph node invasion, as well as location. Accordingly, we found that SFRP2 methylation was higher in the colon compared to the rectum, as observed by Bagci et al. [39]. Therefore, several clinical outcomes should be taken into account when evaluating the diagnostic and prognostic factors of the SFRP2 gene in CRC. Given the variability in SFRP2 methylation status observed across studies and tissues, it is crucial to recognize that this variation may be attributed to tissue type and population differences. Therefore, considering multiple confounding biomarkers is imperative to prevent misinterpretation.’ This addition underscores the importance of considering various factors when interpreting SFRP2 methylation data.

SFRP2 is a member of the SFRP family of proteins involved in the regulation of the Wnt signaling pathway as an antagonist of Wnt proteins. The Wnt pathway is critical in many biological processes, including cancer. Specifically, SFRP2 has a complex role in the development and progression of CRC [40]. Since SFRP2 is primarily expressed in stromal cells and silenced in tumors, as scRNASeq analysis showed, this phenomenon can elucidate its role as a tumor suppressor gene. In our study, we found that SFRP2 was silenced in HCT116, and AZA treatment restored SFRP2 expression in a dose-dependent manner. It is worth to note that DNMT inhibition (through AZA) leads to hypomethylation of CpG sites, resulting in increased accessibility of the gene promoters to transcription factors and other regulatory proteins. This increased accessibility can lead to enhanced gene expression of previously silenced genes. Ectopic rhSFRP2 decreased cell proliferation and migration, probably through decreasing AXIN2 expression, but not for CTNNB1, TIAM1 and VEGF, mainly due to specific action of SFRP2 protein as a Wnt antagonist. AXIN2 act as negative regulators of the Wnt signaling pathway by promoting the phosphorylation and degradation of β-catenin. Several studies have demonstrated similar findings in different cancer contexts, such as in MG63 cells (osteosarcoma) [41], SW480 cells (colorectal cell line) [42], At-T20 cells (mouse pituitary tumor cells) [43], as well as CCLP-1 and QBC939 cells (Cholangiocarcinoma).

This effect was mainly observed through the increase in the expression of phosphorylated-β-catenin and phosphorylated-GSK3β in stem cells of apical papilla [44], and melanocytes [45], as well as decreased AXIN2 expression [44], as observed in our study. Another study found that treatment with rhSFRP2 in MG63 cells increased RhoA activity, which decreased the effects induced by Wnt5A [41]. Indeed, rhSFRP2 decreased RhoA activity induced by Wnt5A in both U251 and T98MG cells (glioblastoma) [46]. Furthermore, a study conducted by Kim et al. (2018) showed that SFRP2 can act through the canonical Wnt/β-catenin pathway in an independent manner in iPSC-derived MSCs [47], suggesting new roles in the Wnt pathway. Overall, the effects of SFRP2 overexpression in HCT116 cells suggest that this protein may have a tumor-suppressive function in colorectal cancer. rhSFRP2 could act as therapeutic target, through its ability to attenuate the Wnt pathway. Given the predominant silencing of SFRP2 observed in colorectal tumors, indicative of its role as a tumor suppressor gene, the rhSFRP2 model provides valuable insights into its impact on pivotal pathways. Specifically, our findings shed light on SFRP2's modulation of the Wnt pathway through its interaction with the AXIN2 gene, presenting novel therapeutic possibilities. However, the role of SFRP2 in cancer cells is multifaceted and likely contingent upon the specific cellular environment and context. Our study highlights several biomedical applications of SFRP2 in the context of CRC. SFRP2 methylation patterns can be used in early detection, disease progression assessment, and therapy monitoring. Further insights into disease mechanisms can be gained by deepening the characterization of SFRP2 in mice and other in vitro models, thereby uncovering the therapeutic potential of rhSFRP2. Further studies are needed to fully elucidate the mechanisms underlying these effects and to determine the potential clinical implications of SFRP2 as a therapeutic target in cancer therapy [48, 49].

There are several limitations that worth to mention related this study. Firstly, this research is based on an observational and prospective study with limited sample size, that need to be validated in future research. Also, longitudinal studies are needed, particularly in standardizing methods for measuring SFRP2 methylation and understanding its variable roles in cancer. Further molecular and mechanistic research on SFRP2 is needed to provide a comprehensive understanding on its role in colorectal carcinogenesis. Additionally, the use of specific SFRP2 promoter DNA pyrosequencing and Infinium Human Methylation 450 K bead chip technology offers several advantages, including high sensitivity and specificity in detecting DNA methylation patterns, when compared to other methods such as Methyl-Specific PCR. However, techniques like bisulfite-DNA sequencing provide extensive coverage, capturing CpG sites not included in this study. It is important to acknowledge potential limitations such as the requirement for specialized equipment and expertise, as well as the necessity for rigorous result validation.

Conclusion

Studies have shown that SFRP2 plays a complex role in the development and progression of CRC. In our study, we found that SFRP2 methylation in whole blood is a potential biomarker for several cancer outcomes, while SFRP2 methylation in tumor tissue is a promising predictor for therapy response, suggesting its potential use in clinical practice. In a CRC cell model, SFRP2 was mainly found to be suppressed, and its expression was restored by DNA demethylation, indicating epigenetic silencing. Treatment with rhSFRP2 was able to decrease cell proliferation and migration, suggesting a role as a tumor suppressor gene, although previous steps in other biological models are needed. However, the exact mechanisms by which SFRP2 influences CRC development and progression are still not fully understood. Nonetheless, our findings suggest that SFRP2 may be a promising target for the development of novel therapeutic approaches for CRC.

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available due the participants were anonymized but are available from the corresponding author on reasonable request.

References

  1. Sung H, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.

    Article  CAS  PubMed  Google Scholar 

  2. Jiang H, Zhou S, Li G. Novel biomarkers used for early diagnosis and tyrosine kinase inhibitors as targeted therapies in colorectal cancer. Front Pharmacol. 2023. https://doi.org/10.3389/fphar.2023.1189799.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Petrelli F, et al. Prognostic survival associated with left-sided vs right-sided colon cancer. JAMA Oncol. 2017;3(2):211. https://doi.org/10.1001/jamaoncol.2016.4227.

    Article  PubMed  Google Scholar 

  4. You JS, Jones PA. Cancer genetics and epigenetics: two sides of the same coin? Cancer Cell. 2012;22(1):9–20. https://doi.org/10.1016/j.ccr.2012.06.008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Sadikovic B, Al-Romaih K, Squire J, Zielenska M. Cause and consequences of genetic and epigenetic alterations in human cancer. Curr Genomics. 2008;9(6):394–408. https://doi.org/10.2174/138920208785699580.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Savio AJ, Mrkonjic M, Lemire M, Gallinger S, Knight JA, Bapat B. The dynamic DNA methylation landscape of the mutL homolog 1 shore is altered by mlh1-93g>a polymorphism in normal tissues and colorectal cancer. Clin Epigenetics. 2017;9(1):26. https://doi.org/10.1186/s13148-017-0326-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Liu J, et al. Wnt/β-catenin signalling: function, biological mechanisms, and therapeutic opportunities. Signal Transduct Target Ther. 2022;7(1):3. https://doi.org/10.1038/s41392-021-00762-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Yu JN, et al. Association between SFRP promoter hypermethylation and different types of cancer: a systematic review and meta-analysis. Oncol Lett. 2019;18(4):3481. https://doi.org/10.3892/OL.2019.10709.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Shariatpanahi AM, Yassi M, Nouraie M, Sahebkar A, Tabrizi FV, Kerachian MA. The importance of stool DNA methylation in colorectal cancer diagnosis: a meta-analysis. PLoS ONE. 2018. https://doi.org/10.1371/JOURNAL.PONE.0200735.

    Article  Google Scholar 

  10. Tang D, Liu J, Wang DR, Yu HF, Li YK, Zhang JQ. Diagnostic and prognostic value of the methylation status of secreted frizzled-related protein 2 in colorectal cancer. Clin Invest Med. 2011. https://doi.org/10.25011/CIM.V34I1.15105.

    Article  PubMed  Google Scholar 

  11. Cabrera-Mulero A, et al. Novel SFRP2 DNA methylation profile following neoadjuvant therapy in colorectal cancer patients with different grades of BMI. J Clin Med. 2019;8(7):1041. https://doi.org/10.3390/jcm8071041.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Carr NJ, Robin LH. WHO Classification of tumors of the digestive system. 4th ed. Geneva: World Health Organization; 2010.

    Google Scholar 

  13. Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972. https://doi.org/10.1177/107424840501000106.

    Article  PubMed  Google Scholar 

  14. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and β-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985. https://doi.org/10.1007/BF00280883.

    Article  PubMed  Google Scholar 

  15. Izquierdo AG, et al. DNA methylome in visceral adipose tissue can discriminate patients with and without colorectal cancer. Epigenetics. 2021. https://doi.org/10.1080/15592294.2021.1950991.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

  17. “Protocol Mycoplasma Detection Test.”.https://projects.iq.harvard.edu/files/hlalab/files/mycoplasm-test_hla.pdf. Accessed: 23 Nov 2022

  18. Molla Kazemiha V, et al. PCR-based detection and eradication of mycoplasmal infections from various mammalian cell lines: a local experience. Cytotechnology. 2009;61(3):117. https://doi.org/10.1007/S10616-010-9252-6.

    Article  CAS  PubMed  Google Scholar 

  19. de Ieso ML and Pei JV. An accurate and cost-effective alternative method for measuring cell migration with the circular wound closure assay. Biosci Rep. 2018;38(5).10.1042/BSR20180698

  20. Aryee MJ, et al. Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014. https://doi.org/10.1093/bioinformatics/btu049.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Tian Y, et al. ChAMP: updated methylation analysis pipeline for illumina beadchips. Bioinformatics. 2017. https://doi.org/10.1093/bioinformatics/btx513.

    Article  PubMed  PubMed Central  Google Scholar 

  22. “The Cancer Genome Atlas - Citing TCGA - NCI.” https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga/using-tcga/citing-tcga. Accessed: 28 Nov 2022

  23. Colaprico A, et al. TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 2016;44(8):e71–e71. https://doi.org/10.1093/NAR/GKV1507.

    Article  PubMed  Google Scholar 

  24. Tang G, Cho M, Wang X. OncoDB: an interactive online database for analysis of gene expression and viral infection in cancer. Nucleic Acids Res. 2022;50(D1):D1334–9. https://doi.org/10.1093/NAR/GKAB970.

    Article  CAS  PubMed  Google Scholar 

  25. Lee HO, et al. Lineage-dependent gene expression programs influence the immune landscape of colorectal cancer. Nat Genet. 2020;526:594–603. https://doi.org/10.1038/s41588-020-0636-z.

    Article  CAS  Google Scholar 

  26. RStudio Team. RStudio: Integrated Development for R. RStudio. Boston: RStudio Inc; 2015.

    Google Scholar 

  27. Sui C, Ma J, Chen Q, Yang Y. The variation trends of SFRP2 methylation of tissue, feces, and blood detection in colorectal cancer development. Eur J Cancer Prev. 2016;25(4):288–98. https://doi.org/10.1097/CEJ.0000000000000185.

    Article  CAS  PubMed  Google Scholar 

  28. Lyberopoulou A, et al. Identification of methylation profiles of cancer-related genes in circulating tumor cells population. Anticancer Res. 2017;37(3):1105–12. https://doi.org/10.21873/ANTICANRES.11423.

    Article  CAS  PubMed  Google Scholar 

  29. Li H, et al. Performance of a methylight assay for methylated sfrp2 dna detection in colorectal cancer tissue and serum. Int J Biol Markers. 2019;34(1):54–9. https://doi.org/10.1177/1724600818820536.

    Article  CAS  PubMed  Google Scholar 

  30. Barták BK, et al. Blood collection and cell-free dna isolation methods influence the sensitivity of liquid biopsy analysis for colorectal cancer detection. Pathol Oncol Res. 2019;25(3):915–23. https://doi.org/10.1007/S12253-018-0382-Z.

    Article  PubMed  Google Scholar 

  31. Barták BK, et al. A liquid biopsy-based approach for monitoring treatment response in post-operative colorectal cancer patients. Int J Mol Sci. 2022;23:7. https://doi.org/10.3390/IJMS23073774/S1.

    Article  Google Scholar 

  32. Zhao G, et al. A novel plasma based early colorectal cancer screening assay base on methylated SDC2 and SFRP2. Clin Chim Acta. 2020;503:84–9. https://doi.org/10.1016/J.CCA.2020.01.010.

    Article  CAS  PubMed  Google Scholar 

  33. Zhang X, et al. Combined detection of plasma GATA5 and SFRP2 methylation is a valid noninvasive biomarker for colorectal cancer and adenomas. World J Gastroenterol: WJG. 2015;21(9):2629. https://doi.org/10.3748/WJG.V21.I9.2629.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Yang Q, Huang T, Ye G, Wang B, Zhang X. Methylation of SFRP2 gene as a promising noninvasive biomarker using feces in colorectal cancer diagnosis: a systematic meta-analysis. Sci Rep. 2016;6(1):1–11. https://doi.org/10.1038/srep33339.

    Article  CAS  Google Scholar 

  35. Hu H, et al. Hypermethylated promoters of secreted frizzled-related protein genes are associated with colorectal cancer. Pathol Oncol Res. 2019;25(2):567–75. https://doi.org/10.1007/s12253-018-0505-6.

    Article  CAS  PubMed  Google Scholar 

  36. Liu X, et al. DNA methylation of SFRP1, SFRP2, and WIF1 and prognosis of postoperative colorectal cancer patients. BMC Cancer. 2019. https://doi.org/10.1186/s12885-019-6436-0.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Bhangu JS, et al. Circulating free methylated tumor dna markers for sensitive assessment of tumor burden and early response monitoring in patients receiving systemic chemotherapy for colorectal cancer liver metastasis. Ann Surg. 2018;268(5):894–902. https://doi.org/10.1097/SLA.0000000000002901.

    Article  PubMed  Google Scholar 

  38. Kumbrink J, et al. Serial analysis of gene mutations and gene expression during first-line chemotherapy against metastatic colorectal cancer: identification of potentially actionable targets within the multicenter prospective biomarker study reveal. Cancers (Basel). 2022;14:15. https://doi.org/10.3390/CANCERS14153631/S1.

    Article  Google Scholar 

  39. Bagci B, Sari M, Karadayi K, Turan M, Ozdemir O, Bagci G. KRAS, BRAF oncogene mutations and tissue specific promoter hypermethylation of tumor suppressor SFRP2, DAPK1, MGMT, HIC1 and p16 genes in colorectal cancer patients. Cancer Biomark. 2016;17(2):133–43. https://doi.org/10.3233/CBM-160624.

    Article  CAS  PubMed  Google Scholar 

  40. Liu Y, Zhou Q, Zhou D, Huang C, Meng X, Li J. Secreted frizzled-related protein 2-mediated cancer events: friend or foe? Pharmacol Rep. 2017;69(3):403–8. https://doi.org/10.1016/J.PHAREP.2017.01.001.

    Article  CAS  PubMed  Google Scholar 

  41. Zhang A, Yan T, Wang K, Huang Z, Liu J. PI3Kα isoform-dependent activation of RhoA regulates Wnt5a-induced osteosarcoma cell migration. Cancer Cell Int. 2017;17(1):27. https://doi.org/10.1186/S12935-017-0396-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Bai Z, Wu Q, Zhang C, Chen J, Cao L. Effects of YAP1 and SFRP2 overexpression on the biological behavior of colorectal cancer cells and their molecular mechanisms. J Gastrointest Oncol. 2021;12(4):1601. https://doi.org/10.21037/JGO-21-418.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Ren J, et al. Decreased expression of SFRP2 promotes development of the pituitary corticotroph adenoma by upregulating Wnt signaling. Int J Oncol. 2018;52(6):1934. https://doi.org/10.3892/IJO.2018.4355.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Jin L, et al. SFRP2 enhances the osteogenic differentiation of apical papilla stem cells by antagonizing the canonical WNT pathway. Cell Mol Biol Lett. 2017;22(1):1–14. https://doi.org/10.1186/S11658-017-0044-2/FIGURES/6.

    Article  Google Scholar 

  45. Kim M, Han JH, Kim JH, Park TJ, Kang HY. Secreted frizzled-related protein 2 (sfrp2) functions as a melanogenic stimulator; the role of SFRP2 in UV-induced hyperpigmentary disorders. J Investig Dermatol. 2016;136(1):236–44. https://doi.org/10.1038/JID.2015.365.

    Article  CAS  PubMed  Google Scholar 

  46. Liu G, et al. Daam1 activates RhoA to regulate Wnt5a-induced glioblastoma cell invasion. Oncol Rep. 2018;39(2):465. https://doi.org/10.3892/OR.2017.6124.

    Article  CAS  PubMed  Google Scholar 

  47. Kim H, et al. Oncogenic role of SFRP2 in p53-mutant osteosarcoma development via autocrine and paracrine mechanism. Proc Natl Acad Sci U S A. 2018;115(47):E11128–37. https://doi.org/10.1073/PNAS.1814044115/-/DCSUPPLEMENTAL.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Khalilov R, Bakishzade A, Nasibova A. Future prospects of biomaterials in nanomedicine. Adv Biol Earth Sci. 2024;9:5–10. https://doi.org/10.62476/abes.9s5.

    Article  Google Scholar 

  49. Rosic G, Selakovic D, Omarova S. Cancer signaling, cell/gene therapy, diagnosis and role of nanobiomaterials. Adv Biol Earth Sci. 2024;9:11–34. https://doi.org/10.62476/abes9s11.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Diana Garcia from IISLaFe for her technical support.

Funding

This study was supported by “Centro de Investigacion Biomédica en Red Fisiopatología de la Obesidad y Nutricion”, which is an initiative of the “Instituto de Salud Carlos III” (ISCIII) of Spain, financed by the European Regional Development Fund, “A way to make Europe”/ “Investing in your future” (CB06/03), a grant from ISCIII (PI18/01399, PI21/00633) and a grant from Consejeria Universidad, Investigacion e Innovacion Junta de Andalucia (PY20-01270). HB is supported by a predoctoral fellowship (“Plan Propio IBIMA 2020 A.1 Contratos predoctorales”, Ref.: predoc20_002) and by a “Sara Borrell” postdoctoral contract (CD22/00053) from the Instituto de Salud Carlos III—Madrid (Spain), “Financiado por la Unión Europea – NextGenerationEU” y mediante el Plan de Recuperación, Transformación y Resiliencia. LAGF was supported by a “Sara Borrell” postdoctoral contract (CD21/000131) from the Instituto de Salud Carlos III—Madrid (Spain). M.M.G. was the recipient of the Nicolas Monardes Programme from the “Servicio Andaluz de Salud, Junta de Andalucia”, Spain (RC-0001–2018 and C-0029–2014).

Author information

Authors and Affiliations

Authors

Contributions

HB, IAS and MMG designed research studies. JS and MMG supervised the study. HB, JP, LAGF conducted experiments. MOC and MRF conducted follow-up of patients with CRC. HB, JP provided statistics analysis. HB and JP acquired data. HB, JS and MMG wrote the manuscript. HB, JS and LAGF and MMG analyzed data. All authors revised the manuscript.

Corresponding authors

Correspondence to Isabel Arranz or Manuel Macias-Gonzalez.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all patients and subjects and reviewed and approved by the Ethics Committees of the “Virgen de la Victoria” University Hospital (Málaga, Spain) (Registration Number: 0311/PI7). The study was conducted in accordance with the guidelines laid down in the Declaration of Helsinki.

Consent for publication

All participants signed a informed consent.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

12967_2024_5329_MOESM1_ESM.pdf

Supplementary material 1: Supplementary Figure 1. We measured the normalized A) SFRP2 expression in blood sample from healthy participants (N=28) and patients with CRC (N=27). B) Kaplan-Meier curve comparing the median of the SFRP2 promoter methylation as low and high methylation. C) We measured the promoter methylation of SFRP2 gene from healthy participants and patients with CRC in early stage (I+II) and late stage (III+IV). D) SFRP2 expression in adipose tissue healthy participants (N=54) and patients with CRC (N=64) and F) SFRP2 methylation in adipose tissue healthy participants (N=56) and patients with CRC (N=41), Kaplan-Meier curve comparing the median of the E) SFRP2 expression in whole blood as low and high expression, SFRP2 expression in adipose tissue as low and high expression and G) promoter SFRP2 methylation in adipose tissue as low and high expression. The significance of differences is evaluated with the Log-rank test. Asterisks indicate significant differences between the groups according to the Mann Whitney test (*p<0.05, **p<0.01, ***p<0.001).

12967_2024_5329_MOESM2_ESM.pdf

Supplementary material 2 A) Kaplan-Meier curve comparing the median of the SFRP2 global β methylation as low and high methylation in A) our population, C) TCGA-COAD and D) TCGA-READ, as well as by studying B) promoter methylation. The significance of differences is evaluated with the Log-rank test. E) Promoter SFRP2 methylation in in tumor by comparing colon (N=30) vs rectum (N=51). Asterisks indicate significant differences between the groups according to the Mann Whitney test (*p<0.05, **p<0.01, ***p<0.001).

12967_2024_5329_MOESM3_ESM.docx

Supplementary material 3: Table 1. Colorectal cancer outcome prediction using SFRP2 methylation as predictive biomarker performed by ROC curve analysis

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

Boughanem, H., pilo, J., García-Flores, L.A. et al. Identification of epigenetic silencing of the SFRP2 gene in colorectal cancer as a clinical biomarker and molecular significance. J Transl Med 22, 509 (2024). https://doi.org/10.1186/s12967-024-05329-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-024-05329-x

Keywords