PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2011; 6(2): e16636.
Published online 2011 February 24. doi:  10.1371/journal.pone.0016636
PMCID: PMC3044708

Comparative Profiling of Primary Colorectal Carcinomas and Liver Metastases Identifies LEF1 as a Prognostic Biomarker

John Luk, Editor

Abstract

Purpose

We sought to identify genes of clinical significance to predict survival and the risk for colorectal liver metastasis (CLM), the most common site of metastasis from colorectal cancer (CRC).

Patients and Methods

We profiled gene expression in 31 specimens from primary CRC and 32 unmatched specimens of CLM, and performed Significance Analysis of Microarrays (SAM) to identify genes differentially expressed between these two groups. To characterize the clinical relevance of two highly-ranked differentially-expressed genes, we analyzed the expression of secreted phosphoprotein 1 (SPP1 or osteopontin) and lymphoid enhancer factor-1 (LEF1) by immunohistochemistry using a tissue microarray (TMA) representing an independent set of 154 patients with primary CRC.

Results

Supervised analysis using SAM identified 963 genes with significantly higher expression in CLM compared to primary CRC, with a false discovery rate of <0.5%. TMA analysis showed SPP1 and LEF1 protein overexpression in 60% and 44% of CRC cases, respectively. Subsequent occurrence of CLM was significantly correlated with the overexpression of LEF1 (chi-square p = 0.042), but not SPP1 (p = 0.14). Kaplan Meier analysis revealed significantly worse survival in patients with overexpression of LEF1 (p<0.01), but not SPP1 (p = 0.11). Both univariate and multivariate analyses identified stage (p<0.0001) and LEF1 overexpression (p<0.05) as important prognostic markers, but not tumor grade or SPP1.

Conclusion

Among genes differentially expressed between CLM and primary CRC, we demonstrate overexpression of LEF1 in primary CRC to be a prognostic factor for poor survival and increased risk for liver metastasis.

Introduction

Despite considerable progress in the diagnosis and treatment of colorectal cancer (CRC) over the last few decades leading to a significant decline in cancer-related mortality,[1], [2] CRC remains a major public health problem throughout the world. In the United States, CRC is the third most common cancer and is also the third leading cause of cancer death in men and women combined.[3]

Worldwide it represents the third most common cancer and second most common cause of cancer-related death.[4] Once metastasis has occurred in CRC, a complete cure of the disease is unlikely. Therefore, there is a need for better understanding of the molecular mechanisms underlying the metastatic phenotype that may provide information leading to the development of drugs to control or prevent metastatic disease.[1]

Colorectal liver metastasis (CLM), occurring in about 60% of CRC patients during the course of their treatment, is the most common distant metastasis from CRC. Several clinical prognostic factors, such as lymph node status and size of the primary tumor, have been identified for CLM.[5], [6] However, little is known about the prognostic significance of molecular markers for CLM.

Recent development and application of human genome and high-throughput technologies, such as DNA microarrays, allows us to simultaneously examine thousands of genes, leading to a much better understanding of carcinogenesis - a great step toward individualized personal medicine.[7] Published studies on CRC gene expression profiling have mainly examined normal vs. tumor tissues or different stages of CRC,[8], [9] or treatment outcomes by the differences in gene expression profiling.[10], [11] To identify molecular markers of clinical significance, we used DNA microarrays to compare the gene-expression profiles of specimens from primary CRC and specimens from CLM. We reasoned that genes upregulated in metastasis might also be relatively overexpressed in a subset of clinically-aggressive primary CRC. After identifying upregulated CLM-signature genes, we used tissue microarrays (TMAs) to study the protein expression of selected signature genes in an independent cohort of primary CRC, to correlate their expression with clinical significance and outcome.

Results

Identification of gene signatures distinguishing CLM from primary CRC by expression profiling

