Next Article in Journal
Pan-Genomic Sequencing Reveals Actionable CDKN2A/2B Deletions and Kataegis in Anaplastic Thyroid Carcinoma
Next Article in Special Issue
Identification of Prognostic Markers of Gynecologic Cancers Utilizing Patient-Derived Xenograft Mouse Models
Previous Article in Journal
Upper Limits of Downstaging for Hepatocellular Carcinoma in Liver Transplantation
Previous Article in Special Issue
Pretreatment Albumin-to-Alkaline Phosphatase Ratio Is a Prognostic Marker in Lung Cancer Patients: A Registry-Based Study of 7077 Lung Cancer Patients
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 24
10.3390/cancers13246338
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clinicopathological and Molecular Prognostic Classifier for Intermediate/High-Risk Clear Cell Renal Cell Carcinoma

by
Fiorella L. Roldán
1,
Juan J. Lozano
2,
Mercedes Ingelmo-Torres
1,
Raquel Carrasco
1,
Esther Díaz
1,
Miguel Ramirez-Backhaus
3,
José Rubio
3,
Oscar Reig
4,
Antonio Alcaraz
1,
Lourdes Mengual
1,5,*,† and
Laura Izquierdo
1,†
1
Department and Laboratory of Urology, Hospital Clínic, Institut d’Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Universitat de Barcelona, 08036 Barcelona, Spain
2
Bioinformatics Platform, Centro de Investigación Biomédica en Red de Enfermedades Hepáticas y Digestivas (CIBERehd), Hospital Clinic, 08036 Barcelona, Spain
3
Department of Urology, Oncologic Institute of Valencia, 46009 Valencia, Spain
4
Translational Genomics and Targeted Therapeutics in Solid Tumors, August Pi i Sunyer Biomedical Research Institute (IDIBAPS) and Medical Oncology Department, Hospital Clínic de Barcelona, 08036 Barcelona, Spain
5
Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, University of Barcelona, 08036 Barcelona, Spain
*
Author to whom correspondence should be addressed.
†
Both senior authors.
Cancers 2021, 13(24), 6338; https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13246338
Submission received: 29 October 2021 / Revised: 13 December 2021 / Accepted: 14 December 2021 / Published: 17 December 2021
(This article belongs to the Collection The Biomarkers for the Diagnosis and Prognosis in Cancer)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

In this report, we identified biomarkers for tumor progression from tissue samples of intermediate/high-risk ccRCC. Using the molecular findings and the clinical data, we developed an improved prognostic model which could help to provide better individualized management recommendations.

Abstract

The probability of tumor progression in intermediate/high-risk clear cell renal cell carcinoma (ccRCC) is highly variable, underlining the lack of predictive accuracy of the current clinicopathological factors. To develop an accurate prognostic classifier for these patients, we analyzed global gene expression patterns in 13 tissue samples from progressive and non-progressive ccRCC using Illumina Hi-seq 4000. Expression levels of 22 selected differentially expressed genes (DEG) were assessed by nCounter analysis in an independent series of 71 ccRCCs. A clinicopathological-molecular model for predicting tumor progression was developed and in silico validated in a total of 202 ccRCC patients using the TCGA cohort. A total of 1202 DEGs were found between progressive and non-progressive intermediate/high-risk ccRCC in RNAseq analysis, and seven of the 22 DEGs selected were validated by nCounter. Expression of HS6ST2, pT stage, tumor size, and ISUP grade were found to be independent prognostic factors for tumor progression. A risk score generated using these variables was able to distinguish patients at higher risk of tumor progression (HR 7.27; p < 0.001), consistent with the results obtained from the TCGA cohort (HR 2.74; p < 0.002). In summary, a combined prognostic algorithm was successfully developed and validated. This model may aid physicians to select high-risk patients for adjuvant therapy.

1. Introduction

Clear cell renal cell carcinoma (ccRCC) is the most frequent renal tumor, accounting for 80–90% of cases, and has the greatest malignant potential of all renal cell carcinoma subtypes. The mainstay treatment for non-metastatic ccRCC is partial/radical nephrectomy. Despite surgical treatment, approximately between 25 to 50% of these patients will develop local relapse or distant metastases during follow-up [1,2,3].
Several prognostic algorithms have been designed to quantify the likelihood of developing disease progression [4,5]. Pathological stage, tumor size, and the Fuhrman/ International Society of Urological Pathology (ISUP) grading system appear to be the most significant prognostic factors [6,7]. Although current risk models including these variables have managed to classify patients into low-, intermediate-, and high-risk of progression, these are insufficient to accurately predict tumor aggressiveness and prognosis at the individual patient-level as tumor biology might not be entirely considered.
Gene expression profiles have been demonstrated to provide valuable prognostic information in several cancer types [8,9], including ccRCC [10,11]. However, as far as we know, none of the proposed classifiers for ccRCC are currently used in clinical practice, nor validated in intermediate/high-risk patients. These patients exhibit a metastatic potential of over 30% [4] and a high mortality rate of between 20–50% at 5 years [12], and yet no adjuvant strategies are recommended by the European Guidelines [13]. Several clinical trials for adjuvant systemic treatment have already been published. However, most of the disappointing results are due to the poor patient selection criteria based only on clinicopathological characteristics. Optimizing the criteria of patient selection has been revealed as fundamental for obtaining positive outcomes. Hence, it is crucial to identify those patients at highest risk of progression in this specific subset of patients, since they may benefit from a better medical management.
Here, we examined gene expression profiles in intermediate/high-risk ccRCC to identify prognostic biomarkers and develop a combined prognostic algorithm, including clinicopathological features and molecular biomarkers, to better predict the potential risk of recurrence after surgery.

2. Results

2.1. Clinical Features of the Cohort

