Routine Molecular Profiling in Both Resectable and Unresectable Pancreatic Adenocarcinoma: Relevance of Cytologic Samples.

BACKGROUND & AIMS: Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive disease, for which it is crucial to promptly detect actionable and prognostic alterations to drive specific therapeutic decisions, regardless of tumor resectability status. Endoscopic ultrasonography-guided fine-needle aspiration (EUS-FNA) is of key importance for PDAC diagnosis and can contribute significantly to tumor molecular profiling. METHODS: Comprehensive genomic profile by targeted next-generation sequencing (NGS) was performed on 2 independent PDAC patient cohorts. Cohort 1 consisted of 77 patients with resectable PDAC for whom the histologic sample at the time of resection was available; for 56 patients cytologic specimens at the time of diagnosis also were obtained by EUS-FNA. Cohort 2 consisted of 20 patients with unresectable PDAC, for whom only the EUS-FNA cytologic sample was available. RESULTS: In cohort 1, a complete concordant mutational profile between the cytologic sample at diagnosis and the corresponding histologic specimen after surgery was observed in 88% of the cases, proving the ability to detect potential clinically relevant alterations in cytologic samples by NGS analysis. Notably, clinically actionable mutations were identified in 20% of patients. In cohort 2, comprehensive mutational profiling was obtained successfully for all samples. Consistent with the findings of cohort 1, KRAS, TP53, CDKN2A, and SMAD4 were the most altered genes. Most importantly, 15% of the patients harbored actionable mutations. CONCLUSIONS: Our findings show the feasibility of an NGS approach using both surgical specimens and cytologic samples. The model proposed in this study can be included successfully in the clinical setting for comprehensive molecular profiling of all PDAC patients irrespective of their surgical eligibility.

Characterization of a novel breast cancer cell line derived from a metastatic bone lesion of a breast cancer patient.

PURPOSE: We aimed to generate and characterize a novel cell line from a breast cancer bone metastasis to better study the progression of the disease. METHODS: The cell line, P7731, was derived from a metastatic bone lesion of a breast cancer patient and assessed for marker expression. P7731 was analyzed for DNA copy number variation, somatic mutations, and gene expression and was compared with the primary tumor. RESULTS: P7731 cells are negative for estrogen receptor alpha (ERα), progesterone receptor (PR), and HER2 (triple-negative); strongly express vimentin (100% of cells positive) and also express cytokeratins 8/18 and 19 but at lower frequencies. Flow cytometry indicates P7731 cells are predominantly CD44+/CD49f+/EpCAM-, consistent with a primitive, mesenchymal-like phenotype. The cell line is tumorigenic in immunocompromised mice. Exome sequencing identified a total of 45 and 76 somatic mutations in the primary tumor and cell line, respectively, of which 32 were identified in both samples and included mutations in known driver genes PIK3CA, TP53, and ARID1A. P7731 retains the DNA copy number alterations present in the matching primary tumor. Homozygous deletions detected in the cell line and in the primary tumor were found in regions containing three known (CDKN2A, CDKN2B, and CDKN1B) and 23 putative tumor suppressor genes. Cell line-specific gene amplification coupled with mRNA expression analysis revealed genes and pathways with potential pro-metastatic functions. CONCLUSION: This novel human breast cancer-bone metastasis cell line will be a useful model to study aspects of breast cancer biology, particularly metastasis-related changes from breast to bone.

Multidimensional phenotyping of breast cancer cell lines to guide preclinical research.

PURPOSE: Cell lines are extremely useful tools in breast cancer research. Their key benefits include a high degree of control over experimental variables and reproducibility. However, the advantages must be balanced against the limitations of modelling such a complex disease in vitro. Informed selection of cell line(s) for a given experiment now requires essential knowledge about molecular and phenotypic context in the culture dish. METHODS: We performed multidimensional profiling of 36 widely used breast cancer cell lines that were cultured under standardised conditions. Flow cytometry and digital immunohistochemistry were used to compare the expression of 14 classical breast cancer biomarkers related to intrinsic molecular profiles and differentiation states: EpCAM, CD24, CD49f, CD44, ER, AR, HER2, EGFR, E-cadherin, p53, vimentin, and cytokeratins 5, 8/18 and 19. RESULTS: This cell-by-cell analysis revealed striking heterogeneity within cultures of individual lines that would be otherwise obscured by analysing cell homogenates, particularly amongst the triple-negative lines. High levels of p53 protein, but not RNA, were associated with somatic mutations (p = 0.008). We also identified new subgroups using the nanoString PanCancer Pathways panel (730 transcripts representing 13 canonical cancer pathways). Unsupervised clustering identified five groups: luminal/HER2, immortalised ('normal'), claudin-low and two basal clusters, distinguished mostly by baseline expression of TGF-beta and PI3-kinase pathway genes. CONCLUSION: These features are compared with other published genotype and phenotype information in a user-friendly reference table to help guide selection of the most appropriate models for in vitro and in vivo studies, and as a framework for classifying new patient-derived cancer cell lines and xenografts.

