Skip to main content
  • Research article
  • Open access
  • Published:

ITGA5 is a prognostic biomarker and correlated with immune infiltration in gastrointestinal tumors

Abstract

Background

Integrin Subunit Alpha 5 (ITGA5), belongs to the integrin alpha chain family, is vital for promoting cancer cell invasion, metastasis. However, the correlation between ITGA5 expression and immune infiltration in gastrointestinal tumors remain unclear.

Methods

The expression level of ITGA5 was detected by Oncomine and Tumor Immune Estimation Resource (TIMER). The association between ITGA5 and prognosis of patients was identified by Kaplan–Meier plotter, Gene Expression Profiling Interactive Analysis 2 (GEPIA2) and PrognoScan. We evaluated the correlation between ITGA5 expression and immune infiltrating level via TIMER. Besides, TIMER, immunohistochemistry (IHC) staining and western blot were used to explore correlations between ITGA5 expression and markers of immune infiltrates cells. Furthermore, we constructed protein-protein interaction (PPI) network and performed functional enrichment by GeneMANIA and Metascape.

Results

ITGA5 was generally overexpressed and correlated with worse prognosis in multiple types of gastrointestinal tumors. In addition, ITGA5 expression level was significantly associated with tumor purity and immune infiltration levels of different immune cells in gastrointestinal tumors. Interestingly, immune markers for monocytes, tumor - associated macrophages (TAMs), macrophages 2 (M2) cells and T-helper 2 (Th2) cells were found to be significantly and positively correlated with ITGA5 expression levels in colon and gastric cancer. Results from IHC staining and western blot further proved that markers of Th2 and M2 cell were significantly increased in gastric cancer patients with high ITGA5 expression levels. Lastly, interaction network and function enrichment analysis revealed ITGA5 was mainly involved in “integrin mediated signaling pathway”, “leukocyte migration”, “cell-substrate adhesion”.

Conclutions

Our study demonstrated that ITGA5 may act as an essential regulator of tumor immune cell infiltration and a valuable prognostic biomarker in gastrointestinal tumors. Additional work is needed to fully elucidate the underlying mechanisms behind these observations.

Peer Review reports

Background

Gastrointestinal tumors are a common malignant tumor affecting individuals worldwide [1]. The high mortality rate associated with gastrointestinal tumors is a result of the tumor invading and metastasizing to other tissues. The increasing resistance observed in gastrointestinal tumors receiving chemotherapy and targeted therapy has led to an increased challenge when treating this disease [2,3,4]. There is extensive evidence that tumor-infiltrating lymphocytes, such as tumor - associated macrophages (TAMs), influence both the prognosis and efficacy of chemotherapy and immunotherapy [5, 6], as well as tumor angiogenesis, progression and metastasis [7]. Therefore, it is crucial to identify novel immune-related therapeutic targets and elucidate the potential mechanisms behind the immune interactions linked to gastrointestinal tumors.

Integrins are a large family of heterodimeric integral membrane proteins that function in cell surface adhesion and signaling pathways. Integrin Subunit Alpha 5 (ITGA5), primarily binds to Integrin Subunit Beta 1 (ITGB1) to form a α5β1 heterodimer [8, 9]. Recently, it has become clear that ITGA5 is essential for cancer proliferation, migration, invasion and metastasis [10,11,12,13]. In addition, ITGA5 was shown to maintain cancer cell stemness and chemotherapy resistance [14, 15]. Interestingly, ITGA5 was identified as a critical factor for promoting the differentiation of human mesenchymal stromal cells in vivo [16]. Recent studies revealed that ITGA5 plays a crucial role in different tumor cell subgroups and cell-cell interactions. Importantly, there is evidence revealing that immune components existing in the tumor microenvironment (TME) comprehensively regulate the biological behavior of the tumor through mutual interactions [17]. Furthermore, immunotherapy, especially the application of immune checkpoint modulators, is an efficient treatment strategy for solid tumors [18, 19].

However, the important role of ITGA5 in gastrointestinal tumors and its relationship with tumor immunity are still unclear. Here, we explored the expression level, prognostic value, protein interaction network and enrichment analysis of ITGA5 in gastrointestinal tumors. More importantly, the correlation between ITGA5 and tumor-infiltrating immune cells and markers in different tumor microenvironments was also analyzed. Our result supported the pro-oncogenic effect of ITGA5 and revealed the overexpression of ITGA5 may have a potential and important relationship with tumor-immune infiltration.

Methods

Clinical tissue specimens

The tissue samples used for immunohistochemistry (IHC) staining and western blot were collected from patients diagnosed with gastric cancer at the First Affiliated Hospital of Anhui Medical University. All patients have already provided written informed consent. This work was approved by the Academic Committee of The First Affiliated Hospital of Anhui Medical University (certification no. Quick-PJ 2020-11-14) and was conducted following the Declaration of Helsinki.

Oncomine analysis