The clinicopathological characteristics of the intermediate/high-risk ccRCC patients and their follow-up information, split by study phases, are listed in Table 1.
Median follow-up of the cohort was 110.76 mo (range 1.1–255.6). During follow-up, 31 (48.43%) patients developed tumor progression. All of them had distant metastasis. Median time to relapse was 12.72 months (range 1–84.84). Twenty-one (67%) of the progressive patients received adjuvant treatment, and 19 (90.5%) remained with stable disease. Thirteen patients (20.3%) died from ccRCC. None of the non-progressive patients received adjuvant treatment. Thirteen patients (20.3%) died from ccRCC.

2.2. Biomarker Discovery Phase

Overall, 1202 genes were identified as differentially expressed between progressive and non-progressive intermediate/high-risk ccRCC patients. Of these, 591 were downregulated and 611 upregulated in progressive compared with non-progressive cases. A heat map based on the most DEGs between the two groups of ccRCC patients is shown in Figure 1A.
Gene set enrichment analysis (GSEA) based on Hallmark, KEGG, and Reactome databases identified that DEGs were positively enriched in pathways, such as the epithelial-mesenchymal transition (EMT), ECM proteoglycans, non-integrin membrane ECM interactions, extracellular matrix organization, MYC targets V1, and PID syndecan 1 pathway, among others (Figure 1B). A complete list of enriched pathways and their target genes is available in Supplementary Materials (Table S1).

2.3. Biomarker Validation Phase

Twenty-two DEGs selected from the previous phase were analyzed by nCounter in an independent cohort of 64 intermediate/high-risk ccRCC samples. Seven (31.8%) genes remained significantly differentially expressed: DUOX1, HS6ST2, KRT20, RCOR2, SFN, SSC4D, and WNT9A. All were upregulated in progressive patients. According to IPA, these seven validated DEGs were enriched in organismal injury and abnormalities, renal and urological disease network, and in malignant genitourinary solid tumor diseases, among others. Significant IPA canonical pathways are depicted in Figure S1. The generated network by GeneMANIA shows that there are no direct interactions between the seven significant DEG, although some of them show co-expression (Figure S2).

2.4. Survival Analyses

Univariate and multivariate Cox regression analysis of the seven validated DEGs and the six clinicopathological variables showed that HS6ST2, pathological stage, tumor size, and ISUP grade were independent prognostic factors of tumor progression (Table 2).
Using IPA we found that HS6ST2 interacts with several molecules involved in different cancer pathways, as illustrated in Figure 2. Moreover, HS6ST2 was found to be associated with formation, activation, sprouting, and tubulation of vascular endothelial cells. Tumor development and angiogenesis are part of the five top organismal processes related to HS6ST2.

2.5. Classifier Development Phase

The RS for disease progression was calculated for each patient according to a mathematical algorithm containing HS6ST2 expression values, pathological stage, ISUP grade, and tumor size. A ROC analysis of this combined gene expression-clinicopathological model was performed, and the optimal Youden’s index cut-off value (0.648) was used to classify patients into high- and low-risk groups for tumor progression (Area Under the Curve [AUC] 0.838). Kaplan Meier curve of the combined classifier generated using the selected threshold was able to discriminate two groups with a significantly different probability of tumor progression (hazard ratio (HR) 7.27; p < 0.0001) (Figure 3A). Remarkably, the predictive combined model outperformed disease progression prediction of clinicopathological variables, as well as of HS6ST2 gene expression (Figure S3).

2.6. Classifier In Silico Validation

We evaluated the predictive capability for disease progression of our combined model in 202 intermediate/high-risk ccRCC samples obtained from the TCGA cohort. A RS for disease progression was computed for each patient. A ROC analysis was performed, and the Youden’s index value (1.106) was used as a cut-off. The AUC value of the model was 0.678, showing greater discriminatory ability than clinicopathological variables and gene expression alone (AUC pT = 0.620; AUC tumor size = 0.669; AUC ISUP grade = 0.556; AUC HS6ST2 = 0.642). The Kaplan Meier curve revealed that patients with a high-risk score have shorter disease-free survival (HR 2.74; p < 0.0002) (Figure 3B). Moreover, we carried out a Kaplan Meier sub-analysis using the RS of our combined model to evaluate its performance in the intermediate and high-risk ccRCC samples separately. The combined classifier was found to be statistically significant in the subset of high-risk ccRCC TCGA cohort (HR 2.40; p = 0.015). (Figure S4).

3. Discussion

