Skip to main content
  • Research article
  • Open access
  • Published:

A functional proteogenomic analysis of endometrioid and clear cell carcinomas using reverse phase protein array and mutation analysis: protein expression is histotype-specific and loss of ARID1A/BAF250a is associated with AKT phosphorylation

Abstract

Background

Ovarian cancer is now recognized as a number of distinct diseases primarily defined by histological subtype. Both clear cell ovarian carcinomas (CCC) and ovarian endometrioid carcinomas (EC) may arise from endometriosis and frequently harbor mutations in the ARID1A tumor suppressor gene. We studied the influence of histological subtype on protein expression with reverse phase protein array (RPPA) and assessed proteomic changes associated with ARID1A mutation/BAF250a expression in EC and CCC.

Methods

Immunohistochemistry (IHC) for BAF250a expression was performed on 127 chemotherapy-naive ovarian carcinomas (33 CCC, 29 EC, and 65 high-grade serous ovarian carcinomas (HGSC)). Whole tumor lysates were prepared from frozen banked tumor samples and profiled by RPPA using 116 antibodies. ARID1A mutations were identified by exome sequencing, and PIK3CA mutations were characterized by MALDI-TOF mass spectrometry. SAM (Significance Analysis of Microarrays) was performed to determine differential protein expression by histological subtype and ARID1A mutation status. Multivariate logistic regression was used to assess the impact of ARID1A mutation status/BAF250a expression on AKT phosphorylation (pAKT). PIK3CA mutation type and PTEN expression were included in the model. BAF250a knockdown was performed in 3 clear cell lines using siRNA to ARID1A.

Results

Marked differences in protein expression were observed that are driven by histotype. Compared to HGSC, SAM identified over 50 proteins that are differentially expressed in CCC and EC. These included PI3K/AKT pathway proteins, those regulating the cell cycle, apoptosis, transcription, and other signaling pathways including steroid hormone signaling. Multivariate models showed that tumors with loss of BAF250a expression showed significantly higher levels of AKT-Thr308 and AKT-Ser473 phosphorylation (p < 0.05). In 31 CCC cases, pAKT was similarly significantly increased in tumors with BAF250a loss on IHC. Knockdown of BAF250a by siRNA in three CCC cell lines wild type for ARID1A showed no increase in either pAKT-Thr308 or pAKT-S473 suggesting that pAKT in tumor tissues is indirectly regulated by BAF250a expression.

Conclusions

Proteomic assessment of CCC and EC demonstrates remarkable differences in protein expression that are dependent on histotype, thereby further characterizing these cancers. AKT phosphorylation is associated with ARID1A/BAF250a deficient tumors, however in ovarian cancers the mechanism remains to be elucidated.

Peer Review reports

Background

There are important clinical and molecular differences between histological subtypes of ovarian cancer [1, 2]. High-grade serous carcinoma (HGSC) is the most common subtype of ovarian cancer representing approximately 70% of cases, whereas clear cell carcinoma (CCC) and endometrioid carcinoma (EC) are less common, occurring at a frequency of approximately 12% and 11% respectively [3]. Whereas HGSC frequently originates from precursor lesions in the fallopian tube [4, 5], both CCC and EC may arise from endometriosis and are associated with a younger age of onset [6–10]. Compared to HGSC, both CCC and EC present with earlier-stage disease, and CCC are less sensitive to common chemotherapy regimens used to treat ovarian cancers [11]. Approximately 50% of CCC and 30% of EC harbor mutations in the chromatin-remodeling gene ARID1A[12, 13]. ARID1A functions as a tumor suppressor gene [14], although there are no proven clinical features or differences in outcomes associated with a reduction or loss of expression of its associated protein (BAF250a) by immunohistochemistry (IHC) [12, 15].

PIK3CA mutations are common in ovarian CCC, endometrial and breast cancers [16–19] and they are frequently co-mutated with ARID1A in ovarian CCC [13]. The impact of these co-mutations on phosphorylation of AKT has been studied in endometrial cancers and it was found that tumors with mutated ARID1A were associated with increased AKT phosphorylation and pathway activation [20]. In this same study, PTEN mutations were also significantly associated with mutations in ARID1A. Further research is needed in order to determine whether these observations in endometrial cancer also pertain to ovarian CCC and EC. Therefore, the purpose of this study was to examine protein expression in ovarian EC and CCC using RPPA and to determine whether expression patterns are specific to histotype. We also sought to determine proteomic patterns of expression based on ARID1A mutation status/BAF250a protein expression in ovarian EC and CCC, and further study the impact of mutation status (ARID1A and PIK3CA) and PTEN expression on AKT phosphorylation.

Methods

Reverse phase protein array

We performed RPPA on whole tumor lysates as previously described [21, 22]. RPPA is suited for functional proteomic assessments of a large number of tumor samples and provides a platform for comparing the relative protein expression between these samples. Using RPPA, 116 antibodies were utilized to assay protein expression of cell surface growth factor receptors, common signaling pathway proteins, steroid hormone receptors, and others proteins involved in proliferation and apoptosis. RPPA provides an estimate of the relative protein expression of each sample in relation to the other samples on the slide. A list of antibodies used for RPPA is included as Additional file 1: Table S1.

Collection of tumor samples and IHC

Tumor samples (n = 127) were obtained from the gynecologic tumor bank at Vancouver General Hospital and the British Columbia Cancer Agency. The research was conducted with approval from the University of British Columbia institutional review board. Tumor samples were collected at the time of primary surgery and snap frozen within 60 minutes after collection. Patients treated by neoadjuvant chemotherapy were excluded from the study. All patient samples were subjected to pathology review (BG) to confirm the histological subtype and site of origin. Clinical data was accessed through the Cheryl Brown Ovarian Cancer Outcomes Unit and is updated on a regular basis. Tumors were studied by IHC as part of a tumor bank tissue microarray (TMA). Methods for preparation, staining, and scoring of the EOC cases for BAF250a have been previously described [12]; tumors with any degree of nuclear staining for BAF250a protein expression considered a positive score for purposes of this study.

