A Comprehensive Pan-Cancer Analysis Reveals GRB7 as a Potential Diagnostic and Prognostic Biomarker
Abstract
Background and aim: Growth factor receptor-bound protein 7 (GRB7) belongs to a group of adaptor proteins characterized by their conserved multidomain structure. These proteins are involved in cellular signaling pathways that regulate cell growth, proliferation, and differentiation. Alterations in GRB7 expression have been linked to multiple human cancers. However, its role as a diagnostic and prognostic marker remains underexplored. This study aimed to assess the diagnostic and prognostic relevance of GRB7 in a comprehensive pan-cancer analysis.
Materials and methods: GRB7 expression across different cancers was evaluated using the Tumor Immune Estimation Resource (TIMER), Gene Expression Profiling Interactive Analysis (GEPIA), and the University of Alabama at Birmingham Cancer Data Analysis Portal (UALCAN). The correlation of GRB7 expression with various clinicopathological parameters was assessed by the UALCAN database. Additionally, the Human Protein Atlas (HPA) (https://www.proteinatlas.org/) was used to illustrate the histology of kidney cancer tissues. The correlation between GRB7 expression and prognosis was explored using the Kaplan-Meier plotter, GEPIA, and UALCAN databases. The TIMER database was used to explore the connection between GRB7 expression in tumor tissues and the infiltration of immune cells. Moreover, genetic alterations of the GRB7 gene were detected by the cBioPortal database. Results were validated by the GEO2R database.
Results: GRB7 expression was significantly upregulated in bladder urothelial carcinoma (BLCA), cervical squamous cell carcinoma (CESC), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), rectum adenocarcinoma (READ), thyroid carcinoma (THCA), and uterine carcinosarcoma (UCEC). Conversely, it was downregulated in kidney chromophobe (KICH) and kidney renal clear cell carcinoma (KIRC) compared to normal tissues (p<0.001). Further analysis confirmed that GRB7 expression in KICH and KIRC was significantly downregulated across various clinicopathological parameters including stage 3 and stage 4 compared to stage 1. It was also significantly downregulated in 61-80 years compared to 41-60 years patients, as confirmed by the immunohistochemistry of kidney tissues. Prognostic analysis revealed that high GRB7 expression was linked to a better prognosis in KIRC and a poorer prognosis in pancreatic adenocarcinoma (PAAD) patients. In KICH, GRB7 expression showed a significant positive correlation with immune infiltration of B cells, CD8+ T cells, and macrophages. In KIRC, GRB7 was positively correlated with immune infiltration of B cells and CD4+ cells. However, in PAAD it was negatively correlated with immune infiltration of macrophages. These findings were validated by gene expression profiling from the Gene Expression Omnibus (GEO) database, confirming a significant GRB7 downregulation in KICH and KIRC and an upregulation in PAAD compared to normal samples.
Conclusion: GRB7 shows potential as a biomarker in both diagnosing and predicting outcomes for various cancers. It may serve as a diagnostic marker for KICH, a prognostic marker for PAAD, and both a diagnostic and prognostic marker for KIRC, making GRB7 a target for future research and therapeutic approaches in oncology.
Article type: Research Article
Keywords: bioinformatics analysis, diagnostic biomarker, grb7, immune infiltration, kich, kirc, paad, pan-cancer, prognostic biomarker
License: Copyright © 2024, Attaelmanan et al. CC BY 4.0 This is an open access article distributed under the terms of the Creative Commons Attribution License CC-BY 4.0., which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Article links: DOI: 10.7759/cureus.74907 | PubMed: 39742197 | PMC: PMC11687406
Relevance: Relevant: mentioned in keywords or abstract
Full text: PDF (9.4 MB)
Introduction
Global cancer data show that cancer is a major global health burden. In 2022, there were over 10 million cancer-related fatalities and roughly 20 million new instances of the disease. According to expectations, there will be 35 million new instances of cancer by 2050, a 77% rise from 2022 [ref. 1].
With the rising incidence and mortality of cancer, understanding its onset, progression, and therapeutic strategies has become crucial. The development of new drugs targeting proteins within essential signaling pathways has accelerated due to advances in genomic research. These targeted inhibitors have demonstrated better efficacy and fewer side effects compared to traditional treatments. Among the promising targets are adaptor proteins, specifically those containing Src-homology 2 (SH2) domains [ref. 2]. Growth factor receptor-bound protein 7 (GRB7) is one such adaptor protein, part of a family that also includes GRB10 and GRB14. This family shares a conserved multidomain structure, consisting of a proline-rich region at the N-terminus, a pleckstrin homology domain, a Ras-associated domain, and a Src-homology 2 (SH2) domain [ref. 3]. Research has demonstrated that these GRB7 domains are an essential component of cellular signaling pathways and have a major impact on controlling cell growth, division, and proliferation [ref. 4]. GRB7 has been linked to growth factor receptors like EGFR and HER2, highlighting its involvement in a range of cellular processes [ref. 5]. It influences the proliferation, survival, and metastasis of cancer cells by interacting with the ErbB family of receptors [ref. 6]. GRB7 controls stress reactions and apoptotic pathways to promote cell survival. Because it knocks down the amounts of BCL-2 family proteins, it increases apoptosis [ref. 4,ref. 7]. GRB7 is associated with anchorage-independent growth and cell cycle progression in terms of proliferation, particularly in malignancies that overexpress ErbB23. High GRB7 expression increases the capacity for cancer cells to migrate, invade, and spread via interacting with different kinases and phospholipids [ref. 4]. GRB7 is the focus of much research in the domains of molecular biology and cancer biology due to its many roles. The GRB7 gene has been found to be dysregulated in various human cancers, indicating its potential as both a biomarker and a therapeutic target [ref. 7,ref. 8]. However, despite these findings, the functional role of GRB7 in cancer remains largely unexplored, especially in the context of its potential importance for molecular-targeted therapies in cancer treatment. Further analysis of the role of GRB7 in cancer could lead to the development of new therapeutic approaches. This work sought to explore the potential of GRB7 in cancer diagnosis and prognosis by performing a thorough pan-cancer investigation using several bioinformatics methods.
Materials and methods
Gene expression analysis
GRB7 expression across different cancers was evaluated using the Tumor Immune Estimation Resource (TIMER) 2.0 (http://timer.cistrome.org/) database. It is an analysis network for tumor immune cell infiltration and can also perform differential gene expression analysis between normal and tumor tissues [ref. 9]. The GRB7 gene’s differential expression across 32 types of tumors obtained from The Cancer Genome Atlas (TCGA) was analyzed using the TIMER database, and the Wilcoxon test was applied to determine the statistical significance of differential expression. Moreover, The Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/) database was used. The GEPIA database integrates data from TCGA and Genotype-Tissue Expression (GTEx). The analysis of the GRB7 gene’s differential expression across 33 types of tumors was performed using a cutoff of 0.05 for the p-value and 1.0 for the |log2(FC)| cutoff [ref. 10]. Further analysis was conducted using the UALCAN database (https://ualcan.path.uab.edu/) to assess the differential GRB7 gene expression in normal and cancer tissues and also to determine GRB7 expression at the protein level across various tumor types, leveraging data from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) datasets (http://ualcan.path.uab.edu/analysis-prot.html), and to explore the relationship between GRB7 expression and clinical pathological parameters of cancer [ref. 11]. Moreover, the Human Protein Atlas (HPA) (https://www.proteinatlas.org/) was employed to illustrate the histological differences in cancer tissues [ref. 12].
Prognostic value of GRB7 in cancers
The prognostic value of the GRB7 gene in various cancers was evaluated using the Kaplan-Meier plotter database (https://kmplot.com/analysis) [ref. 13] and the GEPIA2 database (http://gepia.cancer-pku.cn/) [ref. 10]. The survival analysis module in GEPIA assessed overall survival (OS) and disease-free survival (DFS) based on GRB7 expression in patients with different tumor types from the TCGA database [ref. 4]. Furthermore, the relationship between GRB7 expression and both overall survival and relapse-free survival (RFS) across 21 cancers was derived from the Kaplan-Meier plotter. Hazard ratios (HRs) with 95% confidence intervals (CIs) were computed, and log-rank p-values were used to assess statistical significance. A p-value <0.05 was considered statistically significant for prognostic outcomes. Furthermore, the overall survival was examined using the UALCAN database (http://ualcan.path.uab/index.hrml) [ref. 11].
Immune cell infiltration analysis
TIMER (https://cistrome.shinyapps.io/timer/) is an online website applied to characterize immune infiltration patterns across various cancer types [ref. 9]. This study utilized the TIMER tool to explore the correlation between GRB7 expression and immune cell infiltration. The analysis focused on estimating the abundance of six types of immune cells – B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells across various cancer types.
Genetic alteration analysis
The cBioPortal database (https://www.cbioportal.org/) was used to examine genetic alterations of the GRB7 gene across pan-cancer TCGA datasets [ref. 14]. By selecting the "TCGA Pan-Cancer Atlas Studies" and entering the "GRB7" gene through the query function, the search yielded a total of 10,967 samples from 32 studies.
Validation
Regarding data validation, publicly available datasets from the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/) were utilized [ref. 15]. The GEO2R tool (https://www.ncbi.nlm.nih.gov/geo/geo2r/) was applied to compare different sample groups for identifying differentially expressed genes. The significance of GRB7 in kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), and pancreatic adenocarcinoma (PAAD) was confirmed using several GEO datasets, including GSE11024, GSE53757, and GSE16515. Differential expression was visualized using volcano plots from the bioinformatics.com.cn platform [ref. 16]. Statistical significance was determined based on the criteria p-value <0.05, adjusted p-value <0.05, and |log2(FC)| >1 for detecting differentially expressed genes.
Results
GRB7 gene expression analysis
GRB7 gene expression was assessed in cancerous and normal tissues using TIMER 2.0. The analysis revealed that GRB7 expression was significantly elevated across multiple cancer types, including bladder urothelial carcinoma (BLCA), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), rectum adenocarcinoma (READ), stomach adenocarcinoma (STAD), thyroid carcinoma (THCA), uterine corpus endometrial carcinoma (UCEC) (p<0.001), breast invasive carcinoma (BRCA), cervical squamous cell carcinoma (CESC) (p<0.01), and kidney renal papillary cell carcinoma (KIRP) (p<0.05). In contrast, the expression of GRB7 was significantly reduced in glioblastoma multiforme (GBM) (p<0.01), KICH, and KIRC (p<0.001), as shown in Figure 1A. Additional analysis performed using the GEPIA database showed that GRB7 expression was significantly elevated in BLCA, CESC, CHOL, COAD, OV, READ, testicular germ cell tumors (TGCT), THCA, thymoma (THYM), UCES, and uterine carcinosarcoma (UCS) (Figure 1B). On the other hand, in KICH, KIRC, acute myeloid leukemia (LAML), skin cutaneous melanoma (SKCM), and sarcoma (SARC) the expression levels were lower compared to normal tissues, as shown in Figure 1C.

Due to the lack of normal controls for some tumors in the prior databases, a supplemental analysis was performed using the UALCAN database (Figures 2A–2Q). It indicated that GRB7 expression was significantly elevated in BRCA, CHOL, COAD, KIRP, LIHC, LUAD, and THCA (p<1×10-13), CESC, STAD (p<1×10-11), BLCA, ESCA, head and neck squamous cell carcinoma (HNSC), READ, UCEC (p<1×10-9), and LUSC (p<1×10-5). Conversely, GRB7 was significantly downregulated in KIRC and KICH (p<1×10-9). Whereas the expression of GRB7 in adrenocortical carcinoma (ACC), lymphoid neoplasm diffuse large B cell lymphoma (DLBC), mesothelioma (MESO), pancreatic adenocarcinoma (PAAD), and uveal melanoma (UVM) was found to be insignificant in all three databases.

Moreover, the proteomic expression profiles of the GRB7 gene across various tumors were assessed using the UALCAN database. It was observed that GRB7 proteomic expression showed a significant upregulation in COAD, PAAD, lung cancer, and UCEC samples, while a decrease was noted in KIRC (p<0.001) (Figure 3). By cross-referencing the three databases (GEPIA, TIMER, and UALCAN), GRB7 expression was found to be significantly higher in BLCA, CESC, CHOL, COAD, READ, THCA, and UCEC. Conversely, it was significantly downregulated in KICH and KIRC.

Correlation between clinicopathological parameters and GRB7 expression in KICH and KIRC
Considering the results obtained from the above analysis, KICH and KIRC were selected for further analysis using the UALCAN database. The correlation between clinicopathological parameters (stage, race, gender, and age) and GRB7 expression in KICH was investigated. The gene expression was significantly downregulated in stage 1 (p=1.66×10-8), stage 2 (p=4.52×10-8), stage 3 (p=6.81×10-9), and stage 4 (p=2.55×10-4) (Figure 4A). Furthermore, the gene expression profile showed a clear significance based on race (Figure 4B), with the gene GRB7 downregulated in Caucasians (p=1.50×10-8), and African Americans (p=9.92×10-3). Additionally, the gene expression profile showed high significance in the patient’s gender (Figure 4C), with the gene GRB7 expressed in males and females, (p=1.38×10-8, p=2.77×10-8, respectively). Moreover, GRB7 was significantly downregulated in age groups 21-40 years (p=3.79×10-8), 41-60 years (p=2.12×10-8), and 61-80 years (p=8.12×10-9), as shown in Figure 4D. A detailed analysis of KIRC using the UALCAN database revealed significant alterations in gene expression across various clinicopathological parameters. Gene expression was consistently downregulated in all stages of KIRC compared to normal tissue, with highly significant p-values observed in stage 1 (p=1.62×10-12), stage 2 (p=1.52×10-11), stage 3 (p<1×10-12), stage 4 (p=1.62×10-12). Additionally, notable differences were detected between specific stages, such as stage 1 versus stage 3 (p=2.23×10-4) and stage 1 versus stage 4 (p=9.22×10-4) (Figure 4E). When analyzing different racial groups, significant downregulation in gene expression was found in both Caucasians (p<1×10-12) and African Americans (p=1.19×10-10), while no significant difference was observed in Asians (p=7.04×10-1). Furthermore, a significant difference was noted between Caucasians and African Americans (p=6.43×10-5) (Figure 4F). Gender-specific findings revealed significant downregulation in both males (p=1.62×10-12) and females (p<1×10-16), though there was no substantial difference between genders (Figure 4G). The age-based analysis demonstrated that gene expression was significantly lower across all age groups compared to normal, particularly in those aged 21-40 years (p=2.31×10-7), 41-60 years (p=1.11×10-16), 61-80 years (p<1×10-12), and 81-100 years (p=1.11×10-16). Furthermore, a significant difference was noted between the 41-60 years and 61-80 years age groups (p=1.84×10-2) (Figure 4H). In terms of cancer grade, gene expression was significantly downregulated across all grades compared to normal, with p-values as follows: grade 1 (p=1.20×10-5), grade 2 (p=2.55×10-16), grade 3 (p=1.62×10-12), and grade 4 (p<1×10-12). Additionally, significant differences were observed between specific grades – grade 1 versus grade 3 (p=1.39×10-2), grade 1 versus grade 4 (p=5.12×10-3), grade 2 versus grade 3 (p=2.48×10-2), grade 2 versus grade 4 (p=1.85×10-8), and grade 3 versus grade 4 (p=9.69×10-5) (Figure 4I). Moreover, the histology of kidney cancer tissue has been explored using the Human Protein Atlas showing the tumor elements of KICH and KRIC (Figures 5A–5J). In addition, cancer tissues were also compared between 41-60 years and 61-80 years age groups of patients with kidney cancer (Figures 5E–5J).


Prognostic value of GRB7 in cancers
Patients with cancer were categorized into groups with high and low GRB7 expression levels, and their overall outcomes were analyzed to evaluate the prognostic value of GRB7 (Figures 6A–6J). The GEPIA analysis indicated that higher GRB7 expression was linked with poorer overall survival (OS) in OV and PAAD (Figures 6C, 6D). However, increased GRB7 expression was correlated with better OS in patients with ACC and KIRC (p<0.05) (Figures 6A, 6B). Additionally, elevated GRB7 expression was related to poorer Disease-free survival (DFS) in patients with CHOL, COAD, LGG, and PRAD (Figures 6F, 6G, 6I, 6J). Meanwhile, higher GRB7 expression was associated with improved DFS in patients with ACC and KIRC (p<0.05) (Figure 6E).

The link between GRB7 expression and OS, or relapse-free survival (RFS), was explored using the Kaplan-Meier plotter database (Figures 7A–7X). The analysis showed that increased GRB7 expression was correlated with significantly poorer OS and RFS in OV (OS: HR=1.61, p=0.00037; RFS: HR=1.54, p=0.019), PAAD (OS: HR=2.02, p=0.00068; RFS: HR=3.98, p=0.016), TGCT (OS: HR=582885929.78, p=0.048; RFS: HR=2.31, p=0.026), and UCEC (OS: HR=2.53, p=5.5×10-6; RFS: HR=2.24, p=0.0019). In addition, patients with higher GRB7 expression had significantly poorer OS in LUSC (HR=0.73, p=0.044) and THYM (HR=4.76, p=0.024). However, significantly better OS and RFS were noted in HNSC (OS: HR=0.7, p=0.023; RFS: HR=0.3, p=0.017), and SARC (OS: HR=0.62, p=0.022; RFS: HR=0.55, p=0.014) for individuals with higher GRB7 expression. Furthermore, a significantly better OS was observed in BLCA (HR=0.73, p=0.044), ESCA (HR=0.43, p=0.014), KIRC (HR=0.64, p=0.00046), STAD (HR=0.65, p=0.012) and THCA (HR=0.27, p=0.0055). Whereas in CESC (OS: HR=0.58, p=0.021; RFS: HR=2.88, p=0.0063) and READ (OS: HR=0.41, p=0.032; RFS: HR=9.53, p=0.012), patients with higher GRB7 expression had better OS and poorer RFS.

The relationship between GRB7 expression and overall survival (OS) was further examined using the UALCAN database. This analysis confirmed that higher GRB7 expression was linked with significantly better OS in HNSC (p=0.029), KIRC (p=0.00048), and UVM (p=0.02). The worst OS was noted in LGG (p=0.0024), PAAD (p=0.03), SKCM (p=0.016), THYM (p=0.043), and UCEC (p<0.0001) (Figures 8A–8I). Combining findings from GEPIA, TIMER, and UALCAN, high GRB7 expression was related to poor OS in PAAD, while it was linked to better OS in KIRC. This highlights the potential of GRB7 as a prognostic biomarker in these cancer types.

Correlation analysis of GRB7 expression and immune infiltration of immune cells
The prognosis of tumors is influenced by the interaction of immune cells. In this study, the relationship between GRB7 expression in KICH, KIRC, and PAAD and the infiltration of immune cells, such as CD8+ T cells, CD4+ T cells, B cells, neutrophils, macrophages, and dendritic cells was investigated. In KICH, GRB7 expression was found to have a significantly positively correlated with B cell infiltration (partial. Cor=0.336, p=6.14×10-3), CD8+ cells (partial. Cor=0.343, p=5.21×10-3), and macrophage (partial. Cor=0.311, p=1.71×10-3), as shown in Figure 9A. Similarly, in KIRC (Figure 9B), GRB7 expression has significantly positively correlated with B cells (partial. Cor=0.134, p=4.11×10-3) and CD4+ cells (partial. Cor=0.182, p=8.99×10-5). While in PAAD a significant negative correlation between GRB7 expression and macrophage infiltration was observed (partial. Cor= -0.224, p=3.25×10-3), as shown in Figure 9C.

Genetic alteration analysis
It is well-established that tumorigenesis is strongly linked with genetic alterations [ref. 7]. To examine the genetic alteration of GRB7 in human cancer, a comparative analysis was performed using the cBioPortal database. Among all 10,967 samples queried for 10,953 patients in 32 studies, GRB7 was found to be altered in 463 samples, accounting for approximately 4%. Amplification was identified as the most common type of genetic alteration in GRB7. For ESCA, STAD, BRCA, UCS, and UCEC, the highest frequencies of genetic alterations were observed, with rates of 15.93%, 15%, 11.72%, 8.77%, and 8.7%, respectively. Notably, mutations were more prevalent in SKCM, accounting for 3.83% (17/444 cases), and in UCEC, constituting 3.4% (18/529 cases), as shown in Figure 10A. In 10,967 samples, 113 mutation sites were identified between amino acids 0 and 532, including 88 missense mutations, seven truncating mutations, five splicing mutations, and 13 fusions, as shown in Figure 10C. The X424_splice/A424T mutation was found to be the most frequent. This mutation was situated between the BPS and the SH2 domain (Figure 10B).

Validation
To validate our findings on GRB7 expression, gene expression profiling was performed using data from the Gene Expression Omnibus (GEO) database. A volcano plot was created with a bioinformatics tool (http://www.bioinformatics.com.cn) to visually represent the upregulation or downregulation of GRB7 among the differentially expressed genes (DEGs). In KICH, the study identified differentially expressed genes using the GEO2R tool, which compares sample groups to pinpoint significant differences in gene expression. The GSE11024 dataset was used, consisting of 12 normal samples in the control group and six tumor samples in the tumor group. GRB7 was significantly downregulated (adjusted p-value = 3.25×10-4; |log2(FC)| = -0.8512116). The analysis revealed 10,223 DEGs, of which 6,553 were upregulated and 3,670 were downregulated (Figure 11A). For KIRC, the GSE53757 dataset was used, the control group contained 72 normal samples, and the tumor group contained 72 samples. GRB7 was notably downregulated (adjusted p-value = 2.29×10-17; |log2(FC)| = -0.8176622). The analysis showed 28,064 DEGs, of which 12,974 were upregulated and 15,090 were downregulated (Figure 11B). In PAAD, the GSE16515 dataset was used, it contains 36 tumor samples and 16 normal samples. GRB7 was notably upregulated (adjusted p-value = 0.630; |log2(FC)| = 2.43×10-2). The analysis showed 20,809 DEGs, of which 7,934 were upregulated and 12,875 were downregulated (Figure 11C).

Discussion
The identification of unique markers related to abnormal gene expressions is a critical approach for tumor diagnosis and treatment. GRB7 is located within the 17q12-21 amplicon in proximity to the HER2 gene [ref. 5], it plays a significant role in cellular signaling by interacting with various receptor tyrosine kinases and their downstream partners [ref. 4]. Consequently, GRB7 functions as a pivotal hub, linking diverse potential drivers of oncogenesis to downstream signaling pathways, thus presenting an appealing target for therapeutic intervention. However, a comprehensive bioinformatics analysis of GRB7’s role across pan-cancer tissues using multiple databases has not yet been conducted.
This study examined the expression level of GRB7 in different malignancies using multiple databases, examining its association with gene expression, prognosis, immune infiltration, and gene alterations. The expression of the GRB7 gene was investigated simultaneously in three different databases. Statistically significant expression relative to normal tissue was observed in several cancers, with seven cancers (BLCA, CESC, CHOL, COAD, READ, THCA, and UCEC) showing significant upregulation across all databases. Conversely, GBR7 demonstrated a significant downregulation in expression levels in KICH and KIRC. These findings align with earlier research that identified GRB7 as being upregulated in various tumor types, such as cervical cancer [ref. 6], gastric cancer, bladder cancer, and hepatocellular carcinoma [ref. 17–ref. 19]. Further, the analysis in the UALCAN database confirmed significant downregulation of GRB7 expression in KICH across various clinicopathological parameters, including stage, race, and age, when compared to normal tissue. However, in KIRC, GRB7 expression showed significant variation, being lower in advanced stages compared to earlier stages, and also showed lower expression in Caucasians compared to African Americans when analyzing different racial groups. Additionally, older age groups show lower expression than younger age groups. A similar pattern was observed across cancer grades, as GRB7 expression decreased progressively with advancing grades. This difference in histology has been illustrated between kidney cancer tissues of patients aged 41-60 years and 61-80 years using antibodies CAB005226, CAB005226, and HPA057084. The crucial role of GRB7 in cancer prognosis was illustrated in many studies [ref. 17,ref. 18,ref. 20,ref. 21]. Although HER2 and GRB7 co-amplification was previously reported [ref. 22], another study done by Ramsey et al. suggested that GRB7 plays a unique role in oncogenic pathways independent of HER2, making it a potential standalone biomarker [ref. 23]. In this study, the significant role of GRB7 expression in prognosis was examined using multiple databases (Kaplan-Meier plotter, GEPIA, and UALCAN). Our results revealed that high GRB7 levels were connected to poorer OS in PAAD. This finding aligns with previous studies showing that increased expression of GRB7 protein is associated with regional lymph node metastasis in pancreatic cancer [ref. 24]. However, in KIRC, higher GRB7 expression is correlated with improved OS.
The tumor microenvironment (TME) plays a vital role in both tumor progression and the effectiveness of immunotherapy. It was found that the knockout of GRB7 was associated with increased T cell-mediated cytotoxicity, suggesting the possibility of GRB7 being a potential target for immunotherapy [ref. 25]. In this study, the relationship between GRB7 expression and immune cell infiltration in KICH, KIRC, and PAAD was examined using the TIMER database. For instance, in KICH, a significant positive correlation was observed between GRB7 expression and the infiltration of B cells, CD8+ T cells, and macrophages. While in KIRC, GRB7 expression was also positively correlated with the infiltration of B cells and CD4+ T cells. However, in PAAD, GRB7 expression showed a negative correlation with macrophage infiltration. These results indicate that GRB7 expression is linked to distinct patterns of immune cell infiltration within the TME, which may affect tumor progression and prognosis. Specifically, higher levels of GRB7 may promote the infiltration of B cells, CD8+ T cells, and macrophages while reducing macrophage presence in some cancers, suggesting that GRB7 plays a complex role in shaping the immune landscape across different cancers.
Significant findings regarding GRB7 were identified by analyzing genetic alteration data across various cancers using the cBioPortal platform. GRB7 exhibited genetic alterations in multiple cancer types, with amplification being the predominant alteration detected in ESCA, while mutations were notably prevalent in SKCM and UCEC. However, no significant impact of GRB7 mutations was observed in KICH, KIRC, and PAAD. Additionally, GRB7 gene expression in KICH, KIRC, and PAAD was validated using data from the GEO database and various online resources. The results showed that GRB7 expression was downregulated in KICH and KIRC, while it was upregulated in PAAD.
The primary limitation of the study is its dependence on computational datasets, making it essential to perform wet lab experiments to confirm these computationally derived results. Further experimental validation is needed to determine whether GRB7 is an independent adverse prognostic factor in PAAD and KIRC. In addition to exploring GRB7’s role in modulating immune cell infiltration across the TME in KICH, KIRC, and PAAD to provide critical insights into its role in tumor progression and its potential as a target for immunotherapy-based strategies. These experimental studies will provide valuable insights into the potential of GRB7 as a therapeutic target.
Conclusions
This study highlights that GRB7 was downregulated in KICH, suggesting its potential as a diagnostic biomarker for this cancer. Additionally, it could serve as a prognostic biomarker in PAAD and act as both a diagnostic and prognostic biomarker in KIRC. Moreover, GRB7 gene expression was observed to be associated with immune infiltration, highlighting its potential effect on the tumor microenvironment. The above results pinpoint the importance of GRB7 as a target for future research and therapeutic strategies in oncology.
References
- H Sung, J Ferlay, RL Siegel, M Laversanne, I Soerjomataram, A Jemal, F Bray. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin, 2021. [PubMed]
- SC Pero, L Oligino, RJ Daly. Identification of novel non-phosphorylated ligands, which bind selectively to the SH2 domain of Grb7. J Biol Chem, 2002. [PubMed]
- LJ Holt, K Siddle. Grb10 and Grb14: enigmatic regulators of insulin action – and more?. Biochem J, 2005. [PubMed]
- DC Han, TL Shen, JL Guan. The Grb7 family proteins: structure, interactions with other signaling molecules and potential cellular functions. Oncogene, 2001. [PubMed]
- WW Bivin, O Yergiyev, ML Bunker, JF Silverman, U Krishnamurti. GRB7 expression and correlation with HER2 amplification in invasive breast carcinoma. Appl Immunohistochem Mol Morpho, 2017
- PY Chu, YL Tai, TL Shen. Grb7, a critical mediator of EGFR/ErbB signaling, in cancer development and as a potential therapeutic target. Cells, 2019
- HB Zhao, XF Zhang, XL Jia, HB Wang. Grb7 is over-expressed in cervical cancer and facilitate invasion and inhibit apoptosis in cervical cancer cells. Pathol Res Pract, 2017. [PubMed]
- Y Nadler, AM González, RL Camp, DL Rimm, HM Kluger, Y Kluger. Growth factor receptor-bound protein-7 (Grb7) as a prognostic marker and therapeutic target in breast cancer. Ann Oncol, 2010. [PubMed]
- T Li, J Fan, B Wang. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res, 2017
- Z Tang, B Kang, C Li, T Chen, Z Zhang. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res, 2019
- DS Chandrashekar, B Bashel, SA Balasubramanya, CJ Creighton, I Ponce-Rodriguez, BV Chakravarthi, S Varambally. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia, 2017. [PubMed]
- M Uhlén, L Fagerberg, BM Hallström. Proteomics. Tissue-based map of the human proteome. Science, 2015
- A Lánczky, B Győrffy. Web-based survival analysis tool tailored for medical research (KMplot): development and implementation. J Med Internet Res, 2021
- E Cerami, J Gao, U Dogrusoz. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov, 2012. [PubMed]
- T Barrett, SE Wilhite, P Ledoux. NCBI GEO: archive for functional genomics data sets – update. Nucleic Acids Res, 2013
- D Tang, M Chen, X Huang. SRplot: a free online platform for data visualization and graphing. PLoS One, 2023
- G Zhu, H Cai, Q Xiao, S Zeng, X Jiang, D Liu. GRB7 plays a promoting role in the progression of gastric cancer. BMC Cancer, 2023
- Y Zheng, Y Pei, L Yang. Upregulated GRB7 promotes proliferation and tumorigenesis of bladder cancer via phospho-AKT pathway. Int J Biol Sci, 2020. [PubMed]
- S Itoh, A Taketomi, S Tanaka. Role of growth factor receptor bound protein 7 in hepatocellular carcinoma. Mol Cancer Res, 2007. [PubMed]
- O Giricz, V Calvo, SC Pero, DN Krag, JA Sparano, PA Kenny. GRB7 is required for triple-negative breast cancer cell invasion and survival. Breast Cancer Res Treat, 2012. [PubMed]
- JR Gotovac, DS Liu, MJ Yates. GRB7 is an oncogenic driver and potential therapeutic target in oesophageal adenocarcinoma. J Pathol, 2020. [PubMed]
- PJ Lamy, F Fina, C Bascoul-Mollevi, AC Laberenne, PM Martin, L Ouafik, W Jacot. Quantification and clinical relevance of gene amplification at chromosome 17q12-q21 in human epidermal growth factor receptor 2-amplified breast cancers. Breast Cancer Res, 2011
- B Ramsey, T Bai, A Hanlon Newell. GRB7 protein over-expression and clinical outcome in breast cancer. Breast Cancer Res Treat, 2011. [PubMed]
- S Tanaka, SC Pero, K Taguchi, M Shimada, M Mori, DN Krag, S Arii. Specific peptide ligand for Grb7 signal transduction protein and pancreatic cancer metastasis. J Natl Cancer Inst, 2006. [PubMed]
- L Wen, W Hu, S Hou. GRB7 plays a vital role in promoting the progression and mediating immune evasion of ovarian cancer. Pharmaceuticals (Basel), 2024
