SCG2 is a Prognostic Biomarker Associated With Immune Infiltration and Macrophage Polarization in Colorectal Cancer

Colorectal cancer (CRC) is the second most lethal malignancy around the world. Limited efficacy of immunotherapy creates an urgent need for development of novel treatment targets. Secretogranin II (SCG2) is a member of the chromogranin family of acidic secretory proteins, has a role in tumor microenvironment (TME) of lung adenocarcinoma and bladder cancer. Besides, SCG2 is a stroma-related gene in CRC, its potential function in regulating tumor immune infiltration of CRC needs to be fully elucidated. In this study, we used western blot, real-time PCR, immunofluorescence and public databases to evaluate SCG2 expression levels and distribution. Survival analysis and functional enrichment analysis were performed. We examined TME and tumor infiltrating immune cells using ESTIMATE and CIBERSORT algorithm. The results showed that SCG2 expression was significantly decreased in CRC tumor tissues, and differentially distributed between tumor and adjacent normal tissues. SCG2 was an independent prognostic predictor in CRC. High expression of SCG2 correlated with poor survival and advanced clinical stage in CRC patients. SCG2 might regulate multiple tumor- and immune-related pathways in CRC, influence tumor immunity by regulating infiltration of immune cells and macrophage polarization in CRC.

Introduction

Colorectal cancer (CRC) is the fourth most commonly diagnosed cancer, causing more than ninety thousand deaths worldwide every year (Bray et al., 2018; Siegel et al., 2020). Although diagnosis and treatment have been improved substantially, the prognosis of patients remains poor, especially in those with higher TNM stage (Luebeck et al., 2013; Miller et al., 20162016). Abnormal expression of multiple genes is usually associated with the occurrence of CRC (Meyerson et al., 2010). However, the molecular mechanisms underlying CRC remain unclear. It is urgent to identify novel diagnostic biomarkers and prognostic predictors for this disease.

New treatment modalities have been proposed for CRC, such as immunotherapy (Overman et al., 2018). Recently, studies have shown that the tumor microenvironment (TME) significantly affects CRC progression and therapeutic efficacy (Halama et al., 2011). Monoclonal antibodies against programmed cell death 1 (PD-1), programmed death ligand-1 (PD-L1) and cytotoxic T-lymphocyte–associated antigen 4 (CTLA4) have been proven effective in clinical trials (Rotte, 2019). However, anti-PD-1/PD-L1 therapy is unsatisfactory in metastatic CRC (Le et al., 2015), whereas anti-CTLA-4 treatment showed a limited response rate in advanced CRC (Chung et al., 2010). Therefore, it is of paramount importance to find novel immunotherapy targets of CRC.

Secretogranin II (SCG2) is a member of the tyrosine-sulphated granin family expressed in endocrine, neuroendocrine and neuronal tissues (Troger et al., 2017). It is relevant to secretory vesicle formation and packaging peptide hormones into vesicles (Beuret et al., 2004). SCG2 has an important role in enhancing endothelial cell proliferation, migration, and angiogenesis (Albrecht-Schgoer et al., 2012; Hannon et al., 2018). Studies have revealed that the derived peptides of SCG2, such as secretoneurin (SN) and EM66, are useful markers of neuroendocrine tumors (Guillemot et al., 2012). Besides, SCG2 could predict prognosis in Non-small cell lung cancer as a potentially secreted biomarker (Cury et al., 2019), reflect the TME of lung adenocarcinoma (Wang and Chen, 2020) and bladder cancer (Luo et al., 2020). Recently, Liu et al. (2020) identified SCG2 was a stroma-related gene and predicted poor outcomes in CRC patients. However, the relationship between SCG2 and tumor immunity of CRC is largely unknown, and the mechanisms underlying it remain to be intensively investigated.

This study analyzed SCG2 expression profiles and described its potential prognostic value, multiple biological functions, related signal pathways in CRC. This report’s findings revealed the potential role SCG2 played in regulating tumor immunity and provided a possible mechanism based on bioinformatic analysis.

Materials and Methods

Specimen Collection

We obtained 61 pairs of CRC primary tumor tissues and adjacent tissues from patients who underwent surgery at Zhongnan Hospital of Wuhan University (Wuhan, China) between June 2019 and September 2019. None of the patients received any preoperative therapy such as adjuvant chemotherapy or radiotherapy. Samples of the collected tissues were preserved in RNAlater Stabilization solution (Invitrogen, United States) and stored in the department of Biological Repositories (Zhongnan Hospital of Wuhan University). Our study was approved by Zhongnan Hospital Ethics Committee (ethics code 2018025), and every patient signed informed consent.

RNA Extraction and qRT-PCR

Total RNA was extracted using Trizol reagent (Invitrogen, United States), and cDNA was reverse transcribed from 1ug of purified RNA using PrimeScriptTM RT reagent kit (Toyobo, Osaka). The quantification of mRNA was examined using quantitative real-time polymerase chain reaction (qRT-PCR) on Biorad CFX (Biorad, United States), and gene expression was normalized by GAPDH. The primer sequences were obtained from PrimerBank and synthesized by TSINGKE (Wuhan, China) as follows: GAPDHF-5′CTGGGCTACACTGAGCACC3′, R-5′AAGTGGTCGTTGAGGGCAATG3′, SCG2F-5′ACCAGACCTCAGGTTGGAAAA3′, R-5′ACCAGACCTCAGGTTGGAAAA3′. We used the comparative CT (2-ΔΔCT) to calculate the gene mRNA expression levels, and the experiment was repeated three times with three biological replicates for each treatment.

Protein Extraction and Western Blotting

According to the reagent instructions, total protein in tissues was extracted with RIPA lysis buffer (Boyotime, China), including a protease inhibitor cocktail (Thermo Scientific, United States). We performed western blot using standard methods with the specific antibody, SCG2 (1:1,000, Absin, abs117155), GAPDH (1:10000, Proteintech, 60004-1-Ig).

