Skip to main content
Download PDF

Bioinformatic analysis highlights SNHG6 as a putative prognostic biomarker for kidney renal papillary cell carcinoma

BMC Urology volume 23, Article number: 54 (2023) Cite this article

Abstract

Purpose

Kidney renal papillary cell carcinoma (KIRP) is a highly heterogeneous malignancy and current systemic therapeutic strategies are difficult to achieve a satisfactory outcome for advanced disease. Meanwhile, there is a lack of effective biomarkers to predict the prognosis of KIRP.

Methods

Using TCGA, GTEx, UALCAN, TIMER, TIMER 2.0 and STRING databases, we analyzed the relationship of SNHG6 with KIRP subtypes, tumor-infiltrating immune cells and potential target mRNAs. Based on TCGA data, ROC curves, Kaplan–Meier survival analysis and COX regression analysis were performed to evaluate the diagnostic and prognostic value of SNHG6 in KIRP. Nomogram was used to predict 3- and 5-year disease-specific survival in KIRP patients. In addition, with the help of Genetic ontology and Gene set enrichment analysis, the biological processes and signalling pathways that SNHG6 may be involved in KIRP were initially explored.

Results

In patients with KIRP, SNHG6 was significantly upregulated and associated with a more aggressive subtype (lymph node involvement, pathological stage IV, CIMP phenotype) and poor prognosis. The ROC curve showed good diagnostic efficacy (AUC value: 0.828) and the C-index of the Nomogram for predicting DSS at 3 and 5 years was 0.920 (0.898–0.941). In the immune microenvironment of KIRP, SNHG6 expression levels were negatively correlated with macrophage abundance and positively correlated with cancer-associated fibroblasts. Furthermore, SNHG6 may promote KIRP progression by regulating the expression of molecules such as AURKB, NDC80, UBE2C, NUF2, PTTG1, CENPH, SPC25, CDCA3, CENPM, BIRC5, TROAP, EZH2. Last, GSEA suggests that SNHG6 may be involved in the regulation of the PPAR signalling pathway and the SLIT/ROBO signalling pathway.

Conclusions

Our analysis suggests that a high SNHG6 expression status in KIRP is associated with a poorer prognosis for patients, and also elucidates some potential mechanisms contributing to this poorer outcome. This may provide new insights into the treatment and management of KIRP in the foreseeable future.

Peer Review reports

Introduction

Kidney renal papillary cell carcinoma (KIRP), is the most common type of non-clear cell renal cell carcinoma (non-ccRCC) in the worldwide, accounting for approximately 15–20% of all renal cell carcinomas [1]. The KIRP is histologically divided into two major types, with type I KIRP being more associated with MET mutations, while type II is considered to be more heterogeneous and subdivided into various subtypes based on the genetic and molecular composition of the tumor [2, 3]. At present, systemic treatment strategies for advanced renal cell carcinoma revolve around ccRCC, either with immune checkpoint inhibitors or targeted therapies. However, the relative lack of specific gene mutations in ccRCC, such as VHL and PBRM1, decreases response rates to targeted therapeutic strategies in KIRP patients [4, 5]. Meanwhile, immune checkpoint inhibitors have shown limited activity in the KIRP population [6]. Moreover, the heterogeneity of the tumor itself presents another major challenge for therapeutic strategies. In this context, it is essential to refine the molecular and immunological landscape of KIRP.

Long non-coding RNAs (LncRNAs) are an RNA transcription with over 200 nucleotides and no apparent protein coding capability [7]. Its overexpression, defects or mutations have been found to be associated with a variety of human diseases, including cancer, neurological disorders, cardiovascular diseases and more [8]. Small nucleolar RNA host gene 6 (SNHG6), a lncRNA located in chromosome 8q13.1, has been shown to be linked to poor outcomes in a variety of human cancers [9]. For instance, SNHG6 can promote the progression of colorectal and ovarian clear cell carcinomas by regulating the expression of enhancer of zeste homolog 2 (EZH2) through a ceRNA network [10, 11]. In xenograft mice with lung cancer, SNHG6 can also regulate the differentiation of myeloid-derived suppressor cells by inhibiting protein expression of EZH2 via ubiquitin [12]. Furthermore, SNHG6 can interact with heterogeneous nuclear ribonucleoprotein L (HNRNPL) and polypyrimidine tract binding protein 1 (PTBP1) to promote the progression of hepatocellular carcinoma [13]. Taken together, SNHG6 may be implicated in the regulation of tumor development by multiple pathways and is considered a poor prognostic biomarker for human cancers. Whereas, it is still unclear whether SNHG6 has a place in the tumor microenvironment of KIRP as well.

