Next Article in Journal
Overexpression of p62/IMP2 can Promote Cell Migration in Hepatocellular Carcinoma via Activation of the Wnt/β-Catenin Pathway
Next Article in Special Issue
Prognostic and Predictive Value of PBRM1 in Clear Cell Renal Cell Carcinoma
Previous Article in Journal
PARP1 Inhibition Augments UVB-Mediated Mitochondrial Changes—Implications for UV-Induced DNA Repair and Photocarcinogenesis
Previous Article in Special Issue
GSTO1*CC Genotype (rs4925) Predicts Shorter Survival in Clear Cell Renal Cell Carcinoma Male Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Hypertonicity-Affected Genes Are Differentially Expressed in Clear Cell Renal Cell Carcinoma and Correlate with Cancer-Specific Survival

1
Dr. Margarete Fischer-Bosch Institute of Clinical Pharmacology, 70376 Stuttgart, Germany
2
University of Tübingen, 72074 Tübingen, Germany
3
Department of Medicine, Hematology and Oncology, Martin Luther University Halle-Wittenberg, 06120 Halle (Saale), Germany
4
Max Delbrück Center for Molecular Medicine (MDC), 13125 Berlin, Germany
5
Core Unit DNA–Technologien, Medizinische Fakultät, Universität Leipzig, 04103 Leipzig, Germany
6
Department of Urology, University Hospital Tübingen, 72076 Tübingen, Germany
7
German Cancer Consortium (DKTK) and German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
8
Departments of Clinical Pharmacology, Pharmacy and Biochemistry, University Tübingen, 72076 Tübingen, Germany
9
iFIT Cluster of Excellence EXC 2180, University of Tübingen, 72076 Tübingen, Germany
*
Author to whom correspondence should be addressed.
Cancers 2020, 12(1), 6; https://doi.org/10.3390/cancers12010006
Submission received: 29 November 2019 / Revised: 13 December 2019 / Accepted: 13 December 2019 / Published: 18 December 2019
(This article belongs to the Special Issue Renal Cell Carcinoma)

Abstract

:
The heterogeneity of renal cell carcinoma (RCC) subtypes reflects the cell type of origin in the nephron, with consequences for therapy and prognosis. The transcriptional cues that determine segment-specific gene expression patterns are poorly understood. We recently showed that hypertonicity in the renal medulla regulates nephron-specific gene expression. Here, we analyzed a set of 223 genes, which were identified in the present study by RNA-Seq to be differentially expressed by hypertonicity, for the prediction of cancer-specific survival (CSS). Cluster analyses of these genes showed discrimination between tumor and non-tumor samples of clear cell RCC (ccRCC). Refinement of this gene signature to a four-gene score (OSM score) through statistical analyses enabled prediction of CSS in ccRCC patients of The Cancer Genome Atlas (TCGA) (n = 436) in univariate (HR = 4.1; 95% CI: 2.78−6.07; p = 4.39 × 10−13), and multivariate analyses including primary tumor (T); regional lymph node (N); distant metastasis (M); grading (G)(p = 2.3 × 10−5). The OSM score could be validated in an independent ccRCC study (n = 52) in univariate (HR = 1.29; 95% CI = 1.05–1.59; p = 0.011) and multivariate analyses (p = 0.016). Cell culture experiments using RCC cell lines demonstrated that the expression of the tumor suppressor ELF5 could be restored by hypertonicity. The innovation of our novel gene signature is that these genes are physiologically regulated only by hypertonicity, thereby providing the possibility to be targeted for therapy.

1. Introduction

