- Research article
- Open Access
- Open Peer Review
BRCA1-mutated and basal-like breast cancers have similar aCGH profiles and a high incidence of protein truncating TP53 mutations
BMC Cancervolume 10, Article number: 654 (2010)
Basal-like breast cancers (BLBC) are aggressive breast cancers for which, so far, no targeted therapy is available because they typically lack expression of hormone receptors and HER2. Phenotypic features of BLBCs, such as clinical presentation and early age of onset, resemble those of breast tumors from BRCA1-mutation carriers. The genomic instability of BRCA1-mutated tumors can be effectively targeted with DNA-damaging agents and poly-(ADP-ribose) polymerase 1 (PARP1) inhibitors. Molecular similarities between BLBCs and BRCA1-mutated tumors may therefore provide predictive markers for therapeutic response of BLBCs.
There are several known molecular features characteristic for BRCA1-mutated breast tumors: 1) increased numbers of genomic aberrations, 2) a distinct pattern of genomic aberrations, 3) a high frequency of TP53 mutations and 4) a high incidence of complex, protein-truncating TP53 mutations. We compared the frequency of TP53 mutations and the pattern and amount of genomic aberrations between BRCA1-mutated breast tumors, BLBCs and luminal breast tumors by TP53 gene sequencing and array-based comparative genomics hybridization (aCGH) analysis.
We found that the high incidence of protein truncating TP53 mutations and the pattern and amount of genomic aberrations specific for BRCA1-mutated breast tumors are also characteristic for BLBCs and different from luminal breast tumors.
Complex, protein truncating TP53 mutations in BRCA1-mutated tumors may be a direct consequence of genomic instability caused by BRCA1 loss, therefore, the presence of these types of TP53 mutations in sporadic BLBCs might be a hallmark of BRCAness and a potential biomarker for sensitivity to PARP inhibition. Also, our data suggest that a small subset of genomic regions may be used to identify BRCA1-like BLBCs. BLBCs share molecular features that were previously found to be specific for BRCA1-mutated breast tumors. These features might be useful for the identification of tumors with increased sensitivity to (high-dose or dose-dense) alkylating agents and PARP inhibitors.
Lobules and ducts within the normal human breast are lined with a double layer of epithelial cells: an inner layer of luminal cells and an outer layer of basal/myoepithelial cells that are in direct contact with the basement membrane. Transformation of different mammary epithelial cells results in considerable heterogeneity in breast cancer subtypes, giving rise to differences in clinical presentation, histology and response to therapy . Gene expression profiling has identified five molecular breast tumor subtypes: luminal A, luminal B, normal breast-like, human epidermal growth factor 2 (HER2/ERBB2) positive, and basal-like [2, 3]. In the clinic however, only immunohistochemistry data for estrogen receptor (ER), progesterone receptor (PR) and HER2 status are used to guide treatment choice .
Approximately ~70-80% of all breast tumors are hormone receptor positive  and therefore sensitive to adjuvant endocrine therapy . These tumors classify mostly as luminal A/B breast tumors . The addition of trastuzumab to adjuvant chemotherapy considerably improved the outcome of HER2-positive breast tumors . However, approximately 10-20% of all breast tumors do not express hormone receptors or HER2, and are therefore insensitive to endocrine or trastuzumab treatment. Currently, the only treatment available for these triple-negative breast cancers (TNBCs) is cytotoxic chemotherapy .
The TNBC group as defined by immunohistochemical staining consists for approximately 80% of basal-like breast cancers (BLBCs) as defined by gene expression profiling . BLBCs express luminal (CK19 and CK18) as well as basal cytokeratins (CK5/6, CK17 and CK14), suggesting that these tumors originate from an undifferentiated, dual-lineage stem/progenitor cell type. Furthermore, the TP53 gene is often mutated in BLBCs , and the gene expression profiles of TP53-mutated breast tumors show strong association with the BLBC subgroup . Consequently, BLBCs are aggressive tumors with an expansive growth pattern (pushing margins), a high proliferation rate, high relapse rates and poor survival. Moreover, BLBCs occur more frequently in premenopausal women than in postmenopausal women and are often difficult to detect by mammography or ultrasound .
Although TNBC/BLBC has a relatively poor prognosis in the first five years after diagnosis, approximately 60% of patients - even without adjuvant chemotherapy - do not relapse, and after ~8 years of follow up have a high chance of being cured (reviewed in ). This indicates that within the TNBC/BLBC tumor group there is considerable heterogeneity in tumor behavior. At present however, all TNBC/BLBC patients are treated with cytotoxic adjuvant chemotherapy because there are no clinically useful prognostic and predictive markers to identify patients with aggressive, chemotherapy-sensitive TNBC/BLBC, leading to unnecessary exposure to chemotherapy of a substantial number of patients . In recent years it is becoming clear that phenotypic features of TNBC/BLBC may also apply to the majority of hereditary BRCA1-mutated breast tumors [12, 13]. Since BRCA1 function is required for homology-directed repair of DNA double-strand breaks (DSBs), BRCA1-mutated tumors and BRCA1-like BLBCs are predicted to be sensitive to DSB inducing therapy . Indeed, breast tumors from BRCA1-mutation carriers are sensitive to inhibition of DNA single-strand break (SSB) repair by poly(ADP-ribose) polymerase (PARP) inhibitors  and to chemotherapy that causes DSBs, such as platinum drugs, alkylating agents and topoisomerase I poisons [16, 17]. It will therefore be important to identify features of sporadic BLBCs  that may be useful as predictive biomarkers for response to DSB-inducing chemotherapy or PARP inhibitors. Known molecular features characteristic for BRCA1-mutated breast tumors are 1) a high degree of genomic instability due to homologous recombination (HR) deficiency , 2) a distinct pattern of genomic aberrations [19–22] 3) a high frequency of TP53 mutations and 4) a high incidence of complex, protein-truncating TP53 mutations [10, 23]. In this study, we determined to what extent these molecular characteristics of BRCA1-mutated tumors are present in BLBCs.
Breast tumor groups
To compare molecular characteristics of BRCA1-mutated tumors with BLBCs, we used data from two published tumor sets from the Netherlands Cancer Institute that were sequenced for TP53 and for which aCGH data was generated. The first dataset contains 27 BRCA1-mutated breast tumors and 17 luminal breast tumors (defined by expression profiling) that were previously described by Joosse et al . The luminal tumors from this study were designated luminal-J: J for Joosse. The second dataset from Horlings et al.,  contained 21 non-hereditary BLBCs and 31 luminal breast tumors that were part of a series of 295 breast tumor specimens  which were assigned to breast cancer subgroups according tot their gene expression profiles [2, 3]. For privacy reasons BRCA1-mutation status was not verified in the BLBCs, however, these patients did not have a family history of breast cancer. Luminal tumors from this study were labeled luminal-H, H for Horlings. Six luminal breast tumors were used in both the Joosse study and the Horlings study. TP53 mutation, ER/PR/HER2, CK56 staining pattern and age at diagnosis of all tumors are shown in Table 1.
TP53 mutation analyses
For 21 BRCA1-mutated tumors, and 13 luminal-J tumors exons 2-9 of TP53 were previously sequenced . The abundance of the aberrant base was estimated from the sequence chromatogram from both the forward and reverse sequencing runs. When comparing mutation types found in the tumor groups, the influence of tumor heterogeneity was minimized by only including TP53 mutations that had an estimated abundance of >25% in the tumor DNA .
For all BLBC and luminal-H tumors TP53 exons 2-11 (including +/- 30 bp outside each exon) were sequenced using AB 3730 DNA Analyzer (reference sequence NM_000546).
Frameshift, splice and nonsense mutations and in-frame insertions/deletions are defined as "complex TP53 mutations". All missense mutations found in the BRCA1-mutated and luminal-J tumor groups were classified according to their predicted effect on p53 function as determined by the Sorting Intolerant from Tolerant algorithm SIFT; [27, 28], as used in the IARC TP53 database [29, 30]. Because no matched normal/germ-line DNA was available, some benign germ-line variants may have been identified as deleterious somatic mutations by SIFT. All TP53 missense mutations found in the BLBC and luminal-H tumor groups were classified to be deleterious or non-deleterious according to their predicted effect on TP53 function using "EffectGroup3"  as used in the IARC TP53 database. The 29 most common hotspot mutations (P < 0.001) identified by Walker and colleagues  are referred to as ''hotpot mutations'': K132, C135, P151, V157, R158, Y163, V173, R175, C176, H179, H193, Y205, Y220, Y234, M237, C238, S241, C242, G245, M246, R248, R249, G266, R273, P278, R280, D281, R282, and E285.
Comparisons aCGH data derived from FFPE and fresh frozen tissue
For aCGH procedures of the BRCA1-mutated and luminal-J tumor groups: see Joosse et al., 2009 . For aCGH procedures of the BLBC and luminal-H tumor groups: see Horlings et al., . For the aCGH analyses, a microarray platform containing 3,500 human BAC/PAC clones covering the whole genome with an average spacing of 1 Mb was used . Although the 1 Mb resolution of the BAC platform limits sensitivity for focal changes, the aCGH data is a sound representation of our tumor groups and can be used to find recurrent differences between tumor groups.
When comparing aCGH profiles of tumor DNA isolated from formalin-fixed paraffin-embedded (FFPE) tissue  and fresh-frozen tissue , we noticed that the log2-ratios obtained from the different DNA sources showed different distributions. From the six overlapping samples between the luminal-H and luminal-J tumor groups, as shown in Table 1 we saw that distribution of log2-ratios derived from FFPE samples was consistently wider than the log2-ratios derived from fresh frozen tissue (Additional File 1). Therefore, we transformed both the FFPE (i.e. all log2 ratios from BRCA1-mutated and luminal-J tumors taken together) and the fresh frozen datasets (i.e. all log2 ratios from BLBC and luminal-H tumors taken together) to have a mean of zero and a standard deviation of one before applying KC-SMART and comparative-KC-SMART. This enabled us to construct one luminal tumor group consisting of tumors from both the luminal-H and luminal-J tumor groups to compare the BRCA1-mutated and BLBC-data with. The six tumors that overlapped between the two groups were included only once and taken from the luminal-J tumor group.
KC-SMART and comparative-KC-SMART analysis
(Kernel Convolution - a Statistical Method for Aberrant Region detection) is a computational approach for statistical analysis of non-discretized aCGH data from multiple experiments, and determines which regions are significantly gained or lost relative to randomized data (P < 0.05) . We used KC-SMART to smooth the raw log2-ratios by generating a Kernel Smoothed Estimate (KSE) curves for gains (KSEgains) and losses (KSElosses) separately across a group of tumors.
 detects genomic regions that have a differential aCGH signal between two tumor groups. The comparative-KC-SMART algorithm smoothes raw log2-ratio data from each individual tumor profile by placing Gaussian kernel functions with the height of the log2-ratio at the genomic midposition of each probe (without separating gains and losses, as done for KC-SMART). For each tumor, an aggregated profile is determined by convolution of locally weighted kernel functions. For each genomic position, the KSE values from all tumors that belong to the two tumor groups in the comparison are used to calculate a signal to noise ratio (SNR). We determined a cut-off that defines significant SNR values by generating SNR data using 6000 class-label permutations and calculating the significance threshold corresponding with a False Discovery Rate (FDR) of 0.05. The width of a kernel applied to each data point determines the extent of smoothing and the size of aberrations detected. Smoothing individual tumors with a kernel width of 20 Mb resulted in aCGH profiles that recapitulated raw aCGH data well; therefore, we chose to use this kernel width for all comparisons. R-packages of KC-SMART and comparative-KC-SMART have been submitted to Bioconductor . We used NCBI Build 36 (Hg 18) for these analyses.
