Abstract
Background
The landscape of head and neck cancers has changed with improvements in standard therapy; however, it is necessary to exploit advanced genomic approaches to identify novel diagnostic and prognostic biomarkers for head and neck squamous cell carcinoma (HNSC). ITGA3, ITGA5, and ITGA6, members of the integrin family of proteins, play active roles in cytoskeletal organization and cell migration, proliferation, and survival. However, the expression patterns and prognostic values of ITGA3, ITGA5, and ITGA6 in head and neck squamous cell carcinoma remain unclear.
Material/Methods
Different expression patterns and prognostic values of ITGA3, ITGA5, and ITGA6 were analyzed in patients with HNSC using various databases, including ONCOMINE, GEPIA, TIMER, HPA, Kaplan-Meier Plotter, GEO, and TCGA.
Results
Expression levels of ITGA3, ITGA5, and ITGA6 were substantially increased in patients with HNSC. Additionally, higher expression levels of ITGA3, ITGA5, and ITGA6 were associated with worse overall survival in patients with HNSC, and higher levels of ITGA3 correlated with a worse relapse-free survival.
Conclusions
ITGA3, ITGA5, and ITGA6 are potential diagnostic and prognostic biomarkers for HNSC. In particular, IGTA5 might be used as a significant independent prognostic factor in this cancer.
MeSH Keywords: Head and Neck Neoplasms, Medical Oncology, Prognosis, Transcriptome
Background
Each year approximately 600 000 people worldwide are diagnosed with head and neck squamous cell carcinoma (HNSC), which has a 40% to 50% mortality rate [1]. HNSC occurs in the mucous membranes of the oropharynx, nasopharynx, larynx, oral cavity, and hypopharynx [2]. In addition, there has been no change in the 5-year survival of patients with HNSC for more than 20 years because of the occurrence of distant metastasis and local recurrence [3]. More than 60% of patients with HNSC are diagnosed at an advanced stage (phase III or IV) due to difficulties in early diagnosis, large tumors with major local invasion, and metastases in regional nodes. The landscape of head and neck cancers has changed with improvements in standard therapy, such as minimally invasive, organ-sparing surgical techniques, advances in radiotherapy, and curative multimodal approaches [4]. Although the development of programmed death 1 immune checkpoint inhibitors has an important role in HNSC treatment, these targeted molecular therapies do not considerably benefit patients with HNSC [4,5]. Thus, it is necessary to exploit advanced genomic approaches to find novel diagnostic and prognostic biomarkers for HNSC.
Integrins are heterodimeric transmembrane receptors consisting of α and β subunits, constituting a well-established family of approximately 24 α and β subtypes [6]. As cell adhesion receptors, integrins mediate cell-to-cell and cell-to-extracellular matrix interactions; however, they also play significant roles in cytoskeletal organization and cell migration, proliferation, and survival through their roles in signal transduction [7]. ITGA3 functions as a surface adhesion molecule on the cell membrane [8], contacting extracellular matrix proteins of the laminin family, which are associated with tumors [9]. ITGA5, together with β1, forms a receptor for extracellular fibronectin, which is active in malignant tumor cells and tumor vasculature [6]. ITGA5 has been reported to act as an oncogene in hepatocellular carcinoma [10] and ovarian cancer [11]. Moreover, ITGA6 encodes integrin subunit α6, also known as CD49f, vla-6, or ITGA6B [12], and expression of ITGA6 in breast cancer is considerably decreased after the knockdown of thrombospondin-1, inducing a reduction in adhesion between cancer cells and the extracellular matrix and decreasing cancer angiogenesis [13]. Thus, ITGA3, ITGA5, and ITGA6 are strong biomarker candidates and therapeutic targets for HNSC.
Microarray analysis of gene expression is important in genetic and biomedical research [14]. In this study, we assessed data relating to HNSC using various well-constructed databases and platforms to examine HNSC expression and prognostic analyses based on ITGA3, ITGA5, and ITGA6.
Material and Methods
This research was carried out to examine the expression patterns and prognostic values of ITGA3, ITGA5, and ITGA6 in HNSC using online databases, platforms, and datasets. The study was conducted according to the principles expressed in the Declaration of Helsinki. All the datasets were collected from published literature, and written informed consent was confirmed.
ONCOMINE analysis
ONCOMINE was used as a test for ITGA3, ITGA5, and ITGA6 transcription levels in patients with HNSC [15,16]. The mRNA expression levels of ITGA3, ITGA5, and ITGA6 were compared between patients with different types of cancer and normal controls [17].
Gene expression omnibus analysis
To confirm the expression profiles of ITGA3, ITGA5, and ITGA6 in HNSC, 6 microarray series, GSE2379, GSE6631, GSE29330, GSE53819, GSE58911, and GSE107591, containing tumor and nontumor samples were collected from the Gene Expression Omnibus (GEO) database.
Gene expression profiling interactive analysis
The Gene Expression Profiling Interactive Analysis (GEPIA) database is an interactive service used between The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression data system to examine mRNA levels of ITGA3, ITGA5, and ITGA6 [18]. GEPIA provided associations of ITGA3, ITGA5, and ITGA6 in HNSC.
Human protein atlas analysis
The Human Protein Atlas (HPA) was used to compare the protein expressions of ITGA3, ITGA5, and ITGA6 in patients with different kinds of human cancer with those of normal controls [19].
Tumor IMmune Estimation Resource Analysis (TIMER)
The Tumor IMmune Estimation Resource (TIMER) provides a detailed analysis of immune infiltration for specific cancers [20]. The TIMER database includes 32 cancers and 10 897 samples from TCGA [21]. TIMER was used to evaluate the expression levels of ITGA3, ITGA5, and ITGA6 in different cancers and the correlation of ITGA3, ITGA5, and ITGA6 expression with levels of immune infiltrates [22]. The somatic copy number alterations (SCNA) module in the TIMER database can relate the rate of tumor infiltration of a particular gene to various changes in somatic copy number. SCNA incorporate high amplification (2), arm-level gain (1), normal diploid (0), arm-level deletion (−1), and deep deletion (−2). For each SCNA group, the degree of infiltration was compared to the average using a 2-sided Wilcoxon test.
Functional enrichment analysis
METASCAPE was applied for pathway and process enrichment analysis of ITGA3, ITGA5, and ITGA6 and adjacent genes significantly related to ITGA3, ITGA5, and ITGA6 alterations. The online METASCAPE resource enriched Gene Ontology (GO) terminology for biological process, cellular component, and molecular function categories as well as the Kyoto Genes and Genome Encyclopedia (KEGG). The molecular complex detection (MCODE) algorithm was used to classify strongly connected network components [23].
Kaplan-Meier plotter analysis
The Kaplan-Meier plotter is an interactive service containing microarray gene expression information and survival data from TCGA, GEO, and Cancer Biomedical Informatics Grid [23]. The OS and relapse-free survival (RFS) of patients with HNSC with high and low expressions of ITGA3, ITGA5, and ITGA6 were investigated [24].
TCGA analysis
Level 3 TCGA-HNSC data points were obtained from the UCSC Xena platform [25] and RTCGA package. The RNA-seq gene expression data for HNSC involved 500 cancer samples, of which 483 patients had complete survival data. Univariate and multivariate Cox regression analyses of clinicopathological characteristics, including age, sex, clinical stage, histological grade, and OS time, were performed [17].
Statistical analysis
The mRNA expression of ITGA3, ITGA5, and ITGA6 between HNSC tissues and normal controls was compared using the t test. R software (3.6.1) (https:/www.r-Project.org/) was used for statistical analyses. Data visualization was carried out using the “ggstatsplot” package in R (https://CRAN.R-project.org/package=ggstatsplot). For TCGA-HNSC raw data, Kaplan-Meier curves of OS were drawn by selecting the median expression levels of ITGA3, ITGA5, and ITGA6 as the cut-offs. Significant differences were examined by the log-rank test. Univariate and multivariate survival analyses were conducted using the Cox regression model, and risk factors with P<0.05 according to univariate analysis were selected for multivariate analysis [17].
Results
Expression levels of ITGA3, ITGA5, and ITGA6 in patients with HNSC
We utilized the ONCOMINE database to analyze the mRNA expression levels of ITGA3, ITGA5, and ITGA6 among patients with HNSC. In 13 HNSC datasets, the expression levels of ITGA3, ITGA5, and ITGA6 in most tumor tissues were higher than those in normal tissue (Table 1). For example, in Ginos’s dataset [26], ITGA3 was overexpressed in HNSC with a fold change of 2.207. In Peng’s dataset [27], ITGA3 was overexpressed in oral cavity squamous cell carcinoma with a fold change of 3.147. In Pyeon’s dataset [28], ITGA3 was overexpressed in tongue carcinoma with a fold change of 2.299. According to Talbot’s dataset [29], ITGA3 was overexpressed in tongue squamous cell carcinoma with a fold change of 2.068. In Estilo’s dataset [30], ITGA3 was overexpressed in tongue squamous cell carcinoma with a fold change of 3.789, and ITGA5 with a fold change of 3.027. ITGA5 overexpression in oral cavity squamous cell carcinoma with a fold change of 2.761 in Peng’s dataset [27] and in head and neck squamous cell carcinoma with a fold change of 2.307 in Ginos’s dataset [26] was also detected. For ITGA6 mRNA expression, Pyeon [28] reported a fold change of 2.757 in patients with tongue carcinoma and a fold change of 3.133 in patients with oropharyngeal carcinoma. In Estilo’s dataset [30], ITGA6 was overexpressed in tongue squamous cell carcinoma with a fold change of 6.290. According to Talbot’s dataset [29], ITGA6 was overexpressed in tongue squamous cell carcinoma with a fold change of 4.086. In Ye’s dataset [31], ITGA6 was overexpressed in tongue squamous cell carcinoma with a fold change of 2.958 (Table 1). For validation, we performed a meta-analysis of ITGA3, ITGA5, and ITGA6 expression with thresholds of P value ≤1E-4, fold change ≥2, and the top 10% gene rank in the ONCOMINE database. As illustrated in Figure 1, ITGA3, ITGA5, and ITGA6 were significantly upregulated in HNSC compared with normal tissues (P<0.0001).
Table 1.
Gene ID | Types of HNSC vs. normal | Fold change | P value | t-Test | References |
---|---|---|---|---|---|
ITGA3 | Tongue squamous cell carcinoma vs. normal | 3.789 | 2.94E-12 | 8.728 | Estilo et al., 2009 |
Head and neck squamous cell carcinoma vs. normal | 2.207 | 3.61E-8 | 6.437 | Ginos et al., 2004 | |
Oral cavity squamous cell carcinoma vs. normal | 3.147 | 3.88E-24 | 14.600 | Peng et al., 2011 | |
Tongue carcinoma vs. normal | 2.299 | 1.01E-5 | 5.745 | Pyeon et al., 2007 | |
Tongue squamous cell carcinoma vs. normal | 2.068 | 8.82E-10 | 7.310 | Talbot et al., 2005 | |
ITGA5 | Tongue squamous cell carcinoma vs. normal | 3.027 | 9.33E-7 | 5.339 | Estilo et al., 2009 |
Head and neck squamous cell carcinoma vs. normal | 2.307 | 1.37E-7 | 6.676 | Ginos et al., 2004 | |
Oral cavity squamous cell carcinoma vs. normal | 2.761 | 3.26E-18 | 11.351 | Peng et al., 2011 | |
ITGA6 | Tongue squamous cell carcinoma vs. normal | 6.290 | 2.30E-11 | 8.313 | Estilo et al., 2009 |
Oropharyngeal carcinoma vs. normal | 2.757 | 1.94E-5 | 5.944 | Pyeon et al., 2007 | |
Tongue carcinoma vs. normal | 3.133 | 1.81E-7 | 6.469 | Pyeon et al., 2007 | |
Tongue squamous cell carcinoma vs. normal | 4.086 | 5.58E-12 | 9.438 | Talbot et al., 2005 | |
Tongue squamous cell carcinoma vs. normal | 2.958 | 1.67E-7 | 6.534 | Ye et al., 2008 |
The results of GEPIA analysis indicated that the mRNA levels of ITGA3, ITGA5, and ITGA6 in HNSC were considerably greater than those in normal tissues (Figure 2, P<1E-4).
Moreover, mRNA expression levels of ITGA3, ITGA5, and ITGA6 were analyzed in HNSC and normal tissues using the GEO series (Table 2). The GSE2379, GSE6631, GSE29330, GSE53819, GSE58911, and GSE107591 results showed substantial upregulation in tumor tissues compared with normal tissues (all P<0.05, Figures 3, 4).
Table 2.
GEO series | Contributor(s) | Tumor | Nontumor | Platform |
---|---|---|---|---|
GSE2379 | Cromer A et al., 2004 | 14 | 4 | Affymetrix Human Genome U95A Array |
GSE6631 | Kuriakose MA et al., 2004 | 22 | 22 | Affymetrix Human Genome U95 Version 2 Array |
GSE29330 | Demokan S et al., 2013 | 13 | 5 | Affymetrix Human Genome U133 Plus 2.0 Array |
GSE53819 | Bao YN et al., 2014 | 18 | 18 | Agilent-014850 Whole Human Genome Microarray 4x44K G4112F |
GSE58911 | Lobert S. 2014 | 15 | 15 | Affymetrix Human Gene 1.0 ST Array [transcript (gene) version] |
GSE107591 | Blandino G et al., 2017 | 24 | 23 | Affymetrix Human Gene 1.0 ST Array [transcript (gene) version] |
To further evaluate ITGA3, ITGA5, and ITGA6 expression in patients with HNSC, RNA-seq data for multiple malignancies in TCGA were determined using ONCOMINE and TIMER, with results showing significantly higher expression levels (Supplementary Figure 1, P<0.0001; Supplementary Figure 2, P<0.001).
We also analyzed immunohistochemistry images using HPA to explore the expression of ITGA3, ITGA5, and ITGA6 proteins in HNSC. Compared with normal tissue, HNSC cancer tissues exhibited higher protein expression of all these integrins, as shown in Figure 5.
ITGA3, ITGA5, ITGA6 and neighbor gene network in patients with HNSC
GeneMANIA was employed to perform gene-level correlation analysis of ITGA3, ITGA5, and ITGA6 and neighboring genes (Figure 6A) revealing SPP1, ITGB1, LAMC3, CD151, CD9, RPSA, ITGB3, COL17A1, PLEC, COL18A1, ITGB4, HOXD3, COL4A3, ANGPTL3, THBS2, FHL2, ITGA2, LAMB3, LAMA3, and PMP22 to be closely associated. Moreover, STRING analysis was conducted to identify interactions between ITGA3, ITGA5, and ITGA6 and neighboring genes at the level of protein expression. In this analysis, ITGA5 showed interaction with ITGA6 in coexpression, text mining, and protein homology, and they were both closely associated with CD97, CD151, PLEC, ITGB4, TNN, ITGA2B, and ITGAV (Figure 6B). We evaluated the transcripts per kilobase million association among ITGA3, ITGA5, and ITGA6 in HNSC based on Pearson correlation using GEPIA data (Figure 7) and found positive correlations between ITGA3 and ITGA5 (r=0.62, P<0.001), ITGA3 and ITGA6 (r=0.65, P<0.001), and ITGA5 and ITGA6 (r=0.47, P<0.001).
Expression of ITGA3, ITGA5, and ITGA6 in HNSC is linked to the level of immune infiltration
Sentinel lymph node status and survival are independent predictors of tumor lymph invasion [32], and the tumor purity of clinical samples is a significant factor when evaluating immune infiltration through genomic approaches [33]. In the present study, we investigated associations between the expression of ITGA3, ITGA5, and ITGA6 and immune cell populations using transcriptomic data for various molecular subtypes and the TCGA-HNSC cohort in TIMER (Figure 8). There were statistically significant negative correlations between the expression level of ITGA3 and tumor purity (r=−0.164, P=2.51e-04), B cells (r=−0.226, P=6.60e-07), CD8+ T cells (r=−0.174, P=1.41e-04), and macrophages (r=−0.067, P=1.42e-01) and a statistically significant positive correlations with infiltrating levels of CD4+ T cells (r=0.108, P=1.83e-02) and neutrophils (r=0.224, P=7.41e-07). Furthermore, there was no significant correlation observed between dendritic cells and ITGA3 expression (r=0.098, P=3.13e-02). In addition, the expression level of ITGA5 was significantly negatively correlated with B cells (r=−0.102, P=2.59e-02) and CD8+ T cells (r=−0.102, P=2.59e-02) and significantly positively correlated with infiltrating levels of CD4+ T cells (r=0.208, P=4.10e-06), macrophages (r=0.265, P=3.23e-09), neutrophils (r=0.180, P=7.27e-05), and dendritic cells (r=0.233, P=2.17e-07). However, ITGA5 expression did not have a significant correlation with tumor purity (r=−0.08, P=7.54e-02), whereas ITGA6 expression was significantly negatively correlated with B cells (r=−0.214, P=2.60e-06) and CD8+ T cells (r=−0.253, P=2.36e-08) and significantly positively correlated with infiltrating levels of CD4+ T cells (r=0.117, P=1.01e-02). Nevertheless, no significant correlations between ITGA6 and tumor purity (r=−0.087, P=5.31e-02), macrophages (r=−0.091, P=4.51e-02), neutrophils (r=0.091, P=4.77e-02) or dendritic cells (r=0.015, P=7.45e-01) were found.
Intriguingly, we found that ITGA3 expression levels had a statistically significant correlation with arm-level gain in CD8+ T cells (P=8.69e-05), neutrophils (P=4.92e-06), and dendritic cells (P=0.009). ITGA5 levels also had a statistically significant correlation with arm-level deletion in CD8+ T cells (P=6.52e-05), neutrophils (P=0.003) and dendritic cells (P=0.017), and arm-level gain in B cells (P=0.012). Moreover, ITGA6 expression correlated significantly with arm-level gain in CD8+ T cells (P=5.19e-13), CD4+ T cells (P<0.001), B cells (P=6.89e-05), macrophages (P=0.001), neutrophils (P=1.47e-11), and dendritic cells (P=8.26e-06) (Figure 9).
Functional enrichment analysis of ITGA3, ITGA5, and ITGA6 in patients with HNSC
GO and KEGG analyses in Metascape were carried out for ITGA3, ITGA5, and ITGA6 and neighboring genes. The top 21 GO enrichment items were detected as follows: 9 items in biological process, 6 items in molecular function, and 6 items in cellular component (Figure 10A, 10B and Table 3). ITGA3, ITGA5, and ITGA6 and related genes showed strong enrichment in the biological process categories hemidesmosome assembly, extracellular matrix organization, integrin-mediated signaling pathway, skin morphogenesis, angiogenesis, cell morphogenesis involved in differentiation, receptor internalization, regulation of endocytosis, and positive regulation of the apoptotic process. Molecular function enrichment was mainly found for transcriptional regulation by integrin binding, extracellular matrix binding, neuregulin binding, extracellular matrix structural constituent, sulfur compound binding, and cadherin binding. Among terms in cellular components, the basement membrane, integrin complex, hemidesmosome, collagen trimer, platelet alpha granule, and perinuclear region of cytoplasm were highlighted. The top 4 KEGG pathways for ITGA3, ITGA5, and ITGA6 and the related genes are presented in Figure 10C and 10D and Table 4, including extracellular matrix-receptor interaction, small cell lung cancer, arrhythmogenic right ventricular cardiomyopathy, and protein digestion and absorption.
Table 3.
GO | Category | Description | Count | % | Log10(P) | Log10(q) |
---|---|---|---|---|---|---|
GO:0031581 | GO Biological Processes | Hemidesmosome assembly | 7 | 58.33 | −21.18 | −16.98 |
GO:0030198 | GO Biological Processes | Extracellular matrix organization | 14 | 66.67 | −20.44 | −16.54 |
GO:0007229 | GO Biological Processes | Integrin-mediated signaling pathway | 8 | 66.67 | −16.26 | −13.06 |
GO:0043589 | GO Biological Processes | Skin morphogenesis | 3 | 25.00 | −7.81 | −5.43 |
GO:0001525 | GO Biological Processes | Angiogenesis | 7 | 33.33 | −6.34 | −4.11 |
GO:0000904 | GO Biological Processes | Cell morphogenesis involved in differentiation | 6 | 28.57 | −4.48 | −2.38 |
GO:0031623 | GO Biological Processes | Receptor internalization | 3 | 14.29 | −3.89 | −1.85 |
GO:0030100 | GO Biological Processes | Regulation of endocytosis | 3 | 25.00 | −3.48 | −1.46 |
GO:0043065 | GO Biological Processes | Positive regulation of apoptotic process | 4 | 19.05 | −2.60 | −0.69 |
GO:0005604 | GO Cellular Components | Basement membrane | 9 | 52.94 | −17.36 | −14.07 |
GO:0008305 | GO Cellular Components | Integrin complex | 7 | 41.18 | −16.27 | −13.29 |
GO:0030056 | GO Cellular Components | Hemidesmosome | 4 | 23.53 | −10.93 | −8.90 |
GO:0005581 | GO Cellular Components | Collagen trimer | 3 | 17.65 | −4.53 | −2.82 |
GO:0031091 | GO Cellular Components | Platelet alpha granule | 3 | 17.65 | −4.47 | −2.79 |
GO:0048471 | GO Cellular Components | Perinuclear region of cytoplasm | 5 | 21.74 | −3.28 | −1.73 |
GO:0005178 | GO molecular functions | Integrin binding | 13 | 65.00 | −24.58 | −20.92 |
GO:0050840 | GO molecular functions | Extracellular matrix binding | 7 | 58.33 | −15.78 | −12.81 |
GO:0038132 | GO molecular functions | Neuregulin binding | 3 | 25.00 | −9.03 | −6.54 |
GO:0005201 | GO molecular functions | Extracellular matrix structural constituent | 6 | 30.00 | −8.48 | −6.02 |
GO:1901681 | GO molecular functions | Sulfur compound binding | 3 | 15.00 | −2.94 | −0.99 |
GO:0045296 | GO molecular functions | Cadherin binding | 3 | 15.00 | −2.61 | −0.68 |
Table 4.
GO | Category | Description | Count | % | Log10(P) | Log10(q) |
---|---|---|---|---|---|---|
hsa04512 | KEGG Pathway | ECM-receptor interaction | 12 | 80.00 | −27.34 | −24.84 |
hsa05222 | KEGG Pathway | Small cell lung cancer | 8 | 38.10 | −14.53 | −12.83 |
hsa05412 | KEGG Pathway | Arrhythmogenic right ventricular cardiomyopathy (ARVC) | 7 | 46.67 | −14.01 | −12.40 |
hsa04974 | KEGG Pathway | Protein digestion and absorption | 3 | 20.00 | −4.66 | −3.57 |
Overall, the protein-protein enrichment analysis in Metascape demonstrated the involvement of ITGA3, ITGA5, and ITGA6 in HNSC. Figure 11A and 11B illustrate the network of protein-protein interactions and MCODE components in the gene lists, whereby the most critical components of MCODE are indicated by the interaction network among proteins. The results according to pathway and process enrichment analyses applied separately in MCODE showed that the biological roles of these proteins are mainly related to type I hemidesmosome assembly, and assembly of collagen fibrils and other multimeric structures.
Association of increased mRNA expression of ITGA3, ITGA5, and ITGA6 with improved prognosis in HNSC
We further studied the relevance of ITGA3, ITGA5, and ITGA6 to the survival of patients with HNSC. To explore the association between mRNA levels of ITGA3, ITGA5, and ITGA6 and survival, Kaplan-Meier Plotter tools were applied to support publicly available data sets. Increasing ITGA3, ITGA5, and ITGA6 mRNA levels were significantly linked to OS (P<0.05) (Figure 12) among all patients with HNSC in the Kaplan-Meier curve and log-rank test analyses. For example, survival curves showed that high ITGA3 expression was related to poor RFS (P=0.017), whereas high levels of ITGA5 and ITGA6 mRNA did not appear to be (ITGA3-OS hazard ratio (HR)=1.73, 95% CI=1.30–2.29, P=0.00014; ITGA3-RFS HR=2.40, 95% CI=1.14–5.06, P=0.017; ITGA5-OS HR=1.63, 95% CI=1.25–2.31, P=0.00029; ITGA5-RFS HR=1.95, 95% CI=0.91–4.17, P=0.079; ITGA6-OS HR=1.70, 95% CI=1.30–2.23, P=1e-04; ITGA6-RFS HR=0.74, 95% CI=0.35–1.56, P=0.42). We subsequently performed Cox regression analysis to calculate the prognostic values of ITGA3, ITGA5, and ITGA6 in HNSC based on original TCGA data, which revealed that the high expression of these genes, age >65, and female sex to be associated with worse OS in HNSC (Table 5). In addition, multivariate Cox analysis showed high ITGA5 expression to be on the threshold of significance (HR=1.347; 95% CI [0.988–1.837]; P=0.059), and multivariate Cox analysis adjusted for high ITGA5 expression, age, and sex confirmed high ITGA5 expression as an independent prognostic biomarker in patients with HNSC (HR=1.466; 95% CI [1.108–1.940]; P=0.007; Supplementary Table 1). Conversely, there was a lack of significant findings for other genes with regard to OS in HNSC.
Table 5.
Characteristics | Univariate analysis | Multivariate analysis | ||||
---|---|---|---|---|---|---|
HR | X95.CI | p | HR | X95.CI | P | |
Age: >65 vs. ≤65 | 1.337 | 1.009–1.771 | 0.043 | 1.315 | 0.984–1.758 | 0.064 |
Gender: Male vs. Female | 0.720 | 0.537–0.965 | 0.028 | 0.788 | 0.582–1.067 | 0.123 |
Clinical stage: III/IV vs. I/II | 1.139 | 0.814–1.592 | 0.447 | |||
Histological grade: G3/G4 vs. G1/G2 | 0.931 | 0.685–1.264 | 0.646 | |||
ITGA3 expression: High vs. low | 1.332 | 1.010–1.758 | 0.042 | 1.046 | 0.757–1.446 | 0.784 |
ITGA5 expression: High vs. low | 1.464 | 1.108–1.935 | 0.007 | 1.347 | 0.988–1.837 | 0.059 |
ITGA6 expression: High vs. low | 1.353 | 1.026–1.784 | 0.032 | 1.202 | 0.871–1.658 | 0.264 |
Discussion
ITGA3, ITGA5, and ITGA6 deregulation has been reported in numerous types of cancers. Although the effects of ITGA3, ITGA5, and ITGA6 on the tumorigenesis and prognosis of different cancers have been partially corroborated, advanced bioinformatics analysis of HNSC has yet to be conducted [34]. In fact, the present study is the first to explore the expression and prognostic value on OS and RFS of ITGA3, ITGA5, and ITGA6 in patients with HNSC. We carried out a comprehensive study of the patterns of expression and prognostic value of ITGA3, ITGA5, and ITGA6 in HNSC. We are confident that our findings will provide accessible knowledge that will contribute to an improvement in treatment modeling and prognostic accuracy for HNSC.
Previous research has found that ITGA3 can facilitate cancer development by activating the PI3K-Akt signaling pathway, which promotes the proliferation, migration, and invasion in multiple cancer types [35]. Furthermore, the value of ITGA3 expression during relapse and metastasis in colorectal cancer [36,37], lung cancer [38], prostate cancer [39], and gastric cancer [40] has been determined. ITGA3 can serve as a complementary marker for endometrial cancer [41] and bladder cancer [42]. In the present study, ONCOMINE, GEO, GEPIA,and TIMER data indicated that the mRNA expression level of ITGA3 was markedly higher in patients with HNSC than in normal controls. We also searched the HPA data to validate ITGA3 protein expression, which revealed strongly increased expression in HNSC. Another noteworthy finding was that ITGA3 expression correlated with high levels of immune infiltration in various HNSC. Overall, our results demonstrated that ITGA3 is an important factor for infiltrating immune cells in HNSC and may act as a prognostic biomarker for patients with HNSC. Using the Kaplan-Meier plotter, we observed that upregulated IGTA3 mRNA expression levels are related to worse OS and RFS in all patients with HNSC, indicating that ITGA3 represents an oncogene. Furthermore, we evaluated the connection between the expression level of ITGA3 and OS in HNSC using univariate and multivariate Cox regression analyses of TCGA data, and the results showed a prognostic role for the level of ITGA3 expression.
The knockdown of ITGA5 decreases cell adhesion and promotes apoptosis [43], whereas its upregulation enhances cell adhesion, invasion, and epithelial-mesenchymal transition in colorectal cancer cells [44]. Based on the ONCOMINE, GEO, GEPIA, and TIMER databases, we found that the expression level of ITGA5 was significantly increased in HNSC. Additionally, we examined the protein expression profile of ITGA5 in HPA, which showed levels similar to those of mRNA expression. We speculate that ITGA5 expression is related to different levels of immune infiltration in HNSC; it is also possible that ITGA5 may serve as a prognostic biomarker for HNSC infiltration. Upregulation of ITGA5 mRNA expression correlated with a worse OS but not RFS in patients with HNSC. According to univariate and multivariate Cox regression analyses of TCGA data, ITGA5 mRNA expression might be a significant prognostic factor for OS.
ITGA6 is overexpressed in multiple tumors, promoting tumorigenesis and metastasis [45,46], and expression of ITGA6 is strongly related to the occurrence of intravesical recurrence [47]. Several studies have shown that ITGA6 is a presumptive stem marker [48,49]. Using data from ONCOMINE and GEPIA, the expression level of ITGA6 was significantly upregulated in our study, as confirmed by GEO data and TIMER datasets. Furthermore, we examined the ITGA6 expression profile in HPA, which showed results similar to those for mRNA expression, and ITGA6 can also act as a prognostic biomarker because of its relationship with infiltration of various immune cells in HNSC as revealed by TCGA and TIMER datasets. In addition, we explored the increase in ITGA6 mRNA expression, which correlated significantly with worse OS using the Kaplan-Meier plotter but was not related to RFS. Analysis of ITGA6 expression in TCGA data using univariate and multivariate Cox regression revealed its prognostic value for OS.
A network of ITGA3, ITGA5, and ITGA6 and 20 neighboring genes was generated, which showed mostly enrichment in cancer-related pathways associated with the evolution of multiplex cancers. Our study adds to the growing evidence of the complexity of ITGA3, ITGA5, and ITGA6 and their associated signaling pathways, which provide insights into the rational development of targeted therapy.
Our research has the following limitations that need to be addressed in the future. First, we should explore the mechanisms of ITGA3, ITGA5, and ITGA6 for their development as biomarkers and for prognosis in further experiments. Second, we should validate our research through traditional experiments.
Conclusions
The expression patterns and prognostic values of ITGA3, ITGA5, and ITGA6 were comprehensively studied in patients with HNSC by carrying out a bioinformatics study using a variety of online platforms and data sets. ITGA3, ITGA5, and ITGA6 may serve as diagnostic and prognostic biomarkers for HNSC. In particular, IGTA5 might be used as a significant independent prognostic factor in HNSC. We hope that our results will enrich the diagnostic and therapeutic knowledge of HNSC.
Supplementary Data
Supplementary Table 1.
Characteristics | Univariate analysis | Multivariate analysis | ||||
---|---|---|---|---|---|---|
HR | X95.CI | p | HR | X95.CI | P | |
Age: >65 vs. ≤65 | 1.337 | 1.009–1.771 | 0.043 | 1.299 | 0.974–1.734 | 0.075 |
Gender: Male vs. Female | 0.720 | 0.537–0.965 | 0.028 | 0.78 | 0.578–1.053 | 0.105 |
Clinical stage: III/IV vs. I/II | 1.139 | 0.814–1.592 | 0.447 | |||
Histological grade: G3/G4 vs. G1/G2 | 0.931 | 0.685–1.264 | 0.646 | |||
ITGA5 expression: High vs. low | 1.464 | 1.108–1.935 | 0.007 | 1.466 | 1.108–1.94 | 0.007 |
Acknowledgements
The authors thank Wu Pancheng for his technical support in this research.
Footnotes
Conflict of interest
None.
Source of support: This work was supported by the National Natural Science Foundation of Shandong Provincial (no. ZR2018MH017) and the Key R&D Program Fund Project of Shandong Provincial (no. 2018GSF118001)
References
- 1.Leemans CR, Snijders PJF, Brakenhoff RH. The molecular landscape of head and neck cancer. Nat Rev Cancer. 2018;18(5):269–82. doi: 10.1038/nrc.2018.11. [DOI] [PubMed] [Google Scholar]
- 2.Koshizuka K, Hanazawa T, Kikkawa N, et al. Regulation of ITGA3 by the anti-tumor miR-199 family inhibits cancer cell migration and invasion in head and neck cancer. Cancer Sci. 2017;108(8):1681–92. doi: 10.1111/cas.13298. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 3.Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2016. Cancer J Clin. 2016;66(1):7–30. doi: 10.3322/caac.21332. [DOI] [PubMed] [Google Scholar]
- 4.Chow LQM. Head and neck cancer. New Engl J Med. 2020;382(1):60–72. doi: 10.1056/NEJMra1715715. [DOI] [PubMed] [Google Scholar]
- 5.Bonner JA, Harari PM, Giralt J. Radiotherapy plus cetuximab for locoregionally advanced head and neck cancer: 5-year survival data from a phase 3 randomised trial, and relation between cetuximab-induced rash and survival. Lancet Oncol. 2010;11(1):14. doi: 10.1016/S1470-2045(09)70311-0. [DOI] [PubMed] [Google Scholar]
- 6.Hynes RO. Integrins: Bidirectional, allosteric signaling machines. Cell. 2002;110(6):673–87. doi: 10.1016/s0092-8674(02)00971-6. [DOI] [PubMed] [Google Scholar]
- 7.Schnittert J, Bansal R, Storm G, Prakash J. Integrins in wound healing, fibrosis and tumor stroma: High potential targets for therapeutics and drug delivery. Adv Drug Deliver Rev. 2018;129:37–53. doi: 10.1016/j.addr.2018.01.020. [DOI] [PubMed] [Google Scholar]
- 8.Jiao Y, Li YQ, Liu SY, et al. ITGA3 serves as a diagnostic and prognostic biomarker for pancreatic cancer. Onco Targets Ther. 2019;12:4141–52. doi: 10.2147/OTT.S201675. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 9.Kawataki T, Yamane T, Naganuma H, et al. Laminin isoforms and their integrin receptors in glioma cell migration and invasiveness: Evidence for a role of alpha 5-laminin(s) and alpha 3 beta 1 integrin. Exp Cell Res. 2007;313(18):3819–31. doi: 10.1016/j.yexcr.2007.07.038. [DOI] [PubMed] [Google Scholar]
- 10.Zhao XH, Wu YF, Lv ZP. miR-128 modulates hepatocellular carcinoma by inhibition of ITGA2 and ITGA5 expression. Am J Transl Res. 2015;7(9):1564–73. [Retracted Article; (Am J Transl Res, 2016; 8(4): 1919] [PMC free article] [PubMed] [Google Scholar] [Retracted]
- 11.Gong C, Yang ZY, Wu FH, et al. miR-17 inhibits ovarian cancer cell peritoneal metastasis by targeting ITGA5 and ITGB1. Oncol Rep. 2016;36(4):2177–83. doi: 10.3892/or.2016.4985. [DOI] [PubMed] [Google Scholar]
- 12.Jin H, Ying XL, Que B, et al. N-6-methyladenosine modification of ITGA6 mRNA promotes the development and progression of bladder cancer. Ebiomedicine. 2019;47:195–207. doi: 10.1016/j.ebiom.2019.07.068. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 13.John AS, Rothman VL, Tuszynski GP. Thrombospondin-1 (TSP-1) simulates expression of integrin alpha6 in human breast carcinoma cells: A downstream modulator of TSP-1-induced cellular adhesion. J Oncol. 2010;2010 doi: 10.1155/2010/645376. 645376. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 14.Raghavachari N, Barb J, Yang YQ, et al. A systematic comparison and evaluation of high density exon arrays and RNA-seq technology used to unravel the peripheral blood transcriptome of sickle cell disease. BMC Med Genomics. 2012;5:28. doi: 10.1186/1755-8794-5-28. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 15.Rhodes DR, Yu J, Shanker K, et al. ONCOMINE: A cancer microarray database and integrated data-mining platform. Neoplasia. 2004;6(1):1–6. doi: 10.1016/s1476-5586(04)80047-2. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 16.Rhodes DR, Kalyana-Sundaram S, Mahavisno V, et al. Oncomine 3.0: Genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9(2):166–80. doi: 10.1593/neo.07112. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 17.Wu PC, Wang YY, Wu YJ, et al. Expression and prognostic analyses of ITGA11, ITGB4 and ITGB8 in human non-small cell lung cancer. Peer J. 2019;7:e8299. doi: 10.7717/peerj.8299. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 18.Tang ZF, Li CW, Kang BX, et al. GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102. doi: 10.1093/nar/gkx247. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 19.Uhlen M, Zhang C, Lee S, et al. A pathology atlas of the human cancer transcriptome. Science. 2017;357(6352):eaan2507. doi: 10.1126/science.aan2507. [DOI] [PubMed] [Google Scholar]
- 20.Li TW, Fan JY, Wang BB, et al. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):E108–10. doi: 10.1158/0008-5472.CAN-17-0307. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 21.Li B, Severson E, Pignon JC, et al. Comprehensive analyses of tumor immunity: Implications for cancer immunotherapy. Genome Biol. 2016;17(1):174. doi: 10.1186/s13059-016-1028-7. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 22.Pan JH, Zhou H, Cooper L, et al. LAYN is a prognostic biomarker and correlated with immune infiltrates in gastric and colon cancers. Front Immunol. 2019;10:6. doi: 10.3389/fimmu.2019.00006. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 23.Zhou Q, Zhang F, He Z, Zuo MZ. E2F2/5/8 serve as potential prognostic biomarkers and targets for human ovarian cancer. Front Oncol. 2019;9:161. doi: 10.3389/fonc.2019.00161. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 24.Gyorffy B, Lanczky A, Szallasi Z. Implementing an online tool for genome-wide validation of survival-associated biomarkers in ovarian-cancer using microarray data from 1287 patients. Endocr Relat Cancer. 2012;19(2):197–208. doi: 10.1530/ERC-11-0329. [DOI] [PubMed] [Google Scholar]
- 25.Goldman M, Craft B, Swatloski T, et al. The UCSC Cancer Genomics Browser: Update 2015. Nucleic Acids Res. 2015;43(D1):D812–17. doi: 10.1093/nar/gku1073. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 26.Ginos MA, Page GP, Michalowicz BS, et al. Identification of a gene expression signature associated with recurrent disease in squamous cell carcinoma of the head and neck. Cancer Res. 2004;64(1):55–63. doi: 10.1158/0008-5472.can-03-2144. [DOI] [PubMed] [Google Scholar]
- 27.Peng XC, Gong FM, Zhao YW, et al. Comparative proteomic approach identifies PKM2 and cofilin-1 as potential diagnostic, prognostic and therapeutic targets for pulmonary adenocarcinoma. PLoS One. 2011;6(11):e27309. doi: 10.1371/journal.pone.0027309. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 28.Pyeon D, Newton NA, Lambert PF, et al. Fundamental differences in cell cycle deregulation in human papillomavirus-positive and human papillomavirus-negative head/neck and cervical cancers. Cancer Res. 2007;67(10):4605–19. doi: 10.1158/0008-5472.CAN-06-3619. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 29.Talbot SG, Estilo C, Maghami E, et al. Gene expression profiling allows distinction between primary and metastatic squamous cell carcinomas in the lung. Cancer Res. 2005;65(8):3063–71. doi: 10.1158/0008-5472.CAN-04-1985. [DOI] [PubMed] [Google Scholar]
- 30.Estilo CL, Oc P, Talbot S, et al. Oral tongue cancer gene expression profiling: Identification of novel potential prognosticators by oligonucleotide microarray analysis. BMC Cancer. 2009;9:11. doi: 10.1186/1471-2407-9-11. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 31.Ye H, Yu TW, Temam S, et al. Transcriptomic dissection of tongue squamous cell carcinoma. BMC Genomics. 2008;9:69. doi: 10.1186/1471-2164-9-69. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 32.Azimi F, Scolyer RA, Rumcheva P, et al. Tumor-infiltrating lymphocyte grade is an independent predictor of sentinel lymph node status and survival in patients with cutaneous melanoma. J Clin Oncol. 2012;30(21):2678–83. doi: 10.1200/JCO.2011.37.8539. [DOI] [PubMed] [Google Scholar]
- 33.Yoshihara K, Shahmoradgoli M, Martinez E, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612. doi: 10.1038/ncomms3612. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 34.Son SY, Lee S, Singh D, et al. Comprehensive secondary metabolite profiling toward delineating the solid and submerged-state fermentation of Aspergillus oryzae KCCM 12698. Front Microbiol. 2018;9:1076. doi: 10.3389/fmicb.2018.01076. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 35.Jiao Y, Li Y, Liu S, et al. ITGA3 serves as a diagnostic and prognostic biomarker for pancreatic cancer. Onco Targets Ther. 2019;12:4141–52. doi: 10.2147/OTT.S201675. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 36.Bauer KM, Watts TN, Buechler S, Hummon AB. Proteomic and functional investigation of the colon cancer relapse-associated genes NOX4 and ITGA3. J Proteome Res. 2014;13(11):4910–18. doi: 10.1021/pr500557n. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 37.Linhares MM, Affonso RJ, Viana LD, et al. Genetic and immunohistochemical expression of integrins ITGAV, ITGA6, and ITGA3 as prognostic factor for colorectal cancer: Models for global and disease-free survival. PLoS One. 2015;10(12):e0144333. doi: 10.1371/journal.pone.0144333. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 38.Hsu YC, Chen HY, Yuan SS, et al. Genome-wide analysis of three-way interplay among gene expression, cancer cell invasion and anti-cancer compound sensitivity. BMC Med. 2013;11:106. doi: 10.1186/1741-7015-11-106. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 39.Kurozumi A, Goto Y, Matsushita R, et al. Tumor-suppressive microRNA-223 inhibits cancer cell migration and invasion by targeting ITGA3/ITGB1 signaling in prostate cancer. Cancer Sci. 2016;107(1):84–94. doi: 10.1111/cas.12842. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 40.Scartozzi M, Loretelli C, Bearzi I, et al. Allele polymorphisms of tumor integrins correlate with peritoneal carcinosis capability of gastric cancer cells in radically resected patients. Ann Oncol. 2011;22(4):897–902. doi: 10.1093/annonc/mdq542. [DOI] [PubMed] [Google Scholar]
- 41.Bateman NW, Dubil EA, Wang GS, et al. Race-specific molecular alterations correlate with differential outcomes for black and white endometrioid endometrial cancer patients. Cancer. 2017;123(20):4004–12. doi: 10.1002/cncr.30813. [DOI] [PubMed] [Google Scholar]
- 42.Sakaguchi T, Yoshino H, Yonemori M, et al. Regulation of ITGA3 by the dual-stranded microRNA-199 family as a potential prognostic marker in bladder cancer. Br J Cancer. 2017;116(8):1077–87. doi: 10.1038/bjc.2017.43. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 43.Murillo CA, Rychahou PG, Evers BM. Inhibition of alpha 5 integrin decreases PI3K activation and cell adhesion of human colon cancers. Surgery. 2004;136(2):143–49. doi: 10.1016/j.surg.2004.04.006. [DOI] [PubMed] [Google Scholar]
- 44.Nam EH, Lee Y, Moon B, et al. Twist1 and AP-1 cooperatively upregulate integrin alpha 5 expression to induce invasion and the epithelial-mesenchymal transition. Carcinogenesis. 2015;36(3):327–37. doi: 10.1093/carcin/bgv005. [DOI] [PubMed] [Google Scholar]
- 45.Brooks DLP, Schwab LP, Krutilina R, et al. ITGA6 is directly regulated by hypoxia-inducible factors and enriches for cancer stem cell activity and invasion in metastatic breast cancer models. Mol Cancer. 2016;15:26. doi: 10.1186/s12943-016-0510-x. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 46.Stewart RL, West D, Wang C, et al. Elevated integrin alpha 6 beta 4 expression is associated with venous invasion and decreased overall survival in non-small cell lung cancer. Hum Pathol. 2016;54:174–83. doi: 10.1016/j.humpath.2016.04.003. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 47.Behnsawy HM, Miyake H, Abdalla MA, et al. Expression of integrin proteins in non-muscle-invasive bladder cancer: Significance of intravesical recurrence after transurethral resection. BJU Int. 2011;107(2):240–46. doi: 10.1111/j.1464-410X.2010.09534.x. [DOI] [PubMed] [Google Scholar]
- 48.Bigoni-Ordonez GD, Czarnowski D, Parsons T, et al. Integrin alpha 6 (CD49f), the microenvironment and cancer stem cells. Curr Stem Cell Res Ther. 2019;14(5):428–36. doi: 10.2174/1574888X13666181002151330. [DOI] [PubMed] [Google Scholar]
- 49.Krebsbach PH, Villa-Diaz LG. The role of integrin alpha 6 (CD49f) in stem cells: More than a conserved biomarker. Stem Cells Dev. 2017;26(15):1090–99. doi: 10.1089/scd.2016.0319. [DOI] [PMC free article] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials
Supplementary Table 1.
Characteristics | Univariate analysis | Multivariate analysis | ||||
---|---|---|---|---|---|---|
HR | X95.CI | p | HR | X95.CI | P | |
Age: >65 vs. ≤65 | 1.337 | 1.009–1.771 | 0.043 | 1.299 | 0.974–1.734 | 0.075 |
Gender: Male vs. Female | 0.720 | 0.537–0.965 | 0.028 | 0.78 | 0.578–1.053 | 0.105 |
Clinical stage: III/IV vs. I/II | 1.139 | 0.814–1.592 | 0.447 | |||
Histological grade: G3/G4 vs. G1/G2 | 0.931 | 0.685–1.264 | 0.646 | |||
ITGA5 expression: High vs. low | 1.464 | 1.108–1.935 | 0.007 | 1.466 | 1.108–1.94 | 0.007 |