Stromal gene expression defines poor-prognosis subtypes in colorectal cancer (original) (raw)

Accession codes

Primary accessions

European Nucleotide Archive

Gene Expression Omnibus

Referenced accessions

Gene Expression Omnibus

References

1. De Sousa E Melo, F. et al. Poor-prognosis colon cancer is defined by a molecularly distinct subtype and develops from serrated precursor lesions. Nat. Med. 19, 614–618 (2013).
   Article CAS PubMed Google Scholar
2. Sadanandam, A. et al. A colorectal cancer classification system that associates cellular phenotype and responses to therapy. Nat. Med. 19, 619–625 (2013).
   Article CAS PubMed PubMed Central Google Scholar
3. Marisa, L. et al. Gene expression classification of colon cancer into molecular subtypes: characterization, validation, and prognostic value. PLoS Med. 10, e1001453 (2013).
   Article CAS PubMed PubMed Central Google Scholar
4. Sadanandam, A. et al. Reconciliation of classification systems defining molecular subtypes of colorectal cancer: interrelationships and clinical implications. Cell Cycle 13, 353–357 (2014).
   Article CAS PubMed PubMed Central Google Scholar
5. Nishida, N. et al. Microarray analysis of colorectal cancer stromal tissue reveals upregulation of two oncogenic miRNA clusters. Clin. Cancer Res. 18, 3054–3070 (2012).
   Article CAS PubMed Google Scholar
6. Calon, A. et al. Dependency of colorectal cancer on a TGF-β–driven program in stromal cells for metastasis initiation. Cancer Cell 22, 571–584 (2012).
   Article CAS PubMed PubMed Central Google Scholar
7. Uhlen, M. et al. Towards a knowledge-based Human Protein Atlas. Nat. Biotechnol. 28, 1248–1250 (2010).
   Article CAS PubMed Google Scholar
8. Mouradov, D. et al. Colorectal cancer cell lines are representative models of the main molecular subtypes of primary cancer. Cancer Res. 74, 3238–3247 (2014).
   Article CAS PubMed Google Scholar
9. O'Brien, C.A., Pollett, A., Gallinger, S. & Dick, J.E. A human colon cancer cell capable of initiating tumour growth in immunodeficient mice. Nature 445, 106–110 (2007).
   Article CAS PubMed Google Scholar
10. Ricci-Vitiani, L. et al. Identification and expansion of human colon-cancer-initiating cells. Nature 445, 111–115 (2007).
    Article CAS PubMed Google Scholar
11. Oskarsson, T., Batlle, E. & Massague, J. Metastatic stem cells: sources, niches, and vital pathways. Cell Stem Cell 14, 306–321 (2014).
    Article CAS PubMed PubMed Central Google Scholar
12. Thiery, J.P., Acloque, H., Huang, R.Y. & Nieto, M.A. Epithelial-mesenchymal transitions in development and disease. Cell 139, 871–890 (2009).
    Article CAS PubMed Google Scholar
13. Markowitz, S.D. & Bertagnolli, M.M. Molecular origins of cancer: molecular basis of colorectal cancer. N. Engl. J. Med. 361, 2449–2460 (2009).
    Article CAS PubMed PubMed Central Google Scholar
14. Grady, W.M. & Markowitz, S.D. Genetic and epigenetic alterations in colon cancer. Annu. Rev. Genomics Hum. Genet. 3, 101–128 (2002).
    Article CAS PubMed Google Scholar
15. Markowitz, S. et al. Inactivation of the type II TGF-β receptor in colon cancer cells with microsatellite instability. Science 268, 1336–1338 (1995).
    Article CAS PubMed Google Scholar
16. Jung, P. et al. Isolation and in vitro expansion of human colonic stem cells. Nat. Med. 17, 1225–1227 (2011).
    Article CAS PubMed Google Scholar
17. Sato, T. et al. Long-term expansion of epithelial organoids from human colon, adenoma, adenocarcinoma, and Barrett's epithelium. Gastroenterology 141, 1762–1772 (2011).
    Article CAS PubMed Google Scholar