To survey the differentially expressed genes between CLM and primary CRC (also compared to normal liver, a potential tissue contaminant of CLM), we used cDNA microarrays containing ~19,500 unique genes to profile the gene expression in 31 primary CRC specimens from 30 patients, and 32 unmatched CLM specimens from 31 patients who underwent liver resection. We then performed supervised analysis using SAM (with a Student's t-test metric) and identified 1,186 discriminatory cDNAs (corresponding to 963 unique genes) with significantly higher expression in CLM when compared to primary CRC, and to normal liver (previously profiled,[12]), with a false discovery rate (FDR) of <0.005% (Table S2). The top 35 differentially expressed genes are shown in Fig. 1. The 20 highest-ranking genes were SPP1, CXCR4, GPNMB, LOX, CD53, AIF1, ARHGDIB, SLC12A2, PRG1, SPARC, CD3D, DZIP1, PEG3, FYB, ITM2A, SLA, IGLC2, MGP, LEF1, and MAF. Similar results were obtained using SAM with a non-parametric, Wilcoxin rank-based analysis (Table S3).

Figure 1
Heatmap showing the top 35 ranking genes, based on SAM analysis, with greatest increased expression in CRC with CLM compared to primary CRC or normal liver.

LEF1 protein overexpression in CRC correlates with CLM and overall survival

Primary CRC specimens with relative increased expression of CLM signature genes might exhibit increased metastatic potential. To further study the significance of such potential biomarkers identified via SAM analysis, we examined the protein expression of two highly ranked and biologically plausible signature genes (for which IHC-validated antibodies were also available), SPP1 (linked to metastasis[13]) and LEF1 (lymphoid enhancer factor-1; involved in WNT signaling[14]), in an independent TMA cohort of CRC specimens. Examples of antibody staining by intensity are shown in Fig. 2. Correlation of SPP1 or LEF1 overexpression (IHC scored +2 or +3) with subsequently occurring CLM is shown in Table 1. LEF1 overexpression was significantly correlated with CLM (Fisher's exact test, p = 0.042). No significant correlation was observed between SPP1 overexpression and subsequent CLM.

Figure 2
Immunohistochemical staining and assessment of SPP1 or LEF1 in representative core sections from the CRC tissue microarray.
Table 1
Correlation between liver metastasis and SPP1 or LEF1 overexpressiona.

We used Kaplan-Meier analyses to investigate the impact of SPP1 or LEF1 overexpression on overall survival. LEF1 overexpression was significantly associated with worse survival (log-rank p<0.01; Fig. 3B). In contrast, SPP1 overexpression was not significantly correlated with a worse outcome, though there was a trend toward poor survival (log-rank p = 0.11; Fig. 3A)

Figure 3
Kaplan-Meier overall survival curves based on (A) SPP1 overexpression (log rank p = 0.11) or (B) LEF1 overexpression (log rank p<0.01).

LEF1 is a significant prognostic factor in multivariate analysis

As expected, higher stage was found to be significantly associated with worse overall survival. In the univariate model, the hazard ratios (HRs) for stage 2, 3, and 4 vs. 1 was 1.63, 3.47 (p = 0.02) and 12.63 (p<0.001), respectively; in the multivariate model, the HRs are 1.74, 3.74 (p = 0.01), and 13.74 (p<0.001), respectively. The HRs were also statistically significant for LEF1 overexpression in both the univariate and multivariate models - 1.78 (p<0.05) and 1.66 (p<0.05), respectively. Neither SPP1 overexpression nor tumor grade was significantly associated with survival in either analysis. Table 2 summarizes the HRs and 95% confidence intervals for the variables in both models.

Table 2
Univariate and Multivariate Hazard Ratios From Cox Proportional Hazards Regression Models.

Discussion

In our study, we sought to identify signatures of metastasis embedded in a subset of primary tumors, which might predict clinically-aggressive behavior [15]. Using supervised SAM analysis, we identified 963 unique genes that are significantly overexpressed in CLM vs. primary CRC (and potentially contaminating normal liver tissue). In an independent set of tissue microarrays, we examined two highly-ranked genes (LEF1 and SPP1) as surrogate biomarkers for the CLM signature, and demonstrated that overexpression of LEF1, but not SPP1, in the primary CRC tissues correlates with a statistically significant increased risk of CLM, albeit its sensitivity and specificity in predicting liver metastasis were modest. In addition, independent of tumor stage, overexpression of LEF1, not SPP1, denotes a poor prognosis for survival.

LEF1 was initially identified as a pre-B and T-lymphoid-specific gene encoding a DNA-binding protein of high mobility group (HMG) proteins.[16], [17] It is a member of the T-cell factor/lymphoid-enhancing factor (TCF/LEF) family of transcription factors, which acts through the Wnt signaling pathway[18], [19], [20] to regulate gene expression and coordinate many cellular processes in normal development and tissue homeostasis, and, when deregulated, in colonic tumorogenesis and metastasis. Upon Wnt stimulation, LEF1 or other TCF/LEF-family transcription factors associate with β-catenin, a key cytoplasmic/nuclear mediator of Wnt pathway, and activate Wnt-responsive target genes. In contrast, without Wnt stimulation, glycogen synthase kinase (GSK)-3 (in a complex with APC) constitutively phosphorylates β-catenin, resulting in its proteasome-dependent degradation.[21] Although genetic and epigenetic changes have been documented in several targets throughout the pathway, mutation in either APC or β-catenin appears to be a crucial element in CRC carcinogenesis.[22] The LEF1 gene itself is not normally expressed in the adult intestinal epithelium, but only observed in the embryos while development is in progress. However, its overexpression has been well documented in CRC tumorigenesis[23], and denotes aberrant activation of the Wnt/β-catenin pathway via stabilization of LEF1/β-catenin complex. Indeed, LEF1 is a direct transcriptional target of the LEF1/β-catenin complex[23], indicating a positive feedback loop for Wnt signaling, and suggesting LEF1 might be useful as a prognostic biomarker of Wnt pathway activation. Our findings are consistent with a prior study showing nuclear (active) β-catenin staining to be prognostic in colon cancer,[24] and highlight a role of Wnt signaling in colon cancer progression and liver metastasis.

As the prefix "osteo" suggests, osteopontin (or OPN, also known as secreted phosphoprotein 1 [SPP1], bone sialoprotein I, early T-lymphocyte activation 1) was initially recognized as an important glycosylated, adhesive phosphoprotein in bone.[25] Since then, several lines of evidence have shown its role in controlling tumorigenicity, progression and metastasis via its diverse ability as a cell-matrix mediator to interact with a variety of factors such as cell surface receptors (integrins and CD44), secreted proteases (matrix metalloproteinases and urokinase plasminogen activator), and growth factor/receptor pathways (TGF/EGFR and HGF/Met).[26] Overexpression of SPP1 has been reported in several human cancers, including lung, breast and colon cancers. In a gene-expression profiling study, Agrawal et al identified SPP1 as a lead marker correlating with CRC progression, and strongly expressed in CLM.[27], [28], [29] Most recently, Rohde et al observed that overexpression of SPP1 is indicative of poor survival in CRC and is significantly correlated with CLM.[30] In addition, overexpression of SPP1 correlates with increased immunohistochemical staining of β-catenin and, in an in vivo model, with Wnt activating mutations. These data suggest a crucial role of SPP1 in CRC progression and metastasis likely via molecular cross-talk with the Wnt pathway. Our results partially corroborate with published data in that we observe a trend towards poorer survival and CLM with SPP1 overexpression, though it is not statistically significant. This difference may be due to a smaller sample size and/or shorter follow-up time. Nonetheless, our data suggest that the overexpression of LEF1 is a stronger prognostic factor than SPP1 in correlating overall survival and CLM.

Though our analysis focused primarily on LEF1 and SPP1, other highly-ranked signature genes with increased expression in CLM compared to primary CRC also have biological functions consistent with roles in tumor progression, and might have prognostic utility. For example, CXCR4 (Chemokine (C-X-C motif) receptor 4) has been implicated in breast [31] and colon cancer metastasis.[32] LOX (lysyl oxidase) was shown to be associated with hypoxia where it functions in metastasis and predicts poor outcome in breast cancer.[33] SPARC (secreted protein, acidic, cysteine-rich; osteonectin) was identified among genes that mark and mediate breast cancer metastasis.[34] Further studies are needed to characterize these and other signature genes in CRC progression.

Much of the published data comparing gene expression profiles from primary CRC vs. CLM have reported the differences in up- or down-regulated genes.[35], [36], [37], [38], [39], [40], [41], [42], [43], [44] To the best of our knowledge, our data are the first to identify the expression of LEF1 as a predictor of overall survival, and an indicator of CLM. Comparisons of primary CRC vs. CLM, or of primary CRC associated with or without CLM, have identified gene signatures with relevance to colorectal cancer progression.[27], [30], [35], [36], [37], [38], [39], [40], [41], [43] However, there are minimal overlaps between our CLM signature genes and the top-ranking CLM/progression-associated genes reported in these studies (Table 3). Several issues may explain the discordance among studies, such as: variability in patient cohorts, technical differences (the composition of the microarrays used, study design, statistical methodologies), and variability in the use of independent cohorts of patients to validate candidate prognostic genes.[45]

Table 3
Overlap of CLM genes with top-ranking CLM/progression genes from other published studies.

In conclusion, our study shows that overexpression of LEF1 in primary CRC correlates with a higher risk of CLM and denotes poor overall survival. It is a stronger predictor than SPP1, a marker reported in previous transcriptome studies. High-throughput gene expression profiling technology has revealed new insights into the molecular heterogeneity of CRC and identified new and better molecular markers for risk stratification. This holds promise for personalized medicine and improved targeted therapy. To achieve these goals, further studies are needed to understand the functional roles and clinical implications of LEF1, SPP1 and other signature genes for CLM.

Materials and Methods

Patients and Tumor Specimens

Freshly frozen CRC specimens (from Santa Clara Valley Medical Center) and CLM specimens (from Stanford University Medical Center) were used for gene expression profiling analysis. We used 31 primary CRC specimens from 30 patients (14 males and 16 females; age range: 36 to 80; stage I/II/III/IV = 1/10/11/8, diagnosed between 2000 and 2004), and 32 CLM specimens from 31 patients (16 males and 15 females; age range: 40 to 82). Tissue microarrays were constructed with paraffin blocks from 154 CRC cases (84 males and 70 females; age range: 27 to 92; stage I/II/III/IV = 26/46/38/44) with median follow-up of 3.0 years (range: 0.6–15.3 years). Study protocols were approved by Institutional Review Board (IRB) both at Stanford University Medical Center (IRB # 12473) and Santa Clara Valley Medical Center, a Stanford-affiliated teaching hospital, (IRB #07/28/00-03 and #2/22/2002-04). Individual informed consent was obtained from all participants involved in the gene expression profiling protocols and was waived in the TMA study by IRB due to the use of coded data and retrospective nature of the study. Clinical characteristics of the patient cohorts are summarized in Table S1.

Gene expression profiling

To confirm that the sample was representative of the case, a frozen section from each specimen was first prepared and examined. Tissue was then homogenized in Trizol reagent (Invitrogen, Carlsbad, CA), and total RNA isolated per the manufacturer's protocol. RNA quality was assessed by gel electrophoresis. Gene expression profiling was performed as described previously.[46] Briefly, using microarrays of complementary DNA (cDNA) (manufactured by the Stanford Functional Genomics Facility) containing ~40,000 nonredundant cDNA clones, representing ~19,500 unique UniGene clusters (i.e., genes), we hybridized Cy5-labeled total RNA from the tumor specimens, along with Cy3-labeled universal reference mRNA (pooled from 11 different cancer cell lines). We imaged arrays using an Axon GenePix 4000B scanner (Molecular Devices, Sunnyvale, CA), extracted fluorescence ratios (ratio of the specimen value to the reference value) using the GenePix software, and entered the data into the Stanford Microarray Database[47] for subsequent analysis. The microarray data are accessible from the Gene Expression Omnibus (Accession GSE22834). cDNA microarray expression data from normal liver specimens were previously published.[12]

TMA construction and immunohistochemistry

A tissue arrayer (Beecher Instruments, Sun Prarie, WI) was used to construct a primary CRC tissue microarray as described,[48] comprising of 154 primary colorectal tumors each represented by two 6-mm cores. A 4-µm section was cut from the tissue microarray block, deparaffinized in Citrisolv (Fisher Scientific, Hampton, NH) and hydrated in a graded series of alcohol solutions.

For immunohistochemical staining (IHC), anti-osteopontin (SSP1) mouse monoclonal antibody (Novocatra, Newcastle, UK) and anti-LEF1 (lymphoid enhancer factor-1) rabbit polyclonal antibody (Abcam, Cambridge, MA) were used at 1[ratio]100 and 1[ratio]500 dilutions respectively, and incubated overnight at 4°C. Chromogenic detection was then done using a peroxidase-conjugated secondary antibody and DAB reagents provided with the Envision detection kit (DAKO, Carpinteria, CA). SPP1 expression was scored as positive if distinct cytoplasmic staining was present in more than 10% of tumor cells. LEF1 expression was scored as positive when distinct nuclear staining was present in more than 10% of tumor cells. Weak cytoplasmic staining of LEF1 was not counted. The staining intensity was graded on a semiquantitative score (0, negative; 1+, weak; 2+, moderate; and 3+, strong). Survival analysis was performed in two groups depending on the score: overexpression (2+ and 3+) vs. the remainder (0 or 1+). Immunostains were scored by two pathologists (M.vdR. and YLC) blinded to the clinical data.

Statistical Analysis

For cDNA microarray data, ratios were globally normalized by array and median-centered by gene. We included for analysis the 4,824 cDNAs (corresponding to 3,413 unique genes) that were well-measured (intensity/background >2 in either the test or reference channel) in at least 50% of samples, and variably expressed (>4-fold change from the median) in at least 3 samples. Two-class Significance Analysis of Microarrays (SAM) [49] was used to identify genes that were differentially expressed in CLM compared to primary CRC and normal liver (a potential contaminant of CLM), with statistical significance assessed by a false discovery rate (FDR).

For clinicopathological data, a chi-square test or Fisher's exact test (two-tailed) was used to compare differences in categorical variables across patient groups. Kaplan-Meier methods were used to estimate overall survival. We used both univariate and multivariate Cox proportional hazard regression model to assess the prognostic independence of variables (LEF1, SPP1, stage and tumor grade) for survival. Statistical analyses were performed with the SAS System software, release 9.1.3 (SAS Institute Inc., Cary, NC).

Supporting Information

Table S1

Summary of clinical characteristics of CRC cohorts.

(XLS)

Table S2

Ranked list of significant (FDR<0.005%) SAM genes (Student's t-test) upregulated in CLM vs. primary CRC and normal liver.

(XLS)

Table S3

Ranked list of significant (FDR<0.005%) SAM genes (Wilcoxin rank) upregulated in CLM vs. primary CRC and normal liver.

(XLS)

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: Supported in part by NIH grant CA112016 (JRP); The Valley Foundation (AYL); Lui Hac Minh and C. J. Huang Foundations and the Paul and Tammy Davis Research Fund (MSC and SKS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Markowitz SD, Bertagnolli MM. Molecular Basis of Colorectal Cancer. N Engl J Med. 2009;361:2449–2460. [PMC free article] [PubMed]
2. Jemal A, Ward E, Thun M. Declining Death Rates Reflect Progress against Cancer. PLoS One. 2010;5:e9584. [PMC free article] [PubMed]
3. Jemal A, Siegel R, Ward E, Hao Y, Xu J, et al. Cancer Statistics, 2009. CA Cancer J Clin. 2009;59:225–249. [PubMed]
4. Greenlee RT, Hill-Harmon MB, Murray T, Thun M. Cancer statistics, 2001. CA Cancer J Clin. 2001;51:15–36. [PubMed]
5. Fong Y, Fortner J, Sun RL, Brennan MF, Blumgart LH. Clinical score for predicting recurrence after hepatic resection for metastatic colorectal cancer: analysis of 1001 consecutive cases. Ann Surg. 1999;230:309–318. [PubMed]
6. Jarnagin WR, Conlon K, Bodniewicz J, Dougherty E, DeMatteo RP, et al. A clinical scoring system predicts the yield of diagnostic laparoscopy in patients with potentially resectable hepatic colorectal metastases. Cancer. 2001;91:1121–1128. [PubMed]
7. Pollack JR. A Perspective on DNA Microarrays in Pathology Research and Practice. Am J Pathol. 2007;171:375–385. [PubMed]
8. Notterman DA, Alon U, Sierk AJ, Levine AJ. Transcriptional gene expression profiles of colorectal adenoma, adenocarcinoma, and normal tissue examined by oligonucleotide arrays. Cancer Res. 2001;61:3124–3130. [PubMed]
9. Friederichs J, Rosenberg R, Mages J, Janssen K-P, Maeckl C, et al. Gene expression profiles of different clinical stages of colorectal carcinoma: toward a molecular genetic understanding of tumor progression. Int J Colorectal Dis. 2005;20:391–402. [PubMed]
10. Inokuchi M, Uetake H, Shirota Y, Yamada H, Tajima M, et al. Gene expression of 5-fluorouracil metabolic enzymes in primary colorectal cancer and corresponding liver metastasis. Cancer Chemother Pharmacol. 2004;53:391–396. [PubMed]
11. Matsuyama R, Togo S, Shimizu D, Momiyama N, Ishikawa T, et al. Predicting 5-fluorouracil chemosensitivity of liver metastases from colorectal cancer using primary tumor specimens: three-gene expression model predicts clinical response. Int J Cancer. 2006;119:406–413. [PubMed]
12. Chen X, Cheung ST, So S, Fan ST, Barry C, et al. Gene expression patterns in human liver cancers. Mol Biol Cell. 2002;13:1929–1939. [PMC free article] [PubMed]
13. Rangaswami H, Bulbule A, Kundu GC. Osteopontin: role in cell signaling and cancer progression. Trends in Cell Biology. 2006;16:79–87. [PubMed]
14. Li TW-H, Ting J-HT, Yokoyama NN, Bernstein A, van de Wetering M, et al. Wnt Activation and Alternative Promoter Repression of LEF1 in Colon Cancer. Mol Cell Biol. 2006;26:5284–5299. [PMC free article] [PubMed]
15. Ramaswamy S, Ross KN, Lander ES, Golub TR.Nature Genetics. Nature Publishing Group; 2003. A molecular signature of metastasis in primary solid tumors. [PubMed]
16. Travis A, Amsterdam A, Belanger C, Grosschedl R. LEF-1, a gene encoding a lymphoid-specific protein with an HMG domain, regulates T-cell receptor alpha enhancer function [corrected]. Genes Dev. 1991;5:880–894. [PubMed]
17. Waterman ML, Fischer WH, Jones KA. A thymus-specific member of the HMG protein family regulates the human T cell receptor C alpha enhancer. Genes Dev. 1991;5:656–669. [PubMed]
18. Cadigan KM, Nusse R. Wnt signaling: a common theme in animal development. Genes Dev. 1997;11:3286–3305. [PubMed]
19. Bienz M, Clevers H. Linking colorectal cancer to Wnt signaling. Cell. 2000;103:311–320. [PubMed]
20. Nusse R. WNT targets: repression and activation. Trends in Genetics. 1999;15:1–3. [PubMed]
21. Fuerer C, Nusse R, Ten Berge D. Wnt signalling in development and disease. Max Delbruck Center for Molecular Medicine meeting on Wnt signaling in Development and Disease. EMBO Rep. 2008;9:134–138. [PubMed]
22. Schneikert J, Behrens J. The canonical Wnt signalling pathway and its APC partner in colon cancer development. Gut. 2007;56:417–425. [PMC free article] [PubMed]
23. Hovanes K, Li TWH, Munguia JE, Truong T, Milovanovic T, et al. β-catenin-sensitive isoforms of lymphoid enhancer factor-1 are selectively expressed in colon cancer. Nature Genetics. 2001;28:53. [PubMed]
24. Wong SC, Lo ES, Lee KC, Chan JK, Hsiao WL. Prognostic and diagnostic significance of beta-catenin nuclear immunostaining in colorectal cancer. Clin Cancer Res. 2004;10:1401–1408. [PubMed]
25. Prince C, Oosawa T, Butler W, Tomana M, Bhown A, et al. Isolation, characterization, and biosynthesis of a phosphorylated glycoprotein from rat bone. J Biol Chem. 1987;262:2900–2907. [PubMed]
26. Tuck AB, Chambers AF, Allan AL. Osteopontin overexpression in breast cancer: Knowledge gained and possible implications for clinical management. Journal of Cellular Biochemistry. 2007;102:859–868. [PubMed]
27. Agrawal D, Chen T, Irby R, Quackenbush J, Chambers AF, et al. Osteopontin identified as lead marker of colon cancer progression, using pooled sample expression profiling. J Natl Cancer Inst. 2002;94:513–521. [PubMed]
28. Agrawal D, Chen T, Irby R, Quackenbush J, Chambers AF, et al. Osteopontin identified as colon cancer tumor progression marker. C R Biol. 2003;326:1041–1043. [PubMed]
29. Yeatman TJ, Chambers AF. Osteopontin and colon cancer progression. Clin Exp Metastasis. 2003;20:85–90. [PubMed]
30. Rohde F, Rimkus C, Friederichs J, Rosenberg R, Marthen C, et al. Expression of osteopontin, a target gene of de-regulated Wnt signaling, predicts survival in colon cancer. Int J Cancer. 2007;121:1717–1723. [PubMed]
31. Muller A, Homey B, Soto H, Ge N, Catron D, et al. Involvement of chemokine receptors in breast cancer metastasis. Nature. 2001;410:50–56. [PubMed]
32. Zeelenberg IS, Ruuls-Van Stalle L, Roos E. The chemokine receptor CXCR4 is required for outgrowth of colon carcinoma micrometastases. Cancer Res. 2003;63:3833–3839. [PubMed]
33. Erler JT, Bennewith KL, Nicolau M, Dornhofer N, Kong C, et al. Lysyl oxidase is essential for hypoxia-induced metastasis. Nature. 2006;440:1222–1226. [PubMed]
34. Minn AJ, Gupta GP, Siegel PM, Bos PD, Shu W, et al. Genes that mediate breast cancer metastasis to lung. Nature. 2005;436:518–524. [PMC free article] [PubMed]
35. Ki DH, Jeung HC, Park CH, Kang SH, Lee GY, et al. Whole genome analysis for liver metastasis gene signatures in colorectal cancer. Int J Cancer. 2007;121:2005–2012. [PubMed]
36. Kleivi K, Lind GE, Diep CB, Meling GI, Brandal LT, et al. Gene expression profiles of primary colorectal carcinomas, liver metastases, and carcinomatoses. Mol Cancer. 2007;6:2. [PMC free article] [PubMed]
37. Li M, Lin YM, Hasegawa S, Shimokawa T, Murata K, et al. Genes associated with liver metastasis of colon cancer, identified by genome-wide cDNA microarray. Int J Oncol. 2004;24:305–312. [PubMed]
38. Lin HM, Chatterjee A, Lin YH, Anjomshoaa A, Fukuzawa R, et al. Genome wide expression profiling identifies genes associated with colorectal liver metastasis. Oncol Rep. 2007;17:1541–1549. [PubMed]
39. Nakagawa H, Liyanarachchi S, Davuluri RV, Auer H, Martin EW, Jr, et al. Role of cancer-associated stromal fibroblasts in metastatic colon cancer to the liver and their expression profiles. Oncogene. 2004;23:7366–7377. [PubMed]
40. Tackels-Horne D, Goodman MD, Williams AJ, Wilson DJ, Eskandari T, et al. Identification of differentially expressed genes in hepatocellular carcinoma and metastatic liver tumors by oligonucleotide expression profiling. Cancer. 2001;92:395–405. [PubMed]
41. Takahashi Y, Ishii Y, Nishida Y, Ikarashi M, Nagata T, et al. Detection of aberrations of ubiquitin-conjugating enzyme E2C gene (UBE2C) in advanced colon cancer with liver metastases by DNA microarray and two-color FISH. Cancer Genet Cytogenet. 2006;168:30–35. [PubMed]
42. Takata O, Kawamura YJ, Konishi F, Sasaki J, Kai T, et al. cDNA Array Analysis for Prediction of Hepatic Metastasis of Colorectal Carcinoma. Surgery Today. 2006;36:608–614. [PubMed]
43. Takayama T, Miyanishi K, Hayashi T, Sato Y, Niitsu Y. Colorectal cancer: genetics of development and metastasis. J Gastroenterol. 2006;41:185–192. [PubMed]
44. Yanagawa R, Furukawa Y, Tsunoda T, Kitahara O, Kameyama M, et al. Genome-wide screening of genes showing altered expression in liver metastases of human colorectal cancers by cDNA microarray. Neoplasia. 2001;3:395. [PMC free article] [PubMed]
45. Quackenbush J. Microarray analysis and tumor classification. N Engl J Med. 2006;354:2463–2472. [PubMed]
46. Lapointe J, Li C, Higgins JP, van de Rijn M, Bair E, et al. Gene expression profiling identifies clinically relevant subtypes of prostate cancer. Proc Natl Acad Sci U S A. 2004;101:811–816. [PubMed]
47. Demeter J, Beauheim C, Gollub J, Hernandez-Boussard T, Jin H, et al. The Stanford Microarray Database: implementation of new analysis tools and open source release of software. Nucleic Acids Res. 2007;35:D766–770. [PMC free article] [PubMed]
48. Kononen J, Bubendorf L, Kallioniemi A, Barlund M, Schraml P, et al. Tissue microarrays for high-throughput molecular profiling of tumor specimens. Nat Med. 1998;4:844–847. [PubMed]
49. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98:5116–5121. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science