Immunofluorescence

Tumor samples were immediately fixed with 4% Paraformaldehyde (PFA) for 1 h and dehydrated overnight at 4°C. We used the following primary antibodies for immunofluorescence: mouse anti-SCG2 (1:100, Absin, abs117155) and rat anti-CD68 (1:100, Novusbio, 100–683). Fluorescent secondary antibodies CY3-conjugated anti-mouse SCG2 (1:1,000) and FITC-conjugated anti-rat CD68 antibody (1:1,000) were incubated for 30 min at 25°C. Nuclei were counter stained with DAPI for 10 min. We randomly chose at least five optical fields (20× or 40× magnification) per tumor section for morphometric evaluation.

Data Acquisition

RNA-sequencing expression data and clinicopathological information of 469 CRC patients were downloaded from The Cancer Genome Atlas (TCGA) database. GSE39582 cohort containing microarray gene expression profiles and clinical data was downloaded from the Gene Expression Omnibus (GEO) database. The proteome data of CRC patients was obtained from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) Assay Portal (Whiteaker et al., 2014).

SCG2 Expression Analysis

We identified SCG2 mRNA expression levels in multiple human cancers from the Tumor Immune Estimation Resource (TIMER) database (Li et al., 2017). In TCGA and GEO cohorts, the mRNA expression profiles of SCG2 were visualized, respectively. In the CPTAC cohort, the protein expression levels of SCG2 in 91 pairs of tumor and adjacent tissues were examined. qRT-PCR and WB were utilized to detect RNA and protein levels of SCG2. Finally, immunofluorescence assay was used to detect the expression and distribution of SCG2 protein in tumor and adjacent tissues of CRC patients.

Survival Analysis and Significant Prognostic Marker Analysis

Patients were divided into high-SCG2 expression group and low-SCG2 expression group according to the median expression value of SCG2. We used R packages “survival” and “survminer” to analyze the correlation between SCG2 expression and overall survival (OS) time, disease-free survival (DFS) time. Besides, univariate and multivariate survival analyses were applied to examine the independent prognostic factors in CRC, using the “coxph” function in the “survival” package.

SCG2 Related Genes and Functional Enrichment Analysis

To further study the SCG2-related molecular mechanisms, we performed differential expression genes (DEGs) analysis based on the medium value of SCG2 expression using the package “DESeq2.” Genes were ranked with positive correlation coefficients with SCG2 (Log2FoldChange ≥ 1, p < 0.01). Then Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were conducted using “enrichGo” and “enrichKEGG” functions in the “clusterprofiler” package. We conducted a PPI network based on DEGs by Cytoscape software, and its plug-in “ClueGO” was applied for further function enrichment analysis.

Gene Set Enrichment Analysis

We performed Gene Set Enrichment Analysis (GSEA) to further understand the upregulated pathways in the high SCG2 expression group by GSEA 4.1 software. The number permutations for each gene were set to 1,000, false discovery rate (FDR) < 0.5 and normalized enrichment score (NES) > 1 were set at the cut-off criteria.

Characterization of the TME

Tumor infiltrating immune cells (TIICs) in tumor samples and normal samples of the TCGA cohort were calculated using the CIBERSORT algorithm based on RNA-seq (Newman et al., 2015). The signature TIICs matrix “LM22” was set to 1,000 permutations. Tumor immune single-cell analysis was investigated with the “Dataset” module in Tumor Immune Single Cell Hub (TISCH, http://tisch.comp-genomics.org), which built a Single-cell RNA sequencing (scRNA-seq) atlas of 76 high-quality tumor datasets across 27 cancer types from GEO and ArrayExpress (Sun et al., 2021). Thirteen kinds of major lineage immune cells from GSE146771 cohort were analyzed.

Correlation Analysis of SCG2 and TIICs

We applied the Estimation of STromal and Immune cells in MAlignant Tumours using Expression data (ESTIMATE) algorithm to calculate stromal and immune scores that predicted the levels of infiltrating stromal and immune cells in each patient (Yoshihara et al., 2013) and analyzed the correlations between SCG2 expression and the two scores. We further compared the proportions of twenty-two immune cell phenotypes in SCG2 high or low expression group using CIBERSORT algorithm. The correlations between SCG2 and tumor-associated macrophages (TAMs), M2 gene markers were analyzed to identify the relevance of SCG2 and macrophage polarization. The correlations between SCG2 expression and gene markers of various immune cells including CRC immune checkpoint genes were analyzed by TIMER (Li et al., 2017).

Statistical Analysis

SPSS Statistics 26, GraphPad prism 8, R 4.0.2 were used for statistical analyses. The measurement data were presented as mean ± standard deviation. Statistical differences between two groups were examined by Wilcoxon signed rank test or Mann-Whitney U test. Statistical differences across three or more groups were examined by Kruskal-Wallis H test. Correlations were examined with the Spearman rank correlation test (weakly correlation when R > 0.1, moderate correlation when R > 0.3, strong correlation when R > 0.5). Statistically significant differences were considered when p < 0.05.

Results

Expression Profiles of SCG2 in Human Cancers

To evaluate SCG2 expression in various human cancers, we analyzed TCGA RNA-seq using TIMER database. The SCG2 expression levels were higher in breast invasive carcinoma (BRCA), cholangiocarcinoma (CHOL), head and neck cancer (HNSC), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), lung squamous cell carcinoma (LUSC) compared with adjacent normal tissues. However, SCG2 expression was significantly lower in colon adenocarcinoma (COAD), kidney chromophobe (KICH), prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), stomach adenocarcinoma (STAD), uterine corpus endometrial carcinoma (UCEC) compared with adjacent normal tissues (Figure 1A). These results suggested SCG2 expressed abnormally in various tumor types.

FIGURE 1