18. Bueno, L. et al. Semi-mechanistic modelling of the tumour growth inhibitory effects of LY2157299, a new type I receptor TGF-β kinase antagonist, in mice. Eur. J. Cancer 44, 142–150 (2008).
    Article CAS PubMed Google Scholar
19. Seoane, J., Le, H.V. & Massague, J. Myc suppression of the p21Cip1 Cdk inhibitor influences the outcome of the p53 response to DNA damage. Nature 419, 729–734 (2002).
    Article CAS PubMed Google Scholar
20. Seoane, J. et al. TGFβ influences Myc, Miz-1 and Smad to control the CDK inhibitor p15INK4b. Nat. Cell Biol. 3, 400–408 (2001).
    Article CAS PubMed Google Scholar
21. Giannelli, G., Villa, E. & Lahn, M. Transforming growth factor-β as a therapeutic target in hepatocellular carcinoma. Cancer Res. 74, 1890–1894 (2014).
    Article CAS PubMed Google Scholar
22. Kosinski, C. et al. Gene expression patterns of human colon tops and basal crypts and BMP antagonists as intestinal stem cell niche factors. Proc. Natl. Acad. Sci. USA 104, 15418–15423 (2007).
    Article CAS PubMed PubMed Central Google Scholar
23. Berdiel-Acer, M. et al. A 5-gene classifier from the carcinoma-associated fibroblast transcriptomic profile and clinical outcome in colorectal cancer. Oncotarget. 5, 6437–6452 (2014).
    Article PubMed PubMed Central Google Scholar
24. O'Shannessy, D.J. et al. Influence of tumor microenvironment on prognosis in colorectal cancer: tissue architecture–dependent signature of endosialin (TEM-1) and associated proteins. Oncotarget. 5, 3983–3995 (2014).
    PubMed PubMed Central Google Scholar
25. Francí, C. et al. Snail1 protein in the stroma as a new putative prognosis marker for colon tumours. PLoS ONE 4, e5595 (2009).
    Article PubMed PubMed Central Google Scholar
26. Ngan, C.Y. et al. Quantitative evaluation of vimentin expression in tumour stroma of colorectal cancer. Br. J. Cancer 96, 986–992 (2007).
    Article CAS PubMed PubMed Central Google Scholar
27. Calon, A., Tauriello, D.V. & Batlle, E. TGF-β in CAF-mediated tumor growth and metastasis. Semin. Cancer Biol. 25, 15–22 (2014).
    Article CAS PubMed Google Scholar
28. Barrett, T. & Edgar, R. Gene expression omnibus: microarray data storage, submission, retrieval, and analysis. Methods Enzymol. 411, 352–369 (2006).
    Article CAS PubMed PubMed Central Google Scholar
29. Jorissen, R.N. et al. Metastasis-associated gene expression changes predict poor outcomes in patients with Dukes stage B and C colorectal cancer. Clin. Cancer Res. 15, 7642–7651 (2009).
    Article CAS PubMed PubMed Central Google Scholar
30. De Sousa E Melo, F. et al. Methylation of cancer-stem-cell-associated Wnt target genes predicts poor prognosis in colorectal cancer patients. Cell Stem Cell 9, 476–485 (2011).
    Article CAS PubMed Google Scholar
31. Benjamini, Y., Krieger, A.M. & Yekutieli, D. Adaptive linear step-up procedures that control the false discovery rate. Biometrika 93, 491–507 (2006).
    Article Google Scholar
32. Tibshirani, R., Hastie, T., Narasimhan, B. & Chu, G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc. Natl. Acad. Sci. USA 99, 6567–6572 (2002).
    Article CAS PubMed PubMed Central Google Scholar
33. Merlos-Suárez, A. et al. The intestinal stem cell signature identifies colorectal cancer stem cells and predicts disease relapse. Cell Stem Cell 8, 511–524 (2011).
    Article PubMed Google Scholar
34. Lombardo, Y. et al. Bone morphogenetic protein 4 induces differentiation of colorectal cancer stem cells and increases their response to chemotherapy in mice. Gastroenterology 140, 297–309 (2011).
    Article CAS PubMed Google Scholar
