SPTBN1 abrogates renal clear cell carcinoma progression via glycolysis reprogramming in a GPT2-dependent manner
Journal of Translational Medicine volume 20, Article number: 603 (2022)
Renal clear cell carcinoma (ccRCC) is the most prevalent tumors worldwide. Discovering effective biomarkers is essential to monitor the prognosis and provide alternative clinical options. SPTBN1 is implicated in various cancerous processes. However, its role in ccRCC remains unelucidated. This study intends to explore the biological function and mechanism of SPTBN1 in ccRCC.
Single-cell and bulk RNA-seq, tissue microarray, real-time quantitative PCR, and western blotting were applied to verify the expression and predictive value of SPTBN1 in ccRCC. Gain or loss of functional ccRCC cell line models were constructed, and in vitro and in vivo assays were performed to elucidate its tumorigenic phenotypes. Actinomycin D experiment, RNA immunoprecipitation (RIP), specific inhibitors, and rescue experiments were carried out to define the molecular mechanisms.
SPTBN1 was down-regulated in ccRCC and knockdown of SPTBN1 displayed a remarkably oncogenic role both in vitro and in vivo; while overexpressing SPTBN1 reversed this effect. SPTBN1 mediated ccRCC progression via the pathway of glutamate pyruvate transaminase 2 (GPT2)-dependent glycolysis. The expression of GPT2 was significantly negatively correlated with that of SPTBN1. As an RNA binding protein SPTBN1, regulated the mRNA stability of GPT2.
Our research demonstrated that SPTBN1 is significantly down-regulated in ccRCC. SPTBN1 knockdown promotes ccRCC progression via activating GPT2-dependent glycolysis. SPTBN1 may serve as a therapeutic target for the treatment of ccRCC.
Renal cell carcinoma (RCC), a common urinary malignancy worldwide, is estimated to bring 77,410 new cases in China and 79,000 in USA in 2022 [1,2,3]. Among the three subtypes of RCC , ccRCC makes up about 70% . Although treated with various new diagnostic and surgical strategies, around 1/3 of patients still undergo local recurrence or distant metastasis [6, 7]. Targeted drugs, including tyrosine kinase inhibitors (TKIs), and mTOR have exhibited evident effectiveness as first-line treatment for metastatic RCC [8,9,10]. However, the diagnosis is still difficult, and the survival is poor [11,12,13]. New biomarkers and their molecular mechanisms should be discovered to design better treatment options for ccRCC.
The genes in the Spectrin-family were first discovered in erythrocytes in 1968 by Marchesi and Steersand, including erythroid Spectrin and nonerythroid Spectrin . βII Spectrin (SPTBN1) functions as the most common dynamic intracellular protein of non-erythrocyte Spectrin, a cytoskeletal protein from the F-actin cross-linking protein superfamily presenting in all nucleated cells [15, 16]. As a Spectrin-family gene, SPTBN1 is an essential component of conventional βII Spectrin isoforms. SPTBN1 also serves to maintain cellular structure, function, and cycle [17, 18]. Aberrant SPTBN1 expression drives the progression and metastasis of multiple human malignancies. In breast cancer, SPTBN1 down-regulates miR-21 to suppress the epithelial-mesenchymal transition (EMT), thus inhibiting the proliferation and metastasis . The deletion of SPTBN1 activates abnormal Wnt/β-catenin signaling pathway, thereby manifesting the characteristics of liver cancer stem cells . As a novel tumor suppressor, SPTBN1 expression is dysregulated in hepatocellular carcinoma [20,21,22], pancreatic cancer [23, 24], colorectal cancer , ovarian cancer  and lung cancer . SPTBN1 suppresses tumor progression via Wnt/β-catenin, JAK/STAT3 and TGF-β/SMAD pathways [28, 29]. However, the biological function and underlying mechanism of SPTBN1 in ccRCC have never been elucidated.
Metabolic reprogramming is a hallmark of cancer. ccRCC can be characterized as a “metabolic disease”, due to the excessive metabolic activity during its development [30, 31]. To our knowledge, von Hippel-Lindau (VHL) gene harbors the most frequent mutations in ccRCC, while its absence or alteration can stabilize the activities of hypoxia-inducible factor (HIF) family members , thereby activating a large repertoire of genes involved in glucose uptake and glycolytic metabolism pathways . To obtain the energy required for rapid proliferation, tumor cells consumed glucose and lactate through boosting glycolysis, even under aerobic conditions [34, 35]. Glycolysis process involves various enzymes, including SLC2A1, PKM, HKI, HKII, GPT2, and LDHA, all recognized as potential therapeutic targets [36, 37]. Targeting glycolysis with 2-deoxy-D-glucose can inhibit ccRCC proliferation and sensitize ccRCC cells to TKIs . Likewise, overexpression of PFKFB3 enhances glycolysis to facilitate ccRCC proliferation . These research fruits have highlighted the possibility of targeting glycolysis to treat ccRCC.
In the present study, we aimed to investigate the aberrant expression pattern of SPTBN1 in ccRCC and its tumor-suppressive role, as well as its interaction with downstream GPT2 to suppress tumor-required glycolysis. Our data may value SPTBN1 as a biomarker for the early diagnosis and precise treatment of ccRCC.
Materials and methods
Clinical samples collection
ccRCC and adjacent normal renal samples were obtained by radical nephrectomy from the First Affiliated Hospital of Nanjing Medical University (Jiangsu Province Hospital) between 2005 and 2018. Detailed clinicopathological information of involved patients was listed in Table 3. All diagnoses were confirmed by senior pathologists independently. Informed consent was provided by all patients. The study design and protocol were approved by the ethics committee of the First Affiliated Hospital of Nanjing Medical University (Jiangsu Province Hospital).
Tissue microarray (TMA) and immunohistochemistry (IHC)
IHC was performed as previously described [40, 41]. Briefly, the primary antibodies were diluted as follows: anti-SPTBN1 (1:150, Abcam, USA) and anti-GPT2 (1:200, Proteintech, China). TMA was constructed to analyze a total of 180 ccRCC samples. To assess the degree of protein staining expression level, TMA staining signals were quantified using H-score system  (H-Score = ∑(I × Pi), I = intensity of staining, and Pi = percentage of stained tumor cells, ranging from 0 to 300). The expression level of protein was separately and independently assessed by two experienced pathologists blinded to the clinical outcomes. According to the average expression level of IHC, the SPTBN1 in ccRCC patients was considered as low expression (H-score < 70) and high expression (H-score ≥ 70). The GPT2 was considered as low expression (H-score < 60) and high expression (H-score ≥ 60).
Cell culture and transfection
RCC cell lines (786-O, 769-P, ACHN, Caki-1) and human renal tubular epithelial cell line (HK-2) were purchased from the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China) and cultured in RPMI 1640 (786-O, 769-P), McCoy's 5A (Caki-1), DMEM (ACHN) and DMEM/F12 (HK-2) (Gibco, Thermo Fisher Scientific, USA) containing 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin (Gibco, Thermo Fisher Scientific, USA). The lentiviral vectors containing short hairpin RNAs targeting SPTBN1 (shSPTBN1) and negative control (shNC) were constructed and transfected.2-Deoxy-D-glucose (2-DG; Selleck, China), a glycolysis inhibitor, was applied to inhibit aerobic glycolysis. Cells were transfected with small interfering RNA (siRNA) and overexpression SPTBN1 plasmid (oeSPTBN1) using Lipofectamine 3000 (Invitrogen, Thermo Fisher Scientific, USA).
RNA isolation and real-time quantitative PCR (qRT-PCR)
Total RNA was extracted using TRIzol Reagent (Sigma-Aldrich, Merck, Germany). HiScript III All-in-one RT SuperMix (Vazyme, China) was used for cDNA synthesis according to the manufacturer’s instructions. QRT-PCR was performed with SYBR qPCR Master Mix (Vazyme, China) using StepOne Plus (Applied Biosystems, USA) and LightCycler 480 PCR instrument (Roche Diagnostics, Switzerland). The primers and siRNA Oligo used were listed in Additional file 1: Table S1.
Western blotting (WB) assay
Total proteins were extracted using RIPA buffer supplemented with Halt protease and phosphatase inhibitor cocktail (Thermo Fisher Scientific, USA). Then, total 20 μg of protein was resolved using SDS-PAGE gels (4%-12%) and transferred onto polyvinylidene fluoride (PVDF) membranes (Millipore, USA). The membranes were blocked in TBST containing 5% skim milk powder for 2 h at room temperature. The membranes were incubated with the primary antibodies overnight at 4 °C, then with the secondary antibodies conjugated to horseradish peroxidase (HRP) for 1.5 h at room temperature, respectively. Band signals were detected by chemiluminescence system (Bio-Rad, USA).
Cell proliferation and colony formation assay
Pretreated cells were counted and seeded into 96-well plates. Cell proliferation was measured using the CCK-8 Cell Counting Kit (Vazyme, China). The absorbance was measured at 450 nm with a microplate reader, following a 2-h incubation at 37 °C.
For the colony formation assay, pretreated cells were seeded into 6-well plates and incubated for 10–15 days. Colonies were fixed in 4% paraformaldehyde for 20 min and stained with 0.1% crystal violet for further analysis.
Cell cycle determination
Cells were labeled using PI/RNase Staining Buffer (BD Biosciences, USA) according to the manufacturer’s protocol and followed by flow cytometry (BD Biosciences, USA). The results were further analyzed using Modfit software.
Transwell migration and invasion assay
In the migration assay, pretreated cells were seeded into the upper 24-well Transwell chambers with serum-free medium. Matrigel medium (Corning, USA) was applied for invasion assay. The medium containing 20% FBS was added to the bottom chamber. After incubation at 37 °C for 24 or 36 h, the cells were fixed in 4% paraformaldehyde, stained with 0.1% crystal violet, and captured in five randomly selected fields by a microscope. All the experiments were repeated for at least three times.
Tumor xenograft in vivo experiments
All the mouse experiments were approved by the Institutional Animal Care and Use Committee (IACUC) of Nanjing Medical University (No. IACUC-1908003). Briefly, total 2 × 107 786-O cells with stably knockdown-SPTBN1 (shSPTBN1) and negative control cells were collected and suspended with PBS and Matrigel (1:1, Corning, USA), then subcutaneously injected into 4-week-old female BALB/c nude mice. The tumor formation and volume were measured every other day. The formula of tumor volume was calculated as follows: Tumor volume = length * width2 / 2. Xenograft tumor samples were fixed in 4% paraffin for further research.
RNA stability determination
786-O and Caki-1 cells were treated with 5 μg/mL Actinomycin D (Selleck, China) for 0, 2, 4, 6, 8, or 10 h. Total RNA was harvested for qRT-PCR experiments. The level of GPT2 transcript was normalized to that of β-actin control, and the relative half-life of GPT2 was calculated using GraphPad software.
RNA immunoprecipitation (RIP)
RNA immunoprecipitation (RIP) assay was performed with Magna RIP RNA-Binding Protein Immunoprecipitation Kit (Millipore, USA) according to manufacturer’s protocol. Briefly, a total of 2*107 786-O and Caki-1 cells were lysed in RIP lysis buffer, then incubated with magnetic beads conjugated with anti-SPTBN1 or anti-Rabbit IgG antibodies for 6 h at 4 °C temperature. After centrifugation and washing, the supernatant was treated with Proteinase K. Finally, the immunoprecipitated RNA was purified and extracted using TRIzol Reagent for further qRT-PCR analysis.
Library preparation and transcriptome sequencing
786-O shNC and 786-O shSPTBN1 cells were applied for RNA-sequencing in biological duplicates. For each sample, 1.5 μg of RNA was generated for library preparation using NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, USA). The library preparation was accomplished on an Illumina Hiseq-4000 platform (Allwegene, Beijing, China). Raw data in fastq format were processed for quality control and further analyzed using “DESeq” R package [43, 44] (Benjamini and Hochberg’s algorithm). Genes with adjusted P-value (q-value) < 0.05 were considered as differentially expressed.
Dual-luciferase reporter assay
786-O and Caki-1 ccRCC cells were transfected with pEZX-FR02, a luciferase reporter containing wild-type or mutated GPT2-3’UTR region. The ratios of firefly and Renilla luciferase activities were calculated after 48 h of transfection using Dual-Luciferase Assay System (Promega, USA).
Glucose consumption, lactate secretion, and ATP concentration measurement
Assay kits (Nanjing Jiancheng Bioengineering, China) were used to measure the glucose concentrations, lactic acid production, and ATP level in 786-O and Caki-1 ccRCC cell lines. Concentration was calculated and normalized to cell number.
Bioinformatics analysis of in silico database
TCGA-KIRC cohort, GSE6344, GSE40435, GSE46699, GSE53757, GSE66270, GSE105261, together with single-cell RNA-seq GSE152938 and GSE156632 datasets were employed in this study [45,46,47,48,49,50,51,52,53,54]. Cancer Cell Line Encyclopedia (CCLE), CPTAC and UALCAN databases were also applied [55,56,57]. The proteogenomic expression data in Chinese ccRCC patients were downloaded from the supplementary materials of Qu . For single-cell RNA-seq analysis, R package “Seurat” was used to process data [54, 59,60,61]. The cells expressing > 20% mitochondria-related genes, > 5% of hemoglobin reads, less than 500 genes, or greater than 10,000 genes were filtered out, then we normalized and rescaled the expression matrix. Next, cell clusters were obtained through the UMAP method. Finally, we used the “SingleR” package to annotate the cell clusters .
Statistical analysis was run on SPSS and R software. Data were expressed with the means ± SD. The independent t-test was used to compare continuous variables that exhibited normal distributions. The Wilcoxon test was used to compare the continuous variables that were not normally distributed. Kaplan–Meier survival analysis and Cox hazard regression model were employed to sift out factors predicting the overall survival, disease-specific survival, disease-free interval, and progression-free interval prognostic factors. Correlations analysis was performed by Pearson and Spearman methods. All experiments were repeated at least three times. All statistical tests were two-sided, and P < 0.05 was considered statistically significant. ns: No significant; *:P < 0.05; **:P < 0.01; ***:P < 0.001; ****:P < 0.0001.
Expression profile of Spectrin-family genes in ccRCC
Seven Spectrin-family genes with significantly abnormal expression profiles in ccRCC tissues were screened out from the TCGA-KIRC, along with the clinicopathological parameters of related patients [15, 16, 28, 29]. The flowchart of screening process was illustrated in Fig. 1A. Our results confirmed that SPTA1 and SPTBN5 were up-regulated in ccRCC tissues, while SPTAN1, SPTB, SPTBN1, SPTBN2, and SPTBN4 were down-regulated in ccRCC tissues (Fig. 1B, C). Meanwhile, we analyzed the distribution of Spectrin-family genes using datasets of single-cell RNA-sequencing (Additional file 1: Figure S1A-E). SPTBN1 was mainly enriched in endothelial cells and tissue stem cells, and SPTBN2 in endothelial cells and natural killer cells (Fig. 1D, E). Then, we conducted univariate and multivariate cox regression analyses to assess the prognostic value of Spectrin-family genes. Among these candidate genes, only SPTBN1 and SPTBN2 were significantly associated with overall survival (OS), disease-specific survival (DSS), and progression survival interval (PFI) (Additional file 1: Figure S2A–C). In addition, the correlations between Spectrin-family genes were analyzed (Additional file 1: Figure S3A). Functional enrichment analysis demonstrated that Spectrin-family genes were mainly involved in the interaction between L1 and Ankyrins (Additional file 1: Figure S3B). After targeting small interfering RNA to SPTBN1 and SPTBN2, we identified SPTBN1 as a candidate target in ccRCC among Spectrin-family genes (Fig. 1F). Finally, our pan-cancer analysis discovered that SPTBN1 was significantly down-regulated in the tissue samples of most cancer subtypes, compared to that in adjacent normal samples (Additional file 1: Figure S4A). Based on the Cancer Cell Line Encyclopedia (CCLE) dataset, we observed that SPTBN1 was highly expressed in kidney cancer, compared to that in most of other solid tumors  (Additional file 1: Figure S4B). These findings indicated that SPTBN1 might act as a tumor suppressor role in ccRCC.
SPTBN1 expression was down-regulated in ccRCC
To further investigate the tumor-suppressing role of SPTBN1 in ccRCC, firstly we detected the expression of SPTBN1 in multiple databases. The results showed that SPTBN1 was significantly decreased in both paired and unpaired samples from the TCGA-KIRC cohort (P < 0.001; Fig. 1G). This expression patterns was further validated in external GEO databases, including GSE66270, GSE105261, GSE6344, GSE40435, GSE46699 and GSE53757. (Additional file 1: Figure S5A–F). CPTAC and FUSCC proteome databases also demonstrated that the protein level of SPTBN1 was downregulated in ccRCC tissues (P < 0.001; Fig. 1H).
We further explored the expression of SPTBN1 in clinical ccRCC samples by IHC, qRT-PCR and WB. IHC staining results confirmed that the expression of SPTBN1 was remarkably decreased in tumor tissues (Additional file 1: Figure S5G). Besides, the receiver operating characteristic (ROC) curve analysis showed a high diagnostic value of SPTBN1 for ccRCC (AUC:0.692; 95% CI = 0.644–0.741; Fig. 1I). Moreover, qRT-PCR and WB also demonstrated the same results (NJMU ccRCC cohort 1: N = 53; Fig. 1J, K). In accordance with the findings in TCGA and GEO datasets, SPTBN1 expression was significantly down-regulated in ccRCC cell lines (786-O, 769-P, ACHN, and Caki-1) at both mRNA and protein levels, compared with that in renal tubular normal epithelial cell line HK-2 (P < 0.05; Fig. 1L).
SPTBN1 down-regulation predicted an unfavorable survival of ccRCC patients
We observed that low expression of SPTBN1 was positively correlated with high age, TNM stage, and tumor grade (P < 0.05, Fig. 3A; Table 1). Meanwhile, the analysis based on CPTAC and FUSCC proteome databases revealed low SPTBN1 expression in high-grade ccRCC patients (P < 0.05; Additional file 1: Figure S6A, B). In univariate and multivariate cox regression analyses, SPTBN1 served as an independent predictive factor for the survival of ccRCC patients (Univariate: HR = 0.532, 95% CI = 0.444 − 0.637, P < 0.001; Multivariate: HR = 0.647, 95% CI = 0.489–0.854, P = 0.002; Fig. 2B). We also examined the correlation between SPTBN1 and survival in subgroups of different clinicopathological characteristics, discovering that a low expression of SPTBN1 might predict a poor prognosis (Additional file 1: Figure S7A–D). The logistic regression analysis was adopted to detail the correlativity between SPTBN1 expression and clinicopathological characteristics (Table 2). Moreover, a nomogram based on TNM stage and SPTBN1 was developed to predict the 1-, 3-, and 5-year OS of each ccRCC patient (Fig. 2C). Kaplan–Meier analysis further indicated that a low SPTBN1 expression was correlated with worse overall survival, disease-specific survival, and progression-free survival in the TCGA-KIRC cohort (Fig. 2D). According to SPTBN1 expression in our tissue microarray from NJMU TMA cohort (N = 180), we divided the patients into high-SPTBN1 and low-SPTBN1 expression subgroups (Fig. 3E; Table 3). The results showed that a low SPTBN1 expression was associated with poor overall survival and relapse-free survival (P < 0.01; Fig. 3F). Simultaneously, we observed that SPTBN1 was predominantly down-regulated in ccRCC patients with advanced clinicopathological characteristics, which was consistent with the results of TCGA-based analysis (Table 3). Taken together, SPTBN1 possessed a high prognostic value for ccRCC patients.
SPTBN1 inhibited the proliferation of ccRCC cells and malignant potential in vitro and in vivo
To further investigate the biological function of SPTBN1 in ccRCC, we successfully established SPTBN1-knockdown and SPTBN1-overexpression cell models (786-O and Caki-1) and validated by qRT-PCR and WB (Additional file 1: Figure S8A, B, P < 0.05). Cell counting kit-8 (CCK-8) assay indicated that SPTBN1 knockdown significantly increased cell proliferation ability (Fig. 3A). Colony formation assay was also employed to determine the long-term impact of SPTBN1 on cells proliferation. Colony formation was more evident in SPTBN1 knockdown group than in control group, while SPTBN1 overexpression reversed this change (Fig. 3B). Furthermore, flow cytometry showed that SPTBN1 knockdown increased the proportion of 786-O and Caki-1 cells in the S phase and decreased that in the G0-G1 phase, while SPTBN1 overexpression resulted in opposite effect (Fig. 3C). Additionally, Transwell migration and invasion assay demonstrated that SPTBN1 knockdown turned CCRCC cells more aggressive and metastatic (Additional file 1: Figure S9A).
In addition, to determine the influence of SPTBN1 in vivo, we established subcutaneously xenograft tumor models using nude mice (Fig. 3D). Knocking down SPTBN1 dramatically accelerated tumor growth, as quantified by tumor size and tumor volume (Fig. 3D). Accordingly, IHC experiments further revealed that Ki-67 and PCNA, indicators of proliferation, were notably up-regulated in tumors derived from shSPTBN1 786-O cells, indicating more rapid tumor growth (Fig. 3E). Collectively, SPTBN1 could retard the G1/S progression and suppress the proliferation of ccRCC cells both in vitro and in vivo.
SPTBN1 deficiency reprogrammed glycolysis in ccRCC cells
In order to clarify the underlying mechanism of SPTBN1 involved in RCC cell proliferation, we performed RNA-sequencing on shNC and shSPTBN1 786-O cells. A total of 4491 up-regulated genes and 4179 down-regulated genes were identified (Fig. 4A). Gene ontology (GO) functional enrichment analysis revealed their enrichment in metabolic process-related pathways (Fig. 4B). Subsequently, we performed gene set enrichment analysis (GSEA) using both TCGA-KIRC and our RNA-seq dataset. Interestingly, Reactome Glucose Metabolism and Reactome Glycolysis pathways were significantly associated with the biological function of SPTBN1 (Fig. 4C). To verify this hypothesis, we detected the expression level of rate-limiting enzymes involved in glycolytic process. Unexpectedly, glycolysis-associated genes (PKM, HKII, PFKFB3, ENO2, SLC2A1, PFKP, and LDHA) were dramatically up-regulated after SPTBN1 knockdown, whereas SPTBN1 overexpression reversed this influence (Fig. 4D–G). Consequently, we observed that knockdown of SPTBN1 significantly enhanced lactate concentrations (Fig. 4H), glucose consumption rates (Fig. 4I) and ATP levels (Fig. 4J) in ccRCC cell lines, while SPTBN1 overexpression substantially reversed this increase (Fig. 4H–J).
Ultimately, supplementation with 2 mM 2-deoxy-D-glucose (2-DG), a specific inhibitor of glycolysis, could partially block glycolysis-induced proliferation of ccRCC cells. SPTBN1 knockdown dramatically prevented 2-DG from suppressing 786-O and Caki-1 cells. The half-maximal inhibitory concentration (IC50) data revealed that SPTBN1 knockdown increased the IC50 value and weakened the sensitivity to 2-DG, but SPTBN1 overexpression displayed contradictory effects (Fig. 4K). Taken together, that 2-DG could partially repress the malignant proliferation of SPTBN1-knockdown ccRCC cells (Fig. 4L).
SPTBN1 was negatively correlated with GPT2 expression
We next explored the downstream genes of SPTBN1 in the glycolytic pathway. We screened out the GPT2 gene based on log2FC (Fold Change) and P-value (Fig. 5A). qRT-PCR experiments further confirmed that GPT2 expression was significantly elevated in SPTBN1-knockdown 786-O and Caki-1 cells, while GPT2 mRNA levels fell significantly after SPTBN1 overexpression (Fig. 5B, C). The correlation analysis in the TCGA-KIRC database indicated that the expression levels of SPTBN1 and GPT2 were negatively correlated (Pearson: R = − 0.176, P < 0.001; Spearman: R = − 0.203, P < 0.001; Fig. 5D). In addition, the same results were observed in NJMU ccRCC cohort (N = 53; Pearson: R = − 0.703, P < 0.001; Spearman: R = − 0.668, P < 0.001; Fig. 5E). TMA immunohistochemical data also revealed a negative correlation between SPTBN1 and GPT2 (NJMU ccRCC TMA cohort2: N = 180; Pearson: R = − 0.274, P < 0.001; Spearman: R = − 0.292, P < 0.001; Fig. 5F, G).
GPT2 up-regulation altered the oncogenicity in ccRCC cells
Subsequently, we evaluated the expression and prognostic value of GPT2 in ccRCC patients. qRT-PCR confirmed that GPT2 mRNA was obviously up-regulated in NJMU ccRCC cohort (N = 53, P < 0.001; Fig. 6A). Likewise, GPT2 up-regulation was observed in IHC (Additional file 1: Figure S10A). Meanwhile, qRT-PCR and WB showed that GPT2 expression increased significantly in various types of CCRCC cells, compared with that in HK2 cells (Fig. 6B). The ROC curve analysis showed an AUC of 0.816 (95%CI = 0.770–0.861; Fig. 6C). Kaplan–Meier survival analysis indicated that the prognosis (OS, DSS, and PFI) of patients with high GPT2 expression was significantly worse (OS: P = 0.014; DSS: P = 0.022; PFI: P = 0.003; Fig. 6D), demonstrating the excellent prognostic value of GPT2. According to our tissue microarray from NJMU TMA cohort (N = 180), we divided patients into high- and low-GPT2 expression subgroups. The results showed that high GPT2 expression predicted longer overall survival and relapse-free survival (P < 0.01; Fig. 6E). CCK-8 assay demonstrated that GPT2 knockdown significantly decreased cell proliferation (Fig. 6F). Colony-formation efficiency was reduced after GPT2 depletion (Fig. 6G). GPT2 knockdown significantly suppressed lactate concentration, glucose consumption rate and ATP level in ccRCC cell lines (Fig. 6H, I). In conclusion, GPT2, as a downstream gene of SPTBN1, promoted ccRCC progression via glycolysis.
SPTBN1 regulated the stability of GPT2 mRNA by binding to its 3’-UTR regions
Considering the interaction between SPTBN1 and GPT2 we ascertained above, we suspected whether depletion of SPTBN1 could alter GPT2 expression. The mRNA expression level of SPTBN1 was negatively correlated with that of GPT2, indicating SPTBN1 regulated GPT2 at the transcriptional level. Given that SPTBN1 is known as an RNA-binding protein, we hypothesized that SPTBN1 functions as an RBP binding to the 3’UTR region of GPT2 [64, 65]. Control and SPTBN1-deficient 786-O and Caki-1 cells were stimulated with actinomycin D for 0, 2, 4, 6, 8, and 10 h. Then, mRNA stability was monitored using qRT-PCR. This analysis revealed that SPTBN1-depleted cells predominantly contained more stable GPT2 mRNAs than control cells, whereas SPTBN1 overexpression suppressed GPT2 mRNA level in 786-O and Caki-1 cells (Fig. 7A). Congruently, RIP experiments in 786-O and Caki-1 (with SPTBN1 overexpression) cells further demonstrated that SPTBN1 could bind to GPT2 (Fig. 7B). In the dual-luciferase reporter assay, 786-O and Caki-1 CCRCC cells were transfected with pEZX-FR02 plasmid which contained mutated GPT2 3’-UTR region (GPT2-MUT) and wild type of GPT2 3’-UTR region (GPT2-WT). Compared with mutated GPT2 3’-UTR reporter plasmids, SPTBN1 could significantly reduce the luciferase activity of the wild-type GPT2 3’-UTR, which was in concordance with our hypothesis (Fig. 7C). These findings confirmed that SPTBN1 could function as an RNA-binding protein to regulate the stability of GPT2 mRNA by binding to its 3’-UTR regions.
Silencing GPT2 can reversed the malignant proliferation of ccRCC cells with SPTBN1 deletion
We investigated whether silencing SPTBN1 can reverse the malignant proliferation of ccRCC cells through directly modulating GPT2. CCK-8 assay revealed that GPT2 knockdown attenuated the oncogenic effect of SPTBN1 knockdown (Fig. 7D). Likewise, GPT2 knockdown abolished the ability of SPTBN1 in promoting colony formation (Fig. 7E). It was further validated that GPT2 knockdown could partly reduce ATP concentration, lactate production, and glucose uptake, implying that the glycolysis caused by SPTBN1 depletion in ccRCC cells (Fig. 7F, G). Collectively, SPTBN1 suppressed ccRCC development through degrading GPT2 to reprogram glyolysis (Fig. 8).
The prevalence of RCC has kept increasing worldwide in recent years [1, 3]. ccRCC is a typically heterogeneous disease with various molecular subtypes [66, 67]. scRNA-seq has proven as a potent tool for analyzing the heterogeneity of a tumor and the status of individual cells within tumor microenvironments [68,69,70,71]. In the present research, we for the first-time combined bulk RNA-seq and scRNA-seq to unveil the expression and distribution of Spectrin-family genes in ccRCC. We also teased out SPTBN1 as a novel ccRCC suppressor, as reflected by that its low expression could predict a poor survival and its overexpression could abrogate tumorigenesis.
Metabolic reprogramming is a hallmark as a tumor grows malignant . Malignant cells can activate metabolic pathways, thus providing them nutrition enough to thrive and proliferate, despite that the patient it is taking dietary restrictions [72, 73]. ccRCC is characterized by the dysregulation of metabolic pathways, especially glycolytic pathway . However, the molecular mechanism underlying metabolic reprogramming remains to be further explored. Glycolysis is closely involved in cancer proliferation, metabolism, migration and many another malignant processes [74,75,76]. Here, we found that SPTBN1 suppressed ccRCC progression through blocking glycolysis-related pathways involving PKM1/2, LDHA, HKII, and PFKP.
We further validated the mechanism of SPTBN1 in ccRCC glycolysis. Previous studies have demonstrated that SPTBN1 regulates SOCS3 at the transcriptional level through the TGF-β/Smads2/3 pathway . Lin et al. have proven that SPTBN1 can induce the ubiquitination and degradation of p65 by means of SOCS1. To investigate the mechanism by which SPTBN1 regulates GPT2, we used the online database of “RNAct” and confirmed SPTBN1 as an RNA-binding protein (RBP) [64, 65]. Then, we selected GPT2 as a direct downstream target of SPTBN1. Our research revealed that SPTBN1 could function as an RBP to regulate the mRNA stability of target genes. Based on our RNA-seq data, we found that SPTBN1 could bind to the 3’UTR in GPT2 mRNA to reduce its stability.
Glutamic pyruvate transaminase (GPT) participates in numerous cellular processes, primarily through transamination between alanine and α-ketoglutarate (α-KG) to generate pyruvate and glutamate [77, 78]. Several studies have discovered the implication of GPT2 in glucose metabolism and homeostasis . For instance, GPT2 promotes tumorigenesis of breast cancer by activating sonic Hedgehog signaling , and abrogating GPT2 in triple-negative breast cancer can inhibit tumor growth and promotes autophagy . Our findings elucidated that SPTBN1 interacted with GPT2 in a negative manner in both ccRCC tissues and cell lines, and that GPT2 could promote ccRCC cells proliferation in vitro via boosting glycolysis. These suggest GPT2 as a direct downstream molecule of SPTBN1 in the glycolytic pathway. Silencing GPT2 could partly reverse the malignant proliferation and the activation of glycolysis pathway caused by SPTBN1 deletion, further supporting the role of SPTBN1-GPT2-mediated glycolysis in promoting ccRCC oncogenesis. However, our findings need more ccRCC patient cohorts to verify its feasibility and constancy.
SPTBN1 is significantly down-regulated in ccRCC. SPTBN1 knockdown promotes ccRCC progression via activating GPT2-dependent glycolysis. SPTBN1 may serve as a therapeutic target for the treatment of ccRCC.
Availability of data and materials
All data generated or analyzed during this study were included either in this article methods section or in the Supplemental Information. Other data that support the findings of this study are available from the corresponding author upon reasonable request.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2022. CA Cancer J Clin. 2022;72(1):7–33.
Cohen HT, McGovern FJ. Renal-cell carcinoma. N Engl J Med. 2005;353(23):2477–90.
Xia C, Dong X, Li H, Cao M, Sun D, He S, Yang F, Yan X, Zhang S, Li N, et al. Cancer statistics in China and United States, 2022: profiles, trends, and determinants. Chin Med J (Engl). 2022;135(5):584–90.
Linehan WM, Ricketts CJ. The Cancer Genome Atlas of renal cell carcinoma: findings and clinical implications. Nat Rev Urol. 2019;16(9):539–52.
Jonasch E, Walker CL, Rathmell WK. Clear cell renal cell carcinoma ontogeny and mechanisms of lethality. Nat Rev Nephrol. 2021;17(4):245–61.
Jonasch E, Gao J, Rathmell WK. Renal cell carcinoma. BMJ. 2014;349: g4797.
Staehler M, Rohrmann K, Haseke N, Stief CG, Siebels M. Targeted agents for the treatment of advanced renal cell carcinoma. Curr Drug Targets. 2005;6(7):835–46.
Rini BI, Campbell SC, Escudier B. Renal cell carcinoma. Lancet. 2009;373(9669):1119–32.
Kotecha RR, Motzer RJ, Voss MH. Towards individualized therapy for metastatic renal cell carcinoma. Nat Rev Clin Oncol. 2019;16(10):621–33.
Posadas EM, Limvorasak S, Figlin RA. Targeted therapies for renal cell carcinoma. Nat Rev Nephrol. 2017;13(8):496–511.
Duran I, Lambea J, Maroto P, González-Larriba JL, Flores L, Granados-Principal S, Graupera M, Sáez B, Vivancos A, Casanovas O. Resistance to targeted therapies in renal cancer: the importance of changing the mechanism of action. Target Oncol. 2017;12(1):19–35.
Makhov P, Sohn JA, Serebriiskii IG, Fazliyeva R, Khazak V, Boumber Y, Uzzo RG, Kolenko VM. CRISPR/Cas9 genome-wide loss-of-function screening identifies druggable cellular factors involved in sunitinib resistance in renal cell carcinoma. Br J Cancer. 2020;123(12):1749–56.
Zhao T, Bao Y, Gan X, Wang J, Chen Q, Dai Z, Liu B, Wang A, Sun S, Yang F, et al. DNA methylation-regulated QPCT promotes sunitinib resistance by increasing HRAS stability in renal cell carcinoma. Theranostics. 2019;9(21):6175–90.
Marchesi VT, Steers E. Selective solubilization of a protein component of the red cell membrane. Science. 1968;159(3811):203–4.
Machnicka B, Czogalla A, Hryniewicz-Jankowska A, Bogusławska DM, Grochowalska R, Heger E, Sikorski AF. Spectrins: a structural platform for stabilization and activation of membrane channels, receptors and transporters. Biochim Biophys Acta. 2014;1838(2):620–34.
Zhang R, Zhang C, Zhao Q, Li D. Spectrin: structure, function and disease. Sci China Life Sci. 2013;56(12):1076–85.
Derbala MH, Guo AS, Mohler PJ, Smith SA. The role of βII spectrin in cardiac health and disease. Life Sci. 2018;192:278–85.
Liu Y, Qi J, Chen X, Tang M, Chu C, Zhu W, Li H, Tian C, Yang G, Zhong C, et al. Critical role of spectrin in hearing development and deafness. Sci Adv. 2019;5(4):eaav803.
Wu H, Chen S, Liu C, Li J, Wei X, Jia M, Guo J, Jin J, Meng D, Zhi X. SPTBN1 inhibits growth and epithelial-mesenchymal transition in breast cancer by downregulating miR-21. Eur J Pharmacol. 2021;909: 174401.
Zhi X, Lin L, Yang S, Bhuvaneshwar K, Wang H, Gusev Y, Lee M-H, Kallakury B, Shivapurkar N, Cahn K, et al. βII-Spectrin (SPTBN1) suppresses progression of hepatocellular carcinoma and Wnt signaling by regulation of Wnt inhibitor kallistatin. Hepatology. 2015;61(2):598–612.
Rao S, Yang X, Ohshiro K, Zaidi S, Wang Z, Shetty K, Xiang X, Hassan MI, Mohammad T, Latham PS, et al. β2-spectrin (SPTBN1) as a therapeutic target for diet-induced liver disease and preventing cancer development. Sci Transl Med. 2021;13(624):eabk267.
Lin L, Chen S, Wang H, Gao B, Kallakury B, Bhuvaneshwar K, Cahn K, Gusev Y, Wang X, Wu Y, et al. SPTBN1 inhibits inflammatory responses and hepatocarcinogenesis via the stabilization of SOCS1 and downregulation of p65 in hepatocellular carcinoma. Theranostics. 2021;11(9):4232–50.
Huang T-C, Renuse S, Pinto S, Kumar P, Yang Y, Chaerkady R, Godsey B, Mendell JT, Halushka MK, Civin CI, et al. Identification of miR-145 targets through an integrated omics analysis. Mol Biosyst. 2015;11(1):197–207.
Jiang X, Gillen S, Esposito I, Giese NA, Michalski CW, Friess H, Kleeff J. Reduced expression of the membrane skeleton protein beta1-spectrin (SPTBN1) is associated with worsened prognosis in pancreatic cancer. Histol Histopathol. 2010;25(12):1497–506.
Ying J, Lin C, Wu J, Guo L, Qiu T, Ling Y, Shan L, Zhou H, Zhao D, Wang J, et al. Anaplastic lymphoma kinase rearrangement in digestive tract cancer: implication for targeted therapy in chinese population. PLoS ONE. 2015;10(12): e0144731.
Chen M, Zeng J, Chen S, Li J, Wu H, Dong X, Lei Y, Zhi X, Yao L. SPTBN1 suppresses the progression of epithelial ovarian cancer via SOCS3-mediated blockade of the JAK/STAT3 signaling pathway. Aging (Albany NY). 2020;12(11):10896–911.
Gu F-F, Zhang Y, Liu Y-Y, Hong X-H, Liang J-Y, Tong F, Yang J-S, Liu L. Lung adenocarcinoma harboring concomitant SPTBN1-ALK fusion, c-Met overexpression, and HER-2 amplification with inherent resistance to crizotinib, chemotherapy, and radiotherapy. J Hematol Oncol. 2016;9(1):66.
Yang P, Yang Y, Sun P, Tian Y, Gao F, Wang C, Zong T, Li M, Zhang Y, Yu T, et al. βII spectrin (SPTBN1): biological function and clinical potential in cancer and other diseases. Int J Biol Sci. 2021;17(1):32–49.
Chen S, Li J, Zhou P, Zhi X. SPTBN1 and cancer, which links? J Cell Physiol. 2020;235(1):17–25.
Hu C-J, Wang L-Y, Chodosh LA, Keith B, Simon MC. Differential roles of hypoxia-inducible factor 1alpha (HIF-1alpha) and HIF-2alpha in hypoxic gene regulation. Mol Cell Biol. 2003;23(24):9361–74.
Massari F, Ciccarese C, Santoni M, Brunelli M, Piva F, Modena A, Bimbatti D, Fantinel E, Santini D, Cheng L, et al. Metabolic alterations in renal cell carcinoma. Cancer Treat Rev. 2015;41(9):767–76.
Rathmell WK, Rathmell JC, Linehan WM. Metabolic pathways in kidney cancer: current therapies and future directions. J Clin Oncol. 2018. https://doi.org/10.1200/JCO.2018.79.2309.
Shen C, Kaelin WG. The VHL/HIF axis in clear cell renal carcinoma. Semin Cancer Biol. 2013;23(1):18–25.
Shuch B, Linehan WM, Srinivasan R. Aerobic glycolysis: a novel target in kidney cancer. Expert Rev Anticancer Ther. 2013;13(6):711–9.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Dobbelstein M, Moll U. Targeting tumour-supportive cellular machineries in anticancer drug development. Nat Rev Drug Discov. 2014;13(3):179–96.
Xu F, Guan Y, Xue L, Huang S, Gao K, Yang Z, Chong T. The effect of a novel glycolysis-related gene signature on progression, prognosis and immune microenvironment of renal cell carcinoma. BMC Cancer. 2020;20(1):1207.
Simon AG, Esser LK, Ellinger J, Branchi V, Tolkach Y, Müller S, Ritter M, Kristiansen G, Muders MH, Mayr T, et al. Targeting glycolysis with 2-deoxy-D-glucose sensitizes primary cell cultures of renal cell carcinoma to tyrosine kinase inhibitors. J Cancer Res Clin Oncol. 2020;146(9):2255–65.
Li J, Zhang S, Liao D, Zhang Q, Chen C, Yang X, Jiang D, Pang J. Overexpression of PFKFB3 promotes cell glycolysis and proliferation in renal cell carcinoma. BMC Cancer. 2022;22(1):83.
Miao C, Liang C, Li P, Liu B, Qin C, Yuan H, Liu Y, Zhu J, Cui Y, Xu A, et al. TRIM37 orchestrates renal cell carcinoma progression via histone H2A ubiquitination-dependent manner. J Exp Clin Cancer Res. 2021;40(1):195.
Liang C, Wang S, Qin C, Bao M, Cheng G, Liu B, Shao P, Lv Q, Song N, Hua L, et al. TRIM36, a novel androgen-responsive gene, enhances anti-androgen efficacy against prostate cancer by inhibiting MAPK/ERK signaling pathways. Cell Death Dis. 2018;9(2):155.
Yeo W, Chan SL, Mo FKF, Chu CM, Hui JWY, Tong JHM, Chan AWH, Koh J, Hui EP, Loong H, et al. Phase I/II study of temsirolimus for patients with unresectable Hepatocellular Carcinoma (HCC)-a correlative study to explore potential biomarkers for response. BMC Cancer. 2015;15:395.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Putri GH, Anders S, Pyl PT, Pimanda JE, Zanini F. Analysing high-throughput sequencing data in Python with HTSeq 2.0. Bioinformatics. 2022;38(10):2943–5.
Gumz ML, Zou H, Kreinest PA, Childs AC, Belmonte LS, LeGrand SN, Wu KJ, Luxon BA, Sinha M, Parker AS, et al. Secreted frizzled-related protein 1 loss contributes to tumor phenotype of clear cell renal cell carcinoma. Clin Cancer Res. 2007;13(16):4740–9.
Tun HW, Marlow LA, von Roemeling CA, Cooper SJ, Kreinest P, Wu K, Luxon BA, Sinha M, Anastasiadis PZ, Copland JA. Pathway signature and cellular differentiation in clear cell renal cell carcinoma. PLoS ONE. 2010;5(5): e10696.
Wozniak MB, Le Calvez-Kelm F, Abedi-Ardekani B, Byrnes G, Durand G, Carreira C, Michelon J, Janout V, Holcatova I, Foretova L, et al. Integrative genome-wide gene expression profiling of clear cell renal cell carcinoma in Czech Republic and in the United States. PLoS ONE. 2013;8(3): e57886.
Eckel-Passow JE, Serie DJ, Bot BM, Joseph RW, Cheville JC, Parker AS. ANKS1B is a smoking-related molecular alteration in clear cell renal cell carcinoma. BMC Urol. 2014;14:14.
Eckel-Passow JE, Serie DJ, Bot BM, Joseph RW, Hart SN, Cheville JC, Parker AS. Somatic expression of ENRAGE is associated with obesity status among patients with clear cell renal cell carcinoma. Carcinogenesis. 2014;35(4):822–7.
von Roemeling CA, Radisky DC, Marlow LA, Cooper SJ, Grebe SK, Anastasiadis PZ, Tun HW, Copland JA. Neuronal pentraxin 2 supports clear cell renal cell carcinoma by activating the AMPA-selective glutamate receptor-4. Cancer Res. 2014;74(17):4796–810.
Liep J, Kilic E, Meyer HA, Busch J, Jung K, Rabien A. Cooperative effect of miR-141-3p and miR-145-5p in the regulation of targets in clear cell renal cell carcinoma. PLoS ONE. 2016;11(6): e0157801.
Wotschofsky Z, Gummlich L, Liep J, Stephan C, Kilic E, Jung K, Billaud J-N, Meyer H-A. Integrated microRNA and mRNA signature associated with the transition from the locally confined to the metastasized clear cell renal cell carcinoma exemplified by miR-146-5p. PLoS ONE. 2016;11(2): e0148746.
Nam H-Y, Chandrashekar DS, Kundu A, Shelar S, Kho E-Y, Sonpavde G, Naik G, Ghatalia P, Livi CB, Varambally S, et al. Integrative epigenetic and gene expression analysis of renal tumor progression to metastasis. Mol Cancer Res. 2019;17(1):84–96.
Su C, Lv Y, Lu W, Yu Z, Ye Y, Guo B, Liu D, Yan H, Li T, Zhang Q, et al. Single-cell RNA sequencing in multiple pathologic types of renal cell carcinoma revealed novel potential tumor-specific markers. Front Oncol. 2021;11: 719564.
Ghandi M, Huang FW, Jané-Valbuena J, Kryukov GV, Lo CC, McDonald ER, Barretina J, Gelfand ET, Bielski CM, Li H, et al. Next-generation characterization of the cancer cell line encyclopedia. Nature. 2019;569(7757):503–8.
Chen F, Chandrashekar DS, Varambally S, Creighton CJ. Pan-cancer molecular subtypes revealed by mass-spectrometry-based proteomic characterization of more than 500 human cancers. Nat Commun. 2019;10(1):5679.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, Varambally S. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia (New York, NY). 2017;19(8):649–58.
Qu Y, Feng J, Wu X, Bai L, Xu W, Zhu L, Liu Y, Xu F, Zhang X, Yang G, et al. A proteogenomic analysis of clear cell renal cell carcinoma in a Chinese population. Nat Commun. 2022;13(1):2052.
Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015;33(5):495–502.
Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, Hao Y, Stoeckius M, Smibert P, Satija R. Comprehensive integration of single-cell data. Cell. 2019;177(7):1888-1902.e21.
Liao J, Yu Z, Chen Y, Bao M, Zou C, Zhang H, Liu D, Li T, Zhang Q, Li J, et al. Single-cell RNA sequencing of human kidney. Sci Data. 2020;7(1):4.
Aran D, Looney AP, Liu L, Wu E, Fong V, Hsu A, Chak S, Naikawadi RP, Wolters PJ, Abate AR, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. 2019;20(2):163–72.
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, Wilson CJ, Lehár J, Kryukov GV, Sonkin D, et al. The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483(7391):603–7.
Lang B, Armaos A, Tartaglia GG. RNAct: Protein-RNA interaction predictions for model organisms with supporting experimental data. Nucleic Acids Res. 2019;47(D1):D601–6.
Bellucci M, Agostini F, Masin M, Tartaglia GG. Predicting protein associations with long noncoding RNAs. Nat Methods. 2011;8(6):444–5.
Xu Z, Liu M, Wang J, Liu K, Xu L, Fan D, Zhang H, Hu W, Wei D, Wang J. Single-cell RNA-sequencing analysis reveals MYH9 promotes renal cell carcinoma development and sunitinib resistance via AKT signaling pathway. Cell Death Discov. 2022;8(1):125.
Zhang C, He H, Hu X, Liu A, Huang D, Xu Y, Chen L, Xu D. Development and validation of a metastasis-associated prognostic signature based on single-cell RNA-seq in clear cell renal cell carcinoma. Aging (Albany NY). 2019;11(22):10183–202.
Maynard A, McCoach CE, Rotow JK, Harris L, Haderk F, Kerr DL, Yu EA, Schenk EL, Tan W, Zee A, et al. Therapy-induced evolution of human lung cancer revealed by single-cell RNA sequencing. Cell. 2020;182(5):1232-1251.e22.
Chung W, Eum HH, Lee H-O, Lee K-M, Lee H-B, Kim K-T, Ryu HS, Kim S, Lee JE, Park YH, et al. Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer. Nat Commun. 2017;8:15081.
Darmanis S, Sloan SA, Croote D, Mignardi M, Chernikova S, Samghababi P, Zhang Y, Neff N, Kowarsky M, Caneda C, et al. Single-cell RNA-Seq analysis of infiltrating neoplastic cells at the migrating front of human glioblastoma. Cell Rep. 2017;21(5):1399–410.
Tirosh I, Izar B, Prakadan SM, Wadsworth MH, Treacy D, Trombetta JJ, Rotem A, Rodman C, Lian C, Murphy G, et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science. 2016;352(6282):189–96.
Boroughs LK, DeBerardinis RJ. Metabolic pathways promoting cancer cell survival and growth. Nat Cell Biol. 2015;17(4):351–9.
Martínez-Reyes I, Chandel NS. Cancer metabolism: looking forward. Nat Rev Cancer. 2021;21(10):669–80.
Yu Y, Liang Y, Li D, Wang L, Liang Z, Chen Y, Ma G, Wu H, Jiao W, Niu H. Glucose metabolism involved in PD-L1-mediated immune escape in the malignant kidney tumour microenvironment. Cell Death Discov. 2021;7(1):15.
Feng C, Li Y, Li K, Lyu Y, Zhu W, Jiang H, Wen H. PFKFB4 is overexpressed in clear-cell renal cell carcinoma promoting pentose phosphate pathway that mediates Sunitinib resistance. J Exp Clin Cancer Res. 2021;40(1):308.
Nam H, Kundu A, Karki S, Brinkley GJ, Chandrashekar DS, Kirkman RL, Liu J, Liberti MV, Locasale JW, Mitchell T, et al. The TGF-β/HDAC7 axis suppresses TCA cycle metabolism in renal cancer. JCI Insight. 2021. https://doi.org/10.1172/jci.insight.148438.
Yang R-Z, Blaileanu G, Hansen BC, Shuldiner AR, Gong D-W. cDNA cloning, genomic structure, chromosomal mapping, and functional expression of a novel human alanine aminotransferase. Genomics. 2002;79(3):445–50.
Coloff JL, Murphy JP, Braun CR, Harris IS, Shelton LM, Kami K, Gygi SP, Selfors LM, Brugge JS. Differential glutamate metabolism in proliferating and quiescent mammary epithelial cells. Cell Metab. 2016;23(5):867–80.
Qian K, Zhong S, Xie K, Yu D, Yang R, Gong D-W. Hepatic ALT isoenzymes are elevated in gluconeogenic conditions including diabetes and suppressed by insulin at the protein level. Diabetes Metab Res Rev. 2015;31(6):562–71.
Cao Y, Lin S-H, Wang Y, Chin YE, Kang L, Mi J. Glutamic pyruvate transaminase GPT2 promotes tumorigenesis of breast cancer cells by activating sonic hedgehog signaling. Theranostics. 2017;7(12):3021–33.
Mitra D, Vega-Rubin-de-Celis S, Royla N, Bernhardt S, Wilhelm H, Tarade N, Poschet G, Buettner M, Binenbaum I, Borgoni S, et al. Abrogating GPT2 in triple-negative breast cancer inhibits tumor growth and promotes autophagy. Int J Cancer. 2021;148(8):1993–2009.
We would like to thank the researchers for their contributions. The flowchart was created with BioRender.com.
This research was supported by the National Natural Science Foundation of China (Grant No. 81771640), Project of Nanjing Science and Technology Committee (No. 201605001), “333” Project of Jiangsu Province (No. BRA2018083).
Ethics approval and consent to participate
The studies involving human participants and clinical samples were approved by the First Affiliated Hospital of Nanjing Medical University (No. 2021-SR-430). The patients provided their written informed consent to participate in this study. Animal studies were approved by the Institutional Animal Care and Use Committee (IACUC) of Nanjing Medical University (No. IACUC-1908003).
Consent for publication
The authors declare no competing interests that pertain to this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Single-cell RNA-seq revealed the distribution of Spectrin-family genes. (A). Quality control of single-cell RNA-seq samples (Nine ccRCC sample). The number of gene expressions in each cell, the sum of gene expressions, and the percentage of mitochondrial genes were illustrated. (B). The correlation of number of genes in the cells with the sum of gene expression and the percentage of mitochondrial genes. (C). 3000 hypervariable genes from all the genes shown in red and the top 10 hypervariable genes. (D). Distribution of different clusters in ccRCC tissues obtained by UMAP algorithm. (E). Heatmap showed the results of the cell cluster obtained by cell marker gene annotation were consistent with those obtained by “singleR” package annotation. Figure S2. Univariate and multivariate cox regression to assess the prognosis value of Spectrin-family genes. (A). Univariate and multivariate cox regression between Spectrin-family genes and overall survival (OS). (B). Univariate and multivariate cox regression between Spectrin-family genes and disease-specific survival (DSS). (C). Univariate and multivariate cox regression between Spectrin-family genes and progression survival interval (PFI). Figure S3. Correlation and function enrichment of Spectrin-family genes. (A). The expression correlation analysis of Spectrin-family genes. (B). Functional enrichment demonstrated that Spectrin-family genes were mainly involved in the interaction between L1 and Ankyrins. (*: P<0.05; **: P<0.01). Figure S4. Identification of SPTBN1 expression level in the TCGA and CCLE pan-cancer dataset. (A-B). Pan-cancer expression level of SPTBN1 form TCGA database (A) and CCLE database (B). (ns: no significant; *: P<0.05; **: P<0.01; ***:P<0.001). Figure S5. The expression level of SPTBN1 in GEO datasets and IHC staining. (A-F). The expression level of SPTBN1 in GSE40435 cohort (A), GSE53757 cohort (B), GSE6344 cohort (C), GSE46699 cohort (D), GSE105261 cohort (E) and GSE66270 cohort (F). (G). IHC staining revealed SPTBN1 was down-regulated in ccRCC tissue compared with adjacent normal renal tissue. Figure S6. The protein expression level of SPTBN1 among different clinical grades ccRCC patients. (A). The protein expression level of SPTBN1 among different grade ccRCC patients in CPTAC cohort. (B). The protein expression level of SPTBN1 among different grade ccRCC patients in FUSCC-ccRCC cohort among Chinese ccRCC patients. (**: P<0.05). Figure S7. The correlation between SPTBN1 and survival outcomes in different clinicopathological characteristics subgroups. (A). The correlation between SPTBN1 and survival outcomes among age<=60 and age>60 subgroups. (B). The correlation between SPTBN1 and survival outcomes among female and male subgroups. (C). The correlation between SPTBN1 and survival outcomes among stage I-II and stage III-IV subgroups. (D). The correlation between SPTBN1 and survival outcomes among T1-T2 and T3-T4 subgroups. Figure S8. Validation of SPTBN1 expression level after SPTBN1-knockdown and SPTBN1-overexpressing by qRT-PCR and WB. (A). qRT-PCR and WB assays validated the expression of SPTBN1 after knockdown SPTBN1. (B). qRT-PCR and WB assays validated the expression of SPTBN1 after overexpression SPTBN1. (C). qRT-PCR and WB assays validated the expression of GPT2 after knockdown GPT2. (**: P<0.05). Figure S9. Transwell cell migration and invasion assays after SPTBN1-knockdown and SPTBN1-overexpressing of RCC cells. (**: P<0.05). Figure S10. IHC staining revealed GPT2 was up-regulated in ccRCC tissue compared with adjacent normal renal tissue. Table S1. Oligonucleotide sequences used in this research
About this article
Cite this article
Wu, J., Miao, C., Wang, Y. et al. SPTBN1 abrogates renal clear cell carcinoma progression via glycolysis reprogramming in a GPT2-dependent manner. J Transl Med 20, 603 (2022). https://doi.org/10.1186/s12967-022-03805-w