The site of breast cancer metastases dictates their clonal composition and reversible transcriptomic profile

Optical barcoding reveals enhanced TNFα activation and polyclonality in lung, but not liver, metastases from breast cancer. Intratumoral heterogeneity is a driver of breast cancer progression, but the nature of the clonal interactive network involved in this process remains unclear. Here, we optimized the use of optical barcoding to visualize and characterize 31 cancer subclones in vivo. By mapping the clonal composition of thousands of metastases in two clinically relevant sites, the lungs and liver, we found that metastases were highly polyclonal in lungs but not in the liver. Furthermore, the transcriptome of the subclones varied according to their metastatic niche. We also identified a reversible niche-driven signature that was conserved in lung and liver metastases collected during patient autopsies. Among this signature, we found that the tumor necrosis factor–α pathway was up-regulated in lung compared to liver metastases, and inhibition of this pathway affected metastasis diversity. These results highlight that the cellular and molecular heterogeneity observed in metastases is largely dictated by the tumor microenvironment.

[1]  G. Christofori,et al.  Multi-color clonal tracking reveals intra-stage proliferative heterogeneity during mammary tumor progression , 2020, Oncogene.

[2]  Matthew E. Ritchie,et al.  Targeting triple-negative breast cancers with the Smac-mimetic birinapant , 2020, Cell Death & Differentiation.

[3]  K. Polyak,et al.  Intratumor Heterogeneity: The Rosetta Stone of Therapy Resistance. , 2020, Cancer cell.

[4]  E. Mjolsness,et al.  Transcriptional diversity and bioenergetic shift in human breast cancer metastasis revealed by single-cell RNA sequencing , 2020, Nature Cell Biology.

[5]  C. Curtis,et al.  Multi-cancer analysis of clonality and the timing of systemic spread in paired primary tumors and metastases , 2019, Nature Genetics.

[6]  C. L. Chaffer,et al.  The Complexities of Metastasis , 2019, Cancers.

[7]  A. Teschendorff,et al.  Single-cell landscape in mammary epithelium reveals bipotent-like cells associated with breast cancer risk and outcome , 2019, Communications Biology.

[8]  N. Meyer,et al.  The TNF Paradox in Cancer Progression and Immunotherapy , 2019, Front. Immunol..

[9]  O. Cohen,et al.  Subclonal cooperation drives metastasis by modulating local and systemic immune microenvironments , 2019, Nature Cell Biology.

[10]  S. Tabrizi,et al.  Inhibition of tumour necrosis factor alpha in the R6/2 mouse model of Huntington’s disease by etanercept treatment , 2019, Scientific Reports.

[11]  Maria C. Lecca,et al.  Spatiotemporal regulation of clonogenicity in colorectal cancer xenografts , 2019, Proceedings of the National Academy of Sciences.

[12]  Daniel L. Cameron,et al.  Barcoding reveals complex clonal behavior in patient-derived xenografts of metastatic triple negative breast cancer , 2019, Nature Communications.

[13]  Jeffrey T. Chang,et al.  High-resolution clonal mapping of multi-organ metastasis in triple negative breast cancer , 2018, Nature Communications.

[14]  Christoph Hafemeister,et al.  Comprehensive integration of single cell data , 2018, bioRxiv.

[15]  V. Prassolov,et al.  Cellular Barcoding Identifies Clonal Substitution as a Hallmark of Local Recurrence in a Surgical Model of Head and Neck Squamous Cell Carcinoma. , 2018, Cell reports.

[16]  P. Houghton,et al.  IL-6 and CXCL8 mediate osteosarcoma-lung interactions critical to metastasis. , 2018, JCI insight.

[17]  Joshua M. Dempster,et al.  Genetic and transcriptional evolution alters cancer cell line drug response , 2018, Nature.

[18]  Hong Wang,et al.  Shared-nearest-neighbor-based clustering by fast search and find of density peaks , 2018, Inf. Sci..

