A study exploring critical pathways in clear cell renal cell carcinoma
- Authors:
- Published online on: November 7, 2013 https://doi.org/10.3892/etm.2013.1392
- Pages: 121-130
Abstract
Introduction
Renal cell carcinoma (RCC) is the most lethal type of cancer in the urinary system. It is a heterogeneous group of malignancies, and clear cell, papillary and chromophobe RCC are the major subtypes (1). In addition, clear cell RCC (CCRCC) is the most common type of RCC, constituting >80% of cases, and accounts for 2–3% of all adult malignant neoplasms. There were 58,240 new diagnoses of RCC or cancer of the renal pelvis and 13,040 mortalities due to the disease in 2010 in the United States. Furthermore, RCC accounted for 3% of all malignant tumors in adults and 85% of all primary malignant kidney tumors (2). RCC often presents as a metastatic disease and there are no effective treatments for it. This cancer exhibits resistance to chemotherapy and radiation and <10% of patients suffering from the metastatic disease survive five years subsequent to diagnosis. At present, there are no clinically validated adjuvant therapies for patients with CCRCC who are at a high risk of relapse following surgical resection (3). Overall, 40% of patients with RCC succumb to this disease (4). Certain molecular mechanisms have been revealed in RCC, such as von Hipple-Lindau (VHL) gene mutation. This appears to be responsible for ~60% of the cases of sporadic CCRCC (5). In hereditary papillary renal carcinoma, an oncogene, mesenchymal-epithelial transition factor (MET), has been revealed to be mutated, although the incidence of c-MET mutations is low in sporadic papillary RCC. However, to date, no theory has been able to explain all the cases of CCRCC, and there are no effective tumor markers that may be used for diagnostic, prognostic and treatment purposes in CCRCC. To develop novel therapeutic interventions for CCRCC an improved understanding of the pathogenesis is required, since traditional single-factor analysis for tumor pathogenesis is not effective. In order to obtain an enhanced understanding of the potential ‘cross-talks’ between dysregulated genes and proteins, there has been particular focus on undertaking a more ‘global’ analysis (6). Gene set enrichment analysis (GSEA), based on gene expression profiling, has been performed in CCRCC and has identified certain dysregulated gene sets in patients with G1 and G3 RCC (7).
At present, gene expression profiling between tumor tissue and normal tissue at the RNA level has been used in a previous study of RCC (6). A number of differentially expressed genes have been identified in these studies; however, few studies have been performed based on pathway analysis. While it is not difficult to obtain the gene expression profiles, the extraction of biological insight from such information remains a major challenge (8). Genome analysis provides a large quantity of information concerning molecules that are involved in disease pathogenesis and may be used to interpret diseases. Pathway analysis is a frequently used analytical tool. GSEA, which was introduced by Mootha et al(9), is well known and extensively used in gene set analysis. GSEA is a computational method that identifies pre-defined gene sets that exhibit significant differences in expression between samples from normal individuals and patients. Following its introduction, the methodology of GSEA was subsequently improved by Subramanian et al(8). This modified approach aimed to identify the significant differences in the expression of groups or pathways rather than individual genes. In view of the slight nature of the functions of individual genes, GSEA is more efficient than individual gene analysis for the study of complex diseases (10). Based on the study of a group of associated genes, investigators are likely to have an improved chance of determining the pivotal functional processes responsible for biological phenomena. GSEA in a single dataset is unlikely to identify the significant pathways and there is likely to be a poor concordance. In order to obtain improved data convergence and provide a systematic insight into the pathways altered during CCRCC pathogenesis, standardized microarray preprocessing, GSEA and DAVID functional enrichment analysis were performed in the present study.
Materials and methods
Selection and obtainment of datasets
The Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/) and ArrayExpress (http://www.ebi.ac.uk/arrayexpress/) were searched for gene expression profiling datasets associated with CCRCC, using the keywords ‘renal cell carcinoma’ up to March 25, 2013. The inclusion criteria for the datasets were as follows: i) The dataset was genome-wide; ii) samples included patients with CCRCC and normal controls; iii) complete microarray raw or normalized datasets were available.
In addition, seven public gene expression datasets that concurred with our inclusion criteria were included in our study to assess the transcripts of CCRCC on a genome-wide basis. In datasets GSE15641 and GSE11151, there were other types of pathological samples besides CCRCC and normal tissue; in these instances, only the CCRCC and normal tissue samples were retained. In dataset GSE16449, samples from GSM413237 to GSM413270 had the same probe number of 33998; however, the other 35 samples had a different identical probe number and different data precessing, and were excluded. In dataset GSE12606, the samples comprised three CCRCC tissues, including autologous normal tissue and autologous metastatic tissue; the autologous metastatic samples were excluded. Two platforms, HG-U133A and HG-U133B, were applied in datasets GSE781 and GSE6344, however, only the expression profiles in HG-U133A were included in the present study. The basic information with regard to these datasets, such as the microarray platform, experimental design, probe numbers and sample size, are listed in Table I.
Standardized microarray preprocessing of the datasets
The datasets were processed with an improved software package based on version 2.4.0 of Bioconductor (18) and R version 2.10.1 (19). Dataset GSE16449 comprised normalized data and was retained. Each Affymetrix dataset was computed using the Robust Multichip Averaging (RMA) algorithm in the Affy package (20,21). Genes were excluded from further analysis if it was not possible to map them to any Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway (http://www.genome.jp/kegg/pathway.html). Interquartile ranges (IQRs) were used to describe the level of variability and the genes with an IQR value of <0.5 were excluded. The probe-set with the largest variability was retained when several probe-sets hit one gene.
GSEA for pathways
GSEA was performed separately in each dataset in the Category package (version 2.10.1) (19). Pathways were excluded if <10 genes were identified to be dysregulated. The t-statistic mean of the genes was computed in each pathway. Using a permutation test with 1,000 times, the significantly changed pathways were identified with a P-value of ≤0.05.
Meta-analysis and DAVID functional enrichment analysis
A Student’s t-test comparing the CCRCC and normal samples for every probe was performed independently in all of the six Affymetrix datasets using Excel 2007, and probes with a P-value <0.05 were considered to be significantly different. In all included datasets, with the exception of dataset GSE16449 with the Agilent platform, other datasets were checked with the Affymetrix platform. In order to keep the data uniform, we only do meta-analysis on the datasets with Affymetrix platform. Following the t-test, the χ2 square value of each selected probe was calculated using the formula: χ2=∑(−2 logepi) (22). The probes with χ2 values <0.05 were retained and used to obtain KEGG pathways from DAVID Bioinformatics Resources 6.7 (http://david.abcc.ncifcrf.gov/).
Results
Re-analysis of each dataset to obtain the dysregulated pathways
The seven datasets included 110 CCRCC and 73 normal samples. The common method of GSEA was implemented for all of the seven datasets, respectively. In dataset GSE12606, GSEA did not identify any dysregulated pathways. This result may have been due to the small sample size of GSE12606 (three cases and three controls). For the remaining six datasets, significantly abnormal pathways were identified using GSEA, based on the permutation P-values. The results are described in Table II.
Common pathways identified by GSEA
The results of the GSEA were compared and it was revealed that 17 down- and 12 upregulated pathways coexisted in six datasets. The majority of the downregulated pathways were associated with carbohydrate, lipid and amino acid metabolism, while others were involved in organismal systems or cellular processes. The names of these pathways are listed in Table III. A comparison of the dysregulated pathways, including the up- and downregulated pathways, revealed that certain pathways were upregulated in one or more datasets while being downregulated in another or others.
Dysregulated pathways identified by DAVID functional enrichment analysis
The χ2 value for 3,149 probes <0.05 was calculated using the following formula: χ2=∑(−2 logepi). The IDs of these probes were all used in the DAVID enrichment analysis. A total of 1,023 probes (genes) led to further pathway analysis and 32 dysregulated pathways were identified. The abnormal pathways and the genes involved are listed in Table IV.
![]() | Table IVDysregulated pathways identified by DAVID functional pathway analysis based on the abnormal gene. |
Discussion
CCRCC is a frequently occurring type of cancer that presents numerous challenges. Although a number of preceding studies have attempted to identify the pathogenesis and mechanisms behind this type of cancer, none of the existing theories are able to explain all the cases of CCRCC. Genome-wide expression chips are powerful tools that enable the comprehensive identification of abnormal gene families or pathways present in relevant disease states. Thus, biologically relevant inferences may be reproducible in different studies. For single-gene analysis, examining the same biological condition using different statistical methods and datasets may lead to significant discrepancies (23). Compared with gene expression obtained from different datasets, pathway analysis applied to different datasets yields consistent results and diminishes large discrepancies. Therefore, pathway analysis is able to highlight genes weakly connected to the phenotype that may be difficult to detect in classical univariate statistics.
In our study, GSEA was performed using seven independent publicly available gene expression datasets and the six of these datasets that utilized the Affymetrix platform underwent a meta-analysis and pathway analysis using DAVID functional enrichment analysis. A cross-study based on GSEA was performed to identify critical pathways and to obtain a deeper understand of the common mechanisms involved in CRCC. The results suggested that the majority of the dysregulated pathways were consistent in different studies. Many of the dysregulated pathways identified in the present study have already been indicated to be involved in CCRCC or other types of cancer. In the following paragraphs, a number of the dysregulated pathways and hypotheses regarding the importance of these pathways in CCRCC, based on the functional classification, are discussed.
As shown in Table III, a number of common downregulated pathways were identified by GSEA, including the citrate cycle (TCA cycle) and its associated pathways. Previous studies have demonstrated that the metabolism of a variety of nutrients is dysregulated in certain types of RCC. VHL, MET, folliculin (FLCN), tuberous sclerosis 1 (TSC1), TSC2, fumarate hydratase (FH) and succinate dehydrogenase (SDH) are known as renal cancer genes, and all are involved in pathways that respond to metabolic stress or nutrient stimulation. Individuals who harbor mutations in any of these genes have an increased risk of developing RCC. It has been suggested that RCC may be regarded as a metabolic disease (24). When the current study focused on CCRCC, the TCA cycle pathway and a number of the pathways associated with it were identified as being downregulated. FH and SDH genes (SDHA, SDHC and SDHD) were also identified as having a low expression levels. The FH and SDH genes encode mitochondrial tricarboxylic acid cycle enzymes that are essential in energy metabolism. With deficient FH and SDH expression, the process of the TCA cycle is inhibited and oxidative phosphorylation is hindered. Previous studies have shown that certain FH-deficient kidney tumor cell lines exhibit significantly impaired oxidative phosphorylation (25,26) and are glucose-dependent (25–27). It has been demonstrated that the glucose-mediated generation of cellular reactive oxygen species in an FH-deficient kidney cancer cell line results in the stabilization of hypoxia-inducible factor (HiF) 1-α (27).
Similar to FH, mutations in three of the four SDH genes (SDHB, SDHC and SDHD) have been shown to be involved in familial paraganglioma and familial pheochromocytoma (28,29). Deficient expression of SDH genes severely impairs oxidative phosphorylation and results in SDH-deficient kidney tumors showing glucose dependence. SDH-deficient kidney tumors are likely to be as sensitive to glucose as FH-deficient tumors and targeting the vasculature or glucose transport in these tumors may provide a novel approach to disrupt the fundamental metabolic machinery of these aggressive cancers. Mutations of FH and SDH cause the inactivation of prolyl hydroxylase (PHD), resulting in the stabilization of HiFs and driving the transcriptional activation of genes supporting tumor growth, neovascularization, invasion and metastasis. In the present study, FH and SDH were revealed to have a low expression level in CCRCC. Therefore, we propose that the downregulation of oxidative phosphorylation, the TCA cycle pathway and its associated pathways may have a critical role in the generation and development of CCRCC. The abnormal metabolism of the upstream nutrients of the TCA cycle may therefore be attributed to the blocked TCA cycle.
The tight junction pathway and its associated genes are also indicated to be important in CCRCC. Tight junctions are composed of occludin, claudins and junctional adhesion molecules. Together with adherens junctions and desmosomes, they form the apical junctional complex in epithelial and endothelial cellular sheets. Tight junctions are essential for the tight sealing of the cellular sheets, thus controlling paracellular ion flux and maintaining tissue homeostasis (30). Due to the fact that tight junction proteins are able to recruit signaling proteins (31), these junctions have been assumed to be involved in the regulation of proliferation, differentiation and other cellular functions. There is a generally accepted view that tumorigenesis is accompanied by a disruption of tight junctions. Furthermore, the disruption of tight junctions may be important in the loss of cohesion, invasiveness and lack of differentiation observed in cancer cells. Consistent with this idea, several previous studies have revealed that the expression of occludin and claudins are altered in a number of types of cancer. Occludin was shown to be frequently downregulated in gastrointestinal tumors (32) and a decreased expression of claudin-1 was observed in breast (33,34) and colon cancer (35). The expression pattern of claudins is tissue-specific and it has been demonstrated that claudin-3 (transcript of CLDN3) and/or claudin-4 (transcript of CLDN4 ) are expressed in normal renal tissue (36). To the best of our knowledge, there have been no studies investigating the differential expression of CLDN3 and CLDN4 in normal renal and RCC tissues. When CCRCC tissues were compared with normal renal tissues using GSEA, the tight junction pathway was revealed to be downregulated in CCRCC tissues in different studies. Furthermore, CLDN3, CLDN4, CLDN10, junctional adhesion molecule 2 (JAM2) and JAM3, whose transcripts contribute to the formation of tight junctions, showed decreased expression in CCRCC tissue. Following the review of previous studies, Morin (37) considered that it was possible that the loss of claudins and other tight junction proteins in cancer was responsible for the loss of cell adhesion, an important step in the progression of cancer to metastasis (37). Therefore, the destruction of tight junctions may be important in the development of CCRCC.
A further factor that was considered to be important in the pathogenesis of CRCC was the downregulation of retinol metabolism and its associated genes. It is generally recognized that retinol and its metabolites have significant effects in cell differentiation and apoptosis. In the present study, GSEA revealed that retinol metabolism was downregulated when comparing CCRCC tissues with normal renal tissues, and a number of genes included in this pathway were observed to have a low level of expression, such as alcohol dehydrogenase 1A (ADH1A), ADH1B, ADH1C and cytochrome P450s (CYPs). The inhibition of retinol metabolism results in a lack of its metabolite and, in this case, the blockage of the differentiation and apoptosis of cells. A number of the metabolites are currently widely used in the treatment of several types of cancer. The most notable is the use of all-trans retinoic acid in the treatment of acute promyelocytic leukemia (38). Furthermore, several studies have demonstrated that one of the metabolites of retinol, 13-cis-retinoic acid, showed improved efficacy in the treatment of RCC when combined with other antitumor drugs (39–44). The results of the current study further indicated that the dysregulated metabolism of retinol may be important in CCRCC. In order to further reveal the role of downregulated retinol metabolism in CCRCC, additional studies focusing on these abnormal genes and the affected metabolites are required. The determination of whether retinol or its metabolites exert therapeutic effects in the treatment of CCRCC may be of benefit.
In addition to downregulated pathways, GSEA also identified a number of pathways that were commonly upregulated. Appearing as solid malignant neoplasms in the human body, RCC has been demonstrated to be immunogenic. A number of immune cells have been isolated from RCC, including natural killer (NK) cells, cytotoxic T cells with specificity for autologous tumor cells, helper T cells and dendritic cells that express interleukin (IL)-1 and IL-2 and function as antigen-presenting cells (45–48). When RCC appears as an antigen in the human body, immune activity is induced, leading to a series of enhanced immunization activities. Nine pathways, which were associated with the immune system or immune system diseases and which co-existed in six independent studies, were identified to be upregulated when we compared CCRCC tissue with normal renal tissue using GSEA. Antigen processing and presentation, cytokine-cytokine receptor interaction and its associated pathways, including the chemokine signaling pathway, and the Jak-STAT signaling pathway were revealed to be upregulated in the present study.
When focusing on the chemokine signaling pathway, it was revealed that many genes involved in this pathway were overexpressed. Consistent with a number of previous studies (49,50), CXCR4, a gene with a transcript known as fusin, and its ligand gene CXCL12, whose transcript is known as stromal cell-derived factor-1α (SDF-1α), were demonstrated to be overexpressed in CCRCC in our study. CXCR4 is considered to be the major chemokine receptor expressed on cancer cells (51–53) and its expression is required for tumor metastasis to other organs; CXCR4 activation by CXCL12 induces the migration of neoplastic cells (54). The strong expression of CXCR4 has been indicated to correlate with the metastasis of RCC (55). The expression of CXCR4 and SDF-1α is negatively regulated by VHL (49), i.e. a deficiency in VHL expression results in the strong expression of CXCR4 and CXCL12. It has been proposed that blocking signaling through this chemokine receptor may inhibit micrometastases in RCC and result in improved outcomes (56). Therefore, the in-depth study of the adjustment mechanisms of CXCR4 expression may be beneficial in the development of novel treatments for metastatic RCC.
In mammals, the Jak-STAT pathway is the principal signaling mechanism for numerous cytokines and growth factors. These signals induce proliferation or cell fate determination and are crucial to the proper growth and development of mammalian tissues. Too strong or too weak activity of this signaling pathway leads to severe consequences. One of these cytokines, interferon-α (IFN-α) is widely used in the therapy of several neoplasms, including RCC. IFN-α exerts its effects by activating the Jak-STAT signaling pathway and subsequently inducing target gene transcription (57). Shang et al(58) demonstrated that the defective ability of Jak-STAT in RCC may result in IFN-α resistance (58). In the present study, the Jak-STAT signaling pathway was shown to be upregulated in CCRCC tissue. This was inconsistent with the fact that a substantial proportion of patients with RCC fail to respond to IFN-α treatment. The mechanisms of IFN resistance in RCC have not yet been fully elucidated. Further study for an in-depth understanding of this mechanism may be beneficial for the development of novel treatments.
NK cells are one type of infiltrating immune cell present in RCC and the NK cell-mediated cytotoxicity pathway was shown to exhibit strong activity in CCRCC in the present study. Enhancing CCRCC cell sensitivity to the cytotoxicity of NK cells may be an efficient treatment for CCRCC. Several studies concerning new treatments for RCC based on the toxic effects of NK cells have been performed (59,60).
As malignant neoplasms appear in the human body, CCRCC tissue shows as a self-antigen and induces an autoimmune response. Human leukocyte antigens (HLAs) are key components of the major histocompatibility complex (MHC), which is involved in immune cell signaling processes. Therefore, the overexpression of HLA-B, HLA-C and a number of other HLA alleles, as well as the upregulation of certain immune system disease pathways in CCRCC, may be attributed to the self-antigenicity of CCRCC.
The present study also compared the common dysregulated pathways identified by GSEA with the dysregulated pathways obtained from DAVID enrichment analysis. The TCA cycle pathway, numerous nutrient metabolic pathways and the tight junction pathway were shown to be dysregulated in DAVID functional enrichment analysis, in addition to GSEA, based on the abnormal genes identified by meta-analysis. However,, a number of the dysregulated pathways identified using DAVID functional enrichment analysis were not consistent with the cross-GSEA. In this discussion, only the abnormal pathways co-existing in the six datasets have been regarded as significant and discussed. When we reviewed the design of every dataset included in our study, however, certain datasets had small control samples and it was difficult to identify the significant genes. This indicated that a dataset with a small sample size may lead to bias. We also analyzed the dysregulated pathways co-existing in five of the six datasets. There were 17 down- and 12 upregulated pathways co-existing in five of the CCRCC datasets. When compared with the dysregulated pathways in the meta-analysis, nitrogen metabolism, DNA replication, the peroxisome proliferator-activated receptor (PPAR) signaling pathway, nucleotide excision repair, cell cycle, cell adhesion molecules (CAMs), the Toll-like receptor signaling pathway and Vibrio cholerae infection were reduplicated. This result further indicated that the poor design of the individual dataset may lead to large bias in a cross-GSEA.
Numerous treatments have been developed based on the current understanding of RCC; however, none of the existing therapies are universally valid. The treatment of metastatic RCC remains a significant challenge. More effectual studies are required to reveal the pathogenesis and progress of RCC. The current study demonstrated that a cross-GSEA, based on a number of well-designed datasets, is likely to improve the consistency of the understanding of RCC. Additional cross-GSEAs may further the progress made in this study.
In conclusion, cross-GSEA was used to identify the critical pathways involved in RCC using gene expression datasets of CCRCC. A number of pathways that were associated with RCC or CCRCC were identified as being significant, and the understanding of RCC has been improved. The downregulation of the TCA cycle pathway, retinol metabolism and the tight junction pathway may be critical in the occurrence of CCRCC. In addition, the upregulation of the chemokine signaling pathway and its associated pathways may contribute to the metastasis of CCRCC; however, this is the focus of immunotherapy. A further cross-GSEA based on a number of well-designed datasets may provide a valuable contribution to the progress made in the current study, following the analysis of its results.
Acknowledgements
This study was supported by the National Natural Science Foundation of China (#81060213 and #81272853) and the GuangXi Nature Science Foundation (#2011GXNSFB018100 and #2012GXNSFAA053152).
References
Eble JN, Togashi K, Pisani P, et al: Renal cell carcinoma. Pathology and Genetics of Tumours of the Urinary System and Male Genital Organs. WHO Classification of Tumours. 7. Eble JN, Sauter G, Epstein JI and Sesterhenn IA: IARC Press; Lyon: pp. 12–32. 2004 | |
Jemal A, Siegel R, Xu J and Ward E: Cancer statistics, 2010. CA Cancer J Clin. 60:277–300. 2010. View Article : Google Scholar | |
Jacobsohn KM and Wood CG: Adjuvant therapy for renal cell carcinoma. Semin Oncol. 33:576–582. 2006. View Article : Google Scholar : PubMed/NCBI | |
Janzen NK, Kim HL, Figlin RA and Belldegrun AS: Surveillance after radical or partial nephrectomy for localized renal cell carcinoma and management of recurrent disease. Urol Clin North Am. 30:843–852. 2003. View Article : Google Scholar : PubMed/NCBI | |
Kim WY and Kaelin WG: Role of VHL gene mutation in human cancer. J Clin Oncol. 22:4991–5004. 2004. View Article : Google Scholar : PubMed/NCBI | |
Setlur SR, Royce TE, Sboner A, Mosquera JM, Demichelis F, Hofer MD, Mertz KD, Gerstein M and Rubin MA: Integrative microarray analysis of pathways dysregulated in metastatic prostate cancer. Cancer Res. 67:10296–10303. 2007. View Article : Google Scholar : PubMed/NCBI | |
Maruschke M, Reuter D, Koczan D, Hakenberg OW and Thiesen HJ: Gene expression analysis in clear cell renal cell carcinoma using gene set enrichment analysis for biostatistical management. BJU Int. 108:E29–E35. 2011. View Article : Google Scholar : PubMed/NCBI | |
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES and Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 102:15545–15550. 2005. View Article : Google Scholar : PubMed/NCBI | |
Mootha VK, Lindgren CM, Eriksson KF, et al: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 34:267–273. 2003. View Article : Google Scholar : PubMed/NCBI | |
Shi J and Walker MG: Gene set enrichment analysis (GSEA) for interpreting gene expression profiles. Curr Bioinform. 2:133–137. 2007. View Article : Google Scholar | |
Lenburg ME, Liou LS, Gerry NP, Frampton GM, Cohen HT and Christman MF: Previously unidentified changes in renal cell carcinoma gene expression identified by parametric analysis of microarray data. BMC Cancer. 3:312003. View Article : Google Scholar : PubMed/NCBI | |
Jones J, Otu H, Spentzos D, et al: Gene signatures of progression and metastasis in renal cell cancer. Clin Cancer Res. 11:5730–5739. 2005. View Article : Google Scholar : PubMed/NCBI | |
Wang Y, Roche O, Yan MS, et al: Regulation of endocytosis via the oxygen-sensing pathway. Nat Med. 15:319–324. 2009. View Article : Google Scholar : PubMed/NCBI | |
Yusenko MV, Zubakov D and Kovacs G: Gene expression profiling of chromophobe renal cell carcinomas and renal oncocytomas by Affymetrix GeneChip using pooled and individual tumours. Int J Biol Sci. 5:517–527. 2009. View Article : Google Scholar | |
Gumz ML, Zou H, Kreinest PA, et al: Secreted frizzled-related protein 1 loss contributes to tumor phenotype of clear cell renal cell carcinoma. Clin Cancer Res. 13:4740–4749. 2007. View Article : Google Scholar : PubMed/NCBI | |
Brannon AR, Reddy A, Seiler M, et al: Molecular stratification of clear cell renal cell carcinoma by consensus clustering reveals distinct subtypes and survival patterns. Genes Cancer. 1:152–163. 2010. View Article : Google Scholar | |
Stickel J, Weinzierl AO, Hillen N, Drews O, Schuler MM, Hennenlotter J, Wernet D, Müller CA, Stenzl A, Rammensee HG and Stevanović S: HLA ligand profiles of primary renal cell carcinoma maintained in metastases. Cancer Immunol Immunother. 58:1407–1417. 2009. View Article : Google Scholar : PubMed/NCBI | |
Gentleman RC, Carey VJ, Bates DM, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 5:R802004. View Article : Google Scholar : PubMed/NCBI | |
Gentleman R: Using Categories to Model Genomic Data. 2009, Available from: https://www.bioconductor.org. Accessed April 3, 2013 | |
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U and Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 4:249–264. 2003. View Article : Google Scholar | |
Gautier L, Cope L, Bolstad BM and Irizarry RA: affy - analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 20:307–315. 2004. View Article : Google Scholar : PubMed/NCBI | |
Brown BM: A method for combining non-independent, one-sided tests of significance. Biometrics. 31:987–992. 1975. View Article : Google Scholar | |
Manoli T, Gretz N, Gröne HJ, Kenzelmann M, Eils R and Brors B: Group testing for pathway analysis improves comparability of different microarray datasets. Bioinformatics. 22:2500–2506. 2006. View Article : Google Scholar : PubMed/NCBI | |
Linehan WM, Srinivasan R and Schmidt LS: The genetic basis of kidney cancer: a metabolic disease. Nat Rev Urol. 7:277–285. 2010. View Article : Google Scholar : PubMed/NCBI | |
Yang Y, Valera VA, Padilla-Nash HM, et al: UOK 262 cell line, fumarate hydratase deficient (FH−/FH−) hereditary leiomyomatosis renal cell carcinoma: in vitro and in vivo model of an aberrant energy metabolic pathway in human cancer. Cancer Genet Cytogenet. 196:45–55. 2010.PubMed/NCBI | |
Xie H, Valera VA, Merino MJ, Amato AM, Signoretti S, Linehan WM, Sukhatme VP and Seth P: LDH-A inhibition, a therapeutic strategy for treatment of hereditary leiomyomatosis and renal cell cancer. Mol Cancer Ther. 8:626–635. 2009. View Article : Google Scholar : PubMed/NCBI | |
Sudarshan S, Sourbier C, Kong HS, et al: Fumarate hydratase deficiency in renal cancer induces glycolytic addiction and hypoxia-inducible transcription factor-1alpha stabilization by glucose-dependent generation of reactive oxygen species. Mol Cell Biol. 29:4080–4090. 2009. View Article : Google Scholar | |
Baysal BE, Ferrell RE, Willet-Brozick JE, et al: Mutations in SDHD, a mitochondrial complex II gene, in hereditary paraganglioma. Science. 287:848–851. 2000. View Article : Google Scholar : PubMed/NCBI | |
Vanharanta S, Buchta M, McWhinney SR, et al: Early-onset renal cell carcinoma as a novel extraparaganglial component of SDHB-associated heritable paraganglioma. Am J Hum Genet. 74:153–159. 2004. View Article : Google Scholar : PubMed/NCBI | |
Tsukita S, Furuse M and Itoh M: Multifunctional strands in tight junctions. Nat Rev Mol Cell Biol. 2:285–293. 2001. View Article : Google Scholar : PubMed/NCBI | |
Mitic LL and Anderson JM: Molecular architecture of tight junctions. Annu Rev Physiol. 60:121–142. 1998. View Article : Google Scholar | |
Kimura Y, Shiozaki H, Hirao M, Maeno Y, Doki Y, Inoue M, Monden T, Ando-Akatsuka Y, Furuse M, Tsukita S and Monden M: Expression of occludin, tight-junction-associated protein, in human digestive tract. Am J Pathol. 151:45–54. 1997.PubMed/NCBI | |
Krämer F, White K, Kubbies M, Swisshelm K and Weber BH: Genomic organization of claudin-1 and its assessment in hereditary and sporadic breast cancer. Hum Genet. 107:249–256. 2000.PubMed/NCBI | |
Tokés AM, Kulka J, Paku S, Szik A, Páska C, Novák PK, Szilák L, Kiss A, Bögi K and Schaff Z: Claudin-1, −3 and −4 proteins and mRNA expression in benign and malignant breast lesions: a research study. Breast Cancer Res. 7:R296–R305. 2005. | |
Hewitt KJ, Agarwal R and Morin PJ: The claudin gene family: expression in normal and neoplastic tissues. BMC Camcer. 1862006. View Article : Google Scholar : PubMed/NCBI | |
Long H, Crean CD, Lee WH, Cummings OW and Gabig TG: Expression of Clostridium perfringens enterotoxin receptors claudin-3 and claudin-4 in prostate cancer epithelium. Cancer Res. 61:7878–7881. 2001. | |
Morin PJ: Claudin proteins in human cancer: promising new targets for diagnosis and therapy. Cancer Res. 65:9603–9606. 2005. View Article : Google Scholar : PubMed/NCBI | |
Asou N: 2. All-trans retinoic acid in the treatment of acute promyelocytic leukemia. Intern Med. 46:91–93. 2007. View Article : Google Scholar : PubMed/NCBI | |
Motzer RJ, Schwartz L, Law TM, Murphy BA, Hoffman AD, Albino AP, Vlamis V and Nanus DM: Interferon alfa-2a and 13-cis-retinoic acid in renal cell carcinoma: antitumor activity in a phase II trial and interactions in vitro. J Clin Oncol. 13:1950–1957. 1995.PubMed/NCBI | |
Recchia F, Saggio G, Amiconi G, Di Blasio A, Cesta A, Candeloro G, Necozione S, Fumagalli L and Rea S: Multicenter phase II study of chemo-immunotherapy in the treatment of metastatic renal cell carcinoma. J Immunother. 30:448–454. 2007. View Article : Google Scholar : PubMed/NCBI | |
Boorjian SA, Milowsky MI, Kaplan J, Albert M, Cobham MV, Coll DM, Mongan NP, Shelton G, Petrylak D, Gudas LJ and Nanus DM: Phase 1/2 clinical trial of interferon alpha2b and weekly liposome-encapsulated all-trans retinoic acid in patients with advanced renal cell carcinoma. J Immunother. 30:655–662. 2007. View Article : Google Scholar : PubMed/NCBI | |
Motzer RJ, Murphy BA, Bacik J, Schwartz LH, Nanus DM, Mariani T, Loehrer P, Wilding G, Fairclough DL, Cella D and Mazumdar M: Phase III trial of interferon alfa-2a with or without 13-cis-retinoic acid for patients with advanced renal cell carcinoma. J Clin Oncol. 18:2972–2980. 2000.PubMed/NCBI | |
Aass N, De Mulder PH, Mickisch GH, Mulders P, van Oosterom AT, van Poppel H, Fossa SD, de Prijck L and Sylvester RJ: Randomized phase II/III trial of interferon Alfa-2a with and without 13-cis-retinoic acid in patients with progressive metastatic renal cell carcinoma: the European Organisation for Research and Treatment of Cancer Genito-Urinary Tract Cancer Group (EORTC 30951). J Clin Oncol. 23:4172–4178. 2005. | |
Atzpodien J, Kirchner H, Rebmann U, et al: Interleukin-2/interferon-alpha2a/13-retinoic acid-based chemoimmunotherapy in advanced renal cell carcinoma: results of a prospectively randomised trial of the German Cooperative Renal Carcinoma Chemoimmunotherapy Group (DGCIN). Br J Cancer. 95:463–469. 2006. View Article : Google Scholar | |
Elsässer-Beile U, Grussenmeyer T, Gierschner D, Schmoll B, Schultze-Seemann W, Wetterauer U and Schulte Mönting J: Semiquantitative analysis of Th1 and Th2 cytokine expression in CD3+, CD4+, and CD8+ renal-cell-carcinoma-infiltrating lymphocytes. Cancer Immunol Immunother. 48:204–208. 1999.PubMed/NCBI | |
Finke JH, Rayman P, Edinger M, Tubbs RR, Stanley J, Klein E and Bukowski R: Characterization of a human renal cell carcinoma specific cytotoxic CD8+ T cell line. J Immunother (1991). 11:1–11. 1992. View Article : Google Scholar : PubMed/NCBI | |
Gaudin C, Dietrich PY, Robache S, Guillard M, Escudier B, Lacombe MJ, Kumar A, Triebel F and Caignard A: In vivo local expansion of clonal T cell subpopulations in renal cell carcinoma. Cancer Res. 55:685–690. 1995.PubMed/NCBI | |
Schwaab T, Schned AR, Heaney JA, Cole BF, Atzpodien J, Wittke F and Ernstoff MS: In vivo description of dendritic cells in human renal cell carcinoma. J Urol. 162:567–573. 1999. View Article : Google Scholar : PubMed/NCBI | |
Zagzag D, Krishnamachary B, Yee H, Okuyama H, Chiriboga L, Ali MA, Melamed J and Semenza GL: Stromal cell-derived factor-1alpha and CXCR4 expression in hemangioblastoma and clear cell-renal cell carcinoma: von Hippel-Lindau loss-of-function induces expression of a ligand and its receptor. Cancer Res. 65:6178–6188. 2005. View Article : Google Scholar : PubMed/NCBI | |
Schrader AJ, Lechner O, Templin M, et al: CXCR4/CXCL12 expression and signalling in kidney cancer. Br J Cancer. 86:1250–1256. 2002. View Article : Google Scholar : PubMed/NCBI | |
Balkwill F: The significance of cancer cell expression of the chemokine receptor CXCR4. Semin Cancer Biol. 4:171–179. 2004. View Article : Google Scholar : PubMed/NCBI | |
Kucia M, Jankowski K, Reca R, Wysoczynski M, Bandura L, Allendorf DJ, Zhang J, Ratajczak J and Ratajczak MZ: CXCR4-SDF-1 signalling, locomotion, chemotaxis and adhesion. J Mol Histol. 35:233–245. 2004. View Article : Google Scholar : PubMed/NCBI | |
Zlotnik A: Chemokines in neoplastic progression. Semin Cancer Biol. 14:181–185. 2004. View Article : Google Scholar | |
Burger JA and Kipps TJ: CXCR4: a key receptor in the crosstalk between tumor cells and their microenvironment. Blood. 107:1761–1767. 2006. View Article : Google Scholar : PubMed/NCBI | |
Wang L, Wang L, Yang B, Yang Q, Qiao S, Wang Y and Sun Y: Strong expression of chemokine receptor CXCR4 by renal cell carcinoma cells correlates with metastasis. Clin Exp Metastasis. 26:1049–1054. 2009. View Article : Google Scholar : PubMed/NCBI | |
Reckamp KL, Strieter RM and Figlin RA: Chemokines as therapeutic targets in renal cell carcinoma. Expert Rev Anticancer Ther. 8:887–893. 2008. View Article : Google Scholar : PubMed/NCBI | |
Darnell JE Jr, Kerr IM and Stark GR: Jak-STAT pathways and transcriptional activation in response to IFNs and other extracellular signaling proteins. Science. 264:1415–1421. 1994. View Article : Google Scholar : PubMed/NCBI | |
Shang D, Liu Y, Ito N, Kamoto T and Ogawa O: Defective Jak-Stat activation in renal cell carcinoma is associated with interferon-alpha resistance. Cancer Sci. 98:1259–1264. 2007. View Article : Google Scholar : PubMed/NCBI | |
Lundqvist A, Abrams SI, Schrump DS, Alvarez G, Suffredini D, Berg M and Childs R: Bortezomib and depsipeptide sensitize tumors to tumor necrosis factor-related apoptosis-inducing ligand: a novel method to potentiate natural killer cell tumor cytotoxicity. Cancer Res. 66:7317–7325. 2006. View Article : Google Scholar | |
Wu XX, Zeng Y, Jin XH and Kakehi Y: Enhanced susceptibility of adriamycin-treated human renal cell carcinoma cells to lysis by peripheral blood lymphocytes and tumor infiltrating lymphocytes. Oncol Rep. 18:353–359. 2007. |