Current risk stratification models for non-metastatic ccRCC are based mainly on clinicopathological factors. According to the SSIGN score [5], intermediate and high-risk ccRCC have a disease progression rate of 30 and 70% at 5 years, respectively, and two-thirds of these patients relapse during the first year of follow-up [3]. Several randomized clinical trials (RCT) [14,15] have investigated the benefits of tyrosine kinase inhibitor-based adjuvant therapy in these groups of patients to improve disease progression-free survival (DFS). However, results from these RCT are controversial, as only one of them [14] was positive for its primary endpoint (DFS). The most accepted explanation is that intermediate/high-risk ccRCC involves a heterogeneous group of tumors that do not share neither clinicopathological nor molecular characteristics [16]. Recently, interim data from the first phase-III adjuvant immunotherapy trial [17] has shown an improvement in DFS using an anti-PD-1 agent; however, to date, several adjuvant immunotherapy trials are ongoing, and the dissimilar patient selection criteria will most probably translate into important survival differences [18]. This emphasizes the critical need for more reliable and individualized prognostic markers to select ccRCC patients at the highest risk of recurrence to receive adjuvant treatment or to avoid overtreatment and toxic side effects in those with low progression risk.
The identification of molecular markers and the development of multigene classifiers have already been shown to considerably increase the predictive power of the clinical parameters. In this study, we determined a set of DEGs related to progression in intermediate/high-risk ccRCC using RNAseq. GSEA corroborates that the DEGs identified participate mainly in canonical pathways related to carcinogenetic, cell-cycle, and metabolic processes. To the best of our knowledge, our study is the first to report differential gene expression patterns between progressive and non-progressive intermediate/high-risk ccRCC. Validation of a subset of DEGs in an independent series allowed us to select potential prognostic biomarkers.
Our study indicates that HS6ST2 may serve as a useful prognostic marker among intermediate/high-risk ccRCCs to predict disease progression. Furthermore, we developed and validated, in an independent cohort, a model including HS6ST2 expression levels and three clinicopathological factors which significantly improves the prognostic accuracy upon clinicopathological characteristics for discriminating patients at higher risk of relapse. Our prognostic model is able to classify intermediate/high-risk ccRCC patients into two groups of patients with a different probability of tumor progression which may help physicians to tailor disease management, individualize follow-up, and serve as a tool to improve patient selection for adjuvant treatment in clinical trials. Unfortunately, most of the reported ccRCC gene-based classifiers include all clinical stages of ccRCC; thus, the classifiers’ performance cannot be comparable [10,11]. Only Wu et al. [19], using mRNA expression data from GSE53757 and TCGA databases, evaluated stage III ccRCC prognosis. However, contrary to our study, they analyzed gene expression differences between normal tissue and stage III ccRCC. Furthermore, their model was focused on overall survival, not on disease progression, as in ours.
HS6ST2 (heparan sulfate D-glucosaminyl 6-O-sulfotransferase-2) encodes an enzyme that catalyzes the transfer of sulfate groups in heparan sulfate proteoglycans (HSPGs). HSPGs participate in the regulation of numerous signaling pathways by interacting with various heparin-binding ligands and activating cytokines to influence cell growth, differentiation, adhesion, and migration [20,21]. In vitro and in vivo evidence indicates that HS6ST2 is essential for vascular development [22,23] and plays an important role in tumor angiogenesis due to its interactions with several angiogenic growth factors, including FGF2, VEGF, IL8, and IL6, among others [24,25]. Furthermore, evidence suggests that HS6ST2 is also associated with EMT by interacting with TGF-beta, HIF-1, and estrogen/GPER pathways [24,26,27,28,29], activation of T lymphocytes, and indirectly related with the PD-1, PDL-1 cancer immunotherapy pathway. Interesting, immune-oncology checkpoint inhibitors are being investigated as a potential first line adjuvant therapy after nephrectomy in non-metastatic ccRCC. Numerous studies have described that inhibition of HS6ST2 in tumor cells impairs cell migration, invasion, and tumor angiogenesis and may reverse EMT [20]. Unsurprisingly, reduction of HS6ST2 expression has been investigated as a potential target for future therapies [30].
In line with our findings, overexpression of HS6ST2 has been found to be a poor prognostic factor in several malignant tumors [20,24,30,31,32,33]. Interestingly, HS6ST2 has been recently described as part of a novel 10 glycolysis-related gene signature to predict overall survival in ccRCC [34]. Liep et al. [35] found that overexpression of miR-145-5p and miR-141-3p could inhibit the migration and invasion of RCC cells by decreasing HS6ST2 expression in cellular transfection experiments. Although more studies are required to better comprehend the involvement of this gene in cancer biology, its relationship with tumor development and angiogenesis have important clinical implications. VEGF-target therapies have been the first-line standard of care for metastatic ccRCC and still have a relevant role in metastatic ccRCC as a partner with anti-PD1/L1 immunotherapy in first-line treatment or in monotherapy at progression, due to the natural chemoresistance and radioresistance shown by tumoral renal cells.
Our study has multiple strengths. The first is that we assessed intermediate/high-risk ccRCC patients in a multicenter study. Despite it being well known that these patients have a higher risk of recurrence, our classifier may contribute to determine which of them could most benefit from receiving the new targeted adjuvant treatment. Secondly, RNA was isolated from FFPE ccRCC tissue obtained in routine practice, and the technology used to quantify gene expression is highly cost-effective, facilitating its implementation in the clinical setting. Thirdly, our non-progressive ccRCC patients had a long-term follow-up to avoid misclassification. Lastly, the researchers participating in this study were blinded to all clinical information, and gene expression was matched to clinical data only after all patient cases had been processed.
Nevertheless, we acknowledge some study limitations. First, the retrospective design and the small sample used in the validation phase. Secondly, the biomarker selected may not apply to low-risk tumors and should not be generalizable to all patients with ccRCC. Finally, external validation in the TCGA ccRCC samples revealed that our combined classifier has a good performance predicting disease progression; nevertheless, further experimental validation in larger cohorts is required.

4. Materials and Methods

4.1. Patients

A multicenter study in which a total of 84 non-metastatic intermediate/high-risk ccRCCs, as defined by the Mayo Clinic Stage, Size, Grade, and Necrosis (SSIGN) score [5], who underwent partial or radical nephrectomy between 2000 and 2012 in two different centers (Hospital Clinic of Barcelona, Barcelona, Spain, and Oncologic Institute of Valencia, Valencia, Spain), were retrospectively included.
This study was split into four-stage approach with an initial biomarker discovery phase, a biomarker validation phase, a classifier development phase, and, lastly, the in silico external validation using the TCGA cohort [36,37] (Figure 4). Initial discovery phase included 13 ccRCC patients from Hospital Clinic of Barcelona, six progressive and seven non-progressive. Biomarker validation phase comprises 71 ccRCC patients from the Oncologic Institute of Valencia (IVO). From this, 64 ccRCC patients were eventually included in the study, 31 progressive and 33 non-progressive (Table 1). The remaining seven patients were discarded because of low count values in Nanostring analysis. Tissue samples were obtained under institutional review board-approved protocols (HBC/2016/0333 and 2017-59-BIOBANCOFIVO-15-2017).
All patients were followed-up postoperatively according to the European Urology guidelines [13]. Briefly, CT scans were performed at the 3rd month after surgery, every 6-months for the first 3 years, annually until 5 years, and biannually thereafter. Tumors were considered as progressive when local relapse or distant metastasis were developed during the follow-up. Non-progressive patients had a minimum follow-up of 10 years.