FIGURE 1. Expression analysis of SCG2. (A) The SCG2 expression level in multiple cancer types from TCGA database. (B,C) SCG2 mRNA expression levels in normal and tumor tissues of CRC from TCGA and GEO database. (D) SCG2 protein levels in 96 paired CRC tissues from CPTAC. (E) SCG2 mRNA expression level in tumor and adjacent samples. (F) Quantification of SCG2 protein expression in tumor and adjacent samples. (A–C) Examined by Mann-Whitney U test; (D–F) Examined by Wilcoxon signed rank test; *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.

We then assessed SCG2 mRNA expression levels in TCGA and GEO cohorts and found significantly lower SCG2 expression levels in cancer tissues than normal tissues (Figures 1B,C). SCG2 protein expression was also lower in cancer tissues than adjacent tissues in CPTAC database (Figure 1D). To confirm the expression profile data, we examined mRNA and protein levels in tumor and normal tissues of CRC patients using qRT-PCR and Western blot (Figures 1E,F).

The expression and distribution of SCG2 protein in tumor and adjacent tissues were monitored by immunofluorescence microscopy (Figure 2). In normal tissues, SCG2 was mainly distributed in colonic glands, possibly related to intestinal fluid secretion, whereas in tumor tissues, the intestinal mucosal tissues were destroyed and SCG2 expression decreased.

FIGURE 2

FIGURE 2. Immunofluorescence images of SCG2 (red) and DAPI (blue) in tumor and adjacent tissues of CRC patients. Seven pairs of tumors and adjacent tissues were examined, six of the tumor tissues showed deceased SCG2 expression compared to adjacent normal tissues.

Prognostic Significance of SCG2 Expression in CRC

A univariate Cox survival analysis based on TCGA database indicated that TNM stage (p = 3.6e-04), Invasion depth (p = 0.031), Lymph node metastasis (p = 0.002), distant metastasis (p = 4.96e-07), and SCG2 expression (p = 0.008) were significant risk factors for overall survival (Table 1). Survival analysis using the Kaplan-Meier method showed that higher SCG2 expression was significantly correlated with shorter OS (p = 0.009) and DFS (p < 0.001) in TCGA cohorts (Figures 3A,B). SCG2 mRNA expression was associated with advanced TNM stage, T stage, and lymph invasion (Figures 3C–F). The forest plot based on multivariate Cox regression analysis showed that age (p = 0.002), M stage (p = 2.07e-04), and SCG2 expression (p = 0.016) were independent prognosis factors of CRC patients (Figure 3G). The above results were verified in independent cohort GSE39582 (Table 1; Figures 3H–N). Besides, we used the patients from our hospital to analyze the correlation between SCG2 expression and tumor histological grading, there was no significant difference (Supplementary Figure S1A). Overall, elevated SCG2 mRNA expression was significantly associated with poor prognosis and advanced clinicopathological parameters in CRC patients.

TABLE 1

TABLE 1. Univariate Cox regression of prognostic factors in CRC patients.

FIGURE 3

FIGURE 3. High SCG2 expression is associated with worse survival and advanced clinical stage in CRC patients. (A,B) OS, DFS in the TCGA CRC cohort. (C–F) SCG2 expression was higher in patients of advanced T, N, M, TNM stage in the TCGA CRC cohort. (G) Multivariate Cox analysis of patients in the TCGA CRC cohort showing the HRs of different factors. (H,I) OS, DFS in the GSE39582 CRC cohort. (J–M) SCG2 expression was higher in patients of advanced T, N, M, TNM stage in the GSE39582 CRC cohort. (N) Multivariate Cox analysis of patients in the GSE39582 CRC cohort showing the HRs of different factors. HRs, hazard ratios. (C–E,J–L) Examined by Kruskal-Wallis H test; (F,M) Examined by Mann-Whitney U test; *p < 0.05, **p < 0.01, ***p < 0.001.

Analysis of SCG2 Co-expressed Genes

Genes in co-expression modules are always involved in same biological pathways (Stuart et al., 2003) and have disease predictive value (Yang et al., 2014). We performed genetic difference analysis based on SCG2 medium expression level. The heat map displayed differential expressed genes according to the Pearson correlation (Supplementary Figure S1B). We verified the correlations between SCG2 and MYH11, SYNPO2, DDR2, FABP4, TNS1 in TIMER (all R > 0.5, p < 0.0001, shown in Supplementary Figures S1C–G). Furthermore, high expressions of FABP4, DDR2, and TNS1 were associated with worse overall survival in CRC patients (Supplementary Figures S1H–J), MYH11 and TNS1 were associated with advanced clinical stage (Supplementary Figures S1K,L). It was possible that abnormal expression of SCG2 and these genes involved in CRC progression contributing to a worse prognosis.

To explore the potential molecular function of SCG2 in CRC, we performed KEGG pathway enrichment and GO analyses with the SCG2 co-expressed genes. KEGG enrichment analysis showed that these genes were associated with many signaling pathways (Supplementary Table S1), including known cancer-related pathways such as the PI3K-Akt signaling pathway and ECM-receptor interaction pathway. GO annotations further suggested these genes were associated with many biological processes (Supplementary Table S2)

We then performed PPI network analysis based on SCG2 co-expressed genes in TCGA and GSE39582 cohort (Supplementary Figure S2A). 56 seed and cluster genes in the network were selected for further enrichment analysis (Supplementary Figure S2B). The result indicated that SCG2 and its co-expressed genes could directly or indirectly regulate various immune-related processes, including macrophage polarization (Supplementary Figure S2C).

GSEA Validates SCG2-Related Pathways

