Pancreatic cancer circulating tumour cells express a cell motility gene signature that predicts survival after surgery
© Sergeant et al; licensee BioMed Central Ltd. 2012
Received: 16 July 2012
Accepted: 8 November 2012
Published: 16 November 2012
Most cancer deaths are caused by metastases, resulting from circulating tumor cells (CTC) that detach from the primary cancer and survive in distant organs. The aim of the present study was to develop a CTC gene signature and to assess its prognostic relevance after surgery for pancreatic ductal adenocarcinoma (PDAC).
Negative depletion fluorescence activated cell sorting (FACS) was developed and validated with spiking experiments using cancer cell lines in whole human blood samples. This FACS-based method was used to enrich for CTC from the blood of 10 patients who underwent surgery for PDAC. Total RNA was isolated from 4 subgroup samples, i.e. CTC, haematological cells (G), original tumour (T), and non-tumoural pancreatic control tissue (P). After RNA quality control, samples of 6 patients were eligible for further analysis. Whole genome microarray analysis was performed after double linear amplification of RNA. ‘Ingenuity Pathway Analysis’ software and AmiGO were used for functional data analyses. A CTC gene signature was developed and validated with the nCounter system on expression data of 78 primary PDAC using Cox regression analysis for disease-free (DFS) and overall survival (OS).
Using stringent statistical analysis, we retained 8,152 genes to compare expression profiles of CTC vs. other subgroups, and found 1,059 genes to be differentially expressed. The pathway with the highest expression ratio in CTC was p38 mitogen-activated protein kinase (p38 MAPK) signaling, known to be involved in cancer cell migration. In the p38 MAPK pathway, TGF-β1, cPLA2, and MAX were significantly upregulated. In addition, 9 other genes associated with both p38 MAPK signaling and cell motility were overexpressed in CTC. High co-expression of TGF-β1 and our cell motility panel (≥ 4 out of 9 genes for DFS and ≥ 6 out of 9 genes for OS) in primary PDAC was identified as an independent predictor of DFS (p=0.041, HR (95% CI) = 1.885 (1.025 – 3.559)) and OS (p=0.047, HR (95% CI) = 1.366 (1.004 – 1.861)).
Pancreatic CTC isolated from blood samples using FACS-based negative depletion, express a cell motility gene signature. Expression of this newly defined cell motility gene signature in the primary tumour can predict survival of patients undergoing surgical resection for pancreatic cancer.
Clinical trials.gov NCT00495924
KeywordsCirculating tumour cells Pancreatic ductal adenocarcinoma Gene expression profiling p38 – MAPK signaling Transforming growth factor - β1 Cancer cell migration Cell motility
The vast majority of patients suffering from solid organ malignancies ultimately die from metastases, irrespective of the type of treatment modality . Current cancer treatments are mainly based on primary tumour characterization instead of on characterization of metastases or cancer cells in the blood circulation. Since distant organ metastases are the end-results of circulating tumor cells (CTC), molecular characterization of CTC is promising in the quest for effective therapeutic agents and novel prognostic markers .
Unfortunately, specific isolation of CTC has proved difficult. Identification and isolation of CTC using predefined marker expression has failed to identify certain cancer cell phenotypes, and consequently also their genotypes. Methods to isolate, identify and enable reliable and reproducible molecular characterization of CTC without a priori selection could provide new insights into cancer research and improve management of the disease.
The aim of the present study was to develop a negative depletion strategy to enrich whole blood for CTC without excluding any particular subset of CTC, and to study their gene-expression profiles in order to find novel therapeutic targets and prognostic markers in patients with pancreatic cancer.
Patients and tissues
Clinicopathological characteristics of included patients
AJCC classif. (2002)
Fluorescence activated cell sorting (FACS)
The Panc-1 pancreatic cancer cell line was purchased from the American Type Culture Collection (ATCC; CRL-1469, ABPromochem) . Panc-1 cells were cultured at 37°C in DMEM supplemented with 10% foetal bovine serum, penicillin and streptomycin.
Two different concentrations of Panc-1 cells (5.0 × 104 and 5.0 × 106 cells) were spiked in duplo in 20ml of EDTA-treated blood obtained from healthy volunteers, and gently versed in an Oncoquick® tube. The recovered cell populations (CTC- and G-fractions) were FACS-sorted in CytoRich® red and cytospins were made by spinning the cells at 800 rpm for 10min. Slides were air-dried for 24 hours and a pre-keratin stain (anti-CK8/18, 1:30, Novacastra) was carried out. Cells were counterstained with haematoxylin. An experienced pathologist (TR) semi-quantitatively scored the relative distribution of cell populations.
RNA was isolated from tissue sections using a protocol combining Trizol/chloroform extraction, followed by column chromatography with the RNeasy Mini kit (Qiagen). Briefly, after homogenization of the tissue with 500μl of Trizol® (Invitrogen), 100μl of chloroform was added followed by mixing until a milky solution was obtained.
This mixture was then centrifuged at 12.000 rpm in a cooled benchtop centrifuge for 25 minutes. Subsequently, the supernatant was removed and added to an equal amount of 70% ethanol. After this step the RNeasy Mini kit was followed according to the manufacturer’s recommendations.
RNA from the FACS sorted cell populations was prepared using the RNeasy Micro kit (Qiagen), including a DNase step according to the manufacturer’s recommendations.
RNA concentration and purity were determined spectrophotometrically using the Nanodrop ND-1000 (Nanodrop Technologies) for all samples. RNA integrity of the T and P samples was assessed using a Bioanalyser 2100 (Agilent). Per sample, ideally an amount of 5–10 ng of total RNA spiked with four bacterial RNA transcripts (Affymetrix) was converted and amplified to double-stranded cDNA in a 2-cycle cDNA reverse transcription reaction. Subsequently, the sample was converted to antisense cRNA and labelled with biotin through an in vitro transcription reaction according to the manufacturer’s protocol (Affymetrix). Purified fragmented biotinylated cRNA and hybridization controls (Affymetrix) were mixed, hybridized on Affymetrix HG U133 Plus 2.0 arrays, and subsequently stained and washed in the GeneChip fluidics station 450 (Affymetrix). To assess the raw probe signal intensities, chips were scanned using the GeneChip scanner 3000 (Affymetrix).
Microarray data analysis
R, a free software environment for statistical computing and graphics, was used in combination with the affy and limma libraries in Bioconductor for microarray data analysis [6–9]. To decide whether a signal was significantly above background, the MAS 5.0 algorithm was applied to calculate probe set detection calls. Robust Multichip Average (RMA) was applied to probe sets that had a present detection call in at least 4 out of 6 CCD samples. Using limma, the average expression value for each experimental condition was estimated. Based on these estimates, the contrasts CTC vs. T, CTC vs. P, CTC vs. G, T vs. P, T vs. G, and P vs. G were estimated. A moderated t-statistic (implemented in limma) was used to test whether each contrast was significantly different from 0 Resulting p-values were corrected for multiple testing with Benjamini-Hochberg (BH) to control for false discovery rates . Subsequently, a filter was applied selecting genes that were significantly 2-fold up– AND/OR downregulated using the BH corrected p-values (p < 0.05) in the comparisons CTC vs. T, AND in CTC vs. P, AND CTC vs. G. The resulting filtered data were functionally analyzed with the Ingenuity Pathway Analysis (IPA) application (Ingenuity® Systems, http://www.ingenuity.com). Gene networks were generated within the Ingenuity Pathway Analysis programme, and ordered by a score. This score is a numerical value used to rank networks according to their degree of relevance to the Network Eligible molecules in the dataset. The network score is based on hypergeometric distribution and is calculated with the Ingenuity Pathway Analysis programme using the right-tailed Fisher's Exact Test. The set of genes upregulated in CTC vs. T, AND in CTC vs. P, AND CTC vs. G was further explored for genes associated with cell motility using AmiGO, the GO Consortium’s annotation and ontology tool kit . Gene expression data are publicly available online via: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=zbotdukywssgujm&acc=GSE18670.
Validation of gene signature
Between 2006 and 2010, tissue samples were obtained from patients who had undergone pancreatic resection for PDAC. Samples were stored at −80°C in RNAlater® (Qiagen). To verify that the tumor samples contained more than 70% tumor, hematoxylin and eosin stains were made from each tumour sample and from surrounding control pancreatic tissue to confirm its non-cancerous histology. True cellularity was not determined. From the primary tumour of 143 patients and from surrounding non-tumoural pancreatic (control) tissue of 14 patients, total RNA was extracted using the RNeasy Mini kit according the manufacturer’s instructions. Only samples with an RNA integrity number (RIN) of ≥ 7.0 were used for further analysis, i.e. 78 PDAC-samples (M/F = 40/38, median age 64 y. (range 32–80 y.)) and 6 controls. A limited set of genes was selected based on the results from the microarray gene expression analysis in CTC. Using the nCounter system (Nanostring Technologies, Seattle, WA) gene expression levels were quantified at the VIB Nucleomics Core of the VIB Microarray Facility.
Median overall survival (OS) and disease-free survival (DFS) rates in our patient group were 18.7 months (95% CI: 12.4-25.6 months) and 10.0 months (95% CI: 7.4-12.4 months), respectively. Uni- and multivariable Cox regression analyses were performed to identify predictors of DFS and OS rates after pancretic resection with curative intent. The common closing date was the 1st of July 2011. Variables considered were dichotomized expression values of the genes of interest and clinico-pathological variables: tumour localisation in the pancreas (head/corpus/tail), tumour size (< 20mm, 20-40mm, > 40mm), tumour differentiation grade (pG), pT, pN, pM, extracapsular lymph node involvement (ECLNI), perineural invasion (PNI), lymphovascular invasion (LVI), vascular invasion (VI), resection margin status (R), minimal resection margin > 1mm, adjuvant chemotherapy within the first 3 months after surgery. The cut-off for dichotomisation of gene expression values was the median fold change calculated for every gene. All variables with p-values ≤ 0.10 in univariable analysis were considered for the final multivariable model. A p-value ≤ 0.05 was considered statistically significant.
Similarly, 5.0 × 104 Panc-1 cells (low concentration) were spiked in 20ml of EDTA-treated whole blood. The cytospins’ cell density was greatly reduced owing to the lower concentration of spiked Panc-1 cells and lower recovery of spiked cells. Despite varying cell size, marked cytoplasmatic prekeratin staining in the CTC-fraction could be demonstrated (Figure 3C). On the other hand, the G-fraction was composed almost exclusively of PBMCs similar to the G-fraction after sorting of the high Panc-1 concentration. Estimation of the recovery after spiking, density centrifugation, FACsorting and cytospins was not attempted.
RNA purity and integrity
Similar to enrichment procedures used in our spiking experiments, blood of ten patients (Table 1) was sorted into two fractions (CTC and G). A median number of 1,828 (range: 761 – 33,337) presumed CTC was obtained for RNA extraction. A260/A280 ratios were > 1.80 for all samples. Tissue samples (T and P) were obtained with good RNA integrity numbers (RIN) ≥ 7.5, except for samples T123 and P138, which had a RIN of 7.1 and 6.9, respectively. The RNA quality standards set for microarray analysis were reached in 6 out of 10 patients for all 4 subgroup samples, i.e. for the CTC, G, T and P samples the median starting amount of RNA converted and amplified to antisense cRNA was 0.315ng (range: 0.22 – 0.49), 8.00ng (range: 0.528 – 10), 10ng (range: 5–10) and 10ng (range: 5 – 10), respectively.
Gene expression analysis
Counts of down- and upregulated genes
Corrected p-value < 0.05
Log-ratio < −1
Log-ratio > +1
CTC vs T
CTC vs P
CTC vs G
T vs P
T vs G
P vs G
Molecular network and pathway analysis
To avoid probes referring to the same molecule in the IPA application, a maximum of fold-change values was taken. In order to compare CTC vs. T, AND CTC vs. P, AND CTC vs. G fractions, only genes with a significant 2-fold up– or downregulation were selected (corrected p-values < 0.05). This selection resulted in a total of 1,059 probe sets. Among these 1,059 differentially expressed genes, 738 molecules were eligible to generate gene networks with the Ingenuity Pathway Analysis software, and a total of 572 molecules were eligible for the generation of cellular functions and pathways. Primary functions involved in the top gene network were cancer, cell death, and neurological disease, with transforming growth factor β1 (TGF-β1) as the core molecule (Additional file 2) TGF-β1 was found to be highly upregulated in the CTC fraction compared to the T, P, and G fractions, with a fold change (FC) of 3.18, 3.33, and 2.16, respectively. The three most important biofunctions for the category ‘Diseases and Disorders’ were inflammatory response (6.6E-04 < p-value < 4.81E-02), cancer (1.69E-03 < p-value < 3.92E-02), and genetic disorders (1.69E-03 < p-value < 3.92E-02). The three most important ‘Molecular and Cellular functions’ were cell-to-cell signaling and interaction (7.84E-04 < p-value < 4.98E-02), cellular development (2.47E-03 < p-value < 4.81E-02), and cellular movement (5.22E-03 < p-value < 4.74E-02). P-values were calculated with the Ingenuity Pathway Analysis program using the right-tailed Fisher's Exact Test and represent a measure of its statistical significance with respect to the functions for the dataset and the references of molecules (which define the molecules that may have been functions eligible).
The pathway with the highest expression ratio (10.5%) was p38 mitogen-activated protein kinase (MAPK) signaling. This ratio was calculated as the number of molecules in a given pathway that met cut-off criteria (n=10) divided by the total number of molecules making up that pathway (n=95). In the p38 MAPK pathway TGF-β1, cytosolic phospholipase A2 (cPLA2), and MYC associated factor X (MAX) were significantly upregulated in the CTC fraction as compared to the T, P, and G fractions (p < 0.05). In addition, using the gene ontology application software AmiGO, the list of upregulated genes associated with both p38 MAPK signaling and cell motility was mapped. This resulted in a set of 9 genes, i.e. arachidonate 12-lipoxygenase (ALOX12), autocrine motility factor receptor (AMFR), Rho-guanine nucleotide exchange factor 2 (ARHGEF2), CCL5, engulfment and cell motility protein (ELMO1), signal transducer and activator of transcription 3 (STAT3), Talin-1 (TLN1), Tropomyosin alpha-4 chain (TPM4) and Vinculin (VCL).
Clinically prognostic value of cell motility gene signature
Median (range) fold change PDAC/control tissue (nCounter) for the validation set
10.62 – 15.29
10.90 – 13.07
9.99 – 12.56
7.67 – 11.79
9.73 – 11.54
12.67 – 15.64
0 (0.44) *
0 – 6.45
8.44 – 10.97
5.53 – 9.00
8.76 – 11.64
9.56 – 11.38
4.89 – 8.49
Univariable and multivariable Cox regression analysis for disease-free survival
HR (95% CI)
HR (95% CI)
1.628 (0.898 – 2.822)
0.557 (0.310 – 1.038)
0.597 (0.229 – 1.848)
0.431 (0.158 – 1.372)
1.207 (0.522 – 3.504)
0.899 (0.377 – 2.658)
2.023 (1.155 – 3.724)
2.089 (1.184 – 3.872)
3.454 (1.543 – 6.979)
2.978 (1.313 – 6.123)
≥4 genes or/and TGF-β1
1.903 (1.053 – 3.530)
1.885 (1.025 – 3.559)
Univariable and multivariable Cox regression analysis for overall survival
HR (95% CI)
HR (95% CI)
0.738 (0.232 - 3.267)
.718 (0.220 – 3.220)
1.508 (0.545 – 6.251)
1.191 (0.418 – 5.004)
2.042 (1.132 – 3.922)
1.658 (0.885 – 3.312)
2.986 (1.394 – 5.828)
2.715 (1.234 – 5.499)
1.974 (1.141 – 3.361)
1.870 (1.057 – 3.268)
1.986 (0.947 – 4.879)
1.680 (0.750 – 4.521)
≥ 6 genes or/and TGF-β1
1.734 (0.962 – 3.132)
1.366 (1.004 – 1.861)
On multivariable analysis, high co-expression of TGF-β1 and the cell motility panel (≥ 4 out of 9 genes for DFS and ≥ 6 out of 9 genes for OS) were identified as an independent predictor of both DFS (p=0.041, HR (95% CI) = 1.885 (1.025 – 3.559)) and OS (p=0.047, HR (95% CI) = 1.366 (1.004 – 1.861)) respectively (Table 4 and 5).
Identification of specific genetic and phenotypic changes occurring in CTC could result in a better understanding of the metastatic process and lead to effective therapeutic strategies . To our knowledge, for the first time global gene expression profiles were generated from blood enriched for CTC by a negative depletion technique. Depletion of dead (7AAD staining) and non-epithelial cell types (such as leukocytes (CD45+) and other endothelial/leukocyte progenitor cells (CD34+) results in a more unbiased enrichment of CTC both in spiking experiments and in vivo compared to a priori positive selection using epithelium-specific antibodies.
An important issue in small-sample gene expression microarray analyses is the quality of the final RNA sample. Indeed, all methods requiring additional steps (e.g. tissue biopsies from resection specimens, density centrifugation, FACS) to acquire very small amounts of cells bear the risk of RNA degradation. Moreover, the pancreas is known for its high concentration in RNases, leading to quick RNA degradation. In an effort to minimize this effect, periods of warm ischemia between removal of the surgical specimen and sample processing were kept to a strict minimum. Also, tissue samples were immediately submerged into a RNA stabilizing agent (RNAlater®) and care was taken to store these samples under optimal conditions. Despite these stringent precautions moderate degradation was observed in a substantial number of samples. Nevertheless, moderate RNA degradation probably disturb microarray results only modestly , and the final number of samples with high-quality RNA was sufficient to draw sound conclusions from the present study. Microarray analyses showed differential gene expression profiles between the CTC, T, P and G fractions. In contrast to other reports, we found only 6 genes (C19orf33, ECT2, IL1RN, S100P, SFN, TUBA4A) to be significantly upregulated in the primary PDAC (T) as compared to the non-tumoural pancreatic tissue (P) samples, whilst no genes were found to be downregulated. This finding can be explained by the fact that the ductal epithelial cells, of which PDAC originates, make up only a small proportion of the bulk surrounding pancreas or tumour. Indeed, expression profiles often predominantly reflect the proportion of cell types and obscure more meaningful data. Another explanation may be that we finally retained only 8,152 probe sets out of 46,467 probes initially tested, because Robust Multichip Average (RMA) was applied only to probe sets with a present detection call in at least 4 out of 6 CCD samples.
In our stringent comparison of the various cell fractions (CTC vs. T, AND CTC vs. P, AND CTC vs. G), the top molecular and cellular functions expressed in CTC were cell-to-cell signalling and interaction, cellular development, and cellular movement. The pathway with the highest expression ratio in CTC was p38 MAPK signalling, which is known to be involved in biological processes such as inflammation, apoptosis, and cell differentiation . More recently, the p38 MAPK signalling cascade was also shown to play a pivotal role in cancer cell migration, and could be stimulated by S100A8 . Moreover, potential secretion of TNFα, VEGF-A and TGF-β1 from primary tumours induced S100A8 and S100A9, which in turn stimulated the enrichment of other inflammatory chemo-attractants. These findings suggest that an increased p38 MAPK activity in CTC could occur in response to S100A8 produced not only in homing organs (such as lungs, liver) but also in CTC. Stimulation of p38 MAPK by S100A8 has been shown to induce formation of pseudopodia important to dynamic cell movement in the process of invasion and extravasation . Therefore, we used the gene ontology application software AmiGO at other genes that were overexpressed in CTC involved in the same process and in p38 MAPK signalling. Other genes mediating many of the component processes of invasion include Talin-1 and ELMO1. Both these genes were found to be upregulated in our CTC. Besides these genes, we found another 7 genes, i.e. ALOX12, AMFR, ARHGEF2, CCL5, STAT3, TPM4, and VCL. These 9 genes were termed CTC motility panel for ease of formulation.
In the p38 MAPK pathway TGF-β1, cPLA2, and MAX were significantly upregulated in CTC. TGF-β1 regulates a spectrum of cellular events, including cell proliferation, differentiation, and migration. In addition to the canonical Smad pathway, TGF-β1 can also activate MAPK, phosphatidylinositol 3-kinase (PI3K)/Akt and small GTPases in a cell-specific manner. Activation of ERK or p38 MAPK is required for both TGF-beta-induced epithelial-mesenchymal transition (EMT) and cell migration . Reversible and local activation of TGFβ/Smad-signalling in breast cancer cells has recently been shown to cause a switch from cohesive movement to single cells motility and to promote haematogenous metastasis .
A clear association has been described between the CCL5/CCR5 axis and the directional migration and invasion of human cancer cells. In vitro, this has been shown to involve specific upregulation of matrix-metalloproteinase-9 (MMP-9) . STAT3 is a transcription factor that has been linked with metastasis in different tumourtypes e.g. lung, colon adenocarcinoma, and PDAC . Besides increasing cell motility and invasion, STAT3 is also thought to upregulate the expression of anti-apoptotic and growth-promoting genes and to facilitate the colonisation of the liver by tumour cells [20, 21]. Therefore, STAT3 could be another new target for prevention of metastasis formation by CTC and for the treatment of established metastasis. STAT3 inhibitors are already available and currently undergo testing in various animal cancer models .
Our study provides a valuable source for researchers to explore. Using a negative depletion technique we were capable of isolating CTC, without using a priori selection markers. We have identified several candidate genes in CTC from PDAC, which can be used for various purposes. First, a combination of these genes may distinguish between benign and malignant disease of the pancreas. This was beyond the scope of this study. Also, as these genes represent the molecular signature of CTC in peripheral blood they may be useful for disease monitoring, prediction of survival, and response to therapy. Therefore we determined the expression level of our cell motility gene signature in 78 independent PDAC samples using the nCounter system. This system uses digital technology based on direct multiplexed measurement of gene expression, and offers high levels of sensitivity and precision . We found that the combination of expression of our CTC motility panel and TGF-β1 identified a subgroup of primary PDAC at high risk of early recurrence and worse overall survival. Some of these genes may also be interesting therapeutic targets and need further exploration.However, we should remain cautious as measurement of mRNA levels to infer events at the protein level can be misleading, as discordances may exist between mRNA levels and encoded protein levels.
Pancreatic CTC isolated from blood samples using a FACS-based negative depletion method, express a cell motility gene signature. The expression of this newly defined cell motility gene signature in the primary tumour enables us to predict survival of patients who undergo surgical resection for pancreatic cancer. The expression of this gene signature may be used as a selection tool in patients who are candidates for surgery, not only in pancreatic but also in other solid organ malignancies.
Pancreatic ductal adenocarcinoma
Pancreatic tumoural tissue
Benign surrounding pancreatic tissue
Circulating tumour cells
Haematological and non-viable background.
GS acknowledges support by a doctoral grant from the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT Vlaanderen).
BT acknowledges support by a research grant from the Fund for Scientific Research - Flanders (FWO-Vlaanderen).
Special thanks to François Houtmeyers, Hugo Vankelecom, Joke Allemeersch, Nadine Ectors, Paul Van Hummelen, Jos van Pelt, Petra Windmolders, Raymond Aerts, and Ruth Maes for their technical advice and assistance in sample collection and preparation.
- Hanahan D, Weinberg RA: The hallmarks of cancer. Cell. 2000, 100: 57-70. 10.1016/S0092-8674(00)81683-9.View ArticlePubMed
- Pantel K, Riethdorf S: Pathology: are circulating tumour cells predictive of overall survival?. Nat Rev Clin Oncol. 2009, 6: 190-191. 10.1038/nrclinonc.2009.23.View ArticlePubMed
- Gertler R, Rosenberg R, Fuehrer K, Dahm M, Nekarda H, Siewert JR: Detection of circulating tumour cells in blood using an optimized density gradient centrifugation. Recent Results Cancer Res. 2003, 162: 149-155. 10.1007/978-3-642-59349-9_13.View ArticlePubMed
- Rosenberg R, Gertler R, Friederichs J, Fuehrer K, Dahm M, Phelps R, Thorban S, Nekarda H, Siewert JR: Comparison of two density gradient centrifugation systems for the enrichment of disseminated tumour cells in blood. Cytometry. 2002, 49: 150-158. 10.1002/cyto.10161.View ArticlePubMed
- Lieber M, Mazzetta J, Nelson-Rees W, Kaplan M, Todaro G: Establishment of a continuous tumour-cell line (panc-1) from a human carcinoma of the exocrine pancreas. Int J Cancer. 1975, 15: 741-747. 10.1002/ijc.2910150505.View ArticlePubMed
- Gautier L, Cope L, Bolstad BM, Irizarry RA: affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315. 10.1093/bioinformatics/btg405.View ArticlePubMed
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.PubMed CentralView ArticlePubMed
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-PubMed
- The R Project for Statistical Computing. http://www.r-project.org/,
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995, 57: 289-300.
- Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S: AmiGO: online access to ontology and annotation data. Bioinformatics. 2009, 25: 288-289. 10.1093/bioinformatics/btn615.PubMed CentralView ArticlePubMed
- Albertson DG, Collins C, McCormick F, Gray JW: Chromosome aberrations in solid tumours. Nat Genet. 2003, 34: 369-376. 10.1038/ng1215.View ArticlePubMed
- Schoor O, Weinschenk T, Hennenlotter J, Corvin S, Stenzl A, Rammensee HG, Stevanovic S: Moderate degradation does not preclude microarray analysis of small amounts of RNA. Biotechniques. 2003, 35: 1192-1196. 1198–1201PubMed
- Hiratsuka S, Watanabe A, Aburatani H, Maru Y: Tumour-mediated upregulation of chemoattractants and recruitment of myeloid cells predetermines lung metastasis. Nat Cell Biol. 2006, 8: 1369-1375. 10.1038/ncb1507.View ArticlePubMed
- Ehrchen JM, Sunderkotter C, Foell D, Vogl T, Roth J: The endogenous Toll-like receptor 4 agonist S100A8/S100A9 (calprotectin) as innate amplifier of infection, autoimmunity, and cancer. J Leukoc Biol. 2009, 86: 557-566. 10.1189/jlb.1008647.View ArticlePubMed
- Zuo W, Chen YG: Specific activation of mitogen-activated protein kinase by transforming growth factor-beta receptors in lipid rafts is required for epithelial cell plasticity. Mol Biol Cell. 2009, 20: 1020-1029.PubMed CentralView ArticlePubMed
- Giampieri S, Manning C, Hooper S, Jones L, Hill CS, Sahai E: Localized and reversible TGFbeta signalling switches breast cancer cells from cohesive to single cell motility. Nat Cell Biol. 2009, 11: 1287-1296. 10.1038/ncb1973.PubMed CentralView ArticlePubMed
- Chuang JY, Yang WH, Chen HT, Huang CY, Tan TW, Lin YT, Hsu CJ, Fong YC, Tang CH: CCL5/CCR5 axis promotes the motility of human oral cancer cells. J Cell Physiol. 2009, 220: 418-426. 10.1002/jcp.21783.View ArticlePubMed
- Yu H, Jove R: The STATs of cancer–new molecular targets come of age. Nat Rev Cancer. 2004, 4: 97-105. 10.1038/nrc1275.View ArticlePubMed
- Aggarwal BB, Kunnumakkara AB, Harikumar KB, Gupta SR, Tharakan ST, Koca C, Dey S, Sung B: Signal transducer and activator of transcription-3, inflammation, and cancer: how intimate is the relationship?. Ann N Y Acad Sci. 2009, 1171: 59-76. 10.1111/j.1749-6632.2009.04911.x.PubMed CentralView ArticlePubMed
- Steeg PS: Tumour metastasis: mechanistic insights and clinical challenges. Nat Med. 2006, 12: 895-904. 10.1038/nm1469.View ArticlePubMed
- Geiss GK, Bumgarner RE, Birditt B, Dahl T, Dowidar N, Dunaway DL, Fell HP, Ferree S, George RD, Grogan T, James JJ, Maysuria M, Mitton JD, Oliveri P, Osborn JL, Peng T, Ratcliffe AL, Webster PJ, Davidson EH, Hood L, Dimitrov K: Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat Biotechnol. 2008, 26: 317-325. 10.1038/nbt1385.View ArticlePubMed
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/12/527/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.