We used the MeV program  to cluster tumor aCGH data. Samples and genes are hierarchically clustered using pearson correlation and complete linkage with leaf ordening.
TP53 mutations in BRCA1-mutated, BLBC and luminal breast cancers
We previously found that nearly all BRCA1-mutated breast tumors had deleterious TP53 mutations due to an increased frequency of truncating frameshift, splice, nonsense mutations and in-frame insertions and deletions . BRCA1-mutated tumors and BLBCs are both basal like TNBCs, characteristics that are different from hormone receptor positive luminal breast tumors (Table 1). Therefore, we were motivated to compare TP53 mutation type and frequency found in non-hereditary BLBCs, in BRCA1-mutated tumors and in luminal breast tumors. TP53 mutation types and other tumor characteristics are listed in Table 1. The TP53 gene from the BRCA1-mutated/luminal-J tumors and the BLBC/luminal-H tumors were sequenced in different labs with slightly different methods (see Methods section). At the cost of reducing the power of this analysis we wanted to make sure we did not introduce a methodical bias in our comparisons, therefore, the TP53 mutation data for luminal-J and luminal-H tumors were not combined.
Similar to the 90% (19/21) of BRCA1-mutated tumors, 95% (20/21) of the BLBCs harbored TP53 mutations, significantly more than the 46% (6/14) of the luminal-J and the 26% (8/31) of the luminal-H tumors (p = 7 × 10-3 and p = 5 × 10-7 respectively, two-tailed Fisher's Exact test), Figure 1 Additional File 2. We next compared TP53 mutation types from BRCA1-mutated tumors and non-hereditary BLBCs with luminal-H and luminal-J tumors. We found that 52% (11/21) of the BRCA1-mutated tumors and 57% (12/21) of the BLBCs have complex/truncating TP53 mutations, which is significantly more than the 8% (1/13) and 7% (2/31) of the luminal-J and luminal-H tumors, respectively (p = 1 × 10-2 and p = 8 × 10-5 respectively; two-tailed Fisher's Exact test, Figure 1 Additional File 2), indicating that this feature is common in hereditary BRCA1-mutated breast tumors and non-hereditary BLBCs. Interestingly, the increase in deleterious missense or hotspot TP53 mutations in BRCA1-mutated or BLBCs (8/21) their respective luminal tumor groups was not significant. Deleterious missense mutations: 11/21 BRCA1-mutated vs. 5/13 luminal-J, p = 5 × 10-1, and 8/21 BLBCs vs. 6/31 luminal-H tumors p = 2 × 10-1, two-tailed Fisher's Exact test). Hotspot mutations: 9/21 BRCA1-mutated vs. 3/13 luminal-J, p = 3 × 10-1, and 6/21 BLBCs vs. 6/31 luminal-H tumors p = 5 × 10-1, two-tailed Fisher's Exact test). Together, these data suggest that the increase of TP53 mutations in the BRCA1-mutated tumors and BLBCs is primarily due to the increase in complex/truncating TP53 mutations.
Comparison of aCGH profiles of BRCA1-mutated tumors, BLBCs and luminal breast cancers
To identify DNA copy number aberrations (CNAs) that occur significantly more often in BRCA1-mutated breast tumors than in BLBCs, we analyzed their aCGH with comparative-KC-SMART, a computational method for detection of genomic regions that have a significantly different aCGH signal between two tumor groups. However, the DNA samples used to acquire aCGH profiles for the BLBC and the BRCA1-mutated tumor groups were isolated from fresh frozen tissue and FFPE material respectively, resulting in differences in log2-ratio distribution of the aCGH profiles (Additional File 1). To account for this difference, we normalized the data as explained in the Methods section. We applied comparative-KC-SMART to the normalized aCGH data of the BRCA1-mutated tumor group and the BLBC group (Figure 2a). Comparative-KC-SMART detected small CNAs on chromosomes 5, 7, 8 and 14 that are significantly different between these tumor groups (for regions and cancer-related genes see Additional File 3). BRCA1-specific losses and BLBC-specific gains on chromosomes 5 and 7 flanked each other and seem to be dependent on each other. Interestingly, the BRCA1-specific chromosome 7 loss encompasses EGFR. The BRCA1-specific gain on chromosome 14 peaks at the T-cell receptor alpha (TCRα) locus and encompasses, among other cancer-related genes, poly(ADP-ribose) polymerase 2 (PARP2), and B-cell lymphoma 2 like 2 (BCL2L2). Whether any of the genes located within these CNAs promote survival of BRCA1-deficient cells remains to be established.
Next, we used comparative-KC-SMART to compared CNAs in the BRCA1-mutated tumors and BLBCs vs. the combined luminal tumors (Figure 2b-c). Differential gains and losses are shown in Additional File 4. We compared CNAs between Luminal-J and luminal-H tumors as a control for merging these tumor groups (Figure 2d): as expected, no differences between the two luminal groups were detected.
Gains on chromosomes 1q and 16p, and the loss on chromosome 16q occur more often in the luminal tumors than in the basal-like/BRCA1-mutated tumors (Figure 2b-c, Table 2). The first chromosome 1q gain peaks at 177.31 Mb; the second 1q gain peaks at 202.65 Mb close to MDM4 (202.81 Mb). Similarly, the 16p gain peaks at 15.7 Mb. The chromosome 16q loss consists of two peaks: the first peak maps to 52.19 Mb, at the BRD7, CYLD, and RBL2/p130 genes. The second peak maps to 79.19 Mb, at one of the most active common fragile sites in the human genome, FRA16 D, associated with the WWOX gene (77.24 Mb). Many of the luminal tumors show a co-occurrence of chromosome 1q gain and 16q loss (Additional File 5).
Impact of TP53 mutation on luminal tumors
Because BLBCs and BRCA1-mutated tumors are almost always TP53-mutated, we investigated whether TP53 mutations are associated with specific CNAs in luminal tumors. We stratified the 31 luminal-H tumors by their TP53 mutation status and used comparative-KC-SMART to compare both tumor groups. Although no differences were detected between CNAs from 8 TP53-mutated and 23 TP53 wild-type luminal breast tumors, quantitative differences in KC-SMART profiles can be observed, Figure 3. TP53-mutated tumors have more overall gains on chromosomes 3q, 6p, 20q, 21q and 22q and losses on chromosomes 2q, 3p, 4p, 4q, 13q, 15q and X. On the other hand, TP53 wild-type tumors have a higher incidence of chromosome 1q gain, and 16q loss. However, these data do not provide evidence for an altered profile of TP53-mutated luminal breast tumors, perhaps because the TP53-mutated tumor group is too small for robust statistical analysis.
Analysis of genomic instability in BRCA1-mutated tumors and BLBCs
To determine the amount of CNAs in the different tumor groups, we used KC-SMART to smooth individual aCGH profiles and counted the amount of CNAs exceeding a range of cutoffs for each tumor separately. We found that the median amount of CNAs of BRCA1-mutated tumors was not different from the amount of CNAs found in BLBCs (Figure 4a). In contrast, we found that the median amount of CNAs is significantly greater in the BRCA1-mutated tumors compared with luminal breast tumors between KSE cutoffs 0.02 and 0.1 and between 0.24 and 0.74 (P < 0.01, two sided t-test Figure 4b). Similarly, the median amount of CNAs in the BLBCs was higher compared with luminal tumors for KSE cutoffs between 0.02 and 0.14 and between 0.28 and 0.96 (Figure 4c). We did not detect differences in the median amount of aberrations between luminal-J breast tumors and luminal-H except for KSE cutoffs 0.2 and 0.24 (Figure 4d).
Overlapping gains and losses specific for BLBCs and BRCA1-mutated tumors
We compared gains and losses specific for the BLBC/BRCA1-mutated tumor groups relative to the luminal tumor groups. The BLBC/BRCA1-mutated tumors harbored overlapping differential gains on chromosomes 1p, 2p, 3q, 6p, 6q, 7q, 8q, 10p, 12p, 13q and 19q and losses on chromosomes 3p, 4p, 5q, 10q, 12q, 14q and 15q. The luminal tumor groups contained overlapping gains on chromosomes 1q, 8p and 16p, and an overlapping loss on chromosome 6q, 11q, 13q and 16q (Table 2).
Clearly, the differentially occurring CNAs detected by comparative-KC-SMART are fully dependent on the tumors included in the groups. However, when peaks of recurrent aberrations of two different analyses map closely together, this could point to a region whose gain or loss is relevant for tumorigenesis. Cancer-related genes that map to overlapping differential gains or losses between BRCA1-mutated tumors and BLBCs vs. luminal tumors are shown in Additional File 6.
Unsupervised hierarchical clustering of breast cancers
The overlapping differential gains and losses of BLBCs and BRCA1-mutated tumors may represent regions that discriminate BRCA1-mutated tumors and BLBCs from luminal tumors. To test this possibility, we performed an unsupervised hierarchical clustering analysis. First, we smoothed each tumor profile with KC-SMART to remove experimental noise. Then, for each of the regional aberrations specific for the BLBC/BRCA1-mutated tumors or the luminal tumors (shown in Table 2), we calculated the mean of all KSE values within the region for all tumors. We used the mean KSE values to perform a hierarchical clustering of samples and regions using complete linkage and pearson correlation (Figure 5). The tumors clustered in two branches: 6 luminal tumors and 47 of the 48 BLBCs/BRCA1-mutated tumors clustered in one branch, whereas 1 BRCA1-mutated tumor and 42 of the 48 luminal-H/J tumors clustered in the other branch. Interestingly, the BRCA1 and BLBC cases do not form separate clusters but mix together, meaning that a limited amount of regions can distinguish BLBCs and BRCA1-mutated tumors from luminal tumors. The fact that the luminal-H and luminal-J tumors are mixed, shows that no unwanted biases are introduced by differences in quality of DNA from FFPE vs. fresh-frozen tumor material. As an additional internal control we have used the 6 luminal tumors for which DNA from FFPE and fresh-frozen tumor tissue was both available, and each of the 6 pairs cluster together. Notably, four of the six TP53 wild-type tumors that clustered within the BLBC/BRCA1-mutated branch were luminal. One of the six luminal tumors that clustered within the BLBC/BRCA1-mutated branch also had a complex TP53 mutation, whereas only two of the luminal tumors clustering in the luminal branch had a complex TP53 mutation.
BLBC is an aggressive subgroup of breast cancers for which, so far, no druggable target has been identified. In recent years it has become clear that phenotypic features of a subset of BLBCs resemble those of hereditary BRCA1-mutated breast cancers  and are distinctly different from the more common luminal breast tumors. BRCA1-mutated breast tumors are HR deficient (HRD) and can therefore be targeted with DNA-damaging agents or PARP inhibitors [38, 39]. Because there is evidence that a substantial fraction of BLBCs have HR pathway defects , we set out to determine which molecular characteristics of BRCA1-mutated tumors are common to BLBCs but not to luminal tumors. Since some of these characteristics are likely to be linked to BRCA1-disfunction or the HRD phenotype, they might provide important leads for discovery of new biomarkers or drug targets.
High incidence of protein truncating TP53 mutations in BLBC and BRCA1-mutated breast cancer
We and others have previously shown that BRCA1-mutated breast tumors exhibit an increased frequency of TP53 mutations due to a selective increase in complex TP53 mutations such as frameshift, nonsense and splice mutations or in-frame insertions/deletions [10, 23]. In this study, we found that almost all BRCA1-mutated tumors and BLBCs are TP53 mutated compared with 25-50% of the luminal-tumors. Interestingly, 52.4% of the BRCA1-mutated breast tumors and 57.1% of the BLBCs have complex TP53 mutations, significantly more compared with ~7% luminal tumors. It has been suggested that the increased incidence of complex TP53 mutations in BRCA1-mutated tumors is a direct consequence of the genomic instability resulting from the DNA repair defect induced by BRCA1 loss [10, 23]. Furthermore, the DSB repair defect of BRCA1-deficient tumors might confer strong selection pressure on mutation of TP53 in order to abrogate the p53-dependent DNA damage checkpoint . The high frequency of TP53 mutations in non-hereditary BLBCs might suggest that these tumors are also compromised in homology-directed DSB repair. To test this possibility, it would be interesting to perform functional assays to measure DNA damage response and DNA repair in non-hereditary BLBCs with a (complex) TP53 mutation. It has previously been suggested that TP53 mutations, including complex TP53 mutations, affecting the DNA binding domain of the p53 protein may cause resistance to several different cytotoxic compounds such as anthracyclins, 5FU and mitomycin [42, 43]. However, tumors used in these studies were primarily invasive ductal carcinomas, most of which are likely not compromised in homology directed DSB repair. It is therefore interesting that Silver et al. recently reported a significant association between truncating TP53 mutations and cisplatin response in TNBCs .
It is interesting that, like us, Manie et al., found an increased frequency of complex TP53 mutations in BRCA1-mutated breast tumors; however, they did not detect this feature in BLBCs [10, 23]. It is possible that BLBCs are more effectively identified by gene expression profiling than by immunohistochemical selection of tumors negative for ER/PR/HER2 and positive for CK5/6, CK14 or EGFR, as done by Manie et al. Of note, the TP53 mutation analyses for the BRCA1-mutated breast tumors and the BLBCs described in this manuscript were performed in different labs, thereby reducing the possibility that the detection of a high incidence of complex/truncating TP53 mutations is due to a methodical artifact.
Recurrent CNAs in luminal tumors
Luminal breast tumors are mostly TP53 wild type, or they have a common TP53 hotspot mutation. Interestingly, the chromosome 1q gain and 16q loss (1q/16q) differentiate luminal tumors from BRCA1-mutated/BLBC tumors. Co-occurrence of these aberrations results from an unbalanced translocation event [45, 46] and has been associated with TP53 wild-type status, low tumor grade and good prognosis . Indeed, many of the ER-positive luminal tumors show the 1q/16q co-occurrence. The peak of the chromosome 1q gain of the luminal tumor group maps to 202.75 Mb, with the MDM4 gene at (202.81 Mb), which is a negative regulator of p53 . The 16q loss has two peaks, the first maps near the bromodomain 7 (BRD7) gene, associated with downregulation of p53 , CYLD, loss of which is associated with oncogenesis by activation of NF-ΚB signaling , and the retinoblastoma-like 2 gene (RBL2/p130), involved in G1S cell cycle control and senescence . The second peak on the 16q loss maps to one of the most active common fragile sites in the human genome, FRA16 D (associated with the WWOX gene), which could underlie the translocation process . Together, these data suggest that, whereas development of BLBCs or BRCA1-mutated tumors depends on TP53 mutation, indirect p53 downregulation may be sufficient for luminal tumor development.
aCGH profiles of BRCA1-mutated breast tumors resemble BLBCs
BRCA1-mutated breast tumors are associated with a specific aCGH profile which exhibits features that can be used to identify hereditary breast tumors for which information on BRCA1-mutation is not available [19, 22, 24]. Several previous studies have reported a specific aCGH profile for BLBCs different from other breast cancer subtypes [53, 54]. Interestingly, comparison of CNAs from BRCA1-mutated tumors and BLBCs using comparative-KC-SMART yielded a limited set BRCA1- or BLBC-specific aberrations. The peaks of most gains and losses of BRCA1-mutated tumors and BLBCs co-localized, suggesting a common selection pressure during development of these tumors. Indeed, we found that BRCA1-mutated tumors and BLBCs showed many overlapping CNAs, including the chromosome 3q gain and the chromosome 5q loss. Importantly, clustering on the basis of these regions separated BLBCs and BRCA1-mutated breast tumors from luminal breast tumors.
BRCA1 and TP53 involvement in the BLBC phenotype
Our data show that BRCA1-mutated tumors share molecular characteristics of undifferentiated BLBCs. It has previously been proposed that BRCA1-mutation is associated with BLBC because BRCA1 function has stem cell regulation properties and because loss of BRCA1 impairs DNA damage repair during epithelial cell differentiation . However, it is also possible that a defect in DNA repair mechanisms is primarily harmful in proliferating cells, which are more prone to acquire genetic lesions during cell division. Notably, proliferating cells in the premenopausal mammary gland have been shown to rarely express hormone receptors whereas hormone receptor-positive cells only rarely divide . In line with this, and in contrast to most breast tumors, both BRCA1-mutated tumors and BLBCs often occur in premenopausal breast epithelia . Therefore, we propose that inadequate DNA repair mechanisms result in increased susceptibility to genomic instability in the proliferating hormone receptor-negative cells of the premenopausal mammary gland. Furthermore, we propose that because of this, hormone receptor-negative cells depend heavily on p53-mediated cell cycle arrest and apoptosis to remain untransformed. This line of thought lends great importance to TP53 mutation in TNBC.
Our data suggest that a small subset of genomic regions may be useful for the identification of BRCA1-like BLBCs, which exhibit a high frequency of TP53 mutations, especially protein truncating mutations. These features of basal-like breast cancers might be useful for the identification of tumors with increased sensitivity to (high-dose or dose-dense) alkylating agents and PARP inhibitors. In support of this, it was recently reported that TP53 mutations in non-inflammatory BLBCs are highly predictive of complete response to dose-dense neoadjuvant chemotherapy with epirubicine-cyclophosphamide . Furthermore, a significant positive correlation was found between truncating TP53 mutations and cisplatin response in TNBCs . Together, these and our data support further investigation of (protein truncating) TP53 mutation status as a potential predictor of chemotherapy responsiveness in solid tumors.
Sims AH, Howell A, Howell SJ, Clarke RB: Origins of breast cancer subtypes and therapeutic implications. Nat Clin Pract Oncol. 2007, 4: 516-525. 10.1038/ncponc0908.
Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS, et al: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001, 98: 10869-10874. 10.1073/pnas.191367098.
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al: Molecular portraits of human breast tumours. Nature. 2000, 406: 747-752. 10.1038/35021093.
Harris L, Fritsche H, Mennel R, Norton L, Ravdin P, Taube S, Somerfield MR, Hayes DF, Bast RC: American Society of Clinical Oncology 2007 update of recommendations for the use of tumor markers in breast cancer. J Clin Oncol. 2007, 25: 5287-5312. 10.1200/JCO.2007.14.2364.
Chu KC, Anderson WF, Fritz A, Ries LA, Brawley OW: Frequency distributions of breast cancer characteristics classified by estrogen receptor and progesterone receptor status for eight racial/ethnic groups. Cancer. 2001, 92: 37-45. 10.1002/1097-0142(20010701)92:1<37::AID-CNCR1289>3.0.CO;2-F.
Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: an overview of the randomised trials. Lancet. 2005, 365: 1687-1717. 10.1016/S0140-6736(05)66544-0.
Smith I, Procter M, Gelber RD, Guillaume S, Feyereislova A, Dowsett M, Goldhirsch A, Untch M, Mariani G, Baselga J, et al: 2-year follow-up of trastuzumab after adjuvant chemotherapy in HER2-positive breast cancer: a randomised controlled trial. Lancet. 2007, 369: 29-36. 10.1016/S0140-6736(07)60028-2.
Schneider BP, Winer EP, Foulkes WD, Garber J, Perou CM, Richardson A, Sledge GW, Carey LA: Triple-negative breast cancer: risk factors to potential targets. Clin Cancer Res. 2008, 14: 8010-8018. 10.1158/1078-0432.CCR-08-1208.
Linn SC, Van 't Veer LJ: Clinical relevance of the triple-negative breast cancer concept: Genetic basis and clinical utility of the concept. European journal of cancer (Oxford, England: 1990). 2009, 45: 11-26.
Manie E, Vincent-Salomon A, Lehmann-Che J, Pierron G, Turpin E, Warcoin M, Gruel N, Lebigot I, Sastre-Garau X, Lidereau R, et al: High frequency of TP53 mutation in BRCA1 and sporadic basal-like carcinomas but not in BRCA1 luminal breast tumors. Cancer Res. 2009, 69: 663-671. 10.1158/0008-5472.CAN-08-1560.
Langerod A, Zhao H, Borgan O, Nesland JM, Bukholm IR, Ikdahl T, Karesen R, Borresen-Dale AL, Jeffrey SS: TP53 mutation status and gene expression profiles are powerful prognostic markers of breast cancer. Breast Cancer Res. 2007, 9: R30-10.1186/bcr1675.
Lakhani SR, Reis-Filho JS, Fulford L, Penault-Llorca F, van der Vijver M, Parry S, Bishop T, Benitez J, Rivas C, Bignon YJ, et al: Prediction of BRCA1 status in patients with breast cancer using estrogen receptor and basal phenotype. Clin Cancer Res. 2005, 11: 5175-5180. 10.1158/1078-0432.CCR-04-2424.
Turner NC, Reis-Filho JS: Basal-like breast cancer and the BRCA1 phenotype. Oncogene. 2006, 25: 5846-5853. 10.1038/sj.onc.1209876.
Turner N, Tutt A, Ashworth A: Hallmarks of 'BRCAness' in sporadic cancers. Nat Rev Cancer. 2004, 4: 814-819. 10.1038/nrc1457.
Fong PC, Boss DS, Yap TA, Tutt A, Wu P, Mergui-Roelvink M, Mortimer P, Swaisland H, Lau A, O'Connor MJ, et al: Inhibition of poly(ADP-ribose) polymerase in tumors from BRCA mutation carriers. N Engl J Med. 2009, 361: 123-134. 10.1056/NEJMoa0900212.
Byrski T, Huzarski T, Dent R, Gronwald J, Zuziak D, Cybulski C, Kladny J, Gorski B, Lubinski J, Narod SA: Response to neoadjuvant therapy with cisplatin in BRCA1-positive breast cancer patients. Breast Cancer Res Treat. 2009, 115: 359-363. 10.1007/s10549-008-0128-9.
Kennedy RD, Quinn JE, Mullan PB, Johnston PG, Harkin DP: The role of BRCA1 in the cellular response to chemotherapy. J Natl Cancer Inst. 2004, 96: 1659-1668. 10.1093/jnci/djh312.
Moynahan ME, Cui TY, Jasin M: Homology-directed dna repair, mitomycin-c resistance, and chromosome stability is restored with correction of a Brca1 mutation. Cancer Res. 2001, 61: 4842-4850.
van Beers EH, van Welsem T, Wessels LF, Li Y, Oldenburg RA, Devilee P, Cornelisse CJ, Verhoef S, Hogervorst FB, van't Veer LJ, Nederlof PM: Comparative genomic hybridization profiles in human BRCA1 and BRCA2 breast tumors highlight differential sets of genomic aberrations. Cancer Res. 2005, 65: 822-827.
Tirkkonen M, Johannsson O, Agnarsson BA, Olsson H, Ingvarsson S, Karhu R, Tanner M, Isola J, Barkardottir RB, Borg A, Kallioniemi OP: Distinct somatic genetic changes associated with tumor progression in carriers of BRCA1 and BRCA2 germ-line mutations. Cancer Res. 1997, 57: 1222-1227.
Jonsson G, Naylor TL, Vallon-Christersson J, Staaf J, Huang J, Ward MR, Greshock JD, Luts L, Olsson H, Rahman N, et al: Distinct genomic profiles in hereditary breast tumors identified by array-based comparative genomic hybridization. Cancer Res. 2005, 65: 7612-7621.
Wessels LF, van Welsem T, Hart AA, van't Veer LJ, Reinders MJ, Nederlof PM: Molecular classification of breast carcinomas by comparative genomic hybridization: a specific somatic genetic profile for BRCA1 tumors. Cancer Res. 2002, 62: 7110-7117.
Holstege H, Joosse SA, van Oostrom CT, Nederlof PM, de Vries A, Jonkers J: High incidence of protein-truncating TP53 mutations in BRCA1-related breast cancer. Cancer Res. 2009, 69: 3625-3633. 10.1158/0008-5472.CAN-08-3426.
Joosse SA, van Beers EH, Tielen IH, Horlings H, Peterse JL, Hoogerbrugge N, Ligtenberg MJ, Wessels LF, Axwijk P, Verhoef S, et al: Prediction of BRCA1-association in hereditary non-BRCA1/2 breast carcinomas with array-CGH. Breast Cancer Res Treat. 2009, 116: 479-489. 10.1007/s10549-008-0117-z.
Horlings HM, Lai C, Nuyten DS, Halfwerk H, Kristel P, van Beers E, Joosse SA, Klijn C, Nederlof PM, Reinders MJ, et al: Integration of DNA copy number alterations and prognostic gene expression signatures in breast cancer patients. Clin Cancer Res. 2010, 16: 651-663. 10.1158/1078-0432.CCR-09-0709.
van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, et al: A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002, 347: 1999-2009. 10.1056/NEJMoa021967.
Ng PC, Henikoff S: Predicting deleterious amino acid substitutions. Genome Res. 2001, 11: 863-874. 10.1101/gr.176601.
SIFT Database: Sorting Intolerant From Tolerant. [http://blocks.fhcrc.org/sift/SIFT.html]
Petitjean A, Mathe E, Kato S, Ishioka C, Tavtigian SV, Hainaut P, Olivier M: Impact of mutant p53 functional properties on TP53 mutation patterns and tumor phenotype: lessons from recent developments in the IARC TP53 database. Hum Mutat. 2007, 28: 622-629. 10.1002/humu.20495.
IARC TP53 Database. [http://www-p53.iarc.fr]
Olivier M, Langerod A, Carrieri P, Bergh J, Klaar S, Eyfjord J, Theillet C, Rodriguez C, Lidereau R, Bieche I, et al: The clinical value of somatic TP53 gene mutations in 1,794 patients with breast cancer. Clin Cancer Res. 2006, 12: 1157-1167. 10.1158/1078-0432.CCR-05-1029.
Walker DR, Bond JP, Tarone RE, Harris CC, Makalowski W, Boguski MS, Greenblatt MS: Evolutionary conservation and somatic mutation hotspot maps of p53: correlation with p53 protein structural and functional features. Oncogene. 1999, 18: 211-218. 10.1038/sj.onc.1202298.
Joosse SA, van Beers EH, Nederlof PM: Automated array-CGH optimized for archival formalin-fixed, paraffin-embedded tumor material. BMC Cancer. 2007, 7: 43-10.1186/1471-2407-7-43.
Klijn C, Holstege H, de Ridder J, Liu X, Reinders M, Jonkers J, Wessels L: Identification of cancer genes using a statistical framework for multiexperiment analysis of nondiscretized array CGH data. Nucleic Acids Res. 2008, 36: e13-10.1093/nar/gkm1143.
Holstege H, van Beers E, Velds A, Liu X, Joosse SA, Klarenbeek S, Schut E, Kerkhoven R, Klijn CN, Wessels LF, et al: Cross-species comparison of aCGH data from mouse and human BRCA1- and BRCA2-mutated breast cancers. BMC Cancer. 2010, 10: 455-10.1186/1471-2407-10-455.
Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, Li J, Thiagarajan M, White JA, Quackenbush J: TM4 microarray software suite. Methods Enzymol. 2006, 411: 134-193. 10.1016/S0076-6879(06)11009-5.
Farmer H, McCabe N, Lord CJ, Tutt AN, Johnson DA, Richardson TB, Santarosa M, Dillon KJ, Hickson I, Knights C, et al: Targeting the DNA repair defect in BRCA mutant cells as a therapeutic strategy. Nature. 2005, 434: 917-921. 10.1038/nature03445.
Ashworth A: A synthetic lethal therapeutic approach: poly(ADP) ribose polymerase inhibitors for the treatment of cancers deficient in DNA double-strand break repair. J Clin Oncol. 2008, 26: 3785-3790. 10.1200/JCO.2008.16.0812.
Han J, Haiman C, Niu T, Guo Q, Cox DG, Willett WC, Hankinson SE, Hunter DJ: Genetic variation in DNA repair pathway genes and premenopausal breast cancer risk. Breast Cancer Res Treat. 2009, 115: 613-622. 10.1007/s10549-008-0089-z.
Brugarolas J, Jacks T: Double indemnity: p53, BRCA and cancer. p53 mutation partially rescues developmental arrest in Brca1 and Brca2 null mice, suggesting a role for familial breast cancer genes in DNA damage repair. Nat Med. 1997, 3: 721-722. 10.1038/nm0797-721.
Geisler S, Borresen-Dale AL, Johnsen H, Aas T, Geisler J, Akslen LA, Anker G, Lonning PE: TP53 gene mutations predict the response to neoadjuvant treatment with 5-fluorouracil and mitomycin in locally advanced breast cancer. Clin Cancer Res. 2003, 9: 5582-5588.
Geisler S, Lonning PE, Aas T, Johnsen H, Fluge O, Haugen DF, Lillehaug JR, Akslen LA, Borresen-Dale AL: Influence of TP53 gene alterations and c-erbB-2 expression on the response to treatment with doxorubicin in locally advanced breast cancer. Cancer Res. 2001, 61: 2505-2512.
Silver DP, Richardson AL, Eklund AC, Wang ZC, Szallasi Z, Li Q, Juul N, Leong CO, Calogrias D, Buraimoh A, et al: Efficacy of Neoadjuvant Cisplatin in Triple-Negative Breast Cancer. J Clin Oncol.
Farabegoli F, Hermsen MA, Ceccarelli C, Santini D, Weiss MM, Meijer GA, van Diest PJ: Simultaneous chromosome 1q gain and 16q loss is associated with steroid receptor presence and low proliferation in breast carcinoma. Mod Pathol. 2004, 17: 449-455. 10.1038/modpathol.3800059.
Tsarouha H, Pandis N, Bardi G, Teixeira MR, Andersen JA, Heim S: Karyotypic evolution in breast carcinomas with i(1)(q10) and der(1;16)(q10;p10) as the primary chromosome abnormality. Cancer Genet Cytogenet. 1999, 113: 156-161. 10.1016/S0165-4608(99)00016-3.
Fridlyand J, Snijders AM, Ylstra B, Li H, Olshen A, Segraves R, Dairkee S, Tokuyasu T, Ljung BM, Jain AN, et al: Breast tumor copy number aberration phenotypes and genomic instability. BMC Cancer. 2006, 6: 96-10.1186/1471-2407-6-96.
Toledo F, Krummel KA, Lee CJ, Liu CW, Rodewald LW, Tang M, Wahl GM: A mouse p53 mutant lacking the proline-rich domain rescues Mdm4 deficiency and provides insight into the Mdm2-Mdm4-p53 regulatory network. Cancer Cell. 2006, 9: 273-285. 10.1016/j.ccr.2006.03.014.
Drost J, Mantovani F, Tocco F, Elkon R, Comel A, Holstege H, Kerkhoven R, Jonkers J, Voorhoeve PM, Agami R, Del Sal G: BRD7 is a candidate tumour suppressor gene required for p53 function. Nat Cell Biol. 2010, 12: 380-389. 10.1038/ncb2038.
Brummelkamp TR, Nijman SM, Dirac AM, Bernards R: Loss of the cylindromatosis tumour suppressor inhibits apoptosis by activating NF-kappaB. Nature. 2003, 424: 797-801. 10.1038/nature01811.
Helmbold H, Komm N, Deppert W, Bohn W: Rb2/p130 is the dominating pocket protein in the p53-p21 DNA damage response pathway leading to senescence. Oncogene. 2009, 28: 3456-3467. 10.1038/onc.2009.222.
Durkin SG, Glover TW: Chromosome fragile sites. Annu Rev Genet. 2007, 41: 169-192. 10.1146/annurev.genet.41.042007.165900.
Bergamaschi A, Kim YH, Wang P, Sorlie T, Hernandez-Boussard T, Lonning PE, Tibshirani R, Borresen-Dale AL, Pollack JR: Distinct patterns of DNA copy number alteration are associated with different clinicopathological features and gene-expression subtypes of breast cancer. Genes Chromosomes Cancer. 2006, 45: 1033-1040. 10.1002/gcc.20366.
Adelaide J, Finetti P, Bekhouche I, Repellini L, Geneix J, Sircoulomb F, Charafe-Jauffret E, Cervera N, Desplans J, Parzy D, et al: Integrated profiling of basal and luminal breast cancers. Cancer Res. 2007, 67: 11565-11575. 10.1158/0008-5472.CAN-07-2536.
Foulkes WD: BRCA1 functions as a breast stem cell regulator. J Med Genet. 2004, 41: 1-5. 10.1136/jmg.2003.013805.
Clarke RB, Howell A, Potten CS, Anderson E: Dissociation between steroid receptor expression and cell proliferation in the human breast. Cancer Res. 1997, 57: 4987-4991.
Young SR, Pilarski RT, Donenberg T, Shapiro C, Hammond LS, Miller J, Brooks KA, Cohen S, Tenenholz B, Desai D, et al: The prevalence of BRCA1 mutations among young women with triple-negative breast cancer. BMC Cancer. 2009, 9: 86-10.1186/1471-2407-9-86.
Bertheau P, Turpin E, Rickman DS, Espie M, de Reynies A, Feugeas JP, Plassa LF, Soliman H, Varna M, de Roquancourt A, et al: Exquisite sensitivity of TP53 mutant and basal breast cancers to a dose-dense epirubicin-cyclophosphamide regimen. PLoS Med. 2007, 4: e90-10.1371/journal.pmed.0040090.
Huret JL, Dessen P, Bernheim A: Atlas of Genetics and Cytogenetics in Oncology and Haematology, year 2003. Nucleic Acids Res. 2003, 31: 272-274. 10.1093/nar/gkg126.
Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes. Nat Rev Cancer. 2004, 4: 177-183. 10.1038/nrc1299.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/10/654/prepub
This work was supported by grants from the Dutch Cancer Society (NKI 2002-2635), the Netherlands Organization for Scientific Research (NWO Vidi 917.036.347) and the Susan G. Komen Breast Cancer Foundation (BCTR0403230). The authors thank Dr. Sabine Linn for sharing her expert view on this subject and Dr. Piet Borst and Dr. Hein te Riele for critically reading this manuscript.
The authors declare that they have no competing interests.
HH (first author) designed, performed the analyses and wrote the manuscript, HH (second author) and MJV provided tumor information of the Horlings tumor panel . AV provided help with statistical analyses and performed the bioinformatics. AL and ALB-D sequenced TP53 from the tumor panel described by Horlings et al  and helped with TP53 mutation interpretation, PN provided tumor information from the Joosse tumor panel . JJ supervised the project. All authors read and approved the manuscript.