Clinical Cancer Research CTRC-AACR San Antonio Breast Cancer Symposium Infection and Cancer: Biology, Therapeutics, and Prevention
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Cell Growth & Differentiation

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kosari, F.
Right arrow Articles by Vasmatzis, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kosari, F.
Right arrow Articles by Vasmatzis, G.
Clinical Cancer Research Vol. 11, 5128-5139, July 15, 2005
© 2005 American Association for Cancer Research


Imaging, Diagnosis, Prognosis

Clear Cell Renal Cell Carcinoma: Gene Expression Analyses Identify a Potential Signature for Tumor Aggressiveness

Farhad Kosari1, Alexander S. Parker4, Dagmar Marie Kube1, Christine M. Lohse2, Bradley C. Leibovich3, Michael L. Blute3, John C. Cheville1 and George Vasmatzis1

Authors' Affiliations: Departments of 1 Laboratory Medicine and Pathology, 2 Health Sciences Research, and 3 Urology, Mayo Clinic, Rochester, Minnesota, and 4 Department of Health Sciences Research, Mayo Clinic, Jacksonville, Florida

Requests for reprints: George Vasmatzis or John C. Cheville, Department of Laboratory Medicine and Pathology, 200 First Street SW, Rochester, MN 55905. Phone: 507-266-4617 and 507-284-3867; E-mail: vasmatzis.george{at}mayo.edu.


    Abstract
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
Purpose: The objective of this study was to use gene expression profiling to identify novel biomarkers that are predictive of aggressive behavior in clear cell renal cell carcinoma (CCRCC).

Experimental Design: Candidate genes were discovered using Human Genome U133 Plus 2 Arrays and validated on independent samples by quantitative reverse transcription-PCR (RT-PCR). Both the discovery and the validation cohorts included nonaggressive primary CCRCC, aggressive primary CCRCC, metastatic CCRCC, and nonneoplastic kidney adjacent to tumor.

Results: Aggressive primary and metastatic CCRCC displayed no significant differences in gene expression. In contrast, we identified significant differences in gene expression between nonaggressive and aggressive CCRCC (including metastatic CCRCC). Thirty-four of the 35 transcripts that displayed the most significant differential expression by microarray analysis also displayed significant differential expression in independent validation studies using quantitative RT-PCR (P < 0.001 for 31 candidates and P < 0.005 for the remaining three candidates). Hierarchical clustering of the quantitative RT-PCR data using our candidate markers accurately grouped 88% (23 of 26) of aggressive and metastatic CCRCC samples, 100% (14 of 14) of nonaggressive CCRCC samples, and 100% (15 of 15) of nonneoplastic samples into separate clusters. Finally, we evaluated the ability of protein expression levels of one of our candidate markers (survivin) to predict survival among a cohort of 183 CCRCC patients treated surgically at Mayo Clinic from 1990 to 1992. In multivariate analysis, expression of survivin (BIRC5) was inversely associated with cancer-specific survival (P = 0.017).

Conclusion: We used a combination of genomic profiling and validation by quantitative PCR to identify a panel of candidate biomarkers for determining CCRCC aggressiveness. Our data also indicate that the gene expression alterations that result in aggressive behavior and metastatic potential can be identified in the primary tumor.


Incidence and mortality rates for renal cell carcinoma (RCC) are increasing in the United States and these trends do not seem to be explained by the increased use of abdominal imaging (1). The American Cancer Society estimates that in 2005 there will be >36,000 new cases of RCC and 12,660 deaths attributed to RCC (2). Whereas RCC encompasses a group of at least five unique histologic subtypes, the vast majority (>80%) are classified as clear cell RCC (CCRCC; ref. 3).

The standard of care for patients presenting with localized CCRCC remains surgical excision. Unfortunately, whereas the majority of patients with localized CCRCC will be cured by surgery, ~30% will develop metastases and die following removal of a confined tumor. Indeed, a significant problem in the clinical management of patients presenting with localized CCRCC is the inability to determine tumor aggressiveness and accurately predict which patients are at greater risk of experiencing distant metastases following surgery. Recently, several groups of investigators have combined individual clinical and pathologic predictors of CCRCC aggressiveness into composite multivariate scoring systems that possess a high degree of prognostic ability (47). Whereas these scoring systems have improved our ability to predict CCRCC aggressiveness, they remain surrogate measures of the true underlying molecular alterations that are associated with CCRCC aggressiveness. This limitation underscores the need to identify molecular markers of CCRCC aggressiveness within the tumor that have the potential to not only improve outcome prediction but also inform on new targets for therapies to be used in combination with surgical excision. Unfortunately, to date there are no widely accepted molecular biomarkers for CCRCC aggressiveness.

In recent years, gene expression profiling has emerged as a means of rapidly identifying candidate genes and gene expression patterns that are associated with aggressive tumors (8, 9) including CCRCC (10, 11). To date however, comprehensive gene expression profiling studies of well-characterized CCRCCs with the express purpose of identifying markers of aggressiveness are lacking. Moreover, the studies that do exist often omit analyses of normal and metastatic tumor samples or rely on CCRCC cell lines rather than actual patient tumors (12). More importantly, findings from gene profiling studies of CCRCC aggressiveness are often not validated on an independent cohort of samples using more accurate means of evaluating gene expression (i.e., quantitative PCR).