The kidney’s anatomy and histology consists of different renal cell types located at defined parts of the kidneys, reflected in the complexity of renal function. This is also reflected by the heterogeneity of renal cell carcinoma (RCC) subtypes [1]. The main subtypes are clear cell (ccRCC), papillary (pRCC) and chromophobe (chRCC) renal cell carcinoma [2]. Although several targeted therapies are currently applied, survival rates—especially for metastatic RCC—are still low and innovative treatment strategies are needed [2].
Comprehensive studies carried out by The Cancer Genome Atlas (TCGA) provided further insight into the evolution and origin of RCC, for example, by identifying gene signatures that enable discrimination between the RCC subtypes or define the cell type of origin. Moreover, it was found to be possible to predict clinical outcome in ccRCC patients based on gene expression similarity to the proximal tubule of the nephron, which is the presumed origin of ccRCC [3]. Recently, another study analyzed the impact of different gene expression profiles on RCC ontogeny [1]. The authors were able to identify gene expression programs that were specific for a distinct nephron segment and were also present in the corresponding RCC subtypes. Both studies used data based on nephron-specific gene expression patterns and were able either to improve the prediction of patient outcome or identify gene expression networks defining the origin of RCC. However, as mentioned by Lindgren et al. [1], the transcriptional cues that determine segment-specific gene expression patterns are only partly understood. We have recently shown that the unique hypertonicity in the renal inner medulla regulates kidney and nephron-specific gene expression [4]. The group of Prof. Ian Frew showed that deletion of renal expression of the tumor suppressor von Hippel–Lindau (VHL) protein altered the urine concentration capability in mice [5]. They postulate that the mice cannot build up the hyperosmotic gradient in the kidneys that is necessary for urine concentration. The transcription factor nuclear factor of activated T-cells 5 (NFAT5) is the main transcription factor activated by the hyperosmotic environment, and induces the expression of several genes [6]. A recent study in that Special Issue of Cancers showed that microRNAs that mediate metabolic reprogramming in ccRCC also target NFAT5 [7]. This was also associated with a reduced level of NFAT5 target genes in the ccRCC samples compared to solid normal tissue.
In the present study, we analyzed if the hypertonicity-affected genes were also differentially expressed in ccRCC tumor samples and normal tissue, and if these genes were associated with the clinical outcome of the patients.

2. Results and Discussion