35. Morikawa, K., Walker, S.M., Jessup, J.M. & Fidler, I.J. In vivo selection of highly metastatic cells from surgical specimens of different primary human colon carcinomas implanted into nude mice. Cancer Res. 48, 1943–1948 (1988).
    CAS PubMed Google Scholar
36. Morikawa, K. et al. Influence of organ environment on the growth, selection, and metastasis of human colon carcinoma cells in nude mice. Cancer Res. 48, 6863–6871 (1988).
    CAS PubMed Google Scholar
37. Böttinger, E.P., Jakubczak, J.L., Haines, D.C., Bagnall, K. & Wakefield, L.M. Transgenic mice overexpressing a dominant-negative mutant type II transforming growth factor β receptor show enhanced tumorigenesis in the mammary gland and lung in response to the carcinogen 7,12-dimethylbenz-[_a_]-anthracene. Cancer Res. 57, 5564–5570 (1997).
    PubMed Google Scholar
38. Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 102, 15545–15550 (2005).
    Article CAS PubMed PubMed Central Google Scholar
39. Witten, I.H. et al. Weka: practical machine learning tools and techniques with Java implementations. Proceedings of the ICONIP/ANZIIS/ANNES'99 Workshop on Emerging Knowledge Engineering and Connectionist-Based Information Systems 192–196 (1999).
40. Warren, R.S., Yuan, H., Matli, M.R., Gillett, N.A. & Ferrara, N. Regulation by vascular endothelial growth factor of human colon cancer tumorigenesis in a mouse model of experimental liver metastasis. J. Clin. Invest. 95, 1789–1797 (1995).
    Article CAS PubMed PubMed Central Google Scholar
41. Hu, Y. & Smyth, G.K. ELDA: extreme limiting dilution analysis for comparing depleted and enriched populations in stem cell and other assays. J. Immunol. Methods 347, 70–78 (2009).
    Article CAS PubMed Google Scholar
42. Lander, E.S. et al. Initial sequencing and analysis of the human genome. Nature 409, 860–921 (2001).
    Article CAS PubMed Google Scholar
43. Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
    Article PubMed PubMed Central Google Scholar
44. O'Brien, C.A. et al. ID1 and ID3 regulate the self-renewal capacity of human colon cancer–initiating cells through p21. Cancer Cell 21, 777–792 (2012).
    Article CAS PubMed Google Scholar
45. McKenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).
    Article CAS PubMed PubMed Central Google Scholar
46. DePristo, M.A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498 (2011).
    Article CAS PubMed PubMed Central Google Scholar
47. Pathmanathan, S., Hamilton, E., Atcheson, E. & Timson, D.J. The interaction of IQGAPs with calmodulin-like proteins. Biochem. Soc. Trans. 39, 694–699 (2011).
    Article CAS PubMed Google Scholar
48. Clevers, H. & Batlle, E. EphB/EphrinB receptors and Wnt signaling in colorectal cancer. Cancer Res. 66, 2–5 (2006).
    Article CAS PubMed Google Scholar
49. Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin) 6, 80–92 (2012).
    Article CAS Google Scholar
50. R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2014).
51. Kent, W.J. et al. The human genome browser at UCSC. Genome Res. 12, 996–1006 (2002).
    CAS PubMed PubMed Central Google Scholar
52. Lawrence, M., Gentleman, R. & Carey, V. rtracklayer: an R package for interfacing with genome browsers. Bioinformatics 25, 1841–1842 (2009).
    Article CAS PubMed PubMed Central Google Scholar
53. Wood, S.N. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J. Royal Statist. Soc. (B) 73, 3–36 (2011).
    Article Google Scholar
54. Venkatraman, E.S. & Olshen, A.B. A faster circular binary segmentation algorithm for the analysis of array CGH data. Bioinformatics 23, 657–663 (2007).
    Article CAS PubMed Google Scholar

Download references

Acknowledgements