BAF250a knockdown experiments in clear cell lines

Cell line identities were verified with short tandem repeat (STR) DNA identity testing with the AmpFLSTR® Identifiler® Direct PCR Amplification Kit as per the manufacturer’s protocol (Applied Biosystems®, Life Technologies Corp.) ES2 was obtained from the American Type Culture Collection (ATCC); STR typing to confirm the cell line identity was performed May 12, 2011 and verified against the ATCC STR database. JHOC-5 was obtained from the RIKEN cell bank; STR typing was performed on April 18, 2011 and verified against the RIKEN DSMZ STR Profiler. RMG-1 was obtained from the Health Sciences Research Resources Bank (HSRRB) cell bank; short tandem repeat (STR) DNA identity testing was performed on April 18, 2011 and verified against the HSRRB STR database. Cells were cultured in RPMI with 5-10% FBS, with the exception of the ES-2 cells, which were grown in McCoy’s media. Cells were treated with pooled siRNA (Dharmacon) to ARID1A at 20nM using RNAiMax™ as a transfection reagent (Invitrogen) according to the manufacturer’s recommended protocol. Lysates were prepared 60 hours following transfection using Bicine/CHAPS lysis buffer (20 mM Bicine/0.6% CHAPS, pH 7.6). EGF stimulation was performed using 20 ng/ml EGF for 15–20 minutes as a positive control. Protein levels were assessed by Western blotting (SDS-PAGE). Antibodies used were as follows: BAF250a (Sigma- HPA005456), AKTp473 (Cell Signaling 9271, AKTp308 (Cell Signaling 9275), PDK1 (Cell Signaling 3062), p70S6 Kinase (Santa Cruz sc-8418), PTEN (Cell Signaling 9552), and β-actin (Cell Signaling 3700).

Isoelectric focusing for AKT protein expression

Native capillary isoelectric point focusing was used to assess AKT expression according to recommended protocols using NanoPro™ 2000 (ProteinSimple™, Santa Clara, CA) [23]. G2 premix (pH 5–8) was used for all experiments (ProteinSimple™, # 040–972) Primary antibodies included AKT1 antibody (Santa Cruz, (C-20): sc1618), and ERK 1/2 (ProteinSimple™, Catalog #040-474). Secondary antibodies used were bovine anti-goat (Jackson Labs, Cat# 805-035-180), and goat anti-rabbit (Sigma, #A0545).

DNA sequencing of tumors and mutational analyses

The ARID1A mutations were determined primarily by exome sequencing using next generation technologies and these results were previously described by our group [12]. Common oncogenic mutations were determined by MALDI-TOF mass spectrometry (MassARRAY®, Sequenom Inc.) of primer extension polymerase chain reaction (PCR) products [24]. Mutations of PIK3CA were categorized as kinase domain mutation H1047R, and helical domain mutations (E542K, E545K, and also E545D and Q546R). In addition, there were several other less common mutations included (N345K, R88Q, E110K).

Statistical methods

To assess tumor protein expression unsupervised and supervised cluster analysis was performed using TreeView software (University of Glasgow, Scotland). X-cluster was used to generate heat maps and cluster groups (Stanford). Data was analyzed using SPSS software (Version 20, Chicago, Illinois). Proteins that were differentially expressed in CCC relative to the other subtypes were determined by t-testing and significance analysis of microarray data (SAM) [25]. SAM was used to determine the probability that differences in expression within each group studied are due to false discovery. False discovery rates (FDR) of less 5% were considered significant. Median protein expression levels were used as cut-points for multivariate logistic regression models except for PTEN expression. A cut-point of the lowest 20% of values on RPPA was chosen for PTEN loss based on previous work showing a 28% incidence of PTEN loss in CCC [26], and IHC data from our own center showing a 12% incidence of PTEN loss in 42 cases of CCC/EC.

Results

EOC protein expression is histotype specific

127 ovarian cancers were submitted for RPPA: 33 CCC, 29 EC, and 65 HGSC. Hierarchical clustering of treatment-naïve samples and proteins analyzed by RPPA is shown in Figure 1. It is evident that clustering is driven by histological subtype; CCC and EC form distinct clusters separated by two major HGSC subgroups. Differential expression of proteins was examined by significance analysis of microarray data (SAM) for both CCC and EC compared to HGSC as the reference group. Figure 2 shows the differentially expressed proteins according to histological subtype. Interestingly a number of proteins are down-regulated in both EC and CCC relative to the HGS group. Notable proteins are GAB2, YAP, Cyclin B1, CofilinpSer3, and Caveolin1, 4EBP1pThr37, STAT5, and c-Myc. CCC have lower ERα, AR, and PR expression compared to HGSC. The complete list of protein comparisons by SAM and histotype is included as Additional file 2: Table S2. Other proteins with low expression in CCC include β-Catenin, fibronectin, PTEN, 4EBP1, and stathmin. EC are characterized by lower levels of Cyclin E1, phosphorylated S6, 14-3-3 Zeta, BIM, and STAT3.

Figure 1
figure 1

Hierarchical clustering of samples and proteins analyzed by RPPA. CCC and ECs form distinct clusters separated by two major high-grade serous carcinoma (HGSC) subgroups. Differential expression of proteins was examined by significance analysis of microarray data (SAM analysis) for both CCC and EC compared to HGSC as the reference group.

Figure 2
figure 2

SAM Analysis of differentially expressed proteins by histotype: A. Upregulated and B. Down-regulated Proteins in Clear Cell and Endometrioid Carcinoma by SAM analysis on RPPA data, with High Grade Serous Carcinoma (HGSC) as the reference group. Proteins with FDR* of <5% are shown. Proteins that are over or under expressed in both CCC and EC are shown in the center between the two subtypes. Note: *FDR (False discovery rate): For each analysis (e.g. CCC vs. EC), the percent probability that the difference in protein expression for each comparison is attributable to chance after correcting for multiple parameter testing.

CCC and EC shared four proteins that were up-regulated compared to HGSC. Notably pAKT-Thr308 is higher in both histologies though pAKT-Ser473 was of borderline significance. (FDR = 5.2 and 6.9 respectively). Other up-regulated proteins in both CCC and EC include p21, HSP70, and 4EBP1-pSer65. EC as expected have higher expression of ERα, ERαpSer118 and PR compared to HGSC, while IRS1, CHK1, EGFR, and Cyclin D1 were also elevated. Some of the up-regulated proteins that characterize CCC include α-Catenin, Cyclin E1, HSP27, E-cadherin, p38, p38pThr180Tyr182, SMAD3, and GSK proteins.

We then examined differential protein expression associated with BAF250a IHC score by SAM in the 31 CCC with known ARID1A mutation status. There were only 2 cases of EC with negative BAF250a IHC so the EC were not included in this analysis. The only differentially expressed protein with a low false discovery rate (FDR) was pAKT-Thr308 (FDR < 1%). The top 4 other proteins differentially expressed were as follows: Bcl-2, p27, phosphorylated p38, and pAKT-Ser473. FDR rates by SAM for these 5 proteins varied from 16-20%. AKT phosphorylation levels were higher in BAF250a IHC negative tumors while the other 3 proteins had lower levels.

Investigating the effect of BAF250a/ARID1A and PIK3CAmutation on protein expression

We examined the effect of ARID1A mutation/BAF250a loss on changes in protein expression in the CCC and EC tumors with known ARID1A mutation status, BAF250a expression by IHC, and PIK3CA mutation status. Table 1 shows the characteristics of this study population. Complete sequencing data for ARID1A and PIK3CA mutation status, and BAF250a expression by IHC were available on 90 of the 127 original cases used for proteomic analysis. ARID1A mutations were present in 17 of 31 CCC in this study (55%), and 5 of 24 (21%) cases of EC. Sequenom MassARRAY® analysis showed PIK3CA mutations in 45% of CCC and 46% of EC. Mutations of PIK3CA were categorized as kinase domain (H1047R), helical domain (E542K, E545K, E545D and Q546R), or other (less common) including (N345K, R88Q, and E110K). Fifty-seven percent of tumors (13/23) with ARID1A mutations had PIK3CA mutations versus 38% (12/32) percent of tumors wild type for ARID1A (p = 0.3; Chi-square test). Since previous reports have shown that PIK3CA mutations have differing effects on PIK3CA signaling and AKT phosphorylation, we investigated the role of different PIK3CA mutations in relation to ARID1A mutation or loss by IHC. While 8/11 (73%) PIK3CA helical domain mutants were wild type for ARID1A; kinase/other domain mutants were more common in tumors with ARID1A mutations (10/14; 71%), (p = 0.05; Fisher’s exact test). Figure 3 shows a heat map with samples categorized according to BAF250a IHC result and PIK3CA mutation status. It is evident there are no obvious clustering patterns due either BAF250a or PIK3CA mutation status.

Table 1 Patient/sample characteristics with known ARID1A mutation status according to histotype
Figure 3
figure 3

Heat map with samples categorized according to BAF250a IHC scores and PIK3CA mutation status. No obvious clustering patterns due to either BAF250a by IHC or PIK3CA mutation status are present. Note: capital letters for each sample histotype label (eg. CCC) indicates expression of BAF250a on IHC, small letters (e.g. ccc) indicates loss of BAF250a expression. For PIK3CA mutations, a + sign after the sample label (e.g. CCC+) indicates the presence of a PIK3CA mutation (either helical, kinase, or other) whereas a – sign indicates no mutation.

Loss of BAF250a by IHC, low PTEN expression and activating PIK3CAmutations (particularly in tumors without PTEN loss) are associated with elevated pAKT

Table 2 shows the uni- and multivariate analyses of AKT phosphorylation based on ARID1A mutation status, BAF250a IHC expression, PIK3CA mutation status, and PTEN levels by RPPA. On univariate testing, pAKT-Thr308 expression by RPPA was significantly higher in tumors that were BAF250a negative by IHC (p = 0.007), and those with low PTEN (p = 0.025). In contrast, neither ARID1A mutation status nor the presence of a PIK3CA mutation was associated with a change in pAKT-Thr308 levels. Changes in pAKT-Ser473 were not associated with any of the above factors on univariate testing.

Table 2 Effect of PIK3CA mutations, ARID1A mutations, and BAF250a Expression on Akt phosphorylation

In terms of clarifying the impact of the different PIK3CA mutations on pAKT, the distributions (box plots) of pAKT-Thr308 and pAKT-Ser473 were compared (Additional file 3: Figure S1). By ANOVA, mutation status was not associated with significantly different pAKT levels by RPPA. Helical domain mutants were not associated with an increase in pAKT levels whereas kinase domain and other mutants resulted in slightly higher mean levels of pAKT (both Thr308 and Ser473) that were not statistically significant. Loss of PTEN is another known mechanism for activation of PI3K pathway signaling. When the eleven low PTEN cases were excluded, kinase domain and other less common mutants had significantly higher pAKT-Thr308 (p < 0.05; t-test Bonferroni corrected) than wild type, while helical domain mutants did not. For pAKT-Ser473, only the uncommon mutants had significantly higher levels by RPPA (p < 0.05; t-test Bonferroni corrected) compared to wild type or helical domain mutants.

Multivariate logistic regression showed that pAKT-Thr308 levels by RPPA are significantly higher in tumors with loss of BAF250a expression (p = 0.002), those with low PTEN (p = 0.003), and cancers with activating PIK3CA mutations (p = 0.02). As with the univariate testing, pAKT-Thr308 levels were not associated with ARID1A mutation status. On multivariate analysis, increases in pAKT-Ser473 were only associated with BAF250a expression levels (p = 0.04). Subgroup analysis of AKT phosphorylation in tumors without PIK3CA mutations (18 CCC and 12 EC) showed statistically significant increases in both pAKT-Ser473 and pAKT-Thr308 for 8 cancers lacking BAF250a expression on IHC (p = 0.05 and 0.008 respectively; t-test). Multivariate logistic regression analysis on only the 31 clear cell cancers with known ARID1A mutation status showed both pAKT-Thr308 and pAKT-Ser473 were significantly higher in tumors with BAF250a loss by IHC (p = 0.05 and 0.005 respectively). The EC subgroup was not similarly analyzed, as there were too few ARID1A mutant/BAF250a negative cases for comparison.

siRNA knockdown of BAF250a does not effect AKT phosphorylation in CCC cell lines

In the BAF250a-expressing cell lines RMG1, JHOC5 and ES2, despite good knockdown of BAF250a, changes in AKT phosphorylation (pAKT) or p70S6K, a downstream signaling protein of pAKT were not observed in any of the cell lines tested (Figure 4). Although pAKT-Thr308 was difficult to detect in the RMG1 cell lines, with EGF stimulation all cell lines showed an increase in pAKT except for pAKT-Ser473 in the JHOC5 cells, however changes in pAKT were not altered by BAF250a knockdown. We also observed that PDK1, PTEN levels did not change with BAF250a knockdown in any of the cell lines. Native protein AKT profiles assessed using capillary tube isoelectric point focusing (Figure 4B-D) are sensitive to phosphorylation events [27]. AKT profiles using this technology confirmed western blot results with no change in AKT profiles in two of the cells lines (RMG1 and JHOC5) and possibly a very slight increase in pAKT (isoforms 1 and 2 for the ES2 cell line). Taken together, these findings indicate that BAF250a knockdown has little to no effect in AKT phosphorylation in the cell lines tested.

Figure 4
figure 4

A) Western blot results from siRNA knockdown of BAF250a on cell lines ES2, JHOC5, and RMG1 to clarify the interaction between BAF250a expression and pAKT. Despite good knockdown of BAF250a no change in AKT phosphorylation or levels of p70S6K, a downstream signaling protein of pAKT can be seen in the ES2 and RMG1 cell lines. The baseline levels of pAKT are much higher in the JHOC5 cell line, and there is a suggestion of an increase in pAKT-Thr308 with BAF250a knockdown without obvious similar changes in pAKT-Ser473. PDK1 and PTEN levels did not change with BAF250a knockdown in any of the cell lines. B-D) Native protein AKT profiles using capillary tube isoelectric point focusing. Native AKT profiles are consistent with the western blot result in A, as little change occurs in AKT/pAKT following siRNA mediated BAF250a knockdown.