In contrast to our initial study, where we used microarrays, here we performed RNA-Seq using primary cultured inner medullary collecting duct (IMCD) cells cultivated at 300 or 900 mosmol/kg to identify differentially expressed transcripts affected by hypertonicity (for details see Table S1). We detected significant differences between the two conditions for 355 transcripts (false discovery rate FDR < 0.05; log2 fold change (FC) >3/<−3) and there were matching human transcripts for 284 of these (223 genes) (Figure 1A and Table S1). Hierarchical clustering of the TCGA Kidney Clear Cell Carcinoma (KIRC) samples based RNA-Seq data using the top 223 hypertonicity-affected genes clearly separated the normal non-tumor tissue samples (n = 67) from the tumor samples (n = 449; Figure 1B).
Part of the genes (41) showed a log2 fold change of >3/<−3 between normal and tumor samples. Interestingly, several of the transcripts induced by hypertonicity were suppressed, and transcripts suppressed by hypertonicity were induced in the tumor samples compared to normal samples (Table S1).
The effect of hyper-osmolality on gene expression can be reverted by hypo osmotic switch [4]. For example one of the hypertonicity-induced transcripts (0 fragments per kilobase of transcript per million mapped reads (FPKM) at 300 vs. 75 FPKM at 900 mosmol/kg, see Table S1) was the E74-like ETS transcription factor 5 (ELF5). ELF5 has been described as a tumor suppressor in RCC and is more or less absent in tumor samples (Table S1) [8]. Since ELF5 has an important role as a tumor suppressor, we next asked whether it is possible to induce its expression in a ccRCC cell line by hyperosmolality. To test this, we used the established ccRCC cell line 786-0, and the same cell line that ectopically expresses WT-VHL (786-0-VHL). Both were cultivated either under isotonic (300 mosmol/kg) or for different periods of time under hyperosmotic (600 mosmol/kg) conditions. Indeed, the expression of ELF5 could be induced by cultivation of 786-0 cells under hyperosmotic conditions as shown by PCR or qPCR analyses (Figure 2A,B). Interestingly, the induction of ELF5 expression was higher in VHL+ cells than in VHL-deficient cells.
Our results clearly indicate that it is possible to induce the expression of the tumor suppressor ELF5 in RCC cells only by osmolality without any genetic manipulation. With the hyper-osmolality, we have identified a pathway that could be targeted for future intervention.
In the next step, we analyzed the predictive value of the hypertonicity-related genes for clinical outcome in ccRCC patients using the Cox proportional hazards model. Out of the 223 genes, 111 (49.8%) showed a significant effect (Table 1).
Within the genes that had a significant impact on patients’ survival, hypertonicity-downregulated genes tend to have a negative effect on survival (32 out of 41) while hypertonicity-upregulated genes have equal number of negative (35) and positive (35) effects. The corresponding data with the gene IDs and fold changes are provided in Table S1. This data suggests that the expression of hypertonicity-affected genes can be used to predict cancer-specific survival.
We next selected a minimum set of genes necessary for survival prediction using RNA-Seq data from the TCGA KIRC cohort. We identified 4 (COL1A1, NDUFA4L2, S100A6, MT2A) out of the 223 different genes that were regulated by hypertonicity in rats and subsequently defined our OSM score based on these four genes (Figure S1). Interestingly, all four genes have previously been associated with ccRCC tumorigenesis [9,10,11,12].
Our novel established OSM score based on these four genes was significantly associated with cancer-specific survival (Figure 3A; HR = 4.1; 95% CI: 2.78–6.07; p = 4.39 × 10−13; Cox proportional hazards regression model) in the TCGA cohort.
Multivariate analysis of the score together with clinicopathological parameters (T (primary tumor), N (regional lymph node), M (distant metastasis present at diagnosis), G (grading)) indicated that the score significantly predicted cancer-specific survival (p = 2.3 × 10−5, Table 2).
The independent role for prediction in the multivariate model was proven by analysis of deviance (p = 1.2 × 10−4). To validate these results, the OSM score was calculated in an independent cohort of ccRCC patients (n = 52; for details see [13]) based on expression levels of the selected genes and their model coefficients. Notably, we showed that the OSM score was also significantly associated with cancer-specific survival (Figure 3B; HR = 1.29; 95% CI: 1.05–1.59; p = 0.011). Multivariate analysis confirmed its role in the prediction of cancer-specific survival (p = 0.016) in our validation cohort (Table 1; analysis of deviance p = 0.0215).
A link between loss of VHL and osmolality has also been shown using kidney-specific VHL knock-out mice [5]. The authors observed that the mice had increased diuresis. The same group developed a renal cancer mice model and investigated the gene expression profile in mouse ccRCCs and kidney cortices [14]. Using these gene expression data, we could demonstrate that the mouse ccRCCs and normal kidney cortices could be discriminated based on the osmolality-regulated genes (Figure S2). Our results indicate that VHL function is important for hyper-osmolality-induced gene expression, as seen for ELF5. In a recent manuscript that we have submitted to Cancers we were able to show that the deletion of VHL also reduced the expression of several other hyper-osmolality-induced genes. This implies that VHL is prominently involved in the regulation hyper-osmolality-induced pathways. Since up to 85% of RCC patients harbor loss of VHL function, it is mandatory to identify the underlying cellular and molecular mechanisms.
In summary, our in vitro and in vivo data demonstrate that osmolality is an interesting pathway for the future development of drugs or other interventions in ccRCC which has not been considered so far. Moreover, this is the first report that defines an expression pattern of genes that can not only be used to discriminate between normal vs. tumor tissue and is associated with cancer-specific survival in independent ccRCC cohorts, but have a common physiological mechanism regulating their expression. Thus, targeting osmolality represents a novel interesting option for ccRCC therapy development, and further studies are warranted to identify the functional relevance of hypertonicity-associated pathways in tumor development and proliferation.

3. Materials and Methods

3.1. Primary Renal Cell Culture and RNA-Seq