To further study the SCG2-associated signaling pathways in CRC, GSEA was performed between samples in low and high SCG2 expression groups. GSEA identified a total of 20 considerably upregulated hallmark pathways in the high SCG2 expression group (Supplementary Table S3). The tumor related pathways involved in tumorigenesis, tumor invasion, and metastasis included “Kras signaling up,” “Hedgehog signaling,” “Notch signaling,” “Wnt beta catenin signaling,” “TGF beta signaling,” “Epithelial mesenchymal transition,” and “Angiogenesis.” The inflammatory and immune related pathways included “Complement,” “IL2 STAT5 signaling,” “Inflammatory response,” “IL6 JAK STAT3 signaling,” and “Allograft rejection.” A summary of the enrichment results is shown in Supplementary Figure S3.

Evaluation of TIICs in CRC

TIIC is an integral part of tumor immune microenvironment (Deschoolmeester et al., 2010). In the current study, we estimated the proportion of twenty-two immune cells in the tumor and normal samples of TCGA using the CIBERSORT algorithm (Figures 4A,B). The composition of immune cells in TME varied significantly between both intragroup and intergroup. Additionally, the components of TME at single-cell resolution were investigated using TISCH database. TIICs of GSE146771 cohort and each patient’s cell type proportion were displayed (Figures 4C,D). GSEA showed the enriched upregulated hallmark pathways in different cells (Figure 4E). Interestingly, the pathways enriched in monocyte-macrophages, “Epithelial mesenchymal transition,” “Angiogenesis,” “Complement,” “Inflammatory response” and “IL6 JAK STAT3 signaling,” were also significantly enriched in high SCG2 expression group (Supplementary Figure S3).

FIGURE 4

FIGURE 4. The proportions of TIICs in TME and enrichment analysis by GSEA. (A,B) Composition of 22 immune cells in the tumor and normal tissues of CRC patients in TCGA cohort. (C,D) Components of TME at the single-cell resolution of CRC patients in GSE146771 cohort. (E) Enriched upregulated hallmark pathways in different immune cells.

SCG2 Expression Significantly Correlate With TIICs in CRC

We conducted ESTIMATE analysis to explore the relationship between SCG2 and TME. The results suggested that high SCG2 expression was significantly associated with higher immune score and stromal store in CRC (Figures 5A,B). CIBERSORT analysis showed that the proportions of resting memory CD4 cells, monocytes, activated dendritic cells, and activated mast cells were downregulated in the high SCG2 expression group, whereas M0 macrophages, M2 macrophages, and resting mast cells were significantly upregulated in the high SCG2 group (Figures 5C,D). We further performed correlation analysis between SCG2 expression and various TIICs gene markers (Table 2). Significant correlations were observed between SCG2 expression and T cell exhaustion (immune checkpoint genes), general T cells, CD8 + T cells, CD4 + T cells, Th1, Tfh, Th17, Monocyte, M1 macrophages, M2 macrophages, natural killer (NK) cell, dendritic cells, neutrophils in CRC. The results demonstrated that SCG2 expression was correlated with TME and TIICs in CRC.

FIGURE 5

FIGURE 5. Correlation between SCG2 expression and TIICs. (A) The correlation analysis between SCG2 expression and immune score, stromal score in TCGA cohort. (B) The correlation analysis between SCG2 expression and immune score, stromal score in GSE39582 cohort. (C,D) Twenty-two kinds of tumor-infiltrating immune cells are plotted according to the SCG2 expression level in TCGA and GSE39582 cohort. (A,B) Examined by Spearman’s correlation analysis; (C,D) Examined by Mann-Whitney U test; *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.

TABLE 2

TABLE 2. Correlation analysis between SCG2 and markers of immune cells.

SCG2 Expression Associated With M2 Macrophage Polarization

The expressions of TAM and M2 macrophage marker genes were higher in monocyte-macrophage cells from tumor samples than from normal samples or peripheral blood mononuclear cell (PBMC) (Supplementary Figure S4). It indicated that, in the TME, macrophages tend to polarize towards M2 phenotype and differentiate into TAMs thus performed pro-tumoral functions. CIBERSORT analysis suggested the proportion of M2 macrophages increased in the SCG2 high expression group, but the proportion of M1 macrophages not significantly changed (Figures 5C,D). Additionally, clear correlations existed between SCG2 expression and the marker genes of M2 macrophages and TAMs (Figures 6A–F). The aforementioned results prompted us to hypothesize that SCG2 could promote the polarization of M2 macrophages. To validate this conjecture, we performed confocal immunofluorescence microscopy to confirm the colocalization between SCG2 and macrophage. As shown in Figure 7, SCG2 was colocalized with macrophages in tumor tissues, whereas there was no apparent colocalization in normal tissues. Altogether, the above results suggested that SCG2 might promote tumor infiltrating macrophages polarization toward M2 phenotype and play a pro-tumoral role in CRC.

FIGURE 6

FIGURE 6. Correlation of SCG2 expression with macrophage polarization. (A–C) Correlation analysis between SCG2 expression and TAM markers genes (CCL2, CD68, IL10). (D–F) Correlation analysis between SCG2 expression and M2 Macrophage markers genes (CD163, MS4A4A, VSIG4). Examined by Spearman’s correlation analysis.

FIGURE 7

FIGURE 7. Immunofluorescence images of DAPI (blue), CD68 (green), SCG2 (red), and the merge photos in normal tissues and cancer tissues of CRC patients.

Discussion

CRC is a common health problem and one of the leading causes of cancer death worldwide (Ferlay et al., 2015). Despite numerous studies effort to improve our surgical treatment, radiotherapy, chemotherapy, and immunotherapy over the years, the prognosis of patients with advanced CRC remains poor (Seymour et al., 2007). In the present study, we explored the role of SCG2 in colorectal cancer, revealed its prognostic value, biological functions, associated pathways, and regulation of tumor immunity by analyzing open-access databases comprehensively.