Expression of MAGE-A and NY-ESO-1 cancer/testis antigens is enriched in triple-negative invasive breast cancers.

AIMS: A better understanding of the expression of cancer/testis antigens (CTAs) in breast cancer might enable the identification of new immunotherapy options, especially for triple-negative (TN) tumours, which lack expression of the conventional therapeutic targets oestrogen receptor, progesterone receptor, and human epidermal growth factor receptor 2. The aim of this study was to quantify the expression of MAGE-A and NY-ESO-1 CTAs in breast cancer, and relate this to known clinicopathological parameters. METHODS AND RESULTS: We surveyed MAGE-A and NY-ESO-1 expression in an unselected cohort of 367 breast tumours (of which 65 were TN), with accompanying clinical follow-up data, by using immunohistochemical analysis of tissue microarrays. Relevant to their potential as vaccine targets in breast cancer, MAGE-A was expressed in 13% of cases, and NY-ESO-1 in 3.8%, with the majority of tumours showing fairly homogeneous staining within individual tissue cores (~85% of cases with staining in >75% of tumour cells). Most NY-ESO-1-positive cases also expressed MAGE-A (P = 2.06 × 10-9 ), and both were strongly associated with the TN phenotype (P < 0.0001), with the most proliferative and poorly differentiated cases, in paticular, showing genomic instability. This was characterised by coexpression of c-Kit and TTK, and overexpression of p53. CONCLUSIONS: MAGE-A and NY-ESO-1 are frequently expressed in TN breast cancer (~47% and 17% of TN cases, respectively), suggesting that targeting them could be feasible in this patient group. Expression is reasonably homogeneous in positive cases, suggesting that immunohistochemical analysis of tissue biopsies would be a reliable companion biomarker.

High content screening application for cell-type specific behaviour in heterogeneous primary breast epithelial subpopulations.

BACKGROUND: The complex interaction between multiple cell types and the microenvironment underlies the diverse pathways to carcinogenesis and necessitates sophisticated approaches to in vitro hypotheses testing. The combination of mixed culture format with high content immunofluorescence screening technology provides a powerful platform for observation of cell type specific behavior. METHODS: We have developed a versatile, high-throughput method for assessing cell-type specific responses. In addition to the specificity and sensitivity offered traditionally by immunofluorescent detection in flow cytometry, the 'in-cell' analysis method we describe provides the added benefits of higher throughput and the ability to analyse protein subcellular localisation in situ. Furthermore, elimination of the cell dissociation step allows for more immediate analysis of responses to specific extrinsic stimuli. We applied this method to investigate ionising radiation treatment response in normal breast epithelial cells, measuring growth rate, cell cycle response and double-strand DNA breaks. RESULTS: The 'in-cell' analysis approach elucidated several interesting donor and cell-type specific differences. Notably, in response to ionizing radiation we observed differential expression in luminal and basal-like cells of a member of the APOBEC enzyme family, recently identified as a critical driver of an oncogenic signature. Our findings suggest that this enzyme is active in the normal breast epithelium during DNA damage response. CONCLUSIONS: We demonstrate the practical application of a new method for assessing cell-type specific change in mixed cultures, especially the analysis of normal primary cultures, overcoming a major technical issue of dissecting the response of multiple cell types in a heterogeneous population.

Assessment of ORAI1-mediated basal calcium influx in mammary epithelial cells.