Experiments were approved by a governmental committee on animal welfare (Landesamt für Natur, Umwelt und Verbraucherschutz Nordrhein-Westfalen, Germany) and were performed in accordance with national animal protection guidelines (A 60/1993 and A 67/09).
Primary cultured IMCD cells were prepared as described before [4]. For each group, three biological replicates were used. The groups included cells which had been cultivated at 300 or 900 mosmol/kg for one week. Total RNA was isolated using the mirVana miRNA Isolation Kit (Thermo Scientific, Waltham, MA, USA); 500 ng of total RNA were depleted of ribosomal RNA using the RiboMinus kit (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. Purified RNA was then fragmented by the addition of fragmentation buffer (200 mM Tris acetate, pH 8.2, 500 mM potassium acetate, and 150 mM magnesium acetate) and heating at 94 °C for 3 min in a thermocycler followed by ethanol precipitation with ammonium acetate and GlycoBlue (Thermo Fisher Scientific) as carrier. Fragmented RNA was then reverse transcribed using random hexamer and Superscript III (Thermo Fisher Scientific). The second strand was synthesized using the TargetAmp kit (Epicentre, Madison, WI, USA) according to the manufacturer’s instructions. The final steps of library preparation (e.g., blunt end repair, adapter ligation, adapter fill-in, and amplification) were done according to Meyer and Kircher [15]. The barcoded libraries were purified and quantified using the Library Quantification Kit (Illumina/Universal; KAPA Biosystems, Wilmington, MA, USA) according to the manufacturer’s instructions. A pool of up to 10 libraries was used for cluster generation at a concentration of 10 nM using an Illumina cBot. Sequencing of 2 × 100 bp was performed with an Illumina HiScanSQ sequencer at the sequencing core facility of the IZKF Leipzig (Faculty of Medicine, University Leipzig) using version 3 chemistry and flowcell according to the instructions of the manufacturer. Demultiplexing of raw reads, adapter trimming, and quality filtering were done according to Stokowy et al. [16] using TruSeq (Illumina) adapter sequences.

3.2. 786-0 Renal Cancer Cell Line and Real-Time PCR

The 786-0 and VHL-expressing 786-0-VHL were a kind gift of Prof. Barbara Seliger and were cultivated as described in [17]. For experimental setting, the cell culture medium was adjusted to 600 mosmol/kg by the addition of 100 mM NaCl and 100 mM urea. The cells were cultivated for different time points at 600 mosmol/kg. Total RNA isolation and cDNA synthesis were performed as described previously [4]. Real-time PCR was performed using the SYBR Green PCR Master Mix with the ABI PRISM 7900 Sequence Detection System. All instruments and reagents were purchased from Applied Biosystems (Darmstadt, Germany). Relative gene expression values were evaluated with the 2−ΔΔCt method using GAPDH as reference gene [18]. The primer sequences for ELF5 are ELF5-sense CGT GGA CTG ATC TGT TCA GCA ATG A, ELF5-antisense CAG GGT GGA CTG ATG TCC AGT ATG A and for GAPDH GAPDH-sense CAA GCT CAT TTC CTG GTA TGA C and GAPDH-antisense GTG TGG TGG GGG ACT GAG TGT GG.

3.3. Study Cohorts

Publicly available gene expression data from The Cancer Genome Atlas (TCGA) from a cohort of ccRCC patients (KIRC cohort, n = 449) were used to compare osmolality-induced genes expression between tumor and non-tumor samples. In this data set, 436 patients had both expression and CSS data and were used to develop the osmolality score. Expression data from tumor and non-tumor tissue were downloaded using the Bioconductor TCGAbiolinks package.
The validation cohort consisted of primary tumors with ccRCC histology (n = 52) of patients treated at the Department of Urology, University Hospital Tuebingen, Germany. Details of the study and tissue sample collection were described previously [13]. Transcriptome analyses was performed using the Human Transcriptome Array HTA 2.0 (Affymetrix/Thermo Fisher Scientific, Waltham, MA, USA), as described previously by Büttner et al. [13]. The accession number for genome-wide data at the European Genome-phenome Archive (EGA) (www.ebi.ac.uk/ega/home), which is hosted by the EBI and the CRG, is EGAS00001001176. Cancer-specific survival was used as endpoint in survival analyses of the development cohort (ccRCC KIRC) and the validation cohort, as described previously [13].

3.4. Statistical Analyses

RNA-Seq reads were aligned using bowtie2 and tophat2 to the rat reference genome (rnor6) according to Kim et al. [19]. Rat mRNA-Seq read counts were normalized and tested for differential expression using the Bioconductor edgeR package (v 3.24.3, [20]); 355 rat transcripts showed significant difference between two conditions: 900 and 300 mosmol/kg (Benjamini–Hochberg [21] adjusted p-Value < 0.05; log2FC >3/<−3). Of these, 284 of them had matching human transcripts by gene symbol. For those 284 transcripts (223 genes), mRNA-Seq expression values (FPKM-UQ) of the TCGA-KIRC cohort (449 tumor samples with 67 matching tissue normal samples) were clustered (hierarchical clustering with agglomeration method ward. D2 and Euclidean distance). Clustering proved that selected transcripts expression clearly discriminated between tumor and normal samples. The potential impact of those genes’ expression on TCGA-KIRC patients’ (n = 436) cancer-specific survival (CSS) was tested by building a Cox proportional hazards model on each gene’s expression separately (survival R package v 3.1-7, [21]). We found that 111 genes showed significant effect (Benjamini–Hochberg adjusted p-Value < 0.05). Later, we built a Cox proportional hazards model with lasso penalty based on the expression of the entire set of 223 genes (glmnet R package v 2.0-16, [22]). Four genes had non-zero coefficients according to the model, with minimal cross-validation error. Each TCGA-KIRC patient was assigned a survival score (termed the OSM score) calculated as the weighted sum of the expression of the four selected genes multiplied by the respective model coefficient. Analogously, the score was calculated for our independent RCC cohort (n = 52; [13]) based on selected genes expression (determined by microarray analyses) and the respective model coefficients. The value was multiplied by 100,000 times to avoid infinite hazards ratio. Patient cohorts were recursively partitioned based on the survival score using conditional inference trees [23] with the endpoint CSS. Multivariate survival analysis was performed using Cox proportional hazards regression models. Comparison of Cox models (with and without OSM score) was done using analysis of deviance [24].

4. Conclusions

Our study demonstrates that osmolality is an interesting pathway in ccRCC which has not yet been considered. The expression of hypertonicity-regulated genes is clearly associated with cancer-specific survival in ccRCC. We were also able to induce the expression of potentially tumor-suppressive genes by cultivating ccRCC cell lines under hyper-osmotic conditions.
Thus, targeting osmolality-associated pathways might represent a novel interesting therapeutic option for ccRCC.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/1/6/s1, Figure S1: Selection of differentially expressed transcripts affected by hypertonicity for development of the novel OSM-score, Figure S2: Hierarchical clustering of samples from the normal kidney cortex and mouse ccRCCs based on the hypertonicity affected genes, Table S1: Differentially expressed transcripts affected by hypertonicity and expression differences of corresponding human transcripts in ccRCC tumor and non-tumor tissue of the TCGA KIRC cohort.

Author Contributions

S.K. performed research, analyzed data, and wrote the paper; K.K. analyzed data and wrote the paper; S.W. analyzed data and wrote the paper; J.L. performed research; J.B. contributed samples and wrote the paper; M.S. contributed samples and wrote the paper; E.S. designed research, analyzed data, and wrote the paper; B.E. designed research, performed research, analyzed data, and wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) ED 181/9-1, the Robert Bosch Stiftung (Stuttgart, Germany), the ICEPHA Graduate School Tuebingen-Stuttgart, The German Cancer Consortium (DKTK, Germany), and the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—EXC 2180-390900677.

