
[Cancer Research 65, 2592-2601, April 1, 2005]
© 2005 American Association for Cancer Research
Molecular Biology, Pathobiology, and Genetics |
A Gene Expression Signature for Relapse of Primary Wilms Tumors
Wenliang Li1,2,
Patricia Kessler2,
Herman Yeger3,
Jennifer Alami3,
Anthony E. Reeve4,
Rosemary Heathcott4,
Jane Skeen5 and
Bryan R.G. Williams1,2
1 Department of Genetics, Case Western Reserve University; 2 Department of Cancer Biology, Cleveland Clinic Foundation, Cleveland, Ohio; 3 Department of Pediatric Laboratory Medicine, Hospital for Sick Children, Toronto, Canada; 4 Department of Biochemistry, University of Otago, Dunedin, New Zealand; and 5 Pediatric Haematology/Oncology, Starship Children's Hospital, Auckland, New Zealand
Requests for reprints: Bryan R.G. Williams, Department of Cancer Biology, Lerner Research Institute, Cleveland Clinic Foundation, 9500 Euclid Avenue, Cleveland, OH 44195. Phone: 216-445-9652; Fax: 216-444-3164; E-mail: williab{at}ccf.org.
 |
Abstract
|
|---|
Anaplastic histology and metastasis are each associated with higher relapse and mortality rates in Wilms tumor patients. However, not all anaplastic tumors relapse and some nonanaplastic tumors relapse unexpectedly. To identify more accurate early prognostic indicators, we analyzed expression of 4,900 cancer-related genes in 26 primary Wilms tumors. This analysis revealed that expression of a set of four genes predicts future relapse of primary Wilms tumors with high accuracy, independent of anaplasia. Random permutation testing of this prognostic gene expression signature yielded P = 0.003. Real-time reverse transcription-PCR analysis of the four genes in an independent primary tumor set resulted in correct prediction of future relapse with an accuracy of 92%. One of the four genes in the prognostic signature, CCAAT/enhancer binding protein ß (C/EBPB), is expressed at higher levels in both primary relapsing tumors and metastatic tumors than in primary nonrelapsing tumors. Short interfering RNAmediated down-regulation of C/EBPB expression in WiT49, a cell line derived from a metastatic Wilms tumor, resulted in spontaneous apoptosis. These findings suggest that C/EBPB is a critical survival factor for Wilms tumor cells and that its expression contributes to the prognosis of Wilms tumor patients.
Key Words: Wilms tumor profiling molecular signature relapse C/EBP ß
 |
Introduction
|
|---|
Wilms tumor (nephroblastoma) is the most common malignant abdominal tumor in children (1). Whereas 80% to 85% of Wilms tumor patients survive the disease, there remain 15% to 20% that experience relapse and low survival rates despite aggressive retreatment (2). On the other hand, those patients cured of their disease frequently experience adverse side effects from therapy (3). Therefore, it is important to identify patients at the time of diagnosis who might benefit from intensified initial treatment to prevent future relapse as well as those who might be cured with limited treatment and reduced toxicity. Anaplasia (unfavorable histology) and metastasis (stage) are each associated with higher relapse and mortality rates in Wilms tumor patients (4, 5). Nevertheless, not all anaplastic tumors relapse and some primary tumors with favorable histology relapse unexpectedly (6, 7). In addition, cytogenetic studies have suggested the existence of a genetically distinct subgroup of Wilms tumors that is resistant to chemotherapy, regardless of histology (8). Because neither histology nor stage accurately predicts the course of many Wilms tumors, better prognostic markers are required.
Potential prognostic markers for Wilms tumors have been indicated by studies demonstrating association of loss of heterozygosity of 16q, 11q, 1p, and possibly 22q with adverse outcome (912). In addition, for patients with favorable histology tumors, decreased relapse-free survival times and overall survival rates have been linked to high telomerase activity, gain of chromosome 1q or overexpression of genes on 1q (1315). Notably, strong p53 immunostaining, high levels of full-length TrkB mRNA or low levels of functionally inactive truncated TrkB mRNA have been suggested as negative prognostic factors independent of histology (7, 16).
In this study, we analyzed the expression of 4,900 cancer-related genes in 26 primary Wilms tumors and identified a gene expression pattern that is correlated with future relapse and independent of anaplastic histology. Moreover, one of the prognostic genes, CCAAT/enhancer binding protein ß (C/EBPB), was found to be overexpressed in both primary relapsing and metastatic tumors compared with primary nonrelapsing tumors and to control cell survival in a Wilms tumor cell line model.
 |