SCG2 expression was remarkably decreased in tumor tissues compared with normal tissues of CRC. In normal tissues, our results showed that SCG2 was mainly distributed in colonic glands, which possibly related to intestinal fluid secretion. In contrast, in tumor tissues, the intestinal mucosal tissues were destroyed and total SCG2 expression was decreased. Additionally, SCG2 expression was associated with T, N, M, and TNM stages in CRC patients. Higher expression of SCG2 was significantly correlated to worse OS and DFS. Multivariate cox analysis also showed that SCG2 expression level was an independent prognostic marker in CRC.

To investigate the biological roles of SCG2, we carried out a gene differential expression analysis based on SCG2 median expression level. Then functional enrichment analyses were performed using the genes positively correlated with SCG2. The results suggested that SCG2 was strongly correlated with MYH11, SYNPO2, DDR2, FABP4, and TNS1, which were all involved in tumorigenesis and inflammatory immune response (Gan et al., 2019; Gao et al., 2020; Liu et al., 2020; Matrai et al., 2020; Sun et al., 2020; Tian et al., 2020). Higher expression of these genes might associate with worse survival and higher clinical stage of CRC patients.

Functional enrichment analysis revealed that SCG2 was correlated with several tumorigenesis related pathways such as “PI3K-Akt,” “Wnt,” and “TGF-beta pathways” (Schatoff et al., 2017; Koveitypour et al., 2019). A few pathways related to tumor metastasis included “ECM-receptor interaction” and “focal adhesion” (Gu et al., 2018; Machackova et al., 2020). Furthermore, numerous genes were involved in immune-related processes. “Regulation of actin cytoskeleton” plays a vital role in immune response, tumor migration, and invasion (Schmitz et al., 2000; C. Wickramarachchi et al., 2010); “B cell receptor” is initiated via interaction between B cell receptor and specific antigens and plays a crucial role in immune system (Li et al., 2019); “leukocyte transendothelial migration” tightly regulates inflammation and is associated with the transient crossing of leukocytes through the blood vessel wall (Schimmel et al., 2017). These results indicated that SCG2 might have regulatory roles in tumor progression and immune moderation.

The result of GSEA suggested high SCG2 expression might activate the IL6 JAK STAT3 signaling pathway and IL2–STAT5 signaling pathway in CRC patients. Previous studies have shown that STAT3 was activated in TIICs, including TAMs, amplifying immune suppression, and targeting IL-6 to regulate STAT3 signaling pathway was considered a potential immunotherapy for CRC (Wang and Sun, 2014; Verdeil et al., 2019). STAT5 also played a critical role in tumor immunity, regulated Treg cell’s function and development. Consistent activation of STAT5 was associated with suppressing antitumor immunity and increased tumor proliferation and invasion (Rani and Murphy, 2016). These findings indicated SCG2 might relate to the efficiency of immunotherapy in CRC patients.

The TIICs in TME are critical players in tumor progression, modulate tumor inflammation and metastasis variously (Gonzalez et al., 2018). We quantified each CRC sample’s immune and stromal cell infiltration levels using immune and stromal scores evaluated by ESTIMATE algorithm, respectively. Our results showed the two scores were significantly correlated with SCG2 expression in CRC patients. Some studies have demonstrated that both immune and stromal scores were associated with poor prognosis in CRC patients (Liu et al., 2020; Yuan et al., 2020).

CIBERSORT analysis showed that multiple immune cells had significantly different proportions between SCG2 high and low expression groups. Notably, the proportions of M0 and M2 macrophages in the SCG2 high expression group were considerably higher than those in the SCG2 low expression group, whereas M1 macrophages showed no significant change. Besides, SCG2 expression was moderately or strongly correlated with M2 macrophage marker genes and TAM marker genes but was slightly correlated with M1 macrophage marker genes. Similarly, confocal immunofluorescence microscopy showed SCG2 had an apparent colocalization with macrophages in tumor tissues. The previous study has demonstrated that SN, SCG2 derived peptide, could induce macrophage accumulation and angiogenesis (Liu et al., 2018). These results indicated that high expression of SCG2 might promote M0 macrophages polarize to M2 and eventually differentiate into TAMs, enhancing tumor cell invasion, metastasis, angiogenesis, and inhibiting the anti-tumoral immune surveillance (Kim and Bae, 2016), predict a worse prognosis of CRC patients (Waniczek et al., 2017).

We also found a positive correlation between SCG2 expression and immune checkpoint markers PD-1, PD-L1/2, CTLA-4, TIM-3, TIGIT in CRC patients. The binding of PD-L1/2 to the PD-1 receptor leads to T-cell function’s deactivation, allowing tumor cells to evade immune attacks (Chen and Han, 2015). Upregulation of PD-1 and Tim-3 was associated with the poor prognosis of CRC patients in stage I-III (Kuai et al., 2020), and high expression of TIGIT was associated with advanced TNM stage and poor DFS in CRC patients with mismatch repair deficiency (Zhou et al., 2020). The correlations between SCG2 expression and immune checkpoint genes indicated a potential mechanism of SCG2 regulation on T cell exhaustion and provided a new target for tumor immunotherapy.

However, our study had some limitations. First, we conducted analysis mainly based on samples from TCGA and GEO cohorts, it would be better to confirm the findings in a larger sample size. Second, genetic difference analysis based on RNA-seq might not be precise. Additionally, the mechanisms of SCG2 regulating the infiltration of immune cells and macrophage polarization require further experimental investigation, which would be a future direction for our research.

In conclusion, our study revealed that overexpression of SCG2 predicted poor prognosis and advanced clinical stage in CRC patients. SCG2 was associated with tumor immune cells infiltration, promoted M2 macrophage polarization, and correlated with immune checkpoint expression in CRC. In summary, SCG2 played a critical role in the regulation of tumor immunity and made it a potential biomarker and therapeutic target in CRC.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by the Zhongnan Hospital Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

HoW: Data acquisition, Methodology, Writing -original draft. JY: Software, Writing–review, and editing. YH: Resources. AR: Softwares. HzW: Visualization. ML: Softwares.