Clinical outcomes

Due to the limited cohort size, survival outcomes were not examined. Instead, FIGO stage distribution at presentation was analyzed to determine if this was influenced by BAF250a loss, PIK3CA, or ARID1A mutation status. For this comparison, FIGO stage was classified as low-stage (stage 1/2) or advanced-stage (stage 3/4). Those patients with BAF250a loss were more likely to have advanced stage (5/11 = 45% vs. 9/42 = 21%), however this difference was not statistically significant (p = 0.13). Eighty-four percent (21/25) of patients with a PIK3CA mutation presented with low-stage disease, compared to 64% (18/28; 2 cases with missing stage) without a mutation. Again, this difference was not statistically significant (p = 0.13). Changes in AKT phosphorylation were also not associated with stage.

Discussion

We used functional proteomics Reverse Phase Protein Array (RPPA) to assess differential protein expression in CCC and EC compared to HGSC. RPPA has been validated as a means of outcome prediction in breast and ovarian cancers [28, 29], and is therefore a useful technology. It is apparent that protein expression is clearly related to histotype, as has been previously shown both by IHC classification and RPPA [1, 30]. As expected EC has higher levels of steroid hormone receptor expression (AR, PR) while CCC has lower expression (AR, ER, PR) relative to HGSC. Figure 2 shows that here are more than 50 differentially expressed proteins that characterize EC and CCC from HGSC. It is important to note that our study is descriptive with respect to these protein associations, as they have not been validated by other means.