Acknowledgments

We would like to thank Ian Frew, Clinic of Internal Medicine I, University Clinic Freiburg for his helpful discussion and interpretation of the data. The results shown here are partly based on data generated by the TCGA Research Network. We would like to thank The Cancer Genome Atlas initiative, all tissue donors, and investigators who contributed to the acquisition and analyses of the samples used in this study. Information about TCGA and the investigators and institutions who constitute the TCGA research network can be found at http://cancergenome.nih.gov/.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lindgren, D.; Eriksson, P.; Krawczyk, K.; Nilsson, H.; Hansson, J.; Veerla, S.; Sjölund, J.; Höglund, M.; Johansson, M.E.; Axelson, H. Cell-Type-Specific Gene Programs of the Normal Human Nephron Define Kidney Cancer Subtypes. Cell Rep. 2017, 20, 1476–1489. [Google Scholar] [CrossRef] [Green Version]
  2. Hsieh, J.J.; Purdue, M.P.; Signoretti, S.; Swanton, C.; Albiges, L.; Schmidinger, M.; Heng, D.Y.; Larkin, J.; Ficarra, V. Renal cell carcinoma. Nat. Rev. Dis. Primers 2017, 3, 17009. [Google Scholar] [CrossRef]
  3. Büttner, F.; Winter, S.; Rausch, S.; Reustle, A.; Kruck, S.; Junker, K.; Stenzl, A.; Agaimy, A.; Hartmann, A.; Bedke, J.; et al. Survival Prediction of Clear Cell Renal Cell Carcinoma Based on Gene Expression Similarity to the Proximal Tubule of the Nephron. Eur. Urol. 2015, 68, 1016–1020. [Google Scholar] [CrossRef]
  4. Schulze Blasum, B.; Schroter, R.; Neugebauer, U.; Hofschroer, V.; Pavenstadt, H.; Ciarimboli, G.; Schlatter, E.; Edemir, B. The kidney-specific expression of genes can be modulated by the extracellular osmolality. FASEB J. 2016, 30, 3588–3597. [Google Scholar] [CrossRef] [Green Version]
  5. Schönenberger, D.; Rajski, M.; Harlander, S.; Frew, I.J. Vhl deletion in renal epithelia causes HIF-1α-dependent, HIF-2α-independent angiogenesis and constitutive diuresis. Oncotarget 2016, 7, 60971–60985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Jeon, U.S.; Kim, J.A.; Sheen, M.R.; Kwon, H.M. How tonicity regulates genes: Story of TonEBP transcriptional activator. Acta Physiol. 2006, 187, 241–247. [Google Scholar] [CrossRef] [PubMed]
  7. Bogusławska, J.; Popławski, P.; Alseekh, S.; Koblowska, M.; Iwanicka-Nowicka, R.; Rybicka, B.; Kędzierska, H.; Głuchowska, K.; Hanusek, K.; Tański, Z.; et al. MicroRNA-Mediated Metabolic Reprograming in Renal Cancer. Cancers 2019, 11, 1825, confirmed. [Google Scholar]
  8. Piggin, C.L.; Roden, D.L.; Gallego-Ortega, D.; Lee, H.J.; Oakes, S.R.; Ormandy, C.J. ELF5 isoform expression is tissue-specific and significantly altered in cancer. Breast Cancer Res. 2016, 18, 4. [Google Scholar] [CrossRef] [Green Version]
  9. Ibanez de Caceres, I.; Dulaimi, E.; Hoffman, A.M.; Al-Saleem, T.; Uzzo, R.G.; Cairns, P. Identification of Novel Target Genes by an Epigenetic Reactivation Screen of Renal Cancer. Cancer Res. 2006, 66, 5021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Pal, D.; Sharma, U.; Singh, S.K.; Mandal, A.K.; Prasad, R. Metallothionein gene expression in renal cell carcinoma. Indian J. Urol. IJU J. Urol. Soc. India 2014, 30, 241–244. [Google Scholar]
  11. Wang, L.; Peng, Z.; Wang, K.; Qi, Y.; Yang, Y.; Zhang, Y.; An, X.; Luo, S.; Zheng, J. NDUFA4L2 is associated with clear cell renal cell carcinoma malignancy and is regulated by ELK1. PeerJ 2017, 5, e4065. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Lyu, X.-J.; Li, H.-Z.; Ma, X.; Li, X.-T.; Gao, Y.; Ni, D.; Shen, D.-L.; Gu, L.-Y.; Wang, B.-J.; Zhang, Y.; et al. Elevated S100A6 (Calcyclin) enhances tumorigenesis and suppresses CXCL14-induced apoptosis in clear cell renal cell carcinoma. Oncotarget 2015, 6, 6656–6669. [Google Scholar] [CrossRef] [PubMed]
  13. Büttner, F.; Winter, S.; Rausch, S.; Hennenlotter, J.; Kruck, S.; Stenzl, A.; Scharpf, M.; Fend, F.; Agaimy, A.; Hartmann, A.; et al. Clinical utility of the S3-score for molecular prediction of outcome in non-metastatic and metastatic clear cell renal cell carcinoma. BMC Med. 2018, 16, 108. [Google Scholar] [CrossRef]
  14. Harlander, S.; Schönenberger, D.; Toussaint, N.C.; Prummer, M.; Catalano, A.; Brandt, L.; Moch, H.; Wild, P.J.; Frew, I.J. Combined mutation in Vhl, Trp53 and Rb1 causes clear cell renal cell carcinoma in mice. Nat. Med. 2017, 23, 869–877. [Google Scholar] [CrossRef] [Green Version]
  15. Meyer, M.; Kircher, M. Illumina Sequencing Library Preparation for Highly Multiplexed Target Capture and Sequencing. Cold Spring Harb. Protoc. 2010, 2010, pdb.prot5448. [Google Scholar] [CrossRef]
  16. Stokowy, T.; Eszlinger, M.; Świerniak, M.; Fujarewicz, K.; Jarząb, B.; Paschke, R.; Krohn, K. Analysis options for high-throughput sequencing in miRNA expression profiling. BMC Res. Notes 2014, 7, 144. [Google Scholar] [CrossRef] [Green Version]
  17. Leisz, S.; Schulz, K.; Erb, S.; Oefner, P.; Dettmer, K.; Mougiakakos, D.; Wang, E.; Marincola, F.M.; Stehle, F.; Seliger, B. Distinct von Hippel-Lindau gene and hypoxia-regulated alterations in gene and protein expression patterns of renal cell carcinoma and their effects on metabolism. Oncotarget 2015, 6, 11395–11406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−∆∆Ct Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  19. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [Green Version]
  20. McCarthy, D.J.; Smyth, G.K.; Robinson, M.D. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2009, 26, 139–140. [Google Scholar]
  21. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 1995, 85, 289–300. [Google Scholar] [CrossRef]
  22. Friedman, J.; Hastie, T.; Tibshirani, R. Regularization Paths for Generalized Linear Models via Coordinate Descent. J. Stat. Softw. 2010, 33, 1–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hothorn, T.; Zeileis, A. Partykit: A Modular Toolkit for Recursive Partytioning in R. J. Mach. Learn. Res. 2015, 16, 3905–3909. [Google Scholar]
  24. Therneau, T.M.; Grambsch, P.M. Modeling Survival Data: Extending the Cox Model; Springer: New York, NY, USA, 2000. [Google Scholar]
