Distinct Gene Expression Profiles of Matched Primary and Metastatic Triple-Negative Breast Cancers

Background: Although triple-negative breast cancer (TNBC) is associated with an increased risk of recurrence and metastasis, the molecular mechanisms underlying metastasis in TNBC remain unknown. To identify transcriptional changes and genes regulating metastatic progression in TNBC, we compared the transcriptomic profiles of primary and matched metastatic tumors using massively parallel RNA sequencing. Methods: We performed gene expression profiling using formalin-fixed paraffin-embedded (FFPE) TNBC tissues of patients from two cohorts: the Zurich cohort (n = 31) and the Stavanger cohort (n = 5). Among the 31 patients in the Zurich cohort, 18 had primary TNBC tumors that did not metastasize, and 13 had primary tumors that metastasized (11 paired primary and locoregional recurrences). The Stavanger cohort included five matched primary and metastatic TNBC tumors. Significantly differentially expressed genes (DEGs; absolute fold change ≥2, p < 0.05) were identified and subjected to functional analyses. We investigated if there was any overlap between DEGs from both the cohorts with epithelial-to-mesenchymal-to-amoeboid transition (EMAT) gene signature. xCell was used to estimate relative fractions of 64 immune and stromal cell types in each RNA-seq sample. Results: In the Zurich cohort, we identified 1624 DEGs between primary TNBC tumors and matched metastatic lesions. xCell analysis revealed a significantly higher immune scores for metastatic lesions compared to paired primary tumors in the Zurich cohort. We also found significant upregulation of three MammaPrint signature genes (HRASLS, TGFB3 and RASSF7) in primary tumors that metastasized compared to primary tumors that remained metastasis-free. In the Stavanger cohort, we identified 818 DEGs between primary tumors and matched metastatic lesions. No significant differences in xCell immune scores were observed. We found that 21 and 14 DEGs from Zurich and Stavanger cohort, respectively, overlapped with the EMAT gene signature. In both cohorts, genes belonging to the MMP, FGF, and PDGFR families were upregulated in primary tumors compared to matched metastatic lesions. Conclusions: Our results suggest that distinct gene expression patterns exist between primary TNBCs and matched metastatic tumors. Further studies are warranted to explore whether these discrete expression profiles underlie or result from disease status.