- Open Access
Follistatin-like 1 (FSTL1) is a prognostic biomarker and correlated with immune cell infiltration in gastric cancer
World Journal of Surgical Oncology volume 18, Article number: 324 (2020)
Follistatin-like 1 (FSTL1) plays a central role in the progression of tumor and tumor immunity. However, the effect of FSTL1 on the prognosis and immune infiltration of gastric cancer (GC) remains to be elucidated.
The expression of FSTL1 data was analyzed in Oncomine and TIMER databases. Analyses of clinical parameters and survival data were conducted by Kaplan-Meier plotter and immunohistochemistry. Western blot assay and real-time quantitative PCR (RT-qPCR) were used to analyze protein and mRNA expression, respectively. The correlations between FSTL1 and cancer immune infiltrates were analyzed by Tumor Immune Estimation Resource (TIME), Gene Expression Profiling Interactive Analysis (GEPIA), and LinkedOmics database.
The expression of FSTL1 was significantly higher in GC tissues than in normal tissues, and bioinformatic analysis and immunohistochemistry (IHC) indicated that high FSTL1 expression significantly correlated with poor prognosis in GC. Moreover, FSTL1 was predicted as an independent prognostic factor in GC patients. Bioinformatics analysis results suggested that FSTL1 mainly involved in tumor progression and tumor immunity. And significant correlations were found between FSTL1 expression and immune cell infiltration in GC.
The study effectively revealed useful information about FSTL1 expression, prognostic values, potential functional networks, and impact of tumor immune infiltration in GC. In summary, FSTL1 can be used as a biomarker for prognosis and evaluating immune cell infiltration in GC.
Gastric cancer (GC) is the fourth most common cancer worldwide; despite the rapid development of diagnosis and therapies, GC is still one of the most leading cause of death . Thus, it is urgent to investigate the mechanism of cancer progression, as well as explore the potential biomarkers of diagnosis, prognosis, and therapy in GC. In recent years, immunotherapy has rapidly developed into a promising and effective therapeutic strategy for GC patients ; however, the molecular mechanisms related to immune dysfunction remain unknown. Several studies have uncovered significant correlations between tumor microenvironment (TME) and the mechanisms of immune dysfunction [3,4,5]. Tumor-infiltrating lymphocytes (TILs), including tumor-associated macrophages (TAMs) and tumor-infiltrating neutrophils (TINs), display unique histological features that can distinguish the differentiation of malignancies and often vary between individual patients. In addition, TILs have been proven to play essential roles in the prognosis and treatment of malignancies [6,7,8,9]. Hence, it is important to clarify the immunophenotype of tumor-immune interactions, as this will lead to the identification of novel immunotherapeutic-related targets in various malignancies.
Follistatin-like 1 (FSTL1) is a transmembrane extracellular glycoprotein of the secreted protein acid and belongs to the BM-40/SPARC/osteonectin family. FSTL1 was first cloned in 1988 from the synovial tissues of patients with rheumatoid arthritis after it was detected in the serum and synovial fluid of patients with the disease . Immunological studies suggest that FSTL1 regulates joint inflammation in arthritic diseases, while also modulating allograft tolerance. In addition, FSTL1 acts as an autoantigen in patients with rheumatoid arthritis by inhibiting IL6 and producing IL17 to play cardioprotective roles in several cardiovascular diseases [10, 11]. In previous reports, FSTL1 has been shown to produce varying effects in terms of tumors progression and inhibition. For example, FSTL1 has been found to promote tumor inhibition in breast cancer  and clear-cell renal cell carcinoma ; however, it is a poor prognostic indicator of colorectal cancer , esophageal squamous cell carcinoma , and hepatocellular carcinoma . FSTL1 also affects the progression of multiple tumors by modulating the host immune system, leading to immune dysfunction in many cases. FSTL1-induced immune dysfunction leads to tumor bone metastases  and activation of immune evasion mechanisms in patients with non-small cell lung cancer (NSCLC) through the FSTL1-DIP2A axis . FSTL1 plays a crucial role in tumor-immune and cancer cell progression. However, the effects of FSTL1 on GC prognosis and immune infiltration are still unclear. The potential functions and mechanisms of FSTL1 in cancer progression and immunology need to be discovered.
In this study, we investigated the expression of FSTL1 in GC using bioinformatics analysis, as well as clinical samples, and analyzed the association of FSTL1 expression with different TILs in TME using the Tumor Immune Estimation Resource (TIMER) database and CIBERSORT method. The results demonstrate that FSTL1 holds prognostic value, and explain the potential mechanism and impact of FSTL1 and TME in GC.
FSTL1 mRNA expression was assessed in GC using Oncomine V.4.5 (www.oncomine.org), Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/index.html), and TIMER database (https://cistrome.shinyapps.io/timer/) [19,20,21]. The correlation between FSTL1 expression and survival in GC were analyzed with the Kaplan-Meier plotter . TIMER database was used to investigate FSTL1 expression in multiple cancers and the correction between FSTL1 expression and TME in GC [20, 23]. The GSVA package of R language was used to detect activity changes of immune-related pathway after FSTL1 expression changes , and the sources of immune-related pathway were used in Molecular Signatures Database (http://software.broadinstitute.org/gsea/msigdb) . Median FSTL1 expression level of tumors from The Cancer Genome Atlas (TCGA) database (https://cancergenome.nih.gov) was grouped into high and low expression groups. Immune-related pathway predictions in different levels of FSTL1 expression were performed using immune-related signatures. The screening criteria were p value < 0.01 and false discovery rate (FDR) < 0.01. For the prediction of tumor-infiltrating lymphocytes (TILs), the expression matrices were uploaded to CIBERSORT (https://cibersort.stanford.edu/) . The valid samples were selected using p value < 0.05. Both gene ontology (GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed utilizing the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6. 8[27, 28].
Patients and tissue specimens
All the clinical GC samples were obtained from the First Affiliated Hospital of Nanchang University. Formalin-fixed, paraffin-embedded samples from 240 GC patients were collected from January 2011 to October 2012, and 50 fresh GC tissues and paired adjacent noncancerous tissues were stored in liquid nitrogen before use. Before the operation, the patients had not accepted chemotherapy, radiotherapy, or other treatments. All samples were collected with the consent of the patient, and the study was approved by the Ethics Committee of the First Affiliated Hospital of Nanchang University. All patient specimens and clinical data involved in this study complied with the Declaration of Helsinki.
RNA extraction and real-time quantitative PCR
Total RNA from GC fresh tissues was extracted by TRIzol reagent (Invitrogen, USA) and reverse transcribed by EasyScript First-Strand cDNA Synthesis SuperMix Kit (TransGen Biotech, China). RT-qPCR was evaluated by the StepOnePlus Real-Time PCR System (Applied Biosystems, USA) and Fast Start Universal SYBR Green Master Mix (Takara, Japan). The relative expression of FSTL1, forward (5′-GCCATGACCTGTGACGGAAA-3′); reverse (5′-CAGCGCTGAAGTGGAGAAGA-3′) relative to the internal control (GAPDH) was analyzed using the 2−ΔΔCT method.
The total protein was extracted with RIPA buffer containing protease inhibitor cocktail, separated with 10% SDS-PAGE gel and transferred to polyvinylidene difluoride membrane (Bio-Rad Laboratories, USA). And the membranes were blocked by 10% milk solution for 1 h at room temperature following incubated with primary antibodies against FSTL1 (1:1000, DF12274, Affinity, China) and β-actin (1:1000; #4970; Cell Signaling Technology, USA) for overnight at 4 °C. Then the membranes were washed 3 times in Tri-buffered saline with Tween-20, and incubated with horseradish peroxidase (HRP)-conjugated secondary antibodies (Anti-rabbit IgG, #7074, Cell Signaling Technology, USA) for 1 h at room temperature. After washed 3 times in Tri-buffered saline with Tween-20, the signals were detected with super sensitive regent (Thermo Fisher Scientific, USA).
Formalin-fixed paraffin-embedded tissues sections (4-μm) were degreased in xylene and rehydrated in different concentrations of alcohol and distilled water before antigen retrieval. The sections were blocked in 3% H2O2 at room temperature for 5–10 min to eliminate endogenous peroxidase activity. Then sections were incubated with anti-FSTL1 antibody (1:100, DF12274, Affinity, China) at 4 °C overnight. The slides were washed 2 times in PBS for 5 min each time and incubated with secondary antibodies for 40 min at room temperature, then washed 3 times in PBS. 3,3′-Diaminobenzidine and hematoxylin were used to color these slides. FSTL1 expression level was scored according to the proportion of stained tumor cells and the intensity of the staining . The staining intensity was scored as: 0 (negative), 1 (weak), 2 (moderate), and 3 (strong). The proportion of staining was scored as: 0 (0–5%); 1 (5–25%); 2 (25–50%); 3 (50–75%); and 4 (75–100%). The final FSTL1 expression scores were calculated by multiplying the above two scores. Slides were considered as low or high expression, with final scores of ≤ 6 or > 6, respectively. Two independent pathologists who were blinded observed the results under an optical microscope.
The survival curves of Kaplan-Meier plots and GEPIA databases were displayed by HR with p values or Cox p values from the log-rank test. Gene expression distribution from the TIMER database was displayed using box plots, with statistical significance assessed by the Wilcoxon test, and the correlation of gene expression from the TIMER databases was evaluated with Spearman’s correlation coefficient and estimated using statistical significance. Kaplan-Meier method and log-rank test were used to analyze overall survival curves. Differences in the mRNA expression levels of FSTL1 in fresh GC and matched normal tissues were analyzed using paired t tests. Univariate and multivariate Cox’s proportional hazard models were conducted to analyze the value of each clinic parameter for predicting overall survival. Statistical significance was defined as p < 0.05.
FSTL1 expression and potential prognosis in GC
TIMER and Oncomine databases were used to assess the differential expression of FSTL1 based on cancer type. The majority of data sets from the TIMER and Oncomine database showed that FSTL1 was highly expressed in most malignancies, especially in GC (Fig. 1a, b). In representative Chen Gastric data sets (N = 35) from the Oncomine database , FSTL1 expression was significantly high expression in GC tissues compared with normal tissues, and was elevated in multiple pathological typing, including adenocarcinoma, intestinal carcinoma, diffuse carcinoma, and mixed carcinoma (Fig. 1c, d). The same result was confirmed in GEPIA database, and the expression of FSTL1 increases with the stage of GC (Fig. 1e, f). In addition, we further confirmed that FSTL1 was highly expressed in GC through real-time quantitative PCR (RT-qPCR) and Western blot (Fig. 2a, f). Moreover, the mRNA expression of FSTL1 in the Nx group (P = 0.037) and stage III group (P = 0.018) was significantly higher than N0 group and the I/II group in GC patients, respectively (Fig. 2c, d). But there was no significant change in FSTL1 expression in stage T grouping and differentiation grouping (Fig. 2b, e).
Meanwhile, KM-plotter and GEPIA database survival analysis indicated that high FSTL1 expression indicated poor prognosis in GC patients (Supplementary Figure 1A, B). To confirm the predicted prognostic value of FSTL1 in GC, FSTL1 were analyzed in random selected 240 primary GC tissues and paired adjacent non-tumor tissues by immunohistochemistry. The staining of FSTL1 protein ranged from weak to strong (Fig. 2h), and the result showed that FSTL1 positive staining was increased in GC tissues than adjacent non-tumor tissues (Fig. 2i). The positive expression rate of FSTL1 (43.3%, 104/240) in GC samples was significantly higher than that in adjacent non-tumor tissues (15%, 36/240) (P < 0.001) (Fig. 2g).
Kaplan-Meier analysis and log-rank test were performed to analyze the prognostic value of FSTL1 expression in GC patients. The result showed that patients with positive FSTL1 expression (41.1 ± 3.8 month) had worse overall survival than negative FSTL1 expression (70.9 ± 3.1 month) (P < 0.001) (Fig. 2j).
Correlation of FSTL1 expression and clinical prognosis in GC with different clinicopathological features
To further understand the effect of FSTL1 in the prognosis of GC, we studied the relationship between FSTL1 expression and the clinicopathological features of these cancers using the Kaplan-Meier plotter database and immunohistochemistry assay.
In Kaplan-Meier plotter database, the expression of FSTL1 has a significant correlation with the prognosis of GC patients with various clinicopathological characteristics. It is worth noting that increased FSTL1 expression was associated with worse OS and PFS in stage III–IV and stage N1–3, but not in stage I/II and stage N0 GC patients (Table 1). The result means that FSTL1 expression level can impact the prognosis in the advanced and lymph node metastasis GC patients, which is consistent with the increased expression of FSTL1 in III group and Nx group compared with I/II group and N0 group in Fig. 2c, d.
Next, correlation between FSTL1 expression and clinicopathological was investigated (Table 2). The result showed that expression level of FSTL1 was correlated with tumor size (P < 0.001), lymph node metastasis (P < 0.001), and tumor-node-metastasis (TNM) stage (P = 0.001). No significant correlations between FSTL1 expression and ages (P = 0.48), gender (P = 0.519), depth of invasion (P = 0.147), differentiation (P = 0.151), nerve invasion (P = 0.809), and vascular invasion (P = 0.317) were detected. The results are consistent with the results in the KM-Plotter database.
Univariate and multivariate analyses were conducted to investigate the independent prognostic factors of GC patients (Table 3). Univariate analysis indicated that tumor size (P < 0.001), depth of invasion (P < 0.001), lymph node metastasis (P < 0.001), TNM stage (P < 0.001), vascular invasion (P = 0.003), and FSTL1 expression (P < 0.001) were significantly correlated with OS of GC patients. These factors were subjected to multivariate analysis, which indicated that tumor size (P = 0.006), depth of invasion (P = 0.007), lymph node metastasis (P < 0.001), TNM stage (P < 0.001), and FSTL1 expression (P = 0.002) were independent prognostic factors in GC patients.
Protein-protein interaction of FSTL1 network analysis
Genes co-expressed with FSTL1 were analyzed by LinkedOmics database, which analyzed mRNA sequencing data of 415 GC patients from TCGA. All the genes associated with FSTL1 GC were showed in Fig. 3a. And Fig. 3b revealed the most 50 significant genes positively correlated with FSTL1 in GC. The 50 genes differentially expressed associated with highly FSTL1 expression were primarily enriched in extracellular region/exosome/space/matrix/matrix organization, cell adhesion (Fig. 3c). KEGG pathway analysis found enrichment in pathway in cancer, focal adhesion, PI3K-Akt signaling pathway, and ECM-receptor interaction, which indicated that elevation of FSTL1 was significantly correlated with tumor development and immune response (Fig. 3d).
In addition, we applied GSVA and clustering on a set of immune-specific signatures to further validate the enrichment of FSTL1 expression in the immune-related pathway. As FSTL1 expression increased, most immune-related pathways were activated in GC, CYTOKINE_CYTOKINE_RECEPTOR_INTERACTION, INNATE_IMMUNE_SYSTEM, NATURAL_KILLER_CELL_MEDIATED_CYTOTOXICITY, and others in stomach adenocarcinoma (STAD) (Fig. 4). Hence, FSTL1 plays a crucial role in immune cell infiltration and tumor-immune system interactions in GC.
Correlation between FSTL1 and tumor environment in GC
TILs, which are essential components of the TME, have been shown to reflect the host antitumor immune response and are prognostic factors for cancer [31, 32]. As mentioned, GO analysis revealed that FSTL1 has functions associated with extracellular region/exosome/space/matrix/matrix organization, cell adhesion, while KEGG analysis indicated that FSTL1 involved in the ECM-receptor interaction pathway, PI3K-Akt signaling pathway. It is well-known that remodeling the extracellular matrix (ECM), disordered PI3K-Akt signaling pathway can exert a great influence on the TME [33, 34]. Thus, we analyzed the correlation between FSTL1 expression and TILs in GC using the TIMER database. The result showed that FSTL1 expression was found to be correlated with high immune infiltration of CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells in GC (Fig. 5).
Immune cells are frequently identified based on the expression of cell surface receptors and intracellular markers. To further assess the relationship between FSTL1 and each tumor-infiltrating immune cell, we analyzed the correlation between FSTL1 expression and immune infiltrating cells, including CD8+ T cells, T cells (general), B cells, neutrophils, NK cells and dendritic cells (DCs), monocytes, TAMs, and M1/M2 macrophages, in GC using the TIMER and GEPIA databases. In addition, the functional T cells were also investigated, including T helper type 1 (Th1), Th2, Th17, follicular helper T cell (Tfh), and regulatory T (Treg) cells, as well as exhausted T cells. Even if the correlation was adjusted based on tumor purity, the majority of tumor-infiltrating immune cell markers were positively correlated with FSTL1 expression in GC (Table 4).
Correlation between FSTL1 expression and macrophage polarization in GC
Studies have shown that polarization from the antitumor M1 (classically activated) macrophage to protumor M2 (alternatively activated) macrophage phenotype is correlated with tumor development [35, 36]. Interesting, we found that the markers for monocytes, TAMs, and M2 macrophages showed moderate-to-strong correlations with FSTL1 expression, while the M1 macrophage markers were weakly correlated with FSTL1 expression (Fig. 6a–d, Table 4). To confirm the results with the TIMER database, the correlation between FSTL1 expression and the cell markers for monocytes, TAMs, M1 macrophages, and M2 macrophages in GC was assessed using GEPIA database. Findings from the GEPIA database were consistent with the TIMER database (Supplementary Table 1). Therefore, we used the CIBERSORT method to further explore whether FSTL1 expression was correlated with the polarization of macrophages in GC. The results showed that high FSTL1 expression had a higher ratio of M2 macrophage (Fig. 6e–g).
In this study, we found that (1) FSTL1 was highly expressed and correlated with worse prognosis in GC patients, (2) high expression of FSTL1 was positively correlated with immune infiltration in GC, and (3) high expression of FSTL1 promotes the polarization of M1 macrophages to M2 macrophages in GC.
Studies have showed that FSTL1 was involved in tumor development and tumor immunity by a variety of mechanisms, such as promoting metastasis and chemoresistance through NFκB-BMP signaling pathway in esophageal squamous cell carcinoma , activating the focal adhesion signaling pathway to promotes colorectal cancer , and causing immune dysfunction in lung cancer  and tumor with bone metastasis . From these results, we can speculate that FSTL1 is crucially involved in tumor progression and tumor immunity.
In the present study, we found that FSTL1 expression is upregulated and indicated worse prognosis in GC using bioinformatics. And the results were confirmed by RT-qPCR, Western blot, and immunohistochemistry (IHC) assay. Moreover, the results indicated that high expression level of FSTL1 was correlated with tumor size, lymph node metastasis, and TNM stage, especially. The univariate and multivariate analysis of GC patients suggested that FSTL1 was an independent factor of prognosis in GC. All of the results suggested that FSTL1 may become a potential prognosis predictor in GC patients. To investigate the underlying mechanism, GO and KEGG analyses were performed in DAVID website. Co-expression genes of FSTL1 were enriched in extracellular region/exosome/space/matrix/matrix organization, cell adhesion, and pathway in cancer, focal adhesion, PI3K-Akt signaling pathway, ECM-receptor interaction, which have been proved correlated with tumor development and tumor immunity. Thus, we further investigated the relationship between FSTL1 expression and tumor immunity.
TME, which constitutes an important part of tumor immunity, is a complex environment containing a mixture of nutrients, chemokines, immune cells, and other critical non-cancerous components . The interaction between tumor cells and TME plays a crucial role in tumor growth, progression, chemotherapy, and radiotherapy resistance [39,40,41]. Importantly, tumor-infiltrating immune cells have been found to play roles in tumor growth, invasion, and metastasis in various forms of cancer . Another important aspect of this research is the close connection between FSTL1 expression and tumor-infiltrating immune cells in GC. The result from the GSVA analysis and TIMER database suggests that FSTL1 expression is closely related to immune-related pathways and significantly correlated with CD8+ T cell, CD4+ T cell, macrophages, neutrophils, and DCs in GC (Figs. 4 and 5). Hence, high FSTL1 expression can cause a substantial increase in immune cell tumor infiltration.
Existing research has shown that TAMs, which differentiated from peripheral monocytes, are infiltrating around tumor cells. TAMs secrete a variety of cytokines, which play an important role in tumor occurrence, metastasis, and invasion, under the stimulation of TME. According to the different activation status and function, macrophages can be divided into M1 (classically activated macrophages) macrophage and M2 (alternatively activated macrophages) macrophage. M1 macrophages secrete pro-inflammatory cytokines and chemokines, and present antigens on a full-time basis, participate in the positive immune response, and perform the function of immune surveillance; M2 macrophages have only weak antigen presentation ability and secrete inhibitory cytokines such as IL-10 or TGF-β downregulate the immune response and play an important role in immune regulation [43, 44]. It is well known that polarization from the antitumor M1 (classically activated) macrophage to protumor M2 (alternatively activated) macrophage phenotype is correlated with tumor development, and the proportion of M1 and M2 macrophages in TAM is closely related to the prognosis of tumor. The large number differentiation of M2 macrophage is correlated with worse prognosis, suggesting that the differentiation of macrophages has important guiding significance for the clinical treatment of tumor [36, 45, 46]. To further investigate the relationship between FSTL1 and tumor-infiltrating cells, we analyzed the connection between FSTL1 expression and tumor-infiltrating cell markers. First, we found that FSTL1 expression is significantly correlated with surface markers of monocytes (CD86, CD115), TAMs (CCL2, CD68, IL10), M1 macrophages (INOS, IRF5, COX2), and M2 macrophages (CD163, VSIG4, MS4A4A) (Fig. 6a–d, Table 4, Supplementary Table 1). This indicates that FSTL1 may play an important role in macrophage polarization and affect the progression of GC through macrophage polarization.
In terms of TILs, helper T lymphocytes, including Th1 cells, Th2 cells, Tfh cells, Th17 cells, and Treg cells, are considered to be the major players in tumor immunity [47,48,49]. The surface marker of helper T lymphocytes (Th1 cells, Th2 cells, Tfh cells, Th17 cells, and Treg cells) shows significant correlations with FSTL1 expression (Table 4). Hence, FSTL1 plays a vital role in regulating the potential mechanisms of T cell function in GC. A previous study suggested that high T-bet expression was previously correlated with better DFS and OS in GC patients . Interestingly, FSTL1 expression is positively correlated with T-bet in GC. This indicates that FSTL1 may affect the prognosis of GC by impacting the expression of T-bet. Treg cells have been reported to stimulate GC progression by secreting TGF-β to help the cancer cells escape immune system recognition, and high levels of TGF-β in the TME further stimulate the differentiation and expansion of Treg cells [51, 52]. Previous studies have suggested that Treg cells may be involved in the development and progression of bladder carcinomas and GC [53, 54]. Our results agree that FSTL1 expression is a positive biomarker for Treg cells, including TGF-β in GC. Hence, Treg cells play an essential role in the survival of patients with GC affected by FSTL1.
Furthermore, our results showed that FSTL1 expression is positively related to the infiltration level of DCs and the markers of DCs in GC (Fig. 5, Table 4), which indicated that FSTL1 has strong correlation with DCs infiltration in GC. However, the effects of tumor-infiltrating DCs in GC remain unclear. Additional studies are needed to determine how FSTL1 regulates DCs in GC.
Significant correlations have also been detected between FSTL1 and T cell exhaustion markers in GC (Table 4). Previous studies reported that Treg cells could reversibly inhibit the cytotoxic T cell function of attacking tumor cells . On the other hand, exhausted T cells co-express a wide range of surface inhibitory molecules, such as PD-1, CTLA4, LAG3, and TIM-3 . In general, the higher the number of inhibitory receptors co-expressed by exhausted T cells, the more severe the T cells exhaustion. Although the expression of a single inhibitory receptor does not indicate T cells exhaustion, the co-expression of several inhibitory receptors is an important feature [57, 58]. From our results, the inhibitory receptors, including PD-1, CTLA4, LAG3, and TIM-3, increased in GC with FSTL1 overexpression (Table 4), suggesting that FSTL1 plays a crucial role in the regulation of T cell exhaustion. Hence, FSTL1 is closely correlated with tumor-infiltrating immune cells and indicates that FSTL1 plays a vital role in the immune escape mechanism in patients with GC. Recently, FSTL1 neutralizing antibodies were developed and found to display excellent efficacy in vitro . Also, CTLA4, PD-1, and programmed death ligand-1 (PD-L1) inhibitors have shown promising results in lung cancer, yet the response was unclear in advanced-stage GC [60, 61]. The combined use of FSTL1 neutralizing antibodies with anti-CTLA-4 or anti-PD-1/PD-L1 inhibitors may produce additive or synergistic effects, thus creating a new treatment option. Of course, part of these results have not been proved or confirmed in our specimens, and further investigation should be conducted in the future.
In summary, FSTL1 is commonly expressed and possesses significant prognostic value in GC. High FSTL1 expression is associated with high infiltration levels of CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and DCs in GC, and potentially contributes to the regulation of TAMs, DCs, helper T lymphocytes, and T cell exhaustion. Furthermore, high expression of FSTL1 promotes the polarization of M1 macrophages to M2 macrophages in GC, which may be one of the reasons for the poor prognosis of FSTL1 in GC patients. For these reasons, FSTL1 might be a promising prognostic biomarker for cancer and could provide new ideas for improving tumor immune evasion and immunotherapy in GC.
Availability of data and materials
The data sets during and/or analyzed during the current study are available from the corresponding author on reasonable request
Tumor Immune Estimation Resource
Real-time quantitative PCR
Gene Expression Profiling Interactive Analysis
The Cancer Genome Atlas
Kyoto Encyclopedia of Genes and Genomes
Database for Annotation, Visualization and Integrated Discovery
Bladder urothelial carcinoma
Breast invasive carcinoma
Cervical squamous cell carcinoma and endocervical adenocarcinoma
Lymphoid neoplasm diffuse large B cell lymphoma
Head and neck squamous cell carcinoma
Kidney renal clear cell carcinoma
Kidney renal papillary cell carcinoma
Acute myeloid leukemia
Brain lower grade glioma
Liver hepatocellular carcinoma
Lung squamous cell carcinoma
Ovarian serous cystadenocarcinoma
Pheochromocytoma and paraganglioma
Skin cutaneous melanoma
Testicular germ cell tumors
Uterine corpus endometrial carcinoma
Jemal A, Center MM, DeSantis C, Ward EM. Global patterns of cancer incidence and mortality rates and trends. Cancer Epidemiol Biomark Prev. 2010;19(8):1893–907.
Matsueda S, Graham DY. Immunotherapy in gastric cancer. World J Gastroenterol. 2014;20(7):1657–66.
Finn OJ. Immuno-oncology: understanding the function and dysfunction of the immune system in cancer. Ann Oncol. 2012;23 Suppl 8(Suppl 8):viii6-viii9.
Zarour HM. Reversing T-cell dysfunction and exhaustion in cancer. Clin Cancer Res. 2016;22(8):1856–64.
Ling Z, Shao L, Liu X, Cheng Y, Yan C, Mei Y, et al. Regulatory T cells and plasmacytoid dendritic cells within the tumor microenvironment in gastric cancer are correlated with gastric microbiota Dysbiosis: a preliminary study. Front Immunol. 2019;10:533.
Noy R, Pollard JW. Tumor-associated macrophages: from mechanisms to therapy. Immunity. 2014;41(1):49–61.
Ngambenjawong C, Gustafson HH, Pun SH. Progress in tumor-associated macrophage (TAM)-targeted therapeutics. Adv Drug Deliv Rev. 2017;114:206–21.
Zhou L, Xu L, Chen L, Fu Q, Liu Z, Chang Y, et al. Tumor-infiltrating neutrophils predict benefit from adjuvant chemotherapy in patients with muscle invasive bladder cancer. Oncoimmunology. 2017;6(4):e1293211.
Zhang H, Liu H, Shen Z, Lin C, Wang X, Qin J, et al. Tumor-infiltrating neutrophils is prognostic and predictive for postoperative adjuvant chemotherapy benefit in patients with gastric cancer. Ann Surg. 2018;267(2):311–8.
Tanaka M, Ozaki S, Osakada F, Mori K, Okubo M, Nakao K. Cloning of follistatin-related protein as a novel autoantigen in systemic rheumatic diseases. Int Immunol. 1998;10(9):1305–14.
Le Luduec JB, Condamine T, Louvet C, Thebault P, Heslan JM, Heslan M, et al. An immunomodulatory role for follistatin-like 1 in heart allograft transplantation. Am J Transplant. 2008;8(11):2297–306.
An J, Wang L, Zhao Y, Hao Q, Zhang Y, Zhang J, et al. Effects of FSTL1 on cell proliferation in breast cancer cell line MDAMB231 and its brain metastatic variant MDAMB231BR. Oncol Rep. 2017;38(5):3001–10.
Liu Y, Tan X, Liu W, Chen X, Hou X, Shen D, et al. Follistatin-like protein 1 plays a tumor suppressor role in clear-cell renal cell carcinoma. Chin J Cancer. 2018;37(1):2.
Gu C, Wang X, Long T, Wang X, Zhong Y, Ma Y, et al. FSTL1 interacts with VIM and promotes colorectal cancer metastasis via activating the focal adhesion signalling pathway. Cell Death Dis. 2018;9(6):654.
Lau MC, Ng KY, Wong TL, Tong M, Lee TK, Ming XY, et al. FSTL1 promotes metastasis and chemoresistance in esophageal squamous cell carcinoma through NFkappaB-BMP signaling cross-talk. Cancer Res. 2017;77(21):5886–99.
Yang W, Wu Y, Wang C, Liu Z, Xu M, Zheng X. FSTL1 contributes to tumor progression via attenuating apoptosis in a AKT/GSK-3beta - dependent manner in hepatocellular carcinoma. Cancer Biomarkers. 2017 Jul 19;20(1):75–85.
Kudo-Saito C. FSTL1 promotes bone metastasis by causing immune dysfunction. Oncoimmunology. 2013 Nov 1;2(11):e26528.
Kudo-Saito C, Ishida A, Shouya Y, Teramoto K, Igarashi T, Kon R, et al. Blocking the FSTL1-DIP2A Axis Improves Anti-tumor Immunity. Cell Rep. 2018 Aug 14;24(7):1790–801.
Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007 Feb;9(2):166–80.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune Cells. Cancer Res. 2017;77(21):e108–e10.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102.
Nagy A, Lanczky A, Menyhart O, Gyorffy B. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Sci Rep. 2018;8(1):9227.
Li B, Severson E, Pignon JC, Zhao H, Li T, Novak J, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 2016;17(1):174.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1(6):417–25.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Fang Z, Zhang L, Liao Q, Wang Y, Yu F, Feng M, et al. Regulation of TRIM24 by miR-511 modulates cell proliferation in gastric cancer. J Exp Clin Cancer Res. 2017;36(1):17.
Chen X, Leung SY, Yuen ST, Chu KM, Ji J, Li R, et al. Variation in gene expression patterns in human gastric cancers. Mol Biol Cell. 2003;14(8):3208–15.
Pages F, Galon J, Dieu-Nosjean MC, Tartour E, Sautes-Fridman C, Fridman WH. Immune infiltration in human tumors: a prognostic factor that should not be ignored. Oncogene. 2010;29(8):1093–102.
Zhang D, He W, Wu C, Tan Y, He Y, Xu B, et al. Scoring System for Tumor-Infiltrating Lymphocytes and Its Prognostic Value for Gastric Cancer. Front Immunol. 2019;10:71.
Sherman-Baust CA, Weeraratna AT, Rangel LB, Pizer ES, Cho KR, Schwartz DR, et al. Remodeling of the extracellular matrix through overexpression of collagen VI contributes to cisplatin resistance in ovarian cancer cells. Cancer Cell. 2003;3(4):377–86.
Qiu S, Deng L, Liao X, Nie L, Qi F, Jin K, et al. Tumor-associated macrophages promote bladder tumor growth through PI3K/AKT signal induced by collagen. Cancer Sci. 2019;110(7):2110–8.
Mantovani A, Sozzani S, Locati M, Allavena P, Sica A. Macrophage polarization: tumor-associated macrophages as a paradigm for polarized M2 mononuclear phagocytes. Trends Immunol. 2002;23(11):549–55.
Yang L, Zhang Y. Tumor-associated macrophages: from basic research to clinical application. J Hematol Oncol. 2017;10(1):58.
Kudo-Saito C, Fuwa T, Murakami K, Kawakami Y. Targeting FSTL1 prevents tumor bone metastasis and consequent immune dysfunction. Cancer Res. 2013;73(20):6185–93.
Nagarsheth N, Wicha MS, Zou W. Chemokines in the cancer microenvironment and their relevance in cancer immunotherapy. Nat Rev Immunol. 2017;17(9):559–72.
Barker HE, Paget JT, Khan AA, Harrington KJ. The tumour microenvironment after radiotherapy: mechanisms of resistance and recurrence. Nat Rev Cancer. 2015;15(7):409–25.
Sun Y. Tumor microenvironment and cancer therapy resistance. Cancer Lett. 2016;380(1):205–15.
Yaacoub K, Pedeux R, Tarte K, Guillaudeux T. Role of the tumor microenvironment in regulating apoptosis and cancer progression. Cancer Lett. 2016;378(2):150–9.
Fridman WH, Galon J, Dieu-Nosjean MC, Cremer I, Fisson S, Damotte D, et al. Immune infiltration in human cancer: prognostic significance and disease control. Curr Top Microbiol Immunol. 2011;344:1–24.
Locati M, Mantovani A, Sica A. Macrophage activation and polarization as an adaptive component of innate immunity. Adv Immunol. 2013;120:163–84.
Shapouri-Moghaddam A, Mohammadian S, Vazini H, Taghadosi M, Esmaeili SA, Mardani F, et al. Macrophage plasticity, polarization, and function in health and disease. J Cell Physiol. 2018;233(9):6425–40.
Komohara Y, Fujiwara Y, Ohnishi K, Takeya M. Tumor-associated macrophages: potential therapeutic targets for anti-cancer therapy. Adv Drug Deliv Rev. 2016;99(Pt B):180–5.
Murray PJ. Macrophage Polarization. Annu Rev Physiol. 2017;79:541–66.
Gardner A, Ruffell B. Dendritic cells and cancer immunity. Trends Immunol. 2016;37(12):855–65.
Marshall EA, Ng KW, Kung SH, Conway EM, Martinez VD, Halvorsen EC, et al. Emerging roles of T helper 17 and regulatory T cells in lung cancer progression and metastasis. Mol Cancer. 2016;15(1):67.
Tang M, Diao J, Cattral MS. Molecular mechanisms involved in dendritic cell dysfunction in cancer. Cell Mol Life Sci. 2017;74(5):761–76.
Chen LJ, Zheng X, Shen YP, Zhu YB, Li Q, Chen J, et al. Higher numbers of T-bet(+) intratumoral lymphoid cells correlate with better survival in gastric cancer. Cancer Immunol Immunother. 2013;62(3):553–61.
Li Q, Li Q, Chen J, Liu Y, Zhao X, Tan B, et al. Prevalence of Th17 and Treg cells in gastric cancer patients and its correlation with clinical parameters. Oncol Rep. 2013;30(3):1215–22.
Chen W, Konkel JE. Development of thymic Foxp3(+) regulatory T cells: TGF-beta matters. Eur J Immunol. 2015;45(4):958–65.
Perrone G, Ruffini PA, Catalano V, Spino C, Santini D, Muretto P, et al. Intratumoural FOXP3-positive regulatory T cells are associated with adverse prognosis in radically resected gastric cancer. Eur J Cancer. 2008;44(13):1875–82.
Chi LJ, Lu HT, Li GL, Wang XM, Su Y, Xu WH, et al. Involvement of T helper type 17 and regulatory T cell activity in tumour immunology of bladder carcinoma. Clin Exp Immunol. 2010;161(3):480–9.
Mempel TR, Pittet MJ, Khazaie K, Weninger W, Weissleder R, von Boehmer H, et al. Regulatory T cells reversibly suppress cytotoxic T cell function independent of effector differentiation. Immunity. 2006;25(1):129–41.
Blackburn SD, Shin H, Haining WN, Zou T, Workman CJ, Polley A, et al. Coregulation of CD8+ T cell exhaustion by multiple inhibitory receptors during chronic viral infection. Nat Immunol. 2009;10(1):29–37.
Fuertes Marraco SA, Neubert NJ, Verdeil G, Speiser DE. Inhibitory Receptors Beyond T Cell Exhaustion. Front Immunol. 2015;6:310.
Wherry EJ, Kurachi M. Molecular and cellular insights into T cell exhaustion. Nat Rev Immunol. 2015;15(8):486–99.
Dong Y, Geng Y, Li L, Li X, Yan X, Fang Y, et al. Blocking follistatin-like 1 attenuates bleomycin-induced pulmonary fibrosis in mice. J Exp Med. 2015;212(2):235–52.
Chae YK, Arya A, Iams W, Cruz MR, Chandra S, Choi J, et al. Current landscape and future of dual anti-CTLA4 and PD-1/PD-L1 blockade immunotherapy in cancer; lessons learned from clinical trials with melanoma and non-small cell lung cancer (NSCLC). J Immunother Cancer. 2018;6(1):39.
Xia L, Liu Y, Wang Y. PD-1/PD-L1 Blockade therapy in advanced non-small-cell lung cancer: current status and future directions. Oncologist. 2019;24(Suppl 1):S31–41.
This work was supported in part by grants from the National Natural Science Foundation of China (grant number: 81860427, 81760432, 81660402 and 81660405), the JiangXi Province General Project (grant number: 20164BCD40097, 20171BBG70121, and 2017BBH80027), the Youth Scientific Funds-Youth Fund Project (grant number: 2018ACB21037 and 20171BAB215041), the JiangXi Province Education Fund Project (grant number: 700653002), the department of health of JiangXi Province Project (grant number: 20181041), and the JiangXi Province Postgraduate Special Innovation Fund (grant number: YC2019-B034).
Ethics approval and consent to participate
The study was approved by the First Affiliated Hospital of Nanchang University Ethics Committee, Nanchang, China (No.126), and all patients who underwent hepatectomy in this study signed tissue sample consent for clinical research.
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Supplementary Figure 1. The FSTL1 overall survival curves of GC patients in KM-plotter (A) and GEPIA (B) database. FSTL1, follistatin-like 1; OS, overall survival, HR, hazard ratio.
Additional file 2:
Supplementary Table 1. Correlation analysis between FSTL1 and the related gene markers of monocyte and macrophages by GEPIA.
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.
About this article
Cite this article
Li, L., Huang, S., Yao, Y. et al. Follistatin-like 1 (FSTL1) is a prognostic biomarker and correlated with immune cell infiltration in gastric cancer. World J Surg Onc 18, 324 (2020). https://doi.org/10.1186/s12957-020-02070-9
- Tumor immune infiltrates
- Prognostic biomarker
- Gastric cancer