BACKGROUND: The entry of calcium ions into mammary gland epithelial cells is one of the least well-understood processes in the transport of calcium into milk during lactation. The store-operated calcium entry channel ORAI1, has been suggested as a potential mechanism for the entry of Ca(2+) into mammary gland epithelial cells from the maternal blood supply during lactation. The down regulation of the canonical ORAI1 activator STIM1 during lactation suggests that other known ORAI activators such as STIM2 and SPCA2 may be important during lactation. RESULTS: Differentiation of HC11 mammary gland epithelial cells was associated with enhanced basal Ca(2+) influx. Silencing of Orai1 abolished this enhancement of Ca(2+) influx. Stim2 had a modest effect on Ca(2+) influx in this in vitro model of lactation, whereas Stim1 and Spca2 silencing had no effect. Despite pronounced increases in Spca2 mRNA during lactation there was no change in the generation of the alternative splice product generated by Mist1, which increases during lactation. CONCLUSIONS: These studies support the hypothesis that lactation is associated with a remodelling of Ca(2+) influx and this is associated with enhancement of basal Ca(2+) influx. This enhanced Ca(2+) influx appears to occur through the calcium channel Orai1.

SNORD-host RNA Zfas1 is a regulator of mammary development and a potential marker for breast cancer.

Long noncoding RNAs (lncRNAs) are increasingly recognized to play major regulatory roles in development and disease. To identify novel regulators in breast biology, we identified differentially regulated lncRNAs during mouse mammary development. Among the highest and most differentially expressed was a transcript (Zfas1) antisense to the 5' end of the protein-coding gene Znfx1. In vivo, Zfas1 RNA is localized within the ducts and alveoli of the mammary gland. Zfas1 intronically hosts three previously undescribed C/D box snoRNAs (SNORDs): Snord12, Snord12b, and Snord12c. In contrast to the general assumption that noncoding SNORD-host transcripts function only as vehicles to generate snoRNAs, knockdown of Zfas1 in a mammary epithelial cell line resulted in increased cellular proliferation and differentiation, while not substantially altering the levels of the SNORDs. In support of an independent function, we also found that Zfas1 is extremely stable, with a half-life >16 h. Expression analysis of the SNORDs revealed these were expressed at different levels, likely a result of distinct structures conferring differential stability. While there is relatively low primary sequence conservation between Zfas1 and its syntenic human ortholog ZFAS1, their predicted secondary structures have similar features. Like Zfas1, ZFAS1 is highly expressed in the mammary gland and is down-regulated in breast tumors compared to normal tissue. We propose a functional role for Zfas1/ ZFAS1 in the regulation of alveolar development and epithelial cell differentiation in the mammary gland, which, together with its dysregulation in human breast cancer, suggests ZFAS1 as a putative tumor suppressor gene.

Gene expression profiling of tumour epithelial and stromal compartments during breast cancer progression.

The progression of ductal carcinoma in situ (DCIS) to invasive ductal carcinoma (IDC) marks a critical step in the evolution of breast cancer. There is some evidence to suggest that dynamic interactions between the neoplastic cells and the tumour microenvironment play an important role. Using the whole-genome cDNA-mediated annealing, selection, extension and ligation assay (WG-DASL, Illumina), we performed gene expression profiling on 87 formalin-fixed paraffin-embedded (FFPE) samples from 17 patients consisting of matched IDC, DCIS and three types of stroma: IDC-S (<3 mm from IDC), DCIS-S (<3 mm from DCIS) and breast cancer associated-normal stroma (BC-NS; >10 mm from IDC or DCIS). Differential gene expression analysis was validated by quantitative real time-PCR, immunohistochemistry and immunofluorescence. The expression of several genes was down-regulated in stroma from cancer patients relative to normal stroma from reduction mammoplasties. In contrast, neoplastic epithelium underwent more gene expression changes during progression, including down regulation of SFRP1. In particular, we observed that molecules related to extracellular matrix (ECM) remodelling (e.g. COL11A1, COL5A2 and MMP13) were differentially expressed between DCIS and IDC. COL11A1 was overexpressed in IDC relative to DCIS and was expressed by both the epithelial and stromal compartments but was enriched in invading neoplastic epithelial cells. The contributions of both the epithelial and stromal compartments to the clinically important scenario of progression from DCIS to IDC. Gene expression profiles, we identified differential expression of genes related to ECM remodelling, and specifically the elevated expression of genes such as COL11A1, COL5A2 and MMP13 in epithelial cells of IDC. We propose that these expression changes could be involved in facilitating the transition from in situ disease to invasive cancer and may thus mark a critical point in disease development.

Breast cancer stem cells: treatment resistance and therapeutic opportunities.

