|Cross-platform pathway-based analysis identifies markers of response to the PARP inhibitor olaparib.|
|Jump to Full Text|
|PMID: 22875744 Owner: NLM Status: MEDLINE|
|Poly(ADP-ribose) polymerase (PARP) is an enzyme involved in DNA repair. PARP inhibitors can act as chemosensitizers, or operate on the principle of synthetic lethality when used as single agent. Clinical trials have shown drugs in this class to be promising for BRCA mutation carriers. We postulated that inability to demonstrate response in non-BRCA carriers in which BRCA is inactivated by other mechanisms or with deficiency in homologous recombination for DNA repair is due to lack of molecular markers that define a responding subpopulation. We identified candidate markers for this purpose for olaparib (AstraZeneca) by measuring inhibitory effects of nine concentrations of olaparib in 22 breast cancer cell lines and identifying features in transcriptional and genome copy number profiles that were significantly correlated with response. We emphasized in this discovery process genes involved in DNA repair. We found that the cell lines that were sensitive to olaparib had a significant lower copy number of BRCA1 compared to the resistant cell lines (p value 0.012). In addition, we discovered seven genes from DNA repair pathways whose transcriptional levels were associated with response. These included five genes (BRCA1, MRE11A, NBS1, TDG, and XPA) whose transcript levels were associated with resistance and two genes (CHEK2 and MK2) whose transcript levels were associated with sensitivity. We developed an algorithm to predict response using the seven-gene transcription levels and applied it to 1,846 invasive breast cancer samples from 8 U133A/plus 2 (Affymetrix) data sets and found that 8-21 % of patients would be predicted to be responsive to olaparib. A similar response frequency was predicted in 536 samples analyzed on an Agilent platform. Importantly, tumors predicted to respond were enriched in basal subtype tumors. Our studies support clinical evaluation of the utility of our seven-gene signature as a predictor of response to olaparib.|
|Anneleen Daemen; Denise M Wolf; James E Korkola; Obi L Griffith; Jessica R Frankum; Rachel Brough; Lakshmi R Jakkula; Nicholas J Wang; Rachael Natrajan; Jorge S Reis-Filho; Christopher J Lord; Alan Ashworth; Paul T Spellman; Joe W Gray; Laura J van't Veer|
Related Documents :
|16625354 - Global transcriptional analysis of methanosarcina mazei strain gö1 under different nit...
22466714 - Comparative mapping, genomic structure, and expression analysis of eight pseudo-respons...
22369384 - Effective gene collection from the metatranscriptome of marine microorganisms.
22612724 - Comparative analysis of p2-type remnant prophage loci in xenorhabdus bovienii and xenor...
19262614 - Indoleamine 2,3-dioxygenase attenuates inhibitor development in gene-therapy-treated he...
3014524 - Enhanced expression of the bacterial chloramphenicol acetyltransferase gene in mouse ce...
|Type: Journal Article; Research Support, N.I.H., Extramural; Research Support, Non-U.S. Gov't; Research Support, U.S. Gov't, Non-P.H.S. Date: 2012-08-09|
|Title: Breast cancer research and treatment Volume: 135 ISSN: 1573-7217 ISO Abbreviation: Breast Cancer Res. Treat. Publication Date: 2012 Sep|
|Created Date: 2012-08-28 Completed Date: 2013-01-04 Revised Date: 2013-07-12|
Medline Journal Info:
|Nlm Unique ID: 8111104 Medline TA: Breast Cancer Res Treat Country: Netherlands|
|Languages: eng Pagination: 505-17 Citation Subset: IM|
|Laboratory Medicine, University of California San Francisco, 2340 Sutter Street Box 0808, San Francisco, CA 94115, USA. email@example.com|
|APA/MLA Format Download EndNote Download BibTex|
Breast Neoplasms / drug therapy*
Cell Line, Tumor / drug effects
Cell Survival / drug effects
DNA Repair / genetics
Drug Evaluation, Preclinical
Drug Resistance, Neoplasm
Inhibitory Concentration 50
Neoplasms, Basal Cell / drug therapy*
Phthalazines / pharmacology*
Piperazines / pharmacology*
Poly(ADP-ribose) Polymerases / antagonists & inhibitors*
Tumor Markers, Biological / genetics*, metabolism
|P50 CA 58207/CA/NCI NIH HHS; U54 CA112970/CA/NCI NIH HHS|
|0/AZD 2281; 0/Antineoplastic Agents; 0/Phthalazines; 0/Piperazines; 0/Tumor Markers, Biological; EC 126.96.36.199/PARP1 protein, human; EC 188.8.131.52/Poly(ADP-ribose) Polymerases|
Journal ID (nlm-ta): Breast Cancer Res Treat
Journal ID (iso-abbrev): Breast Cancer Res. Treat
Publisher: Springer US, Boston
© The Author(s) 2012
Received Day: 19 Month: 6 Year: 2012
Accepted Day: 25 Month: 7 Year: 2012
Electronic publication date: Day: 9 Month: 8 Year: 2012
pmc-release publication date: Day: 9 Month: 8 Year: 2012
Print publication date: Month: 9 Year: 2012
Volume: 135 Issue: 2
First Page: 505 Last Page: 517
PubMed Id: 22875744
Publisher Id: 2188
|Cross-platform pathway-based analysis identifies markers of response to the PARP inhibitor olaparib|
|Denise M. Wolf1|
|James E. Korkola2|
|Obi L. Griffith2|
|Jessica R. Frankum3|
|Lakshmi R. Jakkula2|
|Nicholas J. Wang2|
|Jorge S. Reis-Filho3|
|Christopher J. Lord3|
|Paul T. Spellman2|
|Joe W. Gray4|
|Laura J. van’t Veer1||
1Laboratory Medicine, University of California San Francisco, 2340 Sutter Street Box 0808, San Francisco, CA 94115 USA
2Cancer & DNA Damage Responses, Lawrence Berkeley National Laboratories, One Cyclotron Road, Berkeley, CA 94720 USA
3The Breakthrough Breast Cancer Research Centre, The Institute of Cancer Research, 237 Fulham Road, London, SW3 6JB UK
4Department of Biomedical Engineering, Oregon Health and Science University, 3303 SW Bond Avenue, Room #13000, Portland, OR 97239 USA
Several mechanisms have been identified in mammalian cells that function to maintain genome integrity [1, 2]. Mechanisms for the repair of single strand breaks (SSB) include base excision repair (BER), mismatch repair (MMR), and nucleotide excision repair (NER). DNA damage due to double strand breaks (DSB) is repaired via the homologous recombination (HR) pathway or with non-homologous end joining (NHEJ). The Fanconi anemia (FA)/BRCA pathway is involved in repair of DSBs or DNA interstrand cross links, and the DNA damage response (DDR) pathway is a network of DDRs for the regulation of many steps in the DNA repair process .
These repair pathways are frequently deregulated in cancer cells, motivating efforts to develop drugs that are preferentially effective in cells with defective repair. One emerging class of therapeutic agents for this purpose is based on the principle of synthetic lethality. Synthetic lethality is defined as cell death that results from complete inactivation of two genes in cells where inactivation of either alone does not result in death [4, 5]. The first therapy to reach the clinic was based on the concept that synthetic lethality would result from coordinate inactivation of HR repair (due to a genomic aberration) and poly(ADP-ribose) polymerase (PARP). This approach was based on studies showing that cells with HR repair deficiency caused by BRCA1/2 mutation [6, 7], PTEN mutation [8, 9], RAD51D loss of function , and PALB2 loss of function  exhibited synthetic lethality with an inhibitor of one or more PARP family proteins, differentiating cancerous from non-cancerous cells by only targeting cells with complete loss of the former genes. The PARP family consists of 18 PARP domain enzymes , with the most well-studied members being PARP1, PARP2, VPARP (PARP4), tankyrase 1 (TNKS, PARP5a), and TNKS2 (PARP5b). PARP1 is the most ubiquitous member. This protein rapidly binds to sites of damaged DNA to modulate a variety of proteins involved in DNA repair and other cellular processes, making it a key protein in the BER pathway for SSB repair.
PARP is involved in SSB DNA repair and PARP inhibitors cause some of them to be converted into DSBs at replication forks [13, 14]. In HR competent cells, DSBs are repaired so that the cells can survive. However, in HR-deficient cells, DSBs are repaired via the less accurate NHEJ pathway or the single strand annealing subpathway of HR, resulting in chromatid aberrations that usually lead to cell death. These conditions therefore make cells with BRCA mutations or other HR defects [11, 15, 16] preferentially sensitive to (i.e., to show synthetic lethality with) PARP inhibitors. PARP inhibitors also have been proposed as possibly useful for treatment of triple negative breast cancers that exhibit “BRCAness” [7, 17]. BRCAness is defined as the spectrum of phenotypes that some sporadic tumors share with BRCA mutated cancers, reflecting the underlying distinctive DNA repair defects arising from loss of HR; for example, by epigenomic down regulation of BRCA1 and FANCF .
PARP inhibitors in clinical studies for BRCA-associated, triple negative and/or basal-like breast cancer include olaparib (AstraZeneca, London), ABT-888 (also known as Veliparib; Abbott Laboratories, IL), and PF-01367338 (AG014699; Pfizer Inc., NY) [14, 18, 19]. These agents are licensed for monotherapy in DNA repair deficient patients or as chemo-potentiating agents after SSBs are created by common anticancer treatments such as radiotherapy and DNA damaging agents . Results obtained from clinical trials so far, however, appear to vary depending on the specific breast cancer patient population, the specificity of the PARP inhibitor, and the nature of the therapeutic agent used in combination with the PARP inhibitor [20, 21]. A multicenter phase 2 trial showed that olaparib as monotherapy led to objective response rates in 41 % of BRCA1/2 mutation carriers who had previously received several courses of chemotherapy . Results for triple negative breast cancer patients without known BRCA1/2 mutations have been inconsistent. Preclinical studies and phase 1 trials suggested that PARP inhibitors can increase cell death in these patients when combined with paclitaxel , whilst triple negative breast cancer patients largely did not respond to olaparib monotherapy in a phase 2 trial . Thus, our aim in this study was to identify candidate biomarkers that can be tested for their ability to better identify subsets of sporadic cancers with defects in HR-directed repair that will respond to PARP inhibitors.
We focused in this study on olaparib, a small-molecule, reversible, oral inhibitor of both PARP1 and PARP2 . We identified candidate biomarkers associated with response to olaparib by correlating responses to nine concentrations of olaparib in a panel of well-characterized breast cancer cell lines with the transcription levels of genes involved in aspects of DNA repair. Genes tested for correlation with olaparib response included those reported in the literature to be directly relevant to PARP inhibitor response or involved more generally in some aspect of DNA repair (Fig. 1). We applied this signature to primary tumor data to identify the frequency and characteristics of tumors that might be expected to respond to olaparib. These studies set the stage for a clinical test of the sensitivity and specificity of this predictor and indicate known subtypes of breast cancers that might be preferentially sensitive to olaparib.
The sensitivity of a panel of 22 breast cancer cell lines to KU0058948 (olaparib; KuDOS Pharmaceuticals/AstraZeneca) was measured with a growth inhibition assay as described in Supplementary Material and [7, 26]. The following molecular data were collected for the panel: copy number (Affymetrix SNP6), gene expression (Affymetrix U133A, Affymetrix Exon 1.0 ST), transcriptome sequencing (Illumina GAII), methylation (Illumina Methylation27), protein abundance (reverse protein lysate array), and mutation status (COSMIC, ). A detailed description of the availability and preprocessing of all molecular data sets is provided in Supplementary Material and .
The Wilcoxon rank sum test was used to test the association of drug response with individual biomarkers. Drug response was associated with subtype, triple negativity, and mutation status with use of the Fisher’s exact test. Due to the small sample size, a p < 0.05 was deemed significant, whilst a p < 0.1 was considered a trend. Logistic regression (LR) with forward feature selection (fivefold CV) was used to identify candidate biomarkers and was applied to each considered DNA repair pathway separately. The resulting biomarkers were combined into a predictor using a weighted voting algorithm . The Matlab code used for signature development and validation is provided in Supplementary Material. A Chi-square test was used to test for associations of breast cancer subtype with response to olaparib. We refer to Supplementary Material for a detailed description of the statistical methods.
Twenty-two breast cancer cell lines previously profiled for RNA transcript levels were tested for response to nine concentrations of olaparib (see Table 1). These cells mirror many of the transcriptional and genomic characteristics of primary breast tumors and have been used to model responses to a large number of experimental and approved therapeutic compounds [28, 30]. The concentration of olaparib needed to reduce survival to 50 % (SF50) was used as a quantitative measure of sensitivity and ranged from 0.44 nM to 32 µM. The SF50 was not reached for five cell lines at the maximum treatment concentration of 50 µM olaparib. Olaparib response obtained with the growth inhibition assay was not influenced by growth rate assessed as doubling time (Spearman correlation coefficient −0.036, p value 0.874). Figure 2 shows the waterfall plot of SF50 with cell lines ordered from most resistant at the left to most sensitive at the right. Cell lines were divided into a group of 15 resistant and 7 sensitive cell lines, based on an SF50 threshold of 1 µM. Drug response was not significantly associated with breast cancer subtype (p value luminal vs. basal 0.136; Fig. 3), and did not differ between ERBB2-amplified and non-ERBB2-amplified cell lines (p value 1), with transcriptional subtypes assigned to cell lines as previously reported . Four of the seven sensitive cell lines (57 %) were triple negative, compared to 5 of 15 (33 %) resistant cell lines (p value 0.376). Table 1 summarizes characteristics for the 22 cell lines, with SF50, doubling time, transcriptional ER, PR, and ERBB2 status, and the molecular data available for each of them.
We selected candidate molecular features that might be developed as biomarkers for prediction of response to olaparib as those features involved in DNA repair activities that were associated with quantitative response to olaparib in the cell line panel. Molecular features included pretreatment RNA transcript levels, mutation status, copy number variation, and promoter methylation status. Specific genes tested involved aspects of DNA repair listed by Wang and Weaver ; ER, PR and ERBB2 due to the importance of PARP inhibition for triple negative breast cancer [17, 19]; and PARP family members PARP1, PARP2, VPARP, TNKS, and TNKS2. This approach is based on observations that in vitro models showing high sensitivity to PARP inhibitors often have BRCA and PTEN deficiencies [7, 8], copy number variations involving BRCA1 and PARP1 , and/or hypermethylation of the promoter regions of genes BRCA1 and FANCF . Molecular features showing statistically significant associations with SF50 values are summarized in Supplementary Table 1 and illustrated in Fig. 4.
The transcription levels of MRE11A, NBS1, TNKS, TNKS2, XPA, and XRCC5 were significantly lower (p < 0.05; fold-change >2) in the sensitive compared to the resistant cell lines for at least one expression platform (U133A, exon array and RNA-seq), whilst transcription levels for BRCA1, ERCC4, FANCD2, and PR tended to be lower in sensitive lines (p < 0.1). We refer to Supplementary Table 1a for the list of significant associations per platform. PR protein levels measured using reverse phase protein lysate arrays  were also significantly reduced in the sensitive cell lines (p < 0.05). Transcript levels for CHEK2 and MK2 were significantly higher in the sensitive compared to the resistant lines (p < 0.05), with a similar trend for PARP2 and XRCC3 (p < 0.1). Although PARP1 has been shown to be overexpressed in 58 % of invasive breast cancer samples  and upregulated at protein level in 82 % of BRCA1-associated breast cancer samples , there is no consensus on its importance as a biomarker of response to PARP inhibitors [36, 37]. In our cell line panel, expression of PARP1 mRNA levels were not significantly higher in the sensitive lines compared to the resistant lines (median p value 0.277) (Supplementary Table 1a).
The BRCA1-mutated cell lines MDAMB436 and SUM149PT had a trend to be more sensitive to olaparib compared to the wild-type cell lines (p value 0.091) (Supplementary Table 1b). Likewise, cells with reduced BRCA1 copy number were significantly more sensitive to olaparib than cells with normal copy number at this locus (p value 0.012) (Supplementary Table 1c). PTEN loss of function, which was defined as mutation and/or lack of expression, was not significantly associated with olaparib SF50 response (p value 0.145), even though previous studies from our group suggested that PTEN deficiency can cause olaparib sensitivity [8, 9]. Lack of association in the cell line panel could be ascribed to the small sample size and/or to the possibility that the univariate associations do not take into account important multigene effects. Since BRCA1 mutations have been associated with reduced PTEN expression , we tested for association of either BRCA1 mutation or PTEN deficiency with olaparib sensitivity. We found that cell lines with a deficiency in either gene tended to be more sensitive to olaparib than cell lines with functional BRCA1 and PTEN (p value 0.052) (Supplementary Table 1b). No association was found between TP53 mutation status and drug response (p value 0.376).
We used a breast cancer cell line panel comprised luminal, basal, and claudin-low cell lines to develop a multi-transcript predictor of sensitivity to olaparib according to the REMARK recommendations . We limited the predictor to transcript levels to facilitate clinical application. We considered all breast cancer subtypes for the development of the predictor based on a study of RAD51 focus formation in cells responding to a PARP inhibitor. That study showed that 30–40 % of triple negative breast cancers appeared not to have defective HR and therefore might not benefit from a PARP inhibitor whilst ~20 % of non-triple negative breast cancers appeared to have defective HR and therefore might respond to a PARP inhibitor . Thus, we reasoned that a predictor developed using the complete cell line panel might be applicable to the full spectrum of breast cancer covered by the cell line panel. As shown in Fig. 1, the molecular features tested as candidate biomarkers were limited to genes involved in DNA repair pathways BER, NER, MMR, HR/FA, NHEJ, and DDR as defined by Wang and Weaver  and in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database release 55.1 . This led to the selection of 118 genes (see Supplementary Table 2) that were tested for association between transcript levels and response to olaparib. These transcript levels were measured using three different mRNA analysis platforms (Affymetrix U133A arrays, Affymetrix exon arrays and Illumina RNA-seq).
We identified the most important transcripts by applying LR with forward feature selection (fivefold CV) 100 times. Markers significantly associated with olaparib response in over half of the iterations are shown in Table 2. These were further reduced to seven gene transcripts that were significantly associated with olaparib response in at least 2 out of 3 mRNA analysis platforms. Five transcript levels (candidate resistance markers BRCA1, MRE11A, NBS1, TDG, and XPA) were inversely associated with predicted probability of response and two transcript levels (candidate sensitivity markers CHEK2 and MK2) were positively associated with predicted probability of response (see Table 3). BRCA1 is involved in DSB repair via RAD51-mediated HR [42, 43]. CHEK2 is a kinase with signal transduction function in cell-cycle regulation and checkpoint responses , and is involved in the major parallel DDR pathway ATM-CHEK2 . CHEK2 has also been reported as an intermediate-level breast cancer risk gene, regardless of family history [44, 45]. Besides the standard DDR pathways, the cell-cycle checkpoint pathway p38MAPK/MK2 is additionally activated in TP53 mutant cells . MK2 activity is critical for prolonged checkpoint maintenance through a process of posttranscriptional regulation of gene expression . MRE11A and NBS1 are part of the MRN complex, a multifaceted molecular machine for DSB recognition . Finally, TDG is part of the BER pathway, whilst XPA encodes a zinc finger protein that is part of the NER complex.
We combined information on the seven-transcript levels to form a predictive signature using a weighted voting algorithm (Supplementary Material and ). This algorithm assigns a weight and decision boundary to each of the seven genes, based on their expression distribution for the class of sensitive versus resistant cell lines (see Table 3). For this signature to work on external samples, the transcript levels were normalized to the geometric mean of seven control genes, followed by median normalization across the cell lines (see Supplementary Material). The larger the weight for a gene transcript level, the more influence this gene has on predicted probability of response. Positive weights were assigned for sensitivity markers and negative weights were assigned for resistance markers.
We analyzed expression profiles measured for breast cancer patients not treated with PARP inhibitors to understand which patients would have a likelihood of response to olaparib according to our seven-transcript predictor. We used seven U133A and one U133 plus two data sets on 1,846 primary breast tumors with or without metastasis, heterogeneous in treatment and ER/PR/LN status. Our seven-transcript response algorithm predicted that 8–21 % of patients in the eight data sets would be responsive to olaparib (Table 4), using threshold 0.0372 obtained from the cell lines to distinguish sensitive from resistant (see Supplementary Material). The fraction predicted to respond was inversely related to the fraction of ER-positive patients in each data set (Pearson correlation coefficient −0.614, p value 0.1). We also tested the seven-transcript predictor in Agilent mRNA transcript profiles measured for 536 breast invasive carcinoma samples collected by The Cancer Genome Atlas (TCGA) . This required that an Agilent-specific threshold distinguishing sensitive from resistant be established. We accomplished this using a set of Affymetrix and Agilent mRNA transcript profiles measured for 80 I-SPY 1 samples [50, 51]. The Agilent threshold was set so that the fraction of I-SPY 1 samples in the Agilent data set predicted to be sensitive was the same as that predicted to be sensitive using the Affymetrix data (see Supplementary Materials). The fraction of samples predicted to be sensitive in the TCGA data set was 12 % (Table 4). We assessed the transcriptional subtypes of the patient populations predicted to respond to olaparib in 464 samples from GSE25066 and in 528 TCGA tumor samples after exclusion of the normal-like samples. The tumors predicted to respond were enriched in samples classified as basal-like compared to samples classified as luminal A, luminal B or HER2 (p value 0.002 and 2.6 × 10−28 for GSE25066 and TCGA, respectively; Table 5).
In this hypothesis generating study, our overall aim was to use quantitative measurements of response to olaparib in 22 breast cancer cell lines to identify molecular features associated with response as a first step towards development of a molecular signature to predict clinical responses. We limited our search for features associated with olaparib response to copy number, DNA sequence abnormalities or transcription levels for 42 genes suggested in  for their association with DNA repair. Molecular features associated with 15 of these 42 genes were found to be significantly associated or to show a trend of association with olaparib response. Specifically, cell lines that were sensitive to olaparib were enriched in BRCA1 mutations or deletions, PARP1 amplification, reduced expression of BRCA1, ERCC4, FANCD2, MRE11A, NBS1, PR, TNKS, TNKS2, XPA, and XRCC5 and increased expression of CHEK2, MK2, PARP2, and XRCC3.
Since multiple mechanisms may contribute to olaparib sensitivity, we developed a weighted voting signature to combine influences from multiple markers. We included only transcript levels in our algorithm since most molecular features associated with response were apparent at the transcript level. We limited the search space to molecular features of 118 genes from six principal DNA repair pathways in order to increase statistical power. Associations of transcript levels for 118 genes and responses to olaparib for 22 breast cancer cell lines resulted in a seven-gene predictive signature that included five resistance markers (BRCA1, MRE11A, NBS1, TDG, and XPA) and two response markers (CHEK2 and MK2).
The transcript levels of the seven genes in the predictor were consistent with expectations from the literature. Mutations in BRCA1, loss of heterozygosity at the BRCA1 locus and deregulated expression have been described in the literature as potential markers for prediction of response to PARP inhibitors . These studies are consistent with our finding that reduced BRCA1 transcript levels are associated with olaparib sensitivity. PARP1 is required for rapid accumulation of MRE11A at DSB sites. Due to the direct interaction between PARP1 and MRE11A, deficiency in MRE11A has been suggested as a mechanism of sensitizing cells to PARP1 inhibition based on the concept of synthetic lethality . Moreover, a dominant negative mutation in MRE11A in MMR deficient cancers has been shown to sensitize cells to agents causing replication fork stress . These reports are consistent with our finding that reduced MRE11A transcription is associated with olaparib sensitivity. Experimental disruption of the HR pathway protein NBS1 by RNAi has been reported to increase sensitivity to PARP inhibitors . This is consistent with our finding that reduced transcription of NBS1 is associated with olaparib sensitivity. Cells with defective NER have been shown to be hypersensitive to platinum agents, with low XPA protein levels in testis tumor cell lines explaining the low capacity to repair cisplatin-induced DNA damage . PARP inhibitors also enhance lethality in XPA-deficient cells after UV irradiation . Tumor cells with deficiency of the DDR pathway have been suggested to be hypersensitive to PARP inhibitors, with the DNA repair biomarker CHEK1 shown to be overexpressed in BRCA1-like versus non-BRCA1-like triple negative breast cancer . This is consistent with our finding that increased CHEK2 transcription is associated with olaparib sensitivity.
Our seven-gene transcript algorithm suggests that 8–21 % of patients with primary breast cancers may respond to olaparib and that the responsive tumors are enriched in basal-like breast cancers. This represents a hypothesis that can now be tested in clinical trials. Since the signature has not yet been tested clinically, it is inappropriate to use it to select patients for treatment with olaparib or other PARP inhibitors or to use it in any way to manage breast cancer treatment. However, it does present a signature that can be tested in planned translational analyses of ongoing clinical trials of PARP inhibitors and that can be used to determine whether clinical trials are properly sized to detect a response of the magnitude predicted by this signature.
AD is partly supported by a BAEF Fellowship of the Belgian American Educational Foundation for postdoctoral research, OLG is supported by a Fellowship from the Canadian Institutes of Health Research. This study was supported by the Director, Office of Science, Office of Biological & Environmental Research, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231; by the National Institutes of Health, National Cancer Institute Grants P50 CA 58207 Breast SPORE (LJV), U54 CA 112970, SU2C-AACR-DT0409, and the Pfizer Corporation Grant to JWG; Stand Up to Cancer/AACR Breakthrough Breast Cancer, Cancer Research UK and the Breast Cancer Research Foundation. The content of the information does not necessarily reflect the position or the policy of the Government, and no official endorsement should be inferred.
AA and CJL may benefit financially from the development of PARP inhibitors through patents held jointly with AstraZeneca through the Institute of Cancer Research “rewards to inventors” scheme. All other authors declare that they have no conflict of interest.
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
|1..||Rich T,Allen RL,Wyllie AH. Defying death after DNA damageNatureYear: 2000407680577778310.1038/3503771711048728|
|2..||Sancar A,Lindsey-Boltz LA,Unsal-Kacmaz K,Linn S. Molecular mechanisms of mammalian DNA repair and the DNA damage checkpointsAnnu Rev BiochemYear: 200473398510.1146/annurev.biochem.73.011303.07372315189136|
|3..||Ciccia A,Elledge SJ. The DNA damage response: making it safe to play with knivesMol CellYear: 201040217920410.1016/j.molcel.2010.09.01920965415|
|4..||Iglehart JD,Silver DP. Synthetic lethality—a new direction in cancer-drug developmentN Engl J MedYear: 2009361218919110.1056/NEJMe090304419553640|
|5..||Ashworth A,Lord CJ,Reis-Filho JS. Genetic interactions in cancer progression and treatmentCellYear: 20111451303810.1016/j.cell.2011.03.02021458666|
|6..||Bryant HE,Schultz N,Thomas HD,Parker KM,Flower D,Lopez E,Kyle S,Meuth M,Curtin NJ,Helleday T. Specific killing of BRCA2-deficient tumours with inhibitors of poly(ADP-ribose) polymeraseNatureYear: 2005434703591391710.1038/nature0344315829966|
|7..||Farmer H,McCabe N,Lord CJ,Tutt AN,Johnson DA,Richardson TB,Santarosa M,Dillon KJ,Hickson I,Knights C,Martin NM,Jackson SP,Smith GC,Ashworth A. Targeting the DNA repair defect in BRCA mutant cells as a therapeutic strategyNatureYear: 2005434703591792110.1038/nature0344515829967|
|8..||Mendes-Pereira AM,Martin SA,Brough R,McCarthy A,Taylor JR,Kim JS,Waldman T,Lord CJ,Ashworth A. Synthetic lethal targeting of PTEN mutant cells with PARP inhibitorsEMBO Mol MedYear: 200916–731532210.1002/emmm.20090004120049735|
|9..||Dedes KJ, Wetterskog D, Mendes-Pereira AM, Natrajan R, Lambros MB, Geyer FC, Vatcheva R, Savage K, Mackay A, Lord CJ, Ashworth A, Reis-Filho JS (2010) PTEN deficiency in endometrioid endometrial adenocarcinomas predicts sensitivity to PARP inhibitors. Sci Transl Med 2(53):53ra75. doi:10.1126/scitranslmed.3001538|
|10..||Loveday C,Turnbull C,Ramsay E,Hughes D,Ruark E,Frankum JR,Bowden G,Kalmyrzaev B,Warren-Perry M,Snape K,Adlard JW,Barwell J,Berg J,Brady AF,Brewer C,Brice G,Chapman C,Cook J,Davidson R,Donaldson A,Douglas F,Greenhalgh L,Henderson A,Izatt L,Kumar A,Lalloo F,Miedzybrodzka Z,Morrison PJ,Paterson J,Porteous M,Rogers MT,Shanley S,Walker L,Eccles D,Evans DG,Renwick A,Seal S,Lord CJ,Ashworth A,Reis-Filho JS,Antoniou AC,Rahman N. Germline mutations in RAD51D confer susceptibility to ovarian cancerNat GenetYear: 201143987988210.1038/ng.89321822267|
|11..||Buisson R,Dion-Cote AM,Coulombe Y,Launay H,Cai H,Stasiak AZ,Stasiak A,Xia B,Masson JY. Cooperation of breast cancer proteins PALB2 and piccolo BRCA2 in stimulating homologous recombinationNat Struct Mol BiolYear: 201017101247125410.1038/nsmb.191520871615|
|12..||Rouleau M,Patel A,Hendzel MJ,Kaufmann SH,Poirier GG. PARP inhibition: PARP1 and beyondNat Rev CancerYear: 201010429330110.1038/nrc281220200537|
|13..||Caldecott KW. Mammalian single-strand break repair: mechanisms and links with chromatinDNA RepairYear: 20076444345310.1016/j.dnarep.2006.10.00617118715|
|14..||Liang H,Tan A. PARP inhibitorsCurr Breast Cancer RepYear: 20113445410.1007/s12609-010-0036-y|
|15..||McCabe N,Turner NC,Lord CJ,Kluzek K,Bialkowska A,Swift S,Giavara S,O’Connor MJ,Tutt AN,Zdzienicka MZ,Smith GC,Ashworth A. Deficiency in the repair of DNA damage by homologous recombination and sensitivity to poly(ADP-ribose) polymerase inhibitionCancer ResYear: 200666168109811510.1158/0008-5472.CAN-06-014016912188|
|16..||Williamson CT,Muzik H,Turhan AG,Zamo A,O’Connor MJ,Bebb DG,Lees-Miller SP. ATM deficiency sensitizes mantle cell lymphoma cells to poly(ADP-ribose) polymerase-1 inhibitorsMol Cancer TherYear: 20109234735710.1158/1535-7163.MCT-09-087220124459|
|17..||Turner N,Tutt A,Ashworth A. Hallmarks of ‘BRCAness’ in sporadic cancersNat Rev CancerYear: 200441081481910.1038/nrc145715510162|
|18..||Vinayak S,Ford J. PARP inhibitors for the treatment and prevention of breast cancerCurr Breast Cancer RepYear: 2010219019710.1007/s12609-010-0026-022655123|
|19..||Plummer R. Poly(ADP-ribose) polymerase inhibition: a new direction for BRCA and triple-negative breast cancer?Breast Cancer ResYear: 201113421810.1186/bcr287721884642|
|20..||Turner NC,Ashworth A. Biomarkers of PARP inhibitor sensitivityBreast Cancer Res TreatYear: 2011127128328610.1007/s10549-011-1375-821301956|
|21..||Guha M. PARP inhibitors stumble in breast cancerNat BiotechnolYear: 201129537337410.1038/nbt0511-37321552220|
|22..||Tutt A,Robson M,Garber JE,Domchek SM,Audeh MW,Weitzel JN,Friedlander M,Arun B,Loman N,Schmutzler RK,Wardley A,Mitchell G,Earl H,Wickens M,Carmichael J. Oral poly(ADP-ribose) polymerase inhibitor olaparib in patients with BRCA1 or BRCA2 mutations and advanced breast cancer: a proof-of-concept trialLancetYear: 2010376973723524410.1016/S0140-6736(10)60892-620609467|
|23..||Dent R, Lindeman G, Clemons M, Wildiers H, Chan A, McCarthy N, Singer C, Lowe E, Kemsley K, Carmichael J (2010) Safety and efficacy of the oral PARP inhibitor olaparib (AZD2281) in combination with paclitaxel for the 1st or 2nd line treatment of patients with metastatic triple negative breast cancer: results from the safety cohort of a Phase 1/2 multicentre trial. Proc Am Soc Clin Oncol 28(suppl):abstr 1018|
|24..||Gelmon KA,Tischkowitz M,Mackay H,Swenerton K,Robidoux A,Tonkin K,Hirte H,Huntsman D,Clemons M,Gilks B,Yerushalmi R,Macpherson E,Carmichael J,Oza A. Olaparib in patients with recurrent high-grade serous or poorly differentiated ovarian carcinoma or triple-negative breast cancer: a phase 2, multicentre, open-label, non-randomised studyLancet OncolYear: 201112985286110.1016/S1470-2045(11)70214-521862407|
|25..||Fong PC,Boss DS,Yap TA,Tutt A,Wu P,Mergui-Roelvink M,Mortimer P,Swaisland H,Lau A,O’Connor MJ,Ashworth A,Carmichael J,Kaye SB,Schellens JH,Bono JS. Inhibition of poly(ADP-ribose) polymerase in tumors from BRCA mutation carriersN Engl J MedYear: 2009361212313410.1056/NEJMoa090021219553641|
|26..||Edwards SL,Brough R,Lord CJ,Natrajan R,Vatcheva R,Levine DA,Boyd J,Reis-Filho JS,Ashworth A. Resistance to therapy caused by intragenic deletion in BRCA2NatureYear: 200845171821111111510.1038/nature0654818264088|
|27..||Weigelt B,Warne PH,Downward J. PIK3CA mutation, but not PTEN loss of function, determines the sensitivity of breast cancer cells to mTOR inhibitory drugsOncogeneYear: 201130293222323310.1038/onc.2011.4221358673|
|28..||Heiser LM,Sadanandam A,Kuo WL,Benz SC,Goldstein TC,Ng S,Gibb WJ,Wang NJ,Ziyad S,Tong F,Bayani N,Hu Z,Billig JI,Dueregger A,Lewis S,Jakkula L,Korkola JE,Durinck S,Pepin F,Guan Y,Purdom E,Neuvial P,Bengtsson H,Wood KW,Smith PG,Vassilev LT,Hennessy BT,Greshock J,Bachman KE,Hardwicke MA,Park JW,Marton LJ,Wolf DM,Collisson EA,Neve RM,Mills GB,Speed TP,Feiler HS,Wooster RF,Haussler D,Stuart JM,Gray JW,Spellman PT. Subtype and pathway specific responses to anticancer compounds in breast cancerProc Natl Acad Sci USAYear: 201210982724272910.1073/pnas.101885410822003129|
|29..||Moulder S,Yan K,Huang F,Hess KR,Liedtke C,Lin F,Hatzis C,Hortobagyi GN,Symmans WF,Pusztai L. Development of candidate genomic markers to select breast cancer patients for dasatinib therapyMol Cancer TherYear: 2010951120112710.1158/1535-7163.MCT-09-111720423993|
|30..||Neve RM,Chin K,Fridlyand J,Yeh J,Baehner FL,Fevr T,Clark L,Bayani N,Coppe JP,Tong F,Speed T,Spellman PT,DeVries S,Lapuk A,Wang NJ,Kuo WL,Stilwell JL,Pinkel D,Albertson DG,Waldman FM,McCormick F,Dickson RB,Johnson MD,Lippman M,Ethier S,Gazdar A,Gray JW. A collection of breast cancer cell lines for the study of functionally distinct cancer subtypesCancer CellYear: 200610651552710.1016/j.ccr.2006.10.00817157791|
|31..||Wang X,Weaver D. The ups and downs of DNA repair biomarkers for PARP inhibitor therapiesAm J Cancer ResYear: 20111330132721968427|
|32..||Holstege H,Horlings HM,Velds A,Langerod A,Borresen-Dale AL,Vijver MJ,Nederlof PM,Jonkers J. BRCA1-mutated and basal-like breast cancers have similar aCGH profiles and a high incidence of protein truncating TP53 mutationsBMC CancerYear: 20101065410.1186/1471-2407-10-65421118481|
|33..||Tibes R,Qiu Y,Lu Y,Hennessy B,Andreeff M,Mills GB,Kornblau SM. Reverse phase protein array: validation of a novel proteomic technology and utility for analysis of primary leukemia specimens and hematopoietic stem cellsMol Cancer TherYear: 20065102512252110.1158/1535-7163.MCT-06-033417041095|
|34..||Goncalves A,Finetti P,Sabatier R,Gilabert M,Adelaide J,Borg JP,Chaffanet M,Viens P,Birnbaum D,Bertucci F. Poly(ADP-ribose) polymerase-1 mRNA expression in human breast cancer: a meta-analysisBreast Cancer Res TreatYear: 2011127127328110.1007/s10549-010-1199-y21069454|
|35..||Domagala P,Huzarski T,Lubinski J,Gugala K,Domagala W. Immunophenotypic predictive profiling of BRCA1-associated breast cancerVirchows ArchYear: 20114581556410.1007/s00428-010-0988-320941507|
|36..||Cotter M, Pierce A, McGowan P, Madden S, Flanagan L, Quinn C, Evoy D, Crown J, McDermott E, Duffy M (2011) PARP1 in triple-negative breast cancer: expression and therapeutic potential. J Clin Oncol 29 (15_suppl):1061|
|37..||Zaremba T,Ketzer P,Cole M,Coulthard S,Plummer ER,Curtin NJ. Poly(ADP-ribose) polymerase-1 polymorphisms, expression and activity in selected human tumour cell linesBr J CancerYear: 2009101225626210.1038/sj.bjc.660516619568233|
|38..||Saal LH,Gruvberger-Saal SK,Persson C,Lovgren K,Jumppanen M,Staaf J,Jonsson G,Pires MM,Maurer M,Holm K,Koujak S,Subramaniyam S,Vallon-Christersson J,Olsson H,Su T,Memeo L,Ludwig T,Ethier SP,Krogh M,Szabolcs M,Murty VV,Isola J,Hibshoosh H,Parsons R,Borg A. Recurrent gross mutations of the PTEN tumor suppressor gene in breast cancers with deficient DSB repairNat GenetYear: 200840110210710.1038/ng.2007.3918066063|
|39..||McShane LM,Altman DG,Sauerbrei W,Taube SE,Gion M,Clark GM. REporting recommendations for tumor MARKer prognostic studies (REMARK)Breast Cancer Res TreatYear: 2006100222923510.1007/s10549-006-9242-816932852|
|40..||Graeser M,McCarthy A,Lord CJ,Savage K,Hills M,Salter J,Orr N,Parton M,Smith IE,Reis-Filho JS,Dowsett M,Ashworth A,Turner NC. A marker of homologous recombination predicts pathologic complete response to neoadjuvant chemotherapy in primary breast cancerClin Cancer ResYear: 201016246159616810.1158/1078-0432.CCR-10-102720802015|
|41..||Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M (2010) KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res 38(Database issue):D355–360. doi:10.1093/nar/gkp896|
|42..||Gudmundsdottir K,Ashworth A. The roles of BRCA1 and BRCA2 and associated proteins in the maintenance of genomic stabilityOncogeneYear: 200625435864587410.1038/sj.onc.120987416998501|
|43..||Tutt A,Ashworth A. The relationship between the roles of BRCA genes in DNA repair and cancer predispositionTrends Mol MedYear: 200281257157610.1016/S1471-4914(02)02434-612470990|
|44..||CHEK2 Breast Cancer Case-Control ConsortiumCHEK2*1100delC and susceptibility to breast cancer: a collaborative analysis involving 10,860 breast cancer cases and 9,065 controls from 10 studiesAm J Hum GenetYear: 20047461175118210.1086/42125115122511|
|45..||Fletcher O,Johnson N,Santos Silva I,Kilpivaara O,Aittomaki K,Blomqvist C,Nevanlinna H,Wasielewski M,Meijers-Heijerboer H,Broeks A,Schmidt MK,Van’t Veer LJ,Bremer M,Dork T,Chekmariova EV,Sokolenko AP,Imyanitov EN,Hamann U,Rashid MU,Brauch H,Justenhoven C,Ashworth A,Peto J. Family history, genetic testing, and clinical risk prediction: pooled analysis of CHEK2 1100delC in 1,828 bilateral breast cancers and 7,030 controlsCancer Epidemiol Biomarkers PrevYear: 200918123023410.1158/1055-9965.EPI-08-041619124502|
|46..||Reinhardt HC,Aslanian AS,Lees JA,Yaffe MB. p53-deficient cells rely on ATM- and ATR-mediated checkpoint signaling through the p38MAPK/MK2 pathway for survival after DNA damageCancer CellYear: 200711217518910.1016/j.ccr.2006.11.02417292828|
|47..||Reinhardt HC,Hasskamp P,Schmedding I,Morandell S,Vugt MA,Wang X,Linding R,Ong SE,Weaver D,Carr SA,Yaffe MB. DNA damage activates a spatially distinct late cytoplasmic cell-cycle checkpoint network controlled by MK2-mediated RNA stabilizationMol CellYear: 2010401344910.1016/j.molcel.2010.09.01820932473|
|48..||Williams GJ,Lees-Miller SP,Tainer JA. Mre11-Rad50-Nbs1 conformations and the control of sensing, signaling, and effector responses at DNA double-strand breaksDNA RepairYear: 20109121299130610.1016/j.dnarep.2010.10.00121035407|
|49..||The Cancer Genome Atlas Data Portal. http://tcga-data.nci.nih.gov/tcga/tcgaHome2.jsp. Accessed 13 Jan 2012|
|50..||Esserman LJ, Berry DA, Cheang MC, Yau C, Perou CM, Carey L, Demichele A, Gray JW, Conway-Dorsey K, Lenburg ME, Buxton MB, Davis SE, Van’t Veer LJ, Hudis C, Chin K, Wolf D, Krontiras H, Montgomery L, Tripathy D, Lehman C, Liu MC, Olopade OI, Rugo HS, Carpenter JT, Livasy C, Dressler L, Chhieng D, Singh B, Mies C, Rabban J, Chen YY, Giri D, Au A, Hylton N (2011) Chemotherapy response and recurrence-free survival in neoadjuvant breast cancer depends on biomarker profiles: results from the I-SPY 1 TRIAL (CALGB 150007/150012; ACRIN 6657). Breast Cancer Res Treat. doi:10.1007/s10549-011-1895-2. Accessed 13 Jan 2012|
|51..||Hatzis C,Pusztai L,Valero V,Booser DJ,Esserman L,Lluch A,Vidaurre T,Holmes F,Souchon E,Wang H,Martin M,Cotrina J,Gomez H,Hubbard R,Chacon JI,Ferrer-Lozano J,Dyer R,Buxton M,Gong Y,Wu Y,Ibrahim N,Andreopoulou E,Ueno NT,Hunt K,Yang W,Nazario A,DeMichele A,O’Shaughnessy J,Hortobagyi GN,Symmans WF. A genomic predictor of response and survival following taxane-anthracycline chemotherapy for invasive breast cancerJ Am Med AssocYear: 2011305181873188110.1001/jama.2011.593|
|52..||Vilar E,Bartnik CM,Stenzel SL,Raskin L,Ahn J,Moreno V,Mukherjee B,Iniesta MD,Morgan MA,Rennert G,Gruber SB. MRE11 deficiency increases sensitivity to poly(ADP-ribose) polymerase inhibition in microsatellite unstable colorectal cancersCancer ResYear: 20117172632264210.1158/0008-5472.CAN-10-112021300766|
|53..||Wen Q,Scorah J,Phear G,Rodgers G,Rodgers S,Meuth M. A mutant allele of MRE11 found in mismatch repair-deficient tumor cells suppresses the cellular response to DNA replication fork stress in a dominant negative mannerMol Biol CellYear: 20081941693170510.1091/mbc.E07-09-097518256278|
|54..||Koberle B,Masters JR,Hartley JA,Wood RD. Defective repair of cisplatin-induced DNA damage caused by reduced XPA protein in testicular germ cell tumoursCurr BiolYear: 19999527327610.1016/S0960-9822(99)80118-310074455|
|55..||Okano S,Kanno S,Nakajima S,Yasui A. Cellular responses and repair of single-strand breaks introduced by UV damage endonuclease in mammalian cellsJ Biol ChemYear: 200027542326353264110.1074/jbc.M00408520010924509|
|56..||Rodriguez AA,Makris A,Wu MF,Rimawi M,Froehlich A,Dave B,Hilsenbeck SG,Chamness GC,Lewis MT,Dobrolecki LE,Jain D,Sahoo S,Osborne CK,Chang JC. DNA repair signature is associated with anthracycline response in triple negative breast cancer patientsBreast Cancer Res TreatYear: 2010123118919610.1007/s10549-010-0983-z20582464|
[Figure ID: Fig1]
Approach for the development of a predictor of olaparib response in a breast cancer cell line panel with inclusion of prior knowledge of DNA repair pathways. For 22 breast cancer cell lines, growth inhibition assays were used to measure their sensitivity to olaparib, expressed as the surviving fraction at 50 % (SF50). For these cell lines, expression data were obtained with three different platforms (Affymetrix U133A, Affymetrix Exon 1.0 ST, and whole transcriptome shotgun sequencing). The bottom-up approach was used for biomarker selection, incorporating prior knowledge of the principal DNA repair pathways BER, NER, MMR, HR/FA, NHEJ, and DDR. Biomarkers from  were systematically expanded with genes assigned to any of these pathways in the KEGG database, resulting in 118 genes. For each DNA repair pathway and expression data set the most important markers were obtained with LR in combination with forward feature selection, followed by reduction to those selected with consistent pattern of sensitivity for all three platforms
[Figure ID: Fig2]
Waterfall plot of the response to olaparib (expressed as SF50 in µM) for 22 breast cancer cell lines, ordered from most resistant at the left to most sensitive at the right, with bars colored according to subtype (luminal in light grey, basal in black, claudin-low in dark grey, and ERBB2 amplified in white). The threshold of 1 µM used to divide the cell lines into a group of 15 resistant cell lines (indicated with R) and a group of 7 sensitive cell lines (indicated with S) is represented with a horizontal dashed line
[Figure ID: Fig3]
Boxplot of SF50 for the cell lines divided according to breast cancer subtype (9 luminal, 7 claudin-low, and 6 basal lines). No association was found between breast cancer subtype and response to olaparib in the cell line panel (Fisher’s exact test for basal vs. luminal, p value 0.136)
[Figure ID: Fig4]
Overview of individual DNA repair-associated markers that are significantly associated with or do trend towards an association with response to olaparib in the 22 breast cancer cell lines, based on mutation, copy number, and expression data (see Supplementary Table 1 for the complete list of markers). The four boxplots at the top show the association results for BRCA1. The BRCA1-mutated cell lines MDAMB436 and SUM149PT tend to be more sensitive to olaparib compared to the wild-type cell lines (p value 0.091). The sensitive cell lines are also characterized by a significant lower copy number of BRCA1 (p value 0.012) and by BRCA1 down-regulation (RNA-seq, p value 0.055). Cell lines with a deficiency in BRCA1 and/or PTEN tend to be more sensitive to olaparib than cell lines with functional BRCA1 and PTEN (p value 0.052). The boxplots at the bottom show the association for genes NBS1 and XRCC5 that are significantly down-regulated and for genes CHEK2 and MK2 that are significantly up-regulated in the sensitive compared to the resistant cell lines
Overview of the breast cancer cell line panel with response to olaparib expressed as SF50 (µM); ER, PR, and ERBB2 expression with + indicating up-regulation relative to the other cell lines, − down-regulation, and NC no change in expression; and availability of the different molecular data sets indicated with N for unavailability and Y for availability
|Cell line||Olaparib SF50 (µM)||Doubling time (h)||ERa||PRa||ERBB2a||COSMIC||SNP6||RPPA||Methylation||RNA-seq||Exon array||U133A|
|# cell lines||20||21||21||22||20||22||22|
Doubling times were estimated for each cell line from measurements of the number of doublings of untreated cells that occurred in 72 h during the course of assessing responses to 123 therapeutic compounds (Heiser et al. , PNAS 2012)
aFor ER, probe 205225_at on the Affymetrix U133A array was investigated; for PR, probe 208305_at; and for ERBB2 probes 210930_s_at and 216836_s_at
Overview per expression platform of genes from six principal DNA repair pathways that are selected with the LR approach in over half of the iterations
|Biomarker source||Platform||# Genes||Genes selected in >250/500 iterationsa||Avg. AUC (std)b|
|DNA repair biomarkers (Wang and Weaver )||U133A (standard)||11/29||BRCA1, BRCA2, CHEK2, DSS1, MRE11A, NBS1, PALB2, PARP2, PTEN, TP53, XPA||0.793 (0.083)|
|U133A (custom)||7/29||BRCA1, BRCA2, CHEK2, DSS1, NBS1, RAD51, XPA||0.945 (0.059)|
|Exon array||12/29||BRCA2, CHEK2, DSS1, ERCC1, ERCC4, FANCD2, MK2, MRE11A, NBS1, USP11, XPA, XRCC5||0.717 (0.084)|
|RNA-seq||14/29||ATM, BRCA1, DSS1, FANCD2, JTB, MK2, MRE11A, NBS1, PALB2, PARP1, PARP2, XPA, XRCC5, XRCC6||0.715 (0.132)|
|KEGG||U133A (standard)||5/103||DNTT, MUTYH, POLM, RPA2, TOP3B||0.745 (0.075)|
|U133A (custom)||9/103||DNTT, FEN1, MUTYH, NBS1, POLD1, POLM, RAD51, RAD51C, XRCC5||0.725 (0.092)|
|Exon array||4/103||DNTT, MRE11A, TDG, UNG||0.753 (0.083)|
|RNA-seq||5/103||DCLRE1C, FEN1, RPA4, TDG, XRCC5||0.839 (0.054)|
aGenes with consistent pattern of sensitivity for all three platforms (U133A, exon array, RNA-seq) and for both measures of class comparison (mean, median) are shown in bold
bAverage fivefold CV area under the receiver operating characteristics curve (AUC) (standard deviation) across 100 randomizations for a LR model with optimized coefficients and inclusion of the platform-specific genes selected in >1/2 of the iterations
Overview of the seven genes selected for prediction of response to treatment with olaparib based on breast cancer cell line expression data
|Gene symbol||Gene name||Pathway||Entrez gene ID||Marker||Probe||Weight (wg)||Decision boundary (bg)|
|BRCA1||Breast cancer 1, early onset||HR||672||Resistance||204531_s_at||−0.5320||−0.0153|
|CHEK2||CHK2 checkpoint homolog||DDR||11200||Sensitivity||210416_s_at||0.5806||−0.0060|
|MK2||Mitogen-activated protein kinase-activated protein kinase 2||DDR||9261||Sensitivity||201461_s_at||0.0713||0.0031|
|MRE11A||MRE11 meiotic recombination 11 homolog A||DDR/HR||4361||Resistance||205395_s_at||−0.1396||−0.0044|
|XPA||Xeroderma pigmen-tosum, complemen-tation group A||NER||7507||Resistance||205672_at||−0.2335||−0.0126|
The weights and decision boundaries were determined with data from the U133A expression array platform measured for the 22 cell lines used to assess response to olaparib. For each of the five resistance and two sensitivity markers, gene symbol is shown together with gene name, the DNA repair pathway the gene belongs to, entrez gene identifier, corresponding probe set from the Affymetrix U133A array, and weight and decision boundary obtained with the weighted voting algorithm
Prevalence of the seven-gene signature in tumor samples from 9 different studies on primary breast tumors with or without metastasis, heterogeneous in treatment and ER/PR/LN status
|Data set||Platform||# Samples||Characteristics||Treatment||Event rate (%)||# Predicted responders (%)a|
|GSE2034||U133A||286||73.1 % ER+||Untreated||37.4 Distant metastasis||55 (19.2)|
|58 % PR+|
|18.2 % ERBB2+|
|0 % LN+|
|GSE20271||U133A||177||55.7 % ER+||49.2 % FAC;
50.8 % T/FAC
|14.1 pCR||26 (14.7)|
|46.9 % PR+|
|14.2 % ERBB2+|
|GSE23988||U133A||61||52.5 % ER+||FEC/wTx||32.8 pCR||9 (14.8)|
|0 % ERBB2+|
|65.6 % LN+|
|Median tumor size 6 cm (2–17.5)|
|GSE4922||U133A + B||289||86.1 % ER+||37.7 % Systematic adjuvant therapy||35.7 Local/distant recurrence or death||24 (8.3)|
|33.7 % LN+|
|Median tumor size 2 cm (0.2–13)|
|GSE25066||U133A||508||58.9 % ER+||Neoadj. taxane & anthra-cycline-based regimen||19.5 pCR||94 (18.5)|
|69.1 % LN+|
|31.5 % LumA|
|15.3 % LumB|
|37.2 % Basal-like|
|7.3 % HER2-enr|
|8.7 % Normal-like|
|GSE7390||U133A||198||67.7 % ER+||Untreated||31.3 Distant metastasis||33 (16.7)|
|14.1 % ERBB2+|
|0 % LN+|
|Median tumor size 2 cm (0.6–5)|
|GSE11121||U133A||200||78 % ER+||Untreated||23 Distant metastasis||20 (10.0)|
|65 % PR+|
|12.3 % ERBB2+|
|0 % LN+|
|Median tumor size 2 cm (0.1–6.0)|
|GSE5460||U133 plus 2||127||58.3 % ER+||Untreated||–||27 (21.3)|
|23.6 % ERBB2+|
|49.6 % LN+|
|Median tumor size 2.2 cm (0.8–8.5)|
|TCGA||Agilent G4502A||536||44.0 % LumA||Heterogeneous||–||67 (12.5)|
|25.2 % LumB|
|18.5 % Basal-like|
|10.8 % HER2-enr|
|1.5 % Normal-like|
FAC neoadjuvant chemotherapy regimen with 5-fluorouracil, docorubicin, and cyclophosphamide; T/FAC neoadjuvant chemotherapy regimen with paclitaxel and 5-fluorouracil, docorubicin, and cyclophosphamide; FEC/wTx neoadjuvant chemotherapy regimen with four courses of 5-fluorouracil, docorubicin, and cyclophosphamide, followed by four additional courses of weekly docetaxel and capecitabine
aNumber and percentage of patients predicted to respond to treatment with a PARP inhibitor according to the seven-gene predictor with use of threshold 0.0372 for response assignment for Affymetrix data, and threshold 0.174 for Agilent data
Association of breast cancer subtype with predicted response to olaparib in 464 GSE25066 and 528 TCGA tumor samples, after exclusion of the normal-like samples
|Luminal A||120 (75.0)||40 (25.0)||Luminal A||233 (98.7)||3 (1.3)|
|Luminal B||72 (92.3)||6 (7.7)||Luminal B||126 (93.3)||9 (6.7)|
|Basal-like||155 (82.0)||34 (18.0)||Basal-like||54 (54.5)||45 (45.5)|
|HER2-enriched||35 (94.6)||2 (5.4)||HER2-enriched||50 (86.2)||8 (13.8)|
|2.6 × 10−28|
Keywords: Keywords PARP, Olaparib, Breast cancer, Biomarker, Bioinformatics.
Previous Document: Primary empty sella (PES): a review of 175 cases.
Next Document: Lapatinib, trastuzumab or the combination added to preoperative chemotherapy for breast cancer: a me...