[19]  Z. Werb,et al.  Profiling human breast epithelial cells using single cell RNA sequencing identifies cell diversity , 2018, Nature Communications.

[20]  Johan Hartman,et al.  Chemoresistance Evolution in Triple-Negative Breast Cancer Delineated by Single-Cell Sequencing , 2018, Cell.

[21]  David L. Marron,et al.  Integrated RNA and DNA sequencing reveals early drivers of metastatic breast cancer , 2018, The Journal of clinical investigation.

[22]  Wenjing Chen,et al.  Organotropism: new insights into molecular mechanisms of breast cancer metastasis , 2018, npj Precision Oncology.

[23]  A. Shaw,et al.  Tumour heterogeneity and resistance to cancer therapies , 2018, Nature Reviews Clinical Oncology.

[24]  Mary E. Edgerton,et al.  Multiclonal Invasion in Breast Tumors Identified by Topographic Single Cell Sequencing , 2018, Cell.

[25]  Shao-Cong Sun,et al.  NF-κB signaling in inflammation , 2017, Signal Transduction and Targeted Therapy.

[26]  M. Westphal,et al.  Optical Barcoding for Single-Clone Tracking to Study Tumor Heterogeneity. , 2017, Molecular therapy : the journal of the American Society of Gene Therapy.

[27]  Davide Ferrari,et al.  A Spatio-Temporal Model and Inference Tools for Longitudinal Count Data on Multicolor Cell Growth , 2016, The international journal of biostatistics.

[28]  R. Durbin,et al.  Evaluation of GRCh38 and de novo haploid genome assemblies demonstrates the enduring quality of the reference assembly , 2016, bioRxiv.

[29]  Davide Ferrari,et al.  Semisupervised Clustering by Iterative Partition and Regression with Neuroscience Applications , 2016, Comput. Intell. Neurosci..

[30]  James E. Verdone,et al.  Polyclonal breast cancer metastases arise from collective dissemination of keratin 14-expressing tumor cell clusters , 2016, Proceedings of the National Academy of Sciences.

[31]  J. Mesirov,et al.  The Molecular Signatures Database Hallmark Gene Set Collection , 2015 .

[32]  S. Leung,et al.  PDK1-Dependent Metabolic Reprogramming Dictates Metastatic Potential in Breast Cancer. , 2015, Cell metabolism.

[33]  T. Putoczki,et al.  Emerging roles for IL-11 signaling in cancer development and progression: Focus on breast cancer. , 2015, Cytokine & growth factor reviews.

[34]  K. Polyak,et al.  Tumorigenesis: it takes a village , 2015, Nature Reviews Cancer.

[35]  Steffen Dickopf,et al.  A model of breast cancer heterogeneity reveals vascular mimicry as a driver of metastasis , 2015, Nature.

[36]  R. Ponce,et al.  Tumor necrosis factor, tumor necrosis factor inhibition, and cancer risk , 2015, Current medical research and opinion.

[37]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[38]  Martin Hirst,et al.  DNA barcoding reveals diverse growth kinetics of human breast tumour subclones in serially passaged xenografts , 2014, Nature Communications.

[39]  Sridhar Ramaswamy,et al.  Circulating Tumor Cell Clusters Are Oligoclonal Precursors of Breast Cancer Metastasis , 2014, Cell.

[40]  D. Malide,et al.  In vivo clonal tracking of hematopoietic stem and progenitor cells marked by five fluorescent proteins using confocal and multiphoton microscopy. , 2014, Journal of visualized experiments : JoVE.

[41]  N. Navin,et al.  Clonal Evolution in Breast Cancer Revealed by Single Nucleus Genome Sequencing , 2014, Nature.

[42]  Shawn M. Gillespie,et al.  Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma , 2014, Science.

[43]  C. Klein.  Selection and adaptation during metastatic cancer progression , 2013, Nature.

[44]  A. Chenchik,et al.  Measurement of Cancer Cell Growth Heterogeneity through Lentiviral Barcoding Identifies Clonal Dominance as a Characteristic of In Vivo Tumor Engraftment , 2013, PloS one.