The objective of this study was to first do microarray analysis on a group of aggressive and nonaggressive CCRCC tumors to identify candidate markers of aggressiveness and then validate these candidates using quantitative PCR on an independent set of tumors. By comparing gene expression profiles of nonaggressive CCRCC and aggressive CCRCC, a panel of biomarkers with potential prognostic significance was identified. The potential prognostic significance of the set of genes identified by our microarray study was validated using quantitative reverse transcription-PCR (RT-PCR) on an independent cohort of nonneoplastic kidney and well-characterized CCRCC samples. Finally, the prognostic significance of one of these biomarkers was assessed at the protein level using immunoperoxidase staining methods.


    Materials and Methods
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
Patient and tissue selection. After the Mayo Clinic Institutional Review Board and Human Biospecimens Committee approval, fresh snap-frozen CCRCC tumor and nonneoplastic kidney samples were selected from the Mayo Clinic RCC Biospecimens Resource. Nonaggressive and aggressive CCRCC samples were identified as follows. Patients with aggressive CCRCC died of the disease or developed metastases within 4 years following nephrectomy. Two outcome prediction models were used to identify patients with nonaggressive CCRCC due to lack of sufficient follow-up. One of these models, the SSIGN score, uses pathologic features predictive of cancer-specific survival in CCRCC: tumor size, tumor-node-metastasis (TNM) stage, nuclear grade, and tumor necrosis (13). The second model uses the same pathologic features, except the outcome assessed is metastasis-free survival and not cancer-specific survival (14). The pathologic features required for these models are obtained from the nephrectomy specimen. Patients with nonaggressive CCRCC had SSIGN scores ≤3 and consequently had a predicted 5-year cancer-specific survival in excess of 88% and also had a predicted 5-year metastasis-free survival in excess of 84%. In fact, all but one of the patients with nonaggressive CCRCC had expected survivals in excess of 95%, and none of these patients have developed metastases or died of tumor (Table 1).


View this table:
[in this window]
[in a new window]
 
Table 1. Pathologic features for specimens used in the study

 
The following groups were selected for the microarray experiments: 10 tumors from patients with nonaggressive CCRCC and nine tumors from patients with aggressive CCRCC. In addition, nine metastatic CCRCC tumor samples and 12 nonneoplastic kidney samples were studied. The metastatic tumor specimens included four cases that were matched with primary CCRCC. A separate cohort of patient tumor samples was selected for validation by quantitative RT-PCR using the same criteria as used for the microarray experiments. This validation cohort consisted of 14 patients with nonaggressive CCRCC, 17 patients with aggressive CCRCC, and nine metastatic samples. Also included in the validation study were 15 samples of adjacent nonneoplastic tissue from eight cases with nonaggressive tumor and seven cases with aggressive tumor. The predicted metastasis-free and cancer-specific survivals for the nonaggressive tumors and the pathologic features for all cases used in this study are presented in Table 1.

Before all experiments, H&E-stained sections from frozen tissue blocks were reviewed by a urologic pathologist (J.C.C.) to insure appropriate tissue diagnosis as well as quality and quantity of the tumor samples. Frozen tissue sections were also reviewed for pathologic features predictive of tumor aggressive behavior (nuclear grade and necrosis). Tumor blocks were selected to insure that aggressive CCRCC samples were predominantly high grade (nuclear grades 3 and 4) and nonaggressive CCRCC were all low grade (nuclear grades 1 and 2). In tumor blocks, all nonneoplastic tissue was removed from the frozen block before any analyses.

Oligonucleotide microarray experiments. Microarray experiments were done by the Mayo Clinic Microarray Core Facility; 30 mm3 of each tissue were sectioned at 20 or 35 µm, collected in buffer RLT (Qiagen, Valencia, CA) supplemented with ß-mercaptoethanol and homogenized using a PT 1200C (Kinematica AG, Luzerne, Switzerland) rotor/stator homogenizer. Total RNA was isolated using the RNeasy kit (Qiagen) following manufacturer's specifications. Quality and quantity of RNA samples were analyzed by spectrophotometry and Agilent 2100 Bioanalyzer (Agilant Technologies, Palo Alto, CA). Synthesis of biotin-labeled cRNA, hybridization, washing, staining, and scanning were done following manufacture's protocols (Affymetrix Corp., Santa Clara, CA). Microarray experiments were carried out using the Human Genome U133 Plus 2 Arrays.

