- Research article
- Open Access
Relationship of PIK3CA mutation and pathway activity with antiproliferative response to aromatase inhibition
Breast Cancer Research volume 16, Article number: R68 (2014)
PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit α) somatic mutations are the most common genetic alteration in breast cancer (BC). Their prognostic value and that of the phosphatidylinositol 3-kinase (PI3K) pathway in BC remains only partly defined. The effect of PIK3CA mutations and alterations of the PI3K pathway on the antiproliferative response to aromatase inhibitor treatment was determined.
The Sequenom MassARRAY System was used to determine the presence of 20 somatic mutations across the PIK3CA gene in 85 oestrogen receptor–positive (ER+) BC patients treated with 2 weeks of anastrozole before surgery. Whole-genome expression profiles were used to interrogate gene signatures (GSs) associated with the PI3K pathway. Antiproliferative activity was assessed by the change in Ki67 staining between baseline and surgery. Three GSs representing the PI3K pathway were assessed (PIK3CA-GS (Loi), PI3K-GS (Creighton) and PTEN-loss-GS (Saal)).
In our study sample, 29% of tumours presented with either a hotspot (HS, 71%) or a nonhotspot (non-HS, 29%) PIK3CA mutation. Mutations were associated with markers of good prognosis such as progesterone receptor positivity (PgR+) (P = 0.006), low grade (P = 0.028) and luminal A subtype (P = 0.039), with a trend towards significance with degree of ER positivity (P = 0.051) and low levels of Ki67 (P = 0.051). Non-HS mutations were associated with higher PgR (P = 0.014) and ER (P < 0.001) expression than both wild-type (WT) and HS-mutated samples, whereas neither biomarker differed significantly between WT and HS mutations or between HS and non-HS mutations. An inverse correlation was found between the Loi signature and both the Creighton and Saal signatures, and a positive correlation was found between the latter signatures. Lower pretreatment Ki67 levels were observed in mutation compared with WT samples (P = 0.051), which was confirmed in an independent data set. Mutation status did not predict change in Ki67 in response to 2 weeks of anastrozole treatment; there was no significant difference between HS and non-HS mutations in this regard.
PIK3CA mutations are associated with classical markers of good prognosis and signatures of PI3K pathway activity. The presence of a PIK3CA mutation does not preclude a response to neoadjuvant anastrozole treatment.
About 80% of primary breast cancer (BC) presents as oestrogen receptor–positive (ER+) disease . In postmenopausal patients, aromatase inhibition is the most effective means of treatment, reducing the relative risk of recurrence by 20% to 25% from the results achieved with tamoxifen therapy. However, resistance remains a substantial problem and is the focus for the development of new therapeutic strategies targeting putative resistance mechanisms.
To date, although researchers in in vitro studies have reported large numbers of putative mechanisms of resistance, few of these mechanisms have been confirmed as clinically relevant (for example, loss of ER expression, de novo or acquired overexpression of human epidermal growth factor receptor 2 (HER2)) . Detailed molecular characterisation of individual tumours, along with assessment of clinical efficacy, is required in order to identify additional mechanisms and determine their importance. This is particularly relevant, given that multiple new agents with targets that have putative involvement in endocrine resistance are now available .
The phosphatidylinositol 3-kinase (PI3K) pathway is frequently altered in cancer and is a focus of interest in ER + BC. This pathway links the signalling of type I and type II receptor tyrosine kinases (epidermal growth factor (EGFR), HER2 and insulin-like growth factor 1 receptor) to the effector serine/threonine kinase moiety AKT. Complete phosphorylation of AKT results in activation of the mammalian target of rapamycin (mTOR)/regulated associated protein of TOR (Raptor) complex 1 (mTORC1), which regulates cell proliferation and protein synthesis . The phosphatase and tensin homolog (PTEN) tumour suppressor gene, amongst others, is a negative regulator of the PI3K pathway and is often lost in BC , leading to enhanced activation of AKT and downstream partners.
Whilst there are multiple isoforms of the PI3K regulatory and catalytic subunits, single amino acid substitutions in BC are frequently found in the catalytic subunit p110α (PIK3CA), in which the occurrence of mutations is reported to be between 18% and 41% [4–15]. Within PIK3CA, there are three major hotspots (HSs) of mutations concentrated in the helical (E542K and E545K) and kinase (H1047R) domains, accounting for approximately 89% of PIK3CA mutations .
PIK3CA mutations have been reported to be associated with markers of good prognosis, such as high ER, high PgR, smaller size, earlier stage, lower grade [5–7, 13, 14] and better outcome in ER+/HER2- disease, as well as with good prognosis with tamoxifen monotherapy, relapse-free survival, overall survival and BC-specific survival [5, 6, 16]. In contrast, BC patients with PIK3CA mutations and HER2 amplification have been found to exhibit significantly shorter progression-free survival following trastuzumab-based treatment . However, in some cases, study researchers who have reported conflicting results have been unable to verify these associations [10, 13, 16].
Recognition of the frequency of PIK3CA mutations in clinical specimens and their oncogenic nature in vitro has led to the clinical development of numerous inhibitors of the pathway. The appropriate targeting of these agents in postmenopausal patients with ER + disease requires an improved understanding of the importance of PIK3CA mutations and PI3K pathway activation to BC response or resistance to aromatase inhibition. This insight should have a significant impact on the design and interpretation of clinical trials of endocrine therapy with PI3K pathway–targeted inhibitors.
Study of the influence of putative mechanisms and markers of resistance in early disease is not possible at an individual patient level in relation to benefit from adjuvant therapy, owing to the absence of reliable markers of subclinical disease and the inability to characterise the molecular nature of micrometastatic disease. In contrast, presurgical treatment allows the response of primary disease (unaffected by previous therapy) to be characterised by multiple means and to be correlated with the molecular characteristics. The nuclear proliferation marker Ki67 has been validated as an intermediate endpoint of evaluation of benefit from endocrine therapy in multiple studies, and it is therefore helpful for studying both resistance and the impact of targeted therapies alongside endocrine agents .
We therefore conducted a comprehensive assessment of PIK3CA mutation status and gene signatures (GSs) of PI3K pathway activity in tumours from patients with ER + primary BC in a clinical trial of the aromatase inhibitor anastrozole and related these to the drug’s antiproliferative effects and other key biomarkers of prognosis in ER + disease.
All available pretreatment (baseline) core-cut tumour samples were selected for analysis from patients who received the aromatase inhibitor anastrozole alone prior to biopsy in the multicentre 0223 neoadjuvant study (ZD1839IL/0223) of anastrozole with or without EGFR inhibitor gefitinib in ER + primary BC . The original study received approval from an institutional review board at each site and was conducted in accordance with the 1964 Declaration of Helsinki and the International Conference on Harmonization/Good Clinical Practice guidelines. Written informed consent was obtained from each patient before participation. Tumour size was measured by caliper, magnetic resonance imaging, physical palpation or ultrasound as per the preference of each centre. Previously reported scores for ER, PgR, HER2 (positivity (+) or negativity (-) by immunohistochemistry (IHC) and fluorescence in situ hybridisation) and Ki67 levels (percentage of positive cells) were available . ER and PgR expression were assessed by H-score, which incorporates both the staining intensity (scored from 0 to 3) and the percentage of cells stained. For both biomarkers, we used an H-score ≥1 as a cut-off for designation of positive or negative status in accordance with the American Society of Clinical Oncology/College of American Pathologists guidelines for ER and PgR positivity .
Percent residual Ki67 was defined as (2-week Ki67 level ÷ baseline Ki67 level) × 100. This percentage was used as an endpoint for treatment response, both as a continuous variable and dichotomised into responders or nonresponders according to a 2-week/pre-Ki67 residual <65% or > 65%, respectively . This cut-off was a prospectively defined primary endpoint of the original study and has previously been shown to predict clinical response to hormone therapy.
Validation data set
Tumours from a set of 310 pretreatment ER + patients, predominantly HER2-, from 2 neoadjuvant aromatase inhibitor studies were used for confirmation of the relationship of Ki67 data with PIK3CA mutation status .
Gene expression arrays
Whole-genome expression data for 66 pretreatment biopsies and their derivation have been described in detail previously . Intrinsic subtyping was assigned using the 50-gene predictor for luminal A and luminal B subtypes together with the gene expression data .
Multiple (one to three) 10-μm sections were taken from formalin-fixed, paraffin-embedded material, stained with nuclear fast red and needle-dissected for tumour-enriched areas using adjacent haematoxylin and eosin–stained sections as a guide, with the aim of yielding >70% malignant tissue. DNA was extracted using the DNeasy Tissue & Blood Kit (Qiagen, Hilden, Germany) with xylene and ethanol washes to remove paraffin and an overnight incubation in 1 M NaSCN. Nonsynonymous somatic PIK3CA mutations previously described in BC patients were selected from public databases, including COSMIC  and the National Cancer Institute’s somatic mutation database . Mutation detection was conducted by using the Sequenom MassARRAY System (Sequenom, San Diego, CA, USA) as described by others . In brief, with this array, DNA is amplified and treated with shrimp alkaline phosphatase to neutralise unincorporated deoxyribonucleotide triphosphates, the mutation site is extended and the sample is read by matrix-assisted laser desorption/ionisation time-of-flight (MALDI-TOF) mass spectrometry analysis. PCR and extension primers were designed using Sequenom Assay Design Software (Sequenom) and covered 20 mutations (3 HS and 17 low-frequency mutations) using 5 multiplexes (Additional file 1: Table S1). To be considered mutated, the proportion of the mutant allele had to be ≥10% as calculated by the allelotyping function of the Sequenom system.
Droplet digital PCR
HS mutations identified by mass spectrometry–based genotyping were confirmed using a droplet digital PCR (ddPCR) system (Bio-Rad Laboratories, Hercules, CA, USA). Digital PCRs were performed using this system with PIK3CA primers and probes (described in Additional file 2: Table S2) at final concentrations of 900 nM for primers and 250 nM for probes.
PCRs were prepared with 10 μl of master mix in a total volume of 20 μl and partitioned into approximately 14,000 droplets per sample in the droplet generator component of the QX100 ddPCR system according to the manufacturer’s instructions. Emulsified PCRs were run on a 96-well plate on a G-Storm GS4 thermal cycler (G-Storm, Somerton, UK). The plates were incubated at 95°C for 10 minutes, followed by 45 cycles at 95°C for 15 seconds and 63.1°C for 60 seconds (exon 9), or 45 cycles at 95°C for 15 seconds and 61.7°C for 60 seconds (exon 20), and then by a 10-minute incubation at 98°C. The temperature ramp increment was 2.5°C/s for all steps. Plates were read on a QX100 droplet reader with QuantaSoft v22.214.171.124 software (both from Bio-Rad Laboratories) to assess the number of droplets positive for PIK3CA mutant (PIK3CAmu), wild type (WT), both, or neither. At least two negative control wells with no DNA were included in every run.
The concentration of PIK3CA mutant DNA (copies of PIK3CAmu DNA per droplet) was estimated based on the Poisson distribution using the number of PIK3CAmu copies per droplet, MPIK3CAmu = -ln(1 - (nPIK3CAmu/N)), where nPIK3CAmu is the number of droplets positive for the PIK3CA 6-carboxyfluorescein (FAM) probe and N is the total number of droplets. The DNA concentration in the reaction was estimated as follows: MDNAconc = -ln(1 - (nDNAconc/N)), where nDNAconc is the number of droplets positive for PIK3CA FAM probe and PIK3CA VIC probe (Life Technologies, Carlsbad, CA, USA) and N is the total number of droplets. The fraction mutation is MPIK3CAmu/MDNAconc.
Multiple GSs have been developed to represent the PI3K pathway, such as the PIK3CA-GS , PTEN-loss-GS  and PI3K-GS . We reproduced each signature by identifying the genes that were represented on the array (Illumina, Inc, San Diego, CA, USA) and, where multiple probes existed, taking the most variable of those probes. The GS score for each sample was calculated by performing gene centring and then weighted averaging as previously described :
where n is the number of genes in a module, Xi represents the normalised gene expression in the sample and gene-specific weights Wi are equal to +1 or -1 according to the direction of their association with the phenotype in the original publication.
PIK3CA-GS (Loi), PTEN-loss-GS (Saal) and PI3K-GS (Creighton) were also adapted by the removal of proliferation-associated genes (PAGs) to generate three new signatures (PIK3CA-noPAG-GS, PTEN-loss-noPAG-GS and PI3K-noPAG-GS), because the presence of multiple PAGs dominate the GSs and could potentially mask the PI3K pathway–related features of the GSs. These PAGs were identified by three methods: genes functionally associated with cell-cycle progression and cell division according to Gene Ontology annotations , genes showing cell-cycle stage–specific expression [28, 29] and a tumour-based ‘proliferation cluster’ [28, 30–32]. The details of the genes used from each GS are listed in Additional file 3: Table S3.
All statistical analyses were performed in GraphPad Prism 6.0 (GraphPad Software, La Jolla, CA, USA) and R version 3.0.1 (R Foundation for Statistical Computing [http://www.r-project.org/]) or in MS Excel (Microsoft, Redmond, WA, USA). Associations between tumour size (T2 or T3), nodal status (positive or negative for nodal involvement), stage (2, 3 or 4), grade (1, 2 or 3), histology (ductal, lobular or other), HER2 positivity, intrinsic subtypes (luminal A vs luminal B) and PIK3CA mutation status were assessed using Fisher’s exact test. The association between continuous variables, such as ER, PgR, Ki67 and GS scores, with mutation status was evaluated using the nonparametric Mann–Whitney U test or Kruskal–Wallis test with Dunnett’s multiple comparison test. For HS vs non-HS comparisons, if a sample contained both types of mutation, it was excluded from the analysis. Correlations between the GS scores and Ki67 IHC values were calculated using Spearman’s rank correlation analysis. All P-values reported are two-sided, and values <0.05 were considered statistically significant. Cases with missing values were excluded. The number of samples within each analysis can be found in the relevant analysis and in Additional file 4: Table S4.
PIK3CAmutations in ER + breast cancer
DNA from tumours of 126 patients included in the 0223 trial was available, among which 19 samples were not assessed because of a lack of sufficient material. We processed 107 samples using the Sequenom system, which provided us with complete mutational profiles for 85 samples to which the data reported here relate. The demographics of the 85 patients are shown in Table 1.
At least one PIK3CA mutation was detected in 25 tumours (29%). Of those 25, 23 had one mutation, 1 had 3 mutations (exon 7 E453K, exon 9 E542K and E545K) and another had 2 mutations (exon 9 E545K and exon 20 H1047R). Of the 28 mutations, 20 (71%) were in HSs: 14% E542K (n = 4), 32% E545K (n = 9) and 25% H1047R (n = 7) (Figure 1, Additional file 5: Table S5).
The presence of E542K, E545K and H1047L/R mutations was also assessed by ddPCR in 25 samples with mutations, including 19 samples that contained 1 of these mutations, which were revealed by Sequenom analysis. Nineteen of the twenty-one mutations were confirmed by ddPCR. Of the two that were not confirmed, one was below the threshold for sample calling (>25 VIC + droplets) by ddPCR because of low quantity of DNA and the other was discordant. However, as a result of this excellent correlation, the mutations detected with the Sequenom system were carried forward for downstream analyses.
Relationship of PIK3CAmutations with clinicopathological features
The presence of any mutation was associated with low-grade tumours (P = 0.033), but not with size as characterized by T category (T2 vs T3; P = 0.61) (Table 1). Quantitative ER and PgR protein expression (assessed by H-score) was available for 83 of 85 cases. In patients with PIK3CA mutations, the median ER levels were significantly higher: WT = 215.1 (n = 58, IQR = 173.7 to 230.9) and Mu = 230.5 (n = 25, IQR = 214.3 to 251.1) (P = 0.017) (Figure 2A). Amongst the mutation cases, median PgR H-scores were also significantly higher: WT = 90.6 (n = 58, IQR = 2.25 to 205.3) and Mu = 213.8 (n = 25, QR = 109.4 to 275.6) (P = 0.006) (Figure 2B). All except two of the mutation cases were PgR+. Despite this, when PgR was dichotomised as positive or negative, no significant association with mutation status was identified (Table 1).
ER H-score and PgR H-score were also analysed with regard to HS and non-HS mutations. We found that non-HS mutations exhibited significantly higher ER and PgR expression compared to WT and HS samples. For ER, the median expression data were as follows: WT = 215.1 (n = 58, IQR = 173.7 to 230.9), HS = 220.4 (n = 17, IQR = 202.2 to 247.7), non-HS = 249.4 (n = 7, IQR = 235.4 to 255.2), HS vs non-HS/WT = ns, WT vs non-HS (P < 0.001). For PgR expression, the median values were as follows: WT = 90.60 (n = 58, IQR = 2.25 to 205.3), HS = 168.0 (n = 17, IQR = 77.7 to 224.4) and non-HS = 258.3 (n = 7, IQR = 196.2 to 281.9), with a nonsignificant correlation between HS and non-HS/WT and a significant difference between WT and non-HS (P = 0.014). One sample contained both a HS mutation and a non-HS mutation and was therefore excluded from the analysis.
Regarding the relationship between HER2 status and PIK3CA mutations, 69 cases were HER2- and harboured 21 mutations and 13 cases were HER2+ and contained 4 mutations (P = ns) (Table 1).
Intrinsic subtype status according to the PAM50 GS  was available for 66 of 85 samples, among which we classified 22 as luminal A, 14 as luminal B, 17 as normal, 3 as basal-like and 10 as HER2-enriched. There were 11 (50%) of 22 luminal A and 2 (14%) of 14 luminal B mutation cases (P = 0.039) (Table 1). Thus, 11 (85%) of the 13 mutation cases were categorised as luminal A.
Signatures of phosphatidylinositol 3-kinase pathway activity
Three GS scores were calculated for each sample at baseline. The PIK3CA-GS (Loi) was derived from exon 20 mutations and has been shown to negatively correlate with GSs of proliferation, AKT activation and PTEN loss whilst positively correlating with ESR1 expression . The PTEN-loss-GS (Saal) was developed to represent IHC-detectable PTEN loss in BC . Its association with PIK3CA mutation status suggests that the signature integrates various lesions in the PI3K pathway. The PI3K-GS (Creighton) is based on the set of genes in which expression was induced or repressed by PI3K inhibitors . As measurements of the same pathway, it might be expected that the signatures would correlate with one another. However, although the PI3K-GS (Creighton) and PTEN-loss-GS (Saal) were strongly positively correlated with one another (r = 0.657, P < 0.0001), the PIK3CA-GS (Loi) and PI3K-GS (Creighton) showed a weak, nonsignificant trend towards negative correlation (r = -0.218, P = 0.079), and the PIK3CA-GS (Loi) and PTEN-loss-GS (Saal) were minimally inversely correlated (r = -0.292, P = 0.017). The PIK3CA-GS (Loi) has been reported to correlate with ESR1, but we did not see this relationship (r = -0.122, P = 0.330). It did, however, show a slight trend towards correlation with ER protein expression (r = 0.240, P = 0.055).
PIK3CAmutation status and gene signatures representing the phosphatidylinositol 3-kinase pathway
The PIK3CA-GS (Loi) was strongly associated with mutation status, with mutants having higher median GS scores: WT = -0.051 (n = 50, IQR = -0.162 to 0.086) and Mu = 0.140 (n = 16, IQR = 0.062 to 0.254) (P < 0.0001) (Figure 3A). Because this signature was derived using exon 20 mutations, we ran the same analysis excluding the exon 20 mutations and found a significant relationship with non–exon 20 mutations: WT = -0.051 (n = 50, IQR = -0.162 to 0.086) and Mu = 0.143 (n = 11, IQR = 0.056 to 0.240) (P = 0.0003) (data not shown).
The PI3K-GS (Creighton) was not associated with PIK3CA mutation status. Median values were WT = 0.009 (n = 50, IQR = -0.034 to 0.036) and Mu = -0.007 (n = 16, IQR = -0.037 to 0.016) (P = 0.341) (Figure 3B). The PTEN-loss-GS (Saal) was significantly inversely associated with mutation status. Median values were WT = 0.025 (n = 50, IQR = -0.074 to 0.115) and Mu = -0.047 (n = 16, IQR = -0.137 to 0.029) (P = 0.041) (Figure 3C).
When we excluded PAGs from the GSs, the association with mutation status was maintained for the PIK3CA-GS (Loi) and the PI3K-GS (Creighton). For PIK3CA-GS (Loi), the median values were WT = -0.041 (n = 50, IQR = -0.166 to 0.095) and Mu = 0.151 (n = 16, IQR = 0.051 to 0.259) (P < 0.0001) (data not shown). For PI3K-GS (Creighton), the median values were WT = 0.006 (n = 50, IQR = -0.031 to 0.039) and Mu = -0.002 (n = 16, IQR = -0.035 to 0.021) (P = 0.372) (data not shown). However, the PTEN-loss-GS (Saal) association lost its statistical significance. The median values for PTEN-loss-GS (Saal) were WT = -0.002 (n = 50, IQR = -0.049 to 0.089) and Mu = -0.007 (n = 16, IQR = -0.081 to 0.041) (P = 0.503) (data not shown), indicating that this GS is influenced by proliferation genes (27.8% of genes in the GS are PAGs, compared to 12% for PIK3CA-GS and 15% for PI3K-GS).
Association of PIK3CAmutation status and gene signatures with Ki67
There was a borderline significant difference (P = 0.051) in baseline Ki67 between PIK3CA WT and mutation samples, with patients with the mutation having lower levels of proliferation (Table 2). We also analysed baseline Ki67 with regard to HS and non-HS mutations and found no evidence for a distinction to be made (P = 0.144).
There was no significant correlation between baseline Ki67 and the PIK3CA-GS (Loi); however, there was a significant positive correlation with the PI3K-GS (Creighton) and PTEN-loss-GS (Saal) (Table 3). Because the various GSs include PAGs, the influence of these PAGs on the signatures could be considered substantial because proliferation has been strongly associated with the predictive and prognostic utility of gene modules . Therefore, three modified GSs (noPAG version) which had PAGs removed were developed. In the absence of PAGs, the relationships of baseline Ki67 with the GSs were sustained.
Ki67 values after 2 weeks of treatment have been found to be associated with residual risk during endocrine treatment . We found a weak significant relationship of mutation patients’ having lower levels of 2-week Ki67 (P = 0.050) (Table 2), but no association with any of the GSs, with or without PAGs (Table 3).
Ki67 levels were suppressed by median values of 69% and 77% in patients with WT and mutation GSs, respectively, but this difference was not statistically significant (Table 2 and Additional file 6: Figure S1). Similarly, when we divided the patients into good and poor responders, using a cut-off value of 65%, we identified no difference in the proportion of good and poor responders in the WT and mutation groups (P = 0.095) (data not shown). In addition, there was no significant correlation between any of the three GSs and 2-week pre-Ki67% measurements (Table 3). We identified a slight trend towards significance for PIK3CA-noPAG-GS (r = 0.246, P = 0.067) (Table 3 and Additional file 7: Figure S2). There were ten HER2+ samples with mutation status and 2-week/pre-Ki67% data; however, the small sample size precluded any meaningful separate analysis of this subgroup.
Validation of Ki67 findings with regard to mutation status using a publicly available data set
We conducted an analysis of 310 baseline samples from the study by Ellis et al.  for which they had PIK3CA mutation data. In this much larger data set, we found a significant but weak association with baseline Ki67 and discovered that WT samples had higher Ki67 levels than mutation samples, with the following median values: WT = 18.6 (n = 180, IQR = 10.3 to 31.4) and Mu = 14.9 (n = 130, IQR = 7.3 to 26.7) (P = 0.043). This confirms the trend we identified in our smaller cohort. Unlike our data, there was no trend towards mutation samples’ having lower 2-week Ki67 values, with median data as follows: WT = 3.3 (n = 66, IQR = 1.2 to 8.3) and Mu = 3.4 (n = 40, IQR = 1.2 to 7.8) (P = 0.926). In our analysis of Ellis et al.’s data using our methodology for Ki67 change, we also did not identify a significant association, with the following median data: WT = 76% suppression (n = 66, IQR = 46 to 92) and Mu = 72% (n = 40, IQR = 10 to 87) (P = 0.334) .
Both de novo and acquired resistance to endocrine therapy is a major clinical problem in the treatment of ER + BC. In recent years, there has been a concerted effort to develop therapeutic agents targeting kinases within the PI3K pathway. This has grown out of the knowledge that deregulation of this pathway is associated with with increased tumour growth and resistance to therapy . Furthermore, in recent clinical studies, researchers have shown that the combination of everolimus, a targeted mTORC1 inhibitor, with the steroidal aromatase inhibitor exemestane provides greater benefit than exemestane alone in the third-line setting, increasing median progression-free survival from 4.1 to 10.6 months . These data have led to a concerted effort to identify biomarkers of response to PI3K inhibitors in an attempt to identify those patients most likely to gain benefit from this therapy.
Surprisingly, in multiple studies of PIK3CA in BC, researchers have shown PIK3CA to be a highly mutated gene in the luminal subtype that is associated with good clinicopathological markers. Despite this finding, its relationship with response to endocrine therapy remains unclear. In our present study, we aimed to evaluate the relationship of PIK3CA mutations and GSs representing PI3K pathway activation with key biological features of BC and short-term aromatase inhibitor–induced suppression of Ki67, a validated intermediate endpoint for benefit from endocrine therapy [11, 36]. To our knowledge, our present study is, to date, the largest single study of an aromatase inhibitor given in the neoadjuvant setting in which the relationship between response to treatment, PIK3CA mutations and multiple GSs of PI3K has been examined.
Our findings regarding the prevalence of PIK3CA mutations and the associations with lower risk markers in ER + BC confirm previously reported results [5, 7, 9, 13, 14, 37], thus providing further evidence that mutations in this gene seem to play a role in the biology of luminal BC. To complement the mutation data, and for a better overview of the activation status of the PI3K pathway, we investigated a set of previously published GSs. PIK3CA mutation status was associated with both the PIK3CA-GS (Loi) and PTEN-loss-GS (Saal) pathways; however, in each case, there was substantial overlap in the GSs between WT and Mu populations, such that GS could not act as an accurate surrogate for mutation status.
The opposing association of the two signatures representing high pathway activation, PI3K-GS (Creighton) and PTEN-loss-GS (Saal), with mutation status requires careful consideration of the degree to which they represent pathway activation and the population in which they are derived. They are far from equivalent, apparently reflecting different aberrations driving the signatures. The PI3K-GS, derived in vitro, was based on repression or induction as a result of treatment with PI3K inhibitors, whereas the PTEN-loss-GS recapitulates IHC-detectable loss and could be more appropriate for use in studies of basal-like tumours, in which it is associated with high pathway activity . With this in mind, it should also be noted that the GSs are indicative of relative pathway activation rather than absolute.
Although our analysis of PIK3CA mutations was comprehensive, the pathway consists of multiple nodes that could be subject to mutation (or to other aberrations, such as gain or loss of function), such that an alteration at one particular point might not have the same impact as another in terms of altering pathway activity.
Expression of Ki67 at baseline has been found to be associated with a good prognosis. In previous analyses of smaller cohorts, researchers have found that baseline Ki67 was not significantly different between WT and mutation samples overall or grouped by domain [12, 16]. In this study, however, we found a borderline trend for lower levels of Ki67 in mutation samples overall. This novel finding is confirmed by our analysis of a larger sample cohort from a public data set . The median and IQR Ki67 data in both our data set and that of Ellis et al. are extremely similar, suggesting that the weak significance in our data set is most likely due to the small sample size.
In two previous reports, researchers have described the effect of PIK3CA mutation status and response to treatment measured by Ki67 change in small patient cohorts with different treatments [12, 16]. They identified no significant association of mutation with response. We also found no relationship in our larger, homogeneously treated cohort, and, by investigating Ki67% change in a larger public sample set , we confirmed that this result was not a false-negative due to sample size. The lack of a relationship between PIK3CA mutation status and 2-week/pre-Ki67% in all the data sets suggests that the presence of a PIK3CA mutation has limited, if any, impact on the likelihood of benefit from presurgical anastrozole treatment in patients with primary BC.
Interestingly, in a neoadjuvant trial of ER + postmenopausal patients, researchers found that PIK3CA exon 9 mutations, compared with exon 20 mutations, were less responsive to letrozole alone, but that these patients were as sensitive as the overall population to everolimus and letrozole . However, the number of patients in that study was small. A pooled analysis of all published studies would allow greater exploration of the significance of the mutational subgroups because of increased power. The researchers in the recently reported BOLERO-2 and TAMRAD translational studies of ER + metastatic BC patients [39, 40] found that all patients derived benefit from the addition of everolimus to exemestane or tamoxifen, respectively, regardless of their PIK3CA genotype. However, the relevance of PIK3CA mutation to advanced, heavily pretreated disease may differ from primary disease. The early emergence of PIK3CA mutations in BC, as well as their high prevalence and their associations with high ER expression and good prognosis, suggests their functional importance in the disease, but it may be that other aberrations offer predictive potential. For example, patients with independent activation of mTOR via low levels of the tumour suppressor LKB1, an mTOR inhibitor, exhibited greater benefit from everolimus .
The strengths of this study are the comprehensive documentation of the presence of both common and rare PIK3CA mutations (>94% of all PIK3CA mutations recorded in BC) and the use of multiple GSs of the PI3K pathway to catalogue and compare aspects of the PI3K pathway. The use of the presurgical setting allowed the mutations to be linked to response to treatment by using a validated biomarker, in contrast to most other studies in which researchers have considered association with prognosis. We also verified the most important novel data, lower baseline Ki67 in mutation samples, in a large, publicly available data set. The moderate sample size of this study does mean, however, that the ability to determine the importance of the rare mutations is underpowered. The numerous borderline trends could come into significance in a larger data set, as shown by the analysis of the data from Ellis and colleagues.
The results of this presurgical study of 85 breast tumours treated with 2 weeks of anastrozole confirm the high prevalence of PIK3CA mutations and their relationship with markers of good clinical prognosis in ER + disease. The correlation of mutations with GSs of the PI3K pathway suggests that PIK3CA mutants exhibit lower mTORC1 signalling and reduced PI3K activity, despite their oncogenic nature in vitro. The data reported herein show that the presence of a PIK3CA mutation does not seem to attenuate the response to short-term treatment with anastrozole and therefore is not a biomarker of response as measured by Ki67, which is an important finding, given the high frequency of PIK3CA mutations and the widespread use of anastrozole treatment in postmenopausal women.
Epidermal growth factor receptor
Fluorescence in situ hybridisation
Human epidermal growth factor receptor
Matrix-assisted laser/desorption ionisation time of flight
- PIK3CA :
Phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit α
- PTEN :
Phosphatase and tensin homolog
Cancer.net Editorial Board; American Society of Clinical Oncology: Breast Cancer Guide. [http://www.cancer.net/cancer-types/breast-cancer/diagnosis] (accessed 3 July 2014)
Miller WR, Larionov AA: Understanding the mechanisms of aromatase inhibitor resistance. Breast Cancer Res. 2012, 14: 201-10.1186/bcr2931.
Miller TW, Rexer BN, Garrett JT, Arteaga CL: Mutations in the phosphatidylinositol 3-kinase pathway: role in tumor progression and therapeutic implications in breast cancer. Breast Cancer Res. 2011, 13: 224-10.1186/bcr3039.
Pérez-Tenorio G, Alkhori L, Olsson B, Waltersson MA, Nordenskjöld B, Rutqvist LE, Skoog L, Stål O: PIK3CA mutations and PTEN loss correlate with similar prognostic factors and are not mutually exclusive in breast cancer. Clin Cancer Res. 2007, 13: 3577-3584. 10.1158/1078-0432.CCR-06-1609.
Kalinsky K, Jacks LM, Heguy A, Patil S, Drobnjak M, Bhanot UK, Hedvat CV, Traina TA, Solit D, Gerald W, Moynahan ME: PIK3CA mutation associates with improved outcome in breast cancer. Clin Cancer Res. 2009, 15: 5049-5059. 10.1158/1078-0432.CCR-09-0632.
Loi S, Haibe-Kains B, Majjaj S, Lallemand F, Durbecq V, Larsimont D, Gonzalez-Angulo AM, Pusztai L, Symmans WF, Bardelli A, Ellis P, Tutt AN, Gillett CE, Hennessy BT, Mills GB, Phillips WA, Piccart MJ, Speed TP, McArthur GA, Sotiriou C: PIK3CA mutations associated with gene signature of low mTORC1 signaling and better outcomes in estrogen receptor–positive breast cancer. Proc Natl Acad Sci U S A. 2010, 107: 10208-10213. 10.1073/pnas.0907011107.
Cizkova M, Susini A, Vacher S, Cizeron-Clairac G, Andrieu C, Driouch K, Fourme E, Lidereau R, Bièche I: PIK3CA mutation impact on survival in breast cancer patients and in ERα, PR and ERBB2-based subgroups. Breast Cancer Res. 2012, 14: R28-10.1186/bcr3113.
Stemke-Hale K, Gonzalez-Angulo AM, Lluch A, Neve RM, Kuo WL, Davies M, Carey M, Hu Z, Guan Y, Sahin A, Symmans WF, Pusztai L, Nolden LK, Horlings H, Berns K, Hung MC, van de Vijver MJ, Valero V, Gray JW, Bernards R, Mills GB, Hennessy BT: An integrative genomic and proteomic analysis of PIK3CA, PTEN, and AKT mutations in breast cancer. Cancer Res. 2008, 68: 6084-6091. 10.1158/0008-5472.CAN-07-6854.
Ellis MJ, Ding L, Shen D, Luo J, Suman VJ, Wallis JW, Van Tine BA, Hoog J, Goiffon RJ, Goldstein TC, Ng S, Lin L, Crowder R, Snider J, Ballman K, Weber J, Chen K, Koboldt DC, Kandoth C, Schierding WS, McMichael JF, Miller CA, Lu C, Harris CC, McLellan MD, Wendl MC, DeSchryver K, Allred DC, Esserman L, Unzeitig G, et al: Whole-genome analysis informs breast cancer response to aromatase inhibition. Nature. 2012, 486: 353-360.
Barbareschi M, Buttitta F, Felicioni L, Cotrupi S, Barassi F, Del Grammastro M, Ferro A, Dalla Palma P, Galligioni E, Marchetti A: Different prognostic roles of mutations in the helical and kinase domains of the PIK3CA gene in breast carcinomas. Clin Cancer Res. 2007, 13: 6064-6069. 10.1158/1078-0432.CCR-07-0266.
Ellis MJ, Suman VJ, Hoog J, Lin L, Snider J, Prat A, Parker JS, Luo J, DeSchryver K, Allred DC, Esserman LJ, Unzeitig GW, Margenthaler J, Babiera GV, Marcom PK, Guenther JM, Watson MA, Leitch M, Hunt K, Olson JA: Randomized phase II neoadjuvant comparison between letrozole, anastrozole, and exemestane for postmenopausal women with estrogen receptor–rich stage 2 to 3 breast cancer: clinical and biomarker outcomes and predictive value of the baseline PAM50-based intrinsic subtype—ACOSOG Z103. J Clin Oncol. 2011, 29: 2342-2349. 10.1200/JCO.2010.31.6950.
Loi S, Michiels S, Baselga J, Bartlett JM, Singhal SK, Sabine VS, Sims AH, Sahmoud T, Dixon JM, Piccart MJ, Sotiriou C: PIK3CA genotype and a PIK3CA mutation-related gene signature and response to everolimus and letrozole in estrogen receptor positive breast cancer. PLoS One. 2013, 8: e53292-10.1371/journal.pone.0053292.
Loi S, Michiels S, Lambrechts D, Fumagalli D, Claes B, Kellokumpu-Lehtinen PL, Bono P, Kataja V, Piccart MJ, Joensuu H, Sotiriou C: Somatic mutation profiling and associations with prognosis and trastuzumab benefit in early breast cancer. J Natl Cancer Inst. 2013, 105: 960-967. 10.1093/jnci/djt121.
Ramirez-Ardilla A, Helmjr JC, Look MP, Lurkin I, Ruigrok-Ritstier K, van Laere S, Dirix L, Sweep FC, Span PN, Linn SC, Foekens JA, Sleijfer S, Berns EM, Jansen MP: Hotspot mutations in PIK3CA associate with first-line treatment outcome for aromatase inhibitors but not for tamoxifen. Breast Cancer Res Treat. 2013, 139: 39-49. 10.1007/s10549-013-2529-7.
Levine DA, Bogomolniy F, Yee CJ, Lash A, Barakat RR, Borgen PI, Boyd J: Frequent mutation of the PIK3CA gene in ovarian and breast cancers. Clin Cancer Res. 2005, 11: 2875-2878. 10.1158/1078-0432.CCR-04-2142.
Ellis MJ, Lin L, Crowder R, Tao Y, Hoog J, Snider J, Davies S, DeSchryver K, Evans DB, Steinseifer J, Bandaru R, Liu W, Gardner H, Semiglazov V, Watson M, Hunt K, Olson J, Baselga J: Phosphatidyl-inositol-3-kinase α catalytic subunit mutation and response to neoadjuvant endocrine therapy for estrogen receptor positive breast cancer. Breast Cancer Res Treat. 2010, 119: 379-390. 10.1007/s10549-009-0575-y.
Berns K, Horlings HM, Hennessy BT, Madiredjo M, Hijmans EM, Beelen K, Linn SC, Gonzalez-Angulo AM, Stemke-Hale K, Hauptmann M, Beijersbergen RL, Mills GB, van de Vijver MJ, Bernards R: A functional genetic approach identifies the PI3K pathway as a major determinant of trastuzumab resistance in breast cancer. Cancer Cell. 2007, 12: 395-402. 10.1016/j.ccr.2007.08.030.
Dowsett M, Smith IE, Ebbs SR, Dixon JM, Skene A, A’Hern R, Salter J, Detre S, Hills M: Walsh G; IMPACT Trialists Group: Prognostic value of Ki67 expression after short-term presurgical endocrine therapy for primary breast cancer. J Natl Cancer Inst. 2007, 99: 167-170. 10.1093/jnci/djk020.
Smith IE, Walsh G, Skene A, Llombart A, Mayordomo JI, Detre S, Salter J, Clark E, Magill P, Dowsett M: A phase II placebo-controlled trial of neoadjuvant anastrozole alone or with gefitinib in early breast cancer. J Clin Oncol. 2007, 25: 3816-3822. 10.1200/JCO.2006.09.6578.
Hammond ME, Hayes DF, Dowsett M, Allred DC, Hagerty KL, Badve S, Fitzgibbons PL, Francis G, Goldstein NS, Hayes M, Hicks DG, Lester S, Love R, Mangu PB, McShane L, Miller K, Osborne CK, Paik S, Perlmutter J, Rhodes A, Sasano H, Schwartz JN, Sweep FC, Taube S, Torlakovic EE, Valenstein P, Viale G, Visscher D, Wheeler T, Williams RB, et al: American Society of Clinical Oncology/College of American Pathologists guideline recommendations for immunohistochemical testing of estrogen and progesterone receptors in breast cancer. J Clin Oncol. 2010, 28: 2784-2795. 10.1200/JCO.2009.25.6529. A published erratum appears in J Clin Oncol 2010, 28:3543
Dunbier AK, Ghazoui Z, Anderson H, Salter J, Nerurkar A, Osin P, A’Hern R, Miller WR, Smith IE, Dowsett M: Molecular profiling of aromatase inhibitor–treated postmenopausal breast tumors identifies immune-related correlates of resistance. Clin Cancer Res. 2013, 19: 2775-2786. 10.1158/1078-0432.CCR-12-1000.
Dunbier AK, Anderson H, Ghazoui Z, Salter J, Parker JS, Perou CM, Smith IE, Dowsett M: Association between breast cancer subtypes and response to neoadjuvant anastrozole. Steroids. 2011, 76: 736-740. 10.1016/j.steroids.2011.02.025.
COSMIC database. [http://www.sanger.ac.uk/genetics/CGP/cosmic/] (accessed 2 July 2014)
National Cancer Institute, National Genome Research Institute, National Institutes of Health: The Cancer Genome Atlas data portal. [http://tcga-data.nci.nih.gov/docs/somatic_mutations/tcga_mutations.htm] (accessed 2 July 2014)
Shukla N, Ameur N, Yilmaz I, Nafa K, Lau CY, Marchetti A, Borsu L, Barr FG, Ladanyi M: Oncogene mutation profiling of pediatric solid tumors reveals significant subsets of embryonal rhabdomyosarcoma and neuroblastoma with mutated genes in growth signaling pathways. Clin Cancer Res. 2012, 18: 748-757. 10.1158/1078-0432.CCR-11-2056.
Saal LH, Johansson P, Holm K, Gruvberger-Saal SK, She QB, Maurer M, Koujak S, Ferrando AA, Malmström P, Memeo L, Isola J, Bendahl PO, Rosen N, Hibshoosh H, Ringnér M, Borg Å, Parsons R: Poor prognosis in carcinoma is associated with a gene expression signature of aberrant PTEN tumor suppressor pathway activity. Proc Natl Acad Sci U S A. 2007, 104: 7564-7569. 10.1073/pnas.0702507104.
Creighton CJ, Fu X, Hennessy BT, Casa AJ, Zhang Y, Gonzalez-Angulo AM, Lluch A, Gray JW, Brown PH, Hilsenbeck SG, Osborne CK, Mills GB, Lee AV, Schiff R: Proteomic and transcriptomic profiling reveals a link between the PI3K pathway and lower estrogen-receptor (ER) levels and activity in ER+ breast cancer. Breast Cancer Res. 2010, 12: R40-10.1186/bcr2594.
Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, Weinberg RA: An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet. 2008, 40: 499-507. 10.1038/ng.127.
Whitfield ML, Sherlock G, Saldanha AJ, Murray JI, Ball CA, Alexander KE, Matese JC, Perou CM, Hurt MM, Brown PO, Botstein D: Identification of genes periodically expressed in the human cell cycle and their expression in tumors. Mol Biol Cell. 2002, 13: 1977-2000. 10.1091/mbc.02-02-0030..
Ghazoui Z, Buffa FM, Dunbier AK, Anderson H, Dexter T, Detre S, Salter J, Smith IE, Harris AL, Dowsett M: Close and stable relationship between proliferation and a hypoxia metagene in aromatase inhibitor–treated ER-positive breast cancer. Clin Cancer Res. 2011, 17: 3005-3012. 10.1158/1078-0432.CCR-10-1704. A published erratum appears in Clin Cancer Res 2011, 17:4915
Desmedt C, Haibe-Kains B, Wirapati P, Buyse M, Larsimont D, Bontempi G, Delorenzi M, Piccart M, Sotiriou C: Biological processes associated with breast cancer clinical outcome depend on the molecular subtypes. Clin Cancer Res. 2008, 14: 5158-5165. 10.1158/1078-0432.CCR-07-4756.
Hu Z, Fan C, Oh DS, Marron JS, He X, Qaqish BF, Livasy C, Carey LA, Reynolds E, Dressler L, Nobel A, Parker J, Ewend MG, Sawyer LR, Wu J, Liu Y, Nanda R, Tretiakova M, Ruiz Orrico A, Dreher D, Palazzo JP, Perreard L, Nelson E, Mone M, Hansen H, Mullins M, Quackenbush JF, Ellis MJ, Olopade OI, Bernard PS, et al: The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genomics. 2006, 7: 96-10.1186/1471-2164-7-96.
Venet D, Dumont JE, Detours V: Most random gene expression signatures are significantly associated with breast cancer outcome. PLoS Comput Biol. 2011, 7: e1002240-10.1371/journal.pcbi.1002240.
Miller TW, Balko JM, Arteaga CL: Phosphatidylinositol 3-kinase and antiestrogen resistance in breast cancer. J Clin Oncol. 2011, 29: 4452-4461. 10.1200/JCO.2010.34.4879.
Baselga J, Campone M, Piccart M, Burris HA, Rugo HS, Sahmoud T, Noguchi S, Gnant M, Pritchard KI, Lebrun F, Beck JT, Ito Y, Yardley D, Deleu I, Perez A, Bachelot T, Vittori L, Xu Z, Muckopadhyay P, Lebwohl D, Hortobagyi GN: Everolimus in postmenopausal hormone-receptor-positive advanced breast cancer. N Engl J Med. 2012, 366: 520-529. 10.1056/NEJMoa1109653.
Dowsett M, Smith IE, Ebbs SR, Dixon JM, Skene A, Griffith C, Boeddinghaus I, Salter J, Detre S, Hills M, Ashley S, Francis S, Walsh G: IMPACT Trialists: Short-term changes in Ki-67 during neoadjuvant treatment of primary breast cancer with anastrozole or tamoxifen alone or combined correlate with recurrence-free surviva. Clin Cancer Res. 2005, 11: 951s-958s.
The Cancer Genome Atlas Network: Comprehensive molecular portraits of human breast tumours. Nature. 2012, 490: 61-70. 10.1038/nature11412.
Baselga J, Semiglazov V, van Dam P, Manikhas A, Bellet M, Mayordomo J, Campone M, Kubista E, Greil R, Bianchi G, Steinseifer J, Molloy B, Tokaji E, Gardner H, Phillips P, Stumm M, Lane HA, Dixon JM, Jonat W, Rugo HS: Phase II randomized study of neoadjuvant everolimus plus letrozole compared with placebo plus letrozole in patients with estrogen receptor–positive breast cancer. J Clin Oncol. 2009, 27: 2630-2637. 10.1200/JCO.2008.18.8391.
Treilleux I, Arnedos M, Cropet C, Ferrero JM, Abadie-Lacourtoisie S, Spaeth D, Levy C, Legouffe E, Pujade-Lauraine E, Wang Q, Bachelot T: Predictive markers of everolimus efficacy in hormone receptor positive (HR+) metastatic breast cancer (MBC): final results of the TAMRAD trial translational study. J Clin Oncol. 2013, 31: Abstract 510-10.1200/JCO.2012.46.0980.
Hortobagyi GN, Piccart-Gebhart MJ, Rugo HS, Burris HA, Campone M, Noguchi S, Perez AT, Deleu I, Shtivelband M, Provencher L, Masuda N, Dakhill SR, Anderson I, Chen D, Damask A, Huang A, McDonald R, Taran T, Sahmoud T, Baselga J: Correlation of molecular alterations with efficacy of everolimus in hormone receptor–positive, HER2-negative advanced breast cancer: results from BOLERO-2. J Clin Oncol. 2013, 31: Abstract LBA509-
We acknowledge past and present members of the Academic Department of Biochemistry, Royal Marsden Hospital, London, for sample acquisition, DNA extraction (particularly Helen Anderson) and IHC measurements. Thanks also to Maryou Lambros who provided technical guidance. We acknowledge support from the National Institute for Health Research RM/ICR Biomedical Research Centre.
The authors declare that they have no competing interests.
ELK extracted samples, ran the Sequenom analysis, analysed the data and wrote the manuscript. CS analysed the data and wrote the manuscript. QG provided the gene signatures. IGM and NT ran the ddPCR analysis and interpreted the data. LAM and MD conceived the study, participated in its design and coordination, and helped to draft the manuscript. IS designed the clinical study for 0223. All authors read and approved the final manuscript.
Elena López-Knowles, Corrinne V Segal contributed equally to this work.
Electronic supplementary material
Additional file 3: Table S3: Genes comprising the signatures of Creighton et al., Loi et al. and Saal et al. The coefficient describes whether the gene contributes positively or negatively to the gene signature (GS) score. “X” denotes which GS the gene was used in. The red font indicates that the gene was not mapped to the Illumina array data. The total counts for the genes that were mapped from the original publication to the Illumina array data are as follows: 176 of 181 for PTEN-loss-GS (Saal), 231 of 236 for PIK3CA-GS (Loi) and 1,733 of 1,793 for PI3K-GS (Creighton). (XLSX 91 KB)
Additional file 6: Figure S1: 2 wk/pre-Ki67% in wild-type (WT) and mutation (Mu) samples. Data are medians with IQRs (n = 73). (JPEG 503 KB)
About this article
Cite this article
López-Knowles, E., Segal, C.V., Gao, Q. et al. Relationship of PIK3CA mutation and pathway activity with antiproliferative response to aromatase inhibition. Breast Cancer Res 16, R68 (2014). https://doi.org/10.1186/bcr3683
- Breast Cancer