The clinical and pathologic heterogeneity of human breast cancer has long been recognized. Now, molecular profiling has enriched our understanding of breast cancer heterogeneity and yielded new prognostic and predictive information. Despite recent therapeutic advances, including the HER2-specific agent, trastuzumab, locoregional and systemic disease recurrence remain an ever-present threat to the health and well being of breast cancer survivors. By definition, disease recurrence originates from residual treatment-resistant cells, which regenerate at least the initial breast cancer phenotype. The discovery of the normal breast stem cell has re-ignited interest in the identity and properties of breast cancer stem-like cells and the relationship of these cells to the repopulating ability of treatment-resistant cells. The cancer stem cell model of breast cancer development contrasts with the clonal evolution model, whereas the mixed model draws on features of both. Although the origin and identity of breast cancer stem-like cells is contentious, treatment-resistant cells survive and propagate only because aberrant and potentially druggable signaling pathways are recruited. As a means to increase the rates of breast cancer cure, several approaches to specific targeting of the treatment-resistant cell population exist and include methods for addressing the problem of radioresistance in particular.

HER3 and downstream pathways are involved in colonization of brain metastases from breast cancer.

INTRODUCTION: Metastases to the brain from breast cancer have a high mortality, and basal-like breast cancers have a propensity for brain metastases. However, the mechanisms that allow cells to colonize the brain are unclear. METHODS: We used morphology, immunohistochemistry, gene expression and somatic mutation profiling to analyze 39 matched pairs of primary breast cancers and brain metastases, 22 unmatched brain metastases of breast cancer, 11 non-breast brain metastases and 6 autopsy cases of patients with breast cancer metastases to multiple sites, including the brain. RESULTS: Most brain metastases were triple negative and basal-like. The brain metastases over-expressed one or more members of the HER family and in particular HER3 was significantly over-expressed relative to matched primary tumors. Brain metastases from breast and other primary sites, and metastases to multiple organs in the autopsied cases, also contained somatic mutations in EGFR, HRAS, KRAS, NRAS or PIK3CA. This paralleled the frequent activation of AKT and MAPK pathways. In particular, activation of the MAPK pathway was increased in the brain metastases compared to the primary tumors. CONCLUSIONS: Deregulated HER family receptors, particularly HER3, and their downstream pathways are implicated in colonization of brain metastasis. The need for HER family receptors to dimerize for activation suggests that tumors may be susceptible to combinations of anti-HER family inhibitors, and may even be effective in the absence of HER2 amplification (that is, in triple negative/basal cancers). However, the presence of activating mutations in PIK3CA, HRAS, KRAS and NRAS suggests the necessity for also specifically targeting downstream molecules.

Identification and functional analysis of novel BRCA1 transcripts, including mouse Brca1-Iris and human pseudo-BRCA1.

Recent characterization of the mammalian transcriptome has confirmed its predicted complexity, with many loci encoding multiple splice variants and pseudogenes. The breast cancer susceptibility gene BRCA1 is a tumour suppressor gene that produces multiple functional transcripts. For example, BRCA1-IRIS is a splice variant of BRCA1, which encodes a protein that is functionally distinct from BRCA1. Here we describe the identification of ten novel Brca1 splice variants including Brca1-Iris, the mouse orthologue of human BRCA1-IRIS. We show that Brca1-Iris is differentially expressed during mammary epithelial differentiation and regulates survival of mammary epithelial cells. Another transcript, Brca1-Delta22, expressed in both mouse and human cells, was found to be defective in transcriptional activation capacity. Finally, we show that the human BRCA1 pseudogene produces a spliced pseudoBRCA1 transcript. The identification of these transcripts has implications for the understanding of the role of BRCA1 in biology and disease and for the interpretation of mouse knockout models.

Non-coding RNAs underlie genetic predisposition to breast cancer.

BACKGROUND: Genetic variants identified through genome-wide association studies (GWAS) are predominantly non-coding and typically attributed to altered regulatory elements such as enhancers and promoters. However, the contribution of non-coding RNAs to complex traits is not clear. RESULTS: Using targeted RNA sequencing, we systematically annotated multi-exonic non-coding RNA (mencRNA) genes transcribed from 1.5-Mb intervals surrounding 139 breast cancer GWAS signals and assessed their contribution to breast cancer risk. We identify more than 4000 mencRNA genes and show their expression distinguishes normal breast tissue from tumors and different breast cancer subtypes. Importantly, breast cancer risk variants, identified through genetic fine-mapping, are significantly enriched in mencRNA exons, but not the promoters or introns. eQTL analyses identify mencRNAs whose expression is associated with risk variants. Furthermore, chromatin interaction data identify hundreds of mencRNA promoters that loop to regions that contain breast cancer risk variants. CONCLUSIONS: We have compiled the largest catalog of breast cancer-associated mencRNAs to date and provide evidence that modulation of mencRNAs by GWAS variants may provide an alternative mechanism underlying complex traits.