Materials and Methods
|
|---|
Cancer cDNA array construction. The cancer cDNA array used in this study contained probes for 9,240 cDNA clones corresponding to 4,900 unigenes or unclustered expressed sequence tags (EST) relevant to cancer or kidney development. These include (a) 1,800 genes implicated in metastasis or cancer development in general from the literature (17, 18) and Affymetrix cancer G110 array (http://www.affymetrix.com/analysis/index.affx); (b) 1,809 genes involved in kidney development (http://golgi.ana.ed.ac.uk/kidhome.html) or showing >3-fold variation of expression among four primary Wilms tumors and fetal kidney in our previous Affymetrix GeneChip experiments;6 (c) 1,400 potential downstream targets of p53 transcriptional regulation; (d) 950 genes (19) with AU-rich elements in their 3' and/or 5' untranslated region; (e) 1,245 cytokine responsive genes, genes for zinc finger proteins, or genes implicated in apoptosis. Construction of the array was similar to that described previously (19).
Sample collection, array experiment, and data processing. Fifty-four Wilms tumor specimens were obtained from the Cleveland Clinic Foundation (Cleveland, OH), the Hospital for Sick Children (Toronto, Canada), and the National Wilms Tumor Study Group (NWTS) under Institutional Review Board approval. All tumors are prechemotherapy. Patients were subjected to surgery first and treated similarly under the NWTS protocol. Among these tumors, 26 were primary tumors whose future relapse information was available, 12 from patients who went on to relapse within 4 years (termed "relapsing"), 14 from patients who remained relapse-free after their initial diagnosis for a period of at least 4 years (termed "nonrelapsing"). Pooled RNA from 13- to 24-week-old human fetal kidneys (Stratagene, La Jolla, CA) was used as a common reference in the cDNA array experiments. GeneSpring software (Silicongenetics, Redwood City, CA) was used for data analyses. Intensities below 100 were raised to 100 and nonlinear normalization (LOWESS) was done for each array to determine the abundance of each transcript in a tumor sample as a fold change with respect to its abundance in the human fetal kidney sample (common reference).
Relapse prediction. From the 26 primary tumors evaluable with respect to relapse, 18 tumor samples were randomly selected to form the training set, and the remaining eight tumors formed the testing set. Upon filtering out rarely expressed features (cDNAs) in the training set, signal-to-noise ratios, S2N = [Meanrelapsing Meannonrelapsing] / [SDrelapsing + SDnonrelapsing] (20) were calculated for remaining features (Supplementary Information). Welch t test was also used to evaluate the significance for each feature to be associated with relapse. These features were then rank ordered on the basis of their |S2N|. The top 1 feature (cDNA clone) was first selected and evaluated for its overall power to correctly classify the training samples using the "leave-one-out" cross-validation with the k-nearest neighbor (k-NN) algorithm (21). This performance evaluation procedure was repeated when one to two features were added each time from the top of the candidate list to optimize the number of genes in the molecular signature. An optimal set of genes obtained from the training samples through this cross-validation was evaluated by hierarchical clustering and then used to predict the testing samples.
To test the validity of the molecular signature, two other classification methods (discriminant analysis and classification tree model (S plus 6.0, http://www.insightful.com) were applied. To further assess the statistical significance of the molecular signature, 1,500 permutation of labels of the 26 samples were done, and for each permutated data set, a similar number of genes with the greatest |S2N| values were used in the discriminant analysis of the permuted data set. The cross-validation error rates and Hotelling's t test F statistic were recorded to generate the histogram plots and calculate the number of incidences and probability for a permutated data set to have model whose cross-validation error rate or Hotelling's t test F statistic matched or exceeded the molecular signature.
Real-time reverse transcription-PCR of CCL2 and vascular endothelial growth factor. cDNA reverse transcribed from 2.5 µg of total RNA using random primers was used as the template for amplification of CCL2 and vascular endothelial growth factor (VEGF) transcript sequences using SYBR Green Master Mix (Perkin-Elmer, Wellesley, MA) on an ABI 7700 system. Primers for CCL2 (forward primer 5'-TCGCCTCCAGCATGAAAGTC-3' and reverse primer 5'-ATGAAGGTGGCTGCTATGAGC-3') and VEGF (forward primer 5'-TACCTCCACCATGCCAAGTG-3' and reverse primer 5'-TCGTGATGATTCTGCCCTCC-3') were shown to give a single band of expected size on an agarose gel. Amplification of 18s rRNA using Taqman Master Mix (Perkin-Elmer) was used as a normalization control.
CT values = CTgene CT18s rRNA were calculated for CCL2 and VEGF for each primary relapsing and nonrelapsing tumor sample. The Welch t test was used to assess the significant level of
CT difference between primary relapsing and nonrelapsing tumors.
Immunohistochemical staining. Immunohistochemical staining was done on an independent set of samples using an indirect immunoperoxidase labeling method on paraffin sections as described previously (22). Rabbit anti-VEGF (A-20, 1:100) and mouse monoclonal anti-C/EBP ß (H7, 1:100) from Santa Cruz were used. H&E staining was used to verify histopathogical details and identify the most representative regions of each tumor.
Classification of tumors based on real-time reverse transcription-PCR of the four signature genes. Total RNA (1.5 µg) was reverse transcribed into cDNA using oligo-dT primer. The PCR primer sequences are as follows: C/EBP ß (forward primer 5'-GCCAAGAAGACCGTGGACA-3' and reverse primer 5'-GGCGATGTTGTTGCGCTC-3'), cDNA CF542255 (forward primer 5'-TTAGACTGGGATGAACACAGTTCTG-3' and reverse primer 5'-GAGCAATGTCTCAAAGGCTTACCT-3'), p21 (forward primer 5'-TGGAGACTCTCAGGGTCGAAA-3' and reverse primer 5'-CGGCGTTTGGAGTGGTAGAA-3'), H4FG (forward primer 5'-ACAAAACCGGCTATTCGCC-3' and reverse primer 5'-AAGACCGGAAATGCGCTTG-3'), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH forward primer 5'-TGGGCTACACTGAGCACCAG-3' and reverse primer 5'-CAGCGTCAAAGGTGGAGGAG-3'). Before real-time PCR analysis, each primer set was tested to ensure a single PCR product of expected size. Real-time PCR was done using SYBR Green Master mix on an ABI 7700 system. GAPDH expression was used as a normalization control.
CT values = CTgene CTGAPDH were calculated these genes in each primary relapsing and nonrelapsing tumor sample. Discriminant analysis was used to classify the tumor samples and assess the misclassification error rates (S plus 6.0).
Generation of short interfering RNAs. Three short interfering RNA (siRNA) sequences (21 nt) were selected from three positions of C/EBPB cDNA sequence and checked to ensure the lack of repetitive elements or homologies to other genes. The three siRNA target sequences and their starting positions on NM_005194 were siCB-1, AAATCCATGGAAGTGGCCAAC, 259 (5'); siCB-2, AAGTGGCCAACTTCTACTACG, 269 (5'); and siCB-3, AAACCAACCGCACATGCAGAT, 1417 (3' untranslated region). siRNAs were generated using the siRNA construction kit from Ambion according to the manufacture's instructions. Control siRNAs against GAPGH (siGAPDH) and SUMO-1 (siSUMO-1) were also generated. The siRNA target sequence for SUMO-1 (AAGGTGAATATATTAAACTCA) was from Akio Kishi et al. (23). The siRNA templates for GAPDH were provided in the construction kit.
Short interfering RNA transfection and immunoblotting. siCB-1, siCB-2, siCB-3, siSUMO-1, siGAPDH were transfected into WiT49 cells (30-40% confluent) at a final concentration of 25 nmol/L in Opti-MEM medium using oligofectAMINE as transfection reagent according to the manufacturer's protocol (Invitrogen, Carlsbad, CA). Mock transfection (oligofectAMINE reagent only) was done as a negative control. Forty-eight hours post-transfection, detached and adherent WiT49 cells were lysed as described (24). Extracts containing equal amounts of protein as determined by Bradford assay (Bio-Rad, Hercules, CA) were resolved by SDS-PAGE (10% gel) and subjected to immunoblotting. Antibodies used were H-7 (1:500, Santa Cruz Biotechnology, Santa Cruz, CA) for human C/EBPB, 9542 (1:1,000, Cell Signaling, Beverly, MA) for poly(ADP-ribose) polymerase (PARP), p56 (1:2,000, provided by Dr. Ganes Sen; ref. 25), 18-2306 (1:1,000, Zymed, South San Francisco, CA) for SUMO-1, 14AA (1:1,000, Santa Cruz Biotechnology) for C/EBPA, and A-5441 (1:5,000, Sigma, St. Louis, MO) for ß-actin.
 |
Results and Discussion
|
|---|
Prediction of Wilms tumor relapse based upon gene expression signature. We generated gene expression profiles for 26 primary Wilms tumors for which we had information regarding future relapse. Of these, 12 tumors were from patients who went on to relapse within 4 years (termed "relapsing"), seven of which had favorable histology. The remaining 14 were from patients who remained relapse-free after their initial diagnosis for a period of at least 4 years (termed "nonrelapsing"), all having favorable histology. Relapse of Wilms tumor typically occurs 2 to 3 years after original diagnosis. We did not observe a correlation between relapse and stage of these primary tumors (Fisher's exact test P = 0.93, Supplementary Table 1). To identify genes whose expression in the primary tumor is associated with its tendency towards future relapse, we randomly selected 18 tumors from the 26 primary tumors to form a training set. The remaining eight tumors constituted the test set. In the training set, we used both the signal-to-noise (S2N) ratio matrix {[S2N=(Meanrelapsing Meannonrelapsing)/(SDrelapsing + SDnonrelapsing)]; ref. 20} and the Welch t test to search for signature genes associated with relapse. Among the 7,462 clones that remained in the data set after removal of rarely expressed clones (Supplementary Information), 199 clones (182 unique genes or ESTs) were found to be strongly associated with relapse (absolute S2N ratio (|S2N|) > 0.45 and P < 0.04). Expression of multiple clones for the same genes/clusters correlated nicely (Supplementary Table 2). The correlation coefficient is R2 = 0.985. The 199 clones were ordered by decreasing |S2N| ratio. We then optimized the molecular signature by sequentially adding one to two genes from the rank-ordered list and evaluating their overall power for correct classification in the leave-one-out cross-validation with the k-nearest neighbor (k-NN) algorithm. We stopped this process when an optimal gene set consisting of a minimal number of genes yet having maximal accuracy in cross-validation of the training set was reached. The top four genes in the rank-ordered list (C/EBPB, cDNA CF542255, p21, and H4FG) were capable of distinguishing between the primary relapsing and nonrelapsing tumors in the training set with 100% accuracy. Classification of the 18 training samples by the four-gene signature was confirmed by unsupervised hierarchical clustering analysis, where the primary relapsing tumors were clearly separated from primary nonrelapsing tumors (Fig. 1A). The trained model, based on the expression patterns of C/EBPB, cDNA CF542255, p21, and H4FG in the 18 training samples, was then used to predict the relapse status of the eight tumors in the test set using the k-NN algorithm. All eight samples were correctly classified as relapsing or nonrelapsing and this again was confirmed by unsupervised hierarchical clustering (Fig. 1B). We tested the performance of the four-gene signature in the 18-sample training set and eight-sample test set with two other classification methods, discriminant analysis and classification tree model, both of which yielded similarly high accuracy (data not shown). To ensure that the predictive power of the four-gene signature was not biased by the particular composition of the training and testing sets, we randomly formed several other training sets (range, 13-19 samples) from the 26 samples with different relapsing/nonrelapsing ratios and used the expression patterns of the four genes in the training sets to predict the corresponding testing sets. The average accuracies in cross-validation of the training sets and prediction of the testing sets were 98.2% and 94.1%, respectively.

View larger version (38K):
[in this window]
[in a new window]
|
Figure 1. Hierarchical clustering of the 18 training samples (A) and eight testing samples (B) using the four genes. Sample (column) and gene (row). Color scale, expression level for a gene in a sample relative to normal fetal kidney. Sample trees (x-axis) and gene trees (y-axis). Bottom, sample types: relapsing (blue) or nonrelapsing (red). Signal-to-noise ratio S2N = (Meanrelapsing Meannonrelapsing) / (SDrelapsing + SDnonrelapsing).
|
|
To further assess the statistical significance of the four-gene prognostic model, we did 1,500 permutations of labels of the 26 samples, and for each permuted data set, similar numbers of genes with the greatest |S2N| values were used in discriminant analysis. Only 5 of the 1,500 permutations yielded models whose cross-validation error rates or Hotelling's t test F statistics matched or exceeded those of the four-gene model. Thus, the likelihood of matching the success of the observed four-gene model by chance was estimated at P = 0.003. The four-gene signature that we have identified seems to clearly distinguish relapsing tumors from nonrelapsing tumors and show predictive power.
In general, C/EBPB and H4FG (HIST1H4C) were up-regulated in relapsing tumors compared with nonrelapsing tumors, whereas p21 and cDNA CF542255 were down-regulated. It is important to note that it is the expression pattern of all four genes taken together that provides the predictive power described above. Expression of any single gene is not absolutely correlated with prognosis. The up-regulation of C/EBPB in tumors that are destined to relapse is intriguing given that C/EBPB is a critical component of a Ras-dependent transformation pathway (26). Moreover, transcriptional activity of C/EBPB leads to induction of a number of prometastatic genes, including hepatocyte growth factor HGF, IL-6, and IL-8, which in turn induce C/EBPB (27, 28). As further support for its potential role in tumor progression in general, increased expression of C/EBPB has been observed in solid cancers other than Wilms tumor, such as breast, ovarian, and colorectal cancer (29). In tumors that go on to relapse, the activities of C/EBPB are combined with an elevated proliferative rate as indicated by increased histone H4FG expression and decreased levels of the cyclin-dependent kinase inhibitor p21. The contribution of cDNA CF542255 to the progression of Wilms tumor is unclear since its protein product and biological function remain to be determined.
The complete list of 199 clones (182 genes) that are differentially expressed between primary relapsing and primary nonrelapsing Wilms tumors includes many additional genes implicated in processes relevant to tumor progression such as proliferation, metastasis, angiogenesis, apoptosis, transcriptional regulation, and signal transduction. Examples of these genes are those expressed higher in relapsing tumors: VEGF, STAT1, IL-8, integrin ß 4 binding protein, coagulation factor VII (F7), collagen type VIII
2, retinoblastoma binding protein 8, hepatocyte growth factor, early growth response 2, neuromedin B, Ras homologue enriched in brain (RHEB2), C-C motif chemokine receptor 1 (CCR1), and C-C motif chemokine ligand 2 (CCL2); and those expressed lower in relapsing tumors: IGF2R, ZNF23, ZNF33a, caspase 6 (CASP6), mothers against decapentaplegic homologue SMAD3, E2F transcription factor 4 (E2F4), transforming acidic coiled coil containing protein 1, active BCR-related gene, 90-kDa heat shock protein, BCL2-associated X protein (BAX), damage-specific DNA binding protein 2, and suppressor of cytokine signaling 5 (see Supplementary Table 2 for a complete list). Together with our expression data, previously published functional information for a number of these genes suggests potential mechanisms for Wilms tumor progression. For example, our data show that induction of a number of proangiogenic/metastatic factors is associated with Wilms tumor relapse. These factors include CCL2, CCR1, IL-8, VEGF, F7, HGF and of course, C/EBPB. Interestingly, we recently identified significant coexpression of HGF and its receptor, c-Met, in Wilms tumors (22). Furthermore, abnormal expression of adhesion proteins regulated by HGF/c-Met signaling was found to correlate with Wilms tumor invasiveness and metastasis (30). It is also likely that some of the relapse-associated gene products cooperate to promote angiogeneis/metastsis. C/EBPB and HGF, VEGF, and HGF induce each other, forming positive feedback loops. In addition, both F7 and HGF positively regulate the expression of VEGF and IL-8 (31, 32). Whereas proangiogenic/metastasis genes are up-regulated in relapsing Wilms tumors, a number of genes that encode proapoptotic or growth suppressive factors (such as CASP6, BAX, IGF2R, E2F4, and p21) are down-regulated. Studies examining gene expression in a number of other types of cancer have also identified many of the genes on our Wilms tumor relapse-associated list (29, 3335). Thus, our data add support to the idea that a common molecular mechanism underlies metastasis of many different human solid cancers (36).
Cytoplasmic immunopositivity for p53, which indicates a deficiency in its activities, has been suggested as a poor prognostic marker for Wilms tumors, independent of histology (7). It is possible that the gene expression differences between nonrelapsing and relapsing tumors are linked to inactivation of p53 pathway. This is supported by the down-regulation of p21 we detected in relapsing tumors. There are little overlap between the 199 clones we defined and a partial list of 80 cDNA clones differentially expressed between nine nonrelapsing tumors with favorable histology and four relapsing tumors with favorable histology described in a recent report (37). This is probably due to incompleteness of their gene list, different cDNA array platforms and their focus on tumors with favorable histology. Although insights into the aggressiveness of favorable histology tumors may be provided in their study, restricting gene expression profiling to favorable histology tumors can only identify prognostic markers that are supplementary to histology.
Confirmation of microarray gene expression patterns. The protein expression pattern for C/EBPB was examined by immunohistochemical staining on an independent set of samples. Four cases of primary relapsing, three cases of metastatic and four cases of primary nonrelapsing Wilms tumors were examined. C/EBPB protein staining was detected in primary relapsing (2 of 4) and metastatic tumors (2 of 3) but not in primary nonrelapsing tumors (0 of 4; Fig. 2A). Positive C/EBPB staining was nuclear, concurring with its function as a transcription factor. The expression levels of VEGF, another gene expressed higher in primary relapsing tumor than in nonrelapsing tumors, was also analyzed. Its staining was much stronger in primary relapsing tumors (4 of 4) and marginally stronger in metastatic cases (2 of 3) than in nonrelapsing cases (1 of 4). Positive VEGF staining is cytoplasmic and staining of blood vessels serves as an internal positive control. Thus, even in this relatively small sample set, the protein level differences of C/EBPB and VEGF in primary nonrelapsing, primary relapsing, and metastatic Wilms tumors are notable.

View larger version (63K):
[in this window]
[in a new window]
|
Figure 2. Confirmation of cDNA array results by immunohistochemistry and real-time PCR. A, representative immunohistochemical staining photomicrographs for C/EBPB and VEGF in independent primary nonrelapsing, primary relapsing, and metastatic Wilms tumor specimens. B and C, CCL2 and VEGF expression in primary relapsing and nonrelapsing Wilms tumors as assayed by real-time PCR. Relative expression levels (folds, y-axis) to one arbitrary chosen sample (T42) for CCL2 (B) and VEGF (C).
|
|
We also arbitrarily selected two genes down the list (in addition to the top four genes, see next) to validate the cDNA array expression patterns using real-time reverse transcription-PCR (RT-PCR). Similar to VEGF, CCL2 is also expressed at a higher levels in primary relapsing tumors than in primary nonrelapsing tumors based on the cDNA array data. When analyzed by real-time PCR, expression of CCL2 and VEGF were significantly higher in primary relapsing tumors (n = 12) than in primary nonrelapsing tumors (n = 11; Welch's t test P = 0.0067 and 4.26-fold for CCL2, and P = 0.021, 2.95-fold for VEGF; Fig. 2B and C). Both the immunohistochemical staining and real-time RT-PCR results were in agreement with the cDNA array data.
Application of the four-gene prognostic signature to an independent sample set using real-time RT-PCR. To determine whether the four-gene prognostic model might be of practical use in stratification of patients, we designed real-time RT-PCR assays to quantify expression of the four genes. Gene expression as determined by these assays was used to classify an independent set of 12 primary Wilms tumors (6 relapsing and 6 nonrelapsing). Eleven of the 12 samples (92%) were correctly classified with respect to future relapse using the RT-PCR results in either discriminant analysis or tree model classification. When the RT-PCR results were visualized in a three-dimensional scatter plot using the three most significant genes from cDNA array data (C/EBPB, cDNA CF 542255, and p21), six primary nonrelapsing tumors were separated from five primary relapsing tumors, with one primary relapsing tumor behaving as outlier (Fig. 3). Importantly, among these 12 tumors, only one was known to have unfavorable histology. Thus, this sample set provides further support of the improved accuracy of the four-gene signature as a prognosticator compared with histology. These results not only confirm our cDNA array data on an independent sample set using a different approach but also indicate that the four-gene prognostic model may be applied in practice using a simpler technique like RT-PCR.
Identification of genes associated with Wilms tumor metastasis. In addition to identifying genes associated with Wilms tumor relapse, we also mined our array data for genes that are differentially expressed in both metastatic tumors and primary tumors that display future metastasis compared with primary tumors that do not metastasize. Among the 12 primary relapsing tumors in our original sample set, one later relapsed only in the original tumor bed, whereas 11 later relapsed outside of the primary site, mostly in the lung or abdomen. Relapse of primary Wilms tumors typically occurs in the lung or abdomen (38). Because several studies (including this work) have shown that the metastatic potential of human cancers is embodied in the gene expression profiles of the primary tumors (3941), it is presumed that some relevant gene expression changes necessary to initiate the metastatic process occur in primary metastasizing tumors. Some of these gene expression changes might be preserved in the final metastasis, which would suggest roles in maintenance of the phenotype as well. Thus, comparison of both primary metastasizing tumors and metastatic tumors to primary nonmetastasizing tumors may identify gene expression changes involved in both initiation and maintenance of the metastatic phenotype. Accordingly, we identified 95 clones that were persistently up-regulated or down-regulated in both the 11 primary metastasizing tumors and six metastatic (relapsed) tumors compared with the 14 primary nonmetastasizing (nonrelapsing) tumors (|S2N| > 0.4 and Welch t test two-sided P < 0.05).
A subset of the 95 clones are shown in Table 1, including clones whose expression changes at least 50% with a Welch t test two-sided P < 0.05 in both the metastasizing/nonmetastasizing and metastatic/nonmetastasizing comparisons (see Supplementary Table 3 for a complete list). We found that C/EBPB, one of our relapse signature genes, is also among the genes up-regulated in both primary metastasizing tumors and metastatic tumors compared with primary nonmetastasizing tumors. This finding is consistent with other studies implicating C/EBPB in tumor metastasis as described above. Increased expression of C/EBPB in metastatic Wilms tumors was confirmed at the protein level by immunohistochemical staining (Fig. 2A).
View this table:
[in this window]
[in a new window]
|
Table 1. Genes persistently up-regulated or down-regulated in both primary metastasizing and metastatic tumors compared with primary nonmetastasizing tumors
|
|
To discover enriched biological themes within the list of 95 clones associated with Wilms tumor metastasis, we used the EASE program (http://david.niaid.nih.gov/david/ease.htm; 42) and noticed that five general categories were significantly overrepresented: calcineurin signaling, cell cycle regulation, protein/nucleus import, metabolism, and apoptosis regulation (Supplementary Table 4). It was noted that biological themes relating to calcineurin signaling ranked highest among the descriptive categories. Three calcineurin isoforms, including calcineurin A
, Aß, and B type I, were all up-regulated in both primary metastasizing and metastatic tumors compared with primary nonmetastasizing tumors (Table 1; Supplementary Table 3). Calcineurin activates NFAT (nuclear factor of activated T cells), which mediates VEGF signaling and results in transcriptional activation of proangiogenic genes such as tissue factor and cyclooxygenase 2 (43, 44). This result suggests that calcineurin-NFAT signaling plays an important roles in Wilms tumor metastasis.
C/EBPB is required for survival of Wilms tumor cells. C/EBPB is an essential member of the molecular signature for relapse of primary Wilms tumor that we identified through microarray expression analysis. Three distinct polypeptides are produced from the C/EBPB gene through use of different translation initiation sites. The LAP-1 (46 kDa) and LAP-2 (42 kDa) isoforms differ at the NH2 terminus, whereas translation from a third downstream in-frame AUG results in a 20-kDa LIP isoform (ref. 45; Fig. 4A). Previous studies suggest that C/EBPB may be involved in tumor development, with elevated levels of C/EBPB being found in human breast (reviewed in ref. 45), ovarian (46), colorectal (29), renal cell (47), and other cancers. Furthermore, in a mouse model of carcinogen-induced skin tumorigenesis involving oncogenic Ras signaling, C/EBPB-null mice were resistant to tumor development (26). Here the role of C/EBPB in Wilms tumor cell growth and survival was investigated by siRNA-mediated suppression of C/EBPB expression in WiT49, a cell line derived from a metastatic Wilms tumor (48). Three different siRNA target sequences specific for the C/EBPB transcript were designed (Fig. 4A). As controls, siRNAs against GAPDH and SUMO-1 were also generated.

View larger version (49K):
[in this window]
[in a new window]
|
Figure 4. A, schematics of C/EBPB cDNA sequence, three protein isoforms, and siRNA target positions. B, effects of siRNAs on the protein expression of C/EBPBA and the three protein isoforms of C/EBPB. Immunoblotting was performed 48 hours after siRNA transfection.
|
|
As measured by immunoblotting, all three isoforms of C/EBPB protein were significantly down-regulated following transfection of WiT49 cells with C/EBPB siRNAs, with siCB-2 being the most potent (Fig. 4B). Notably, levels of the closely related protein CCAAT/enhancer binding protein
(C/EBPA) were not reduced by siRNAs directed against C/EBPB. Rather, C/EBPA protein levels increased as C/EBPB protein expression was down-regulated by C/EBPB-specific siRNAs. This observation is in agreement with some other circumstances where C/EBPA expression decreases when C/EBPB expression increases or vice versa (4951). Transfection of siRNA directed against the SUMO-1 transcript resulted in reduced levels of SUMO-1 protein, but not C/EBPB protein (Fig. 5). These observations show the specificity of gene targeting by the different siRNAs.

View larger version (48K):
[in this window]
[in a new window]
|
Figure 5. Down-regulation of C/EBPB protein by siRNAs correlates with an increase in PARP cleavage in Wit49 cells. Immunoblotting was done with antibodies against human C/EBPB (H-7, Santa Cruz Biotechnology), PARP (9542, Cell Signaling), p56 (25), SUMO-1 (18-2306, Zymed), and ß-actin (A-5441, Sigma).
|
|
It has been reported by our lab and others that siRNAs can activate double-stranded RNA (dsRNA) responsive pathways in cells (5254). Indeed, as shown in Fig. 5, transfection of the different siRNAs targeting C/EBPB, GAPDH, and SUMO-1 led to induction of the classic dsRNA-responsive gene product p56 (25). Moreover, the GAPDH and SUMO-1 siRNAs induced C/EBPB transcript and protein expression (Figs. 4 and 5), suggesting that C/EBPB is also a dsRNA-responsive gene. Importantly, an increase in spontaneous apoptosis was observed in C/EBPB siRNA-treated cells as shown by proteolytic cleavage of PARP. The amount of PARP cleavage correlated with the extent of down-regulation of C/EBPB proteins (Fig. 5). However, unlike those targeting C/EBPB, siRNAs directed against the GAPDH and SUMO-1 transcripts did not induce PARP cleavage in WiT49 cells (Fig. 5). These results indicate that down-regulation of C/EBPB specifically leads to spontaneous apoptosis of the tumor cells, suggesting that continued expression of C/EBPB is critical survival factor for Wilms tumor cells.
In summary, we have identified a four-gene molecular signature that predicts future relapse of primary Wilms tumors with high accuracy. Real-time RT-PCR of the four genes showed that this prognostic model, if replicated in larger tumor sets, has the potential for clinical application in the stratification of Wilms tumor patients at diagnosis. It is important to note that lacking samples with anaplasia that are nonrelapsing in our array analyses may limit the generalizability of this molecular signature, although one of the tumors stained by immunohistochemistry was anaplastic and nonrelapsing and it showed no expression of C/EBPB proteins similar to other nonrelapsing tumors with favorable histology (Supplementary Table 1). Furthermore, our observations that C/EBPB contributes to both patient prognosis and survival of cultured Wilms tumor cells support the hypothesis that C/EBPB is involved in the development and/or progression of Wilms tumor. Taken together with previous studies implicating C/EBPB in cell growth and tumorigenesis, our findings suggest that C/EBPB may be a physiologically relevant therapeutic target for cancer treatment.
 |
Acknowledgments
|
|---|
Grant support: NIH grant R01-CA89279 (B.R.G. Williams), Health Research Council of New Zealand (A.E. Reeve), Cancer Society of NZ (A.E. Reeve), and Lottery Grants Board (A.E. Reeve).
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.
We thank the National Wilms Tumor Study Group for providing Wilms tumor specimens; Drs. Alexandru Almasan, Patricia Stanhope-Baker, Serpil Erzurum, and Shaoqi Rao for critical comments on the article; and Dr. Anthony Sadler for a fruitful discussion.
 |
Footnotes
|
|---|
Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org).
6 Unpublished data. 
Received 4/30/04.
Revised 1/19/05.
Accepted 1/28/05.
 |
References
|
|---|
- Breslow N, Olshan A, Beckwith JB, Green DM. Epidemiology of Wilms tumor. Med Pediatr Oncol 1993;21:17281.[Medline]
- Jones KP, Grundy PE, Coppes MJ. Recent advances in the genetics of childhood renal cancers: a report of the 3rd International Conference on the molecular and clinical genetics of childhood renal tumors, together with the Mitchell Ross symposium on anaplastic and other high risk embryonal tumors of childhood, 8-10th April 1999, Wistar Institute, Philadelphia, PA. Med Pediatr Oncol 2000;35:12630.[CrossRef][Medline]
- Egeler RM, Wolff JE, Anderson RA, Coppes MJ. Long-term complications and post-treatment follow-up of patients with Wilms' tumor. Semin Urol Oncol 1999;17:5561.[Medline]
- Bonadio JF, Storer B, Norkool P, Farewell VT, Beckwith JB, D'Angio GJ. Anaplastic Wilms' tumor: clinical and pathologic studies. J Clin Oncol 1985;3:51320.[Abstract]
- Beckwith JB, Palmer NF. Histopathology and prognosis of Wilms tumors: results from the First National Wilms' Tumor Study. Cancer 1978;41:193748.[CrossRef][Medline]
- Gow KW, Murphy JJ. Cytogenetic and histologic findings in Wilms' tumor. J Pediatr Surg 2002;37:8237.[CrossRef][Medline]
- Lahoti C, Thorner P, Malkin D, Yeger H. Immunohistochemical detection of p53 in Wilms' tumors correlates with unfavorable outcome. Am J Pathol 1996;148:157789.[Abstract]
- Stock C, Ambros IM, Lion T, et al. Genetic changes of two Wilms tumors with anaplasia and a review of the literature suggesting a marker profile for therapy resistance. Cancer Genet Cytogenet 2002;135:12838.[CrossRef][Medline]
- Grundy P, Telzerow P, Moksness J, Breslow NE. Clinicopathologic correlates of loss of heterozygosity in Wilm's tumor: a preliminary analysis. Med Pediatr Oncol 1996;27:42933.[CrossRef][Medline]
- Grundy PE, Telzerow PE, Breslow N, Moksness J, Huff V, Paterson MC. Loss of heterozygosity for chromosomes 16q and 1p in Wilms' tumors predicts an adverse outcome. Cancer Res 1994;54:23313.[Abstract/Free Full Text]
- Klamt B, Schulze M, Thate C, et al. Allele loss in Wilms tumors of chromosome arms 11q, 16q, and 22q correlate with clinicopathological parameters. Genes Chromosomes Cancer 1998;22:28794.[CrossRef][Medline]
- Grundy RG, Pritchard J, Scambler P, Cowell JK. Loss of heterozygosity on chromosome 16 in sporadic Wilms' tumour. Br J Cancer 1998;78:11817.[Medline]
- Dome JS, Chung S, Bergemann T, et al. High telomerase reverse transcriptase (hTERT) messenger RNA level correlates with tumor recurrence in patients with favorable histology Wilms' tumor. Cancer Res 1999;59:43017.[Abstract/Free Full Text]
- Hing S, Lu YJ, Summersgill B, et al. Gain of 1q is associated with adverse outcome in favorable histology Wilms' tumors. Am J Pathol 2001;158:3938.[Abstract/Free Full Text]
- Lu YJ, Hing S, Williams R, Pinkerton R, Shipley J, Pritchard-Jones K. Chromosome 1q expression profiling and relapse in Wilms' tumour. Lancet 2002;360:3856.[CrossRef][Medline]
- Eggert A, Grotzer MA, Ikegaki N, et al. Expression of the neurotrophin receptor TrkB is associated with unfavorable outcome in Wilms' tumor. J Clin Oncol 2001;19:68996.[Abstract/Free Full Text]
- Hippo Y, Yashiro M, Ishii M, et al. Differential gene expression profiles of scirrhous gastric cancer cells with high metastatic potential to peritoneum or lymph nodes. Cancer Res 2001;61:88995.[Abstract/Free Full Text]
- Scholl FA, Betts DR, Niggli FK, Schafer BW. Molecular features of a human rhabdomyosarcoma cell line with spontaneous metastatic progression. Br J Cancer 2000;82:123945.[CrossRef][Medline]
- Frevel MA, Bakheet T, Silva AM, Hissong JG, Khabar KS, Williams BR. p38 mitogen-activated protein kinase-dependent and -independent signaling of mRNA stability of AU-rich element-containing transcripts. Mol Cell Biol 2003;23:42536.[Abstract/Free Full Text]
- Golub TR, Slonim DK, Tamayo P, et al. Molecular classification of cancer: class discovery and class pre-diction by gene expression monitoring. Science 1999;286:5317.[Abstract/Free Full Text]
- Nutt CL, Mani DR, Betensky RA, et al. Gene expression-based classification of malignant gliomas correlates better with survival than histological classification. Cancer Res 2003;63:16027.[Abstract/Free Full Text]
- Alami J, Williams BR, Yeger H. Expression and localization of HGF and met in Wilms' tumours. J Pathol 2002;196:7684.[CrossRef][Medline]
- Kishi A, Nakamura T, Nishio Y, Maegawa H, Kashiwagi A. Sumoylation of Pdx1 is associated with its nuclear localization and insulin gene activation. Am J Physiol Endocrinol Metab 2003;284:E83040.[Abstract/Free Full Text]
- Goh KC, Haque SJ, Williams BR. p38 MAP kinase is required for STAT1 serine phosphorylation and transcriptional activation induced by interferons. EMBO J 1999;18:56018.[CrossRef][Medline]
- Guo J, Peters KL, Sen GC. Induction of the human protein P56 by interferon, double-stranded RNA, or virus infection. Virology 2000;267:20919.[CrossRef][Medline]
- Zhu S, Yoon K, Sterneck E, Johnson PF, Smart RC. CCAAT/enhancer binding protein-ß is a mediator of keratinocyte survival and skin tumorigenesis involving oncogenic Ras signaling. Proc Natl Acad Sci U S A 2002;99:20712.[Abstract/Free Full Text]
- Cho MK, Kim SG. Hepatocyte growth factor activates CCAAT enhancer binding protein and cell replication via PI3-kinase pathway. Hepatology 2003;37:68695.[CrossRef][Medline]
- Jiang JG, Zarnegar R. A novel transcriptional regulatory region within the core promoter of the hepatocyte growth factor gene is responsible for its inducibility by cytokines via the C/EBP family of transcription factors. Mol Cell Biol 1997;17:575870.[Abstract]
- Rask K, Thorn M, Ponten F, et al. Increased expression of the transcription factors CCAAT-enhancer binding protein-ß (C/EBß) and C/EBzeta (CHOP) correlate with invasiveness of human colorectal cancer. Int J Cancer 2000;86:33743.[CrossRef][Medline]
- Alami J, Williams BR, Yeger H. Differential expression of E-cadherin and ß catenin in primary and metastatic Wilms's tumours. Mol Pathol 2003;56:21825.[Abstract/Free Full Text]
- Dong G, Chen Z, Li ZY, Yeh NT, Bancroft CC, Van Waes C. Hepatocyte growth factor/scatter factor-induced activation of MEK and PI3K signal pathways contributes to expression of proangiogenic cytokines interleukin-8 and vascular endothelial growth factor in head and neck squamous cell carcinoma. Cancer Res 2001;61:59118.[Abstract/Free Full Text]
- Wang X, Gjernes E, Prydz H. Factor VIIa induces tissue factor-dependent up-regulation of interleukin-8 in a human keratinocyte line. J Biol Chem 2002;277:236206.[Abstract/Free Full Text]
- Yao R, Wang Y, Lubet RA, You M. Differentially expressed genes associated with mouse lung tumor progression. Oncogene 2002;21:581421.[CrossRef][Medline]
- Vaarala MH, Porvari K, Kyllonen A, Vihko P. Differentially expressed genes in two LNCaP prostate cancer cell lines reflecting changes during prostate cancer progression. Lab Invest 2000;80:125968.[Medline]
- Ugurel S, Rappl G, Tilgen W, Reinhold U. Increased serum concentration of angiogenic factors in malignant melanoma patients correlates with tumor progression and survival. J Clin Oncol 2001;19:57783.[Abstract/Free Full Text]
- Ramaswamy S, Ross KN, Lander ES, Golub TR. A molecular signature of metastasis in primary solid tumors. Nat Genet 2003;33:4954.[CrossRef][Medline]
- Takahashi M, Yang XJ, Lavery TT, et al. Gene expression profiling of favorable histology Wilms tumors and its correlation with clinical features. Cancer Res 2002;62:6598605.[Abstract/Free Full Text]
- Miser JS, Tournade MF. The management of relapsed Wilms tumor. Hematol Oncol Clin North Am 1995;9:1287302.[Medline]
- Ramaswamy S, Ross KN, Lander ES, Golub TR. A molecular signature of metastasis in primary solid tumors. Nat Genet 2003;33:4954.
- van't Veer LJ, Dai H, van de Vijver MJ, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature 2002;415:5306.[CrossRef][Medline]
- van de Vijver MJ, He YD, van't Veer LJ, et al. A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med 2002;347:19992009.[Abstract/Free Full Text]
- Dennis G Jr, Sherman BT, Hosack DA, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol 2003;4:P3.[CrossRef][Medline]
- Armesilla AL, Lorenzo E, Gomez del Arco P, Martinez-Martinez S, Alfranca A, Redondo JM. Vascular endothelial growth factor activates nuclear factor of activated T cells in human endothelial cells: a role for tissue factor gene expression. Mol Cell Biol 1999;19:203243.[Abstract/Free Full Text]
- Hernandez GL, Volpert OV, Iniguez MA, et al. Selective inhibition of vascular endothelial growth factor-mediated angiogenesis by cyclosporin A: roles of the nuclear factor of activated T cells and cyclooxygenase 2. J Exp Med 2001;193:60720.[Abstract/Free Full Text]
- Zahnow CA. CCAAT/enhancer binding proteins in normal mammary development and breast cancer. Breast Cancer Res 2002;4:11321.[CrossRef][Medline]
- Sundfeldt K, Ivarsson K, Carlsson M, et al. The expression of CCAAT/enhancer binding protein (C/EBP) in the human ovary in vivo: specific increase in C/EBPß during epithelial tumour progression. Br J Cancer 1999;79:12408.[CrossRef][Medline]
- Oya M, Horiguchi A, Mizuno R, Marumo K, Murai M. Increased activation of CCAAT/enhancer binding protein-ß correlates with the invasiveness of renal cell carcinoma. Clin Cancer Res 2003;9:10217.[Abstract/Free Full Text]
- Alami J, Williams BR, Yeger H. Derivation and characterization of a Wilms' tumour cell line, WiT 49. Int J Cancer 2003;107:36574.[CrossRef][Medline]
- Gilpin DA, Hsieh CC, Kuninger DT, Herndon DN, Papaconstantinou J. Effect of thermal injury on the expression of transcription factors that regulate acute phase response genes: the response of C/EBP
, C/EBPß, and C/EBP
to thermal injury. Surgery 1996;119:67483.[CrossRef][Medline]
- Alam T, An MR, Papaconstantinou J. Differential expression of three C/EBP isoforms in multiple tissues during the acute phase response. J Biol Chem 1992;267:50214.[Abstract/Free Full Text]
- Skrtic S, Ekberg S, Wallenius V, Enerback S, Hedin L, Jansson JO. Changes in expression of CCAAT/enhancer binding protein
(C/EBP
) and C/EBPß in rat liver after partial hepatectomy but not after treatment with cyproterone acetate. J Hepatol 1997;27:90311.[CrossRef][Medline]
- Gitlin L, Karelsky S, Andino R. Short interfering RNA confers intracellular antiviral immunity in human cells. Nature 2002;418:4304.[CrossRef][Medline]
- Bridge AJ, Pebernard S, Ducraux A, Nicoulaz AL, Iggo R. Induction of an interferon response by RNAi vectors in mammalian cells. Nat Genet 2003;34:2634.[CrossRef][Medline]
- Sledz CA, Holko M, de Veer MJ, Silverman RH, Williams BR. Activation of the interferon system by short-interfering RNAs. Nat Cell Biol 2003;5:8349.[CrossRef][Medline]
This article has been cited by other articles:

|
 |

|
 |
 
P. Gade, S. K. Roy, H. Li, S. C. Nallar, and D. V. Kalvakolanu
Critical Role for Transcription Factor C/EBP-{beta} in Regulating the Expression of Death-Associated Protein Kinase 1
Mol. Cell. Biol.,
April 15, 2008;
28(8):
2528 - 2548.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
A. D Singh, K. Sisley, Y. Xu, J. Li, P. Faber, S. J Plummer, H. S Mudhar, I. G Rennie, P. M Kessler, G. Casey, et al.
Reduced expression of autotaxin predicts survival in uveal melanoma
Br. J. Ophthalmol.,
October 1, 2007;
91(10):
1385 - 1392.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. Delettre, V. J. Yuste, R. S. Moubarak, M. Bras, J.-C. Lesbordes-Brion, S. Petres, J. Bellalou, and S. A. Susin
AIFsh, a Novel Apoptosis-inducing Factor (AIF) Pro-apoptotic Isoform with Potential Pathological Relevance in Human Cancer
J. Biol. Chem.,
March 10, 2006;
281(10):
6413 - 6427.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. L. Metzger and J. S. Dome
Current Therapy for Wilms' Tumor
Oncologist,
November 1, 2005;
10(10):
815 - 826.
[Abstract]
[Full Text]
[PDF]
|
 |
|