Figure 1. (A) Differentially expressed transcripts affected by hypertonicity in primary cultured inner medullary collecting duct (IMCD) cells either cultivated at 300 or 900 mosmol/kg. In total, 355 transcripts were differentially expressed (with a cut off log2 fold change of >3 and <−3) and there were matched human transcripts for 284 of those. Of those, 110 transcripts were downregulated and 174 transcripts were upregulated by hypertonicity. (B) Hierarchical clustering of samples from The Cancer Genome Atlas (TCGA) Kidney Clear Cell Carcinoma (KIRC) cohort based on the hypertonicity-affected genes. The expression levels of the top 223 regulated genes were extracted from the TCGA KIRC cohort, and hierarchical clustering was performed. This set of genes was able to clearly separate clear cell renal cell carcinoma (ccRCC) samples (dark green) from the normal tissue samples (light green).
Figure 1. (A) Differentially expressed transcripts affected by hypertonicity in primary cultured inner medullary collecting duct (IMCD) cells either cultivated at 300 or 900 mosmol/kg. In total, 355 transcripts were differentially expressed (with a cut off log2 fold change of >3 and <−3) and there were matched human transcripts for 284 of those. Of those, 110 transcripts were downregulated and 174 transcripts were upregulated by hypertonicity. (B) Hierarchical clustering of samples from The Cancer Genome Atlas (TCGA) Kidney Clear Cell Carcinoma (KIRC) cohort based on the hypertonicity-affected genes. The expression levels of the top 223 regulated genes were extracted from the TCGA KIRC cohort, and hierarchical clustering was performed. This set of genes was able to clearly separate clear cell renal cell carcinoma (ccRCC) samples (dark green) from the normal tissue samples (light green).
Cancers 12 00006 g001
Figure 2. E74-like ETS transcription factor 5 (ELF5) expression in 786-0 cells is induced by environmental hypertonicity. von Hippel–Lindau (VHL)-deficient 786-0 and VHL-expressing (VHL+) 786-0 cells were cultivated either in normal medium (300) or for 24, 48, and 72 h in 600 mosmol/kg medium. (A) The expression of ELF5 was analyzed by PCR. The osmolality was increased by the addition of 100 mM NaCl and 100 mM urea. The expression of ELF5 was VHL dependent. (B) The expression of ELF5 was quantified by real-time PCR in VHL-deficient 786-0 and VHL-expressing 786-0 VHL+ cells (N > 3, p-Value < 0.05 compared to control (CTL) using one-way ANOVA are marked by *, p-Value < 0.01 are marked by **). For more details about the PCR product of ELF5, please view the Supplementary Materials.
Figure 2. E74-like ETS transcription factor 5 (ELF5) expression in 786-0 cells is induced by environmental hypertonicity. von Hippel–Lindau (VHL)-deficient 786-0 and VHL-expressing (VHL+) 786-0 cells were cultivated either in normal medium (300) or for 24, 48, and 72 h in 600 mosmol/kg medium. (A) The expression of ELF5 was analyzed by PCR. The osmolality was increased by the addition of 100 mM NaCl and 100 mM urea. The expression of ELF5 was VHL dependent. (B) The expression of ELF5 was quantified by real-time PCR in VHL-deficient 786-0 and VHL-expressing 786-0 VHL+ cells (N > 3, p-Value < 0.05 compared to control (CTL) using one-way ANOVA are marked by *, p-Value < 0.01 are marked by **). For more details about the PCR product of ELF5, please view the Supplementary Materials.
Cancers 12 00006 g002
Figure 3. The OSM score could predict patient survival. (A) Kaplan–Meier plot indicating that hypertonicity-affected genes using the four selected genes (OSM score) can be applied for the prediction of patients’ cancer-specific survival in the TCGA KIRC cohort. (B) Kaplan–Meier plot indicating that the OSM score can be applied for the prediction of patients’ cancer-specific survival in the ccRCC validation cohort.
Figure 3. The OSM score could predict patient survival. (A) Kaplan–Meier plot indicating that hypertonicity-affected genes using the four selected genes (OSM score) can be applied for the prediction of patients’ cancer-specific survival in the TCGA KIRC cohort. (B) Kaplan–Meier plot indicating that the OSM score can be applied for the prediction of patients’ cancer-specific survival in the ccRCC validation cohort.
Cancers 12 00006 g003
Table 1. Number of hypertonicity-affected genes and their impact on patient survival.
Table 1. Number of hypertonicity-affected genes and their impact on patient survival.
Effect on Cancer Specific SurvivalDownregulatedUpregulated
hazardous3235
indifferent4666
favorable935
Table 2. Multivariate Cox regression for cancer-specific survival in the TCGA cohort (n = 409) and the validation cohort (n = 51).
Table 2. Multivariate Cox regression for cancer-specific survival in the TCGA cohort (n = 409) and the validation cohort (n = 51).
Multivariate AnalysesVariableLevelp-Value
(Wald Test)
HR (95% CI)
Including T, N, M, G and OSM scoreOSM score 2.35 × 10−52.6 (1.67–3.69)
(TCGA cohort)Primary tumorT1/T2 1
T3/T45.05 × 10−21.7 (1–2.86)
Lymph nodesN0 1
N16.15 × 10−22.71 (0.95–7.72)
NX2.4 × 10−20.6 (0.38–0.93)
Distant metastasisM0 1
M11.55 × 10−125.62 (3.48–9.07)
GradeG1/G2 1
G3/G46.69 × 10−32.35 (1.27–4.37)
Including T, N, M, G and OSM scoreOSM score 1.58 × 10−21.34 (1.06–1.7)
(Validation cohort)Primary tumorT1/T2 1
T3/T42.54 × 10−12.63 (0.5–13.9)
Lymph nodesN0 1
N1/N21.27 × 10−10.33 (0.08–1.36)
Distant metastasisM0 1
M12.78 × 10−541 (7.22–233.06)
Fuhrman gradeG1/G2 1
G3/G49.94 × 10−11 (0.29–3.43)
Abbreviations: CI, confidence interval; HR, hazard ratio; Ref., reference level; T, primary tumor; N, regional lymph node; M, distant metastasis present at diagnosis; G, grading. Cases with grading information “GX” or metastasis status “MX” were excluded from multivariate analyses. OSM scores were determined based on gene expression data measured by RNA-Seq (TCGA) or microarray analyses (validation cohort).

