
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
Authors' Affiliations: Departments of 1 Pathology and Laboratory Medicine and 2 Public Health, Division of Biostatistics and Epidemiology, Weill Medical College, Cornell University, New York, New York
Requests for reprints: Yao-Tseng Chen, Department of Pathology, Weill Medical College, Cornell University, Room C-458A, 1300 York Avenue, New York, NY 10021. Phone: 212-746-6472; Fax: 1-212-746-4483; E-mail: ytchen{at}med.cornell.edu.
| Abstract |
|---|
|
|
|---|
Experimental Design: The expression of four genes in renal tumor was evaluated by quantitative reverse transcription-PCR: carbonic anhydrase IX (CA9), methylacyl-CoA racemase (AMACR), parvalbumin (PVALB), and chloride channel kb (CLCNKB). Thirty-one fresh-frozen and 63 formalin-fixed, paraffin-embedded tumor specimens were analyzed.
Results: CA9 expression was highest in clear cell carcinoma and lowest in chromophobe RCC and in oncocytoma. AMACR expression was highest in papillary RCC, and CLCNKB was highest in chromophobe RCC/oncocytoma. PVALB was highest in chromophobe RCC, variable in oncocytoma, and low in clear cell and papillary types. Similar findings were observed in fresh-frozen and formalin-fixed specimens. The mRNA expression ratios among these genes (i.e., CA9/AMACR and AMACR/CLCNKB ratios) further accentuate the gene expression differences among these tumors, and a molecular diagnostic algorithm was established. This algorithm accurately classified the 31 fresh-frozen tumors into 14 clear cell, 5 papillary, 6 chromophobe, and 6 oncocytomas. In the formalin-fixed group, the molecular criteria accurately classified the cases into 15 clear cell, 16 papillary, and 32 in the chromophobe/oncocytoma group but could only separate some, but not all, oncocytomas from chromophobe RCC.
Conclusions: RNA expression ratios based on the four-gene panel can accurately classify subtypes of RCC as well as help distinguish some oncocytomas from chromophobe RCC.
5% are of the chromophobe type (2). The distinction of these different types is clinically important, as papillary and chromophobe types have a better prognosis than the clear cell type (3). For instance, Thoenes et al. (4) reported a 5-year overall survival rate of 92% for chromophobe in contrast to 62% for clear cell type. Similarly, a 5-year survival rate of 82% to 90% has also been reported for papillary renal carcinoma (5, 6). Although the histologic diagnosis of these different subtypes of RCC can be made unequivocally in most cases, selected cases can be difficult, particularly the chromophobe type. Moreover, the differential diagnosis between RCC and renal oncocytoma can pose a diagnostic dilemma. Oncocytoma, comprising 3% to 5% of the renal tumors, is generally considered a benign neoplasm (7). However, some cases of oncocytoma can show atypical histologic features, leading to their being interpreted as RCC, either as chromophobe type or as RCC, unclassified, for fear of underdiagnosing a malignant neoplasm. To help resolve these morphologic dilemmas in surgical pathology, attempts have been made to segregate RCC types by immunohistochemical staining, with variable results (e.g., refs. 810).
An alternative to immunohistochemical staining would be to distinguish these tumors at the DNA or RNA level. Despite their overlapping histologic features, different RCC subtypes are biologically distinct. This is evident from their characteristic cytogenetic abnormalities, the most common ones being the chromosome 3p aberrations in clear cell type, polysomes in papillary RCC, most often trisomy 17 and trisomy or tetrasomy 7, and multiple chromosomal losses in chromophobe type (6). These biological differences imply that these tumors should bear distinguishable molecular signatures that could be useful diagnostic tools. This notion is supported by recent gene profiling studies (1115). In compilation, these studies indicated that clear cell, papillary, and chromophobe RCC could be reliably separated by cDNA microarray technique and hierarchical clustering analysis. However, chromophobe RCC and oncocytoma were found to closely clustered, reflecting their common origin from the collecting duct (6).
Although current evidence indicates that accurate classification of RCC can be achieved by DNA microarray analysis, the complexity and the cost of the assay exclude its use as a routine diagnostic test at present. In this study, we show that by quantitative reverse transcription-PCR (qRT-PCR) analysis of four genes, carbonic anhydrase IX (CA9), methylacyl-CoA racemase (AMACR), chloride channel kb (CLCNKB), and parvalbumin (PVALB) the histologic subtypes of the RCCs could be accurately predicted. In addition, oncocytomas can be separated from clear cell and papillary types in all cases and distinguished from chromophobe RCC in some cases.
| Materials and Methods |
|---|
|
|
|---|
In silico analysis to identify candidate genes for quantitative reverse transcription-PCR and primer design. To identify the most promising candidates from the selected gene list, the expression profile of each gene in normal tissues was evaluated using a combination of the Serial Analysis of Gene Expression Anatomic Reviewer and its Virtual Northern tool (http://cgap.nci.nih.gov/SAGE/AnatomicViewer) and database searches using BLASTN (http://www.ncbi.nlm.nih.gov/BLAST). The objective of the analysis was to identify genes with limited expression in normal tissues, focusing on genes with exclusive or predominant expression in kidney. Once such a gene is identified, the exon-intron structure of the gene was defined, and trans-intronic primers for conventional as well as real-time RT-PCR were designed using Primer3 software (http://frodo.wi.mit.edu/primer3/primer3_code.html).
Conventional PCR primers were custom synthesized by Qiagen (Valencia, CA), whereas the qRT-PCR primers and probes were synthesized by Applied Biosystems (Foster City, CA) with 5'-FAM and 3'-TAMRA labeling. The primer sequences are listed in Table 1. Primer-probe premix for glyceraldehyde-3-phosphate dehydrogenase (GAPDH) endogenous control was purchased from Applied Biosystems.
|
30 mg) were used for each sample, with the final RNA yield in the range of 10 to 60 µg based on A260 nm by spectrophotometer. For extraction of RNA from paraffin-embedded tissue, the Optimum FFPE RNA isolation kit (Ambion, Austin, TX) was used using materials derived from four 8-µm sections. The quantity of RNA isolated ranged from 0.14 to 9.70 µg. Conventional reverse transcription-PCR. Total RNA (2 µg) isolated from fresh-frozen tissues was used for a 20 µL reverse transcription reaction. Synthesized cDNA (2 µL) was then used per 25 µL PCR assay. The PCR was set up using AmpliTaq Gold and its corresponding 10x PCR mix (Applied Biosystems) in an ABI 9600, with 35 cycles of amplification, each cycle consisting of 15 seconds at 94°C, 1 minute at 60°C, and 1 minute at 72°C. The PCR products were visualized by 1.2% agarose gel electrophoresis and ethidium bromide staining.
Quantitative reverse transcription-PCR. qRT-PCR was done using an ABI PRISM 7000 Sequence Detection System. Total RNA (up to 2 µg) was used per 20 µL reverse transcription reaction, and cDNA (2 µL) was then used for each 25 µL PCR, corresponding to 100 ng total RNA. Forty-five two-step cycles of amplification were done, each cycle consisting of 15 seconds at 95°C and 1 minute at 60°C. Control amplification for housekeeping gene GAPDH was done in all samples. All fresh-frozen specimens included in the final analysis have Ct values differing by less than two cycles, indicating similar cDNA quality and quantity. One case of papillary carcinoma had much poorer GAPDH amplification possibly due to inclusion of necrotic areas in the tissue used for RNA extraction. This case was excluded from further study. The GAPDH amplification results were more variable for paraffin-embedded samples, but all samples successfully amplified and were included in the final analysis.
Analysis of quantitative reverse transcription-PCR data. Following the amplification, the same threshold was set for analyzing all experiments to compare Ct values derived from different experiments. The mean Ct values from each sample were normalized against the corresponding GAPDH Ct values calculated as (Ctexperimental gene CtGAPDH). All data in Results represented normalized Ct values.
For comparing the expression levels of individual gene in various types of tumor, the normalized Ct values were plotted using box-whiskers plots according to each tumor type. These plots showed the median and the range of the Ct values as well as Ct values corresponding to 25% and 75% of the specimens. As such, these plots allowed more quantitative comparisons between groups than other graphic tools, such as scatter plots.
In addition to evaluating individual gene expression difference among tumor types, the relative expression between paired genes was found to be useful predictors of the tumor type (see Results). For this comparison, the relative expression of the paired genes was calculated as (raw CtGene A raw CtGene B) and expressed as
Ct (gene A gene B) or (gene A gene B) for short. This
Ct value reflected directly the ratio between mRNA levels of the two genes, which equaled 2
Ct if these two genes were amplified with similar efficiencies. For example, a tumor with (CA9 AMACR) value of 4.0 would have CA9 mRNA level 16-fold (24) of the AMACR mRNA level at the same specimen, as CA9 and AMACR did amplify with similar efficiency.
Statistical analysis and modeling. Regression tree method (16, 17) was used to create the decision rules for classification of renal tumors based on the expression ratios of genes chosen as predictors (see Results), with the goal being to provide highest classification accuracy using the known histologic diagnosis as the "gold standard."
Once the model was established, grouped cross-validation is used to estimate the prediction accuracy and model fitting accuracy. In this analysis, the data are divided into 12 similar size groups. Sequentially, each of the 12 groups is treated as the validation sample and the other 11 groups are treated as training sample. The model is fitted based on data from the 11 groups and the prediction was made for the 12th group. This process was repeated 20 times to compute the average misclassification rate in the training sample and the validation sample.
| Results |
|---|
|
|
|---|
|
|
CA9 was strongly expressed in clear cell RCC, with moderate expression in papillary carcinoma, and low expression in chromophobe RCC and in oncocytoma. A representative qRT-PCR amplification plot is shown in Fig. 2. Shown in normalized Ct (see Materials and Methods), for which a lower Ct value indicates a higher expression, the median Ct for clear cell RCC was 5.612 (range, 4.627-8.022) in contrast to 19.410 for chromophobe RCC and oncocytoma combined (range, 17.109-22.334) and 12.042 for papillary renal carcinoma (Table 3; Fig. 3). As the fold difference at levels of mRNA expression equal 2
Ct in this experiment, the Ct range difference indicated that the lowest CA9 mRNA expression in clear cell type is >512 (29)fold that of the highest level observed in chromophobe/oncocytoma, conclusively distinguishing clear cell carcinoma from chromophobe RCC/oncocytoma. In contrast, chromophobe RCC and oncocytoma showed very similar levels of CA9 expression and were essentially indistinguishable.
|
|
|
PVALB expression was similar to CLCNKB in that the expression was highest in chromophobe (median Ct, 2.978) and lowest in clear cell type (median Ct, 17.787). Papillary RCC also showed low PVALB mRNA expression (median Ct, 15.596). Intriguingly, oncocytomas showed a broad range of PVALB expression ranging from 4.777 to 14.024, indicating that expression of this gene in oncocytoma is variable, in some cases much lower than in chromophobe RCC (range, 1.463-3.478). Because chromophobe RCC is always associated with strong PVALB expression, the finding of low PVALB expression is thus supportive of the diagnosis of oncocytoma, if chromophobe RCC is the only alternative diagnosis being considered.
AMACR expression was highest in papillary RCC (median Ct, 3.312). All other categories showed only minor differences in their AMACR levels, with the median for clear cell RCC, chromophobe RCC, and oncocytoma being 6.625, 6.697, and 8.706, respectively.
DEFB1 expression was evaluated in 25 of the 33 cases. It showed highest expression in chromophobe (median Ct, 4.907; range, 0.831 to 5.869) and oncocytoma (median Ct, 3.164; range, 1.251-10.868), with lower expression in papillary type (median Ct, 5.777; range, 3.565-11.549) and clear cell type (median Ct, 9.348; range, 5.293-15.465). The overlapping ranges in these four groups indicate that there is substantial variation in DEFB1 mRNA expression within individual types. Therefore, this gene was judged as having an expression profile similar to CLCNKB and PVALB but less consistent, and DEFB1 was not analyzed in subsequent experiments.
Gene expression ratios highlight differences among tumor types. From the distribution of the Ct values in Fig. 3, it is clear that a very high-level mRNA expression of CA9, AMACR, or CLCNKB/PVALB is essentially diagnostic of clear cell, papillary, and chromophobe RCC/oncocytoma, respectively. The inverse expression profile seen in some of these genes in different tumor types further suggested that evaluation of the expression ratios between paired genes could be more powerful in distinguishing these tumor types. For instance, because clear cell carcinomas express higher CA9 than papillary RCC and because papillary RCC express higher AMACR than clear cell RCC, the expression difference of these two genes, expressed as
Ct (CA9 AMACR; see Materials and Methods), could be anticipated to be more reliable than either gene alone in separating these two tumor types. Similarly,
Ct (CLCNKB AMACR) should help distinguish papillary RCC from chromophobe RCC/oncocytoma. This notion was proven correct, and Fig. 4 showed that
Ct (CA9 AMACR) could distinguish clear cell RC from all three other categories, whereas
Ct (CLCNKB AMACR) could distinguish papillary RCC from chromophobe RCC and oncocytoma.
|
In comparison to PCR data using fresh tissues, RNA from formalin-fixed, paraffin-embedded amplified poorer, with Ct values obtained
4 to 8 cycles behind what was observed in fresh tissues. However, the effects of this poor RNA quality were eliminated by the process of normalization with GAPDH endogenous control, and the values of the normalized Ct were in general similar to those described above for fresh tissues. However, for low abundant mRNA species, such as CA9 in chromophobe/oncocytoma group, no amplification was seen at the end of 45 cycles and an accurate numerical Ct could not be assigned. For such cases, a Ct value of 45 was used for the calculation of gene ratios (
Ct).
The qRT-PCR data on all 63 cases (excluding the 3 unclassified RCC) were summarized in Table 4, and the critical gene ratios were plotted in Fig. 5. Comparing Tables 3 and 4, it is clear that the gene expression profiles were highly similar between fresh-frozen and formalin-fixed tissues. However, two notable differences were observed. One was that the ranges of mRNA level distribution tend to be wider in the formalin-fixed, paraffin-embedded group for the same tumor category, likely attributed from the greater variation in the RNA quality of formalin-fixed tissues. The second difference observed was the PVALB expression in the chromophobe/oncocytoma group. In the fresh sample group, the chromophobe RCC showed higher level of PVALB expression than oncocytomas, and a Ct cutoff at 4.20 (Figs. 3 and 4), for instance, would clearly separate cases in these two categories. In the formalin-fixed group, all cases of the chromophobe RCC showed strong PVALB expression at levels similar to those seen in the fresh samples. However, the oncocytomas showed a broader range of PVALB expression in this group of 13 cases (median 3.302, range 2.814 to 11.018), significantly overlapping with the range seen in chromophobe RCC (2.815 to 3.673).
|
|
Ct (CA9 AMACR),
Ct (CLCNKB AMACR), and Ct (PVALB) were selected as predictors, and a diagnostic algorithm was generated using regression tree method and data from all 94 cases (31 fresh and 63 formalin-fixed). The algorithm was first built to classify the cases into clear cell, papillary, and chromophobe/oncocytoma groups based on
Ct (CA9 AMACR) and
Ct (CLCNKB AMACR). An attempt was then made to separate chromophobe RCC and oncocytoma based on the PVALB expression. The resulting regression tree and the decision rules were shown in Fig. 6.
|
For the 31 fresh-frozen cases, the chromophobe RCC and oncocytoma could be accurately classified using a Ct (PVALB) cutoff of 4.20. Among the 13 oncocytomas in the formalin-fixed group, however, only 6 cases would have been correctly classified. The other 7 cases showed higher level of PVALB expression, indistinguishable from chromophobe RCC. Therefore, we conclude that only oncocytomas with low PVALB expression (Ct > 4.20) can be reliably classified as oncocytoma, whereas cases with high PVALB expression can be either chromophobe RCC or oncocytoma.
This algorithm was then used to evaluate the five unclassifiable RCC in our series, two in the fresh-frozen group and three in the formalin-fixed group. All five cases were found to fall into the chromophobe/oncocytoma category. Their PVALB Ct values were 6.757, 12.915, 1.604, 3.673, and 9.604, respectively. These values predicted that three of the five cases were oncocytoma, whereas two cases would be classified as chromophobe RCC versus oncocytoma, not further classifiable. These predictions were fully compatible with the histologic features of these cases.
| Discussion |
|---|
|
|
|---|
From these studies, clear cell carcinoma, papillary carcinoma, and chromophobe carcinoma were found to have discrete microarray patterns. On the other hand, chromophobe RCC has clustered with oncocytoma and cannot be distinguished (11, 13, 14). Moreover, the technical complexity and the cost of DNA microarray prevent its use as a routine diagnostic technique at present, and alternative methodologies to translate the microarray-derived knowledge to daily practice are needed. One strategy is to identify or produce antibodies for the gene product(s) of interest, thus converting the DNA/RNAbased analysis to an immunohistochemistry-based assay that can be routinely done. This was indeed attempted in RCC, and antibodies against DEFB1, PVALB, and vimentin were tested, showing differential staining patterns in different types of RCC (10). A second approach, shown in the present study, is to evaluate the mRNA expression of a few selected genes in tumors by qRT-PCR. This qRT-PCRbased approach has several advantages over the antibody approach. One major advantage is that no antibody is needed. As many of the target genes identified through microarray encoded proteins that had no commercially available antibodies, production of new antibodies would be required. These new antibodies will need to have high sensitivity and specificity on paraffin sections to be useful diagnostic reagents, and this goal is not easily achievable for most new antibodies. Second, unlike qRT-PCR, antibody staining is at best semiquantitative, likely unreliable to reflect the quantitative differences detected by the microarray. The expression level of AMACR, for example, is consistently higher in papillary RCC than in others, but only by several-fold by qRT-PCR. Not unexpectedly, a proportion of clear cell RCC showed positive staining with anti-racemase antibody P504s, although papillary RCC did show preferential staining (24).3 The last but probably the most important advantage of the qRT-PCR strategy is its capability to quantitatively evaluate the expression of more than one gene in a cumulative fashion as illustrated by the gene ratio approach in our study. As expression of multiple genes are often coordinated, the ability of qRT-PCR to evaluate the cumulative effects of multiple gene expression and compare them between specimens is a distinct advantage. This concept was the basis of the recently developed assay for predicting prognosis in breast cancer (25, 26) and in lymphoma (27) in which a panel of genes were analyzed by qRT-PCR, resulting in a calculated index that was used to predict the prognosis and/or likelihood of recurrence.
This qRT-PCRbased approach certainly has its disadvantages in comparison with immunohistochemical assays. It is still technically more demanding and cannot be readily done in current pathology laboratories. However, this will likely change in the near future. For example, user-friendly machines using qRT-PCR technique are being developed for the detection of cytokeratin-positive cells in the sentinel lymph node biopsy specimens of breast cancer (28) and these machines can be easily modified to accommodate the analysis described in our study. Another major disadvantage is the contamination of normal tissues in the RNA preparations. Because tissue microdissection is unlikely a feasible approach for a routine diagnostic test, we had taken two measures to circumvent this drawback: one was to obtain sections on which adjacent nonneoplastic tissue is either absent or could be easily removed manually and the other was to carefully choose the gene targets, limiting the analysis to genes that are expressed exclusively or predominantly in tumor tissues but not in stromal and/or adjacent normal tissues. The choice of CA9, AMACR, CLCNKB, and PVALB reflected such considerations, and the nonneoplastic cells did not seem to adversely affect the data obtained.
CA9 is only expressed in pancreaticobiliary epithelium normally, and no expression was found in normal kidney (29). In contrast, in malignant neoplasm, CA9 is a marker of tumor hypoxia (30), found to be expressed in carcinomas of the kidney, lung, colon, and head and neck squamous cell carcinoma (3133). Among renal cancer, significant expression was found in the clear cell type (34, 35), and immunotherapeutic monoclonal antibody trials targeting CA9 are currently ongoing (36). Two antibodies against CA9 have been described [i.e., G250 (36) and M75 (37)]. Only M75 works on paraffin sections, but this antibody is not commercially available. Our current study confirmed high level CA9 expression in clear cell carcinoma, and a good commercial anti-CA9 antibody would be highly useful.
AMACR is a racemase that is normally expressed in several normal tissues, including liver, kidney, and salivary gland (38). Using anti-AMACR antibody P504s, renal tubular epithelial cells were found to express AMACR (38). For malignancies, AMACR is best recognized for its expression in prostate cancer (39). RCC and hepatocellular carcinoma have also been found to express AMACR (38). In kidney tumors, AMACR protein expression is most consistently found in papillary carcinoma but has also been observed in some cases of clear cell carcinoma (24), correlating with our observation of AMACR mRNA expression in these tumors.
PVALB, a high-affinity calcium binding protein, is expressed in several tissues, notably brain and kidney (40, 41). In kidney, the expression is limited to the distal nephron, particularly in the early distal convoluted tubules (42). Reflecting this tissue origin, PVALB mRNA expression has been found in chromophobe RCC and oncocytomas (10, 34, 43). It is of interest that we found chromophobe RCC to be more consistent in this expression, whereas some oncocytomas showed lower levels of PVALB. This quantitative difference may not be detectable by antibody staining. Rather, as suggested by Young et al. (10), positive PVALB immunostaining maybe a marker for chromophobe RCC and oncocytoma collectively, distinguishing them from clear cell and papillary RCC.
CLCNKB encodes one of two kidney-specific chloride channels (CLCNKA and CLCNKB), with expression localized to glomeruli, proximal tubules, and collecting ducts (44). Microarray data have suggested higher expression of this gene in chromophobe RCC/oncocytoma than in clear cell type, and this was confirmed by the present study. Given its preferential expression in selective kidney tumors, production of antibodies against this protein would be worthy of pursuing.
Based on the expression of these four genes, we were able to establish a molecular diagnostic model that is highly accurate in classifying renal cancer. As can be expected, chromophobe RCC and oncocytoma remain more difficult to separate. The ability of PVALB expression to predict at least some of the oncocytomas suggests the potential of this approach, and microarray studies to identify additional distinguishing markers between these two tumor types are ongoing.
| Footnotes |
|---|
3 Y.T. Chen, J.J. Tu, S.K. Tickoo, unpublished data. ![]()
Received 3/22/05; revised 6/ 9/05; accepted 6/21/05.
| References |
|---|
|
|
|---|
-methylacyl-CoA racemase (P504s) in various malignant neoplasms and normal tissues: a study of 761 cases. Hum Pathol 2003;34:7926.[Medline]
-methylacyl CoA racemase (P504S). Am J Clin Pathol 2004;122:27589.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
S. Rohan, J. J. Tu, J. Kao, P. Mukherjee, F. Campagne, X. K. Zhou, E. Hyjek, M. A. Alonso, and Y.-T. Chen Gene Expression Profiling Separates Chromophobe Renal Cell Carcinoma from Oncocytoma and Identifies Vesicular Transport and Cell Junction Proteins as Differentially Expressed Genes Clin. Cancer Res., December 1, 2006; 12(23): 6937 - 6945. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. J. Camp, J. M. Farnham, and L. A. Cannon-Albright Localization of a Prostate Cancer Predisposition Gene to an 880-kb Region on Chromosome 22q12.3 in Utah High-Risk Pedigrees. Cancer Res., October 15, 2006; 66(20): 10205 - 10212. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 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 |