Funding

This work was supported by the Hubei Provincial National Natural Science Foundation of China (Project 2018CFB446) and the Health Commission of Hubei Province Scientific Research (Project WJ 2019H025).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We acknowledge the public database TIMER, TCGA, GEO, CPTAC, TISCH for free use.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.795133/full#supplementary-material

References

Albrecht-Schgoer, K., Schgoer, W., Holfeld, J., Theurl, M., Wiedemann, D., Steger, C., et al. (2012). The Angiogenic Factor Secretoneurin Induces Coronary Angiogenesis in a Model of Myocardial Infarction by Stimulation of Vascular Endothelial Growth Factor Signaling in Endothelial Cells. Circulation 126 (21), 2491–2501. doi:10.1161/circulationaha.111.076950

PubMed Abstract | CrossRef Full Text | Google Scholar

Beuret, N., Stettler, H., Renold, A., Rutishauser, J., and Spiess, M. (2004). Expression of Regulated Secretory Proteins Is Sufficient to Generate Granule-like Structures in Constitutively Secreting Cells. J. Biol. Chem. 279 (19), 20242–20249. doi:10.1074/jbc.m310613200

CrossRef Full Text | Google Scholar

Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA: A Cancer J. Clinicians 68 (6), 394–424. doi:10.3322/caac.21492

CrossRef Full Text | Google Scholar

Chen, L., and Han, X. (2015). Anti-PD-1/PD-L1 Therapy of Human Cancer: Past, Present, and Future. J. Clin. Invest. 125 (9), 3384–3391. doi:10.1172/jci80011

CrossRef Full Text | Google Scholar

Chung, K. Y., Gore, I., Fong, L., Venook, A., Beck, S. B., Dorazio, P., et al. (2010). Phase II Study of the Anti-cytotoxic T-Lymphocyte-Associated Antigen 4 Monoclonal Antibody, Tremelimumab, in Patients with Refractory Metastatic Colorectal Cancer. Jco 28 (21), 3485–3490. doi:10.1200/jco.2010.28.3994

CrossRef Full Text | Google Scholar

Cury, S. S., de Moraes, D., Freire, P. P., de Oliveira, G., Marques, D. V. P., Fernandez, G. J., et al. (2019). Tumor Transcriptome Reveals High Expression of IL-8 in Non-small Cell Lung Cancer Patients with Low Pectoralis Muscle Area and Reduced Survival. Cancers (Basel) 11 (9), 1251. doi:10.3390/cancers11091251

PubMed Abstract | CrossRef Full Text | Google Scholar

C. Wickramarachchi, D., Theofilopoulos, A. N., and Kono, D. H. (2010). Immune Pathology Associated with Altered Actin Cytoskeleton Regulation. Autoimmunity 43 (1), 64–75. doi:10.3109/08916930903374634

PubMed Abstract | CrossRef Full Text | Google Scholar

Deschoolmeester, V., Baay, M., Van Marck, E., Weyler, J., Vermeulen, P., Lardon, F., et al. (2010). Tumor Infiltrating Lymphocytes: an Intriguing Player in the Survival of Colorectal Cancer Patients. BMC Immunol. 11, 19. doi:10.1186/1471-2172-11-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferlay, J., Soerjomataram, I., Dikshit, R., Eser, S., Mathers, C., Rebelo, M., et al. (2015). Cancer Incidence and Mortality Worldwide: Sources, Methods and Major Patterns in GLOBOCAN 2012. Int. J. Cancer 136 (5), E359–E386. doi:10.1002/ijc.29210

CrossRef Full Text | Google Scholar

Gan, L., Camarena, V., Mustafi, S., and Wang, G. (2019). Vitamin C Inhibits Triple-Negative Breast Cancer Metastasis by Affecting the Expression of YAP1 and Synaptopodin 2. Nutrients 11 (12), 2997. doi:10.3390/nu11122997

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, J., Zhang, H.-P., Sun, Y.-H., Guo, W.-Z., Li, J., Tang, H.-W., et al. (2020). Synaptopodin-2 Promotes Hepatocellular Carcinoma Metastasis via Calcineurin-Induced Nuclear-Cytoplasmic Translocation. Cancer Lett. 482, 8–18. doi:10.1016/j.canlet.2020.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez, H., Hagerling, C., and Werb, Z. (2018). Roles of the Immune System in Cancer: from Tumor Initiation to Metastatic Progression. Genes Dev. 32 (19-20), 1267–1284. doi:10.1101/gad.314617.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, C., Wang, X., Long, T., Wang, X., Zhong, Y., Ma, Y., et al. (2018). FSTL1 Interacts with VIM and Promotes Colorectal Cancer Metastasis via Activating the Focal Adhesion Signalling Pathway. Cell Death Dis. 9 (6), 654. doi:10.1038/s41419-018-0695-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Guillemot, J., Thouënnon, E., Guérin, M., Vallet-Erdtmann, V., Ravni, A., Montéro-Hadjadje, M., et al. (2012). Differential Expression and Processing of Secretogranin II in Relation to the Status of Pheochromocytoma: Implications for the Production of the Tumoral Marker EM66. J. Mol. Endocrinol. 48 (2), 115–127. doi:10.1530/jme-11-0077

CrossRef Full Text | Google Scholar

Halama, N., Michel, S., Kloor, M., Zoernig, I., Benner, A., Spille, A., et al. (2011). Localization and Density of Immune Cells in the Invasive Margin of Human Colorectal Cancer Liver Metastases Are Prognostic for Response to Chemotherapy. Cancer Res. 71 (17), 5670–5677. doi:10.1158/0008-5472.can-11-0268

PubMed Abstract | CrossRef Full Text | Google Scholar