[45]  Thomas Boudier,et al.  TANGO: a generic tool for high-throughput 3D image analysis for studying nuclear organization , 2013, Bioinform..

[46]  Wei Shi,et al.  featureCounts: an efficient general purpose program for assigning sequence reads to genomic features , 2013, Bioinform..

[47]  Ellen T. Gelfand,et al.  The Genotype-Tissue Expression (GTEx) project , 2013, Nature Genetics.

[48]  K. Pogoda,et al.  Analysis of pattern, time and risk factors influencing recurrence in triple-negative breast cancer patients , 2013, Medical Oncology.

[49]  Guangchuang Yu,et al.  clusterProfiler: an R package for comparing biological themes among gene clusters. , 2012, Omics : a journal of integrative biology.

[50]  K. Polyak,et al.  Intra-tumour heterogeneity: a looking glass for cancer? , 2012, Nature Reviews Cancer.

[51]  Michael Thomaschewski,et al.  RGB marking with lentiviral vectors for multicolor clonal cell tracking , 2012, Nature Protocols.

[52]  Caroline Dive,et al.  Clinical significance and molecular characteristics of circulating tumor cells and circulating tumor microemboli in patients with small-cell lung cancer. , 2012, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[53]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[54]  Michael Thomaschewski,et al.  RGB marking facilitates multicolor clonal cell tracking , 2011, Nature Medicine.

[55]  L. Matrisian,et al.  Rapid Extravasation and Establishment of Breast Cancer Micrometastases in the Liver Microenvironment , 2010, Molecular Cancer Research.

[56]  Larry Norton,et al.  Tumor Self-Seeding by Circulating Cancer Cells , 2009, Cell.

[57]  Davis J. McCarthy,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[58]  B. Fehse,et al.  A multicolor panel of novel lentiviral "gene ontology" (LeGO) vectors for functional gene analysis. , 2008, Molecular Therapy.

[59]  K. Hess,et al.  Response to neoadjuvant therapy and long-term survival in patients with triple-negative breast cancer. , 2008, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[60]  Vishva M. Dixit,et al.  IAP Antagonists Induce Autoubiquitination of c-IAPs, NF-κB Activation, and TNFα-Dependent Apoptosis , 2007, Cell.

[61]  J. Minna,et al.  Autocrine TNFalpha signaling renders human cancer cells susceptible to Smac-mimetic-induced apoptosis. , 2007, Cancer cell.

[62]  W. Gerald,et al.  Distinct organ-specific metastatic potential of individual breast cancer cells and primary tumors. , 2005, The Journal of clinical investigation.

[63]  Meng Yang,et al.  Determination of clonality of metastasis by cell-specific color-coded fluorescent-protein imaging. , 2003, Cancer research.

[64]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[65]  B. Beutler,et al.  A tumor necrosis factor (TNF) receptor-IgG heavy chain chimeric protein as a bivalent antagonist of TNF activity , 1991, The Journal of experimental medicine.

[66]  C. Vallejo,et al.  Bone‐Only Versus Visceral‐Only Metastatic Pattern in Breast Cancer: Analysis of 150 Patients A GOCS Study , 1990, American journal of clinical oncology.

[67]  L. Liotta,et al.  The significance of hematogenous tumor cell clumps in the metastatic process. , 1976, Cancer research.

[68]  J. Mesirov,et al.  The Molecular Signatures Database (MSigDB) hallmark gene set collection. , 2015, Cell systems.

[69]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[70]  Gábor Csárdi,et al.  The igraph software package for complex network research , 2006 .

[71]  W. Gerald,et al.  Identifying site-specific metastasis genes and functions. , 2005, Cold Spring Harbor symposia on quantitative biology.

[72]  G. Moore,et al.  The comparative size and structure of tumor cells and clumps in the blood, bone marrow, and tumor imprints , 1960, Cancer.