Here, based on the TCGA database, we were the first to analyze the expression levels of SNHG6 in KIRP, the association of SNHG6 with KIRP subtypes and its predictive value for KIRP. Subsequently, in the context of these assumptions holding true, we further explored the relationship between SNHG6 and the abundance of tumor-infiltrating immune cells in the tumor microenvironment, potential mRNAs downstream of SNHG6, and the biological processes and signalling pathways that may be involved in SNHG6 in KIRP.

Methods

Data access

Sequencing data and corresponding clinicopathological information for 288 KIRP samples and 32 adjacent normal samples were downloaded from the TCGA database (https://portal.gdc.cancer.gov/) (14). 136 SNHG6 expression data in normal kidney tissue samples were obtained from the GTEx database (https://www.gtexportal.org/) to compensate for the relative lack of normal samples when performing differential gene expression analysis [15].

UALCAN database

UALCAN (http://ualcan.path.uab.edu/) is a comprehensive, publicly available web resource that encompasses multiple cancer OMICS data (including TCGA, MET500, CPTAC and CBTTC) [16]. We used it to comprehensively analyze the relationship between SNHG6 expression and age, gender, race, pathological stage and histological subtype of KIRP patients.

Nomogram construction and evaluation

First, 282 patients with KIRP were included in univariate and multivariate COX regression analyses to evaluate the impact of clinicopathological information (age, gender, BMI, race, smoking, pathological stage, and SNHG6 expression levels) on patients' disease-specific survival (DSS). Subsequently, indicators that were statistically significant in the multivariate COX analysis were included in the construction of the Nomogram to predict patients' DSS at 3 and 5 years. The calibration curve drawn by the "rms" R package is used to evaluate the agreement between the predicted and actual values.

TIMER databases

TIMER is a comprehensive resource containing information on the immune infiltration of multiple cancer types and is currently available in TIMER (https://cistrome.shinyapps.io/timer/) and TIMER 2.0 versions (http://timer.cistrome.org/)[17, 18]. We used these tools to explore the correlation of SNHG6 with tumor-infiltrating immune cells in KIRP; for statistically significant cell types, Kaplan–Meier curves were further plotted to analyze the impact of these correlations on overall survival.

Target mRNA screening and protein–protein interaction (PPI)

On the one hand, we performed gene correlation analysis for SNHG6, and mRNAs with spearman's correlation coefficient > 0.3 were selected. On the other hand, we screened the differentially expressed genes (DEGs) in KIRP based on logFC > 1.5 and plotted the volcano plot. subsequently, the two sets of screened genes were analyzed by intersection analysis via a Venn diagram and the target molecules were used in the next step of PPI and genetic ontology (GO) analysis. PPI networks for the target mRNAs were constructed in the STRING database (https://cn.string-db.org/) and next imported into Cytoscape (version 3.9.1) to seek out more critical molecules.

GO and gene set enrichment analysis (GSEA)

GO analysis was used to evaluate the biological processes, molecular functions, and cellular composition of the target mRNAs. GSEA explored the potential signalling pathways of SNHG6 in KIRP by gathering differential genes between the high and low expression of SNHG6 groups. 1000 operations were performed per analysis, with normalized enrichment score > 1, false discovery rate < 0.25 and a nominal P value < 0.05 being considered statistically significant. These were performed using the 'clusterprofiler' of the R package.

Statistical analyses

RNAseq data obtained from TCGA and GTEx databases were analyzed using the R package (version 3.6.3). Wilcoxon signed rank test, Wilcoxon rank sum test and Kruskal–Wallis test were used for intergroup comparison and analysis of the relationship between SNHG6 and clinicopathological features. ROC curves and Kaplan–Meier analysis were used to assess the diagnostic and prognostic predictive value of SNHG6 for KIRP. Spearman correlation analysis was used to evaluate the correlation between SNHG6 expression levels with tumor-infiltrating immune cells and mRNA. p < 0.05 was considered statistically significant.

Results

Pan-cancer analysis of SNHG6

By integrating genomic data from the TCGA database and GTEx database, we found that SNHG6 was highly expressed in a variety of human cancers, including glioma, renal cancer, gastric cancer, colorectal cancer, head and neck squamous cell carcinoma, lung squamous carcinoma and melanoma, etc. (see Fig. 1A for details). Subsequently, we comprehensively analyzed the effects of SNHG6 on overall survival (OS), DSS, disease-free interval (DFI), and progression-free interval (PFI) of cancer patients, and the results showed that upregulation of SNHG6 may be detrimental to the prognosis of KIRP patients. (Fig. 1B–E). Therefore, KIRP was chosen as the target cancer for this study.

Fig. 1
figure 1

Expression of SNHG6 in pan-cancer and the effect of SNHG6 upregulation on the prognosis of human cancers. A SNHG6 expression is upregulated in a variety of cancer tissues. Effect of upregulation of SNHG6 expression on B OS, C DSS, D DFI and E PFI in human cancer patients. *, p < 0.05; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001. SNHG6, Small nucleolar RNA host gene 6; KIRP, kidney renal papillary cell carcinoma; OS, overall survival; DSS, Disease-specific survival; DFI, disease-free interval; PFI, progression-free interval

Full size image

Expression levels of SNHG6 in KIRP and correlation with KIRP subtypes

As shown in Fig. 2A, SNHG6 was upregulated in KIRP tissues with or without the addition of data from the GTEx database. In the meantime, this phenomenon could be similarly observed in the 32 paired samples. Supported by this evidence, we further compared the expression levels of SNHG6 in different subgroups of KIRP. The results showed that the expression of SNHG6 was not significantly correlated with the age, gender and racial differences of the patients, while it was significantly upregulated in patients with lymph node metastasis and pathological stage IV (p < 0.01) (Fig. 2B). The CpG island methylator phenotype (CIMP) has been found to be strongly associated with poor prognosis in KIRP patients, and therefore SNHG6 expression levels in the CIMP group have been explored [19]. As expected, the expression of SNHG6 was significantly higher in the CIMP group than in classical type I and type II KIRP patients.

Fig. 2
figure 2

Expression of SNHG6 in KIRP and the relationship between SNHG6 expression and KIRP subtypes. A SNHG6 is highly expressed in KIRP tissues (based on TCGA and GTEx databases). B High SNHG6 expression was observed in patients with lymph node involvement, pathological stage IV, CIMP phenotype; and there was no significant correlation with gender, age and race. **, p < 0.01; ***, p < 0.001; ns, not significant. SNHG6, Small nucleolar RNA host gene 6; KIRP, kidney renal papillary cell carcinoma

Full size image

The diagnostic and prognostic utility of SNHG6 in KIRP

ROC curves were used to assess the diagnostic efficacy of SNHG6 in KIRP. As shown in Fig. 3A, the AUC value for SNHG6 was 0.828 when differentiating KIRP tissue from normal tissue; for stage I-II and III-IV patients, the AUC values were 0.819 and 0.850, respectively. Kaplan–Meier curves visualised the impact of SNHG6 expression on oncological outcomes in KIRP patients; the high SNHG6 expression group was associated with worse PFI, DSS and OS in KIRP patients with hazard ratios (HR) of 2.28, 3.59 and 2.02 respectively (Fig. 3B) (p < 0.05). Subsequently, we implemented univariate and multivariate COX analyses to look for potential predictors of DSS, and pathological stage. SNHG6 expression were found to have non-negligible predictive value (Table 1). Accordingly, these two metrics were incorporated into the construction of the Nomogram to predict disease-specific survival at 3 and 5 years in KIRP patients (Fig. 3C). Encouragingly, the C-index was 0.920 (0.898–0.941) and the good agreement between predicted and actual values was also objectively demonstrated by the calibration curve in Fig. 3D.

Fig. 3
figure 3

Diagnostic and prognostic value analysis of SNHG6. A ROC curves of SNHG6 expression in normal and tumor tissues. B Kaplan–Meier analysis showed that KIRP patients with high SNHG6 expression had worse OS, DSS and PFI. C Nomogram for predicting DSS at 3 and 5 years in KIRP patients. D Calibration curve to evaluate Nomogram consistency. SNHG6, Small nucleolar RNA host gene 6; KIRP, kidney renal papillary cell carcinoma; OS, overall survival; DSS, disease-specific survival; PFI, progression-free interval

Full size image
Table 1 Univariate and Multivariate COX regression analysis for DSS
Full size table

The relationship between SNHG6 and tumor-infiltrating immune cells

Using the TIMER database, we found that there was no significant correlation between SNHG6 and B lymphocytes, T lymphocytes and dendritic cells in the tumor microenvironment of KIRP, but was positively correlated with neutrophils and cancer associated fibroblasts, and negatively correlated with macrophages (Fig. 4A). To more visually demonstrate the significance of this correlation on the survival prognosis in KIRP patients, we combined those two variables to plot the Kaplan–Meier curves (Fig. 4B). Interestingly, lower macrophage abundance in the tumor microenvironment under high expression state of SNHG6 led to a deterioration of OS in KIRP patients. At the same time, we found that patients in the SNHG6 high expression plus high abundance of cancer-associated fibroblasts group had a worse OS than any other subgroup. As a result, the pro-oncogenic effect of SNHG6 on KIRP may be partially related to tumor-infiltrating immune cells, especially macrophages and cancer associated fibroblasts.

Fig. 4
figure 4

Analysis of the correlation between SNHG6 and tumor-infiltrating immune cells and the impact of this correlation on OS in KIRP patients. A There was no significant correlation between SNHG6 and B lymphocytes, T lymphocytes and dendritic cells, but a positive correlation with neutrophils and cancer associated fibroblasts, and a negative correlation with macrophages. B High cancer-associated fibroblast abundance and low macrophage abundance worsened OS in KIRP patients. SNHG6, Small nucleolar RNA host gene 6; KIRP, kidney renal papillary cell carcinoma; OS, overall survival

Full size image

mRNAs potentially regulated by SNHG6

Following the screening criteria, we found that SNHG6 was closely associated with 1506 mRNAs in KIRP, and the top 20 molecules have been shown via the heat map in Fig. 5A. On the other hand, differentially expressed mRNAs in KIRP were screened by the volcano map, where 1500 mRNAs were up-regulated (Fig. 5B). Venn diagram further identified overlapping genes in both groups (Fig. 5C). These 88 target molecules were then constructed into a PPI network through the STRING database and hub genes were screened based on the centrality of the nodes (Fig. 5D). Furthermore, Kaplan–Meier survival analysis was performed for mRNAs (AURKB, NDC80, UBE2C, NUF2, PTTG1, CENPH, SPC25, CDCA3, CENPM, BIRC5, TROAP, EZH2) with more edges to explore their impact on OS in KIRP patients. Intriguingly, the upregulation of almost all these molecules (except NDC80) in KIRP worsened the OS of the patients (Fig. 6A–L). Apparently, the adverse effect of these mRNAs on the survival prognosis of KIRP is in line with the trend of SNHG6.

Fig. 5
figure 5

Screening of target mRNAs associated with SNHG6 A The top 20 mRNAs positively correlated with SNHG6. B A volcano map showing DEGs in the KIRP. C An intersection analysis of SNHG6-related genes and DEGs. D PPI network built with STRING database and Cytoscape. SNHG6, Small nucleolar RNA host gene 6; KIRP, kidney renal papillary cell carcinoma; PPI, Protein–protein interaction; DEGs: differentially expressed genes

Full size image
Fig. 6
figure 6

Kaplan–Meier analysis exploring the impact of hub gene expression on OS in KIRP patients. A AURKB. B NDC80. C UBE2C. D NUF2. E PTTG1. F CENPH. G SPC25. H CDCA3. I CENPM. J BIRC5. K TROAP. L EZH2. SNHG6, Small nucleolar RNA host gene 6; KIRP, kidney renal papillary cell carcinoma

Full size image

Exploration of the potential mechanisms of SNHG6 in KIRP

Based on the 88 target mRNAs obtained from the previous screening, we conducted GO analysis and found that SNHG6 may be involved in the following biological processes, including: T cell activation, regulation of peptidase activity, regulation of DNA metabolic process, negative regulation of cell cycle process, regulation of lymphocyte activation, etc. Meanwhile, it is involved in cellular components such as the chromosome centromeric region, kinetochore and proteasome core complex, and its potential molecular function is mainly to influence enzyme activity (Fig. 7A). In addition, GSEA was used to explore the potential signalling pathway of SNHG6 in KIRP. As shown in Figs. 7B, the KEGG pathways that were significantly enriched included the ribosome, proximal tubule bicarbonate reclamation, pyruvate metabolism, citrate cycle tca cycle, and peroxisome proliferator-activated receptor (PPAR) signaling pathway [20,21,22]. The REACTOME pathway includes translation, Roundabout (ROBO) signalling receptor transmission, RRNA processing, regulation of short-lived non-coding transcripts (SLITs) and ROBOs expression and initiation of eukaryotic translation.

Fig. 7
figure 7

GO analysis and GSEA. A GO analysis: biological processes, molecular functions and cellular composition. B Enrichment plots by GSEA. GO, Gene Ontology. GSEA, gene set enrichment analysis

Full size image

Discussion

KIRP is the second most common subtype of renal cell carcinomas and is characterized by the presence of basophilic or eosinophilic cells in papillary or tubular structures [23]. KIRP has long been considered a low-grade malignancy; and in patients with localized renal cell carcinomas, the papillary subtype is an independent prognostic factor for improved oncological outcome [24]. Conversely, a meta-analysis showed no significant improvement in the prognosis of patients with metastatic KIRP compared to patients with metastatic ccRCC, with patients with type II KIRP having an even worse outcome [25]. Furthermore, KIRP is a highly heterogeneous tumor and its prognosis is significantly associated with histological subtypes and genetic phenotypes. In this context, the treatment and management of advanced KIRP undoubtedly faces enormous challenges. Accordingly, more promising biomarkers or therapeutic targets need to be explored in the current work.

Human SNHG6 is the housekeeping gene of the 5′ TOP family, consisting of five transcripts, SNHG6-201 to SNHG6-205 [26]. High expression of SNHG6 has been found in a variety of cancer tissues (including colorectal cancer, hepatocellular carcinoma, breast cancer, ccRCC, etc.) and can predict a poor prognosis [11, 27,28,29]. In the present study, these results were also largely validated; in addition to this, we found that SNHG6 was also significantly upregulated in KIRP and strongly associated with unfavorable oncological outcomes. Patients with KIRP have a higher incidence of lymph node involvement compared to the metastatic pattern of ccRCC, which makes lymph node involvement more critical in the pathological grading system of KIRP [30]. Interestingly, our results showed that SNHG6 expression was positively correlated with lymph node metastatic lesions. Furthermore, we found that SNHG6 was significantly upregulated in KIRP patients with stage IV pathological stage and those carrying the CIMP phenotype. Not surprisingly, these results suggest that SNHG6 expression is associated with the more aggressive KIRP subtype.

To further explore the potential of SNHG6 as a biomarker for KIRP, we have comprehensively assessed its diagnostic and prognostic predictive value with the help of ROC curves and Kaplan–Meier analysis. As expected, SNHG6 showed good diagnostic performance and its up-regulated status predicted worse PFI, DSS and OS; this is in line with the findings of previous studies [9]. Subsequently, based on the results of the COX analysis, we constructed a Nomogram for predicting 3-year DSS and 5-year DSS in KIRP patients and showed satisfactory predictive efficacy. Consequently, based on the above data, there is reason to believe that SNHG6 is an independent predictor of poor prognosis in KIRP patients and may contribute in some way to the progression of KIRP.

Renal cell carcinoma has long been recognized as an immunogenic tumor and a variety of immune cells in the tumor microenvironment can be involved in mediating the anti-tumor immune response [31]. The presence of tumor-associated macrophages has been independently shown to reduce the risk of death in patients with KIRP [32]. Cancer associated fibroblasts have also been shown to play multiple roles in tumor development by secreting growth factors, inflammatory ligands and extracellular matrix proteins that promote cancer cell proliferation and immune rejection [33].According to the TIMER database, we found that SNHG6 expression in KIRP was negatively correlated with macrophage abundance and positively correlated with cancer associated fibroblasts. Furthermore, survival analyses further confirmed that this correlation have deleterious effects on OS in KIRP patients (Fig. 4B). Therefore, we speculated that SNHG6 may influence disease progression to some extent by participating in the regulation of the immune microenvironment of KIRP.

There is growing evidence that LncRNAs can be involved in pre-mRNA alternative splicing as transcriptional regulators in the nucleus on the one hand, and in post-transcriptional regulation (mRNA stability, mRNA translation, protein stability and ceRNA networks) in the cytoplasm on the other hand [34]. After qualifying the expression level of mRNAs with strong SNHG6 correlation in KIRP, we obtained 88 target molecules. Subsequently, we performed GO analysis using these 88 genes and the results showed that SNHG6 may be involved in a range of biological processes including: activation and proliferation of immune cells, regulation of peptidase activity and regulation of DNA metabolic process, etc. Encouragingly, these results seem to corroborate the previously mentioned role of SNHG6 in the regulation of the tumor immune microenvironment. Additionally, we performed PPI network construction using 88 target molecules and identified several hub genes such as AURKB, NDC80, UBE2C, NUF2, PTTG1, CENPH, SPC25, CDCA3, CENPM, BIRC5, TROAP, EZH2, etc. Also, upregulation of the expression of these molecules has been found to predict a poor prognosis for KIRP patients. Intriguingly, previous studies have shown that these molecules can contribute to the progression of a variety of cancers [35,36,37,38,39,40,41,42], with PTTG1 and CDCA3 being found to be prognostic biomarkers for KIRP [43, 44]. On the other hand, AURKB, BIRC5 and SPC25 were found to be significantly upregulated in CIMP-positive ccRCC tissues and were associated with poor prognosis [45]. EZH2 is considered to be a new target for cancer therapy and has been shown to be a downstream target gene regulated by SNHG6 in a variety of malignancies [11, 42, 46, 47]. Thus, given this favorable evidence, it seems reasonable to presume that SNHG6 may influence the biological behavior of KIRP by regulating the expression of these target genes.

Finally, we explored the potential mechanisms of SNHG6 in KIRP, and GSEA showed that SNHG6 is closely associated with the PPAR signalling pathway, Robo signalling receptor transmission and the regulation of Slit and Robo expression. To our knowledge, the PPAR signalling pathway can function pleiotropically in cancer, and PPARα antagonists can be involved in the regulation of multiple reprogrammed metabolic pathways and attenuate tumor growth in renal cell carcinomas [48, 49]. SLITs are a series of secreted proteins that regulate angiogenesis, inflammatory cell chemotaxis, tumor cell migration and metastasis by binding to ROBO receptors [50, 51]. However, GSEA can only provide researchers with preliminary evidence and specific signalling pathways will need to be further explored in subsequent studies.

It has to be mentioned that although this study systematically analyses the potential association of SNHG6 with the malignant phenotype of KIRP, it still has some shortcomings. For example, the data in this study were derived from public databases and lacked validation from the clinical sample in our study centre. Secondly, the exact mechanism of SNHG6 in the tumourigenesis development of KIRP remains unclear and needs to be refined in subsequent in vivo and in vitro experiments.

Conclusion

Here, SNHG6 expression was found to be significantly upregulated in KIRP tissues and associated with a more aggressive KIRP subtype and poorer prognosis. Besides, SNHG6 may interfere with anti-tumor immune responses by affecting macrophages and cancer-associated fibroblasts in the tumor microenvironment, and may also promote KIRP progression by regulating the expression of molecules such as EZH2. Meanwhile, the PPAR signalling pathway and SLIT/ROBO signalling pathway may be specific signalling pathways for SNHG6 in KIRP.

Availability of data and materials

RNA-seq data for KIRP patients are available from both the TCGA database (https://portal.gdc.cancer.gov/) and the GTEx database (https://www.gtexportal.org/).

References

  1. de Vries-Brilland M, McDermott DF, Suárez C, Powles T, Gross-Goupil M, Ravaud A, et al. Checkpoint inhibitors in metastatic papillary renal cell carcinoma. Cancer Treat Rev. 2021;99: 102228.

    Article  PubMed  Google Scholar 

  2. Akhtar M, Al-Bozom IA, Al HT. Papillary renal cell carcinoma (PRCC): an update. Adv Anat Pathol. 2019;26(2):124–32.

    Article  PubMed  Google Scholar 

  3. Rhoades Smith KE, Bilen MA. A review of papillary renal cell carcinoma and MET inhibitors. Kidney Cancer. 2019;3(3):151–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Cancer Genome Atlas Research Network. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature. 2013;499(7456):43–9.

    Article  Google Scholar 

  5. Connor Wells J, Donskov F, Fraccon AP, Pasini F, Bjarnason GA, Beuselinck B, et al. Characterizing the outcomes of metastatic papillary renal cell carcinoma. Cancer Med. 2017;6(5):902–9.

    Article  PubMed  PubMed Central  Google Scholar 

  6. de Vries-Brilland M, Gross-Goupil M, Seegers V, Boughalem E, Beuselinck B, Thibault C, et al. Are immune checkpoint inhibitors a valid option for papillary renal cell carcinoma? A multicentre retrospective study. Eur J Cancer. 2020;136:76–83.

    Article  PubMed  Google Scholar 

  7. Quinn JJ, Chang HY. Unique features of long non-coding RNA biogenesis and function. Nat Rev Genet. 2016;17(1):47–62.

    Article  CAS  PubMed  Google Scholar 

  8. Esteller M. Non-coding RNAs in human disease. Nat Rev Genet. 2011;12(12):861–74.

    Article  CAS  PubMed  Google Scholar 

  9. Shen H, Mo Q, Xu X, Liu B. The prognostic value of lncRNA SNHG6 in cancer patients. Cancer Cell Int. 2020;20:286.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Wu Y, Deng Y, Guo Q, Zhu J, Cao L, Guo X, et al. Long non-coding RNA SNHG6 promotes cell proliferation and migration through sponging miR-4465 in ovarian clear cell carcinoma. J Cell Mol Med. 2019;23(8):5025–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Xu M, Chen X, Lin K, Zeng K, Liu X, Xu X, et al. lncRNA SNHG6 regulates EZH2 expression by sponging miR-26a/b and miR-214 in colorectal cancer. J Hematol Oncol. 2019;12(1):3.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Lu W, Cao F, Feng L, Song G, Chang Y, Chu Y, et al. LncRNA Snhg6 regulates the differentiation of MDSCs by regulating the ubiquitination of EZH2. J Hematol Oncol. 2021;14(1):196.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wang H, Ma P, Liu P, Guo D, Liu Z, Zhang Z. lncRNA SNHG6 promotes hepatocellular carcinoma progression by interacting with HNRNPL/PTBP1 to facilitate SETD7/LZTFL1 mRNA destabilization. Cancer Lett. 2021;520:121–31.

    Article  CAS  PubMed  Google Scholar 

  14. Tomczak K, Czerwińska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge. Contemp Oncol (Pozn). 2015;19(1a):A68-77.

    PubMed  Google Scholar 

  15. GTEx Consortium. Human genomics. The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science. 2015;348(6235):648–60.

    Article  Google Scholar 

  16. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Can Res. 2017;77(21):e108–10.

    Article  CAS  Google Scholar 

  18. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Linehan WM, Ricketts CJ. The cancer genome atlas of renal cell carcinoma: findings and clinical implications. Nat Rev Urol. 2019;16(9):539–52.

    Article  CAS  PubMed  Google Scholar 

  20. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci Publ Protein Soc. 2019;28(11):1947–51.

    Article  CAS  Google Scholar 

  22. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.

    Article  CAS  PubMed  Google Scholar 

  23. Mendhiratta N, Muraki P, Sisk AE Jr, Shuch B. Papillary renal cell carcinoma: review. Urol Oncol. 2021;39(6):327–37.

    Article  PubMed  Google Scholar 

  24. Nguyen DP, Vertosick EA, Corradi RB, Vilaseca A, Benfante NE, Touijer KA, et al. Histological subtype of renal cell carcinoma significantly affects survival in the era of partial nephrectomy. Urol Oncol. 2016;34(6):259.e1-8.

    Article  PubMed  Google Scholar 

  25. Deng J, Li L, Xia H, Guo J, Wu X, Yang X, et al. A comparison of the prognosis of papillary and clear cell renal cell carcinoma: evidence from a meta-analysis. Medicine (Baltimore). 2019;98(27): e16309.

    Article  CAS  PubMed  Google Scholar 

  26. Wang HS, Zhang W, Zhu HL, Li QP, Miao L. Long noncoding RNA SNHG6 mainly functions as a competing endogenous RNA in human tumors. Cancer Cell Int. 2020;20:219.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Chen K, Hou Y, Liao R, Li Y, Yang H, Gong J. LncRNA SNHG6 promotes G1/S-phase transition in hepatocellular carcinoma by impairing miR-204-5p-mediated inhibition of E2F1. Oncogene. 2021;40(18):3217–30.

    Article  CAS  PubMed  Google Scholar 

  28. Lv Y, Lv X, Yang H, Qi X, Wang X, Li C, et al. LncRNA SNHG6/miR-125b-5p/BMPR1B axis: a new therapeutic target for triple-negative breast cancer. Front Oncol. 2021;11: 678474.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zhao P, Deng Y, Wu Y, Guo Q, Zhou L, Yang X, et al. Long noncoding RNA SNHG6 promotes carcinogenesis by enhancing YBX1-mediated translation of HIF1α in clear cell renal cell carcinoma. Faseb j. 2021;35(2): e21160.

    Article  CAS  PubMed  Google Scholar 

  30. Margulis V, Tamboli P, Matin SF, Swanson DA, Wood CG. Analysis of clinicopathologic predictors of oncologic outcome provides insight into the natural history of surgically managed papillary renal cell carcinoma. Cancer. 2008;112(7):1480–8.

    Article  PubMed  Google Scholar 

  31. Díaz-Montero CM, Rini BI, Finke JH. The immunology of renal cell carcinoma. Nat Rev Nephrol. 2020;16(12):721–35.

    Article  PubMed  Google Scholar 

  32. Hutterer GC, Pichler M, Chromecki TF, Strini KA, Klatte T, Pummer K, et al. Tumour-associated macrophages might represent a favourable prognostic indicator in patients with papillary renal cell carcinoma. Histopathology. 2013;63(3):309–15.

    Article  PubMed  Google Scholar 

  33. Biffi G, Tuveson DA. Diversity and biology of cancer-associated fibroblasts. Physiol Rev. 2021;101(1):147–76.

    Article  CAS  PubMed  Google Scholar 

  34. Gao N, Li Y, Li J, Gao Z, Yang Z, Li Y, et al. Long non-coding RNAs: the regulatory mechanisms, research strategies, and future directions in cancers. Front Oncol. 2020;10: 598817.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Zhang X, Wang F, Wang Z, Yang X, Yu H, Si S, et al. ALKBH5 promotes the proliferation of renal cell carcinoma by regulating AURKB expression in an m(6)A-dependent manner. Ann Transl Med. 2020;8(10):646.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Qu Y, Li J, Cai Q, Liu B. Hec1/Ndc80 is overexpressed in human gastric cancer and regulates cell growth. J Gastroenterol. 2014;49(3):408–18.

    Article  CAS  PubMed  Google Scholar 

  37. Chen Z, Wang L. The clinical significance of UBE2C gene in progression of renal cell carcinoma. Eur J Histochem. 2021. https://doi.org/10.4081/ejh.2021.3196.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Xie X, Lin J, Fan X, Zhong Y, Chen Y, Liu K, et al. LncRNA CDKN2B-AS1 stabilized by IGF2BP3 drives the malignancy of renal clear cell carcinoma through epigenetically activating NUF2 transcription. Cell Death Dis. 2021;12(2):201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Zhang H, Li W, Gu W, Yan Y, Yao X, Zheng J. MALAT1 accelerates the development and progression of renal cell carcinoma by decreasing the expression of miR-203 and promoting the expression of BIRC5. Cell Prolif. 2019;52(5): e12640.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Xiao Y, Najeeb RM, Ma D, Yang K, Zhong Q, Liu Q. Upregulation of CENPM promotes hepatocarcinogenesis through mutiple mechanisms. J Exp Clin Cancer Res. 2019;38(1):458.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Gao B, Wang L, Zhang Y, Zhang N, Han M, Liu H, et al. MiR-532-3p suppresses cell viability, migration and invasion of clear cell renal cell carcinoma through targeting TROAP. Cell Cycle. 2021;20(16):1578–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Duan R, Du W, Guo W. EZH2: a novel target for cancer treatment. J Hematol Oncol. 2020;13(1):104.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Tian X, Xu WH, Xu FJ, Li H, Anwaier A, Wang HK, et al. Identification of prognostic biomarkers in papillary renal cell carcinoma and PTTG1 may serve as a biomarker for predicting immunotherapy response. Ann Med. 2022;54(1):211–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Li H, Li M, Yang C, Guo F, Deng S, Li L, et al. Prognostic value of CDCA3 in kidney renal papillary cell carcinoma. Aging. 2021;13(23):25466–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Arai E, Gotoh M, Tian Y, Sakamoto H, Ono M, Matsuda A, et al. Alterations of the spindle checkpoint pathway in clinicopathologically aggressive CpG island methylator phenotype clear cell renal cell carcinomas. Int J Cancer. 2015;137(11):2589–606.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Li Y, Li D, Zhao M, Huang S, Zhang Q, Lin H, et al. Long noncoding RNA SNHG6 regulates p21 expression via activation of the JNK pathway and regulation of EZH2 in gastric cancer cells. Life Sci. 2018;208:295–304.

    Article  CAS  PubMed  Google Scholar 

  47. Pu FF, Shi DY, Chen T, Liu YX, Zhong BL, Zhang ZC, et al. SP1-induced long non-coding RNA SNHG6 facilitates the carcinogenesis of chondrosarcoma through inhibiting KLF6 by recruiting EZH2. Cell Death Dis. 2021;12(1):59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Abu Aboud O, Donohoe D, Bultman S, Fitch M, Riiff T, Hellerstein M, et al. PPARα inhibition modulates multiple reprogrammed metabolic pathways in kidney cancer and attenuates tumor growth. Am J Physiol Cell Physiol. 2015;308(11):C890–8.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Chang WH, Lai AG. The pan-cancer mutational landscape of the PPAR pathway reveals universal patterns of dysregulated metabolism and interactions with tumor immunity and hypoxia. Ann N Y Acad Sci. 2019;1448(1):65–82.

    Article  CAS  PubMed  Google Scholar 

  50. Tong M, Jun T, Nie Y, Hao J, Fan D. The role of the slit/robo signaling pathway. J Cancer. 2019;10(12):2694–705.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Jiang Z, Liang G, Xiao Y, Qin T, Chen X, Wu E, et al. Targeting the SLIT/ROBO pathway in tumor progression: molecular mechanisms and therapeutic perspectives. Ther Adv Med Oncol. 2019;11:1758835919855238.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by the Science and Technology Plan Project of Jiangxi Provincial Health and Health Commission (202130209).

Author information

Author notes

  1. Yifu Liu, Xiaofeng Cheng, Ping Xi made equal contributions to the study.

Authors and Affiliations

  1. Department of Urology, The First Affiliated Hospital of Nanchang University, Nanchang, 330006, Jiangxi, China

    Yifu Liu, Xiaofeng Cheng, Ping Xi, Zhicheng Zhang, Ting Sun & Binbin Gong

  2. Jiangxi Institute of Urology, Nanchang, 330006, Jiangxi, China

    Yifu Liu, Xiaofeng Cheng, Ping Xi, Zhicheng Zhang, Ting Sun & Binbin Gong

Authors
  1. Yifu Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiaofeng Cheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ping Xi
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Zhicheng Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ting Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Binbin Gong
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YFL, XFC and PX wrote the manuscript. YFL and BBG reviewed the data and manuscript. ZCZ and TS organized and analyzed data. BBG edited the manuscript and supervised the study.

Corresponding author

Correspondence to Binbin Gong.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

All authors declare that they have no competing interests related to this manuscript.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, Y., Cheng, X., Xi, P. et al. Bioinformatic analysis highlights SNHG6 as a putative prognostic biomarker for kidney renal papillary cell carcinoma. BMC Urol 23, 54 (2023). https://doi.org/10.1186/s12894-023-01218-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12894-023-01218-5

Keywords

BMC Urology

ISSN: 1471-2490

Contact us