Comprehensive analysis of prognostic value and immune infiltration of kindlin family members in non-small cell lung cancer
BMC Medical Genomics volume 14, Article number: 119 (2021)
Kindlin Family Members have been reported to be aberrantly expressed in various human cancer types and involved in tumorigenesis, tumor progression, and chemoresistance. However, their roles in non-small cell lung cancer (NSCLC) remain poorly elucidated.
We analyzed the prognostic value and immune infiltration of Kindlins in NSCLC through Oncomine, GEPIA, UALCAN, CCLE, Kaplan‑Meier plotter, cBioPortal, TIMER, GeneMANIA, STRING, and DAVID database. Additionally, the mRNA expression levels of Kindlins were verified in 30 paired NSCLC tissues and NSCLC cell lines by real-time PCR.
The expression level of FERMT1 was remarkably increased in NSCLC tissues and NSCLC cell lines, while FERMT2 and FERMT3 were reduced. Kindlins expressions were associated with individual cancer stages and nodal metastasis. We also found that higher expression level of FERMT1 was obviously correlated with worse overall survival (OS) in patients with NSCLC, while higher FERMT2 was strongly associated with better overall survival (OS) and first progression (FP). Additionally, the expression of FERMT2 and FERMT3 were obviously correlated with the immune infiltration of diverse immune cells. Functional enrichment analysis has shown that Kindlins may be significantly correlated with intracellular signal transduction, ATP binding and the PI3K-Akt signaling pathway in NSCLC.
The research provides a new perspective on the distinct roles of Kindlins in NSCLC and likely has important implications for future novel biomarkers and therapeutic targets in NSCLC.
Lung cancer is the most common cancer and the leading cause of cancer mortality, causing approximately 1,760,000 deaths worldwide each year . Non-small cell lung cancer (NSCLC) accounts for approximately 80–85% of lung cancers, including lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC) and lung large cell carcinoma (LCC) . Although gradual improvements in diagnosis and treatment therapies, the 5-year overall survival rate of NSCLC remains poor [3, 4]. Therefore, it is necessary to explore the mechanism underlying the tumorigenesis and progression of NSCLC and to search for novel biomarkers with high sensitivity and specificity.
The Kindlin family members are newly discovered focal adhesion proteins consisting of three members (FERMT1, FERMT2 and FERMT3) that share a conserved FERM domain-containing three subdomains (F1, F2 and F3) with an inserted pleckstrin homology domain. Remarkably, Kindlins may function as tumor promoters or suppressors depending on the cancer type. For instance, FERMT2 severs as a tumor-promoting role in gastric cancer , pancreatic cancer [6, 7] and breast cancer [8, 9], but as a tumor-suppressive role in epithelial ovarian cancer  and colorectal cancer . Weinstein et al. reported a higher mRNA expression of FERMT1 in lung cancer tissues . Besides, Zhan et al. also found that FERMT1 and FERMT2 had different expressions in lung cancer. FERMT1 was highly expressed in NSCLC, especially in LUSC. On the other hand, FERMT2 was highly expressed in LCC and weakly expressed in LUAD and LUSC . Additionally, Djaafri et al. found that FERMT3 expression was decreased in lung cancer . However, the series of research only reported the expression levels of Kindlins, and few studies have been conducted on the prognosis and mechanism of Kindlins.
Rapid advances in online platforms and various databases have contributed to the widely used bioinformatics analysis in the field of cancer research. Up to date, bioinformatics analysis has not been adequately used to investigate the roles of Kindlins in NSCLC. Hence, the present study investigated the role of individual kindlins in NSCLC using RT-PCR combined with bioinformatics. In this study, we aim to investigate the prognostic value and immune infiltration of Kindlins in NSCLC, providing new clues for early diagnosis, prognostic judgments and individualized treatments for NSCLC patients.
Materials and methods
Expression database of Kindlins in NSCLC
Oncomine database (http://www.oncomine.org) is a comprehensive cancer microarray database based on 86,733 samples from 715 databases for gene transcriptome analysis in different cancers . The screening and data mining conditions included: “Gene: FERMT1, FERMT2 and FERMT3”; “Analysis Type: Cancer versus Normal Analysis”; “Cancer Type: Lung Cancer”; “Data Type: mRNA”; “P value < 0.05”; “Fold change:2”; “Gene rank: top 10%”. GEPIA database (http://gepia.cancer-pku.cn/) is a recently developed interactive online platform, applying for analyzing the RNA sequencing (RNA-Seq) expression based on over 9000 tumors from the TCGA and 8000 normal samples from the GTEx . GEPIA was used to analyze Kindlins expression and tumor stage in LUAD and LUSC. UALCAN is a user-friendly web that provides publicly available cancer transcriptome data based on TCGA database . UALCAN was used to confirm the association between mRNA expressions of Kindlins in NSCLC and clinicopathological parameters. CCLE (www.broadinstitute.org/ccle) provides the public information on gene expression, chromosome copy number and mutation profile of 947 human cancer cells . In our research, we mainly use it to verify the expression levels of FERMT1, FERMT2 and FERMT3 in different kinds of cancer cell lines.
Kaplan Meier plotter (http://kmplot.com/analysis) aims to discover and validate survival biomarkers based on a meta-analysis from 11 k samples from 20 different cancer types [19, 20]. In the study, we used the KM plotter database from a set of 1926 lung cancer samples to assess Kindlins’ prognostic values, including overall survival (OS), first progression (FP) and post-progression survival (PPS). Moreover, we evaluated the associations of the Kindlins with various clinical parameters of NSCLC, including the histology, clinical stages, gender and smoking history.
cBioPortal (http://cbioportal.org) is an open website used to explore, visualize and analyze multilayer cancer genome data from over 5,000 tumors from 20 different cancer studies [21, 22]. According to cBioPortal's online instructions, the Kindlins gene alterations information in different cancer types was obtained, including genetic mutations, gene fusions, gene amplifications, deep deletions and multiple alterations.
TIMER (https://cistrome.shinyapps.io/timer/) is a comprehensive database for analysis of the tumor-infiltrating immune cells from 32 cancer types . In our research, Spearman correlation was used in the gene module to explore the correlation between the expression of Kindlins and immune infiltration, including tumor purity and six types of immune cells (B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells).
GeneMANIA and STRING database
GeneMANIA (http://www.genemania.org) is a useful and flexible prediction server, displaying a functional interaction network to explore the association between genes and data sets . In the current study, we used GeneMANIA to analyze the relationships between Kindlins concerning the co-expression, co-localization, physical interactions, pathway, genetic interactions, prediction, and shared protein domains.
STRING (https://string-db.org/) is an online database, predicting protein–protein interactions network in terms of direct (physical) and indirect (functional) associations . In this study, we constructed the protein–protein interactions network of Kindlins using and the selection criteria: “organism: Homo sapiens”; “the minimum required interaction score > 0.4”; “the max number of interactors: 20”.
DAVID database (https://david.ncifcrf.gov) is a comprehensive bioinformatics enrichment platform, providing integrative and systematic annotations of biological functions from a series of genes/proteins [26, 27]. We applied the DAVID database (version 6.8) for gene ontology (GO) terms analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis of Kindlins and their related proteins. GO terms covers three aspects: biological processes, cellular components and molecular functions. The P Value < 0.05 was set as a criterion and regarded as significant enrichment.
Lung tissue samples
In this study, between July 2019 and December 2019, thirty pairs of NSCLC tissue and adjacent normal tissues were collected from the Fuzhou Pulmonary Hospital, China. These tissues were used to detect the expression level of Kindlins mRNA by quantitative real-time PCR (RT-PCR). The collection and use of the samples were approved by the ethics committee of the Second Affiliated Hospital of Fujian Medical University. The approval number is 2019 (ethical research review)-207.
The NSCLC cell lines (A549, SPCA-1 and H1299) and one normal cell line (BEAS‐2B) were obtained from Procell life science and Technology Co., Ltd (Wuhan, China). Human bronchial epithelial cells (BEAS-2B) were obtained from FenghBio Co., Ltd. (Changsha, China). All the cell lines were cultured at RPMI-1640 medium (GIBCO, Los Angeles, CA, USA) with 10% fetal bovine serum (FBS, Gibco) and grown in a humidified incubator at 37 °C (5% CO2) environment.
Total RNA in tissues or cells was extracted using TRIzol reagent (Invitrogen) and cDNA was synthesized using Primescript RT Reagent (Takara Bio Inc., Japan). PCR amplification was performed using the SYBR Green PCR kit (Takara Bio Inc., Japan) in a 7500 PCR system (Thermo Fisher Scientific), and GAPDH was used as an endogenous control. The relative quantification analysis was performed using the comparative CT method. The following PCR primers were used:
FERMT1 forward, 5′-TTGAAGATGGTGAGGTTGCGAGTC-3′
FERMT1 reverse, 5′-GGGTTGGCTGAATGCGAGGATG-3′
FERMT2 forward, 5′-TGGCTCTGGACGGGATAAGGATG-3′
FERMT2 reverse, 5′-TTTGTGCTGAGGGGTGAACTGAAG-3′
FERMT3 forward, 5′-ACTGCACCGAGGAGGAGATGATG-3′
FERMT3 reverse, 5′-CCTTGAGGTTGAGCTGCTGAATGG-3′
GAPDH forward, 5′-CTCCTGCACCACCAACTGCTTAG-3′
GAPDH reverse, 5′-GACGCCTGCTTCACCACCTTC-3′
RT-PCR was performed in triplicate. The data were analyzed using the GraphPad Prism (version 8.0). Student’s t-test was used to compare the expression of Kindlins mRNA between NSCLC tissue samples tissues and adjacent normal tissues. P < 0.05 was considered a statistically significant difference.
Differential expression of Kindlins in NSCLC patients and cell lines
We first used ONCOMINE database to explore the expression of Kindlins in NSCLC. Multiple datasets showed that the mRNA expression level of FERMT1 was significantly increased in NSCLC tissues, while FERMT2 and FERMT3 were reduced in NSCLC versus normal tissues (Fig. 1a and Table 1,). These database include the Hou’s dataset , Su dataset , Okayama dataset , Wachi dataset , Selamat’s dataset , Landi’s dataset , and Stearman’s dataset . Moreover, Meta-analysis of Kindlins genes expression in NSCLC studies from Oncomine databases was consistent with the above results (Fig. 1b). Then, we used the GEPIA and UALCAN dataset to further confirm these findings. The results indicated that the expression level of FERMT1 was overexpressed in LUAD and LUSC tissues than in normal tissues, while FERMT2 and FERMT3 were decreased inversely (Fig. 1c–e). We also used CCLE to explore expression levels of Kindlins in NSCLC cell lines. As presented in Fig. 2a, the expression level of FERMT1 was higher in NSCLC cell lines and FERMT3 was downregulated. These results were consistent with those from Oncomine, GEPIA and UALCAN dataset. Different from these results, FERMT2 expression was increased in NSCLC, which requires further examination.
To verify the above results, we detect the mRNA expression level of Kindlins in 30 pairs of NSCLC tissues and NSCLC cell lines using RT-PCR. As expected, the results indicated that the FERMT1 expression was upregulated, while FERMT2 and FERMT3 were down-regulated in the tissue samples from 30 NSCLC cases (P < 0.05, Fig. 2b). Furthermore, the expression level of FERMT1 was highly expressed in NSCLC cell lines (A549, SPCA-1 and H1299) compared with normal cell lines (BEAS‐2B) (P < 0.05, Fig. 2c). For FERMT2 and FERMT3, the mRNA expression levels were reduced in NSCLC cell lines (P < 0.05, Fig. 2c). In conclusion, FERMT1 was significantly up-regulated in NSCLC patients and NSCLC cell lines, and FERMT2 and FERMT3 were significantly down-regulated.
Association between the Kindlins expression and clinicopathological parameters in NSCLC
We then further investigated the associations between the Kindlins mRNA expression and tumor stages in NSCLC using GEPIA and UALCAN datasets. As presented in Fig. 3a, the expression levels of FERMT1 and FERMT2 were remarkably distinct in different tumor stages of LUAD and LUSC, while there was no significance between FERMT3 and different tumor stages. We then explored the expression of Kindlins in distinct individual clinicopathological stages and nodal metastasis status of LUAD and LUSC using UALCAN dataset. It was shown that FERMT1 was highly expressed compared to normal tissues in every clinicopathological stage of LUAD and LUSC, while FERMT2 and FERMT3 genes were downregulated (Fig. 3b). Besides, the expression levels of Kindlins were significantly associated with nodal metastasis status (Fig. 3c). However, the expression levels of the Kindlin family members in different clinicopathological stages and nodal metastasis status were not different. These data suggested that FERMT1, FERMT2 and FERMT3 might play important parts in the tumorigenesis and progression of NSCLC.
Prognostic analysis of Kindlins in patients with NSCLC
Furthermore, we investigated the prognostic significance of Kindlins expression in NSCLC using Kaplan–Meier Plotter. We found that higher FERMT1 expression was significantly correlated with shorter OS in patients with NSCLC (HR = 1.28, log-rank P = 0.00013), while higher FERMT2 expression was strongly associated with better long-term OS (HR = 0.76, log-rank P = 1.4e−05) and FP (HR = 0.82, log-rank P = 0.04) (Fig. 4). Besides, patients with high FERMT3 mRNA levels had longer OS, but the difference was not significant (HR = 0.86, log-rank P = 0.078). Moreover, the expression of FERMT2 was correlated with favorable OS in NSCLC patients with stage 1, suggesting that FERMT2 might act as a prognostic role in early-stage NSCLC (S1). Taken together, the findings point out that FERMT1 expression might be a danger factor, whereas FERMT2 expression might be a protective factor for the prognosis in NSCLC.
Alterations of Kindlins in patients with NSCLC
We then explored the gene alterations of Kindlins in NSCLC using the cBioPortal database. The Kindlins gene alterations were analyzed in 31 cancer studies, which included 10,931 samples. The results showed that gene alterations in Kindlins were present in different types of NSCLC, including LUSC and LUAD, compared with other cancer types (Fig. 5a). Furthermore, alterations frequencies and alterations types of Kindlins were determined in 3897 patients/ 181 samples of 11 NSCLC studies (Fig. 5b). The results indicated that 5.45% of the 3025 cases of NSCLC had mutations, amplifications, deep deletions and multiple changes in Kindlins, the frequencies were 2.25% (68 cases), 2.64% (80 cases), 0.53% (16 cases) and 0.03% (1 case), respectively, and gene amplifications and mutations accounted for the majority. As shown in Fig. 5c, the percentages of gene alterations in individual Kindlins in NSCLC at a range of 1.4–2.5% (FERMT1, 2%; FERMT2, 2.5%; FERMT3, 1.4%;). In addition, it was found that gene amplifications made up the majority of gene alterations for FERMT1 and FERMT2, whereas gene mutations made up the majority of genetic alterations for FERMT3 (Fig. 5c).
Correlation between Kindlins expression and immune infiltration in NSCLC
Since immune cells are closely related to tumor proliferation and progression, we analyzed the correlation between kindlin family members and immune infiltration in LUAD and LUSC using TIMER. We found that FERMT3 expression is significantly negatively associated with tumor purity and remarkably positive correlated with infiltrating levels of B cell, CD8+ T cell, CD4+ T cell, Macrophage, Neutrophil, and Dendritic cell in LUAD and LUSC as presented in Fig. 6c. Additionally, FERMT2 expression showed a similar relationship with tumor purity and the infiltration of CD8+ T cell, CD4+ T cell, Macrophage, Neutrophil, and Dendritic cell in LUAD and LUSC (Fig. 6b). However, the correlation strengths between FERMT1 expression and the six types of immune cells were all relatively weaker (Fig. 6a).
Interaction network of Kindlins at the gene and protein levels
Moreover, the gene–gene interaction network of Kindlins is generated through GeneMANIA. As presented in Fig. 7a, three central nodes implying Kindlins were surrounded by 20 nodes implying genes that were closely associated with Kindlins concerning physical interactions, co-expression, predictions, co-localization, pathway, genetic interactions and shared protein domains. The top five genes showing the greatest connections with the Kindlins contained FBLIM1 (filamin binding LIM protein1), PARVB (Parvin beta), TRIM15 (tripartite motif containing15), DUOXA1 (dual oxidase maturation factor 1) and VCL (vinculin), among which FBLIM1 was linked with FERMT2 for pathway, physical interactions and predictions. Moreover, the functional analysis suggested that these genes were strongly correlated with focal adhesion, integrin activation, cell junction organization and cell-substrate adherens junction.
Using STRING, we also conducted a protein–protein interaction network of the seed genes (FERMT1, FERMT2, and FERMT3) and 20 functional partners with the highest confidence scores (score > 0.900, Fig. 7b). The top five proteins concerning the greatest associations with the Kindlins contained FBLIM1 (filamin binding LIM protein1), CDC5L (cell division cycle 5-like protein), ILK (integrin-linked protein kinase), PICALM (phosphatidylinositol-binding clathrin assembly protein) and ITGB3 (Integrin beta-3). The results indicated that the main biological processes involved the PPI network were platelet degranulation, cell-substrate adhesion, cell adhesion, substrate adhesion-dependent cell spreading, and cell–matrix adhesion, and the main KEGG pathways were focal adhesion, platelet activation, complement and coagulation cascades, amoebiasis and ECM-receptor interaction.
Functional enrichment analysis of Kindlins in NSCLC
Proteomic data associated with differentially expressed Kindlins in NSCLC were harvested by Reverse Phase Protein Array (RPPA) from cBioportal (S2). These data suggest that CTLA4, SQSTM1, CDKN1B, GATA3, FASN, PARK7, BAX, RICTOR, NAPSA, PDK1, MACC1, PEA15, BCL2, GSK3A, GSK3B, TSC2, STAT3, AKT1, AKT2, AKT3, TSC1, AXL, PDCD1, MAPK1, MSH2, RPS6KB1, BID, MTOR, EIF4EBP1, BAP1, ETS1, KEAP1, XRCC5, ERBB3, TP63, CASP3, PXN, FOXO3 and NOTCH (P value < 0.05) were primarily correlated with genomic alterations of Kindlins in NSCLC. Furthermore, enrichment analysis of GO terms and KEGG pathways of these differential proteins was performed using DAVID (S3-6). The top 5 results of GO terms analysis with the terms of biological process analysis were intracellular signal transduction, cell surface receptor signaling pathway, cellular response to chemical stimulus, programmed cell death and cell death (Fig. 8a). For cellular components, Kindlins were significantly enriched in the cytosol, nucleoplasm, mitochondrion, actin cytoskeleton and actin filament (Fig. 8b). Molecular function mainly included ATP binding, protein serine/threonine kinase activity, kinase activity, transcription factor activity and sequence-specific DNA binding (Fig. 8c). Moreover, KEGG pathway enrichment analysis showed that, based on gene count and P value, the PI3K-Akt signaling pathway is regarded to be the most significant for Kindlins in NSCLC. Additionally, KEGG pathway enrichment included pathways in cancer, HIF-1 signaling pathway, insulin signaling pathway and mTOR signaling pathway (Fig. 8d). To sum up, the above results indicated that Kindlins and their related differential proteins were closely linked to tumor-related signaling pathways in NSCLC.
Kindlins belong to the newly discovered focal adhesion proteins. Previous studies have shown that Kindlins were aberrantly expressed in various human cancer types and these dysregulations have involvement in proliferation, epithelial-mesenchymal transition, migration, invasion, and metastasis [35, 36]. However, the complex roles and significance of Kindlins in NSCLC remain unclear. This study for the first time explored the expression patterns, prognostic value, immune infiltration and potential functions of Kindlins in NSCLC. We found Kindlins could play a tumor-promoting or tumor-suppressive role in NSCLC development. Compared with normal tissue, FERMT1 was highly expressed and correlated a worse prognosis in patients with NSCLC, while FERMT2 had an opposite effect. Besides, the expression level of FERMT3 was reduced in patients with NSCLC and higher FERMT3 mRNA levels had longer OS, although the prognosis was not statistically significant. Additionally, a genetic alteration rate of Kindlins (5.9%) was found in patients with NSCLC. These results indicated that Kindlins play critical roles in NSCLC.
FERMT1, also known as Kindlin-1 and URP1, is primarily expressed in epithelial cells . Prior studies have revealed that FERMT1 is overexpressed and acts as a tumor promoter in various cancer types, including breast [38, 39], colon , pancreatic , and hepatocellular  cancer. Mechanistically, several interesting linkages have indicated that FERMT1 may be involved in tumor occurrence and development by regulating the transforming growth factor-beta (TGFβ) signaling [39, 40, 43]. It was reported that FERMT1 was overexpressed in NSCLC and inhibited tumor epithelial-mesenchymal transition, growth and invasion [12, 13]. Adversely, Sin et al. reported that FERMT1 expression was related to a poor prognosis in lung and breast adenocarcinoma, and FERMT1 may be a strategy for inhibiting metastasis . A similar tumorigenic effect of FERMT1 in NSCLC was also shown in our study. The expression of FERMT1 was remarkably increased in patients with NSCLC and NSCLC cell lines, and FERMT1 was significantly associated with tumor grade and nodal metastasis. Moreover, the higher FERMT1 expression was significantly related to poor OS in NSCLC. These results suggested that FERMT1 might act as a potential therapeutic target in NSCLC.
FERMT2, also known as Kindlin-2 and MIG-2, is widely expressed but not in blood cells . The potential role of FERMT2 in various types of cancer has been extensively explored. It has been suggested that FERMT2 acts either as a tumor promoter or suppressor in different cancers. Recent research has demonstrated that FERMT2 is overexpressed in many tumor types, including gastric , breast [8, 9], kidney , glioblastoma , pancreatic [6, 7], bladder , prostate , hepatocellular  and esophageal  cancer, and associated with tumor progression and poor prognosis. In addition, studies have shown that FERMT2 can down-regulate and inhibit tumor cell functions in epithelial ovarian cancer , colorectal cancer  and mesenchymal cancer cell . In lung cancer, the expression pattern of FERMT2 in previous studies is controversial. Weinstein et al. found that the expression of FERMT2 did not observe a significant change in lung cancer . Zhan et al. reported that FERMT2 was overexpressed in LCC, weakly expressed in LUAD and LUSC and promoted tumor invasion and growth in NSCLC . The present study showed that FERMT2 was markedly decreased in NSCLC tissues and NSCLC cell lines. The previous result  was inconsistent with ours, which may be caused by the small sample size and different detection methods in different studies. Additionally, higher FERMT2 expression was strongly associated with long-term OS and FP in NSCLC, and FERMT2 was correlated with tumor grade and nodal metastasis. Therefore, FERMT2 may be used as a prognostic biomarker for NSCLC.
FERMT3, also known as Kindlin-3 and URP2, is mainly expressed in hematopoietic cells . It was reported that FERMT3 was down-regulated and plays a tumor suppressor in many solid tumors including breast cancer, melanoma, as well as lung cancer . Nevertheless, its function as a cancer suppressor or promoter is controversial. Previous research showed that FERMT3 was overexpressed in breast cancer and enhances cancer invasion and metastasis . Moreover, Lu et al. found that FERMT3 was obviously upregulated in glioma, and knockdown of FERMT3 effectively suppressed glioma cell proliferation and chemoresistance . Our study was consistent with the above that FERMT3 expression was significantly decreased in NSCLC, while the expression was not obviously associated with different tumor stages and OS in NSCLC patients. Hence, the prognostic significance of FERMT3 in NSCLC requires further investigation.
The tumor microenvironment (TME) serves an important role in tumor proliferation and progression and is considered as a potential predictor of clinical outcome and immunotherapy reactivity . Our study found that the expression of FERMT2 and FERMT3 was remarkably correlated with the immune infiltration of six immune cells (B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells) in LUAD and LUSC, suggesting that FERMT2 and FERMT3 may be potential biomarkers of NSCLC for immune checkpoint blockade therapy.
Mechanically, GeneMANIA and STRING database showed that Kindlins and their interacting genes and proteins are strongly correlated with integrin activation, cell junction, focal adhesion and cell adhesion. This suggests that Kindlins may be involved in the tumor progression by regulating integrin activation, cell junction, and focal adhesion. Enrichment analysis indicated that the functions of Kindlins and their related differential proteins are mainly associated with the intracellular signal transduction, cell surface receptor signaling pathway, cellular response to chemical stimulus, programmed cell death and cell death in NSCLC. These findings were consistent with previous studies that Kindlins serve important roles as regulators of integrin inside-out signaling . Furthermore, enrichment analysis suggested that Kindlins were primarily correlated with the PI3K-Akt signaling pathway, pathways in cancer, HIF-1 signaling pathway, insulin signaling pathway and mTOR signaling pathway. These data suggest that Kindlins may modulate tumor development and progression by regulating these signaling pathways in NSCLC, emphasizing their potential as anti-NSCLC therapeutic targets.
There were some limitations to this study that need to be noted. First, our study was primarily based on the transcriptional level, whereas protein expressions need further investigation. Additionally, this study didn't carry out experiments on the biological mechanisms of Kindlins in NSCLC. To address these issues, further in vitro and in vivo studies need to validate these findings.
In summary, we analyzed the prognostic value and immune infiltration of individual Kindlins in NSCLC through RT-PCR combined with bioinformatics analysis. The expression level of FERMT1 was significantly elevated in NSCLC tissues and NSCLC cell lines, and the expression levels of FERMT2 and FERMT3 were reduced. We also found that FERMT1 expression was significantly associated with short OS in NSCLC patients, while higher FERMT2 expression was strongly associated with better OS and FP. In addition, FERMT2 and FERMT3 were significantly correlated with immune infiltration. Collectively, the research provides a new perspective on the distinct roles of Kindlins in NSCLC and likely has important implications for future effective therapeutic targets in NSCLC.
Availability of data and materials
All methods were performed in accordance with the relevant guidelines and regulations; Direct web links of datasets about; Oncomine: http://www.oncomine.org; GEPIA database: http://gepia.cancer-pku.cn/; UALCAN: http://ualcan.path.uab.edu; CCLE: www.broadinstitute.org/ccle; Kaplan Meier plotter: http://kmplot.com/analysis; cBioPortal:http://cbioportal.org; TIMER: https://cistrome.shinyapps.io/timer/; GeneMANIA: http://www.genemania.org; STRING: https://string-db.org/; DAVID: https://david.ncifcrf.gov.
Non-small cell lung cancer
Lung squamous cell carcinoma
Lung large cell carcinoma
Gene expression profiling interactive analysis
CBio cancer genomics portal
Search tool for the retrieval of interacting genes
The database for annotation visualization and integrated discovery
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424. https://doi.org/10.3322/caac.21492.
Osmani L, Askin F, Gabrielson E, Li QK. Current WHO guidelines and the critical role of immunohistochemical markers in the subclassification of non-small cell lung carcinoma (NSCLC): moving from targeted therapy to immunotherapy. Semin Cancer Biol. 2018;52:103–9. https://doi.org/10.1016/j.semcancer.2017.11.019.
Duma N, Santana-davila R, Molina JR. Non-small cell lung cancer: epidemiology, screening, diagnosis, and treatment. Mayo Clin Proc. 2019;94:1623–40. https://doi.org/10.1016/j.mayocp.2019.01.013.
Ettinger DS, Wood DE, Aggarwal C, Aisner DL. NCCN guidelines insights: non-small cell lung cancer, version 1.2020. J Natl Compr Cancer Netw. 2019;17:1464–72. https://doi.org/10.6004/jnccn.2019.0059.
Shen Z, Ye Y, Kauttu T, Seppänen H, Vainionpää S, Wang S, et al. Novel focal adhesion protein kindlin-2 promotes the invasion of gastric cancer cells through phosphorylation of integrin β1 and β3. J Surg Oncol. 2013;108:106–12. https://doi.org/10.1002/jso.23353.
Mahawithitwong P, Ohuchida K, Ikenaga N, Fujita H, Zhao M, Kozono S, et al. Kindlin-2 expression in peritumoral stroma is associated with poor prognosis in pancreatic ductal adenocarcinoma. Pancreas. 2013;42:663–9. https://doi.org/10.1097/MPA.0b013e318279bd66.
Yoshida N, Masamune A, Hamada S, Kikuta K, Takikawa T, Motoi F, et al. Kindlin-2 in pancreatic stellate cells promotes the progression of pancreatic cancer. Cancer Lett. 2017;390:103–14. https://doi.org/10.1016/j.canlet.2017.01.008.
Guo B, Gao J, Zhan J, Zhang H. Kindlin-2 interacts with and stabilizes EGFR and is required for EGF-induced breast cancer cell migration. Cancer Lett. 2015;361:271–81. https://doi.org/10.1016/j.canlet.2015.03.011.
Xue X, Li J, Wan W, Shi X, Zheng Y. Kindlin-2 could influence breast nodule elasticity and improve lymph node metastasis in invasive breast cancer. Sci Rep. 2017;7:1–10. https://doi.org/10.1038/s41598-017-07075-1.
Ren C, Du J, Xi C, Yu Y, Hu A, Zhan J, et al. Kindlin-2 inhibits serous epithelial ovarian cancer peritoneal dissemination and predicts patient outcomes. Biochem Biophys Res Commun. 2014;446:187–94. https://doi.org/10.1016/j.bbrc.2014.02.087.
Ren Y, Jin H, Xue Z, Xu Q, Wang S, Zhao G, et al. Kindlin-2 inhibited the growth and migration of colorectal cancer cells. Tumor Biol. 2015;36:4107–14. https://doi.org/10.1007/s13277-015-3044-8.
Weinstein EJ, Bourner M, Head R, Zakeri H, Bauer C, Mazzarella R. URP1: a member of a novel family of PH and FERM domain-containing membrane-associated proteins is significantly over-expressed in lung and colon carcinomas. Biochim Biophys Acta. 2003;1637:207–16. https://doi.org/10.1016/S0925-4439(03)00035-8.
Zhan J, Zhu X, Guo Y, Wang Y, Wang Y, Qiang G, et al. Opposite Role of Kindlin-1 and Kindlin-2 in Lung Cancers. PLoS ONE. 2012;7:1–12. https://doi.org/10.1371/journal.pone.0050313.
Djaafri I, Khayati F, Menashi S, Tost J, Idbaih A, Setterblad N, et al. A novel tumor suppressor function of Kindlin-3 in solid cancer. Oncotarget. 2014;5:8970–85. https://doi.org/10.18632/oncotarget.2125.
Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9:166–80. https://doi.org/10.1593/neo.07112.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45:W98-102. https://doi.org/10.1093/nar/gkx247.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia (United States). 2017;19:649–58. https://doi.org/10.1016/j.neo.2017.05.002.
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. Erratum: Addendum: The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity (Nature (2012) 483 7391 (603–607)). Nature. 2019;565:603–7.
Gyorffy B, Surowiak P, Budczies J, Lánczky A. Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer. PLoS ONE. 2013. https://doi.org/10.1371/journal.pone.0082241.
Nagy Á, Lánczky A, Menyhárt O, Gyorffy B. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Sci Rep. 2018;8:1–9. https://doi.org/10.1038/s41598-018-27521-y.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4. https://doi.org/10.1158/2159-8290.CD-12-0095.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013. https://doi.org/10.1126/scisignal.2004088.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Physiol Behav. 2017;176:139–48.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(SUPPL. 2):214–20. https://doi.org/10.1093/nar/gkq537.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–13. https://doi.org/10.1093/nar/gky1131.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13. https://doi.org/10.1093/nar/gkn923.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. https://doi.org/10.1038/nprot.2008.211.
Hou J, Aerts J, den Hamer B, van Ijcken W, den Bakker M, Riegman P, et al. Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS ONE. 2010;5:e10312. https://doi.org/10.1371/journal.pone.0010312.
Su LJ, Chang CW, Wu YC, Chen KC, Lin CJ, Liang SC, et al. Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme. BMC Genomics. 2007;8:1–12. https://doi.org/10.1186/1471-2164-8-140.
Okayama H, Kohno T, Ishii Y, Shimada Y, Shiraishi K, Iwakawa R, et al. Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 2012;72:100–11. https://doi.org/10.1158/0008-5472.CAN-11-1403.
Wachi S, Yoneda K, Wu R. Interactome-transcriptome analysis reveals the high centrality of genes differentially expressed in lung cancer tissues. Bioinformatics. 2005;21:4205–8. https://doi.org/10.1093/bioinformatics/bti688.
Selamat SA, Chung BS, Girard L, Zhang W, Zhang Y, Campan M, et al. Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. Genome Res. 2012;22:1197–211. https://doi.org/10.1101/gr.132662.111.
Landi MT, Dracheva T, Rotunno M, Figueroa JD, Liu H, Dasgupta A, et al. Gene expression signature of cigarette smoking and its role in lung adenocarcinoma development and survival. PLoS ONE. 2008. https://doi.org/10.1371/journal.pone.0001651.
Stearman RS, Dwyer-Nield L, Zerbe L, Blaine SA, Chan Z, Bunn PA, et al. Analysis of orthologous gene expression between human pulmonary adenocarcinoma and a carcinogen-induced murine model. Am J Pathol. 2005;167:1763–75. https://doi.org/10.1016/S0002-9440(10)61257-6.
Plow EF, Das M, Katarzyna, Bialkowska, Sossey-Alaoui K. Of Kindlins and cancer. Discov. 2016;4:1–15. doi:https://doi.org/10.15190/d.2016.6.Of.
Zhan J, Zhang H. Kindlins: Roles in development and cancer progression. Int J Biochem Cell Biol. 2018;98:93–103. https://doi.org/10.1016/j.biocel.2018.03.008.
Ussar S, Wang HV, Linder S, Fässler R, Moser M. The Kindlins: subcellular localization and expression during murine development. Exp Cell Res. 2006;312:3142–51. https://doi.org/10.1016/j.yexcr.2006.06.030.
Sarvi S, Patel H, Li J, Dodd GL, Creedon H, Muir M, et al. Kindlin-1 promotes pulmonary breast cancer metastasis. Cancer Res. 2018;78:1484–96. https://doi.org/10.1158/0008-5472.CAN-17-1518.
Sin S, Bonin F, Petit V, Meseure D, Lallemand F, Biche I, et al. Role of the focal adhesion protein kindlin-1 in breast cancer growth and lung metastasis. J Natl Cancer Inst. 2011;103:1323–37. https://doi.org/10.1093/jnci/djr290.
Kong J, Du J, Wang Y, Yang M, Gao J, Wei X, et al. Focal adhesion molecule Kindlin-1 mediates activation of TGF-β signaling by interacting with TGF-βRI, SARA and Smad3 in colorectal cancer cells. Oncotarget. 2016;7:76224–37. https://doi.org/10.18632/oncotarget.12779.
Mahawithitwong P, Ohuchida K, Ikenaga N, Fujita H, Zhao M, Kozono S, et al. Kindlin-1 expression is involved in migration and invasion of pancreatic cancer. Int J Oncol. 2013;42:1360–6. https://doi.org/10.3892/ijo.2013.1838.
Ma HX, Shu QH, Pan JJ, Liu D, Xu GL, Li JS, et al. Expression of Kindlin-1 in human hepatocellular carcinoma and its prognostic significance. Tumor Biol. 2015;36:4235–41. https://doi.org/10.1007/s13277-015-3060-8.
Kloeker S, Major MB, Calderwood DA, Ginsberg MH, Jones DA, Beckerle MC. The kindler syndrome protein is regulated by transforming growth factor-β and involved in integrin-mediated adhesion. J Biol Chem. 2004;279:6824–33. https://doi.org/10.1074/jbc.M307978200.
Yan M, Zhang L, Wu Y, Gao L, Yang W, Li J, et al. Increased expression of kindlin-2 is correlated with hematogenous metastasis and poor prognosis in patients with clear cell renal cell carcinoma. FEBS Open Bio. 2016;6:660–5. https://doi.org/10.1002/2211-5463.12063.
Ou Y, Zhao Z, Zhang W, Wu Q, Wu C, Liu X, et al. Kindlin-2 interacts with β-catenin and YB-1 to enhance EGFR transcription during glioma progression. Oncotarget. 2016;7:74872–85. https://doi.org/10.18632/oncotarget.12439.
Wu J, Yu C, Cai L, Lu Y, Jiang L, Liu C, et al. Effects of increased Kindlin-2 expression in bladder cancer stromal fibroblasts. Oncotarget. 2017;8:50692–703. https://doi.org/10.18632/oncotarget.17021.
Yang J, Pan T, Yang H, Wang T, Liu W, Liu B, et al. Kindlin-2 promotes invasiveness of prostate cancer cells via NF-κB-dependent upregulation of matrix metalloproteinases. Gene. 2016;576:571–6. https://doi.org/10.1016/j.gene.2015.11.005.
Lin J, Lin W, Ye Y, Wang L, Chen X, Zang S, et al. Kindlin-2 promotes hepatocellular carcinoma invasion and metastasis by increasing Wnt/β-catenin signaling. J Exp Clin Cancer Res. 2017;36:1–14.
Wang P, Zhan J, Song J, Wang Y, Fang W, Liu Z, et al. Differential expression of Kindlin-1 and Kindlin-2 correlates with esophageal cancer progression and epidemiology. Sci China Life Sci. 2017;60:1214–22. https://doi.org/10.1007/s11427-016-9044-5.
Shi X, Wu C. A suppressive role of mitogen inducible gene-2 in mesenchymal cancer cell invasion. Mol Cancer Res. 2008;6:715–24. https://doi.org/10.1158/1541-7786.MCR-07-2026.
Fagerholm SC, Lek HS, Morrison VL. Kindlin-3 in the immune system. Am J Clin Exp Immunol. 2014;3:37–42.
Sossey-Alaoui K, Pluskota E, Davuluri G, Bialkowska K, Das M, Szpak D, et al. Kindlin-3 enhances breast cancer progression and metastasis by activating Twist-mediated angiogenesis. FASEB J. 2014;28:2260–71. https://doi.org/10.1096/fj.13-244004.
Lu C, Cui C, Liu B, Zou S, Song H, Tian H, et al. FERMT3 contributes to glioblastoma cell proliferation and chemoresistance to temozolomide through integrin mediated Wnt signaling. Neurosci Lett. 2017;657:77–83. https://doi.org/10.1016/j.neulet.2017.07.057.
Papait A, Stefani FR, Cargnoni A, Magatti M, Parolini O, Silini AR. The multifaceted roles of MSCs in the tumor microenvironment: interactions with immune cells and exploitation for therapy. Front Cell Dev Biol. 2020;8:447.
Rognoni E, Ruppert R, Fässler R. The kindlin family: Functions, signaling properties and implications for human disease. J Cell Sci. 2016;129:17–27. https://doi.org/10.1242/jcs.161190.
The present study was supported by The National Key Research and Development Program of China (Grant Number 2016YFC1304000) and Startup Fund for scientific research of Fujian Medical University (Grant Number 2019QH2038).
Ethics approval and consent to participate
The study protocol was approved by the ethics committee of the Second Affiliated Hospital of Fujian Medical University (Quanzhou, China; approval number. 2019-207). Written informed consent was obtained from all participants or their guardians.
Consent for publication
The authors declare that there are no competing interests associated with the study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Table S1: Association of the expression of Kindlins with OS in NSCLC patients with different clinical parameters
. Proteomic data
. GO Biological Process
. GO Cellular Component
. GO Molecular Function
. KEGG pathway
About this article
Cite this article
Su, X., Liu, N., Wu, W. et al. Comprehensive analysis of prognostic value and immune infiltration of kindlin family members in non-small cell lung cancer. BMC Med Genomics 14, 119 (2021). https://doi.org/10.1186/s12920-021-00967-2
- Non-small cell lung cancer