Abstract
Triple-negative breast cancer (TNBC) has been subdivided into six distinct subgroups: basal-like 1 (BL1), basal-like 2 (BL2), mesenchymal (M), mesenchymal stem–like (MSL), immunomodulatory (IM), and luminal androgen receptor (LAR). We recently identified a subgroup of TNBC with loss of the tumor suppressor PTEN and five specific microRNAs that exhibits exceedingly poor clinical outcome and contains TP53 mutation, RB1 loss and high MYC and WNT signalling. Here, show that these PTEN-low/miRNA-low lesions cluster with BL1 TNBC. These tumors exhibited high RhoA signalling and were significantly stratified on the basis of PTEN-low/RhoA-signalling-high with hazard ratios (HRs) of 8.2 (P = 0.0009) and 4.87 (P = 0.033) in training and test cohorts, respectively. For BL2 TNBC, we identified AKT1 copy gain/high mRNA expression as surrogate for poor prognosis (HR = 3.9; P = 0.02 and HR = 6.1; P = 0.0032). In IM, programmed cell death 1 (PD1) was elevated and predictive of poor prognosis (HR = 5.3; P = 0.01 and HR = 3.5; P < 0.004). Additional alterations, albeit without prognostic power, characterized each subtype including high E2F2 and TGFβ signalling and CXCL8 expression in BL2, high IFNα and IFNγ signalling and CTLA4 expression in IM, and high EGFR signalling in MSL, and may be targeted for therapy. This study identified PTEN-low/RhoA-signalling-high, and high AKT1 and PD1 expression as potent prognostications for BL1, BL2 and IM subtypes with survival differences of over 14, 2.75 and 10.5 years, respectively. This intrinsic heterogeneity could be exploited to prioritize patients for precision medicine.
Similar content being viewed by others
Introduction
Breast Cancer (BC) is pathologically classified as oestrogen-positive (ER+), HER2/ERBB2/NEU-positive (HER2+) and triple negative (TNBC) subtypes1,2,3. The latter group represents ~15% of all BC cases but has poor prognosis and affects young women with a dramatically higher incidence in African and African-American women4,5. TNBC can be sub-divided into 6 subtypes: basal-like (BL1 and BL2), mesenchymal (M), mesenchymal stem-like (MSL), immunomodulatory (IM), and luminal androgen receptor (LAR), as well as an unspecified group (UNS)3. An alternative classification divides TNBC into BL1 and BL2, M and LAR6. Whether these subgroups can be further stratified, for example through multi-omic approaches7, is largely unknown.
Approximately 20% of TNBC patients respond well to standard therapy (tumor resection, radiation and cytotoxic chemotherapy), but the rest develop lethal metastatic disease. A recent clinical trial has demonstrated an encouraging response of some TNBC patients to the immune-checkpoint blockade (ICB) inhibitor, Atezolizumab (anti-PD-L1), in combination with conventional chemotherapy8 (Reviewed in9). Yet, most patients, including those with high expression of PD-L1, succumbed to the disease. The identification of TNBC subtype that are likely to respond to ICB therapy is of great interest.
Germ line mutations in the tumor suppressor genes BRCA1 and BRCA2 lead to basal-like BC10,11,12; and many TNBC with intact BRCA1/2 are nonetheless classified as BRCAness lesions13. PARP inhibitors are synthetic lethal with BRCA1/2 mutant/BRCAness TNBC and have been approved for therapy14. However, emergence of clones that resist PARP inhibition through multiple mechanisms is a major clinical problem15. Genomic analysis of sporadic TNBC identified combined loss of RB1 (mutation/deletion; promoter/mRNA silencing) plus TP53 (mutation/deletion) in as many as 28–40% of cases16,17,18,19. Additional alterations include PTEN loss or PIK3CA mutation, and enhanced EGFR, WNT and MYC signaling. RB1, PTEN and TP53 are also the most frequent drivers of metastasis in diverse types of solid human cancers including breast cancer20.
We recently demonstrated that inactivation of Pten in the mouse mammary gland induces mammary tumors that fail to induce secondary tumors after orthotopic injection into recipient mice21. An exception was a relatively rare group of tumors resembling basal-like BC, which could be efficiently propagated in recipient mice. These transplantable mouse tumors exhibited low expression of the tumor suppressor microRNA-145, raising the question of whether in breast cancer patients, PTEN-deficiency cooperates with miR-145 loss and/or other microRNAs to define an aggressive subgroup of TNBCs. We found that low expression of four of the following microRNAs: hsa-miR-145, hsa-miR-4324, hsa-miR-125b, hsa-miR-381 and hsa-miR136 in cooperation with PTEN loss marks highly aggressive TNBCs. These PTEN-low/miRs-low TNBCs exhibit TP53 mutation (not deletion), loss of RB1 signature, and high MYC, WNT and PI3K signalling.
Here we sought to determine whether these PTEN-low/miRs-low TNBCs are spread among different TNBC subtypes or concentrate within one particular subgroup. We found that these tumors cluster almost exclusively with the BL1 subtype and represent lesions with the poorest clinical outcome. Further analysis revealed that this aggressive subgroup of BL1 TNBC exhibit unique signalling attributes such as high RhoA pathway activity. Furthermore, PTEN-low/Rho-signalling high BL1 TNBC patients had the worst prognosis. This observation indicates that the broad survival curve of BL1 patients (and by extension, other subtypes) is not merely stochastic in nature but a reflection of intrinsic heterogeneity that can be molecularly identified and potentially exploited therapeutically. These findings prompted further analysis of other TNBC subtypes and the identification of AKT1 copy number gain and high expression as predictor of poor prognosis in BL2, and high Programmed Cell Death 1 (PD1) expression as predictor of poor prognosis in IM TNBC. This and additional analysis presented herein provide evidence that patients within each TNBC subtype can be molecularly stratified and prioritized for targeted therapy.
Results
Pten-low/miR-low TNBCs cluster with BL1 lesions with exceedingly poor clinical outcome
We first classified a publically available dataset of 205 TNBCs with mRNA and miRNA expression and clinical data into 6 (BL1, BL2, M, MSL, IM, LAR) or 4 (BL1 and BL2, M and LAR) subtypes as described by Lehmann et al.3,6, (Fig. 1A). TNBCs that did not cluster with any of these subtypes were designated UNS (unspecified). In addition, “PAM50 + cla” classification identified most TNBCs as Basal-like or Claudin-low (cla) lesions, as well as some Normal-like and HER2+ tumors. We then asked how PTEN-low and PTEN-low/miRNA-low (group “a”) TNBCs clustered with respect to these TNBC subtypes. Nine of the eleven PTEN-low/miRNA-low TNBCs grouped together with BL1. The remaining two tumors clustered with BL2 and UNS, though the one that classified with UNS (r = 0.5454; p < 0.0001) by correlation to TNBCType centroids was very close to BL1 (r = 0.5188; p < 0.0001).
Genomic and transcriptomic analysis has stratified breast cancer into 10 intrinsic subtypes with IntClust10 comprising most basal-like breast cancers16,22. Consistent with this, all PTEN-low/miRNA-low TNBCs that were classified as BL1 also clustered with IntClust10 (Fig. 1A). Notably, the UNS group showed similar PAM50 and IntClust10 clustering as BL1. In addition, while BL-1 group “a” had poor outcome, other pathological parameters were similar across the different TNBC subtypes (black boxes in Fig. 1A).
Kaplan-Meier survival curves for the six Lehmann TNBC-subtypes as well as UNS tumors are shown in Fig. 1B. BL1 TNBCs in this dataset are second most aggressive after BL2. When BL1 tumors were divided according to PTEN and miRNA status, the PTEN-low/miRNA-low BL1 subgroup (group “a”) exhibited the worst clinical outcome compared to all other BL1 or other TNBC subtypes including BL2, with hazard ratios (HRs) of 4.2 and 4.37, respectively (P < 0.0001; Fig. 1C,D). Using a validation cohort of 44 TNBC, PTEN-low/miRNA-low TNBCs also clustered with BL1 and stratified this subgroup into poor and relatively good clinical outcome; HR for the PTEN-low/miRNA-low subgroup compared with the other BL1 lesions, or relative to all other TNBCs was 3.77 (P = 0.071) and 3.07 (P = 0.014), respectively (Fig. 1E). Specifically, the median latency was 32 and 31.4 months for PTEN-low/miRNA-low BL1 patients versus 204 months or undefined period (i.e. % survival not crossing the 50% mark) for all other BL1 patients in the two cohorts.
We previously showed that PTEN-low/miRNA-low TNBCs exhibited high WNT signalling and loss of RB1 signature23. We therefore asked whether these two classifiers were also able to stratify BL1 patients. Using two different signatures for RB1 loss19,23,24,25, PTEN-low/RB-loss and PTEN-low/β-catenin-high TNBCs almost completely overlapped with the PTEN-low/miRNA-low subgroup (Fig. 1A). Moreover, these classifiers significantly stratified BL1 patients into poor and good clinical outcomes on both test and validation cohorts with high hazard ratios of 3.3 and 4.7 for PTEN-low/β-catenin pathway high, and 3.8 to 9.1 for PTEN-low/RB-signature-low (Fig. 1F). Thus, in two independent cohorts, PTEN-low/miRNA-low TNBCs cluster with BL1 and represent patients with extremely poor prognosis. Moreover, these patients can be identified using PTEN-low and either β-catenin high or RB1-loss signatures.
Unique signalling pathways in different TNBC subtypes; high RhoA signalling predicts poor clinical outcome in BL1 patients
The aforementioned results showing that BL1 TNBC patients can be molecularly stratified into good-poor prognosis groups raised the question of whether other TNBC subtypes might also be clinically stratified. To identify Lehmann TNBC subtype-specific alterations, we first determined activity for twenty-five pathways/signatures using established conditions25,26 (Fig. 2A). Consistent with our previous study, BL1 and particularly PTEN-low/miRNA-low BL1 lesions exhibited RB1-loss signature and high WNT and MYC signalling. Within BL1, the MYC signature alone sufficed to stratify this subtype into poor and moderate clinical outcome with HRs of ~2.1 in the test (205 patient) and validation (190 patient) cohorts (P = 0.027; P = 0.036; Fig. 2B,C). The validation cohort (190 patients) included 146 TNBCs from METABRIC (excluding the 205 TNBC samples in the test cohort) and the 44 TNBCs in cohort GSE22220 (both Illumina platforms).
Several signalling pathways were elevated in specific TNBC subtypes (Fig. 2D–H). Thus, E2F2 and TGFβ pathways were specifically elevated in BL2; IFNα and IFNγ pathways were high in IM; and EGFR signaling was prominent in MLS relative to other subtypes. With the exception of IFNα in one clinical cohort, none of these pathways had any prognostication value in the respective TNBC subtype (Supplemental Fig. 1). Nonetheless, these pathways represent specific markers and potential targets for therapy (see Discussion).
Importantly, RhoA signalling was significantly elevated in PTEN-low/miRNA-low (group “a”) TNBC, compared with PTEN-low/miRNA-high (group “b”), PTEN-low/miRNA-others (group “c”) or PTEN-high (Fig. 3A), or in BL1 (Fig. 3B) vs most other subtypes in two different cohorts. Moreover, high RhoA signalling identified TNBCs with exceedingly poor clinical outcome (Fig. 3A,B). High RhoA signalling was as efficient as or an even better predictor than RB1-loss in identifying TNBC with poor prognosis. Moreover, in combination with PTEN status, PTEN-low/RhoA-high BL1 TNBC showed the poorest clinical outcome compared with all other BL1 lesions with HRs of 8.2 (P = 0.0009) and 4.87 (P = 0.025) in the two clinical cohorts (Fig. 3C). Median latency was 32 and 42.7 months for PTEN-low/RhoA-high BL1 relative to 204 or undefined period for all other BL1 patients in the two cohorts – a difference of over 14 years. Thus, RhoA signalling represents another targetable pathway that is elevated in BL1 and, together with PTEN-low status, predictive of poor outcome.
Unique spectrum of mutations and gains/losses in different TNBC subtypes
To ask whether other TNBC subtypes can be molecularly stratified into tumors with distinct clinical outcomes, we analyzed the spectrum of mutations and copy number alterations (CNAs) in each subtype. We identified 135 genes with at least one mutation in 185 TNBCs (Supplemental Fig. 2). Only TP53 and PIK3CA were frequently mutated (Supplemental Fig. 3A,B). TNBC samples showing no TP53 mutation likely harbor p53 deletion or amplification of HDM2, the E3 ligase for this tumor suppressor27. Indeed, pathway activity for p53 showed loss of this tumor suppressor signaling in nearly all TNBCs, including all BL1, BL2 and IM (Fig. 2A). In accordance, TP53 mutation alone did not predict poor prognosis (not shown). In contrast, TNBCs with PIK3CA mutation showed significantly worst prognosis relative to TNBCs with wild-type PIK3CA in two different cohorts (Supplemental Fig. 3C). Furthermore, TNBCs selected on the basis of combined PIK3CA plus TP53 mutations showed similar poor prognosis as those selected on the basis of PIK3CA mutation alone (Supplemental Fig. 3D).
CNA analysis identified alterations using 22544 probes in all TNBC subtypes with the exception of MSLs (arrow in Supplemental Fig. 4). In search for specific regions that are gained or lost in each TNBC subtype using ANOVA-test-selected 2744 CNA probes (Supplemental Fig. 5), we identified common gains in chromosome 16p13.3 in BL2. This region contains 169 genes, of which, 12 are shown in Supplemental Fig. 6. These include 3-Phosphoinositide Dependent Protein Kinase 1 (PDPK1 also known as PDK1), which phosphorylates and activates a subgroup of AGC protein kinases like protein kinase B (PKBα/AKT1, PKBβ/AKT2, PKBγ/AKT3) downstream of PI3K. Chr16p13.3 gain is implicated in prostate cancer progression via induction of PDPK128. However, although PDPK1 mRNA levels are significantly higher in BL2 compared to other subtypes, its expression did not correlate with poor prognosis in all TNBCs or specifically within BL2 (Supplemental Fig. 6). Likewise, expression of other genes and microRNAs such as miR3176, located within the 16p13.3 amplicon, was also elevated in BL2, but did not predict clinical outcome.
High copy number alteration (CNA) and mRNA expression of AKT1 predict poor clinical outcome in BL2 TNBC subtype
Compilation of CNA, exome and whole genome sequencing analysis identified 94 common oncogenic alterations in breast cancer18. Having found no obvious large CNA associated with poor prognosis, we asked whether any of the 94 common oncogenic alterations in breast cancer are associated with CNA in specific TNBC subtypes. Twelve such genes were identified including MYC, GATA3, AXIN1, FBXW7, AKT1, MAP3K1, APC, PIK3R1, BULB1B, PBRM1 and PTEN (Fig. 4A) from the ANOVA-test-selected CNA probes. Of these, only AKT1 copy number gain was associated with poor clinical outcome in BL2 with HR = 3.89 (P = 0.02; Fig. 4C; Supplemental Fig. 7). AKT1 is located on chromosome 14q32.3, a region of gain in multiple BL2 TNBC lesions (Fig. 4B). CNAs in neighbouring genes such as CDCA4, GPR132 and SIVA also predicted poor clinical outcome in BL2 (Supplemental Fig. 8A–C).
We next determined the effect of AKT1 mRNA expression on clinical outcome in BL2 TNBC. Using the 205 TNBC cohort with 14 BL2 tumors, high AKT1 expression identified patients with extremely poor prognosis compared with AKT1-low BL2 samples with HR = 3.93 (P = 0.022; Fig. 4D). The validation cohort with 190 TNBCs contained 10 BL2 samples that could be segregated on the basis of AKT1 expression with HR = 6.1 (P = 0.032; Fig. 4E). The median latency in the 14 BL2 and 10 BL2 cohorts was 38 and 43.4 months for high AKT1 expression versus 73.65 months and undefined period for low AKT1 expression. Thus, AKT1 copy number gain and high mRNA expression can be used to detect a subgroup of BL2 TNBCs with exceedingly poor prognosis. These tumors are likely sensitive to AKT1 inhibitors.
CDCA4 showed a similar profile as AKT1 in BL2 with significant HR for CNA, but its mRNA level stratified BL2 patients only in one of two cohorts (Supplemental Fig. 8A). Notably, expression of AKT1 and CDCA4 was highest in LAR and BL1, respectively, but they showed significant prognostic power only in BL2, in which copy number gain occurs. Several other genes in the chr14q32.3 amplicon in BL2 such as CLBA1 also predicted poor clinical outcome in the test cohort but only showed a trend toward significance in the validation cohort (Supplemental Fig. 8B,C). Thus, high copy number and mRNA expression of AKT1 are the only alterations that consistently showed significant prognostic power in all BL2 cohorts (CNA, mRNA).
Expression analysis of immune-checkpoint blockade and inflammation genes identifies PD1 as prognostic marker for IM TNBC
As noted, some TNBC patients respond well to the anti-PD-L1 immune-checkpoint blockade (ICB) inhibitor, Atezolizumab, in combination with conventional chemotherapy8, yet, most patients, including those with PD-L1+ tumors, succumb to the disease. Here we asked whether ICB proteins or other inflammatory chemokines and cytokines are elevated in specific TNBC subtypes. We first examined their expression in PTEN-low/miRNA-low (group “a”) BL1 TNBCs versus all other TNBC, ER+ and HER2+ breast cancer samples (Supplemental Fig. 9A). The MYC pathway, PD-L1 (CD274) and CD47 were elevated in TNBC, particularly in PTEN-low/miRNA-low TNBCs (group “a”), relative to HER2+ or ER+ breast cancer in most cohorts (Supplemental Fig. 9B–D). Notably, MYC transcriptionally regulates expression of PD-L1 (CD274) and CD4729. CCL7 and CSF2 exhibited a significant increase in expression in group “a” vs TNBC, ER+ and HER2+ breast cancer samples, whereas CXCL8, ICAM1, CCL8, IGFBP3, CXCL1, CCL13, CCL2 and IL6 were elevated in group “a” relative to ER+ and HER2+ lesions (Supplemental Fig. 9E).
Within TNBC subtypes, multiple ICB and inflammatory genes exhibited subtype specificity (Fig. 5A). The ICB genes PDCD1 (PD1) and CTLA4 were highly expressed in most IM TNBC samples. Importantly, PDCD1 expression predicted poor clinical outcome in two independent cohorts with HRs of 5.27 (P = 0.011) (Fig. 5B) and 3.55 (P = 0.0037) (Fig. 5C). Median latency for high PDCD1 IM patients in these cohorts was 124 and 104 months compared with undefined period or 241.5 months for IM patients expressing low level of PDCD1, respectively – a difference of over 10 years. In contrast, PD1 expression had little prognostic power in all TNBC lesions indicating that other factors in the immune microenvironment likely impact outcome in the other subtypes (see Discussion). High CTLA4 expression showed HR of 4.19 but the P value only showed a trend towards significance (P = 0.072) possibly due to a small number of patients in the CTLA4-negative arm (Supplemental Fig. 10A). In addition, CCL2 was highly expressed in most MSL, whereas, CXCL8 was elevated in all BL2 but with no prognostic value (Supplemental Fig. 10B,C).
Discussion
The development of effective therapy for TNBC is hampered in part by the high level of heterogeneity in these tumors. Indeed, TNBC patients are pathologically defined in negative terms by their low expression of ER, PR and HER2, not positively via specific markers. BRCA1/2 mutations and BrCAness are a hallmark of familial basal-like breast cancer, but sporadic TNBCs are poorly defined in the clinic. This low information state has been improved by molecular classifications that initially defined basal and claudin-low TNBC and more recently six different TNBC subtypes: basal-like 1 (BL1), basal-like 2 (BL2), mesenchymal (M), mesenchymal stem–like (MSL), immunomodulatory (IM), and luminal androgen receptor (LAR)3. Each of these TNBC subgroups shows a wide range of survival over 20 years from initial detection, suggesting either that (i) actual death may be stochastic or affected by extrinsic factors such as overall health status and treatment regimens (e.g. successful/timely removal of primary tumor prior to dissemination), or (ii) each subtype is inherently heterogeneous and driven by different oncogenic alterations within each subtype. Here we provide evidence for the latter scenario; we show that BL1, BL2 and IM TNBC can be stratified on the basis of specific oncogenic alterations into poor and relatively moderate subgroups. High risk BL1 tumors show multiple alterations including high WNT, MYC, and RhoA signalling and loss of RB1 and PTEN plus p53 mutation rather than deletion; high risk BL2 lesions exhibit AKT1 copy number gain and high AKT1 mRNA expression, whereas high risk IM patients show high mRNA expression of the ICB gene PDCD1 (PD1).
The difference in survival between TNBCs with these alterations compared with those without them is remarkable. For BL1, the average median difference is <32 vs >204 months for group “a” or ~36 vs >204 for Pten-low/RhoA signaling high (a greater than 14 year difference); for BL2 it is <41 vs >74 months for AKT1 high (>2.75 year difference); and for IM, 114 vs >242 months for high PD1 (>10.67 year difference). Thus, the wide range of survival in each subtype of TNBC patients is driven by intrinsic heterogeneity and specific oncogenic alterations. Upon further validation of these findings, these patients can be identified and prioritized for therapy.
The oncogenic alterations that stratify BL1, BL2 and IM offer potential targets for therapeutic intervention. Thus, high risk BL1 may be treated with inhibitors of the PI3K pathway, WNT, MYC, RhoA and/or RB1 downstream effectors30; BL2 with AKT1 inhibitors, and IM with PD1/PD-L1 based therapy. Further analysis may identify oncogenic alterations that can stratify M, MSL and LAR TNBC subtypes. We also found specific alterations in each subtype that did not stratify these patients into different risk groups, but may serve as excellent targets for therapy. Moreover, these genes are candidates for the creation of new genetically engineered mouse models for testing TNBC stratification and aggressiveness, as well as novel therapies.
Genetically engineered mouse models for each specific TNBC subtype would greatly facilitate not only the analysis of subtype specific vulnerabilities but also provide immune-competent hosts to test potential ICB-based therapeutics. BL1 TNBC is modeled by targeted deletion of RB1 or PTEN together with p53 deletion or mutation19,30,31,32. Given the combined loss of all three tumor-suppressor genes in aggressive BL1 (this study), the generation of triple Rb/Pten/p53 triple mutant mouse model in the mammary gland is warranted. Notably, recent preclinical analysis in several laboratories point to the CHK1/WEE1/CDC25/Aurora Kinase pathway as a target for RB1/PTEN/TP53-deficient basal and mesenchymal TNBC subtypes24,33.
For BL2, a logical model would involve Akt1/PKBα overexpression on a p53 mutant background. Notably, expression of activated PKB/Akt1 in mammary epithelium provides a survival signal that cooperates with a mutant middle T antigen, defective in PI3K signaling, to promote tumor growth but does not promote metastatic progression34. Similarly, p53-R172H mutation increased mammary tumors in transgenic myristylated-AKT mice, but did not promote full tumor formation35. This activated AKT1 allele induces senescence, a barrier for tumorigenesis. Perhaps a BL2 TNBC model would require p53 mutation plus high wild-type AKT1/PKB expression either with other genes on the 14q32.3 amplicon (through chromosome engineering) or in combination with other oncogenic alterations that occur in BL2.
The tumor immune microenvironment is controlled by complex interaction with tumor cells and can be exploited therapeutically36,37,38. Indeed, one mechanism by which tumors evade the immune system involves suppression of cytotoxic T cells through expression of Programmed Death-Ligand 1 (PD-L1; cluster of differentiation 274, CD274). PD-L1 binds Programmed cell death protein 1 (PD-1, CD279) on the surface of cytotoxic CD8+ T cells, shutting down their ability to kill tumor cells. High expression of PD-L1 on tumor cells or in the TME, and expression of PD-1 receptor on immune cells, inversely correlates with survival of BC patients39. PD-L1 also inhibits CD4 T helper cells and promotes survival of regulatory anti-inflammatory, suppressive T regulatory (Treg) cells. As noted, a phase III trial on TNBC patients with PD-L1 positive tumors revealed a median OS of 27 months following anti-PD-L1 (Atezolizumab) plus nab-Paclitaxel treatment, compared with 15.5 months for nab-Paclitaxel only therapy8. While these are seminal results, the response was not durable and many patients, including those with high PD-L1 expression, showed only marginal response, and ultimately succumbed to their disease, albeit at a slower pace. This raises the question of whether specific subgroups of TNBC patients may respond better to anti-PD-L1 therapy.
A recent analysis of immune microenvironment of TNBC identified four groups on the basis of CD8+ T cells in the tumor core, stroma and epithelial compartments as well as expression of other immune determinants such as PD-L1, uncovering extensive complexity40. Here, we found that high expression of the PD-L1 receptor, PD1, alone identifies IM TNBC patients at high risk who succumb, on average, over 10 years sooner than patients with PD1-low tumors. The effect was not seen when all TNBC patients were considered, indicating that different immune factors impact survival in other TNBC subtypes. Our results suggest that IM patients with high PD1 expression may benefit from anti-PD1 or PD-L1 therapy. Interestingly, TGFβ was shown to attenuate tumour response to PD-L1 blockade by trapping cytotoxic T cells in the stroma, excluding them from the tumor parenchyma41,42. In this regard, our observation of high TGFβ signalling in BL2 TNBC points to these patients as excellent candidates for anti-PD-L1 plus anti-TGFβ combination therapy.
In conclusions, we report that TNBC subtypes, in particular BL1, BL2 and IM, can be further stratified on the basis of intrinsic oncogenic alterations. The effect of stratification is considerable with survival differences of more than 14, 2.75 and 10.5 years between high and low risk BL1, BL2 and IM patients, respectively. These poorest outcome subgroups within each subtype can now be identified and prioritized for aggressive, personalized therapy. The oncogenic alterations used to stratify these TNBC subtypes can be further used as targets for therapeutic interventions.
Methods
Datasets
Training set included 1302 BC/205 TNBCs from public mRNA dataset EGAD00010000434 and miRNA dataset EGAD00010000438 were obtained from EGAS0000000012243; CNA dataset was obtained from EGAS0000100175344 in which, 185 TNBCs have mutation data. In the validation set, a cohort containing 207 BC/44 TNBC from mRNA dataset GSE22219 and miRNA dataset GSE22216 were obtained from GSE2222045 and used to verify PTEN/miRNA subtypes. In addition, another validation cohort containing 146 TNBCs with mRNA and mutation data was obtained from METABRIC (https://www.cbioportal.org/study/summary?id=brca_metabric). mRNA data from the integrated 190 TNBC in the two validation cohorts were used to confirm TNBC subtypes expression in the training cohort.
TNBC subtypes and pathway signatures
PTEN/miRNA groups23, TNBCType-63, TNBCType-46, PAM50 + Claudin2 and IntCluster22 subtypes were classified according to the methods in each original publication. Pathway activity of 25 oncogenic and tumor suppressor gene signatures25,26 and RB1-loss signature24 were calculated as described.
Statistical analysis
Statistical analysis was performed using Prism Software (GraphPad Software, La Jolla, CA, USA). Differences in gene expression, pathway and signature activity, CNA, and mutation status between different groups were calculated using ANOVA and student t-test. For Kaplan-Meier survival analysis log-rank (Mantel-Cox) test was used to calculate P values and hazard ratios (HRs), and to compare survival curves. P value < 0.05 was considered significant and marked as *. P < 0.01, P < 0.001 and P < 0.0001 were marked as **, *** and ****, respectively. P > 0.05 was considered not significant (ns).
Data availability
All data generated and/or analyzed during this study are referenced or included in this published article.
References
Prat, A. et al. Molecular characterization of basal-like and non-basal-like triple-negative breast cancer. The oncologist 18, 123–133, https://doi.org/10.1634/theoncologist.2012-0397 (2013).
Prat, A. et al. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res 12, R68 (2010).
Lehmann, B. D. et al. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J Clin Invest 121, 2750–2767, https://doi.org/10.1172/JCI45014 (2011).
Millikan, R. C. et al. Epidemiology of basal-like breast cancer. Breast Cancer Res Treat 109, 123–139 (2008).
Stark, A. et al. African ancestry and higher prevalence of triple-negative breast cancer: findings from an international study. Cancer 116, 4926–4932, https://doi.org/10.1002/cncr.25276 (2010).
Lehmann, B. D. et al. Refinement of Triple-Negative Breast Cancer Molecular Subtypes: Implications for Neoadjuvant Chemotherapy Selection. PLoS One 11, e0157368, https://doi.org/10.1371/journal.pone.0157368 (2016).
Golubnitschaja, O. et al. Multi-omic approach decodes paradoxes of the triple-negative breast cancer: lessons for predictive, preventive and personalised medicine. Amino Acids 50, 383–395, https://doi.org/10.1007/s00726-017-2524-0 (2018).
Schmid, P. et al. Atezolizumab and Nab-Paclitaxel in Advanced Triple-Negative Breast Cancer. N Engl J Med, https://doi.org/10.1056/NEJMoa1809615 (2018).
Garrido-Castro, A. C., Lin, N. U. & Polyak, K. Insights into Molecular Classifications of Triple-Negative. Breast Cancer: Improving Patient Selection for Treatment. Cancer Discov 9, 176–198, https://doi.org/10.1158/2159-8290.CD-18-1177 (2019).
Turner, N. C. & Reis-Filho, J. S. Basal-like breast cancer and the BRCA1 phenotype. Oncogene 25, 5846–5853 (2006).
Jung, J. et al. Association between basal-like phenotype and BRCA1/2 germline mutations in Korean breast cancer patients. Curr Oncol 23, 298–303, https://doi.org/10.3747/co.23.3054 (2016).
Lim, E. et al. Aberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers. Nat Med 15, 907–913 (2009).
Turner, N., Tutt, A. & Ashworth, A. Hallmarks of ‘BRCAness’ in sporadic cancers. Nat Rev Cancer 4, 814–819, https://doi.org/10.1038/nrc1457 (2004).
Geenen, J. J. J., Linn, S. C., Beijnen, J. H. & Schellens, J. H. M. PARP Inhibitors in the Treatment of Triple-Negative Breast Cancer. Clin Pharmacokinet 57, 427–437, https://doi.org/10.1007/s40262-017-0587-4 (2018).
Lord, C. J., Tutt, A. N. & Ashworth, A. Synthetic lethality and cancer therapy: lessons learned from the development of PARP inhibitors. Annu Rev Med 66, 455–470, https://doi.org/10.1146/annurev-med-050913-022545 (2015).
Curtis, C. et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature 486, 346–352, https://doi.org/10.1038/nature10983 (2012).
Koboldt, D. C. et al. Comprehensive molecular portraits of human breast tumours. Nature 490, 61–70, https://doi.org/10.1038/nature11412 (2012).
Nik-Zainal, S. et al. Landscape of somatic mutations in 560 breast cancer whole-genome sequences. Nature 534, 47–54, https://doi.org/10.1038/nature17676 (2016).
Jones, R. A. et al. RB1 deficiency in triple-negative breast cancer induces mitochondrial protein translation. J Clin Invest 126, 3739–3757, https://doi.org/10.1172/JCI81568 (2016).
Robinson, D. R. et al. Integrative clinical genomics of metastatic cancer. Nature, https://doi.org/10.1038/nature23306 (2017).
Wang, S. et al. microRNA-143/145 loss induces Ras signaling to promote aggressive Pten-deficient basal-like breast cancer. JCI Insight 2, https://doi.org/10.1172/jci.insight.93313 (2017).
Rueda, O. M. et al. Dynamics of breast-cancer relapse reveal late-recurring ER-positive genomic subgroups. Nature 567, 399–404, https://doi.org/10.1038/s41586-019-1007-8 (2019).
Wang, D. Y., Gendoo, D. M. A., Ben-David, Y., Woodgett, J. R. & Zacksenhaus, E. A subgroup of microRNAs defines PTEN-deficient, triple-negative breast cancer patients with poorest prognosis and alterations in RB1, MYC, and Wnt signaling. Breast Cancer Res 21, 18, https://doi.org/10.1186/s13058-019-1098-z (2019).
Liu, J. C. et al. Identification of CDC25 as a Common Therapeutic Target for Triple-Negative Breast Cancer. Cell Rep 23, 112–126, https://doi.org/10.1016/j.celrep.2018.03.039 (2018).
Andrechek, E. R. et al. Genetic heterogeneity of Myc-induced mammary tumors reflecting diverse phenotypes including metastatic potential. Proc Natl Acad Sci USA 106, 16387–16392, https://doi.org/10.1073/pnas.0901250106 (2009).
Gatza, M. L. et al. A pathway-based classification of human breast cancer. Proc Natl Acad Sci USA 107, 6994–6999 (2010).
Cancer Genome Atlas, N. Comprehensive molecular portraits of human breast tumours. Nature 490, 61–70, https://doi.org/10.1038/nature11412 (2012).
Choucair, K. A. et al. The 16p13.3 (PDPK1) Genomic Gain in Prostate Cancer: A Potential Role in Disease Progression. Transl Oncol 5, 453–460, https://doi.org/10.1593/tlo.12286 (2012).
Casey, S. C. et al. MYC regulates the antitumor immune response through CD47 and PD-L1. Science 352, 227–231, https://doi.org/10.1126/science.aac9935 (2016).
Liu, J. C. et al. Combined deletion of Pten and p53 in mammary epithelium accelerates triple-negative breast cancer with dependency on eEF2K. EMBO Mol Med 6, 1542–1560, https://doi.org/10.15252/emmm.201404402 (2014).
Jiang, Z. et al. Rb deletion in mouse mammary progenitors induces luminal-B or basal-like/EMT tumor subtypes depending on p53 status. J Clin Invest 120, 3296–3309, https://doi.org/10.1172/JCI41490 (2010).
Wang, S., Liu, J. C., Kim, D., Datti, A. & Zacksenhaus, E. Targeted Pten deletion plus p53-R270H mutation in mouse mammary epithelium induces aggressive claudin-low and basal-like breast cancer. Breast Cancer Res 18, 9, https://doi.org/10.1186/s13058-015-0668-y (2016).
Knudsen, E. S. & Zacksenhaus, E. The vulnerability of RB loss in breast cancer: Targeting a void in cell cycle control. Oncotarget 9, 30940–30941, https://doi.org/10.18632/oncotarget.25797 (2018).
Hutchinson, J., Jin, J., Cardiff, R. D., Woodgett, J. R. & Muller, W. J. Activation of Akt (protein kinase B) in mammary epithelium provides a critical cell survival signal required for tumor progression. Mol Cell Biol 21, 2203–2212, https://doi.org/10.1128/MCB.21.6.2203-2212.2001 (2001).
Blanco-Aparicio, C. et al. Exploring the gain of function contribution of AKT to mammary tumorigenesis in mouse models. PLoS One 5, e9305, https://doi.org/10.1371/journal.pone.0009305 (2010).
Sanmamed, M. F. & Chen, L. A Paradigm Shift in Cancer Immunotherapy: From Enhancement to Normalization. Cell 175, 313–326, https://doi.org/10.1016/j.cell.2018.09.035 (2018).
Wellenstein, M. D. & de Visser, K. E. Cancer-Cell-Intrinsic Mechanisms Shaping the Tumor Immune Landscape. Immunity 48, 399–416, https://doi.org/10.1016/j.immuni.2018.03.004 (2018).
Sun, C., Mezzadra, R. & Schumacher, T. N. Regulation and Function of the PD-L1 Checkpoint. Immunity 48, 434–452, https://doi.org/10.1016/j.immuni.2018.03.014 (2018).
Sun, S. et al. PD-1(+) immune cell infiltration inversely correlates with survival of operable breast cancer patients. Cancer immunology, immunotherapy: CII 63, 395–406, https://doi.org/10.1007/s00262-014-1519-x (2014).
Gruosso, T. et al. Spatially distinct tumor immune microenvironments stratify triple-negative breast cancers. J Clin Invest 129, 1785–1800, https://doi.org/10.1172/JCI96313 (2019).
Mariathasan, S. et al. TGFbeta attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature 554, 544–548, https://doi.org/10.1038/nature25501 (2018).
Tauriello, D. V. F. et al. TGFbeta drives immune evasion in genetically reconstituted colon cancer metastasis. Nature 554, 538–543, https://doi.org/10.1038/nature25492 (2018).
Dvinge, H. et al. The shaping and functional consequences of the microRNA landscape in breast cancer. Nature 497, 378–382, https://doi.org/10.1038/nature12108 (2013).
Pereira, B. et al. The somatic mutation profiles of 2,433 breast cancers refines their genomic and transcriptomic landscapes. Nat Commun 7, https://doi.org/10.1038/ncomms11479 (2016).
Buffa, F. M. et al. microRNA-associated progression pathways and potential therapeutic targets identified by integrated mRNA and microRNA expression profiling in breast cancer. Cancer Res 71, 5635–5645, https://doi.org/10.1158/0008-5472.CAN-11-0489 (2011).
Acknowledgements
This work was funded by a research grant from the Canadian Cancer Society/Canadian Breast Cancer Foundation to EZ, and a Terry Fox Foundation Program Grant to JRW and EZ.
Author information
Authors and Affiliations
Contributions
D.Y.W. performed all the bioinformatic analysis; D.Y.W. and E.Z. conceived the idea; D.Y.W., J.Z., Y.B.D., J.R.W. and E.Z. reviewed the data; all authors made substantial contributions to writing and revising the manuscript, and final approval.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wang, DY., Jiang, Z., Ben-David, Y. et al. Molecular stratification within triple-negative breast cancer subtypes. Sci Rep 9, 19107 (2019). https://doi.org/10.1038/s41598-019-55710-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-55710-w
This article is cited by
-
Cilengitide sensitivity is predicted by overall integrin expression in breast cancer
Breast Cancer Research (2024)
-
Mesenchymal–epithelial transition and AXL inhibitor TP-0903 sensitise triple-negative breast cancer cells to the antimalarial compound, artesunate
Scientific Reports (2024)
-
Understanding the chemistry & pharmacology of antibody–drug conjugates in triple-negative breast cancer with special reference to exatecan derivatives
Medical Oncology (2024)
-
CDK4/6 inhibitors and the pRB-E2F1 axis suppress PVR and PD-L1 expression in triple-negative breast cancer
Oncogenesis (2023)
-
The potential roles of lncRNA TINCR in triple negative breast cancer
Molecular Biology Reports (2023)