4.2. Tissue Specimens and RNA Isolation

Formalin-fixed paraffin-embedded (FFPE) tissue blocks were reviewed. RNA was isolated from FFPE specimens (total thickness 80 µm) using the RecoverAll™ Total Nucleic Acid Isolation kit for FFPE (Ambion, Inc., Austin, TX, USA) following manufacturer’s instructions. RNA was quantified by spectrophotometric analysis at 260 nm (NanoDrop Technologies, Wilmington, DE, USA).

4.3. Biomarker Discovery Phase—RNA Sequencing

4.3.1. Library Preparation and Sequencing Method

RNA from 13 selected ccRCC samples was processed for library preparation using the TruSeq® RNA Access Library Preparation Kit (Illumina, San Diego, CA, USA) that allows generating libraries starting from degraded RNA. Briefly, cDNA strands were synthetized from input RNA in order to be adaptor-tagged, labeled, and amplified. cDNA was then pooled and enriched by a double step of probe hybridization. The enriched targets were captured by streptavidin labeled beads, cleaned up, and amplified to obtain the final multiplexed libraries. The libraries were then sequenced on an Illumina HiSeq® 4000 platform (Illumina, San Diego, CA, USA).

4.3.2. Read Alignment and Differential Gene Expression Analysis

Paired-end RNA-Seq FASTQ files were trimmed from a 3′ end to a fixed length based on the Phred quality score (trimmed if score fell below 20, with a minimum read length of 25) [38]. Trimmed RNA-seq reads were aligned to the GRCh38 reference genome with STAR [39], and gene counts were determined using quantMode GeneCounts. Trimmed reads were then aligned using STAR. We used limma-voom transformation and cyclic-loess to normalize the non-biological variability. An assessment of differential expression between groups was evaluated using moderated t-statistics [40]. Significant DEGs in progressive relative to non-progressive patients were identified based on an adjusted p value of <0.05. The heatmap and statistical analyses were performed using R statistical package (v3.3.2). Gene set enrichment analysis (GSEA) was performed using GSEA2-2.2.0 software for testing specific gene sets based on Hallmark, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Reactome pathway databases [41]. RNAseq files and clinical information were deposited into Gene Expression Omnibus (GEO) with the accession number GSE175648.

4.4. Biomarker Validation Phase

4.4.1. nCounter Elements Gene Expression Analysis

The 22 differentially expressed genes (DEGs) with a higher fold change were selected from the discovery phase and validated in an independent cohort of 71 ccRCC patients. nCounter Elements Analysis System (NanoString Technologies, Seattle, WA, USA) was used to quantify gene expression from these 22 targets and two housekeeping (PPIA and TBP) genes [42]. Briefly, probes were hybridized with the Elements TagSet and 250 ng of RNA samples for 21 h at 67 °C in a Thermal Cycler. Thereafter, samples were purified in the nCounter Prep Station to remove unligated probes. Expression data were collected using the nCounter Digital analyzer. Those counts with values ≤ 10 were excluded from the analysis, and those genes (n = 14) expressed in less than 80% of the samples were discarded for further examinations. Furthermore, those samples (n = 7) with less than 80% of expressed genes were excluded from further analysis. Empty input values for each gene were inputted to the minimum count value for that gene. Gene expression normalization was performed using the PPIA housekeeping gene since TBP was discarded because of its lack of expression in more than 20% of samples.

4.4.2. Survival Analysis

Cox stepwise regression analysis was used on the established clinical prognostic factors of ccRCC (pT stage, pN status, ISUP grade, tumor size, necrosis, vascular, and perirenal fat invasion) and the seven DEGs to investigate their influence on tumor progression. Statistical significance was defined at a p value of 0.05.

4.5. Classifier Development Phase

After establishing the multivariate model, a risk score (RS) for the variables of the model was calculated for each patient. RS was subjected to a Receiver Operating Characteristics (ROC) curve analysis to choose the most appropriate threshold for predicting tumor progression. Thereafter, Kaplan-Meier curves were generated using the selected cut-off point and compared according to the log-rank test.

4.6. Pathway Enrichment Analysis

Ingenuity Pathway Analysis (IPA) software was used to identify interactions and networks between significant DEGs, possible altered canonical pathways, regulators, diseases, and functions based on direct/indirect and experimental targets.

4.7. Classifier In Silico Validation