Current and future approach to the pathologist's assessment for targeted therapy in breast cancer.

Breast cancer is a common disease in the population. Contrary to public perception, it is a heterogeneous disease with varying morphology, prognosis and response to therapy. The pathological analysis is at the heart of information provided to surgeons and oncologists to plan further management. The pathologist is increasingly asked to test for biomarkers that provide prognostic and predictive information to direct treatment. Staining cancers for ER, PgR and HER2 has become routine and it is likely that addition of other biomarkers including 'basal markers', VEGF and growth factor receptors such as HER1 (EGFR) will soon follow. Microarray based genomic, transcription and proteomic methods are changing our classification systems and identifying novel targets for the development of new therapeutics. It is important for pathologists to appreciate and embrace the new developments as they will impact on daily clinical practice and require accurate assessment of biomarkers to determine treatment options as part of multidisciplinary teams.

Localization of plasma membrane and secretory calcium pumps in the mammary gland.

Until recently the mechanism for the enrichment of milk with calcium was thought to be almost entirely via the secretory pathway. However, recent studies suggest that a plasma membrane calcium ATPase, PMCA2, is the primary mechanism for calcium transport into milk, highlighting a major role for apical calcium transport. We compared the expression of the recently identified secretory calcium ATPase, SPCA2, and SPCA1, in the mouse mammary gland during development. SPCA2 levels increased over 35-fold during lactation with expression localized to luminal secretory cells, while SPCA1 increased only a modest 2-fold and was expressed throughout the cells of the mammary gland. We also observed major differences in the localization of PMCA2 and PMCA1. Our studies highlight the likely specific roles of PMCA2 and SPCA2 in lactation and indicate that calcium transport into milk is a complex interplay between apical and secretory pathways.

Targeted disruption of Brca1 in restricted compartments of the mouse mammary epithelia.

Tumours arising in BRCA1 mutation carriers have a characteristic phenotype, the molecular and cellular basis of which is unknown. To address the hypothesis that this phenotype reflects a role for BRCA1 in either in the basal or the stem cell compartments of the mammary epithelia, we have targeted its disruption to K14 and K6a expressing cells of the mouse. Unlike MMTV and WAP driven conditional knockout models of Brca1, these two models did not result in any observable changes in the mammary gland. Our results suggest that BRCA1-associated tumours arise either in K14 and K6a negative basal cells of the mammary gland, or possibly from transdifferentiation of luminal epithelia.

Regulation of BRCA1 messenger RNA stability in human epithelial cell lines and during cell cycle progression.

The mechanisms regulating expression of breast cancer 1 (BRCA1) are not fully characterised. By studying regulation of endogenous BRCA1 in human epithelial cell lines and during cell cycle progression, we provide evidence to suggest BRCA1 is regulated post-transcriptionally at the level of messenger RNA stability. We also show that RNA-binding proteins associate with an AU-rich, cis-active sequence of the BRCA1 3' untranslated region in a cell cycle-dependent manner. Our data identify a new post-transcriptional regulatory axis and a novel mechanism for modulating the levels of BRCA1 protein, with possible implications for understanding the mechanisms underlying BRCA1 repression in breast cancer.

The calcium pump plasma membrane Ca(2+)-ATPase 2 (PMCA2) regulates breast cancer cell proliferation and sensitivity to doxorubicin.

Regulation of Ca(2+) transport is vital in physiological processes, including lactation, proliferation and apoptosis. The plasmalemmal Ca(2+) pump isoform 2 (PMCA2) a calcium ion efflux pump, was the first protein identified to be crucial in the transport of Ca(2+) ions into milk during lactation in mice. In these studies we show that PMCA2 is also expressed in human epithelia undergoing lactational remodeling and also report strong PMCA2 staining on apical membranes of luminal epithelia in approximately 9% of human breast cancers we assessed. Membrane protein expression was not significantly associated with grade or hormone receptor status. However, PMCA2 mRNA levels were enriched in Basal breast cancers where it was positively correlated with survival. Silencing of PMCA2 reduced MDA-MB-231 breast cancer cell proliferation, whereas silencing of the related isoforms PMCA1 and PMCA4 had no effect. PMCA2 silencing also sensitized MDA-MB-231 cells to the cytotoxic agent doxorubicin. Targeting PMCA2 alone or in combination with cytotoxic therapy may be worthy of investigation as a therapeutic strategy in breast cancer. PMCA2 mRNA levels are also a potential tool in identifying poor responders to therapy in women with Basal breast cancer.