Hannon, P. R., Duffy, D. M., Rosewell, K. L., Brännström, M., Akin, J. W., and Curry, T. E. (2018). Ovulatory Induction of SCG2 in Human, Nonhuman Primate, and Rodent Granulosa Cells Stimulates Ovarian Angiogenesis. Endocrinology 159 (6), 2447–2458. doi:10.1210/en.2018-00020

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J., and Bae, J.-S. (2016). Tumor-Associated Macrophages and Neutrophils in Tumor Microenvironment. Mediators Inflamm. 2016, 6058147. doi:10.1155/2016/6058147

PubMed Abstract | CrossRef Full Text | Google Scholar

Koveitypour, Z., Panahi, F., Vakilian, M., Peymani, M., Seyed Forootan, F., Nasr Esfahani, M. H., et al. (2019). Signaling Pathways Involved in Colorectal Cancer Progression. Cell Biosci. 9, 97. doi:10.1186/s13578-019-0361-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuai, W., Xu, X., Yan, J., Zhao, W., Li, Y., Wang, B., et al. (2020). Prognostic Impact of PD-1 and Tim-3 Expression in Tumor Tissue in Stage I-III Colorectal Cancer. Biomed. Res. Int. 2020, 5294043. doi:10.1155/2020/5294043

PubMed Abstract | CrossRef Full Text | Google Scholar

Le, D. T., Uram, J. N., Wang, H., Bartlett, B. R., Kemberling, H., Eyring, A. D., et al. (2015). PD-1 Blockade in Tumors with Mismatch-Repair Deficiency. N. Engl. J. Med. 372 (26), 2509–2520. doi:10.1056/NEJMoa1500596

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Yin, W., Jing, Y., Kang, D., Yang, L., Cheng, J., et al. (2019). The Coordination between B Cell Receptor Signaling and the Actin Cytoskeleton during B Cell Activation. Front. Immunol. 9, 3096. doi:10.3389/fimmu.2018.03096

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.can-17-0307

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J. W., Yu, F., Tan, Y.-F., Huo, J.-P., Liu, Z., Wang, X.-J., et al. (2020). Profiling of Tumor Microenvironment Components Identifies Five Stroma-Related Genes with Prognostic Implications in Colorectal Cancer. Cancer Biother. Radiopharm. doi:10.1089/cbr.2020.4118

CrossRef Full Text | Google Scholar

Liu, W., Wang, F., Zhao, M., Fan, Y., Cai, W., and Luo, M. (2018). The Neuropeptide Secretoneurin Exerts a Direct Effect on Arteriogenesis In Vivo and In Vitro. Anat. Rec. 301 (11), 1917–1927. doi:10.1002/ar.23929

PubMed Abstract | CrossRef Full Text | Google Scholar

Luebeck, E. G., Curtius, K., Jeon, J., and Hazelton, W. D. (2013). Impact of Tumor Progression on Cancer Incidence Curves. Cancer Res. 73 (3), 1086–1096. doi:10.1158/0008-5472.can-12-2198

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, Y., Chen, L., Zhou, Q., Xiong, Y., Wang, G., Liu, X., et al. (2020). Identification of a Prognostic Gene Signature Based on an Immunogenomic Landscape Analysis of Bladder Cancer. J. Cel Mol Med. 24 (22), 13370–13382. doi:10.1111/jcmm.15960

CrossRef Full Text | Google Scholar

Machackova, T., Vychytilova-Faltejskova, P., Souckova, K., Trachtova, K., Brchnelova, D., Svoboda, M., et al. (2020). MiR-215-5p Reduces Liver Metastasis in an Experimental Model of Colorectal Cancer through Regulation of ECM-Receptor Interactions and Focal Adhesion. Cancers (Basel) 12 (12), 3518. doi:10.3390/cancers12123518

PubMed Abstract | CrossRef Full Text | Google Scholar

Matrai, C., Motanagh, S., Mirabelli, S., Ma, L., He, B., Chapman-Davis, E., et al. (2020). Molecular Profiles of Mixed Endometrial Carcinoma. Am. J. Surg. Pathol. 44 (8), 1104–1111. doi:10.1097/PAS.0000000000001519

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyerson, M., Gabriel, S., and Getz, G. (2010). Advances in Understanding Cancer Genomes through Second-Generation Sequencing. Nat. Rev. Genet. 11 (10), 685–696. doi:10.1038/nrg2841

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, K. D., Siegel, R. L., Lin, C. C., Mariotto, A. B., Kramer, J. L., Rowland, J. H., et al. (20162016). Cancer Treatment and Survivorship Statistics, 2016. CA: A Cancer J. Clinicians 66 (4), 271–289. doi:10.3322/caac.21349

CrossRef Full Text | Google Scholar

Newman, A. M., Liu, C. L., Green, M. R., Gentles, A. J., Feng, W., Xu, Y., et al. (2015). Robust Enumeration of Cell Subsets from Tissue Expression Profiles. Nat. Methods 12 (5), 453–457. doi:10.1038/nmeth.3337

PubMed Abstract | CrossRef Full Text | Google Scholar

Overman, M. J., Lonardi, S., Wong, K. Y. M., Lenz, H.-J., Gelsomino, F., Aglietta, M., et al. (2018). Durable Clinical Benefit with Nivolumab Plus Ipilimumab in DNA Mismatch Repair-Deficient/Microsatellite Instability-High Metastatic Colorectal Cancer. Jco 36 (8), 773–779. doi:10.1200/jco.2017.76.9901

PubMed Abstract | CrossRef Full Text | Google Scholar

Rani, A., and Murphy, J. J. (2016). STAT5 in Cancer and Immunity. J. Interferon Cytokine Res. 36 (4), 226–237. doi:10.1089/jir.2015.0054

PubMed Abstract | CrossRef Full Text | Google Scholar

Rotte, A. (2019). Combination of CTLA-4 and PD-1 Blockers for Treatment of Cancer. J. Exp. Clin. Cancer Res. 38 (1), 255. doi:10.1186/s13046-019-1259-z