This study provides evidence that BAF250a expression identifies a subgroup of CCC and EC with higher AKT-Thr308 and AKT-Ser473 phosphorylation. Based on RPPA expression levels, tumors lacking BAF250a by IHC show higher levels of AKT phosphorylation independent of PIK3CA mutation and PTEN loss. Similarly in the cohort of CCC, multivariate analysis showed that increased pAKT (pAKT-Thr308 and pAKT-Ser473) was associated with BAF250a expression loss. As previously shown in endometrial cancers, increased pAKT was also found in tumors with ARID1A/BAF250a loss when both PIK3CA mutation or PTEN loss were absent [20].

We previously reported that ARID1A mutation status correlates with BAF250a loss by IHC [12]. In this previous report 73% of CCC with an ARID1A mutation showed BAF250a loss by IHC, while the correlation was less (50%) in EC. This current study contained a subset of cases from the original report, and found that loss of BAF250a protein by IHC showed a stronger association with pAKT changes than actual ARID1A mutation status. Thus, the assessment of BAF250a by IHC may be preferable to sequencing ARID1A in tumor samples to identify associations with AKT signaling. Mutations in ARID1A are often heterozygous [13] and in such cases it is unknown to what extent BAF250a expression/function is impaired and therefore still detectable by IHC. This could also be true for some somatic mutations. Conversely, IHC can detect epigenetic loss of function events such as methylation that may be as biologically deleterious as homozygous nonsense mutations. These same considerations have previously been described in relation to PTEN assessment. IHC appears to be preferable to sequencing for assessment of PTEN function in endometrial cancers, though a much greater proportion of cases (44%) show PTEN loss by IHC in the presence of a normal PTEN sequence [31]. Further research should be done to confirm the findings of this study and as our sample size is limited due to the fact that CCC and EC are relatively uncommon.

