Extensive intratumoral heterogeneity (ITH) is believed to contribute to therapeutic failure and tumor recurrence, as treatment-resistant cell clones can survive and expand. However, little is known about ITH in triple-negative breast cancer (TNBC) because of the limited number of single-cell sequencing studies on TNBC. In this study, we explored ITH in TNBC by evaluating gene expression-derived and imaging-derived multi-region differences within the same tumor. We obtained tissue specimens from 10 TNBC patients and conducted RNA sequencing analysis of 2–4 regions per tumor. We developed a novel analysis framework to dissect and characterize different types of variability: between-patients (inter-tumoral heterogeneity), between-patients across regions (inter-tumoral and region heterogeneity), and within-patient, between-regions (regional intratumoral heterogeneity). We performed a Bayesian changepoint analysis to assess and classify regional variability as low (convergent) versus high (divergent) within each patient feature (TNBC and PAM50 subtypes, immune, stroma, tumor counts and tumor infiltrating lymphocytes). Gene expression signatures were categorized into three types of variability: between-patients (108 genes), between-patients across regions (183 genes), and within-patients, between-regions (778 genes). Based on the between-patient gene signature, we identified two distinct patient clusters that differed in menopausal status. Significant intratumoral divergence was observed for PAM50 classification, tumor cell counts, and tumor-infiltrating T cell abundance. Other features examined showed a representation of both divergent and convergent results. Lymph node stage was significantly associated with divergent tumors. Our results show extensive intertumoral heterogeneity and regional ITH in gene expression and image-derived features in TNBC. Our findings also raise concerns regarding gene expression based TNBC subtyping. Future studies are warranted to elucidate the role of regional heterogeneity in TNBC as a driver of treatment resistance.

Purpose: Triple-negative breast cancer (TNBC) is a molecularly complex and heterogeneous breast cancer subtype with distinct biological features and clinical behavior. Although TNBC is associated with an increased risk of metastasis and recurrence, the molecular mechanisms underlying TNBC metastasis remain unclear. We performed whole-exome sequencing (WES) analysis of primary TNBC and paired recurrent tumors to investigate the genetic profile of TNBC. Methods: Genomic DNA extracted from 35 formalin-fixed paraffin-embedded tissue samples from 26 TNBC patients was subjected to WES. Of these, 15 were primary tumors that did not have recurrence, and 11 were primary tumors that had recurrence (nine paired primary and recurrent tumors). Tumors were analyzed for single-nucleotide variants and insertions/deletions. Results: The tumor mutational burden (TMB) was 7.6 variants/megabase in primary tumors that recurred (n = 9); 8.2 variants/megabase in corresponding recurrent tumors (n = 9); and 7.3 variants/megabase in primary tumors that did not recur (n = 15). MUC3A was the most frequently mutated gene in all groups. Mutations in MAP3K1 and MUC16 were more common in our dataset. No alterations in PI3KCA were detected in our dataset. Conclusions: We found similar mutational profiles between primary and paired recurrent tumors, suggesting that genomic features may be retained during local recurrence.