- Research article
- Open Access
Next generation sequencing of triple negative breast cancer to find predictors for chemotherapy response
Breast Cancer Research volume 17, Article number: 134 (2015)
In triple negative breast cancers (TNBC) the initial response to chemotherapy is often favorable, but relapse and chemotherapy resistance frequently occur in advanced disease. Hence there is an urgent need for targeted treatments in this breast cancer subtype. In the current study we deep sequenced DNA of tumors prior to chemotherapy to search for predictors of response or resistance.
Next generation sequencing (NGS) was performed for 1,977 genes involved in tumorigenesis. DNA from 56 pre-treatment TNBC-biopsies was sequenced, as well as matched normal DNA. Following their tumor biopsy, patients started neoadjuvant chemotherapy with doxorubicin and cyclophosphamide. We studied associations between genetic alterations and three clinical variables: chemotherapy response, relapse-free survival and BRCA proficiency.
The mutations observed were diverse and few recurrent mutations were detected. Most mutations were in TP53, TTN, and PIK3CA (55 %, 14 %, and 9 %, respectively). The mutation rates were similar between responders and non-responders (average mutation rate 9 vs 8 mutations). No recurrent mutations were associated with chemotherapy response or relapse. Interestingly, PIK3CA mutations were exclusively observed in patients proficient for BRCA1. Samples with a relapse had a higher copy number alteration rate, and amplifications of TTK and TP53BP2 were associated with a poor chemotherapy response.
In this homogenous cohort of TNBCs few recurrent mutations were found. However, PIK3CA mutations were associated with BRCA proficiency, which can have clinical consequences in the near future.
Neoadjuvant chemotherapy (NAC) is the standard of care for locally advanced breast cancer. For triple negative breast cancer, which is characterized by the absence of expression of estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2), it is currently the only available systemic treatment option. Approximately half of all TNBC reaches a pathological complete remission (pCR) upon NAC. However, tumors without a pCR often relapse. When distant metastases develop, prognosis is poor, with a median post recurrence survival of only one year [1, 2]. As it is impossible to predict which patients will respond and which ones will not, predictive biomarkers are urgently needed. The same is true for targeted treatments, particularly for non-responding tumors.
With the introduction of next generation sequencing (NGS), tumor collections have been sequenced to categorize mutations and to find new treatment targets [3–5]. Many new mutations and several important signaling pathways have been identified. However, one of the main findings of the sequencing efforts is the marked heterogeneity within breast cancer. This heterogeneity is highlighted by the fact that recurrent mutations are infrequent. No two tumors have the same molecular make up. Consequently, subsequent studies attempting to identify broad mechanisms of chemotherapy resistance are not likely to be successful. There are many different routes to tumorigenesis, and each tumor has its own sequence of tumor-causing molecular events. The emergence of drug resistance may be driven by additional mutations and add a further level of complexity. Clinical data are not available for most studies, precluding the integration of molecular and clinical data. In addition, most existing molecular studies have been undertaken on a broad mix of breast cancer samples, with widely varying clinical stage or molecular subtype. A more homogeneous study population with well-characterized clinical data may be required to find clinical predictive factors.
To find mutations predicting chemotherapy benefit in TNBC, we performed NGS in a clinically well-characterized cohort. We sequenced an NGS panel covering 1,977 genes involved in tumorigenesis and treatment. We sequenced pretreatment biopsies from patients with a good response and patients with a poor response, to identify chemotherapy predictive markers. The purpose of this study was twofold: 1) to find chemotherapy response-predictive biomarkers and 2) to identify new treatment targets.
Patients and tumor biopsies
Pre-treatment biopsies and peripheral blood samples were obtained from 56 patients with untreated, primary triple negative breast cancer. All patients had received neoadjuvant treatment at the Netherlands Cancer Institute between 2004 and 2013 as part of ongoing clinical trials, or were treated off protocol according to the standard arms of one of these studies  (NCT00448266, NCT01057069). The ethical committee of the Netherlands Cancer Institute approved the studies and all patients gave informed consent (reference numbers of Ethical approval: PTC06.1725/N06IAA and PTC09.2716/M09TNM). All tumors were either at least 3 cm in size, or the presence of axillary lymph node metastases had been established by fine needle aspiration or pre-treatment sentinel node biopsy. Biopsies were taken using a 14-G core needle under ultrasound guidance. After collection, specimens were snap-frozen in liquid nitrogen and stored at −80 °C. Each patient had two or three biopsies taken to ensure that enough tumor material was available for both diagnostic and research purposes.
All patients started NAC with three courses of 2-weekly administrations of cyclophosphamide and doxorubicin (ddAC). The selection of the chemotherapy regimen for the fourth to the sixth courses depended on the specific clinical trial and could consist of continuation of ddAC in the case of a favorable magnetic resonance imaging (MRI) response (for criteria see previous publication ), randomization to intensive carboplatin-based chemotherapy or a switch to a combination of capecitabine and docetaxel in the case of an unfavorable MRI response. All patients subsequently underwent surgery, either breast conserving or mastectomy. The clinical results of these strategies have not been published because the studies are ongoing (NCT00448266, NCT01057069).
Triple negative status was defined by the absence of ER and PR expression and no amplification of HER2. ER (Roche Diagnostics cat. no. 5278406001 (Roche, Basel, Switserland))- and PR (Roche Diagnostics cat. no. 5278392001)-negativity was defined as immunohistochemical (IHC) staining of fewer than 10 % of tumor cell nuclei. None of the samples in this cohort had between 1 % and 10 % positive tumor cell nuclei, so all samples fulfilled the ASCO/CAP criteria for triple negativity. Negativity for HER2 (Roche Diagnostics cat. no. 5278368001) amplification was present if IHC staining was graded as 0 or 1+. In the case of 2+ or 3+ staining, chromogenic in situ hybridization was performed to determine HER2 amplification (gene copy number ≥6 per tumor cell). Chemotherapy response was assessed by microscopic examination of the surgery resection specimen. The complete absence of any invasive tumor cells in both the breast and the lymph nodes was considered as a pCR. All other responses were assigned to the no-pCR group. An experienced breast cancer pathologist (JW) reviewed all pathology slides.
The series was well-characterized for BRCA function. Briefly, germ line DNA was isolated from peripheral blood lymphocytes of affected patients. We used mutation-scanning methods. The protein truncation test was used for exon 11 of BRCA1 and exons 10 and 11 of BRCA2. The remaining exons were tested using denaturing gradient gel electrophoresis. Confirmation of aberrant samples was done by Sanger sequencing . In addition, multiplex ligation-dependent probe amplification (MLPA) was performed using MLPA kit P087 (BRCA1) to detect large genomic deletions or duplications in the genes.
Hypermethylation of the BRCA1 promoter was determined using methylation-specific MLPA analysis, according to the manufacturer’s protocol (MLPA assay ME005-custom, MRC-Holland, Amsterdam, The Netherlands). For normalization and analysis the Coffalyzer program was used (MRC-Holland). According to the manufacturer’s protocol, we used a cutoff of 20 % to call a sample methylated . Employing this cutoff, methylated samples have very low BRCA1 gene expression [10, 11], which was also true for this series (Additional file 1: Figure S1).
DNA was isolated from frozen biopsies with a minimal tumor percentage of 50 %. Tumor cell percentage was assessed by an experienced breast pathologist (JW). Isolation was performed with (Qiagen, Venlo, the Netherlands) DNA mini kit. Matched normal DNA was obtained from peripheral blood and extracted by DNAzol and purified with Qiagen DNeasy kit. Samples were interrogated by a designed “Cancer mini-genome” consisting of 1,977 cancer genes . This set comprises genes involved in DNA repair, cell cycle, apoptosis, epigenetic modification (methylation, acetylation), angiogenesis (vascular endothelial growth factor (VEGF) pathway), and genes from the epidermal growth factor receptor (EGFR), mammalian target of rapamycin (mTOR), insulin, TP53, transforming growth factor (TGF)-beta, Notch, Wnt and hedgehog pathways, known tumor suppressor and proto-oncogenes, and all additional kinases not present in any of the above mentioned processes. Genes are listed in Additional file 2: Table S1. Barcoded fragment libraries were generated from 300–600 ng of isolated DNA from tumor and reference samples as previously described . Pools of libraries were enriched for 1,977 cancer-related genes using (SureSelect technology, Agilent, Santa Clara, California, US). Enriched libraries were sequenced to achieve an average coverage of at least 150× on a (SOLiD 5500xl, Applied Biosystems, Waltham, Massachusetts, US) instrument according to the manufacturers’ protocol. Sequencing statistics are available in Additional file 3: Table S2.
Sequence data analysis
Sequence reads were mapped on the human reference genome version 19 (GRCh37), using BWA (−c –l 25 –k 2 –n 10)  and variant calling was done using a custom pipeline identifying variants with at least 10× coverage, a 20 % allele frequency, and multiple (> = 2) occurrences in the seed (the first 25-bp most accurately mapped part of the read) and support from independent reads (> = 3). Validation of our custom variant-calling pipeline is described in Nijman et al. . Putative somatic variants were identified by subsequently genotyping all variant positions in the raw datasets of both the tumor and reference sample using samtools mpileup, to ensure the absence of variant alleles in the reference sample; only positions showing less than 5 % variant alleles in the reference sample were considered somatic. Next, somatic variants with a minor allele frequency >0.01 in non-Catalogue of Somatic Mutations in Cancer (COSMIC) populations in dbSNP were filtered out, as were variants seen > = 1× homozygously or > = 3× heterozygously in reference samples in our in-house database and variants with an allele frequency <25 % if the coverage was <50×. Finally, manual curation was performed to remove obvious noise in poorly mapped regions. For the analysis of association with clinical data and the pathway analysis, we included only putative driver mutations defined as either 1) mutations with a major effect on protein expression, defined as a) truncating mutations, indels, splice acceptors/donors, initiator codon variants or b) non-synonymous mutations predicted to be damaging according to PolyPhen  and/or SIFT , or 2) known COSMIC mutations.
Copy number analysis was done by comparing depth of coverage in reference and tumor samples using the robust z score according to the method of Iglewicz and Hoaglin . Briefly, the median and median absolute deviation of coverage per exon of all reference samples, normalized on the total number of reads per sample, was calculated to determine the coverage distribution. The z -scores (i.e., the difference between the coverage of the tumor sample and the median coverage of the reference pool of normal samples, multiplied by a correction factor of 0.6745, and divided by the median absolute deviation of the reference pool) were calculated for all exons of each tumor sample. The z score of the total gene was calculated by taking the average of the z scores per exon. As we noticed an effect of GC-content on the exon-coverage profiles in some, but not all samples, we used unsupervised hierarchical clustering of the normalized coverage per exon of all samples. Two clear clusters were formed, but 12 samples (including both tumor and reference samples) could not be placed in either cluster. Exon coverage of the tumor samples was then compared to the reference samples within the same cluster. Genes with z scores > = 10 were considered high-level amplifications; a z score > = 3 indicates copy number gain and < = −3 copy number loss.
Sanger sequencing for validation of PIK3CA mutations
Samples were analyzed by Sanger Sequencing using a BigDye Terminator Cycle Sequencing Kit and ABI 3730 capillary sequencer (Life Technologies/Applied Biosystems, Waltham, Massachusetts, US).
Statistics: association of sequencing data with clinical information
Statistical analyses were performed with R/Bioconductor . We focused on three clinical variables: 1) pCR, 2) relapse and 3) BRCA proficiency. We encoded each of the clinical variables into a binary variable and then studied associations at the gene and at the pathway levels: association between mutations in an individual gene and the given clinical variable was tested using Fisher’s exact test of the presence/absence of a mutation in a sample and the presence/absence of the clinical variable. The pathway analysis was based on a list of 29 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG), including at least three genes targeted in the sequencing experiments. We considered a pathway to be mutated if at least one of its genes was found mutated. We then performed the Fisher’s exact test of the presence/absence of a mutation in a given pathway and the presence/absence of the clinical variable in the various samples. All tests were corrected for multiple testing using the Benjamini-Hochberg method. In all analyses, associations were considered significant when the adjusted p value was <0.05.
Description of the dataset
We performed deep sequencing on a set of fifty-six triple negative pre-treatment biopsies, scheduled to receive NAC. The pathological complete remission rate was 45 % (Table 1). With a median follow up of 2.5 years, nine patients had a relapse (16 %). Twelve patients had a BRCA1 or BRCA2 germline mutation, and seventeen biopsies showed BRCA1 promoter methylation. In total, BRCA function (based on either a BRCA1/2 mutation or BRCA1 promoter methylation) was impaired in at least 29 tumors (52 %). This number could be higher as BRCA1 germline screening was not performed for five cases, as these cases did not qualify for routine diagnostic BRCA1 and BRCA2 germline testing (based on age, family history). Between 1 and 36 somatic mutations per sample were detected (10 mutations on average; Additional file 4: Table S3).
Few recurrent mutations in TNBC
The somatic mutation rates were similar for responders and non-responders (average mutation rate of nine versus eight mutations per tumor, p = 0.32). The mutation spectrum was broad with few recurrent mutations in the dataset. Figure 1 shows the mutations occurring in three or more samples. Mutations in TP53, TTN and PIK3CA were the most frequently observed alterations in this study (55 %, 14 %, and 9 %, respectively). Potentially dubious genes were extracted from Lawrence et al.  and indicated as such in the mutation data (Additional file 4: Table S3).
Association of mutations with clinical variables
There were no putative driver mutations significantly associated with either chemotherapy response (pCR rate) or with relapse (Additional file 5: Table S4). An explanation could be the heterogeneous nature of TNBCs and lack of recurrent mutations in TNBCs.
Next, we assessed whether any putative driver mutations were associated with BRCA deficiency, either resulting from a germline mutation or from promoter methylation. PIK3CA mutations were exclusively observed in patients with functional BRCA genes (5/27 (19 %)) versus 0/29 in the patients with a BRCA mutation or BRCA1 promoter methylation (p = 0.02, adjusted p value (p-adj) = 0.11) (Table 2). Tumors with BRCA impairment may develop via different routes than tumors with functioning BRCA genes, and the PI3K pathway may play a role in the latter. This finding has also been recently observed in an independent dataset . All PIK3CA mutations were confirmed with Sanger sequencing.
We next investigated whether mutations in different genes would cluster in some known pathways. Using KEGG, we derived a pathway mutation status based on its gene members (see “Methods”) and looked for association with the clinical variables (Table 3; Additional file 6: Table S5 and Additional file 7: Table S6). Calcium signaling was inversely associated with a pCR (p = 0.03), as was Notch signaling (p = 0.03). BRCA mutated or methylated tumors were associated with mutations in the Wnt signaling pathway (p = 0.06). However, after correcting for multiple testing, none of the associations remained statistically significant. No pathways were associated with relapse.
Copy number analysis
Fifty out of fifty-six tumors had NGS data suitable for copy number analysis. We identified high-level, focal amplifications of several growth factor receptors such as EGFR and FGFR2, and amplification of AIM1 in three samples, MET amplification in three samples and CCNE1 and PIK3CA in single cases (Additional file 8: Figure S2). Multiple focal, high-level amplifications were visible in several samples, affecting up to five different chromosomes per patient, while other samples showed no such events. We also detected recurrent gains of cyclins (CCND1, CCND2), MYC and FGFR1. Chromosome 1q and 8q (including MYC) showed gains in several samples, both events were frequently detected in the same sample.
Patients with a relapse or without a pCR in the breast had a higher number of genes with copy number alterations (CNA) than patients without a relapse or with a pCR in the breast (p = 0.0091 and p = 0.023, respectively). Specific amplifications and deletions were associated with response, relapse and BRCA proficiency. Most interesting is that TTK and TP53BP2 amplifications were associated with a poor response (p = 0.02 and p = 0.06, respectively). However, after correction for multiple testing, the p values were no longer significant. Table 4 gives an overview of the CNAs related to the clinical phenotypes. TTK is involved in mitotic checkpoint, and over-expression was recently reported to be associated with TNBC. TTK might be an attractive therapeutic target . TTK and TP53BP2 amplifications were also observed in Nimblegen copy number data of the same samples (Additional file 9: Figure S3).
In this study we deep-sequenced a cohort of 56 pre-treatment breast cancer biopsies to find chemotherapy predictive biomarkers and new treatment targets. We found that TNBC was very heterogeneous with few recurrent mutations. As a result, we were not able to find biomarkers associated with chemotherapy response and relapse. We found that tumors without BRCA impairment (through a mutation or through BRCA1-promoter methylation) frequently had PIK3CA mutations.
The heterogeneous nature of TNBC has already been described by others [3, 5, 23] and makes it difficult to find predictive biomarkers. There may well be hundreds of different mutations associated with resistance to standard chemotherapy. This fact would explain why overall comparison of groups of sensitive versus resistant tumors does not lead to candidate markers. The same problem exists for gene expression studies: due to several factors, among which the highly variable makeup of tumors may be central, it has proven to be extremely difficult to find predictive profiles . Different genes may explain resistance in subgroups of tumors as we have shown for HER2-negative tumors in our gene expression studies [25, 26]. A possibility to analyze NGS mutation data in a different manner is to investigate logical combinations of the mutation status of several genes and how such a combination relates to a clinical outcome variable of interest. For instance, mutation in gene A and not in gene B could be predictive of response. However, the number of possible combinations is so large that such an analysis is not practical with datasets of the size presented here. In this study we have used external knowledge (KEGG pathways) to focus on the combination of some specific subsets of genes. It would be interesting to investigate combinations of genes beyond these annotated pathways.
Previous research by our group and by others has shown that a proportion of TNBC shares similarities with BRCA1-mutated breast cancers [10, 27–29]. This phenotype is often called “BRCAness”. As BRCA1 and BRCA2 are involved in the repair of DNA double-strand breaks (DSBs), a process called homologous recombination, dysfunctional BRCA could render a tumor hypersensitive to drugs that (ultimately) induce DNA DSBs during the DNA replication phase, such as the PARP inhibitors or alkylating chemotherapy. In this study we observed that 22 % of the tumors had a BRCA mutation and 30 % had BRCA1 promoter methylation, indicating that BRCA inactivation is observed in over half of the tumors. PIK3CA mutations were only observed in the tumors proficient for BRCA. Our group has been extensively studying tumors with BRCA1 promoter methylation or a BRCA1-like genomic profile. We found that these BRCA1-like tumors resemble BRCA1 mutation carriers in multiple respects, i.e., age, grade and good chemotherapy response . Most importantly those BRCA1-like tumors have an excellent response to DSB-inducing chemotherapy . The association between non-BRCA1-like tumors and PIK3CA mutations was also seen in an independent study . Other studies found that PIK3CA mutations are more frequent in HER2+ and luminal breast cancer; it might be that we defined a subgroup of TNBCs where PIK3CA mutations are also frequent, i.e., the non-BRCA subgroups. Another consequence of this finding could be that some non-BRCA-related tumors could be susceptible to PIK3CA inhibitors.
As far as we know this is the first deep sequencing study comparing well-responding and poorly-responding breast tumors in the neoadjuvant setting to identify biomarkers. Balko et al.  performed NGS on TNBC after neoadjuvant treatment. This study focused on a cohort of cancers with macroscopic residual disease after NAC, which is associated with poor survival. When we compare our results with the Balko study we see the following. The Balko study detected a higher frequency of several targetable alterations when comparing data with The Cancer Genome Atlas (TCGA) Breast . These alterations included MCL1, PTEN, and JAK2. As expected, the percentage of these mutations in our study was more comparable to the Balko data than to TCGA (Additional file 10: Table S7). Whereas the Balko dataset contained only samples that had been taken after neoadjuvant chemotherapy, the TCGA cohort contains all clinical breast cancer stages. However, the Balko series only consists of tumors with a poor response to chemotherapy response, whereas our series consists of both responding and non-responding tumors. We also detected a high percentage of specific copy number alterations, like DDR2, PYGO2 and SCYL3. As the Balko study sequenced only 182 oncogenes, those genes were not in their sequencing panel.
A limitation of this study is the small number of samples with relapse. This fact limits the power and reduces the likelihood of finding associations with relapse. However, if there were strong associations with relapse they would have been picked up.
Although our dataset is small with 56 tumors, and precludes us from making any firm conclusions about markers associated with neoadjuvant chemotherapy response in TNBC, we can draw several potentially important conclusions. First, the mutation spectrum remains diverse even in a carefully selected and untreated group of patients with TNBC from the same clinical stage (locally advanced breast cancer, i.e., larger than 3 cm and/or positive lymph nodes). In addition, all samples were from the same institute and the laboratory procedures were thoroughly monitored: tumor percentage was 50 % or more, a matching normal DNA sample was available for every tumor, and pathology assessment of the whole series had been performed by a single specialized breast pathologist (JW). Our results strongly suggest that every tumor has its unique molecular makeup. Second, we found that PIK3CA mutations were specifically associated with BRCA1 proficient tumors, a finding that was also observed in another study . This finding may have clinical consequences. This study shows that it is increasingly important to have large series of homogenous sample series, with clinically well-annotated data.
Two important conclusions can be drawn from these data. First, TNBC is heterogeneous. Even in a homogenous, clinically well-annotated sample series, few recurrent mutations can be identified. Second, PIK3CA mutations were associated with BRCA proficient tumors. After validation in larger series, triple negative tumors with PIK3CA mutations can be candidates for agents targeting these oncogenic pathways. To identify recurrent mutations associated with chemotherapy or relapse, larger sample series are urgently needed.
copy number alteration
Catalogue of Somatic Mutations in Cancer
dose dense cyclophosphamide, doxorubicin
double strand break
epidermal growth factor receptor
human epidermal growth factor receptor 2
Kyoto Encyclopedia of Genes and Genomes
multiplex ligation-dependent probe amplification
mammalian target of rapamycin
next generation sequencing
pathological complete remission
transforming growth factor
triple negative breast cancer
vascular endothelial growth factor
Guarneri V, Broglio K, Kau SW, Cristofanilli M, Buzdar AU, Valero V, et al. Prognostic value of pathologic complete response after primary chemotherapy in relation to hormone receptor status and other factors. J Clin Oncol. 2006;24:1037–44.
Liedtke C, Mazouni C, Hess KR, Andre F, Tordai A, Mejia JA, et al. Response to neoadjuvant therapy and long-term survival in patients with triple-negative breast cancer. J Clin Oncol. 2008;26:1275–81.
The Cancer Genome Atlas Project. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490:61–70.
Stephens PJ, Tarpey PS, Davies H, Van Loo P, Greenman C, Wedge DC, et al. The landscape of cancer genes and mutational processes in breast cancer. Nature. 2012;486:400–4.
Shah SP, Roth A, Goya R, Oloumi A, Ha G, Zhao Y, et al. The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012;486:395–9.
Hannemann J, Oosterkamp HM, Bosch CA, Velds A, Wessels LF, Loo C, et al. Changes in gene expression associated with response to neoadjuvant chemotherapy in breast cancer. J Clin Oncol. 2005;23:3331–42.
Loo CE, Straver ME, Rodenhuis S, Muller SH, Wesseling J, Vrancken Peeters MJ, et al. Magnetic resonance imaging response monitoring of breast cancer during neoadjuvant chemotherapy: relevance of breast cancer subtype. J Clin Oncol. 2011;29:660–6.
van der Hout AH, van den Ouweland AM, van der Luijt RB, Gille HJ, Bodmer D, Bruggenwirth H, et al. A DGGE system for comprehensive mutation screening of BRCA1 and BRCA2: application in a Dutch cancer clinic setting. Hum Mutat. 2006;27:654–66.
MLPA company website www.mlpa.com.
Joosse SA, Brandwijk KI, Mulder L, Wesseling J, Hannemann J, Nederlof PM. Genomic signature of BRCA1 deficiency in sporadic basal-like breast tumors. Genes Chromosomes Cancer. 2011;50:71–81.
Lips EH, Mulder L, Hannemann J, Laddach N, Vrancken Peeters MT, van de Vijver MJ, et al. Indicators of homologous recombination deficiency in breast cancer and association with response to neoadjuvant chemotherapy. Ann Oncol. 2011;22:870–6.
Vermaat JS, Nijman IJ, Koudijs MJ, Gerritse FL, Scherer SJ, Mokry M, et al. Primary colorectal cancers and their subsequent hepatic metastases are genetically different: implications for selection of patients for targeted treatment. Clin Cancer Res. 2012;18:688–99.
Harakalova M, Mokry M, Hrdlickova B, Renkens I, Duran K, van Roekel H, et al. Multiplexed array-based and in-solution genomic enrichment for flexible and cost-effective targeted next-generation sequencing. Nat Protoc. 2011;6:1870–86.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
Nijman IJ, van Montfrans JM, Hoogstraat M, Boes ML, van de Corp L, Renner ED, et al. Targeted next-generation sequencing: a novel diagnostic tool for primary immunodeficiencies. J Allergy Clin Immunol. 2014;133:529–34.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.
Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4:1073–81.
Iglewicz B, Hoaglin DC. How to detect and handle outliers. Milwaukee (WI): ASQC Quality Press; 1993.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
Lawrence MS, Stojanov P, Polak P, Kryukov GV, Cibulskis K, Sivachenko A, et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature. 2013;499:214–8.
Severson TM, Peeters J, Majewski I, Michaut M, Bosma A, Schouten PC, et al. BRCA1-like signature in triple negative breast cancer: Molecular and clinical characterization reveals subgroups with therapeutic potential. Mol Oncol. 2015;9:1528-38.
Maire V, Baldeyron C, Richardson M, Tesson B, Vincent-Salomon A, Gravier E, et al. TTK/hMPS1 is an attractive therapeutic target for triple-negative breast cancer. PLoS One. 2013;8:e63712.
Wheler JJ, Parker BA, Lee JJ, Atkins JT, Janku F, Tsimberidou AM, et al. Unique molecular signatures as a hallmark of patients with metastatic breast cancer: implications for current treatment paradigms. Oncotarget. 2014;5:2349–54.
Borst P, Wessels L. Do predictive signatures really predict response to cancer chemotherapy? Cell Cycle. 2010;9:4836–40.
de Ronde JJ, Lips EH, Mulder L, Vincent AD, Wesseling J, Nieuwland M, et al. SERPINA6, BEX1, AGTR1, SLC26A3, and LAPTM4B are markers of resistance to neoadjuvant chemotherapy in HER2-negative breast cancer. Breast Cancer Res Treat. 2013;137:213–23.
de Ronde JJ, Rigaill G, Rottenberg S, Rodenhuis S, Wessels LF. Identifying subgroup markers in heterogeneous populations. Nucleic Acids Res. 2013;41:e200.
Turner N, Tutt A, Ashworth A. Hallmarks of ‘BRCAness’ in sporadic cancers. Nat Rev Cancer. 2004;4:814–9.
Silver DP, Richardson AL, Eklund AC, Wang ZC, Szallasi Z, Li Q, et al. Efficacy of neoadjuvant Cisplatin in triple-negative breast cancer. J Clin Oncol. 2010;28:1145–53.
Lips EH, Mulder L, Oonk A, van der Kolk LE, Hogervorst FB, Imholz AL, et al. Triple-negative breast cancer: BRCAness and concordance of clinical features with BRCA1-mutation carriers. Br J Cancer. 2013;108:2172–7.
Vollebergh MA, Lips EH, Nederlof PM, Wessels LF, Schmidt MK, van Beers EH, et al. An aCGH classifier derived from BRCA1-mutated breast cancer and benefit of high-dose platinum-based chemotherapy in HER2-negative breast cancer patients. Ann Oncol. 2011;22:1561–70.
Balko JM, Giltnane JM, Wang K, Schwarz LJ, Young CD, Cook RS, et al. Molecular profiling of the residual disease of triple-negative breast cancers after neoadjuvant chemotherapy identifies actionable therapeutic targets. Cancer Discov. 2014;4:232–45.
This research was supported by the Center for Translational Molecular Medicine (project Breast CARE (030–104)) and the Center for Personalized Cancer Treatment, which is a collaboration between all academic cancer centers and several teaching hospitals in the Netherlands and supported by Barcode for Life. We would like to acknowledge the NKI-AVL Core Facility Molecular Pathology & Biobanking (CFMPB) for supplying NKI-AVL Biobank material and/or laboratory support.
The authors have declared no conflicts of interest.
EHL, MM, SR and LFAW were responsible for the study design. EHL and MM coordinated the study. LM, NJMB, MJK performed the laboratory experiments, under the supervision of EC, PMN and JW. Data and statistical analysis were performed by MM and MH under the supervision of LFAW. EHL, MM, EC, EEV, RB, PMN, JW, SR and LFAW took part in data interpretation. EHL and MM drafted the manuscript with significant help of all authors. All authors read and approved the final manuscript.
Esther H. Lips and Magali Michaut contributed equally to this work.
BRCA1 gene expression, measured by RT-PCR, for BRCA1-methylated, mutated and non-mutated/non-methylated triple negative breast cancer (TNBC) samples. This figure shows the association between BRCA proficiency/deficiency and BRCA1 gene expression, as measured by RT-PCR. Samples with either BRCA1 promoter methylation or a BRCA1 mutation all had low gene expression values. (DOCX 134 kb)
1977 Cancer Gene Panel. This file lists all 1,977 genes in the Cancer Gene Panel, used in this study. (XLS 111 kb)
Overview of sequencing statistics. (XLSX 48 kb)
List of all mutations. This file shows all mutations passing thresholds as explained in “Methods”. The column “confirmed” shows mutations that have been confirmed by either Sanger Sequencing or deep whole exome sequencing in another research project (unpublished data). The column “included in analysis” shows the mutations that passed our threshold for putative driver mutations, based on effect, SIFT, polyphen scores and COSMIC (see “Methods”). The column “dubious” shows the mutations that might be dubious based on gene length, expression or replication time . (XLSX 64 kb)
Association between mutations and chemotherapy response or relapse. This table shows the top five genes associated with either response or relapse. No significant associations were found. (DOCX 18 kb)
Association between relapse and pathway mutations. This table shows the top five pathways associated with relapse. No significant associations were found. (DOCX 14 kb)
Specific mutations within affected pathways. This file shows all pathways with KEGG IDs, and mutations affecting the different pathways. (XLSX 10 kb)
Heatmap of all copy number aberrations in the dataset. Samples are represented by rows, genomic regions by the columns sorted on genomic order. Red identifies a high level gain or amplification; blue a loss. (PNG 130 kb)
Validation of TTK, TP53BP2 and genome-wide copy number profiles with Nimblegen 135 K aCGH arrays. aCGH and exome data, for the samples for which aCGH data were available (three cases for TTK, five cases for TP53BP2). Arrow indicates respectively the TTK and TP53BP2 locus. From page 3 on we show all 14 samples for which we have both exome and aCGH data, irrespective of TTK and TP53BP2 status. In these rainbow plots, the samples containing a TTK or TP53BP2 gain are indicated with an arrow at the respective genomic location. (PDF 1639 kb)
Comparison of top alterations listed by Balko et al. (reference ) (A) and our study (B) for TCGA , Balko et al.  and our study. Part A shows the top alterations as listed by Balko (reference ) in our study and TCGA (reference ), part B shows the top alterations in our study. (DOCX 17 kb)
About this article
Cite this article
Lips, E.H., Michaut, M., Hoogstraat, M. et al. Next generation sequencing of triple negative breast cancer to find predictors for chemotherapy response. Breast Cancer Res 17, 134 (2015) doi:10.1186/s13058-015-0642-8
- Next Generation Sequencing
- Triple Negative Breast Cancer
- PIK3CA Mutation
- Copy Number Alteration
- Chemotherapy Response