The Cancer Genome Atlas (TCGA) gene expression dataset was obtained from portal (https://firebrowse.org) (accessed on 8 October 2021) [36], and TCGA clinical data were downloaded from the portal (https://0-www-sciencedirect-com.brum.beds.ac.uk/science/article/pii/S0092867418302290?via%3Dihub#app2) (accessed on 8 October 2021) [37]. After selecting samples matching our selection criteria and excluding patients without survival status or missing clinical data, a total of 202 ccRCCs were selected. Level 3 RNAseq expression data from ccRCC samples by the Illumina HiSeq2000 RNA sequencing platform, and corresponding clinical data were downloaded from the aforementioned websites.

5. Conclusions

A novel prognostic algorithm based on gene expression and clinicopathological factors was successfully developed and validated in the TCGA cohort. This model is able to identify a subset of intermediate/high-risk ccRCC patients with a higher risk of tumor progression and may contribute to identify those patients that would most benefit from closer surveillance and adjuvant treatment.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/cancers13246338/s1, Figure S1. Significant canonical pathways and overlapping network for the seven genes validated, Figure S2. Gene-gene interaction network for the 7 significant differentially expressed genes, Figure S3. Performance of the combined classifier. ROC curves for predicting tumor progression after surgery according to the combined classifier, gene expression, and clinicopathological variables, Figure S4. Kaplan Meier survival analyses for tumor progression according to the combined classifier (A) in high-risk TCGA Cohort and (B) in the intermediate-risk TCGA cohort, Table S1: Enriched pathways and their target genes from the DEG in the RNAseq analysis.

Author Contributions

Conceptualization, A.A., L.M. and L.I.; Methodology, M.R.-B., J.R., O.R. and F.L.R.; Software and validation, F.L.R. and J.J.L.; Formal analysis, F.L.R., L.I., L.M. and A.A.; Investigation, F.L.R., L.I. and L.M.; Resources E.D., M.I.-T. and R.C.; data curation J.J.L. and F.L.R.; Writing—original draft preparation, F.L.R.; Writing—review and editing, all authors contributed to review and editing the writing; Visualization, F.L.R. and L.M.; Supervision, L.M., L.I. and A.A.; Project administration, E.D., M.I.-T. and R.C.; Funding acquisition, F.L.R. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by an Emili Letang grant from Hospital Clínic de Barcelona to FLR/2019.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Hospital clinic of Barcelona (protocol code HBC/2016/0333 and date of approval was 24 January 2019) and Oncologic Institute of Valencia (protocol code 2017-59-BIOBANCOFIVO-15-2017 and date of approval was 1 May 2017).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available in the article and supplementary material. RNAseq files and clinical information were deposited into Gene Expression Omnibus (GEO) with the accession number GSE175648. Further details can be obtained on request from the corresponding author.

Acknowledgments

The authors thank Macrogen, Inc., for RNAseq analysis and Débora Martínez from Medical Oncology Department of Hospital Clinic of Barcelona for technical support with the nCounter platform. We are indebted to the IDIBAPS Biobank for sample and fecha procurement. We thank Azucena Salas and Lluís Revilla for IPA software support. This work was developed at the building Centre de Recerca Biomèdica Cellex, Barcelona.

Conflicts of Interest

The authors have no conflict of interest to disclose.

References

  1. Ghatalia, P.; Gordetsky, J.; Kuo, F.; Dulaimi, E.; Cai, K.Q.; Devarajan, K.; Bae, S.; Naik, G.; Chan, T.A.; Uzzo, R.; et al. Prognostic impact of immune gene expression signature and tumor infiltrating immune cells in localized clear cell renal cell carcinoma. J. Immunother. Cancer 2019, 7, 139. [Google Scholar] [CrossRef] [Green Version]
  2. Sharma, T.; Tajzler, C.; Kapoor, A. Is There a Role for Adjuvant Therapy After Surgery in “High Risk for Recurrence” Kidney Cancer? An Update on Current Concepts. Curr. Oncol. 2018, 25, 444–453. [Google Scholar] [CrossRef] [Green Version]
  3. Janzen, N.K.; Kim, H.L.; Figlin, R.A.; Belldegrun, A.S. Surveillance after radical or partial nephrectomy for localized renal cell carcinoma and management of recurrent disease. Urol. Clin. N. Am. 2003, 30, 843–852. [Google Scholar] [CrossRef]
  4. Leibovich, B.C.; Blute, M.L.; Cheville, J.C.; Lohse, C.M.; Frank, I.; Kwon, E.D.; Weaver, A.L.; Parker, A.S.; Zincke, H. Prediction of progression after radical nephrectomy for patients with clear cell renal cell carcinoma: A stratification tool for prospective clinical trials. Cancer 2003, 97, 1663–1671. [Google Scholar] [CrossRef]
  5. Frank, I.G.O.R.; Blute, M.L.; Cheville, J.C.; Lohse, C.M.; Weaver, A.L.; Zincke, H. An outcome prediction model for patients with clear cell renal cell carcinoma treated with radical nephrectomy based on tumor stage, size, grade and necrosis: The SSIGN score. J. Urol. 2002, 168, 2395–2400. [Google Scholar] [CrossRef]
  6. Ficarra, V.; Schips, L.; Guillè, F.; Li, G.; De La Taille, A.; Galetti, T.P.; Cindolo, L.; Novara, G.; Zigeuner, R.E.; Bratti, E.; et al. Multiinstitutional European validation of the 2002 TNM staging system in conventional and papillary localized renal cell carcinoma. Cancer 2005, 104, 968–974. [Google Scholar] [CrossRef] [PubMed]
  7. Fuhrman, S.A.; Lasky, L.C.; Limas, C. Prognostic significance of morphologic parameters in renal cell carcinoma. Am. J. Surg. Pathol. 1982, 6, 655–664. [Google Scholar] [CrossRef] [PubMed]
  8. Cuzick, J.; Swanson, G.P.; Fisher, G.; Brothman, A.R.; Berney, D.M.; Reid, J.E.; Mesher, D.; Speights, V.O.; Stankiewicz, E.; Foster, C.S.; et al. Prognostic value of an RNA expression signature derived from cell cycle proliferation genes in patients with prostate cancer: A retrospective study. Lancet Oncol. 2011, 12, 245–255. [Google Scholar] [CrossRef] [Green Version]
  9. Kwa, M.; Makris, A.; Esteva, F.J. Clinical utility of gene-expression signatures in early stage breast cancer. Nat. Rev. Clin. Oncol. 2017, 14, 595–610. [Google Scholar] [CrossRef] [PubMed]
  10. Rini, B.; Goddard, A.; Knezevic, D.; Maddala, T.; Zhou, M.; Aydin, H.; Campbell, S.; Elson, P.; Koscielny, S.; Lopatin, M.; et al. A 16-gene assay to predict recurrence after surgery in localised renal cell carcinoma: Development and validation studies. Lancet Oncol. 2015, 16, 676–685. [Google Scholar] [CrossRef]
  11. Brooks, S.A.; Brannon, A.R.; Parker, J.S.; Fisher, J.C.; Sen, O.; Kattan, M.W.; Hakimi, A.A.; Hsieh, J.J.; Choueiri, T.K.; Tamboli, P.; et al. ClearCode34: A Prognostic Risk Predictor for Localized Clear Cell Renal Cell Carcinoma. Eur. Urol. 2014, 66, 77–84. [Google Scholar] [CrossRef] [Green Version]
  12. Patard, J.-J.; Kim, H.L.; Lam, J.S.; Dorey, F.J.; Pantuck, A.J.; Zisman, A.; Ficarra, V.; Han, K.-R.; Cindolo, L.; De La Taille, A.; et al. Use of the University of California Los Angeles Integrated Staging System to Predict Survival in Renal Cell Carcinoma: An International Multicenter Study. J. Clin. Oncol. 2004, 22, 3316–3322. [Google Scholar] [CrossRef]
  13. Ljungberg, B.; Albiges, L.; Bensalah, K.; Bex, A.; Giles, R.H.; Hora, M. EAU Guidelines. Edn. presented at the EAU Annual Congress Amsterdam 2020. Eur. Urol. 2020, 67, 913–924. [Google Scholar] [CrossRef] [PubMed]
  14. Ravaud, A.; Motzer, R.J.; Pandha, H.S.; George, D.J.; Pantuck, A.J.; Patel, A.; Chang, Y.-H.; Escudier, B.; Donskov, F.; Magheli, A.; et al. Adjuvant Sunitinib in High-Risk Renal-Cell Carcinoma after Nephrectomy. N. Engl. J. Med. 2016, 375, 2246–2254. [Google Scholar] [CrossRef] [PubMed]
  15. Motzer, R.J.; Haas, N.B.; Donskov, F.; Gross-Goupil, M.; Varlamov, S.; Kopyltsov, E.; Lee, J.L.; Melichar, B.; Rini, B.I.; Choueiri, T.K.; et al. Randomized Phase III Trial of Adjuvant Pazopanib Versus Placebo after Nephrectomy in Patients with Localized or Locally Advanced Renal Cell Carcinoma. J. Clin. Oncol. 2017, 35, 3916–3923. [Google Scholar] [CrossRef]
  16. Ricketts, C.J.; De Cubas, A.A.; Fan, H.; Smith, C.C.; Lang, M.; Reznik, E.; Bowlby, R.; Gibb, E.A.; Akbani, R.; Beroukhim, R.; et al. The Cancer Genome Atlas Comprehensive Molecular Characterization of Renal Cell Carcinoma. Cell Rep. 2018, 23, 313–326.e5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. ASCO 2021: Results from the KEYNOTE-564 Study of Adjuvant Pembrolizumab in Resected Clear-Cell Renal Cell Carcinoma (RCC). Available online: https://ascopost.com/issues/june-25-2021/keynote-564-adjuvant-pembrolizumab-extends-disease-free-survival-in-high-risk-renal-cell-carcinoma/ (accessed on 12 August 2021).
  18. Palumbo, C.; Mazzone, E.; Mistretta, F.A.; Knipper, S.; Perrotte, P.; Shariat, S.F.; Saad, F.; Kapoor, A.; Lattouf, J.-B.; Simeone, C.; et al. A Plea for Optimizing Selection in Current Adjuvant Immunotherapy Trials for High-Risk Nonmetastatic Renal Cell Carcinoma According to Expected Cancer-specific Mortality. Clin. Genitourin. Cancer 2020, 18, 314–321. [Google Scholar] [CrossRef]
  19. Wu, J.; Jin, S.; Gu, W.; Wan, F.; Zhang, H.; Shi, G.; Qu, Y.; Ye, D. Construction and Validation of a 9-Gene Signature for Predicting Prognosis in Stage III Clear Cell Renal Cell Carcinoma. Front. Oncol. 2019, 9, 152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Jin, Y.; He, J.; Du, J.; Zhang, R.-X.; Yao, H.-B.; Shao, Q.-S. Overexpression of HS6ST2 is associated with poor prognosis in patients with gastric cancer. Oncol. Lett. 2017, 14, 6191–6197. [Google Scholar] [CrossRef]
  21. Sarrazin, S.; Lamanna, W.C.; Esko, J.D. Heparan sulfate proteoglycans. Cold Spring Harb. Perspect. Biol. 2011, 3, a004952. [Google Scholar] [CrossRef] [Green Version]
  22. Lundin, L.; Larsson, H.; Kreuger, J.; Kanda, S.; Lindahl, U.; Salmivirta, M.; Claesson-Welsh, L. Selectively Desulfated Heparin Inhibits Fibroblast Growth Factor-induced Mitogenicity and Angiogenesis. J. Biol. Chem. 2000, 275, 24653–24660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chen, E.; Stringer, S.E.; Rusch, M.A.; Selleck, S.B.; Ekker, S.C. A unique role for 6-O sulfation modification in zebrafish vascular development. Dev. Biol. 2005, 284, 364–376. [Google Scholar] [CrossRef] [PubMed]
  24. Song, K.; Li, Q.; Peng, Y.-B.; Li, J.; Ding, K.; Chen, L.-J.; Shao, C.-H.; Zhang, L.-J.; Li, P. Silencing of hHS6ST2 inhibits progression of pancreatic cancer through inhibition of Notch signalling. Biochem. J. 2011, 436, 271–282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Cole, C.L.; Rushton, G.; Jayson, G.C.; Avizienyte, E. Ovarian Cancer Cell Heparan Sulfate 6-O-Sulfotransferases Regulate an Angiogenic Program Induced by Heparin-Binding Epidermal Growth Factor (EGF)-Like Growth Factor/EGF Receptor Signaling. J. Biol. Chem. 2014, 289, 10488–10501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Pupo, M.; Pisano, A.; Abonante, S.; Maggiolini, M.; Musti, A.M. GPER activates Notch signaling in breast cancer cells and cancer-associated fibroblasts (CAFs). Int. J. Biochem. Cell Biol. 2014, 46, 56–67. [Google Scholar] [CrossRef]
  27. Sahlgren, C.; Gustafsson, M.V.; Jin, S.; Poellinger, L.; Lendahl, U. Notch signaling mediates hypoxia-induced tumor cell migration and invasion. Proc. Natl. Acad. Sci. USA 2008, 105, 6392–6397. [Google Scholar] [CrossRef] [Green Version]
  28. Zhou, J.; Jain, S.; Azad, A.K.; Xu, X.; Yu, H.C.; Xu, Z.; Godbout, R.; Fu, Y. Notch and TGFβ form a positive regulatory loop and regulate EMT in epithelial ovarian cancer cells. Cell. Signal. 2016, 28, 838–849. [Google Scholar] [CrossRef]
  29. De Francesco, E.M.; Maggiolini, M.; Musti, A.M. Crosstalk between notch, HIF-1α and GPER in breast cancer EMT. Int. J. Mol. Sci. 2018, 19, 2011. [Google Scholar] [CrossRef] [Green Version]
  30. Pollari, S.; Käkönen, R.S.; Mohammad, K.S.; Rissanen, J.P.; Halleen, J.M.; Wärri, A.; Nissinen, L.; Pihlavisto, M.; Marjamäki, A.; Perälä, M.; et al. Heparin-like Polysaccharides Reduce Osteolytic Bone Destruction and Tumor Growth in a Mouse Model of Breast Cancer Bone Metastasis. Mol. Cancer Res. 2012, 10, 597–604. [Google Scholar] [CrossRef] [Green Version]
  31. Hatabe, S.; Kimura, H.; Arao, T.; Kato, H.; Hayashi, H.; Nagai, T.; Matsumoto, K.; De Velasco, M.; Fujita, Y.; Yamanouchi, G.; et al. Overexpression of heparan sulfate 6-O-sulfotransferase-2 in colorectal cancer. Mol. Clin. Oncol. 2013, 1, 845–850. [Google Scholar] [CrossRef] [Green Version]
  32. Backen, A.C.; Cole, C.L.; Lau, S.C.; Clamp, A.R.; McVey, R.; Gallagher, J.T.; Jayson, G.C. Heparan sulphate synthetic and editing enzymes in ovarian cancer. Br. J. Cancer 2007, 96, 1544–1548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Schulten, H.-J.; Al-Mansouri, Z.; Baghallab, I.; Bagatian, N.; Subhi, O.; Karim, S.; Al-Aradati, H.; Al-Mutawa, A.; Johary, A.; Meccawy, A.A.; et al. Comparison of microarray expression profiles between follicular variant of papillary thyroid carcinomas and follicular adenomas of the thyroid. BMC Genom. 2015, 16, S7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Xing, Q.; Zeng, T.; Liu, S.; Cheng, H.; Ma, L.; Wang, Y. A novel 10 glycolysis-related genes signature could predict overall survival for clear cell renal cell carcinoma. BMC Cancer 2021, 21, 381. [Google Scholar] [CrossRef]
  35. Liep, J.; Kilic, E.; Meyer, H.A.; 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, e0157801. [Google Scholar] [CrossRef] [PubMed]
  36. FireBrowse. Available online: http://firebrowse.org/ (accessed on 8 October 2021).
  37. Liu, J.; Lichtenberg, T.; Hoadley, K.A.; Poisson, L.M.; Lazar, A.J.; Cherniack, A.D.; Kovatich, A.J.; Benz, C.C.; Levine, D.A.; Lee, A.V.; et al. An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell 2018, 173, 400–416.e11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Babraham Bioinformatics-Trim Galore! Available online: https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/ (accessed on 15 May 2019).
  39. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  40. Ritchie, M.E.; Phipson, B.; Wu, D.I.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  41. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  42. Bjerregaard, H.; Pedersen, S.; Kristensen, S.R.; Marcussen, N. Reference Genes for Gene Expression Analysis by Real-time Reverse Transcription Polymerase Chain Reaction of Renal Cell Carcinoma. Diagn. Mol. Pathol. 2011, 20, 212–217. [Google Scholar] [CrossRef] [PubMed]
Cancers 13 06338 g001 550
Figure 1. Differentially expressed genes in the discovery phase. (A) Heat map displaying the 50 most DEGs between progressive and non-progressive intermediate/high-risk ccRCC patients. Red pixels correspond to up-regulated genes, whereas green pixels indicate down-regulated genes. (B) GSEA shows positive correlation of DEGs in pathways involved in tumor progression. Abbreviations: DEGs, Differentially expressed genes. GSEA, Gene set enrichment analysis.
Figure 1. Differentially expressed genes in the discovery phase. (A) Heat map displaying the 50 most DEGs between progressive and non-progressive intermediate/high-risk ccRCC patients. Red pixels correspond to up-regulated genes, whereas green pixels indicate down-regulated genes. (B) GSEA shows positive correlation of DEGs in pathways involved in tumor progression. Abbreviations: DEGs, Differentially expressed genes. GSEA, Gene set enrichment analysis.
Cancers 13 06338 g001
Cancers 13 06338 g002 550
Figure 2. Disease and functions related to the validated genes and HS6ST2 significant molecular interactions and involved pathways. Abbreviations: CP, Canonical pathways.
Figure 2. Disease and functions related to the validated genes and HS6ST2 significant molecular interactions and involved pathways. Abbreviations: CP, Canonical pathways.
Cancers 13 06338 g002
Cancers 13 06338 g003 550
Figure 3. Performance of the combined classifier. Kaplan Meier survival analyses for tumor progression according to the combined classifier (A) in our cohort and (B) in the TCGA cohort.
Figure 3. Performance of the combined classifier. Kaplan Meier survival analyses for tumor progression according to the combined classifier (A) in our cohort and (B) in the TCGA cohort.
Cancers 13 06338 g003
Cancers 13 06338 g004 550
Figure 4. Study outline. Tissue samples were obtained from a total of 77 patients with intermediate/high-risk ccRCC. Samples were split into a biomarker discovery (N = 13) and validation (N = 64) phase. Transcripts differentially expressed between progressive and non-progressive intermediate/high-risk ccRCCs were first identified in the discovery phase using RNAseq. Twenty-two DEGs were selected for validation in an independent set of 64 tissue samples using nCounter Elements (Nanostring). A prognostic model was generated using gene expression and clinical data in the classifier development phase. Finally, the prognostic model was in silico validated using a TCGA cohort. Abbreviations: ccRCC, Clear cell renal cell carcinoma, DEGs, Differentially expressed genes.
Figure 4. Study outline. Tissue samples were obtained from a total of 77 patients with intermediate/high-risk ccRCC. Samples were split into a biomarker discovery (N = 13) and validation (N = 64) phase. Transcripts differentially expressed between progressive and non-progressive intermediate/high-risk ccRCCs were first identified in the discovery phase using RNAseq. Twenty-two DEGs were selected for validation in an independent set of 64 tissue samples using nCounter Elements (Nanostring). A prognostic model was generated using gene expression and clinical data in the classifier development phase. Finally, the prognostic model was in silico validated using a TCGA cohort. Abbreviations: ccRCC, Clear cell renal cell carcinoma, DEGs, Differentially expressed genes.
Cancers 13 06338 g004
Table
Table 1. Demographic and pathological characteristics of enrolled patients.
Table 1. Demographic and pathological characteristics of enrolled patients.
Clinicopathological
Characteristics		Discovery Phase Hospital Clinic Barcelona (n = 13)Validation Phase Institute Valenciano of Oncology (n = 64)
Gender
Male9 (69.2)49 (76.6)
Female4 (30.8)15 (23.4)
Age at diagnosis
(yr)		54.85 (36–81)58.6 (35–87)
Pathological
tumor size (cm)		8.2 (2.5–14)8 (3.1–24)
ISUP
ISUP 1-4 (6.3)
ISUP 22 (15.4)20 (31.3)
ISUP 36 (46.2)31 (48.4)
ISUP 45 (38.4)9 (14)
Tumor stage
pT15 (38.4)9 (14)
pT25 (38.4)16 (25)
pT32 (15.4)36 (56.3)
pT41 (7.8)3 (4.7)
N stage
N0/x11 (84.6)58 (90.6)
N12 (15.4)6 (9.4)
Perirenal fat invasion3 (23.1)40 (62.5)
Vascular invasion2 (15.4)10 (15.6)
Necrosis1 (7.8)24 (37.5)
SSIGN score
Intermediate risk7 (53.8)39 (60.9)
High risk6 (46.2)25 (39.1)
( ) Range or %.
Table
Table 2. Univariate and multivariate Cox regression analysis of statistically significant genetic and clinical variables in the validation set.
Table 2. Univariate and multivariate Cox regression analysis of statistically significant genetic and clinical variables in the validation set.
GenesUnivariateMultivariate
p95% CIHRp95% CIHR
DUOX1<0.0011.404–2.7511.965
HS6ST20.0012.084–14.1575.432<0.0012.710–14.8806.35
KRT20<0.0011.419–2.9262.037
RCOR2<0.0012.417–16.6276.340
SFN<0.0011.806–5.6473.194
SSC4D0.0011.982–14.8315.422
WNT9A0.0021.600–7.4753.459
pT Stage<0.0011.610–5.3622.9390.0161.150–4.0902.17
Tumor size<0.0011.068–1.2461.1540.0181.020–1.2301.12
ISUP0.0451.010–2.5331.5990.0211.100–3.3701.93
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Roldán, F.L.; Lozano, J.J.; Ingelmo-Torres, M.; Carrasco, R.; Díaz, E.; Ramirez-Backhaus, M.; Rubio, J.; Reig, O.; Alcaraz, A.; Mengual, L.; et al. Clinicopathological and Molecular Prognostic Classifier for Intermediate/High-Risk Clear Cell Renal Cell Carcinoma. Cancers 2021, 13, 6338. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13246338

AMA Style

Roldán FL, Lozano JJ, Ingelmo-Torres M, Carrasco R, Díaz E, Ramirez-Backhaus M, Rubio J, Reig O, Alcaraz A, Mengual L, et al. Clinicopathological and Molecular Prognostic Classifier for Intermediate/High-Risk Clear Cell Renal Cell Carcinoma. Cancers. 2021; 13(24):6338. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13246338

Chicago/Turabian Style

Roldán, Fiorella L., Juan J. Lozano, Mercedes Ingelmo-Torres, Raquel Carrasco, Esther Díaz, Miguel Ramirez-Backhaus, José Rubio, Oscar Reig, Antonio Alcaraz, Lourdes Mengual, and et al. 2021. "Clinicopathological and Molecular Prognostic Classifier for Intermediate/High-Risk Clear Cell Renal Cell Carcinoma" Cancers 13, no. 24: 6338. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13246338

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop