Functional and genomic characterisation of a xenograft model system for the study of metastasis in triple-negative breast cancer

Triple-negative breast cancer (TNBC) represents 10-20% of all human ductal adenocarcinomas and has a poor prognosis relative to other subtypes. Hence, new molecular targets for therapeutic intervention are necessary. Analyses of panels of human or mouse cancer lines derived from the same individual that differ in their cellular phenotypes but not in genetic background have been instrumental in defining the molecular players that drive the various hallmarks of cancer. To determine the molecular regulators of metastasis in TNBC, we completed a rigorous in vitro and in vivo characterisation of four populations of the MDA-MB-231 human breast cancer line ranging in aggressiveness from non-metastatic to spontaneously metastatic to lung, liver, spleen and lymph node. Single nucleotide polymorphism (SNP) array analyses and genome-wide mRNA expression profiles of tumour cells isolated from orthotopic mammary xenografts were compared between the four lines to define both cell autonomous pathways and genes associated with metastatic proclivity. Gene set enrichment analysis (GSEA) demonstrated an unexpected association between both ribosome biogenesis and mRNA metabolism and metastatic capacity. Differentially expressed genes or families of related genes were allocated to one of four categories, associated with either metastatic initiation (e.g. CTSC, ENG, BMP2), metastatic virulence (e.g. ADAMTS1, TIE1), metastatic suppression (e.g. CST1, CST2, CST4, CST6, SCNNA1, BMP4) or metastatic avirulence (e.g. CD74). Collectively, this model system based on MDA-MB-231 cells should be useful for the assessment of gene function in the metastatic cascade and also for the testing of novel experimental therapeutics for the treatment of TNBC.