Expression of ITGA5 mRNA levels were analyzed to compare normal and tumor tissues in various cancer types using the Oncomine database (https://www.oncomine.org) [20, 21]. The filter condition was set to a P value < 0.001, fold change > 1.5, gene rank: 10%, data type: mRNA.

PrognoScan analysis

PrognoScan database (http://www.abren.net/PrognoScan) was used to identify the correlation between ITGA5 expression and prognosis in various cancers across a large collection of publicly available microarray datasets [22]. The threshold was adjusted to a Cox P-value < 0.05.

Kaplan-Meier plotter analysis

Kaplan-Meier plotter (https://kmplot.com) possesses gene expression data and survival information associated with 1065 gastric cancer patients [23]. To evaluate differential expression of ITGA5 mRNA levels and the effects on patient progression, Kaplan–Meier survival curves for overall survival (OS) and disease or progression free survival (DFS or PFS) were generated for patients containing low and high ITGA5 mRNA expression levels. In addition, we used the Kaplan-Meier Plotter to evaluate the relationship between the expression of ITGA5 and the clinicopathological characteristics of gastric cancer patients. Results contained a hazard ratio (HR) with a 95% confidence interval (CI) and log rank P value.

TIMER analysis

Tumor Immune Estimation Resource (TIMER) database (https://cistrome.shinyapps.io/timer) is a publicly available comprehensive resource for systematic analysis of tumor immune-infiltrates [24]. It includes 10,897 samples across 32 different cancer types from The Cancer Genome Atlas (TCGA). The correlation between ITGA5 expression levels and immune infiltrates including B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages and dendritic cells, were explored using the gene module in various cancer types. The correlation between genes markers of tumor-infiltrating immune cells and ITGA5 expression was analyzed using the correlation module. Gene markers included CD8+ T cells, general T cells, B cells, monocytes, TAMs, M1 macrophages, M2 macrophages, neutrophils, natural killer (NK) cells, dendritic cells (DCs), T-helper 1 (Th1) cells, T-helper 2 (Th2) cells, follicular helper T (Tfh) cells, T-helper 17 (Th17) cells, Tregs and exhausted T cells. ITGA5 expression was plotted on the X-axis and the expression of related marker genes were represented as gene symbols on the Y-axis. Correlation coefficients were estimated using Spearman’s correlation method. Gene expression levels were shown as log2 RSEM.

GEPIA2 analysis

Gene expression Profiling Interactive Analysis 2 (GEPIA2) (http://gepia2.cancer-pku.cn) is an updated and enhanced version of GEPIA that supports 198,619 isoforms and 84 cancer subtypes [25, 26]. In our study, the “Expression DIY” and “Survival” module of GEPIA2 were used to determine the correlation between ITGA5 expression and prognosis in pancreatic and esophageal cancers.

Immunohistochemistry staining analysis

Collected tissue specimens were formalin fixed and embedded with paraffin. Tissue sections (4 μm thickness) were used in IHC staining analyses as previously described [27]. The intensity of immunostaining and the percentage of positive cells were used to determine the immunoreactivity to proteins. The density grading of intensity was analyzed as follows: 0:no staining, 1:weak staining; 2: medium staining, and 3: strong staining. The proportion of positive cells was further classified as 0: < 1%; 1: 1–30%, 2: 30–70%, and 3:> 70%. Semi-quantitative assessment and final immunostaining scores were determined by multiplying the proportion score and intensity scores. Scores 0, 1, 2 and 3 indicated low expression levels, while scores 4, 6 and 9 indicated high expression levels. IHC results were evaluated by two pathologists using a blind test and differences were resolved through consensus. Antibody information is listed in Additional file: Table S1.

Western blot analysis

All the collected tissues were cut and ground into a homogenate. Using RIPA buffer (Radio-Immunoprecipitation assay buffer, Beyotime) to lyse aggregated proteins in the cells. In addition, the concentrations were detected with the help of the BCA Protein Assay Kit (Beyotime). Segregation of equivalent quantities of proteins was made using the SDS-PAGE, followed by transferring to PVDF membranes (Millipore, Boston, MA, USA). After that, the membrane was incubated with the targeted primary antibody overnight at 4 °C. After washing, incubation of the membranes was performed with the corresponding secondary antibody for 2 h at the room temperature. GAPDH was used as an internal loading control. The findings were analyzed with the help of the ECL detection system (Pierce Biotech, Rockford, IL, USA). Antibody information is listed in Additional file: Table S2.

GeneMANIA analysis

GeneMANIA (http://genemania.org/) is used to construct a protein-protein interaction (PPI) network and analyze the function of interactive genes [28]. This online tool uses bioinformatic methods to display a list of interacting genes, including gene co-expression, physical interactions, gene co-localization, gene enrichment analysis and website prediction. GeneMANIA was used to construct the PPI network for ITGA5.

Metascape analysis

Metascape (https://metascape.org/gp/index.html) is a gene function annotation website with obvious advantages of wide coverage and fast updating [29]. It integrates multiple authoritative data resources such as Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, UniProt and DrugBank to complete thorough pathway enrichment and biological process annotation. Genes interacting with ITGA5 were put into Metascape to perform GO and KEGG annotation.

Statistical analysis

All statistical analyses were performed using SPSS 21.0 software (IBM Corp; Armonk, NY). Survival curves and relative results generated from PrognoScan, Kaplan-Meier Plotter and GEPIA2 databases were shown with HR and P or Cox P values from a log-rank test. Correlation analyses in TIMER were evaluated using Spearman’s correlation. Associations between ITGA5 expression and IHC results in patients were evaluated using Pearson’s χ2 test. Student’s t-test was used to determine the significance for the western blot experiments. P-value of < 0.05 was considered as statistically significant.

Results

Pan-cancer analysis of ITGA5 expression levels

To determine ITGA5 mRNA expression levels in both normal and tumor tissues, Oncomine was interrogated to analyze ITGA5 expression in various cancer types. The expression levels of ITAG5 were greater in colorectal, esophageal, gastric and pancreatic cancers relative to normal tissues (Fig. 1a). In addition, the transcriptional levels for ITGA5 were analyzed using RNA-seq data for multiple malignancies in the TCGA with the use of TIMER. Results revealed significant differences in ITGA5 expression levels when comparing normal and tumor tissues. Hepatocellular carcinoma (LIHC) and gastric adenocarcinoma (STAD) were significantly higher compared with normal tissue levels (Fig. 1b).

Fig. 1
figure 1

ITGA5 expression in different types of human cancers. a Different expression of ITGA5 between tumors and normal tissues in Oncomine database. The filter condition was set as: p value < 0.0001, fold change > 2, gene rank: 10%, data type: mRNA. b Different expression of ITGA5 between tumors and normal tissues in TIMER database (*P < 0.05, **P < 0.01, ***P < 0.001)

Correlation between ITGA5 expression levels and patient prognosis

To further identify the prognostic potential of ITGA5 in various cancers, clinical data from TCGA were used to explore the prognostic value of ITGA5 expression level using GEPIA2. High expression of ITGA5 was marginally correlated with poor prognosis in pancreatic cancer (OS HR = 1.6, P = 0.12; DFS HR = 2.1, P = 0.015) (Fig. 2a-b) and esophageal cancer (OS HR = 1.4, P = 0.37; DFS HR = 1.2, P = 0.52) (Fig. 2c-d). Subsequently, the Kaplan-Meier plotter database was used to determine the prognostic value of ITGA5. Significant differences were observed between ITGA5 expression and prognosis in gastric cancer (OS HR = 2.4, 95% CI = 1.96 to 2.97, P < 1e-16; DFS HR = 2.69, 95% CI = 2.04 to 3.54, P = 2.2e-13) (Fig. 2e-f) and live cancer (OS HR = 1.73, 95% CI = 1.22 to 2.44, P = 0.0017; DFS HR = 1.23, 95% CI = 0.88 to 1.72, P = 0.23) patients (Fig. 2g-h). In addition, the PrognoScan database was used to investigates the relationship between ITGA5 mRNA levels and the survival of cancer patients using high-throughput analysis and detailed clinical prognosis data. ITGA5 expression was found to significantly impact prognosis in colorectal. Specifically, the cohorts (GSE17536) [30] included 177 samples at different stages of colorectal cancer and showed a remarkable association between high ITGA5 expression levels and poor prognosis (OS HR = 1.66, 95% CI = 1.16 to 2.37, Cox P = 0.005; DFS HR = 2.97, 95% CI = 1.70 to 5.21, Cox P = 0.0004) (Figs. 2i–j) in colorectal cancer.

Fig. 2
figure 2

The prognostic value of ITGA5 expression in gastrointestinal cancers. a-d Correlation between ITGA5 expression and prognosis of pancreatic cancer and esophageal cancer in GEPIA2; e-h Correlation between ITGA5 expression and prognosis of gastric cancer and liver cancer in Kaplan-Meier Plotter; i-j Correlation between ITGA5 expression and prognosis of colorectal cancer in PrognoScan. OS, overall survival; DFS, disease free survival; HR: Hazard ratio

Interestingly, we also found that ITGA5 also has a significant impact on the prognosis of patients in some non-gastrointestinal tumors. High expression of ITGA5 was also associated with poor prognosis in lung cancer (OS HR = 1.6, 95% CI = 1.4 to 1.81, P = 5.9e-13; DFS HR = 1.75, 95% CI = 1.44 to 2.12, P = 1.2e-08) and breast cancer (OS HR = 1.21, 95% CI = 0.96 to 1.52, P = 0.11; DFS HR = 1.28, 95% CI = 1.14 to 1.44, P = 4.1e-05) (Additional file: Fig. S1). Similarly, the cohorts (GSE26712) [31] showed that high ITGA5 expression levels were correlated with poor prognosis (OS HR = 1.58, 95% CI = 1.06 to 2.36, Cox P = 0.025; DFS HR = 1.54, 95% CI = 1.06 to 2.24, Cox P = 0.025) in ovarian cancers (Additional file: Fig. S1).

Relationship between ITGA5 expression and immune infiltration levels in gastrointestinal cancers

Recent work has shown that gastrointestinal tumors exhibit extensive immune infiltration characteristics [32,33,34]. Therefore, to better understand the role of ITGA5 in gastrointestinal tumors, the relationship between ITGA5 expression and immune infiltration in gastrointestinal tumors was investigated using the TIMER database. ITGA5 expression was found to be significantly correlated with tumor purity in colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), STAD and rectum adenocarcinoma (READ) (P < 0.05), CD8+ cell in COAD, LIHC and pancreatic adenocarcinoma (PAAD), CD4 + T cells in COAD, PAAD, LIHC, READ and STAD (P < 0.05) (Fig. 3a-f). Furthermore, ITGA5 expression was significantly associated with macrophages and dendritic cells in COAD, ESCA, READ, LIHC, PAAD and STAD (P < 0.05) (Fig. 3a-f). ITGA5 expression was also significantly associated with neutrophils in COAD, LIHC, READ, PAAD and STAD (P < 0.05) (Fig. 3a-f).

Fig. 3
figure 3

Correlation of ITGA5 expression with immune infiltration level in Gastrointestinal tumors. a Colon adenocarcinoma (COAD); b Esophageal carcinoma (ESCA); c Liver hepatocellular carcinoma (LIHC); d Pancreatic adenocarcinoma (PAAD); e Rectum adenocarcinoma (READ); f Stomach adenocarcinoma (STAD)

Specifically, ITGA5 expression levels showed significantly positive correlations with infiltrating levels of CD8 + T cells (r = 0.248, P = 4.18e-07), CD4+ T cells (r = 0.423, P = 6.67e-19), macrophages (r = 0.545, P = 1.20e-32), neutrophils (r = 0.590, P = 4.92e-39) and dendritic cells (r = 0.556, P = 5.48e-34) in COAD (Fig. 3a). Similarly, positive correlations were identified with infiltrating levels of CD4+ T cells (r = 0.426, P = 1.52e-17), macrophages (r = 0.447, P = 1.42e-19), neutrophils (r = 0.189, P = 2.47e-04) and dendritic cells (r = 0.28, P = 4.19e-08) in STAD (Fig. 3f). These results revealed that ITGA5 expression is closely associated with immune infiltration in gastrointestinal tumors.

Correlation between ITGA5 expression levels and immune markers

Gastric cancer and colorectal cancer are the most common gastrointestinal-related malignancies. Therefore, COAD and STAD were chosen to further investigate the relationship between ITGA5 expression and immune marker genes of diverse immune cells using TIMER databases. The immune cells analyzed included CD8+ T cells, general T cells, different functional T cells, B cells, monocytes, TAMs, M1 macrophages and M2 macrophages, neutrophils, NK cells and dendritic cells. After adjusting correlation by purity, ITGA5 expression was found to be significantly correlated with immune markers of most immune cells in STAD and COAD (Table 1 and Fig. 4a-h). Immune markers for monocytes, TAMs, M2 and Th2 phenotypes were found to be significantly and positively correlated with ITGA5 expression levels in COAD and STAD (P < 0.05; Figs. 4a-h).

Table 1 Correlation analysis between ITGA5 and markers of immune cells in TIMER database
Fig. 4
figure 4

Correlation analysis between ITGA5 expression and immune marker set in STAD and COAD. a Scatterplots of correlations between ITGA5 expression and gene markers of Th2 cells in COAD; b Scatterplots of correlations between ITGA5 expression and gene markers of Monocyte in COAD; c Scatterplots of correlations between ITGA5 expression and gene markers of Th2 cells in STAD; d Scatterplots of correlations between ITGA5 expression and gene markers of Monocyte in STAD; e Scatterplots of correlations between ITGA5 expression and gene markers of TAM in COAD; f Scatterplots of correlations between ITGA5 expression and gene markers of M2 in COAD; g Scatterplots of correlations between ITGA5 expression and gene markers of TAM in STAD; h Scatterplots of correlations between ITGA5 expression and gene markers of M2 cells in STAD

High ITGA5 expression affects the prognosis of gastric cancer patients exhibiting lymph node metastasis

To better understand the relationship between ITGA5 expression levels and clinicopathological features, further research was focused on gastric cancer. According to Kaplan-Meier Plotter database analyses, overexpression of ITGA5 is strongly correlated with the deterioration of OS and PFS based on gender, differentiation and Lauren classification (P < 0.05) (Table 2). In addition, stage I through IV cancer patients with high expression levels of ITGA5 showed worse OS (P < 0.05) (Table 2). Similarly, high expression levels of ITGA5 were associated with worse PFS in stage II through IV gastric cancer patients (P < 0.05) but not in stage I gastric cancer patients (Table 2). However, there was no strong correlation observed between ITGA5 expression levels and patient prognosis in stage N0 (OS HR = 1.92, P = 0.13; PFS HR = 1.95, P = 0.11), mixed Lauren classification (OS HR = 1.88, P = 0.22; PFS HR = 0.48, P = 0.16) or poor differentiation (OS HR = 0.75, P = 0.16; PFS HR = 1.39, P = 0.25) (Table 2). Moreover, high ITGA5 levels express the highest HR value of N1 with OS and PFS among four N categories (Table 2). These results indicate that the expression levels of ITGA5 are related to lymph node metastasis in gastric cancer patients.

Table 2 Correlation of ITGA5 and clinicopathological factors in STAD from Kaplan-Meier Plotter database

ITGA5 protein expression and M2 and Th2 immune marker genes in STAD

To confirm the results obtained from publicly available databases, IHC were performed on adjacent normal and STAD tumor tissues obtained from 40 gastric patients. Compared with adjacent normal tissues, ITGA5, CD163, STAT6 and GATA3 levels were greater in most tumor tissues (Fig. 5a). In addition, these results revealed that the expression levels of CD163 (χ2 = 8.750, P = 0.003), STAT6 (χ2 = 8.174, P = 0.004) and GATA3 (χ2 = 5.079, P = 0.024) in the ITGA5 high-expression group were significantly higher than what was observed in the ITGA5 low-expression group (Table 3). Among the patients with higher ITGA5 expression levels, 12 patients showed higher CD163, STAT6 and GATA3 levels at the same time. In addition, we attempted to detect the proteins levels of ITGA5, CD163, STAT6 and GATA3 in patient tissues by western blot. As Fig. 5b, c shows, the expression levels in tumor tissues were indeed higher than that in adjacent normal tissues.

Fig. 5
figure 5

IHC and western blot analysis between ITGA5 and different markers of Th2 and M2 cells in gastric cancer tissues and adjacent normal tissues. a Representative IHC images of CD163, ITGA5, GATA3 and STAT6 in gastric cancer tissues and adjacent normal tissues; b Representative western blot images of CD163, ITGA5, GATA3 and STAT6 in 6 paired gastric cancer tissues (T) and adjacent normal tissues (N) (The full-length, original blots are presented in Additional file: Fig. S2); c Relative protein density of CD163, ITGA5, GATA3 and STAT6 in 6 paired gastric cancer tissues and adjacent normal tissues

Table 3 The protein expression level of ITGA5 and immune marker genes of M2 and Th2 cells in STAD

ITGA5 PPI network and enrichment analyses

PPI network analysis revealed interactions between ITGA5 and specific genes. As shown by GeneMANIA, genes interacting with ITGA5 included SPP1, ITGB1, ITGB3, COL18A1, HOXD3, ANGPTL3, CD9, FBN1, VEGFD, ITGA2B, ITGA4, FN1, FLT4, TLN1, ACVRL1, FUBP1, THBS4, ACTN1, ANGPT2 and ITGA8 (Fig. 6a).

Fig. 6
figure 6

PPI network of ITGA5 and functional enrichment analysis. a PPI network of ITGA5 in GeneMANIA, different colors of the network edge indicate the bioinformatics methods applied: physical interactions, co-expression, predicted, co-localization, pathway, genetic interactions, and shared protein domains; b A heat map of GO and KEGG analysis of ITGA5 and it’s 20 most interactive genes, orange represents the enrichment terms colored by p values; c,d Interactive network of the top enrichment terms colored by cluster ID, different colors represent various enrichment pathways of ITGA5 correlated genes

To further predict enriched function and pathway information of ITGA5 interacting genes, GO function and KEGG pathway analyses were performed using Metascape. According to GO term analysis, genes interacting with ITGA5 in the PPI network are mainly related with cell-substrate adhesion, regulation of leukocyte migration, response to wounding and myeloid cell differentiation. KEGG pathway analysis showed these genes are enriched in focal adhesion, PID avb3 integrin pathway and PID integrin5 pathway. Based on these results, ITGA5 and its interacting proteins were shown to play an important role in the integrin mediated signaling pathway, leukocyte migration, cell-substrate adhesion and other essential biological processes (Fig. 6b). In addition, Cytoscape was used to determine the relationship for the enriched terms and to build a network diagram (Fig. 6c and d).

Discussion

ITGA5 belongs to the integrin alpha chain family and interacts with ITGB1 to form a heterodimeric integral membrane protein - integrin α5β1 complex that performs diverse biological functions. Despite not knowing all of the functions performed by ITGA5, multiple studies have demonstrated that increased expression of ITGA5 is associated with poor prognosis in multiple tumors types, such as triple negative breast cancer [14], lung cancer [35], hepatocellular carcinoma [36] and ovarian cancer [37]. Here, public databases were used to determine ITGA5 expression levels in different cancer types. Compared to normal tissues, ITGA5 was generally overexpressed in different cancers including multiple types of gastrointestinal tumors. We explored whether the expression of ITGA5 was correlated with patient prognosis in these different tumors. Results revealed that high expression levels of ITGA5 were significantly correlated with worse OS and DFS in four types of gastrointestinal tumors, including colorectal, pancreatic, gastric and liver cancers. Although the results did not show statistical significance, esophageal cancer patients with high ITGA5 expression levels still showed a trend of high survival risk. To understand the clinical significance of ITGA5 in gastric cancer, the prognostic values for ITGA5 expression levels with various clinicopathological factors were determined. High expression of ITGA5 was correlated with poor prognosis in gastric cancer stages II-IV, T2-T4, N1-N3 and M0-M1. These results demonstrated that ITGA5 has potential as a prognostic biomarker in gastrointestinal tumors.

Both innate and adaptive immune cells regulate the biological behavior of tumors and the response to treatment through direct contact or signal transduction mechanisms [17]. We therefore explored the relationship between ITGA5 expression levels and immune infiltration levels. Our results indicated that ITGA5 expression was significantly and negatively correlated with immune purity in four types of gastrointestinal tumors, including COAD, ECAD, STAD and READ. In STAD and COAD, significant and positive relationships were observed between ITGA5 expression and infiltration levels for CD4 + T cells, macrophages, neutrophils and dendritic cells. In addition, the relationship between ITGA5 expression levels and immune cell markers implicates underlying regulation in STAD and COAD. First, gene markers for Th2 cells (GATA3, STAT6, STAT5A and IL13) showed a positive correlation with ITGA5 expression. The M2 macrophage markers (CD163, VSIG4 and MS4A4A) showed a moderate or strong correlation. In addition, we also observed a significant correlation between ITGA expression and monocytes and TAMs in COAD and STAD. These results revealed a potential regulatory role for ITGA5 in the polarization of macrophages. TAMs, derived from mononuclear cells, are induced to differentiate into M2 cells by cytokines such as IL-4, IL-10 and IL-13, which are secreted by Th2 cells [38,39,40]. Moreover, M2 cells exhibit protumoral activity by promoting genetic instability, angiogenesis, stem cell nurturing and local immunosuppression [41, 42]. Subsequently, IHC and western blot were performed on adjacent normal and tumor tissues extracted from gastric cancer patients to find that Th2 cell markers such as STAT6, GATA3 and the M2 cell polarization marker CD163 were significantly increased in patients with high ITGA5 expression levels. This suggested that ITGA5 may play an important role in promoting the recruitment and activation of monocytes, TAMs, Th2 cells and M2 cells in the tumor microenvironment.

Emerging evidence provides potential mechanisms to explain the correlation between ITGA5 expression levels and immune infiltration. Integrin α5β1 usually acts as a receptor for fibronectin to perform biological functions. Previous studies mainly focused on the role of ITGA5 in tumor cells, but recent studies found that ITGA5 is also expressed in cancer-associated fibroblasts [43], TAMs [44] and chimeric antigen receptor expressing T cells [45]. Based on previous research, cancer-associated fibroblasts, which are regulators of immune cell recruitment and function, affect both innate and adaptive immune responses [46, 47]. Similarly, ITGA5 is expressed on the surface of TAMs and chimeric antigen receptor expressing T cells and can directly regulate the recruitment and alternative activation of immune cells. Interestingly, further work has shown that decreased expression of ITGA5 down-regulates pancreatic stellate cells differentiation into cancer-associated fibroblasts in pancreatic ductal cancer [15]. Importantly, pancreatic stellate cell are seen as the main source of cancer-associated fibroblasts [48]. The study presented here revealed an important role for ITGA5 in the differentiation and maturation of cells. These results partly explain the effects of ITGA5 expression on immune cell infiltration in malignant tumors and provides insight for further work. ITGA5 and its interacting proteins play an important role in integrin-mediated signaling pathways, leukocyte migration and cell-matrix adhesion as shown by the PPI network that was constructed as well as gene functional enrichment analysis. ITGA5 may act as an “anchor” for cell positioning, promoting the aggregation, adhesion and migration of partial immune cells and changing components of the tumor microenvironment. This may result in immune escape and suppression, ultimately accelerating tumor development and metastasis. In the future, the possible role and mechanisms behind ITAG5 in tumor immunity should be explored.

Conclusions

In summary, our study demonstrated that ITGA5 may be an essential regulator of tumor immune cell infiltration and a valuable prognostic biomarker in gastrointestinal tumors. Additional work is needed to fully elucidate the underlying mechanisms behind these observations.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ITGA5:

Integrin subunit alpha 5

TAMs:

Tumor - associated macrophages

ITGB1:

Integrin subunit beta 1

TME:

Tumor microenvironment

PPI:

Protein-protein interaction

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

OS:

Overall survival

DFS:

Disease free survival

PFS:

Progression free survival

HR:

Hazard ratio

CI:

Confidence interval

TCGA:

The cancer genome atlas

NK:

Natural killer

DCs:

Dendritic cells

Th1:

T-helper 1

Th2:

T-helper 2

Tfh:

Follicular helper T

Th17:

T-helper 17

GEPIA2:

Gene expression profiling interactive analysis 2

IHC:

Immunohistochemistry

LIHC:

Liver hepatocellular carcinoma

STAD:

Stomach adenocarcinoma

COAD:

Colon adenocarcinoma

ESCA:

Esophageal carcinoma

READ:

Rectum adenocarcinoma

PAAD:

Pancreatic adenocarcinoma

References

  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  PubMed  Google Scholar 

  2. Van Cutsem E, Sagaert X, Topal B, Haustermans K, Prenen H. Gastric cancer. Lancet (London, England). 2016;388(10060):2654–64.

    Article  CAS  Google Scholar 

  3. Becht E, de Reyniès A, Giraldo NA, Pilati C, Buttard B, Lacroix L, et al. Immune and stromal classification of colorectal Cancer is associated with molecular subtypes and relevant for precision immunotherapy. Clin Cancer Res. 2016;22(16):4057–66.

    Article  CAS  PubMed  Google Scholar 

  4. Myint ZW, Goel G. Role of modern immunotherapy in gastrointestinal malignancies: a review of current clinical progress. J Hematol Oncol. 2017;10(1):86.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Waniczek D, Lorenc Z, Śnietura M, Wesecki M, Kopec A, Muc-Wierzgoń M. Tumor-associated macrophages and regulatory T cells infiltration and the clinical outcome in colorectal Cancer. Arch Immunol Ther Exp. 2017;65(5):445–54.

    Article  CAS  Google Scholar 

  6. 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.

    Article  PubMed  Google Scholar 

  7. Spranger S, Bao R, Gajewski TF. Melanoma-intrinsic β-catenin signalling prevents anti-tumour immunity. Nature. 2015;523(7559):231–5.

    Article  CAS  PubMed  Google Scholar 

  8. Morgan MR, Byron A, Humphries MJ, Bass MD. Giving off mixed signals--distinct functions of alpha5beta1 and alphavbeta3 integrins in regulating cell behaviour. IUBMB Life. 2009;61(7):731–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Marelli UK, Rechenmacher F, Sobahi TR, Mas-Moruno C, Kessler H. Tumor targeting via integrin ligands. Front Oncol. 2013;3:222.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Tani N, Higashiyama S, Kawaguchi N, Madarame J, Ota I, Ito Y, et al. Expression level of integrin alpha 5 on tumour cells affects the rate of metastasis to the kidney. Br J Cancer. 2003;88(2):327–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Sawada K, Mitra AK, Radjabi AR, Bhaskar V, Kistner EO, Tretiakova M, et al. Loss of E-cadherin promotes ovarian cancer metastasis via alpha 5-integrin, which is a therapeutic target. Cancer Res. 2008;68(7):2329–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Qian F, Zhang ZC, Wu XF, Li YP, Xu Q. Interaction between integrin alpha (5) and fibronectin is required for metastasis of B16F10 melanoma cells. Biochem Biophys Res Commun. 2005;333(4):1269–75.

    Article  CAS  PubMed  Google Scholar 

  13. Deng Y, Wan Q, Yan W. Integrin α5/ITGA5 promotes the proliferation, migration, invasion and progression of Oral squamous carcinoma by epithelial-Mesenchymal transition. Cancer Manag Res. 2019;11:9609–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Xiao Y, Li Y, Tao H, Humphries B, Li A, Jiang Y, et al. Integrin α5 down-regulation by miR-205 suppresses triple negative breast cancer stemness and metastasis by inhibiting the Src/Vav2/Rac1 pathway. Cancer Lett. 2018;433:199–209.

    Article  CAS  PubMed  Google Scholar 

  15. Kuninty PR, Bansal R, De Geus SWL, Mardhian DF, Schnittert J, van Baarlen J, et al. ITGA5 inhibition in pancreatic stellate cells attenuates desmoplasia and potentiates efficacy of chemotherapy in pancreatic cancer. Sci Advances. 2019;5(9):eaax2770.

    Article  CAS  Google Scholar 

  16. Hamidouche Z, Fromigué O, Ringe J, Häupl T, Vaudin P, Pagès JC, et al. Priming integrin alpha5 promotes human mesenchymal stromal cell osteoblast differentiation and osteogenesis. Proc Natl Acad Sci U S A. 2009;106(44):18587–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wu T, Dai Y. Tumor microenvironment and therapeutic response. Cancer Lett. 2017;387:61–8.

    Article  CAS  PubMed  Google Scholar 

  18. Valsecchi ME. Combined Nivolumab and Ipilimumab or Monotherapy in untreated melanoma. N Engl J Med. 2015;373(13):1270.

    Article  PubMed  Google Scholar 

  19. Borghaei H, Paz-Ares L, Horn L, Spigel DR, Steins M, Ready NE, et al. Nivolumab versus Docetaxel in advanced nonsquamous non-small-cell lung Cancer. N Engl J Med. 2015;373(17):1627–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. 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 (New York, NY). 2007;9(2):166–80.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  22. Mizuno H, Kitada K, Nakai K, Sarai A. PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med Genet. 2009;2:18.

    Google Scholar 

  23. Lánczky A, Nagy Á, Bottai G, Munkácsy G, Szabó A, Santarpia L, et al. miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients. Breast Cancer Res Treat. 2016;160(3):439–46.

    Article  PubMed  CAS  Google Scholar 

  24. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47(W1):W556–w60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Wang G, Zhong WC, Bi YH, Tao SY, Zhu H, Zhu HX, et al. The prognosis of Peroxiredoxin family in breast Cancer. Cancer Manag Res. 2019;11:9685–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(Web Server issue):W214–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Smith JJ, Deane NG, Wu F, Merchant NB, Zhang B, Jiang A, et al. Experimentally derived metastasis gene expression profile predicts recurrence and death in patients with colon cancer. Gastroenterology. 2010;138(3):958–68.

    Article  CAS  PubMed  Google Scholar 

  31. Bonome T, Levine DA, Shih J, Randonovich M, Pise-Masison CA, Bogomolniy F, et al. A gene signature predicting for survival in suboptimally debulked patients with ovarian cancer. Cancer Res. 2008;68(13):5478–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Bijlsma MF, Sadanandam A, Tan P, Vermeulen L. Molecular subtypes in cancers of the gastrointestinal tract. Nat Rev Gastroenterol Hepatol. 2017;14(6):333–42.

    Article  CAS  PubMed  Google Scholar 

  33. Leone RD, Powell JD. Metabolism of immune cells in cancer. Nat Rev Cancer. 2020;20(9):516–31.

  34. Solinas C, Pusole G, Demurtas L, Puzzoni M, Mascia R, Morgan G, et al. Tumor infiltrating lymphocytes in gastrointestinal tumors: controversies and future clinical implications. Crit Rev Oncol Hematol. 2017;110:106–16.

    Article  PubMed  Google Scholar 

  35. Zheng W, Jiang C, Li R. Integrin and gene network analysis reveals that ITGA5 and ITGB1 are prognostic in non-small-cell lung cancer. OncoTargets Therap. 2016;9:2317–27.

    Article  CAS  Google Scholar 

  36. Zhao X, Wu Y, Lv Z. miR-128 modulates hepatocellular carcinoma by inhibition of ITGA2 and ITGA5 expression. Am J Transl Res. 2015;7(9):1564–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Gong C, Yang Z, Wu F, Han L, Liu Y, Gong W. miR-17 inhibits ovarian cancer cell peritoneal metastasis by targeting ITGA5 and ITGB1. Oncol Rep. 2016;36(4):2177–83.

    Article  CAS  PubMed  Google Scholar 

  38. Chen D, Xie J, Fiskesund R, Dong W, Liang X, Lv J, et al. Chloroquine modulates antitumor immune response by resetting tumor-associated macrophages toward M1 phenotype. Nat Commun. 2018;9(1):873.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Xiang W, Shi R, Kang X, Zhang X, Chen P, Zhang L, et al. Monoacylglycerol lipase regulates cannabinoid receptor 2-dependent macrophage activation and cancer progression. Nat Commun. 2018;9(1):2574.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Chen Y, Song Y, Du W, Gong L, Chang H, Zou Z. Tumor-associated macrophages: an accomplice in solid tumor progression. J Biomed Sci. 2019;26(1):78.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Genard G, Wera AC, Huart C, Le Calve B, Penninckx S, Fattaccioli A, et al. Proton irradiation orchestrates macrophage reprogramming through NFκB signaling. Cell Death Dis. 2018;9(7):728.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Najafi M, Hashemi Goradel N, Farhood B, Salehi E, Nashtaei MS, Khanlarkhani N, et al. Macrophage polarity in cancer: a review. J Cell Biochem. 2019;120(3):2756–65.

    Article  CAS  PubMed  Google Scholar 

  43. Lu L, Xie R, Wei R, Cai C, Bi D, Yin D, et al. Integrin α5 subunit is required for the tumor supportive role of fibroblasts in colorectal adenocarcinoma and serves as a potential stroma prognostic marker. Mol Oncol. 2019;13(12):2697–714.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Xu M, Zhang S, Jia L, Wang S, Liu J, Ma X, et al. E-M, an engineered Endostatin with high ATPase activity, Inhibits the Recruitment and Alternative Activation of Macrophages in Non-small Cell Lung Cancer. Front Pharmacol. 2017;8:532.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Guha P, Cunetta M, Somasundar P, Espat NJ, Junghans RP, Katz SC. Frontline science: functionally impaired geriatric CAR-T cells rescued by increased α5β1 integrin expression. J Leukoc Biol. 2017;102(2):201–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Harper J, Sainson RC. Regulation of the anti-tumour immune response by cancer-associated fibroblasts. Semin Cancer Biol. 2014;25:69–77.

    Article  CAS  PubMed  Google Scholar 

  47. Jiang H, Hegde S, DeNardo DG. Tumor-associated fibrosis as a regulator of tumor immunity and response to immunotherapy. Cancer Immunol Immunotherap. 2017;66(8):1037–48.

    Article  CAS  Google Scholar 

  48. Schnittert J, Bansal R, Prakash J. Targeting pancreatic stellate cells in Cancer. Trends Cancer. 2019;5(2):128–42.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by the National Natural Science Foundation of China (No. 81572350 to A. M. X), and the funder contributed to the collection and analysis of data, carrying out the experiment, writing the manuscript. In addition, this study was also supported by the project funding for the Construction of National and Provincial Key Clinical Specialty (No. Z155080000004 to the Department of General Surgery, The Fourth Affiliated Hospital of Anhui Medical University), this funder bodies (A. M. X. and G.W) played a key role in study design, experimental operation, data analysis and drafting and revising the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

H.Z and G.W conceived and designed the experiments. H.Z, G.W and A.M.X extracted the data and performed the analysis, H.Z and G.W performed the experiments, H.Z, H.X.Z and A.M.X wrote the manuscript. H.Z and A.M.X revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Aman Xu.

Ethics declarations

Ethics approval and consent to participate

This work was approved by the Academic Committee of The First Affiliated Hospital of Anhui Medical University and was conducted following the Declaration of Helsinki (certification no. Quick-PJ 2020-11-14). All patients provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declared no conflicts of interest in this work.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1

: Table S1. The information of antibody for immunohistochemistry. Table S2. The information of antibody for western blot.

Additional file 2

: Figure S1. The prognostic value of ITGA5 expression in non-gastrointestinal cancers. (A-D) Correlation between ITGA5 expression and prognosis of lung cancer and breast cancer in Kaplan-Meier Plotter; (E-F) Correlation between ITGA5 expression and prognosis of ovarian cancer in PrognoScan. OS, overall survival; DFS, disease free survival. Figure S2. The full-length, original blots for CD163 (A), ITGA5 (B), GATA3 (C) and STAT6 (D) in 6 paired gastric cancer tissues (T) and adjacent normal tissues (N). (the order of samples from left to right are: marker, N1, T1, N2, T2, N3, T3, marker, N4, T4, N5, T5, N6, T6, marker).

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

Zhu, H., Wang, G., Zhu, H. et al. ITGA5 is a prognostic biomarker and correlated with immune infiltration in gastrointestinal tumors. BMC Cancer 21, 269 (2021). https://doi.org/10.1186/s12885-021-07996-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-021-07996-1

Keywords