Unlike previously reported observations in endometrial cancers cell lines [20], we have been unable to demonstrate a mechanism by which AKT phosphorylation (pAKT) occurs in tumors using cell line models (RMG1, ES2, and JHOC5). Despite good knockdown of BAF250a in these lines, we could not demonstrate clear increases in pAKT, or downstream signaling. Our findings suggest cell line models may not accurately reflect the signaling changes in tumor samples. Other possible explanations include pAKT modulation by tumor/stromal interactions, cell lineage specificity, or the association of pAKT with other mutations such as PIK3CA regulatory domain mutations as they have also been shown to alter pAKT [31].

PIK3CA mutations are frequent in EC and CCC, as well as endometrial and breast cancers [16–19]. There is evidence indicating different PIK3CA mutations have differing effects on PIK3CA signaling and AKT phosphorylation, and therefore may determine prognosis. In breast cancer patients, a recent study suggests patients with PIK3CA mutations may actually have a more favorable prognosis [32], and that tumors with PIK3CA mutations irrespective of cancer type might have improved response rates to PI3K-directed therapies. In cell lines, helical domain mutants are not consistently associated with increases in pAKT and appear to act through alternate signaling mechanisms potentially involving SGK3 [33]. Similarly, we were not able to identify obvious increases in pAKT signaling in EOC patients with helical domain PIK3CA mutants.

This study provides important insights relating to protein expression particularly in ovarian CCC and EC. Because CCC and EC are relatively uncommon compared to high grade serous ovarian cancer (HGSC), it is difficult to gather large numbers of patient samples in order to have the power to detect multiple pathway aberrations. We were unable to show downstream pathway effects resulting from pAKT. Furthermore, selection and validation of antibodies for RPPA may potentially bias results so further study of a larger population of patients is important to confirm and clarify the associations we have found between AKT phosphorylation and BAF250a loss. The identification of downstream signaling aberrations also provides supporting evidence of the biological significance of pAKT. It will be helpful in future studies to also screen for other mutations such as regulatory domain mutations of PIK3CA as it may be that ARID1A loss is simply associated with other cellular events that increase pAKT rather than being a direct cause.

It is interesting to note that aside from pAKT, proteomic aberrations resulting from ARID1A/BAF250a loss in this study are otherwise absent. This finding supports clinical observations to date those patients with ARID1A mutant ovarian clear cell endometrial carcinomas do not have a phenotype associated with differing outcomes or treatment responses [12, 15, 34]. Recently however a report in patients with bladder cancer found that ARID1A/BAF250a deficient tumors are associated with a more aggressive phenotype [35]. These findings do not negate the importance of studying signaling pathway aberrations as this information may lead to finding new therapeutic targets. Furthermore, much work is being done to better understand the mechanism of action of novel AKT/PI3K directed therapies [36]. CCC is less sensitive to chemotherapy [37] and the addition of other chemotherapy agents such as irinotecan has been disappointing [38]. While histotype designation is key to classifying ovarian cancers, more research is also needed to find markers that identify patient subgroups who will benefit from novel therapies.

Conclusions

Using functional proteomics this study identified over 50 proteins that are differentially expressed when CCC and EC are compared to HGSC. Interestingly, functional protein aberrations resulting from the common occurrence of ARID1A mutation/BAF250a loss are limited. Only increases in AKT phosphorylation were found to be associated with BAF250a loss though our ability to detect changes in other PI3K pathway proteins may be limited based on sample size. Knockdown experiments in clear cell cancer lines failed to show a direct effect on AKT phosphorylation suggesting an indirect/alternative mechanism for AKT phosphorylation in CCC. Further study is needed to determine the mechanism for the association between BAF250a loss and pAKT in CCC and EC, and explore the diagnostic and therapeutic implications of the described changes in protein expression relating to histotype.

Note: For full protein names see antibody list (Additional file 1: Table S1).

Authors’ information

KCW was supported by a Dr. Nelly Auersperg Studentship award from the Michael Smith Foundation for Health Research (MSFHR). YW is a recipient of postdoctoral fellowship awards from Canadian Institute of Health Research (CIHR) and MSFHR. GBM and BTH were supported by NCI P50CA083639 and CCSG P30 CA16672. GBM is also supported by the Stand Up to Cancer Dream Team Translational Research Grant, a Program of the Entertainment Industry Foundation (SU2C-AACR-DT0209). BTH is additionally supported by a Career Development Award (CDA) from the Conquer Cancer Foundation of the American Society of Clinical Oncology (ASCO) and a Translational Research Award (TRA) from the Health Research Board (HRB) of Ireland and Science Foundation Ireland (SFI). DGH is supported by the MSFHR and the Dr. Chew Wei Memorial Professorship. MSC was the recipient of a Myer-Levy Fellowship and GBM, BTH, and MSC, were supported by the Kleberg Center for Molecular Markers.

Abbreviations

ATCC:

American type culture collection

CCC:

Ovarian clear cell carcinoma

EC:

Ovarian endometrioid carcinoma

FDR:

False discovery rate

HSRRB:

Health sciences research resources bank

HGSC:

High-grade serous ovarian cancer

IHC:

Immunohistochemistry

MALDI-TOF MS:

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry

pAKT:

AKT phosphorylation

PCR:

Polymerase chain reaction

RPPA:

Reverse phase protein array

SAM:

Significance analysis of microarrays

siRNA:

Small interfering RNA

STR:

Short tandem repeat

TMA:

Tissue microarray.

References

  1. Kobel M, Kalloger SE, Boyd N, McKinney S, Mehl E, Palmer C, Leung S, Bowen NJ, Ionescu DN, Rajput A, Prentice LM, Miller D, Santos J, Swenerton K, Gilks CB, Huntsman D: Ovarian carcinoma subtypes are different diseases: implications for biomarker studies. PLoS Med. 2008, 5 (12): e232-10.1371/journal.pmed.0050232.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Thigpen T, du Bois A, McAlpine J, DiSaia P, Fujiwara K, Hoskins W, Kristensen G, Mannel R, Markman M, Pfisterer J, Quinn M, Reed N, Swart AM, Berek J, Colombo N, Freyer G, Gallardo D, Plante M, Poveda A, Rubinstein L, Bacon M, Kitchener H, Stuart GC, Gynecologic Cancer InterGroup: First-line therapy in ovarian cancer trials. Int J Gynecol Cancer. 2011, 21 (4): 756-762. 10.1097/IGC.0b013e31821ce75d.

    Article  PubMed  Google Scholar 

  3. Kobel M, Kalloger SE, Huntsman DG, Santos JL, Swenerton KD, Seidman JD, Gilks CB, Cheryl Brown Ovarian Cancer Outcomes Unit of the British Columbia Cancer Agency VBC: Differences in tumor type in low-stage versus high-stage ovarian carcinomas. Int J Gynecol Pathol. 2010, 29 (3): 203-211. 10.1097/PGP.0b013e3181c042b6.

    Article  PubMed  Google Scholar 

  4. Crum CP, Drapkin R, Miron A, Ince TA, Muto M, Kindelberger DW, Lee Y: The distal fallopian tube: a new model for pelvic serous carcinogenesis. Curr Opin Obstet Gynecol. 2007, 19 (1): 3-9. 10.1097/GCO.0b013e328011a21f.

    Article  PubMed  Google Scholar 

  5. Kuhn E, Kurman RJ, Vang R, Sehdev AS, Han G, Soslow R, Wang TL: Shih IeM: TP53 mutations in serous tubal intraepithelial carcinoma and concurrent pelvic high-grade serous carcinoma–evidence supporting the clonal relationship of the two lesions. J Pathol. 2012, 226 (3): 421-426. 10.1002/path.3023.

    Article  CAS  PubMed  Google Scholar 

  6. LaGrenade A, Silverberg SG: Ovarian tumors associated with atypical endometriosis. Hum Pathol. 1988, 19 (9): 1080-1084. 10.1016/S0046-8177(88)80090-X.

    Article  CAS  PubMed  Google Scholar 

  7. Ogawa S, Kaku T, Amada S, Kobayashi H, Hirakawa T, Ariyoshi K, Kamura T, Nakano H: Ovarian endometriosis associated with ovarian carcinoma: a clinicopathological and immunohistochemical study. Gynecol Oncol. 2000, 77 (2): 298-304. 10.1006/gyno.2000.5765.

    Article  CAS  PubMed  Google Scholar 

  8. Vigano P, Somigliana E, Chiodo I, Abbiati A, Vercellini P: Molecular mechanisms and biological plausibility underlying the malignant transformation of endometriosis: a critical analysis. Hum Reprod Update. 2006, 12 (1): 77-89.

    Article  CAS  PubMed  Google Scholar 

  9. Munksgaard PS, Blaakaer J: The association between endometriosis and gynecological cancers and breast cancer: a review of epidemiological data. Gynecol Oncol. 2011, 123 (1): 157-163. 10.1016/j.ygyno.2011.06.017.

    Article  PubMed  Google Scholar 

  10. Pearce CL, Templeman C, Rossing MA, Lee A, Near AM, Webb PM, Nagle CM, Doherty JA, Cushing-Haugen KL, Wicklund KG, Chang-Claude J, Hein R, Lurie G, Wilkens LR, Carney ME, Goodman MT, Moysich K, Kjaer SK, Hogdall E, Jensen A, Goode EL, Fridley BL, Larson MC, Schildkraut JM, Palmieri RT, Cramer DW, Terry KL, Vitonis AF, Titus LJ, Ziogas A, et al: Association between endometriosis and risk of histological subtypes of ovarian cancer: a pooled analysis of case–control studies. Lancet Oncol. 2012, 13 (4): 385-394. 10.1016/S1470-2045(11)70404-1.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Anglesio MS, Carey MS, Kobel M, Mackay H, Huntsman DG: Clear cell carcinoma of the ovary: a report from the first ovarian clear cell symposium, june 24th, 2010. Gynecol Oncol. 2011, 121 (2): 407-415. 10.1016/j.ygyno.2011.01.005.

    Article  PubMed  Google Scholar 

  12. Wiegand KC, Shah SP, Al-Agha OM, Zhao Y, Tse K, Zeng T, Senz J, McConechy MK, Anglesio MS, Kalloger SE, Yang W, Heravi-Moussavi A, Giuliany R, Chow C, Fee J, Zayed A, Prentice L, Melnyk N, Turashvili G, Delaney AD, Madore J, Yip S, McPherson AW, Ha G, Bell L, Fereday S, Tam A, Galletta L, Tonin PN, Provencher D, et al: ARID1A mutations in endometriosis-associated ovarian carcinomas. N Engl J Med. 2010, 363 (16): 1532-1543. 10.1056/NEJMoa1008433.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Jones S, Wang TL, Shih I, Mao TL, Nakayama K, Roden R, Glas R, Slamon D, Diaz LA, Vogelstein B, Kinzler KW, Velculescu VE, Papadopoulos N: Frequent mutations of chromatin remodeling gene ARID1A in ovarian clear cell carcinoma. Sci (New York, NY). 2010, 330 (6001): 228-231. 10.1126/science.1196333.

    Article  CAS  Google Scholar 

  14. Guan B, Wang TL, Shih Ie M: ARID1A, a factor that promotes formation of SWI/SNF-mediated chromatin remodeling, is a tumor suppressor in gynecologic cancers. Cancer Res. 2011, 71 (21): 6718-6727. 10.1158/0008-5472.CAN-11-1562.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Lowery WJ, Schildkraut JM, Akushevich L, Bentley R, Marks JR, Huntsman D, Berchuck A: Loss of ARID1A-associated protein expression is a frequent event in clear cell and endometrioid ovarian cancers. Int J Gynecol Cancer. 2012, 22 (1): 9-14. 10.1097/IGC.0b013e318231f140.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Kuo KT, Mao TL, Jones S, Veras E, Ayhan A, Wang TL, Glas R, Slamon D, Velculescu VE, Kuman RJ, Shih IeM: Frequent activating mutations of PIK3CA in ovarian clear cell carcinoma. Am J Pathol. 2009, 174 (5): 1597-1601. 10.2353/ajpath.2009.081000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Campbell IG, Russell SE, Choong DY, Montgomery KG, Ciavarella ML, Hooi CS, Cristiano BE, Pearson RB, Phillips WA: Mutation of the PIK3CA gene in ovarian and breast cancer. Cancer Res. 2004, 64 (21): 7678-7681. 10.1158/0008-5472.CAN-04-2933.

    Article  CAS  PubMed  Google Scholar 

  18. Gonzalez-Angulo AM, Ferrer-Lozano J, Stemke-Hale K, Sahin A, Liu S, Barrera JA, Burgues O, Lluch AM, Chen H, Hortobagyi GN, Mills GB, Meric-Bernstam F: PI3K pathway mutations and PTEN levels in primary and metastatic breast cancer. Mol Cancer Ther. 2011, 10 (6): 1093-1101. 10.1158/1535-7163.MCT-10-1089.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. McConechy MK, Ding J, Cheang MC, Wiegand KC, Senz J, Tone AA, Yang W, Prentice LM, Tse K, Zeng T, McDonald H, Schmidt AP, Mutch DG, McAlpine JN, Hirst M, Shah SP, Lee CH, Goodfellow PJ, Gilks CB, Huntsman DG: Use of mutation profiles to refine the classification of endometrial carcinomas. J Pathol. 2012, 228 (1): 20-30.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Liang H, Cheung LW, Li J, Ju Z, Yu S, Stemke-Hale K, Dogruluk T, Lu Y, Liu X, Gu C, Guo W, Scherer SE, Carter H, Westin SN, Dyer MD, Verhaak RG, Zhang F, Karchin R, Liu CG, Lu KH, Broaddus RR, Scott KL, Hennessy BT, Mills GB: Whole-exome sequencing combined with functional genomics reveals novel candidate driver cancer genes in endometrial cancer. Genome Res. 2012, 22 (11): 2120-2129. 10.1101/gr.137596.112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Hennessy BT, Smith DL, Ram PT, Lu Y, Mills GB: Exploiting the PI3K/AKT pathway for cancer drug discovery. Nat Rev Drug Discov. 2005, 4 (12): 988-1004. 10.1038/nrd1902.

    Article  CAS  PubMed  Google Scholar 

  22. 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 cells. Mol Cancer Ther. 2006, 5 (10): 2512-2521. 10.1158/1535-7163.MCT-06-0334.

    Article  CAS  PubMed  Google Scholar 

  23. Anderson CL, Wang Y, Rustandi RR: Applications of imaged capillary isoelectric focussing technique in development of biopharmaceutical glycoprotein-based products. Electrophoresis. 2012, 33 (11): 1538-1544. 10.1002/elps.201100611.

    Article  CAS  PubMed  Google Scholar 

  24. 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 (15): 6084-6091. 10.1158/0008-5472.CAN-07-6854.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Ho CM, Lin MC, Huang SH, Huang CJ, Lai HC, Chien TY, Chang SF: PTEN promoter methylation and LOH of 10q22-23 locus in PTEN expression of ovarian clear cell adenocarcinomas. Gynecol Oncol. 2009, 112 (2): 307-313. 10.1016/j.ygyno.2008.09.040.

    Article  CAS  PubMed  Google Scholar 

  27. Iacovides DC, Johnson AB, Wang N, Boddapatti S, Korkola J, Gray JW: Identification and quantification of AKT isoforms and phosphoforms in breast cancer using a novel nanofluidic immunoassay. Mol Cell Proteomics. 2013, 12 (11): 3210-3220. 10.1074/mcp.M112.023119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Kamel D, Brady B, Tabchy A, Mills GB, Hennessy B: Proteomic classification of breast cancer. Curr Drug Targets. 2012, 13 (12): 1495-1509. 10.2174/138945012803530080.

    Article  CAS  PubMed  Google Scholar 

  29. Yang JY, Yoshihara K, Tanaka K, Hatae M, Masuzaki H, Itamochi H, Takano M, Ushijima K, Tanyi JL, Coukos G, Lu Y, Mills GB, Verhaak RG: Predicting time to ovarian carcinoma recurrence using protein markers. J Clin Invest. 2013, 123 (9): 3740-3750.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Carey MS, Agarwal R, Gilks B, Swenerton K, Kalloger S, Santos J, Ju Z, Lu Y, Zhang F, Coombes KR, Miller D, Huntsman D, Mills GB, Hennessy BT: Functional proteomic analysis of advanced serous ovarian cancer using reverse phase protein array: TGF-beta pathway signaling indicates response to primary chemotherapy. Clin Cancer Res. 2010, 16 (10): 2852-2860. 10.1158/1078-0432.CCR-09-2502.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Cheung LW, Hennessy BT, Li J, Yu S, Myers AP, Djordjevic B, Lu Y, Stemke-Hale K, Dyer MD, Zhang F, Ju Z, Cantley LC, Scherer SE, Liang H, Lu KH, Broaddus RR, Mills GB: High frequency of PIK3R1 and PIK3R2 mutations in endometrial cancer elucidates a novel mechanism for regulation of PTEN protein stability. Cancer Discov. 2011, 1 (2): 170-185. 10.1158/2159-8290.CD-11-0039.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Cizkova M, Susini A, Vacher S, Cizeron-Clairac G, Andrieu C, Driouch K, Fourme E, Lidereau R, Bieche I: PIK3CA mutation impact on survival in breast cancer patients and in ERalpha, PR and ERBB2-based subgroups. Breast Cancer Res. 2012, 14 (1): R28-10.1186/bcr3113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Vasudevan KM, Barbie DA, Davies MA, Rabinovsky R, McNear CJ, Kim JJ, Hennessy BT, Tseng H, Pochanard P, Kim SY, Dunn IF, Schinzel AC, Sandy P, Hoersch S, Sheng Q, Gupta PB, Boehm JS, Reiling JH, Silver S, Lu Y, Stemke-Hale K, Dutta B, Joy C, Sahin AA, Gonzalez-Angulo AM, Lluch A, Rameh LE, Jacks T, Root DE, Lander ES, Lluch A, et al: AKT-independent signaling downstream of oncogenic PIK3CA mutations in human cancer. Cancer Cell. 2009, 16 (1): 21-32. 10.1016/j.ccr.2009.04.012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Fadare O, Gwin K, Desouki MM, Crispens MA, Jones HW, Khabele D, Liang SX, Zheng W, Mohammed K, Hecht JL, Parkash V: The clinicopathologic significance of BAF250a (AIRD1A) expression in clear cell carcinomas of the endometrium. Mol Pathol. 2013, 26 (8): 1101-1110. 10.1038/modpathol.2013.35.

    Article  CAS  Google Scholar 

  35. Balbas-Martinez C, Rodriguez-Pinilla M, Casanova A, Dominguez O, Pisano DG, Gomez G, Lloreta J, Lorente JA, Malats N, Real FX: ARID1A alterations are associated with FGFR3-wild type, poor prognosis, urothelial bladder tumors. PLoS One. 2013, 8 (5): e62483-10.1371/journal.pone.0062483.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Agarwal R, Carey M, Hennessy B, Mills GB: PI3K pathway-directed therapeutic strategies in cancer. Curr Opin Investig Drugs. 2010, 11: 615-628.

    CAS  PubMed  Google Scholar 

  37. Mackay HJ, Brady MF, Oza AM, Reuss A, Pujade-Lauraine E, Swart AM, Siddiqui N, Colombo N, Bookman MA, Pfisterer J, du Bois A, Gynecologic Cancer InterGroup: Prognostic relevance of uncommon ovarian histology in women with stage III/IV epithelial ovarian cancer. Int J Gynecol Cancer. 2010, 20 (6): 945-952. 10.1111/IGC.0b013e3181dd0110.

    Article  PubMed  Google Scholar 

  38. Takakura S, Takano M, Takahashi F, Saito T, Aoki D, Inaba N, Noda K, Sugiyama T, Ochiai K, Japanese Gynecologic Oncology G: Randomized phase II trial of paclitaxel plus carboplatin therapy versus irinotecan plus cisplatin therapy as first-line chemotherapy for clear cell adenocarcinoma of the ovary: a JGOG study. Int J Gynecol Cancer. 2010, 20 (2): 240-247. 10.1111/IGC.0b013e3181cafb47.

    Article  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

