- Research
- Open access
- Published:
CYFIP2 serves as a prognostic biomarker and correlates with tumor immune microenvironment in human cancers
European Journal of Medical Research volume 28, Article number: 364 (2023)
Abstract
Background
The mechanisms whereby CYFIP2 acts in tumor development and drives immune infiltration have been poorly explored. Thus, this study aimed to identifying the role of CYFIP2 in tumors and immune response.
Methods
In this study, we first explored expression patterns, diagnostic role and prognostic value of CYFIP2 in cancers, particularly in lung adenocarcinoma (LUAD). Then, we performed functional enrichment, genetic alterations, DNA methylation analysis, and immune cell infiltration analysis of CYFIP2 to uncover its potential mechanisms involved in immune microenvironment.
Results
We found that CYFIP2 significantly differentially expressed in different tumors including LUAD compared with normal tissues. Furthermore, CYFIP2 was found to be significantly correlated with clinical parameters in LUAD. According to the diagnostic and survival analysis, CYFIP2 may be employed as a potential diagnostic and prognostic biomarker. Moreover, genetic alterations revealed that mutation of CYFIP2 was the main types of alterations in different cancers. DNA methylation analysis indicated that CYFIP2 mRNA expression correlated with hypomethylation. Afterwards, functional enrichment analysis uncovered that CYFIP2 was involved in tumor-associated and immune-related pathways. Immune infiltration analysis indicated that CYFIP2 was significantly correlated with immune cells infiltration. In particular, CYFIP2 was strongly linked with immune microenvironment scores. Additionally, CYFIP2 exhibited a significant relationship with immune regulators and immune-related genes including chemokines, chemokines receptors, and MHC genes.
Conclusion
Our results suggested that CYFIP2 may serve as a prognostic cancer biomarker for determining prognosis and might be a promising therapeutic strategy for tumor immunotherapy.
Introduction
The occurrence and development of malignancies may be affected by many factors, such as genetics, living environment and habits. At present, the efficacy of tumor treatment is still far from satisfactory, probably because the occurrence of tumor is linked to genetic variants [1]. In-depth knowledge of the full range of functions of a gene is meaningful to help us better understand the mechanism of tumorigenesis and promote personalized treatment [2]. With the rapid development of basic research and clinical medicine, people’s understanding of the etiology and pathophysiological processes of diseases has been deepened, and more and more biomarkers involved in tumorigenesis, development and influencing prognosis have been discovered one after another, and play an important role in early diagnosis of tumors, evaluation of therapeutic efficacy and prediction of prognosis [3]. Sensitive and specific biomarkers are especially necessary in clinical precision medicine [4]. Biomarkers can also be used as potential targets for drug design. The discovery of powerful biomarkers and new therapeutic targets is critical to the development of the next generation of precision medicine. During past years, the improvement of high-throughput sequencing based cancer atlas initiatives and omics technologies has brought cancer research into a new era [5]. Computational or bioinformatics approaches have made it more efficient to investigate the ability and function of target molecules as well as their interactions, thereby contributing to a comprehensive understanding of cancer initiation and progression [6].
Lung adenocarcinoma (LUAD) is now the most common histological subtype of primary lung cancer accounting for greater than 40% of cases [7]. Despite achievements in the development of new approaches in the treatment of LUAD, unfortunately, LUAD is still one of the most aggressive and rapidly fatal tumor types with overall survival less than 5 years [8]. Thus, it is important to find stable and reliable biomarkers for LUAD to identify patients with poor prognosis and provide more aggressive treatment [9].
Study identified Cytoplasmic FMR Interacting Protein 2 (CYFIP2) as a p53-inducible gene which was a direct p53 target gene that may be part of a redundant network of genes responsible for p53-dependent apoptosis [10]. During the past decades, accumulating evidence has demonstrated the role of CYFIP2 in different diseases. For obesity and related metabolic disorders, data emphasized the potential of CYFIP2 as a pharmacotherapeutic target for treating obesity and other metabolic disorder [11]. Moreover, CYFIP protein family, including CYFIP1 and CYFIP2, is involved in neural development and maturation as well as in different neural disorders, such as intellectual disabilities, autistic spectrum disorders, and Alzheimer’s disease [12]. Recent studies have shown that CYFIP2 is also closely related to tumor initiation and progression. Fox example, NUAK family kinase 2 (NUAK2) may regulate CYFIP2 expression to promote cervical cancer cell proliferation, migration, invasion and epithelial-to-mesenchymal transition [13]. Circulating CYFIP2 was demonstrated to be significantly upregulated in gastric cancer tissues and cell lines and high expression of circulating CYFIP2 was associated with metastasis and poor prognosis of gastric cancer patients [14]. Studies also indicated that CYFIP2 may take part in the vascular endothelial growth factor receptor signaling pathway and the purine ribonucleoside diphosphate metabolic process [15]. However, the comprehensive function and prognostic landscape of CYFIP2 involved in human cancers are still not fully understood. Thus, understanding its mechanisms may be helpful to enhance diagnosis and treatment and may also promote its potential clinical practice.
Our study, for the first time, focused on the comprehensive landscape of CYFIP2 function based on a multi-omics analysis. We first revealed the expression patterns including its expression in normal tissues, cancer tissues and intracellular localization. Then, we evaluated the clinical utility of CYFIP2 by studying its associations with different clinicopathological characteristics and potential role in the diagnosis and prognosis of different cancers. Moreover, we investigated the potential molecular mechanisms of CYFIP2 in the pathogenesis or clinical prognosis of different cancers, especially in LUAD, by performing co-expression gene analysis, protein–protein iteration network analysis, functional enrichment analysis, genetic alterations and DNA methylation analysis. In addition, we uncovered the landscape of CYFIP2 expression in tumor immunity and microenvironment through the correlation analysis of CYFIP2 expression with immune cell infiltration, immune regulators, and immune-related genes (chemokines, chemokines receptors, and MHC genes). The schematic pipeline for this study is plotted at Fig. 1.
Materials and methods
Expression analysis of CYFIP2
The expression of CYFIP2 in normal tissues was analyzed based on GTEx database [16]. To assess the differences of CYFIP2 in protein expression level, the IHC image was downloaded and analyzed from the HPA database [17]. The intracellular localization of CYFIP2 was also identified by the HPA database. The expression differences of CYFIP2 between cancer and normal tissues in different cancer types were accomplished by the TIMER database [18]. We input CYFIP2 in the “Gene_DE”module of TIMER2 web and observed the expression difference of CYFIP2 between cancer and normal tissues for the different tumors of the TCGA project [19]. We further compared the expression difference of CYFIP2 in LUAD with normal tissues and its paired adjacent tissues.
Correlations between CYFIP2 and clinical features
Correlations between CYFIP2 and different clinical features such as stage, age, and gender in TCGA-LUAD cohort were analyzed using the Kruskal–Wallis test and visualized by violin plots. Moreover, correlations between CYFIP2 and different clinicopathological characteristics in LUAD were analyzed using the Chi-squared test or Fisher’s exact test and presented in the baseline datasheet.
Diagnostic analysis and survival analysis of CYFIP2
The diagnostic significance of CYFIP2 in LUAD was estimated using RNA sequencing data from TCGA and illustrated by receiver operating characteristic (ROC) curves using the “pROC” R package. The correlation between CYFIP2 expression and survival in pan-cancer was analyzed by Kaplan–Meier Plotter [20]. The relationship of CYFIP2 expression with overall survival (OS) was analyzed in LUAD, esophageal squamous cell carcinoma (ESCC), bladder carcinoma (BC), Kidney renal clear cell carcinoma (KIRC), Kidney renal papillary cell carcinoma (KIPP), pancreatic ductal adenocarcinoma (PDAC), lung squamous cell carcinoma (LUSC), head and neck squamous cell carcinoma (HNSC), sarcoma, thymoma. Hazard ratios (HRs) with 95% confidence intervals (CIs) and log-rank P-values were calculated. Moreover, multivariable risk analysis with CYFIP2 expression and other variables was performed to explore independent factors.
Genetic alterations and DNA methylation analysis of CYFIP2
The cBioPortal for Cancer Genomic provides a web resource for exploring, visualizing, and analyzing multidimensional cancer genomics data [21]. Alteration of CYFIP2 status in cancer patients and mutation site were acquired from this database. The UALCAN portal, an easy to use, interactive web-portal to perform to in-depth analyses of TCGA gene expression data, allows us to analyze the promoter methylation level of CYFIP2 in normal and primary tumor tissues [22].
Co-expression genes analysis
The LinkedOmics database was applied to screen the co-expression genes of CYFIP2 based on the TCGA-LUAD cohort [23]. The top 50 genes significantly related to CYFIP2 were identified and presented in a heat map and volcano plot. Then, the survival maps evaluating the relationships between top 50 genes most positively and negatively associated with CYFIP2 and overall survival (OS)/disease-free survival (DFS) in LUAD were generated.
PPI analysis and functional enrichment analysis of CYFIP2
The STRING database aims to integrate all known and predicted associations between proteins, including both physical interactions as well as functional associations [24]. In this study, we used STRING to search interacting genes and construct PPI networks. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of interacting genes of CYFIP2 were performed by the “ClusterProfiler” package and visualized by the “ggplot2” package [25, 26]. Moreover, significant network module of CYFIP2 PPI network was identified. Survival maps of the relationship between module genes expression and OS/DFS in LUAD were established. In addition, the correlation analyses between the expression of CYFIP2 and these genes in LUAD were performed by Spearman’s correlation analysis and visualized by scatter plot.
Correlation analysis of CYFIP2 expression and immune cell infiltration
TISIDB is an online web integrated multiple types of data resources in tumor immunology [27]. In this study, we performed TISIDB analysis to determine the expression of CYFIP2 and tumor-infiltrating lymphocytes (TILs) across human cancers. Based on the gene expression profile, the relative abundance of TILs was inferred using gene set variation analysis. The correlations between CYFIP2 and TILs were measured by Spearman’s test. The relationships between CYFIP2 expression and immune infiltration were also determined using the TIMER tool.
Correlation analysis of CYFIP2 and tumor microenvironment
The tumor microenvironment (TME) is composed of a variety of cells including stromal cells and immune cells, which is significant to efficiently identify appropriate patients for immunotherapies [28]. Correlation analysis of CYFIP2 expression and immune microenvironment scores (stromal score, immune score, and estimate score) were explored with R package “estimate” and presented by scatter plot.
Correlation analysis of CYFIP2 and immune related genes
The immune regulators and related genes including chemokines, chemokine receptors and MHC genes may play an important role in the tumor initiation and immunity. The correlation analyses of immune regulators and related genes with CYFIP2 were assessed by Spearman’s correlation analysis and visualized by heatmap and radar map. P-value < 0.05 was regarded as the significant threshold.
Results
Data information
The pan-cancer analysis was carried out in 32 different cancer types. The sample number for each cancer type is listed in Table 1. The TCGA-LUAD cohort was divided into two groups including low- and high- CYFIP2 expression level based on the median expression value. The detailed clinicopathological characteristics are presented at Table 2.
The expression patterns of CYFIP2
To understand the expression and functions of CYFIP2 in different cancers, we should first know its expression in normal tissues. Therefore, the expression level of CYFIP2 in 29 types of normal tissues was first analyzed on the basis of the GTEx database. The results indicated that CYFIP2 expression varied across different types of normal tissues with an expression above medium level of expression in lung tissues (Fig. 2A).
Then, we evaluated the expression of CYFIP2 in tumor tissues using the HPA database. The representative IHC images of CYFIP2 protein expression of LUAD tissues with low-, medium- and strong staining are presented in Fig. 2B–D.
The intracellular localization of CYFIP2 was evaluated in HEK293, HeLa and U2OS cells by immunofluorescent staining of microtubules, endoplasmic reticulum (ER), and nucleus according to the HPA database. Displayed in Fig. 2E, CYFIP2 was found to be localized on the endoplasmic reticulum, plasma membrane and cytosol.
CYFIP2 is linked with clinical features in different cancer types
We then applied the TIMER2 approach to analyze human CYFIP2 expression levels in different cancer types based on TCGA data. As shown in Fig. 3A, the expression levels of CYFIP2 in the tumor tissues of BLCA (P < 0.001), ESCA (P < 0.001), KICH (P < 0.001), KIRC (P < 0.001), KIRP (P < 0.001), LUAD (P < 0.01), LUSC (P < 0.001), and UCEC (P < 0.05) were significantly lower than the corresponding control tissues. In contrast, the expression levels of CYFIP2 in the tumor tissues of BRCA (P < 0.001), CHOL (P < 0.001), PRAD (P < 0.001) and THCA (P < 0.001) were significantly higher than the corresponding control tissues.
We further explored the clinical significance of CYFIP2 expression in LUAD. The results showed that CYFIP2 was markedly downregulated in LUAD tumor tissue compared with normal tissue (Fig. 3B) and paired adjacent normal tissues (Fig. 3C). The diagnostic potential of CYFIP2 in LUAD was then evaluated. ROC analysis indicated that the AUC of CYFIP2 expression in diagnosing LUAD was 0.656 (95% CI 0.606–0.707) based on TCGA data (Fig. 3D). Interestingly, when combining TCGA and GTEx data (Fig. 3E), the diagnostic power of CYFIP2 in LUAD became more significant with an AUC of 0.766 (95% CI 0.735–0.797). Furthermore, we found that CYFIP2 expression was correlated with some clinicopathological characteristics such as T stages, age, gender, and smoking status while no significant association was detected with pathologic stages (Fig. 3F–J).
CYFIP2 is associated with survival outcomes in different cancer types
We used Kaplan–Meier Plotter, which is mainly based on Affymetrix microarray information from TCGA, to assess CYFIP2-related OS in different cancers (Fig. 4A–H). High expression of CYFIP2 was identified as a favorable prognostic factor in LUAD (HR = 0.61, 95% CI 0.45–0.84, log-rank P < 0.01), thymoma (HR = 0.06, 95% CI 0.01–0.47, log-rank P < 0.01), sarcoma (HR = 0.54, 95% CI 0.36–0.82, log-rank P < 0.01), PDAC (HR = 0.43, 95% CI 0.28–0.67, log-rank P < 0.01), clear cell RCC (HR = 0.34, 95% CI 0.25–0.46, log-rank P < 0.01), BC (HR = 0.73, 95% CI 0.54–0.98, log-rank P = 0.035). For UCEC and ESCC, CYFIP2 was found to have an unfavorable effect on OS (HR = 2.74, 95% CI 1.76–4.26, log-rank P < 0.01; HR = 3.21, 95% CI 1.17–8.79, log-rank P = 0.018).
We further explored the potential relationship between CYFIP2 expression and several clinical features of patients with LUAD (Fig. 4I). The results of multivariate regression analyses showed that high expression of CYFIP2 played a benefited role in patients with LUAD (HR = 0.668, 95% CI 0.449–0.994, P = 0.047). Some characteristics also played a detrimental role in patients with LUAD, such as T3&T4 stage (HR = 1.883, 95% CI 1.124–3.153, P = 0.016), lymph node-positive disease (HR = 1.899, 95% CI 1.294–2.787, P = 0.001), progressive disease (PD) & stable disease (SD) (HR = 2.760, 95% CI 1.841–4.138, P < 0.001).
CYFIP2 is mutated with missense in different cancer types
We performed comparative analysis of CYFIP2 using cBioPortal database to figure out genomic mutation of CYFIP2 in different cancers (Fig. 5A). The genetic alteration profiling of CYFIP2 showed that its mutation is one of the most important single factors for alteration in UCEC, SKCM, STAD, UCS, LUAD, ACC, COAD, BLCA, and LUSC. In addition, CYFIP2 amplification frequencies are the highest in KIRC, CHOL, PAAD, and PCPG. The types, sites and case numbers of the CYFIP2 genetic alteration are further presented in Fig. 5B. Missense mutation of CYFIP2 is the main type of genetic alteration.
CYFIP2 mRNA expression correlates with hypomethylation
DNA methylation has been recognized as one of the most important ways in the regulation of gene transcription, and thus we supposed that dysregulation of CYFIP2 expression may be regulated in this way. Accordingly, as shown in Fig. 5C, the results revealed that gene methylation negatively correlated with CYFIP2 gene expression (P < 0.05) in TCGA-LUAD cohort. In addition, the promoter methylation of CYFIP2 in tumor tissues of TCGA-LUAD was significantly lower than that of normal tissues adjacent to the cancer regardless of different clinical stages, N stages, gender, age and smoking status (Fig. 5D–H).
CYFIP2 has a synergistic effect with its co-expressed genes
For gaining the in-depth knowledge of CYFIP2 biological function in LUAD, we used the LinkedOmics database to study the correlated significant genes based on TCGA-LUAD cohort. The whole genes significantly related CYFIP2 are identified and presented in the volcano plot (Fig. 6A). As shown in the heat map (Fig. 6B, C), the top 50 marked genes were positively and negatively associated with CYFIP2. We further explored the prognostic role of the marked genes positively and negatively related to CYFIP2. As plotted by the survival maps (Fig. 6D–G), most of the top 50 positively correlated genes had favorable protective HRs for both OS and DFS in LUAD, while nearly all the negatively correlated genes had an unfavorable protective HR for both OS and DFS, suggesting these co-expressed genes highly owned the probability of becoming potential biomarkers in LUAD. These results indicated that CYFIP2 may have a synergistic effect with its co-expressed genes in the prognosis of LUAD.
CYFIP2 is implicated in the regulation of signaling pathways
CYFIP2 may also interact closely with some genes in the initiation and progression of different cancers. STRING database was used to construct an integrated PPI network of CYFIP2-interacted genes and performed functional annotations through GO and KEGG analyses. Figure 7A shows a PPI network of CYFIP2 and its top 50 interacted genes. As shown in Fig. 7B, CYFIP2-interacting genes were highly involved in a series of LUAD and immunity related pathways such as Ras signaling pathway, PI3K/AKT signaling pathway, VEGF signaling pathway, chemokine signaling pathway, EGFR tyrosine kinase inhibitor resistance, B cell receptor signaling pathway, T cell receptor signaling pathway, HIF-1 signaling pathway, and non-small cell lung cancer.
Consistent with the KEGG enrichment results, the CYFIP2-interacting genes were highly enriched in a series of important GO functional annotations involved in tumor initiation and progression. In the BP category (Fig. 7C), these genes were significantly enriched in actin filament organization, regulation of actin filament-based process, and actin filament polymerization. The cell leading edge, actin cytoskeleton, lamellipodium, cell-substrate junction, focal adhesion, SCAR complex, neuron spine, dendritic spine, site of polarized growth, and growth cone were the top ten most significantly enriched items in the CC category (Fig. 7D). In terms of MF, the GO terms were highly involved in processes including actin binding, small GTPase binding, and GTPase binding (Fig. 7E).
Furthermore, a significant module with top ten interacting genes was identified (Fig. 8A). We further evaluated the prognostic value of these genes in LUAD. The survival maps indicated that the top ten CYFIP2-interacting genes may also have a predictive role in the OS and DFS of LUAD (Fig. 8B, C). Scatter plots by Spearman’s correlation analysis further presented the correlations between CYFIP2 expression and the top ten CYFIP2-interacting genes based on the TCGA-LUADA cohort (Fig. 8D–M). In LUAD, the expression of CYFIP2 was negatively correlated with the expression of ABI1 (r = − 0.17, P < 0.001), ABI2 (r = − 0.12, P = 0.011), and WASF1 (r = − 0.21, P < 0.001), while the expression of CYFIP2 was positively correlated with the expression of NCKAP1L (r = 0.23, P < 0.001) and WASF3 (r = 0.3, P < 0.001).
CYFIP2 is strongly associated with the infiltration of immune cells
Since CYFIP2 participated in regulating immune-related pathways, we further analyzed the correlations between CYFIP2 expression and 28 types of TILs by the TISIDB database. Figure 9A displays the relationships between expression of CYFIP2 and 28 types of TILs across different cancers. Figure 9B presents that the expression of CYFIP2 was correlated with abundance of act_B cells (r = 0.277, P = 1.8e–10), imm_B cells (r = 0.236, P = 6.17e–8), macrophage cells (r = 0.176, P = 5.83e–5), Th1 cells (r = 0.131, P = 2.83e–3), Tfh cells (r = 0.137, P = 1.83e–8), CD8 + cells (r = 0.207, P = 2.19e-6), pDC cells (r = 0.157, P = 3.4e–4), CD4 + cells (r = 0.152, P = 5.16e–4) and Th17 cells (r = 0.199, P = 5.64e–6).
Using the TIMER database, we further evaluated the relationships between CYFIP2 expression and six main types of tumor infiltrating immune cells of LUAD. As shown in Fig. 9C, CYFIP2 expression had significantly negative associations with tumor purity (r = − 0.149, P = 9.29e-04) and significantly positive correlations with infiltrating levels of B cell (r = 0.388, P = 8.01e–19), CD8 + T cell (r = 0.12, P = 7.93e–03), CD4 + T cell (r = 0.305, P = 6.78e–12), macrophage (r = 0.075, P = 1.00e–1), neutrophil (r = 0.106, P = 1.94e–02), and dendritic cell (r = 0.218, P = 1.22e–06).
Moreover, different algorithms, including CIBERSORT, CIBERSORT-ABS, XCELL, TIMER, EPIC, QUANTISEQ, and MCPCOUNTER, were employed to further investigate the correlations between the above six tumor infiltrating immune cells and CYFIP2 expression in different cancers. As shown in Fig. 10, the results revealed that CYFIP2 expression was positively correlated with B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils, and dendritic cells in most cancer types based on all or most algorithms, especially in LUAD. These results highly suggested that CYFIP2 played a vital part in regulating the infiltration levels of different immune cells.
CYFIP2 is highly involved in the tumor immune microenvironment
For tumor microenvironment analysis, the R package “estimate” was applied to calculate the immune score, stromal score and estimate score based on the CYFIP2 expression level. The expression correlations between CYFIP2 and immune microenvironment scores in different cancers were plotted by heatmap in Fig. 11A. The results revealed that the stromal score, immune score and estimate score of ACC, GBM, KIRC, LGG, SARC, and UCS were negatively associated with CYFIP2 expression levels. In addition, the stromal score, immune score and estimate score of CESC, ESCA, HNSC, LUAD, LUSC, PAAD, SKCM, STAD, and UVM were positively associated with CYFIP2 expression levels. As shown in Fig. 11B-D, all tumor immune microenvironment scores of LUAD were significantly related to the CYFIP2 expression level.
CYFIP2 is significantly related to immune regulators and related genes
As shown in Fig. 12A, B, CYFIP2 was positively associated with most immune stimulators and immune inhibitors in different cancers. Figure 12C, D illustrates the detailed expression correlations between CYFIP2 and immune stimulators/inhibitors. In addition, CYFIP2 expression was also highly related to almost most chemokines, chemokine receptors, and MHC genes in different cancers (Fig. 13A-C). The detailed expression correlations between CYFIP2 and chemokines/chemokine receptors/MHC genes were presented at Fig. 13D–F. These results suggested that CYFIP2 might regulate immune cell infiltration via these genes.
Discussion
CYFIP2, a p53-inducible gene, has been demonstrated by previous studies to be involved in some tumors. Jiao et al. reported CYFIP2 was a novel p53-mediated pro-apoptotic protein whose expression was decreased in gastric cancer and they showed CYFIP2 knockdown promoted proliferation and colony formation, and inhibited apoptosis in gastric cancer cells [29]. Chang’s study revealed that mutational dynamics and genetic variation, as well as aberrant DNA repair, tumor cell cycle control and apoptotic pathways were associated with CYFIP2 in endometrial cancer in the Taiwanese population [30]. Previous studies also showed that CYFIP2 was downregulated in ccRCC patients with an unfavorable survival outcome and might be a potential promising prognostic biomarker associated with immune infiltration, metabolism, and epithelial–mesenchymal transition process in ccRCC [31]. However, the associations between CYFIP2 and different cancers remain inclusive. Besides, the specific mechanisms of CYFIP2 involved in cancer initiation and progression remain to be unclear. The present study first comprehensively uncovered the landscape of CYFIP2 regarding its expression pattern, clinical relevance, genetic alterations, DNA methylation and immune cell infiltration.
We first explored the expression patterns in different levels including normal tissues, cell lines and tumor tissues. Our study demonstrated that CYFIP2 was aberrantly expressed in different cancers. The difference in CYFIP2 expression levels in different tumor types may reflect distinct underlying functions and mechanisms. In addition, CYFIP2 was identified to have the potential for diagnosis of LUAD. We further found that overexpression of CYFIP2 generally predicted good prognosis for patients with tumors, such as LUAD, thymoma, sarcoma, PDAC, clear cell RCC, and BC. In contrast, its highly expression was correlated with poor prognosis in UCEC and ESCC. In addition, a high level of CYFIP2 expression was shown to be related to a favorable prognosis in LUAD with T1&T2 stage, N0 stage, CR & PR diseases. These results suggest that CYFIP2 may serve as a potential biomarker for diagnosis and predicting the prognosis of tumor patients.
Genetic alterations in the key genes have been demonstrated to be highly involved in the initiation and progression of cancer through regulating important signaling pathways that are indispensable for the fundamental intracellular functions including cell growth, survival, invasion, and metastasis [32]. The genetic alteration profiling of CYFIP2 revealed that its mutation is the most important single factors. The impact of mutant subtypes on prognosis and predictive value in LUAD is becoming increasingly understood [33]. Knowledge of genetic alterations will provide more useful information for the diagnosis, treatment and prognosis of cancer.
DNA methylation is a major form of epigenetic modification and is involved in numerous fundamental biological activities through post-transcriptional regulation of gene expression [34]. Aberrant patterns of DNA methylation can lead to the pathogenesis of many diseases including various cancers [35]. Hypermethylation within the promoter region may contribute to silencing or inactivation of tumor suppressor genes in cancer cells [36]. In addition, recent evidence has indicated that DNA methylation has also important roles in the immune system and influence tumor immunotherapy response as DNA methylation is essential for the proper function of B-cell differentiation and mature cells [37]. The present study showed that DNA methylation of CYFIP2 was downregulated in LUAD, which is consistent with the fact that CYFIP2 had a protective role for LUAD patients.
In complex diseases, a specific gene generally plays a coordinating role through its co-expressed and interacting genes. Our results also indicated that CYFIP2 may have a synergistic effect with its co-expressed genes in the prognosis of LUAD. Moreover, our KEGG analysis revealed that CYFIP2-interacting genes were strongly involved in a series of tumor-related and immune-related signaling pathways. For example, CYFIP2-interacting genes were remarkedly linked with non-small cell lung cancer, which was consistent with our present study that CYFIP2 took part in the initiation and progression of LUAD. CYFIP2-interacting genes were also strongly related to Ras signaling, which is a well-studied pathway that is involved in cell proliferation and growth, cell survival and apoptosis, metabolism, and motility and plays an important role in tumor growth and progression [38]. VEGF signaling pathway and EGFR tyrosine kinase inhibitor resistance are two important pathways for the targeted therapy of LUAD [39, 40]. The PI3K/Akt/mTOR pathway has been one of the most altered molecular pathways implicated in both tumorigenesis and the progression of LUAD [41]. Recent evidence also revealed that the PI3K/AKT/mTOR pathway can play a vital role in regulating PD-L1 expression, the abnormal activation of which can contribute to the increased PD-L1 protein translation [42]. In addition, CYFIP2-interacting genes were significantly associated with some immune-related pathways such as chemokine signaling pathway, B cell receptor signaling pathway, and T cell receptor signaling pathway, which play a vital role of tumor immunity and immunotherapy response [43,44,45]. These results confirmed that CYFIP2 played an important role in the tumor development and immune microenvironment.
The immune cells within the TME play important roles in tumorigenesis [46]. It has been known that these tumors associated immune cells may possess tumor-antagonizing or tumor-promoting functions [47]. Tumor-infiltrating immune cells as well as other immune components within the TME can shape an immunosuppressive TME that enables cancer cells to evade surveillance of the immune system [48]. Immune cells are important constituents of the tumor stroma. Studies suggested that the innate immune cells as well as adaptive immune cells contribute to tumor progression when present in the TME [47]. Our study showed the correlation analysis between CYFIP2 expression and immune cell infiltration in pan-cancer. We found that several tumor infiltrating immune cells (Act_B cells, Imm_B cells, macrophage cells, Th1 cells, Tfh cells, CD8 + cells, pDC cells, CD4 + cells and Th17 cells) were correlated with the expression of CYFIP2 in cancers using TISIDB. We further found that positive correlation was indicated between CYFIP2 expression and B cell, CD8 + T cell, CD4 + T cell, macrophage, neutrophil, and dendritic cell in LUAD using TIMER. This finding suggests that there is a potential correlation between CYFIP2 and immune infiltration in LUAD.
Besides that, the CYFIP2 expression presented a significantly correlation with most immune regulators, chemokines, chemokine receptors, and MHC genes. Emerging evidence has supported that immune regulators play a crucial part in establishing an immunosuppressive TME, causing NK cell exhausting and tumor immune escaping [49]. Immune checkpoints blockade can reduce immune escape of tumor cells and limits tumor growth. Thus, targeting immune checkpoints has become one of the most important ways of immunotherapy [50]. Chemokines have been recognized as small, secreted proteins with important roles in directing the migration of immune cells, which is essential for initiating and providing an effective anti-tumor immune response [51]. Each chemokine can activate several different receptors, which also play a predominant part in the coordination of cell trafficking in the immune response processes [52]. MHC genes are clusters of genes encoding molecules that are primarily associated with innate and adaptive immune responses and other molecular and biological processes [53]. These results elucidated the underlying mechanisms of CYFIP2 in tumor immunity and its influence on immune related genes.
Cellular immunotherapy is a hot field in medicine today. From cancer to systemic lupus erythematosus, cellular immunotherapy has shown promise. With clinical progress, cellular immunotherapy is expected to play an even greater role in the clinic. In this study, we suggested that CYFIP2 might be a promising therapeutic strategy for tumor immunotherapy. However, the immunotherapeutic mechanism of CYFIP2 have not been well elucidated. Shabani et al. made the following points that CYFIP2 participated in the vascular endothelial growth factor receptor signaling pathway and the purine ribonucleoside diphosphate metabolic process, respectively [15]. Moreover, Zhao et al. found that CYFIP2 was closely related to T-cell CD8 + , T-cell CD4 + and neutrophils and was highly associated with tumor mutation burden and microsatellite instability in various cancers [54]. In addition, Tong et al. convinced that CYFIP2 was significantly associated with CD4 + cells, CD8 + cells and a series of immune markers via metabolic pathways and epithelial–mesenchymal transition in clear cell renal cell carcinoma [31]. Our study also indicated that CYFIP2 potentially contributed to regulation of tumor-associated and immune-related pathways and regulated immune-related genes. In a word, possible mechanisms through which CYFIP2 might influence tumor immunity are worthy of in-depth study.
There are several limitations in our study. First, most of the present studies were limited to bioinformatics analysis, further study is required to validate these results. Second, the diagnostic and prognostic value of CYFIP2 were evaluated based on TCGA cohorts, which have not been validated by other datasets or clinical data from real-world. In addition, despite the finding indicating that CYFIP2 expression correlated with both immune cell infiltration and patient survival in cancers, we could not prove that CYFIP2 affects patient survival through immune infiltration.
Conclusions
Taken together, our present study indicated that CYFIP2 was aberrantly expressed in various tumors and highly involved in tumor-associated and immune-related pathways. The results also indicated that CYFIP2 could be employed as potential diagnostic and prognostic indicator for numerous tumors. In addition, the abnormal expression of CYFIP2 was linked with immune cells infiltration and may affect tumor immunity by acting on immune regulators and immune-related genes, which may provide a foundation for more precise and personalized immunotherapy in future. However, the mechanisms by which CYFIP2 may affect tumor immunity have not been thoroughly investigated by biological experiments. Thus, further experimental and clinical studies in future are required to promote practical application in immunotherapy and prognosis prediction.
Availability of data and materials
Data are available from the corresponding author upon reasonable request.
Abbreviations
- CYFIP2:
-
Cytoplasmic FMR Interacting Protein 2
- LUAD:
-
Lung adenocarcinoma
- ESCC:
-
Esophageal squamous cell carcinoma
- BC:
-
Bladder carcinoma
- KIRC:
-
Kidney renal clear cell carcinoma
- KIPP:
-
Kidney renal papillary cell carcinoma
- PDAC:
-
Pancreatic ductal adenocarcinoma
- LUSC:
-
Lung squamous cell carcinoma
- HNSC:
-
Head and neck squamous cell carcinoma
- HR:
-
Hazard ratios
- CI:
-
Confidence intervals
References
Barnoud T, Parris JLD, Murphy ME. Common genetic variants in the TP53 pathway and their impact on cancer. J Mol Cell Biol. 2019;11(7):578–85.
Vadapalli S, Abdelhalim H, Zeeshan S, Ahmed Z. Artificial intelligence and machine learning approaches using gene expression and variant data for personalized medicine. Brief Bioinform. 2022. https://doi.org/10.1093/bib/bbac191.
Aveta A, Cilio S, Contieri R, Spena G, Napolitano L, Manfredi C, et al. Urinary MicroRNAs as biomarkers of urological cancers: a systematic review. Int J Mol Sci. 2023. https://doi.org/10.3390/ijms241310846.
Peng Q, Zhang X, Min M, Zou L, Shen P, Zhu Y. The clinical role of microRNA-21 as a promising biomarker in the diagnosis and prognosis of colorectal cancer: a systematic review and meta-analysis. Oncotarget. 2017;8(27):44893–909.
He X, Liu X, Zuo F, Shi H, Jing J. Artificial intelligence-based multi-omics analysis fuels cancer precision medicine. Semin Cancer Biol. 2023;88:187–200.
Lin Y, Qian F, Shen L, Chen F, Chen J, Shen B. Computer-aided biomarker discovery for precision medicine: data resources, models and applications. Brief Bioinform. 2019;20(3):952–75.
Inamura K. Clinicopathological characteristics and mutations driving development of early lung adenocarcinoma: tumor initiation and progression. Int J Mol Sci. 2018. https://doi.org/10.3390/ijms19041259.
Soltis AR, Bateman NW, Liu J, Nguyen T, Franks TJ, Zhang X, et al. Proteogenomic analysis of lung adenocarcinoma reveals tumor heterogeneity, survival determinants, and therapeutically relevant pathways. Cell Rep Med. 2022;3(11): 100819.
Seguin L, Durandy M, Feral CC. Lung adenocarcinoma tumor origin: a guide for personalized medicine. Cancers. 2022. https://doi.org/10.3390/cancers14071759.
Jackson RS 2nd, Cho YJ, Stein S, Liang P. CYFIP2, a direct p53 target, is leptomycin-B sensitive. Cell Cycle. 2007;6(1):95–103.
Manigandan S, Yun JW. Loss of cytoplasmic FMR1-interacting protein 2 (CYFIP2) induces browning in 3T3-L1 adipocytes via repression of GABA-BR and activation of mTORC1. J Cell Biochem. 2022;123(5):863–77.
Biembengut IV, Silva ILZ, Souza T, Shigunov P. Cytoplasmic FMR1 interacting protein (CYFIP) family members and their function in neural development and disorders. Mol Biol Rep. 2021;48(8):6131–43.
Li Y, Song X, Liu L, Yue L. NUAK2 silencing inhibits the proliferation, migration and epithelial-to-mesenchymal transition of cervical cancer cells via upregulating CYFIP2. Mol Med Rep. 2021. https://doi.org/10.3892/mmr.2021.12457.
Lin J, Liao S, Li E, Liu Z, Zheng R, Wu X, et al. circCYFIP2 acts as a sponge of miR-1205 and affects the expression of its target gene E2F1 to regulate gastric cancer metastasis. Mol Ther Nucleic Acids. 2020;21:121–32.
Shabani S, Khayer N, Motalebzadeh J, Majidi Zadeh T, Mahjoubi F. Characterization of pathways involved in colorectal cancer using real-time RT-PCR gene expression data. Gastroenterol Hepatol Bed Bench. 2021;14(2):123–31.
Consortium GT. The genotype-tissue expression (GTEx) project. Nat Genet. 2013;45(6):580–5.
Ponten F, Schwenk JM, Asplund A, Edqvist PH. The human protein atlas as a proteomic resource for biomarker discovery. J Intern Med. 2011;270(5):428–46.
Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(1):509–14.
Hutter C, Zenklusen JC. The cancer genome atlas: creating lasting value beyond its data. Cell. 2018;173(2):283–5.
Lanczky A, Gyorffy B. Web-Based survival analysis tool tailored for medical research (KMplot): development and implementation. J Med Internet Res. 2021;23(7): e27633.
Brunner M, Mullen L, Jauk F, Oliver J, Cayol F, Minata J, et al. Automatic integration of clinical and genetic data using cBioPortal. Stud Health Technol Inform. 2022;290:799–803.
Chandrashekar DS, Karthikeyan SK, Korla PK, Patel H, Shovon AR, Athar M, et al. UALCAN: An update to the integrated cancer data analysis platform. Neoplasia. 2022;25:18–27.
Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63.
Szklarczyk D, Kirsch R, Koutrouli M, Nastou K, Mehryary F, Hachilif R, et al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2023;51(D1):D638–46.
The Gene Ontology C. The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47(1):330–8.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.
Ru B, Wong CN, Tong Y, Zhong JY, Zhong SSW, Wu WC, et al. TISIDB: an integrated repository portal for tumor-immune system interactions. Bioinformatics. 2019;35(20):4200–2.
Goenka A, Khan F, Verma B, Sinha P, Dmello CC, Jogalekar MP, et al. Tumor microenvironment signaling and therapeutics in cancer progression. Cancer Commun. 2023. https://doi.org/10.1002/cac2.12416.
Jiao S, Li N, Cai S, Guo H, Wen Y. Inhibition of CYFIP2 promotes gastric cancer cell proliferation and chemoresistance to 5-fluorouracil through activation of the Akt signaling pathway. Oncol Lett. 2017;13(4):2133–40.
Chang YS, Huang HD, Yeh KT, Chang JG. Identification of novel mutations in endometrial cancer patients by whole-exome sequencing. Int J Oncol. 2017;50(5):1778–84.
Tong J, Meng X, Lv Q, Yuan H, Li W, Xiao W, et al. The downregulation of prognosis- and immune infiltration-related gene CYFIP2 serves as a novel target in ccRCC. Int J Gen Med. 2021;14:6587–99.
Fois SS, Paliogiannis P, Zinellu A, Fois AG, Cossu A, Palmieri G. Molecular epidemiology of the main druggable genetic alterations in non-small cell lung cancer. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms22020612.
Castellanos E, Feld E, Horn L. Driven by mutations: the predictive value of mutation subtype in EGFR-Mutated non-small cell lung cancer. J Thorac Oncol. 2017;12(4):612–23.
Nishiyama A, Nakanishi M. Navigating the DNA methylation landscape of cancer. Trends Genet. 2021;37(11):1012–27.
Oliver J, Garcia-Aranda M, Chaves P, Alba E, Cobo-Dols M, Onieva JL, et al. Emerging noninvasive methylation biomarkers of cancer prognosis and drug response prediction. Semin Cancer Biol. 2022;83:584–95.
Ehrlich M. DNA hypermethylation in disease: mechanisms and clinical relevance. Epigenetics. 2019;14(12):1141–63.
Cao X, Geng Q, Fan D, Wang Q, Wang X, Zhang M, et al. m(6)A methylation: a process reshaping the tumour immune microenvironment and regulating immune evasion. Mol Cancer. 2023;22(1):42.
Punekar SR, Velcheti V, Neel BG, Wong KK. The current state of the art and future trends in RAS-targeted cancer therapies. Nat Rev Clin Oncol. 2022;19(10):637–55.
Zhao Y, Guo S, Deng J, Shen J, Du F, Wu X, et al. VEGF/VEGFR-Targeted therapy and immunotherapy in non-small cell lung cancer: targeting the tumor microenvironment. Int J Biol Sci. 2022;18(9):3845–58.
Wang Q, Zeng A, Zhu M, Song L. Dual inhibition of EGFR-VEGF: An effective approach to the treatment of advanced non-small cell lung cancer with EGFR mutation (Review). Int J Oncol. 2023. https://doi.org/10.3892/ijo.2023.5474.
Sanaei MJ, Razi S, Pourbagheri-Sigaroodi A, Bashash D. The PI3K/Akt/mTOR pathway in lung cancer; oncogenic alterations, therapeutic opportunities, challenges, and a glance at the application of nanoparticles. Transl Oncol. 2022;18: 101364.
Quan Z, Yang Y, Zheng H, Zhan Y, Luo J, Ning Y, et al. Clinical implications of the interaction between PD-1/PD-L1 and PI3K/AKT/mTOR pathway in progression and treatment of non-small cell lung cancer. J Cancer. 2022;13(13):3434–43.
Gupta M, Chandan K, Sarwat M. Natural products and their derivatives as immune check point inhibitors: targeting cytokine/chemokine signalling in cancer. Semin Cancer Biol. 2022;86(Pt 2):214–32.
Lin A, Fang J, Cheng Q, Liu Z, Luo P, Zhang J. B Cell receptor signaling pathway mutation as prognosis predictor of immune checkpoint inhibitors in lung adenocarcinoma by bioinformatic analysis. J Inflamm Res. 2022;15:5541–55.
Morton LT, Wachsmann TLA, Meeuwsen MH, Wouters AK, Remst DFG, van Loenen MM, et al. T cell receptor engineering of primary NK cells to therapeutically target tumors and tumor immune evasion. J Immunother Cancer. 2022. https://doi.org/10.1136/jitc-2021-003715.
de Visser KE, Joyce JA. The evolving tumor microenvironment: from cancer initiation to metastatic outgrowth. Cancer Cell. 2023;41(3):374–403.
Arner EN, Rathmell JC. Metabolic programming and immune suppression in the tumor microenvironment. Cancer Cell. 2023;41(3):421–33.
Li HX, Wang SQ, Lian ZX, Deng SL, Yu K. Relationship between tumor infiltrating immune cells and tumor metastasis and its prognostic value in cancer. Cells. 2022. https://doi.org/10.3390/cells12010064.
Guo Z, Zhang R, Yang AG, Zheng G. Diversity of immune checkpoints in cancer immunotherapy. Front Immunol. 2023;14:1121285.
Ghaedrahmati F, Esmaeil N, Abbaspour M. Targeting immune checkpoints: how to use natural killer cells for fighting against solid tumors. Cancer Commun (Lond). 2023;43(2):177–213.
Ozga AJ, Chow MT, Luster AD. Chemokines and the immune response to cancer. Immunity. 2021;54(5):859–74.
Griffith JW, Sokol CL, Luster AD. Chemokines and chemokine receptors: positioning cells for host defense and immunity. Annu Rev Immunol. 2014;32:659–702.
Sachini N, Papamatheakis J. NF-Y and the immune response: dissecting the complex regulation of MHC genes. Biochim Biophys Acta Gene Regul Mech. 2017;1860(5):537–42.
Zhao Z, He S, Yu X, Lai X, Tang S, Mariya ME, et al. Analysis and experimental validation of rheumatoid arthritis innate immunity gene CYFIP2 and Pan-Cancer. Front Immunol. 2022;13: 954848.
Acknowledgements
Not applicable.
Funding
This work was supported by National Natural Science Foundation of China grants (82003219), Scientific Research Program for Young Talents of China National Nuclear Corporation, Leader Project of Clinical Technology Application Research for Jiangsu Gerontology (LR2021023), Postgraduate Research & Practice Innovation Program of Jiangsu Province (SJCX23_1674), Key Medical Discipline Construction Unit of Jiangsu Province for the 14th Five-year plan (JSDW202236), Suzhou Radiotherapy Clinical Medical Center (Szlcyxzx202103), Nuclear Technology Medical Application Program for Key Talents of The Second Affiliated Hospital of Soochow University (XKTJ-HRC20210013), and Project of State Key Laboratory of Radiation Medicine and Protection, Soochow University (GZK1202021, GZK12023016).
Author information
Authors and Affiliations
Contributions
QP, BR, KX and WL performed the computational analysis, and wrote the manuscript. YY and XG took part in the data analysis. MA polished the language of the manuscript. YZ and YT conceived of the study, and took part in its design and coordination. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not Applicable.
Competing interests
The authors declare that they have no competing interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Peng, Q., Ren, B., Xin, K. et al. CYFIP2 serves as a prognostic biomarker and correlates with tumor immune microenvironment in human cancers. Eur J Med Res 28, 364 (2023). https://doi.org/10.1186/s40001-023-01366-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40001-023-01366-2