We thank G. Stassi (University of Palermo) for providing PDO7 and PDO8, L. Wakefield (US National Cancer Institute) for providing the plasmid encoding DNR, I. Joval for assistance in mounting the figures, M. Virtudes Cespedes and R. Mangues (IIB Sant Pau) for logistic support with CRC samples, and all members of the Batlle laboratory for support and discussions We are grateful for the excellent assistance of the IRB Barcelona core facilities for Histology, Functional Genomics and Advanced Digital Micropscopy. D.V.F.T. holds a Juan de la Cierva postdoctoral fellowship, from the Spanish Ministry of Economy and Competitiveness, and E.L. holds a fellowship from Fundación Olga Torres and Asociación Española contra el Cáncer (AECC). This work has been supported by grants from the Doctor Josef Steiner Foundation, AECC, Red Temática de Investigación Cooperativa en Cáncer, Instituto de Salud Carlos III (RTICC:RD12/0036/0024) and grant SAF2011-27068, the latter two from the Spanish Ministry of Economy and Competitiveness, and by ‘Xarxa de Bancs de Tumors’ sponsored by Pla Director d'Oncologia de Catalunya (XBTC).

Author information

Author notes

1. Elisa Espinet & David Rossell
   Present address: Present addresses: Division of Stem Cells and Cancer, German Cancer Research Center (DKFZ), Heidelberg, Germany (E.E.) and Department of Statistics, University of Warwick, Warwick, UK (D.R.).,

Authors and Affiliations

1. Institute for Research in Biomedicine (IRB) Barcelona, Barcelona, Spain
   Alexandre Calon, Enza Lonardo, Antonio Berenguer-Llergo, Elisa Espinet, Xavier Hernando-Momblona, Marta Sevillano, Sergio Palomo-Ponce, Daniele V F Tauriello, Daniel Byrom, Carme Cortina, Clara Morral, Carles Barceló, Sebastien Tosi, Antoni Riera, Camille Stephan-Otto Attolini, David Rossell, Elena Sancho & Eduard Batlle
2. Department of Pathology, Hospital del Mar, Barcelona, Spain
   Mar Iglesias
3. Cancer Research Program, Hospital del Mar Research Institute (IMIM), Barcelona, Spain
   Mar Iglesias
4. Universitat Autónoma de Barcelona, Barcelona, Spain
   Mar Iglesias
5. Institució Catalana de Recerca i Estudis Avançats (iCREA), Barcelona, Spain
   Eduard Batlle

Authors

1. Alexandre Calon
   You can also search for this author inPubMed Google Scholar
2. Enza Lonardo
   You can also search for this author inPubMed Google Scholar
3. Antonio Berenguer-Llergo
   You can also search for this author inPubMed Google Scholar
4. Elisa Espinet
   You can also search for this author inPubMed Google Scholar
5. Xavier Hernando-Momblona
   You can also search for this author inPubMed Google Scholar
6. Mar Iglesias
   You can also search for this author inPubMed Google Scholar
7. Marta Sevillano
   You can also search for this author inPubMed Google Scholar
8. Sergio Palomo-Ponce
   You can also search for this author inPubMed Google Scholar
9. Daniele V F Tauriello
   You can also search for this author inPubMed Google Scholar
10. Daniel Byrom
    You can also search for this author inPubMed Google Scholar
11. Carme Cortina
    You can also search for this author inPubMed Google Scholar
12. Clara Morral
    You can also search for this author inPubMed Google Scholar
13. Carles Barceló
    You can also search for this author inPubMed Google Scholar
14. Sebastien Tosi
    You can also search for this author inPubMed Google Scholar
15. Antoni Riera
    You can also search for this author inPubMed Google Scholar
16. Camille Stephan-Otto Attolini
    You can also search for this author inPubMed Google Scholar
17. David Rossell
    You can also search for this author inPubMed Google Scholar
18. Elena Sancho
    You can also search for this author inPubMed Google Scholar
19. Eduard Batlle
    You can also search for this author inPubMed Google Scholar

Contributions

A.C., E.L. and E.E. designed, planned and performed experiments and analyzed the results. X.H.-M. and S.P.-P. provided crucial assistance with in vivo experiments. M.S. performed immunohistochemistry. M.I. scored the tumor microarray. S.T. developed an algorithm to quantify xenograft images. D.V.F.T., C.C., C.B. and C.M. performed experiments and/or analyzed results. D.V.F.T., D.B. and A.R. synthesized the LY2157299 inhibitor. A.B.-L., C.S.-O.A. and D.R. designed and performed biostatistical analyses. E.B. conceptualized and supervised the project, analyzed results and wrote the manuscript, with the assistance of E.S.