Emerging role of long non-coding RNA SOX2OT in SOX2 regulation in breast cancer.

The transcription factor SOX2 is essential for maintaining pluripotency in a variety of stem cells. It has important functions during embryonic development, is involved in cancer stem cell maintenance, and is often deregulated in cancer. The mechanism of SOX2 regulation has yet to be clarified, but the SOX2 gene lies in an intron of a long multi-exon non-coding RNA called SOX2 overlapping transcript (SOX2OT). Here, we show that the expression of SOX2 and SOX2OT is concordant in breast cancer, differentially expressed in estrogen receptor positive and negative breast cancer samples and that both are up-regulated in suspension culture conditions that favor growth of stem cell phenotypes. Importantly, ectopic expression of SOX2OT led to an almost 20-fold increase in SOX2 expression, together with a reduced proliferation and increased breast cancer cell anchorage-independent growth. We propose that SOX2OT plays a key role in the induction and/or maintenance of SOX2 expression in breast cancer.

Kinome profiling reveals breast cancer heterogeneity and identifies targeted therapeutic opportunities for triple negative breast cancer.

Our understanding of breast cancer heterogeneity at the protein level is limited despite proteins being the ultimate effectors of cellular functions. We investigated the heterogeneity of breast cancer (41 primary tumors and 15 breast cancer cell lines) at the protein and phosphoprotein levels to identify activated oncogenic pathways and developing targeted therapeutic strategies. Heterogeneity was observed not only across histological subtypes, but also within subtypes. Tumors of the Triple negative breast cancer (TNBC) subtype distributed across four different clusters where one cluster (cluster ii) showed high deregulation of many proteins and phosphoproteins. The majority of TNBC cell lines, particularly mesenchymal lines, resembled the cluster ii TNBC tumors. Indeed, TNBC cell lines were more sensitive than non-TNBC cell lines when treated with targeted inhibitors selected based on upregulated pathways in cluster ii. In line with the enrichment of the upregulated pathways with onco-clients of Hsp90, we found synergy in combining Hsp90 inhibitors with several kinase inhibitors, particularly Erk5 inhibitors. The combination of Erk5 and Hsp90 inhibitors was effective in vitro and in vivo against TNBC leading to upregulation of pro-apoptotic effectors. Our studies contribute to proteomic profiling and improve our understanding of TNBC heterogeneity to provide therapeutic opportunities for this disease.

Thrombospondin-4 expression is activated during the stromal response to invasive breast cancer.

The thromobospondins are a family of extracellular glycoproteins that are activated during tissue remodeling processes such as embryogenesis, wound healing and cancer. Thrombospondin-4 (THBS4) is known to have roles in cellular migration, adhesion and attachment, as well as proliferation in different contexts. Data to support a role in cancer biology is increasing, including for gastrointestinal and prostate tumours. Here, using a combination of immunohistochemistry, immunofluorescence and analysis of publicly available genomic and expression data, we present the first study describing the pattern of expression of THBS4 in normal breast and breast cancer. THBS4 was located to the basement membrane of large ducts and vessels in normal breast tissue, but was absent from epithelium and extracellular matrix. There was a significant induction in expression in cancer-associated stroma relative to normal stroma (P = 0.0033), neoplastic epithelium (P < 0.0001) and normal epithelium (P < 0.0001). There was no difference in stromal expression of THBS4 between invasive ductal carcinomas (IDC) and invasive lobular carcinomas (ILC). The THBS4 mRNA levels were variable yet were generally highest in tumours typically rich in stromal content (ILC, ER positive low grade IDC; luminal A and normal-like subtypes). Genomic alterations of the THBS4 gene (somatic mutations and gene copy number) are rare suggesting this dramatic activation in expression is most likely dynamically regulated through the interaction between invading tumour cells and stromal fibroblasts in the local microenvironment. In summary, THBS4 expression in breast cancer-associated extracellular matrix contributes to the activated stromal response exhibited during tumour progression and this may facilitate invasion of tumour cells.