The authors would like to thank Yang Li, and Jie Li, Sylvia Bamford, Sarah Maines-Bandiera, and Clara Salamanca for their valuable technical assistance, and the CCSG RPPA core facility at the MD Anderson Cancer Center for their assistance with the RPPA.

Funding

This work was supported by grants from the British Columbia (BC) Cancer Foundation, the Vancouver General Hospital (VGH)–University of British Columbia Hospital Foundation (to the OvCaRe ovarian cancer research team in Vancouver), the Carraressi Foundation and the Canadian Institute of Health Research operating grant 106654.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mark S Carey.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Concept and design of study, formulation of research questions: MSC, KCW, BTH, GBM, DGH. Development/application of research methodologies: BG, KS-H, MM, YW, ZJ. YL, BTH, DGH, GBM. Acquisition of data: YL, KS-H, KCW, YW, MSC. Analysis and interpretation of data: FZ, MSC, KCW, SL, SEK, ZJ, BTH. Writing, review and/or revision of the manuscript: MSC, KCW, BG, BTH, GBM, DGH. Study supervision: MSC, DGH, BTH, GBM. All authors approved the final version of the manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Wiegand, K.C., Hennessy, B.T., Leung, S. et al. A functional proteogenomic analysis of endometrioid and clear cell carcinomas using reverse phase protein array and mutation analysis: protein expression is histotype-specific and loss of ARID1A/BAF250a is associated with AKT phosphorylation. BMC Cancer 14, 120 (2014). https://doi.org/10.1186/1471-2407-14-120

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2407-14-120

Keywords