Corresponding author

Correspondence toEduard Batlle.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Integrated supplementary information

Supplementary Figure 1 High levels of genes characteristic of CAFs identify poor-prognosis patients across multiple subtypes in the De Sousa E Melo et al. molecular classification.

(a) Clustering analysis of the 146-probeset signature used to classify CRC patients into subtypes by De Sousa E Melo et al. in the GSE14333 cohort. We allowed unsupervised hierarchical clustering of the probesets using correlation distance and the Ward agglomeration method, whereas we enforced the classification of patients into subtypes. Data show normalized, centered and scaled Affymetrix probeset intensities on a log2 scale. The HR lane represents the hazard ratio for the corresponding genes on a log2 scale. Genes significantly upregulated in the microdissected epithelial or stromal compartment are depicted in the epithelial (LCM) and stroma (LCM) lanes, respectively. Genes specifically upregulated in epithelial, endothelial, leukocyte or FAP cell populations are represented in the EpCAM, CD31, CD45 and FAP lanes, respectively. Patients with high average expression (z scores >0) of the CAF-enriched gene cluster are marked in green (GP-HiC). (b) The smooth estimate of HR (+1 s.d.) (increase in recurrence risk for one standard deviation) shows higher risk of relapse for patients within the good-prognosis subtypes (CCS1 + CCS2) presenting higher average expression of the CAF cluster (GP-HiC group). Red dashed lines indicate 95% confidence bands. HR values for the linear continuous model and corresponding P values are indicated. (c) The Kaplan-Meier curve shows the recurrence-free survival of patients within good-prognosis cancer subtypes and presenting low expression levels of the CAF cluster gene set (blue; GP-LoC) or patients within good-prognosis cancer subtypes that present high expression levels of the CAF cluster gene set (yellow; GP-HiC), both compared to the stem cell–like poor-prognosis subtype (red; CCS3). HR and P values are indicated.

Supplementary Figure 2 High levels of genes characteristic of CAFs in C4 poor-prognosis patients in the Marisa et al. molecular classification.

(a) Clustering analysis of the 1,459-probeset signature used to classify CRC patients into subtypes by Marisa et al. in the GSE39582 cohort. We allowed unsupervised hierarchical clustering of the probesets using correlation distance and the Ward agglomeration method, whereas we enforced the classification of patients into subtypes. Data show normalized, centered and scaled Affymetrix probeset intensities on a log2 scale. The HR lane represents the hazard ratio for the corresponding genes on a log2 scale. Genes significantly upregulated in the microdissected epithelial or stromal compartment are depicted in the epithelial (LCM) and stroma (LCM) lanes, respectively. Genes specifically upregulated in epithelial, endothelial, leukocyte or FAP cell populations are represented in the EpCAM, CD31, CD45 and FAP lanes, respectively. Patients with high average expression (z scores >0) of the CAF-enriched gene cluster are marked in green (GP-HiC). (b) The smooth estimate of HR (+1 s.d.) (increase in recurrence risk for one standard deviation) shows higher risk of relapse for patients within the good-prognosis subtypes (C1 + C2 + C3 + C5 + C6) presenting higher average expression of the CAF cluster (GP-HiC group). Red dashed lines indicate 95% confidence bands. HR values for the linear continuous model and corresponding P values are indicated. (c) The Kaplan-Meier curve shows the recurrence-free survival of patients within good-prognosis cancer subtypes and presenting low expression levels of the CAF cluster gene set (blue; GP-LoC) or patients within good-prognosis cancer subtypes that present high expression levels of the CAF cluster gene set (yellow; GP-HiC), both compared to the C4 poor-prognosis subtype (red). HR and P values are indicated.

Supplementary Figure 3 Expression patterns of markers of poor prognosis.

(a) The pie chart shows the expression pattern distribution (percentage) of clinically relevant antibodies in individuals with CRC within the Human Protein Atlas database. Stromal-specific antibodies account for approximately 31% of this distribution, including the newly identified CALD1 and POSTN proteins (red). Antibodies to FAP and IGFBP7 belong to the 62% of antibodies that stain both the stromal and epithelial compartments (yellow). We found that 3% of antibodies are epithelial specific (blue), whereas the rest were not detectable in either fraction (4%; gray slice). (b) _z_-score mean of CALD1, FAP, POSTN and IGFBP7 in the cell populations from disaggregated primary CRC samples. ***P < 0.001.