CrossRef Full Text | Google Scholar

Schatoff, E. M., Leach, B. I., and Dow, L. E. (2017). Wnt Signaling and Colorectal Cancer. Curr. Colorectal Cancer Rep. 13 (2), 101–110. doi:10.1007/s11888-017-0354-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Schimmel, L., Heemskerk, N., and van Buul, J. D. (2017). Leukocyte Transendothelial Migration: A Local Affair. Small GTPases 8 (1), 1–15. doi:10.1080/21541248.2016.1197872

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmitz, A. A. P., Govek, E.-E., Böttner, B., and Van Aelst, L. (2000). Rho GTPases: Signaling, Migration, and Invasion. Exp. Cel Res. 261 (1), 1–12. doi:10.1006/excr.2000.5049

CrossRef Full Text | Google Scholar

Seymour, M. T., Maughan, T. S., Ledermann, J. A., Topham, C., James, R., Gwyther, S. J., et al. (2007). Different Strategies of Sequential and Combination Chemotherapy for Patients with Poor Prognosis Advanced Colorectal Cancer (MRC FOCUS): a Randomised Controlled Trial. The Lancet 370 (9582), 143–152. doi:10.1016/s0140-6736(07)61087-3

CrossRef Full Text | Google Scholar

Siegel, R. L., Miller, K. D., Goding Sauer, A., Fedewa, S. A., Butterly, L. F., Anderson, J. C., et al. (2020). Colorectal Cancer Statistics, 2020. CA A. Cancer J. Clin. 70 (3), 145–164. doi:10.3322/caac.21601

CrossRef Full Text | Google Scholar

Stuart, J. M., Segal, E., Koller, D., and Kim, S. K. (2003). A Gene-Coexpression Network for Global Discovery of Conserved Genetic Modules. Science 302 (5643), 249–255. doi:10.1126/science.1087447

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, D., Wang, J., Han, Y., Dong, X., Ge, J., Zheng, R., et al. (2021). TISCH: a Comprehensive Web Resource Enabling Interactive Single-Cell Transcriptome Visualization of Tumor Microenvironment. Nucleic Acids Res. 49 (D1), D1420–D1430. doi:10.1093/nar/gkaa1020

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y. L., Zhang, Y., Guo, Y. C., Yang, Z. H., and Xu, Y. C. (2020). A Prognostic Model Based on Six Metabolism-Related Genes in Colorectal Cancer. Biomed. Res. Int. 2020, 5974350. doi:10.1155/2020/5974350

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, W., Zhang, W., Zhang, Y., Zhu, T., Hua, Y., Li, H., et al. (2020). FABP4 Promotes Invasion and Metastasis of colon Cancer by Regulating Fatty Acid Transport. Cancer Cel Int. 20, 512. doi:10.1186/s12935-020-01582-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Troger, J., Theurl, M., Kirchmair, R., Pasqua, T., Tota, B., Angelone, T., et al. (2017). Granin-derived Peptides. Prog. Neurobiol. 154, 37–61. doi:10.1016/j.pneurobio.2017.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Verdeil, G., Lawrence, T., Schmitt-Verhulst, A. M., and Auphan-Anezin, N. (2019). Targeting STAT3 and STAT5 in Tumor-Associated Immune Cells to Improve Immunotherapy. Cancers (Basel) 11 (12), 1832. doi:10.3390/cancers11121832

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S.-W., and Sun, Y.-M. (2014). The IL-6/JAK/STAT3 Pathway: Potential Therapeutic Strategies in Treating Colorectal Cancer. Int. J. Oncol. 44 (4), 1032–1040. doi:10.3892/ijo.2014.2259

CrossRef Full Text | Google Scholar

Wang, Z., and Chen, X. (2020). Establishment and Validation of an Immune-Associated Signature in Lung Adenocarcinoma. Int. Immunopharmacology 88, 106867. doi:10.1016/j.intimp.2020.106867

PubMed Abstract | CrossRef Full Text | Google Scholar

Waniczek, D., Lorenc, Z., Śnietura, M., Wesecki, M., Kopec, A., and Muc-Wierzgoń, M. (2017). Tumor-Associated Macrophages and Regulatory T Cells Infiltration and the Clinical Outcome in Colorectal Cancer. Arch. Immunol. Ther. Exp. 65 (5), 445–454. doi:10.1007/s00005-017-0463-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Whiteaker, J. R., au, fnm., Halusa, G. N., Hoofnagle, A. N., Sharma, V., MacLean, B., et al. (2014). CPTAC Assay Portal: a Repository of Targeted Proteomic Assays. Nat. Methods 11 (7), 703–704. doi:10.1038/nmeth.3002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Han, L., Yuan, Y., Li, J., Hei, N., and Liang, H. (2014). Gene Co-expression Network Analysis Reveals Common System-Level Properties of Prognostic Genes across Cancer Types. Nat. Commun. 5, 3231. doi:10.1038/ncomms4231

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshihara, K., Shahmoradgoli, M., Martínez, E., Vegesna, R., Kim, H., Torres-Garcia, W., et al. (2013). Inferring Tumour Purity and Stromal and Immune Cell Admixture from Expression Data. Nat. Commun. 4, 2612. doi:10.1038/ncomms3612

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, W., Cai, W., Huang, X., and Peng, S. (2020). Prognostic Value of Immune Scores in the Microenvironment of Colorectal Cancer. Oncol. Lett. 20 (5), 1. doi:10.3892/ol.2020.12119

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Ding, X., Li, H., Yang, C., Ma, Z., Xu, G., et al. (2020). Upregulation of TIGIT and PD-1 in Colorectal Cancer with Mismatch-Repair Deficiency. Immunological Invest. 50, 338–355. doi:10.1080/08820139.2020.1758130

PubMed Abstract | CrossRef Full Text | Google Scholar