Microarray data analysis. Affymetrix microarray analysis software GCOS was used to process scanned chip images. The software generates a cell intensity file for each chip, which contains a single intensity value for each probe cell (.CEL file). dChip 1.3 (http://biosun1.harvard.edu/complab/dchip/) was used to calculate model-based expression index after data from all chips were normalized against an array with median overall intensity using invariant set method (15). The model-based expression index was calculated using perfect match/mismatch model with outlier detection and correction, and the calculated expression values were log2 transformed.

To identify differentially expressed genes in aggressive and nonaggressive cases, three algorithms were used. First, using the dChip program, probe sets with a difference of >2.0 on the log scale (>4.0-fold change) between the average expression levels of the aggressive and nonaggressive cases and a P < 0.001 were identified (129 genes, list 1). To estimate the number of false positives in this list, the cases were randomly assigned to two groups 1,000 times and the same criteria were applied to identify differentially expressed genes. The median false discovery rate (FDR) by this process was 0.8% (one gene) and a 90th percentile of 3.9% (five genes). Second, expression values of probe sets (11,715) determined by the dChip program to be most variable across the aggressive and nonaggressive cases were exported to GeneCluster 2.0 (Whitehead/MIT for Genome Research) to identify 120 probe sets with highest signal to noise ratios (list 2). The signal to noise ratio estimate, also called the discriminate score (10), is computed as SNR = (µ1 µ2) / ({sigma}1 + {sigma}2), where µ and {sigma} refer to the mean and SD, respectively. A high SNR typically suggests that the expression level of a gene displays a much larger variation between the two groups compared with the variation within each group. Finally, probe set expression levels (54,605) from dChip were imported to the prediction analysis of microarray (PAM) algorithm to identify 117 genes that best distinguish aggressive and nonaggressive cases (list 3). PAM uses the "shrunken centroid" approach to reduce the effects of "noisy" genes (16). The threshold for shrinking the centroids was set at 3.5, producing an overall misclassification error of ~11% (3 of 28). Probe sets common to the three lists were identified. From this list, candidates with >30% absent calls in the group determined to overexpress the gene were discarded. Finally, the redundant probe sets representing a gene were removed. The final list (Table 2) included 35 probe sets. We used this list for supervised clustering in the dChip program using the centroid linkage method and Euclidean distance metric (Fig. 2).


View this table:
[in this window]
[in a new window]
 
Table 2. Candidate biomarkers predictive of CCRCC aggressive behavior

 


View larger version (63K):
[in this window]
[in a new window]
 
Fig. 2. Heat map depicting the expression levels of the 35 genes selected using three algorithms. High (red) and low (blue) expression levels according to the scale at the bottom of the heat map. Up-regulated (red bars) and down-regulated (blue bars) genes in primary tumors with nonaggressive CCRCC compared with the aggressive primaries and metastatic CCRCC, respectively. Dendogram (top) of the supervised clustering results based on the 35 selected genes. Gene names and the corresponding Affymetrix probe sets are on the right side of the heat map. Legend colors are as defined in Fig. 1.

 
Quantitative reverse transcription-PCR. Validation experiments were done using tissue obtained from an independent cohort from the RCC Biospecimens Resource. Total RNA isolation and DNase treatment were carried out using RNeasy Mini kit and RNase-Free DNase Set (Qiagen) following manufacturer's specifications. RNA integrity was assessed using the Agilent 2100 Bioanalyzer. One hundred and sixty nanograms of total RNA as measured by spectrophotometry (Nanodrop, Wilmington, DE) were used in reverse transcription using Superscript III reverse transcriptase enzyme (Invitrogen, Carlsbad, CA) following manufacturer's protocol.

Quantitative RT-PCR experiments were done on an ABI 7900 HT system (Applied Biosystems, Foster City, CA). For each primer set, the optimum primer concentration (typically 0.15 µmol/L final concentration) was determined and standard curves were generated using four to five dilutions of a pooled cDNA sample from the validation cohort. Typical standard curve included 4, 1, 0.25, 0.0625, 0.0156, and 0 ng (no template control) of total RNA equivalents of cDNA. To confirm that the amplification occurred on the target sequences, the amplicons were analyzed by gel electrophoresis and the dissociation curves were examined for the presence of a single sharp peak at the melting temperature of the amplicon. The expression level of each gene was normalized by karyopherin {alpha} 6 (KPNA6) as {Delta}CT = CT-KPNA6CT-gene, where CT is the threshold cycle in the quantitative PCR experiment (see Results). To select the most significantly differentially expressed genes (Fig. 4), we used the z score from the Mann-Whitney test (http://faculty.vassar.edu/lowry/utest.html).



View larger version (40K):
[in this window]
[in a new window]
 
Fig. 4. Quantitative RT-PCR validations of selected candidate biomarkers. Left, 10 genes with most significantly down-regulated expression in aggressive and metastatic CCRCC compared with nonaggressive CCRCC. Right, three genes with most significantly up-regulated expression in aggressive and metastatic CCRCC compared with nonaggressive CCRCC. Color designations are as defined in Fig. 3.

 
Clustering analysis of the quantitative reverse transcription-PCR data. Expression levels of genes were measured by quantitative PCR, normalized to KPNA6 (see Results), and imported into the CLUSTER program. In the CLUSTER program, gene expression levels were mean centered and scaled such that for each gene, the sum of the squares of the values across all samples was set to 1. Next, genes and samples were clustered using the correlation similarity metric and the centroid linkage clustering method. Finally, the TREEVIEW program was used to visualize the results (Fig. 5).



View larger version (65K):
[in this window]
[in a new window]
 
Fig. 5. Quantitative RT-PCR experimental data on the 55 sample validation cohort visualized by the TREEVIEW program. Gene names are listed on the right of the heat map. The dendogram displays the clustering of the cohort by the CLUSTER program. In the map, red and green indicate expression levels higher and lower than the mean expression, respectively. The color scheme for the dendogram labels is as defined in Fig. 1.

 
Survivin immunoperoxidase staining and analysis. We examined one of the genes, survivin (BIRC5) that was up-regulated in aggressive CCRCC for its association with cancer-specific survival at the protein level. Immunohistochemical analyses were done on formalin-fixed paraffin-embedded tissue blocks from a large cohort of surgically treated patients with CCRCC with long follow-up. We identified 252 patients treated with open radical nephrectomy or nephron-sparing surgery for CCRCC between 1990 and 1992 from the Mayo Clinic Nephrectomy Registry. Of these, 183 (72.6%) patients had adequate tissue specimens available for immunohistochemical staining and quantification. There were no statistically significant differences in tumor size, TNM classification, nuclear grade, or the presence of tumor necrosis between patients with and without available tissue. Furthermore, there was no statistically significant difference in cancer-specific survival between these two groups (P = 0.804, log-rank test).

Five-micrometer sections were prepared from formalin-fixed paraffin-embedded tissue blocks containing the highest grade. Immunoperoxidase staining with anti-survivin antibody (DAKO, Carpinteria, CA; 1:100 dilution) was done using standard techniques. Slides were reviewed (J.C.C.), and the area of greatest staining was circled (diameter, 1-2 mm) for analysis. Each slide was then scanned using the Bacus Slide Scanner (Bacus Laboratories, Inc., Lombard, IL). The BLISS system (Bremson Company, Lenexa, KS) was used to digitally capture images at 480 x 752 pixel resolution and at 40x magnification. The entire slide is composed of multiple tiles to create a mosaic or composite picture. Computer assisted analysis was done by using the IHCScore software (Bacus Laboratories). This software allows for analysis measurements of (a) total area, (b) immunohistochemistry area, (c) total immunohistochemistry stain density, (d) average immunohistochemistry stain density, and (e) immunohistochemistry percentage area (referenced to the total area). For this study, the all filter was selected to capture the staining intensity of anti-survivin antibody. A series of images were viewed to subjectively create an immunohistochemistry and reference threshold. Each slide was segmented to capture the entire area of interest. Thresholds were adjusted for each slide and data were exported to Excel spreadsheet.

Cancer-specific survival was estimated using the Kaplan-Meier method. The duration of follow-up was calculated from the date of surgery to the date of death or last follow-up. Scatter plots of survivin expression versus the difference in observed survival and the survival expected form a Cox proportional hazards regression model (formally known as the Martingale residuals) were used to identify a potential cut point for survivin expression for illustration purposes. The associations of survivin expression with death from CCRCC were evaluated using Cox proportional hazards regression models univariately and after adjusting for pathologic features important in outcome prediction among patients with CCRCC including primary tumor size, TNM classification, nuclear grade, and the presence of tumor necrosis. Survival analyses were done using the SAS software package (SAS Institute, Cary, NC).


    Results
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
Clustering of cases based on the overall gene expression profiles. The overall gene expression profile from the microarray experiments was studied to determine if CCRCC cases could be separated and appropriately classified into aggressive and nonaggressive categories. Genes with variable expression across the samples were selected for unsupervised clustering (Fig. 1A). This analysis identified two major clades. One clade included all of the nonneoplastic cases from patients with nonaggressive and aggressive CCRCC, and the other clade included all the cases of CCRCC. This indicated that the gene expression profile common to all CCRCC is significantly different from the expression profile in nonneoplastic renal tissue. The clade containing the CCRCC cases consisted of two smaller clades. One clade included only the tumors from patients with aggressive CCRCC and the metastatic tumor samples. The other clade included all tumor samples from patients with nonaggressive CCRCC, three cases from the aggressive CCRCC, and two metastatic tumor samples. This distribution of the cohort suggests that gene expression profiles can stratify the majority of patients into the appropriate categories according to tumor behavior.



View larger version (22K):
[in this window]
[in a new window]
 
Fig. 1. A, unsupervised clustering of the 40 cases from the microarray data. Genes (n = 1,787) that were present in at least 50% of the cases and had expression levels that varied by at least 1.2 SD of log intensity unit were used. The clade on the left consists of normal samples (green legends) exclusively, the clade on the right includes two smaller clusters; a cluster on the left consisting of primary tumors in patients with aggressive CCRCC (red legend) and metastatic tumor samples (pink legend), and a cluster on the right consisting mostly of nonaggressive CCRCC primaries (blue legend). B, unsupervised clustering of the nonneoplastic tissues. Genes (n = 1,273) that were present in at least 50% of the cases and had expression levels that varied by at least 1.0 SD of log intensity unit were used. Dark and light green legends indicate nonneoplastic tissues adjacent to aggressive and nonaggressive primaries, respectively. C, unsupervised clustering of the aggressive CCRCC primary and metastasis cases. Genes (n = 1,568) that were present in at least 50% of the cases and had expression levels that varied by at least 1.2 SD of log intensity unit were used.

 
Comparison of the expression profiles of the nonneoplastic tissues adjacent to aggressive and nonaggressive primary tumors. We examined if the expression profile of the nonneoplastic kidney is associated with the aggressive behavior of CCRCC. In the overall unsupervised clustering plot (Fig. 1A), nonneoplastic tissues clustered together irrespective of the behavior of the matched tumor. Examination of the gene expression profile of nonneoplastic renal tissue was done excluding tumor cases. Again, we noted that the four matched nonneoplastic kidney samples from patients with nonaggressive CCRCC were interspersed among the eight nonneoplastic samples from patients with aggressive CCRCC, and there was no significant difference in expression profiles among these groups of nonneoplastic tissue (Fig. 1B). Although the number of nonneoplastic cases adjacent to nonaggressive CCRCC tumor in this study is limited, these analyses suggest that the gene expression in the nonneoplastic kidney is not associated with the behavior of the matched CCRCC.

Comparison of the expression profiles of the aggressive clear cell renal cell carcinoma primary tumors and the clear cell renal cell carcinoma metastatic cases. Expression profiles were examined to determine if they could discriminate between aggressive primary CCRCC and metastatic CCRCC. In the unsupervised clustering analysis, aggressive primary and metastatic CCRCC did not separate into unique clades (Fig. 1A). To insure that the clustering pattern was not influenced by the expression profiles of the nonneoplastic and nonaggressive CCRCC, expression profiles from the aggressive primary CCRCC and metastatic samples were analyzed separately. Here again, it was noted that the aggressive CCRCC cases were interspersed among the metastatic tumor samples (Fig. 1C). We compared the expression profiles of the two groups using the dChip and PAM algorithms. By the dChip algorithm, the number of differentially expressed genes between the two groups was comparable with the number of differentially expressed genes found by randomly assigning the metastatic samples and the aggressive primary CCRCC to two groups. The median FDR was ~100% and the 90th percentile FDR was 300% to 400%, depending on the significance criteria. We used PAM to identify a group of genes that could be used for classification of aggressive primary and CCRCC metastasis cases (see Materials and Methods). The average misclassification error with any possible threshold for "shrinking centroids" was 40% to 60%. These analyses suggested that no set of genes could separate the metastatic and aggressive primary CCRCC into two distinct groups.

Comparison of expression profiles of nonaggressive clear cell renal cell carcinoma against aggressive and metastatic clear cell renal cell carcinoma. Because the aggressive primary and metastatic CCRCC samples showed similar expression profiles, they were grouped together and compared with nonaggressive CCRCC. This increased the statistical power for identification of significantly differentially expressed genes.

To identify probe sets that are most relevant to CCRCC behavior, we used the signal to noise selection criteria and the PAM algorithm in addition to the fold change and P value criteria provided by the dChip software (see Materials and Methods). In each case, comparisons were made for the gene expression values in the primary tumors from nonaggressive CCRCC versus aggressive primary tumors and metastatic CCRCC samples. The top 120 to 130 candidate prognostic biomarkers were selected using the three statistical algorithms. We identified 129 probe sets that displayed a fold change of at least 4.0 and P < 0.001 (median FDR = 0.8% and 90 percentile FDR = 3.9%) by dChip, 125 probe sets with highest signal to noise ratio by GeneCluster, and 120 probe sets by PAM after the centroids were "shrunken" by a factor of 3.75. Probe sets common in the three lists that also had a present (P) call by the dChip algorithm in at least 70% of the cases were selected. Finally, multiple probe sets representing the same gene were discarded so that the listing would represent unique individual genes. The final candidate list included 35 probe sets corresponding to 35 unique transcripts (Table 2). The majority of the candidate biomarkers identified by our analysis (29 of 35, ~83%) displayed down regulation of expression in the aggressive CCRCC compared with the nonaggressive CCRCC.

With this set of 35 differentially expressed genes, hierarchical clustering of the 28 CCRCC tissues was repeated. This analysis produced two major subgroups. One subgroup contained 17 of 18 (94%) of the aggressive CCRCC and metastatic CCRCC samples. The other subgroup included all 10 (100%) tissues from the nonaggressive CCRCCs and the remaining case of aggressive CCRCC (Fig. 2).

Validation by quantitative reverse transcription-PCR. The results from the microarray experiments were validated by determining the expression levels of the 35 candidate biomarkers in an independent cohort of CCRCC and nonneoplastic samples using quantitative RT-PCR. Compared with the microarray technology, quantitative RT-PCR provides a much wider dynamic range (5-6 orders of magnitude) and thus a more accurate measurement of the relative gene expression values.

We first used the microarray data to identify genes that could be used for normalization of expression levels determined by quantitative PCR. Two genes, eukaryotic translation elongation factor 1 {alpha} 1 (EEF1A1) and karyopherin {alpha} 6 (KPNA6), were selected from the five genes with the lowest SDs across all samples in our microarray expression data. In addition, two common housekeeping genes, ß-2-microglobin (B2M) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH), were examined. The expression levels of all four genes were measured by quantitative RT-PCR (Fig. 3). As expected from the microarray data, expression of GAPDH was most variable across the samples followed by B2M. On the other hand, KPNA6 expression was least variable across all samples. More importantly, expression levels of GAPDH and B2M were lower in nonneoplastic kidney than in the RCC cases (P < 0.001 for both genes). In contrast, KPNA6 expression was not statistically different among the CCRCC and nonneoplastic tissues. Furthermore, expression levels of KPNA6 were comparable with those of most of the candidate biomarkers and on average 10- to 20-fold (approximately four cycles in a quantitative RT-PCR experiment) lower than the expression levels of GAPDH. Thus, KPNA6 was selected for normalization of the quantitative RT-PCR data.



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 3. Unnormalized (raw) data depicting the expression values of the four candidate normalization genes across the 55 sample validation cohort. Nonneoplastic cases adjacent to nonaggressive (light green) and aggressive primaries (dark green). Nonaggressive CCRCC primaries (blue), aggressive CCRCC primaries (red), and metastatic cases (pink). GAPDH and B2M display the largest SD and have higher expression in CCRCC tissues than in nonneoplastic samples. KPNA6 expression levels display the lowest variation and do not show differential expression between the tumors and the nonneoplastic cases.

 
All of the candidate biomarkers, except interleukin 8 (IL-8), displayed significant differential expression by quantitative RT-PCR (P < 0.001 for 31 candidates and P < 0.005 for the remaining three candidates), as predicted by our microarray analysis. In the microarray experiments, IL-8 expression was up-regulated in aggressive primary and metastatic CCRCC relative to nonaggressive CCRCC. In the validation cohort, the up-regulation of IL-8 in aggressive primary and metastatic CCRCC was marginal (P <0.055). Figure 4 (left) illustrates expression levels of 10 candidate biomarkers that were most significantly down-regulated in aggressive primary and metastatic CCRCC compared with nonaggressive CCRCC by the Mann-Whitney test. Figure 4 (right) illustrates three candidate biomarkers that showed the highest level of up-regulation in aggressive and metastatic CCRCC compared with nonaggressive CCRCC. Of note, every cycle difference in these experiments represents an ~2-fold difference in expression. For ECRG4, for example, we detected a difference of more than four cycles in the mean expression levels between nonaggressive CCRCC and aggressive CCRCC, indicating an ~16-fold difference in expression levels between the two groups. Hierarchical clustering of the quantitative RT-PCR data confirmed that the 35 genes selected from the microarrays had prognostic significance for CCRCC. As Fig. 5 shows, two main subgroups were identified in the validation cohort. One subgroup included 23 of the 26 aggressive and metastatic CCRCC cases (88%) and none of the nonaggressive CCRCC cases. The other main subgroup included two further clusters, one containing all 15 (100%) of the nonneoplastic tissues and the other containing all 14 (100%) of the nonaggressive cases and the remaining three cases of aggressive primaries.

Survivin immunoperoxidase staining and analysis. At last follow-up, 113 of the 183 patients with adequate tissue specimens available for immunoperoxidase analysis had died, including 63 patients who died of CCRCC at a median 2.4 years following surgery (range, 0-13). Among the 70 patients who were still alive, the median duration of follow-up was 12.5 years (range, 1-15). The median survivin expression was 0.87 (mean, 2.37; range, 0.03-35.80). The staining was predominantly nuclear (Fig. 6A and B). A scatterplot of survivin expression versus the expected risk of death for each patient suggested that a cut point value of 2 would be most appropriate (17). There were 58 (31.7%) patients with survivin expression levels of ≥2. Univariately, the risk ratio for a 1-unit increase in survivin expression was 1.14 (95% confidence interval, 1.10-1.19; P < 0.001). In our categorical analysis, patients with survivin expression levels of ≥2 were over four times more likely to die from CCRCC compared with patients with survivin expression of <2 (risk ratio, 4.27; 95% confidence interval, 2.59-7.05; P < 0.001). The association between survivin expression levels and cancer-specific survival is shown in Fig. 7. The estimated cancer-specific survival rates at 1, 5, and 10 years following surgery were 96.7%, 85.3%, and 75.4%, respectively, for patients with survivin expression levels of <2 compared with 71.2%, 47.6%, and 38.0%, respectively, for patients with survivin expression of ≥2. After multivariate adjustment for primary tumor size, TNM classification, nuclear grade and tumor necrosis, the risk ratio for a 1-unit increase in survivin expression was 1.10 (95% confidence interval, 1.04-1.16; P < 0.001), whereas the risk ratio for survivin expression of ≥2 was 2.07 (95% confidence interval, 1.14-3.76; P = 0.017).



View larger version (125K):
[in this window]
[in a new window]
 
Fig. 6. A, nonaggressive CCRCC with survivin score of <2. Survivin score is nuclear in location. B, aggressive CCRCC with high survivin score and numerous positively staining nuclei.

 


View larger version (12K):
[in this window]
[in a new window]
 
Fig. 7. Cancer-specific survival by survivin expression assessed using immunoperoxidase methods for 183 patients with CCRCC. Patients with tumors that had survivin score of ≥2 had a significantly worse cancer-specific survival than patients with survivin score of <2 (P < 0.001).

 

    Discussion
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
This study included gene expression profiling of well-characterized CCRCC primary tumors, adjacent nonneoplastic kidney, and metastatic tumor samples. Unsupervised clustering using data from the oligonucleotide microarray experiments separated the CCRCC samples from nonneoplastic kidney samples. Furthermore, this analysis separated CCRCC cases into their respective categories indicating unique gene expression profiles predictive of tumor aggressiveness. Using genomic profiling on tumor samples from a well-characterized cohort of CCRCC patients, we identified a panel of genes that were differentially expressed between patients with aggressive and nonaggressive CCRCC. Additionally, primary CCRCC with aggressive behavior did not exhibit a significantly different gene expression profile from metastatic CCRCC. This observation suggests that gene expression alterations that result in aggressive behavior and metastatic potential can be identified in the primary tumor. In subsequent analyses, we identified 35 unique transcripts whose expression values differed significantly between nonaggressive and aggressive CCRCC. Validation studies using quantitative RT-PCR on an independent set of CCRCC tissues confirmed the results from oligonucleotide microarray experiments and further supported this set of genes as potential biomarkers for aggressive CCRCC. Finally, one transcript, survivin (BIRC5), was assessed at the protein level in 183 patients treated surgically for CCRCC. Survivin protein expression was found to be significantly associated with decreased cancer-specific survival even after adjusting for pathologic features predictive of outcome (grade, TNM stage, tumor size, and presence of necrosis).

The pursuit to accurately predict outcome among patients with primary CCRCC as well as metastatic CCRCC has resulted in several recently published scoring systems that incorporate clinical and pathologic features to determine prognosis (6, 7, 13, 18). These prediction models have improved our ability to predict outcome in patients with CCRCC, but the incorporated clinical and pathologic features are surrogate measures of the molecular alterations that determine CCRCC aggressiveness. As such, they do not provide insight regarding the biology of CCRCC nor do they provide potential targets for adjuvant therapies. In this study, the use of CCRCC primary tumors as well as metastatic CCRCC samples allowed for the identification of an expression profile indicative of tumor aggressiveness. In addition, we analyzed one of these biomarkers, survivin, using a immunoperoxidase assay, and found that overexpression of the survivin protein was significantly associated with decreased cancer-specific survival in a multivariate model. Survivin is a member of the inhibitor of apoptosis protein family, and its expression both at the mRNA and protein level has been associated with more aggressive tumor behavior in carcinomas of the larynx, liver, prostate, lung, ovary, and stomach (1922). This is the first study to show that survivin protein overexpression is predictive of outcome in patients with CCRCC. Interestingly, it has been shown that inhibition of survivin activity and down-regulation of nuclear factor-{kappa}B, both mediated by VHL protein, are critical to conferring sensitivity of RCC cells to tumor necrosis factor-{alpha} cytotoxicity (23). This suggests that survivin has potential not only as a prognostic biomarker but also as a therapeutic target, although additional confirmatory studies are needed.

A number of other genes that showed increased expression in aggressive as compared with nonaggressive CCRCC are noteworthy. IL-8, a potent chemotactic cytokine for inflammatory cells, exhibited higher expression levels in the aggressive compared with the nonaggressive CCRCCs in our microarray experiments. IL-8 is implicated in the migration of lymphocytes into tumors through an {alpha}-1 integrin-mediated pathway in the extracellular matrix, and studies show that neutralizing antisera specific to IL-8 inhibit tumor-infiltrating lymphocyte migration (24). Although our microarrray experiments identified significantly increased expression of IL-8 in aggressive primary and metastatic CCRCC, this differential expression was marginally significant (P < 0.055) by our quantitative RT-PCR experiments. Another gene, serum amyloid A, has been identified in the serum of CCRCC patients, and elevated serum levels are associated with aggressive CCRCC (25). Serum amyloid A1 and A2 are acute-phase reactants whose expression is regulated in part by IL-1 and IL-6 (2628). Serum amyloid A can be induced in renal tubular epithelial cells, but previous to this study, serum amyloid A mRNA had not been associated with CCRCC aggressive behavior. CKS2, determined to be up-regulated in aggressive CCRCC, has been associated with cancer (up-regulated in metastatic colon cancer; ref. 29), but its function and significance in CCRCC will require further study. Finally, HSP150 was found to be up-regulated in aggressive CCRCC. This gene has not been associated with CCRCC or any cancer to the best of our knowledge. In contrast to a limited number of up-regulated genes in aggressive CCRCC, there were numerous genes that exhibited decreased mRNA levels relative to nonaggressive CCRCC. Several of these genes have been described previously, yet their functional role in CCRCC remains unknown. Esophageal cancer-related gene 4 (ECRG4) has been identified to be down-regulated in squamous cell carcinoma of the esophagus through hypermethylation of the CpG islands (30, 31). In this study, ECRG4 was identified as one of the most significantly down-regulated genes in aggressive tumors and we are investigating the influence of this gene in progression of CCRCC.

At the present time, there is no standard method for the analysis of microarray data. In this report, we used three algorithms that were relevant to our analysis and identified the best candidate biomarkers common to all three of the algorithms. The fact that all of the candidate biomarkers on the list were validated by the quantitative RT-PCR experiments suggests that our approach for the analysis of microarray data was justified. To identify genes for normalization of quantitative RT-PCR results, we searched our microarray data for transcripts that displayed minimum variation in normalized expression values across the samples. The two transcripts selected by this analysis, EEF1A1 and KPNA6, were confirmed by quantitative RT-PCR to have considerably less variation across the 55-sample validation cohort than the commonly used GAPDH and B2M. Furthermore, GAPDH and B2M had significantly higher expression levels in CCRCC samples than in nonneoplastic kidney. Increased expression of GAPDH mRNA in tumor samples is consistent with reports suggesting increased expression of GAPDH protein in renal cell carcinoma to meet the energy demands of the tumor cells following diminished oxidative phosphorylation in the mitochondria (32). Similarly, increased expression of B2M is consistent with reports indicating elevated levels of B2M protein in the serum of patients with RCC (33). Comparing the expression levels of EEF1A1 and KPNA6, we chose KPNA6 for normalization because the expression levels of KPNA6 across the validation samples were more comparable to the expression levels of the selected biomarkers.

Other investigators have examined gene expression in CCRCC using microarrays. Takahashi et al. examined gene expression profiling in 26 cases of CCRCC and reported distinctive gene expression profiles associated with tumor aggressiveness similar to our findings (10). In that study, 17 CCRCC cases were nonaggressive and nine were aggressive. In the gene expression analyses, tumors clustered into defined outcome groups based on 40 transcripts. However, the expression profile was not validated either on the original cohort or an independent set of tumors. Univariate analyses showed a significant correlation of expression profile with outcome as well as histologic tumor grade. However, the number of patients was limited for these types of analyses, and the selection criteria for inclusion into the study were not provided. It is difficult to compare our analyses to that of Takahashi et al. We employed the latest Affymetrix oligonucleotide arrays, which contain more genes than the cDNA arrays used by Takahashi et al. In addition, patient selection and criteria for determining differential gene expression differed between our study and that of Takahashi et al.

The selection of cases in our study warrants comment. CCRCC samples were selected based on a combination of outcome and pathologic features (Table 1). In cases of nonaggressive CCRCC with limited follow-up data, two prediction models were employed to insure that patients considered to have nonaggressive CCRCC had a predicted 5-year metastasis-free survival of at least 84% and 5-year cancer-specific survival of at least 88%. In addition, all frozen tissue blocks were reviewed to insure that nonaggressive tumors were all low grade (nuclear grades 1 and 2). In contrast, patients with CCRCC considered aggressive died of disease or developed metastases within 4 years of diagnosis. In addition, because of the grade heterogeneity found within many CCRCC, aggressive tumors used in the analysis were predominantly nuclear grades 3 and 4. It is possible that this selection process, which included pathologic features, improved our ability to identify significant differences in gene expression and confirms the importance of pathologic review of all tumors analyzed by microarrays and quantitative RT-PCR.

At least one of the transcripts in our list of differentially expressed genes, endomucin (EMCN), is believed to be associated with the nonepithelial renal components. Because bulk tissues were used in this study, we cannot determine if these genes are expressed in nontumor or tumor cells, although tissues were selected to be composed entirely of tumor with limited stromal and inflammatory components. Studies using laser capture microdissection, quantitative RT-PCR, and immunohistochemistry are under way to localize mRNA and protein expression patterns.

In conclusion, our experimental analyses identified a panel of potential biomarkers that are differentially expressed between the aggressive and the nonaggressive CCRCC. In addition, protein expression of one of these genes, survivin, was an independent predictor of patient outcome in a large cohort of patients with surgically treated CCRCC. Further studies are required to determine if expression of the other genes identified in this study provide prognostic information beyond that provided by routine pathologic examination and prognostic scoring systems and algorithms.


    Footnotes
 
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Received 1/11/05; revised 4/ 5/05; accepted 4/27/05.


    References
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 

  1. Chow WH, Devesa SS, Warren JL, Fraumeni JF Jr. Rising incidence of renal cell cancer in the United States. JAMA 1999;281:1628–31.[Abstract/Free Full Text]
  2. Jemal A, Murray T, Ward E, et al. Cancer statistics, 2005. CA Cancer J Clin 2005;55:10–30.[Abstract/Free Full Text]
  3. Cheville JC, Lohse CM, Zincke H, Weaver AL, Blute ML. Comparisons of outcome and prognostic features among histologic subtypes of renal cell carcinoma. Am J Surg Pathol 2003;27:612–24.[CrossRef][Medline]
  4. Sorbellini M, Kattan MW, Snyder ME, et al. A postoperative prognostic nomogram predicting recurrence for patients with conventional clear cell renal cell carcinoma. J Urol 2005;173:48–51.[Medline]
  5. Frank I, Blute ML, Cheville JC, et al. A multifactorial postoperative surveillance model for patients with surgically treated clear cell renal cell carcinoma. J Urol 2003;170:2225–32.[Medline]
  6. Kattan MW, Reuter V, Motzer RJ, Katz J, Russo P. A postoperative prognostic nomogram for renal cell carcinoma. J Urol 2001;166:63–7.[CrossRef][Medline]
  7. Zisman A, Pantuck AJ, Dorey F, et al. Improved prognostication of renal cell carcinoma using an integrated staging system. J Clin Oncol 2001;19:1649–57.[Abstract/Free Full Text]
  8. Alizadeh AA, Ross DT, Perou CM, van de Rijn M. Towards a novel classification of human malignancies based on gene expression patterns. J Pathol 2001;195:41–52.[CrossRef][Medline]
  9. Alizadeh AA, Eisen MB, Davis RE, et al. Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature 2000;403:503–11.[CrossRef][Medline]
  10. Takahashi M, Rhodes DR, Furge KA, et al. Gene expression profiling of clear cell renal cell carcinoma: gene identification and prognostic classification. Proc Natl Acad Sci U S A 2001;98:9754–9.[Abstract/Free Full Text]
  11. Vasselli JR, Shih JH, Iyengar SR, et al. Predicting survival in patients with metastatic kidney cancer by gene-expression profiling in the primary tumor. Proc Natl Acad Sci U S A 2003;100:6958–63.[Abstract/Free Full Text]
  12. Liou LS, Shi T, Duan ZH, et al. Microarray gene expression profiling and analysis in renal cell carcinoma. BMC Urol 2004;4:9.[Medline]
  13. Frank I, Blute ML, Cheville JC, Lohse CM, Weaver AL, 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–400.[CrossRef][Medline]
  14. Leibovich BC, Blute ML, Cheville JC, et al. 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–71.[CrossRef][Medline]
  15. Li C, Wong WH. Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci U S A 2001;98:31–6.[Abstract/Free Full Text]
  16. Tibshirani R, Hastie T, Narasimhan B, Chu G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci U S A 2002;99:6567–72.[Abstract/Free Full Text]
  17. Therneau T, Grambsch P. Modeling survival data: extending the Cox model. Edition 1. Ann Arbor: Springer-Verlag; 2000. p. 87–92.
  18. Leibovich BC, Han KR, Bui MH, et al. Scoring algorithm to predict survival after nephrectomy and immunotherapy in patients with metastatic renal cell carcinoma: a stratification tool for prospective clinical trials. Cancer 2003;98:2566–75.[CrossRef][Medline]
  19. Kren L, Brazdil J, Hermanova M, et al. Prognostic significance of anti-apoptosis proteins survivin and bcl-2 in non-small cell lung carcinomas: a clinicopathologic study of 102 cases. Appl Immunohistochem Mol Morphol 2004;12:44–9.[Medline]
  20. Pizem J, Cor A, Gale N. Survivin expression is a negative prognostic marker in laryngeal squamous cell carcinoma and is associated with p53 accumulation. Histopathology 2004;45:180–6.[Medline]
  21. Shariat SF, Lotan Y, Saboorian H, et al. Survivin expression is associated with features of biologically aggressive prostate carcinoma. Cancer 2004;100:751–7.[CrossRef][Medline]
  22. Miyachi K, Sasaki K, Onodera S, et al. Correlation between survivin mRNA expression and lymph node metastasis in gastric cancer. Gastric Cancer 2003;6:217–24.[CrossRef][Medline]
  23. Qi H, Ohh M. The von Hippel-Lindau tumor suppressor protein sensitizes renal cell carcinoma cells to tumor necrosis factor-induced cytotoxicity by suppressing the nuclear factor-{kappa}B-dependent antiapoptotic pathway. Cancer Res 2003;63:7076–80.[Abstract/Free Full Text]
  24. Ferrero E, Fabbri M, Poggi A, Galati G, Bernasconi S, Zocchi MR. Tumor-driven matrix invasion by infiltrating lymphocytes: involvement of the {alpha}1 integrin I-domain. Eur J Immunol 1998;28:2530–6.[CrossRef][Medline]
  25. Kimura M, Tomita Y, Imai T, et al. Significance of serum amyloid A on the prognosis in patients with renal cell carcinoma. Cancer 2001;92:2072–5.[CrossRef][Medline]
  26. Glojnaric I, Casl MT, Simic D, Lukac J. Serum amyloid A protein (SAA) in colorectal carcinoma. Clin Chem Lab Med 2001;39:129–33.[CrossRef][Medline]
  27. Blay JY, Rossi JF, Wijdenes J, et al. Role of interleukin-6 in the paraneoplastic inflammatory syndrome associated with renal-cell carcinoma. Int J Cancer 1997;72:424–30.[CrossRef][Medline]
  28. Raynes JG, McAdam KP. Serum amyloid A isoforms in inflammation. Scand J Immunol 1991;33:657–66.[Medline]
  29. Li M, Lin YM, Hasegawa S, et al. Genes associated with liver metastasis of colon cancer, identified by genome-wide cDNA microarray. Int J Oncol 2004;24:305–12.[Medline]
  30. Lu J, Liu Z, Xiong M, et al. Gene expression profile changes in initiation and progression of squamous cell carcinoma of esophagus. Int J Cancer 2001;91:288–94.[CrossRef][Medline]
  31. Yue CM, Deng DJ, Bi MX, Guo LP, Lu SH. Expression of ECRG4, a novel esophageal cancer-related gene, downregulated by CpG island hypermethylation in human esophageal squamous cell carcinoma. World J Gastroenterol 2003;9:1174–8.[Medline]
  32. Cuezva JM, Krajewska M, de Heredia ML, et al. The bioenergetic signature of cancer: a marker of tumor progression. Cancer Res 2002;62:6674–81.[Abstract/Free Full Text]
  33. Selli C, Cozzolino F, Carini M, Lenzi R, Vercelli D. Serum ß 2 microglobulin levels in patients with renal cell carcinoma. Urol Res 1984;12:261–3.[Medline]



This article has been cited by other articles:


Home page
Clin. Cancer Res.Home page
A. E. Krambeck, H. Dong, R. H. Thompson, S. M. Kuntz, C. M. Lohse, B. C. Leibovich, M. L. Blute, T. J. Sebo, J. C. Cheville, A. S. Parker, et al.
Survivin and B7-H1 Are Collaborative Predictors of Survival and Represent Potential Therapeutic Targets for Patients with Renal Cell Carcinoma
Clin. Cancer Res., March 15, 2007; 13(6): 1749 - 1756.
[Abstract] [Full Text] [PDF]


Home page
Clin. Cancer Res.Home page
J. Jones and T. A. Libermann
Genomics of Renal Cell Cancer: The Biology Behind and the Therapy Ahead
Clin. Cancer Res., January 15, 2007; 13(2): 685s - 692s.
[Abstract] [Full Text] [PDF]


Home page
Clin. Cancer Res.Home page
J. S. Lam, A. J. Pantuck, A. S. Belldegrun, and R. A. Figlin
Protein Expression Profiles in Renal Cell Carcinoma: Staging, Prognosis, and Patient Selection for Clinical Trials
Clin. Cancer Res., January 15, 2007; 13(2): 703s - 708s.
[Abstract] [Full Text] [PDF]


Home page
Brief BioinformHome page
J. C. Hsu, J. Chang, T. Wang, E. Steingrimsson, M. K. Magnusson, and K. Bergsteinsdottir
Statistically designing microarrays and microarray experiments to enhance sensitivity and specificity
Brief Bioinform, January 1, 2007; 8(1): 22 - 31.
[Abstract] [Full Text] [PDF]


Home page
Clin. Cancer Res.Home page
K.-i. Shioi, A. Komiya, K. Hattori, Y. Huang, F. Sano, T. Murakami, N. Nakaigawa, T. Kishida, Y. Kubota, Y. Nagashima, et al.
Vascular Cell Adhesion Molecule 1 Predicts Cancer-Free Survival in Clear Cell Renal Carcinoma Patients
Clin. Cancer Res., December 15, 2006; 12(24): 7339 - 7346.
[Abstract] [Full Text] [PDF]


Home page
JCOHome page
J. S. Lam, A. Breda, A. S. Belldegrun, and R. A. Figlin
Evolving Principles of Surgical Management and Prognostic Factors for Outcome in Renal Cell Carcinoma
J. Clin. Oncol., December 10, 2006; 24(35): 5565 - 5575.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kosari, F.
Right arrow Articles by Vasmatzis, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kosari, F.
Right arrow Articles by Vasmatzis, G.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Cell Growth & Differentiation