Supplementary Figure 4 Poor-prognosis CRC subtypes are characterized by high levels of TGF-β and TBRS expression.

(a) _z_-score means of TGFB1, TGFB2 and TGFB3 mRNA expression levels in the three GEO data sets in each molecular subtype. The poor-prognosis group (red box) for each data set contains the highest levels of TGFB1 and TGFB3. (b,c) _z_-score means of TBRSs (TGF-β–responsive signatures) in the three data sets in each molecular subtype. The TBRSs are extracted from (b) macrophages (Ma-TBRS) or (c) lymphocytes (T-TBRS) treated with recombinant TGF-β1. The poor-prognosis groups for each data set are depicted in red. P values for pairwise comparisons are in Supplementary Table 7.

Supplementary Figure 5 TGF-β–activated fibroblasts promote tumor initiation.

Colony formation assay (a) and qRT-PCR analysis of cell cycle genes (b) in HT29-M6 and KM12L4a cells treated for 7 d with recombinant TGF-β1 (red bar) compared to untreated (blue bar) CRC cells. qRT-PCR data are normalized to PP1A levels and are presented as fold change versus the control (Con) (n ≥ 3 independent experiments). (c) The Kaplan-Meier curves show the disease-free survival (DFS) of mice after subcutaneous injections of 6 × 103 KM12L4aTGFβ cells alone (n = 8 tumor cell inoculations; black line) or coinoculated with FIB-shCon (n = 8 tumor cell inoculations; red line) or FIB-shTBRI (n = 8 tumor cell inoculations; blue line). FIB, fibroblasts (5 × 104 cells). *P < 0.05. (d) Growth kinetics of HT29-M6 (blue; n = 4 tumor cell inoculations) and HT29-M6TGFβ (red; n = 4 tumor cell inoculations) cells inoculated subcutaneously alone or in combination with FIB (green and yellow lines, respectively). Day 1, day of first detection. Values are means ± s.e.m.

Supplementary Figure 6 Coinjection of CRC cells with wild-type or TGF-β pathway–defective fibroblasts did not modify xenograft growth rates.

(a) Growth kinetics of HT29-M6TGFβ cells coinoculated subcutaneously with FIB-shCon (n = 8 tumor cell inoculations; blue) or FIB-shTBRI (n = 5 tumor cell inoculations; red). (b) Growth kinetics of HT29-M6TGFβ cells coinoculated subcutaneously with FIB-Con (n = 7 tumor cell inoculations; blue) or FIB-DNR (n = 5 tumor cell inoculations; red). Day 1, day of first detection. Values are means ± s.e.m. (c,d) Histology of xenografts with HT29-M6 and HT29-M6TGFβ cells inoculated subcutaneously alone (Con) or in combination with FIB, FIB-shTBRI or FIB-DNR. (Scale bars, 200 μm). (c) We use trichrome Masson staining to mark epithelial cells in red and collagen fibers in blue. (d) EpCAM immunostaining marks epithelial tumor cells. The histological patterns were equivalent in all conditions. (e) We measured the surfaces occupied by stromal and epithelial cells using trichrome Masson staining (Online Methods) in multiple tumors formed after subcutaneous inoculation of HT29-M6TGFβ cells alone or in combination with FIB or FIB-DNR cells. This analysis confirmed no significant differences between the three experimental conditions.

Supplementary Figure 7 Fibroblasts lacking a TGF-β response showed minor capacity to support tumor formation.