Share and Cite

MDPI and ACS Style

Kandabarau, S.; Leiz, J.; Krohn, K.; Winter, S.; Bedke, J.; Schwab, M.; Schaeffeler, E.; Edemir, B. Hypertonicity-Affected Genes Are Differentially Expressed in Clear Cell Renal Cell Carcinoma and Correlate with Cancer-Specific Survival. Cancers 2020, 12, 6. https://doi.org/10.3390/cancers12010006

AMA Style

Kandabarau S, Leiz J, Krohn K, Winter S, Bedke J, Schwab M, Schaeffeler E, Edemir B. Hypertonicity-Affected Genes Are Differentially Expressed in Clear Cell Renal Cell Carcinoma and Correlate with Cancer-Specific Survival. Cancers. 2020; 12(1):6. https://doi.org/10.3390/cancers12010006

Chicago/Turabian Style

Kandabarau, Siarhei, Janna Leiz, Knut Krohn, Stefan Winter, Jens Bedke, Matthias Schwab, Elke Schaeffeler, and Bayram Edemir. 2020. "Hypertonicity-Affected Genes Are Differentially Expressed in Clear Cell Renal Cell Carcinoma and Correlate with Cancer-Specific Survival" Cancers 12, no. 1: 6. https://doi.org/10.3390/cancers12010006

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