(a) Kaplan-Meier curves showing the disease-free survival (DFS) of mice after subcutaneous inoculation of HT29-M6 cells alone (black line; n = 6 tumor cell inoculations) or coinjection with FIB-shCon (red line; n = 6 tumor cell inoculations) or FIB-shTBRI (blue line; n = 6 tumo cells inoculations) cells. FIB, fibroblasts (5 × 104 cells). (b) Growth kinetics of FIB cultured in the presence (red line) or absence (blue line) of recombinant TGF-β1. Cells were treated with TGF-β1 from day 1. Growth was followed for 6 d, and the number of live cells was monitored at the indicated time points. Values are means ± s.d. **P < 0.01. (c) The table shows raw data for the subcutaneous inoculation experiments in the indicated conditions. The columns tumors/n indicate the number of tumor initiation events per number of inoculations. The percentage (%) of successful tumor initiation events and DFS (days ± s.d.) are indicated for each condition. (d) Estimated frequency of tumor-initiating cells (TICs) calculated using ELDA (Online Methods) from the above data (***P < 0.001 for comparisons against TGF-β + FIB); 95% confidence intervals are indicated.

Supplementary Figure 8 Heterogeneity of patient-derived tumor organoids (PDO) growth under TGF-β treatment.

Growth kinetics of tumor organoids cultured in the presence or absence of recombinant TGF-β1 (red line) or LY2157299 (yellow line). The blue line represents untreated cells. Growth was followed for 7 days (D), and tumor organoids were counted at the indicated time points. Values are means ± s.d.

Supplementary Figure 9 TGF-β inhibition impairs tumor initiation.

(a) Kaplan-Meier curves show disease-free survival (DFS) for mice injected subcutaneously with cells from PDO1, 2, 3, 4, 5, 6, 7 or 8, treated (red line; n = 12 tumor cell inoculations) or untreated (blue line; n = 12 tumor cell inoculations) with LY2157299. *P < 0.05; ***P < 0.001. (b) Tumor volume over time for the indicated tumor organoids (from a), measured from the day of appearance, treated (red) or untreated (blue) with LY2157299. Values are means ± s.e.m. (c) Kaplan-Meier curves and growth curves for mice injected subcutaneously with cells from PDO5 induced (ind.) with TGF-β1, treated (yellow) or untreated (blue) with doxycycline (4 × 105 cells; n = 18 tumor cell inoculations). Values are means ± s.e.m.

Supplementary information

Supplementary Text and Figures

Supplementary Figures 1–9 and Supplementary Note. (PDF 1934 kb)

Supplementary Table 1

Genes positively associated to relapse (HR > 1, P < 0.01) and upregulated in poor-prognosis subtypes in each cohort (probe level data). (XLS 1455 kb)

Supplementary Table 2

HR and stromal enrichment of poor-prognosis gene sets across multiple thresholds of gene expression and cutoffs of significance. (XLS 39 kb)

Supplementary Table 3

Reclassification of CRC patients according to the EpCAM+ cluster or excluding the stromal cluster. (XLS 45 kb)

Supplementary Table 4

Expression patterns of genes commonly associated with poor prognosis in different patient cohorts according to the Human Protein Atlas database. (XLS 41 kb)

Supplementary Table 5

Cox proportional hazards multivariate analysis of stromal protein intensities, adjusted by AJCC stage, sex, age, tumor site and treatment. (XLS 30 kb)

Supplementary Table 6

Statistical analysis of IHC data (adjusted for TMA, stage, sex, age, tumor site and treatment). (XLS 21 kb)

Supplementary Table 7

Statistics for group-to-group comparisons for TGF-β and FTBRS average levels in different CRC patient data sets. (XLS 35 kb)

Supplementary Table 8

TGF-β response signatures (TBRSs) in normal stromal cells. (XLS 278 kb)

Supplementary Table 9

TBRS predicts disease relapse in the three CRC patient cohorts independently of the main clinical parameters. (XLS 28 kb)

Supplementary Table 10

Clinicopathological features and summary of genetic alterations in major pathways according to exome sequencing in tumor organoids. (XLS 29 kb)

Rights and permissions

About this article

Cite this article

Calon, A., Lonardo, E., Berenguer-Llergo, A. et al. Stromal gene expression defines poor-prognosis subtypes in colorectal cancer.Nat Genet 47, 320–329 (2015). https://doi.org/10.1038/ng.3225

Download citation

• Received: 24 September 2014
• Accepted: 28 January 2015
• Published: 23 February 2015
• Issue Date: April 2015
• DOI: https://doi.org/10.1038/ng.3225