Figures
Abstract
Background
The PI3K/AKT pathway plays a pivotal role in breast cancer development and maintenance. PIK3CA, encoding the PI3K catalytic subunit, is the oncogene exhibiting a high frequency of gain-of-function mutations leading to PI3K/AKT pathway activation in breast cancer. PIK3CA mutations have been observed in 30% to 40% of ERα-positive breast tumors. However the physiopathological role of PIK3CA mutations in breast tumorigenesis remains largely unclear.
Methodology/Principal Findings
To identify relevant downstream target genes and signaling activated by aberrant PI3K/AKT pathway in breast tumors, we first analyzed gene expression with a pangenomic oligonucleotide microarray in a series of 43 ERα-positive tumors with and without PIK3CA mutations. Genes of interest were then investigated in 249 ERα-positive breast tumors by real-time quantitative RT-PCR. A robust collection of 19 genes was found to be differently expressed in PIK3CA-mutated tumors. PIK3CA mutations were associated with over-expression of several genes involved in the Wnt signaling pathway (WNT5A, TCF7L2, MSX2, TNFRSF11B), regulation of gene transcription (SEC14L2, MSX2, TFAP2B, NRIP3) and metal ion binding (CYP4Z1, CYP4Z2P, SLC40A1, LTF, LIMCH1).
Citation: Cizkova M, Cizeron-Clairac G, Vacher S, Susini A, Andrieu C, Lidereau R, et al. (2010) Gene Expression Profiling Reveals New Aspects of PIK3CA Mutation in ERalpha-Positive Breast Cancer: Major Implication of the Wnt Signaling Pathway. PLoS ONE 5(12): e15647. https://doi.org/10.1371/journal.pone.0015647
Editor: James McCubrey, East Carolina University, United States of America
Received: August 2, 2010; Accepted: November 19, 2010; Published: December 30, 2010
Copyright: © 2010 Cizkova et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by INCA (Institut National du CAncer), Programme Cancéropôles 2003-2007 Ile de France, and Comité Départemental des Hauts-de-Seine of the Ligue Nationale Contre le Cancer. This work is part of the Cartes d'Identité des Tumeurs (CIT) program (http://cit.ligue-cancer.net/index.php/en) launched and directed by the Ligue Nationale Contre le Cancer. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Deregulation of the phosphatidylinositol 3-kinase (PI3K) signaling pathway is frequent in human cancers. Activation of PI3K, which catalyzes inositol lipid phosphorylation to produce phosphatidylinositol-3,4,5-trisphosphate, is one of the most important downstream molecular events following tyrosine kinase receptor activation. Phosphatidylinositol-3,4,5-trisphosphate activates the serine/threonine kinase AKT, which in turn regulates several signaling pathways controlling cell survival, apoptosis, proliferation, motility, and adhesion [1]. PI3K is a heterodimeric enzyme composed of a p110α catalytic subunit encoded by the PIK3CA gene, and a p85 regulatory subunit encoded by the PIK3R1 gene [2].
Gain-of-function mutations in PIK3CA have recently been found in several malignancies, including breast cancer [1], [3], [4]. PIK3CA is frequently mutated at “hotspots” in exons 9 and 20, corresponding to the helical (E542K and E545K) and kinase (H1047R) domains, respectively. P110α carrying a hotspot mutation has oncogenic activity, transforming primary fibroblasts in culture, inducing anchorage-independent cell growth, and causing tumors in animals [5], [6].
After the TP53 suppressor gene, the PIK3CA oncogene is the most frequently mutated gene in human breast cancers (up to 40% of breast tumors) [7], [8]. Activating somatic mutations of other oncogenes (EGFR, KRAS, HRAS, NRAF, BRAF and AKT1) involved in downstream molecular events following tyrosine kinase receptor activation are frequent in several malignancies but rare in breast cancer. Several studies suggest that PIK3CA mutations are more frequent in estrogen receptor alpha (ERα)-positive breast tumors (30–40%) than in ERα-negative breast tumors (10–20%) [7].
The pathological role of these gain-of-function PIK3CA mutations in breast tumors, and particularly in ERα-positive breast tumors, is largely unknown. Better knowledge of PIK3CA mutation impact requires the identification of downstream target genes and signaling pathways activated by aberrant PI3K/AKT signaling. Here, we compared gene expression in PIK3CA-mutated and PIK3CA wild-type ERα-positive breast tumors, using a genome-wide microarray and subsequently real-time quantitative reverse transcriptase-polymerase chain reaction (RT-PCR).
Materials and Methods
Patients and Samples
We analyzed samples of 292 primary unilateral non metastatic ERα-positive postmenopausal breast tumors excised from women at René Huguenin Hospital (Saint-Cloud, France) from 1978 to 2008. Other characteristics of the patients are listed in Table S1. Each patient gave written informed consent and this study was approved by the Local Ethical Committee (Breast Group of René Huguenin Hospital). Immediately after surgery the tumor samples were stored in liquid nitrogen until RNA extraction. The samples analyzed contained more than 70% of tumor cells. ERα status was determined at the protein level by using biochemical methods (Dextran-coated charcoal method until 1988 and enzyme immunoassay thereafter) and was confirmed at mRNA level by real-time RT-PCR. Forty-three samples were used as a microarray and RT-PCR screening set to identify differentially expressed genes. These genes were then validated in the remaining 249 ERα-positive tumors by means of RT-PCR. Control samples consisted of eight specimens of normal breast tissue collected from women undergoing cosmetic breast surgery or adjacent normal breast tissue from breast cancer patients.
RNA extraction
Total RNA was extracted from breast tissue by using the acid-phenol guanidium method, and its quality was determined by agarose gel electrophoresis and ethidium bromide staining. The 18S and 28S RNA bands were visualized under ultraviolet light.
PIK3CA mutation screening
PIK3CA mutation screening was performed on cDNA fragments obtained by RT-PCR amplification of exons 9 and 20 and their flanking exons. Details of the primers and PCR conditions are available on request. The amplified products were sequenced with the BigDye Terminator kit on an ABI Prism 3130 automatic DNA sequencer (Applied Biosystems, Courtabœuf, France). Sequences thus obtained were compared with the corresponding cDNA reference sequence (NM_006218).
Microarray analysis
Microarray experiments used Human Genome U133 Plus 2.0 arrays from Affymetrix, containing 54675 probe sets. Gene chips were hybridized and scanned using standard Affymetrix protocols. Expression data were obtained as CEL files. BRB ArrayTools (version 3.6.0 available on http://linus.nci.nih.gov/BRB-ArrayTools.html) were used to import CEL files with Robust Method Average (RMA) normalization, and to analyze gene expression. A class comparison based on a univariate t test applied to log-normalized data was used to identify genes differentially expressed in breast tumors with and without PIK3CA mutations. Supervised class prediction analysis was implemented with the Prediction Analysis for Microarrays (PAM) algorithm to identify genes required for optimal prediction [9].
The Database for Annotation, Visualization and Integrated Discovery (DAVID, available on http://david.abcc.ncifcrf.gov/) was used to interpret the lists of differentially expressed probes and to identify statistically overrepresented biological function categories of Gene Ontology (GO) and biological pathways, as defined in the Kyoto Encyclopedia of Genes and Genomes (KEGG).
In compliance with the Minimun Information About a Microarray Experiment (MIAME) recommendations, raw data were deposited in the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) under series accession number GSE22035.
Real-time quantitative RT-PCR
RT-PCR was applied to the selected genes, as well as ERα (NM_000125), MKI67 (NM_002417), and TBP (NM_003194; endogenous mRNA control). Primers and PCR conditions are available on request, and the RT-PCR protocol using the SYBR Green Master Mix kit on the ABI Prism 7900 Sequence Detection System (Perkin-Elmer Applied Biosystems, Foster City, CA, USA) is described in detail elsewhere [10]. The relative mRNA expression level of each gene, expressed as the N-fold difference in target gene expression relative to the TBP gene, and termed “Ntarget”, was calculated as Ntarget = 2ΔCtsample. The value of the cycle threshold (ΔCt) of a given sample was determined by subtracting the average Ct value of the target gene from the average Ct value of the TBP gene. The Ntarget values of the samples were subsequently normalized such that the median Ntarget value of the normal breast samples was 1. The relative expression of each gene was characterized by the median and range, and the differences in gene expression between tumors with and without PIK3CA mutations were analysed for significance with the non parametric Mann-Whitney U test.
Clustering analysis
Hierarchical clustering analyses of gene expression and samples were performed using BRB ArrayTools. Classification performance was calculated as overall accuracy, defined as the proportion of correctly classified tumors in each cluster, using Matthews' correlation coefficient (MCC) [11]. This parameter was used to discriminate identical accuracies. The chi-square test was used to determine the statistical significance of the clustering.
Results
Analysis of differentially expressed genes in 43 ERα-positive tumors
Overview of transcriptome changes in PIK3CA-mutated tumors.
To identify PIK3CA mutation-related genes, microarray analysis (Affymetrix U133 Plus 2.0 arrays) was first applied to 43 ERα-positive breast tumors, of which 14 were PIK3CA-mutated and 29 were wild-type (Table S1). We found that 6124 probes were differentially expressed between breast tumors with and without PIK3CA mutations, with P values <0.05. Of these, 2538 probes (1630 unique genes) were up-regulated (Table S2) and 3586 (2672 unique genes) were down-regulated (Table S3). Only 216 up-regulated probes (153 unique genes) and 28 down-regulated probes (18 unique genes) showed at least a 2-fold change (FC).
Gene ontology analysis of differentially expressed genes.
To identify families of genes that might have significant roles related to specific biological or molecular processes, we used the DAVID database to annotate the 6124 probes and categorize them by function. As shown in Table 1, these genes were mainly involved in the regulation of transcription, cell cycling, proliferation, death, adhesion and cytoskeleton organization, and also ion binding and transport, and ATP and RNA binding activity.
The 2672 down-regulated genes were mainly associated with ATP binding, acetylation and ion transport (Table 1). Among the down-regulated genes with FC≥2, no significantly overrepresented GO categories appeared.
Most of the 1630 up-regulated genes were involved in transcriptional regulation (17.3%) (biological process) and ion binding (25.6%) (molecular function) (Table 1). The latter included the metal ion-binding and zinc ion-binding categories (Table 1). As shown in Figure 1A, the 216 probes most strongly up-regulated in PIK3CA-mutated tumors (153 unique genes) belonged mainly to the ion-binding category (35.5%) but also to categories of structural molecule activity (including structural cytoskeleton constituents) (9.3%), transcription regulatory activity (9.3%) and nucleotide binding (including ATP and GTP binding) (7.5%).
Molecular functions were attributed to 107 of the 153 genes using GOTERM_MF_FAT from the DAVID database. Categories with at least three genes are represented in A. Subclassification of the 36 genes belonging to the metal ion-binding category is shown in B. All categories were represented and several genes were common to more than one category. Genes belonging to the metal ion-binding and transcription activity categories are listed in Table 2.
In the ion-binding category, the genes corresponded to genes encoding metal ion-binding proteins in 95% of cases: 28% encoding iron ion-binding and 23% with zinc ion-binding proteins (Figure 1B), pointing to a role of ion-binding proteins, and especially iron ion-binding proteins, in breast cancer with PIK3CA mutations. Interestingly, the genes belonging to the metal ion-binding category (Table 2) included two families of genes that were among the most strongly up-regulated in PIK3CA-mutated tumors. They comprised four genes of cytochrome P450 family 4 (CYP4Z1, CYP4X1, CYP4B1 and the pseudogene (CYP4Z2P) and two solute carrier genes (SLC4A4 and SLC40A1). All these genes, with exception of SLC4A4, are associated with iron ion binding. In addition to these genes, we found on the top of the list lactoferrin (LTF), also known to be involved in iron metabolism. Among the genes encoding zinc ion-binding proteins, three (ANPEP, LIMCH1 and NR2F2) are known to be cancer-related.
Besides NR2F2, five other transcription factors, all known to be involved in tumorigenesis, were identified (Table 2): (a) TFAP2B, a tumor suppressor gene in breast cancer [12], (b) SEC14L2, a gene possibly involved in the antiproliferative effect of vitamin E in cancer [13], (c) ID4, which has been proposed to be involved in breast cancer, inhibiting mammary epithelial cell differentiation and stimulating mammary epithelial cell growth [14], (d) TCF7L2, also named TCF4, a cancer-promoting gene involved in the Wnt signaling pathway [15], and (e) MSX2, a gene implicated in mammary gland and breast cancer development [16], and which is also activated by Wnt signaling [17].
These five transcriptional factors (TFAP2B, SEC14L2, ID4, TCF7L2 and MSX2), as well as ten genes involved in metal ion binding (CYP4Z1, CYP4X1, CYP4B1, CYP4Z2P, SLC4A4, SLC40A1, LTF, ANPEP, LIMCH1 and NR2F2), were selected for validation by RT-PCR.
Pathway analysis of differentially expressed genes.
By applying KEGG pathway analysis to the 6124 probes differentially expressed in PIK3CA-mutated tumors, we identified physiological pathways directly or indirectly associated with PIK3CA mutations. The most significantly overrepresented pathways are shown in Table 1. In addition to signaling pathways in cancer cells, the following five signaling networks were thus identified: MAPK, Calcium, Jak-STAT, Wnt and apoptosis. The Calcium signaling pathway was specifically altered by the down-regulated genes, whereas the Wnt signaling pathway was specifically altered by the up-regulated genes. The same method applied to the 216 probes (153 unique genes) that were up-regulated with FC≥2 also revealed the Wnt signaling pathway (P = 0.015) (data not shown), highlighting the importance of this pathway in PIK3CA-mutated tumors. Five major genes of the Wnt signaling were thus recognized among the 216 probes (Table S2): MSX2 and TCF7L2 (already cited), and WNT5A, VANGL2 and TNFRSF11B/osteoprotegerin. These genes were also selected for RT-PCR validation.
Finally, among the genes up-regulated with FC≥2 (Table S2), we identified PIK3R1, the gene encoding the PI3K regulatory subunit, and two other genes of interest: HMGCS2, a nuclear gene encoding a mitochondrial matrix enzyme involved in ketogenesis and cholesterol synthesis, processes possibly implicated in the etiology or progression of breast cancer [18] and MAPT, a protein enhancing microtubule assembly and stability, that might be involved in taxane resistance [19]. These three genes were added to the RT-PCR validation set.
Two-class prediction analysis of differentially expressed genes.
Two-class prediction analysis with the PAM algorithm was used to identify the group of genes that best characterized PIK3CA-mutated and wild-type tumors and that classified the tumors with the smallest number of predictive features. A threshold of 2.81, that minimized the error, identified 56 differentially expressed probes corresponding to 39 unique genes (Table S4). Thirty-eight of these 39 unique genes were over-expressed in ERα-positive breast tumors with PIK3CA mutations, 16 being up-regulated at least 3-fold, while only one gene (NKAIN1, encoding Na+/K+ ATPase interacting protein) was down-regulated, with a FC of 3.52. Interestingly, two major genes involved in the Wnt signaling pathway were also identified by PAM, namely WNT5A (the most predictive gene) and TCF7L2, further confirming the importance of this pathway in PIK3CA-mutated tumors. The previously selected up-regulated genes were almost all included in the list of the most predictive genes. PAM analysis identified five interesting new genes that were up-regulated with FC≥3, namely VTCN1, TMC5, NTN4, REEP1 and NRIP3, which were added to the RT-PCR validation set.
Among the down-regulated genes, NKAIN1 was selected for RT-PCR validation, along with two other genes known to be involved in cancer biology: TUSC3 and TPD52, that were among the 28 most strongly down-regulated probes (FC≥2) (Table S3) and that were also among the most predictive genes in PAM analysis with a lower FC threshold of 2.5 (data not shown).
Combined analysis of the GO, KEGG and PAM approaches identified 29 most promising genes (26 up-regulated and 3 down-regulated) for RT-PCR validation. The expression status of these genes was first confirmed in the same series of 43 breast tumors (Table 3). Strong positive correlations were observed between the microarray and RT-PCR expression levels of each gene (Spearman's correlation coefficients ranged from 0.69 to 0.97 and were all significant, at P<0.0001; data not shown).
mRNA expression of the 29 genes of interest in 249 ERα-positive breast tumors
Overall expression of the 29 differentially expressed genes.
The expression levels of the 29 genes selected by microarray analysis were then verified by RT-PCR in a large independent cohort of 249 ERα-positive breast tumors, of which 157 were PIK3CA wild-type and 92 were PIK3CA-mutated (Table S1). This PIK3CA mutation frequency of 37% was in keeping with the results of previous studies showing a mutation rate of up to 40% in ERα-positive breast tumors [7], [8]. Almost all the tumors had a single mutation, 44 (47.8%) in exon 9 (helical domain) and 46 (50%) in exon 20 (kinase domain) [7]. Two tumors (2.2%) carried two mutations, located in exons 9 and 20 in one case, and in exon 20 in the second case.
Nineteen (66%) of the 29 selected genes showed significantly different expression between mutated and wild-type tumors in the validation cohort, with a distribution similar to that observed in the screening cohort (Table 4). Among the three down-regulated genes of interest in the screening set, only NKAIN1 was significantly down-regulated in the validation set. Among the 26 up-regulated genes, 18 were also up-regulated in the validation set. With exception of VANGL2, up-regulation of the genes involved in Wnt signaling pathway, namely WNT5A, MSX2, TCF7L2 and TNFRSF11B, was confirmed in the validation set, further emphasizing the important role of the Wnt signaling pathway in PIK3CA-mutated breast cancer. Up-regulation was also confirmed for genes related to breast cancer (MAPT, HMGCS2, NR2F2 and TFAP2B), genes involved in metal ion binding (CYP4Z1, CYP4Z2P, SLC40A1, LTF and LIMCH1) and also NRIP3, NTN4, REEP1, SEC14L2 and TMC5. Deregulation of these genes was not related to ERα status or proliferation since similar expression levels of ERα and MKI67 were observed in PIK3CA-mutated and -non mutated tumors (Table 4). Only 2 of the 29 selected genes showed significantly different expression between PIK3CA exon 9- and exon 20-mutated tumors, namely TFAP2B and NRIP3 (Table S5). Interestingly, TFAP2B was over-expressed in exon 20-mutated tumors and NRIP3 in exon 9-mutated tumors.
Identification of the most discriminatory genes.
PAM prediction analysis was then used to test the ability of each gene to classify the 249 ERα-positive breast tumors according to PIK3CA mutation status. NKAIN1 was the most predictive gene (PAM rank) (Table 4). NKAIN1 was also an essential classifier in supervised hierarchical clustering analysis. Indeed, the 19-gene set including NKAIN1 classified the 249 breast tumors significantly more accurately than the set of 18 up-regulated genes without NKAIN1 (accuracy 59% and 57%, Χ2 test P values of 0.0006 and 0.0141, respectively) (Table S6). Three different minimal sets of 4, 5 and 6 genes, all including both NKAIN1 and CYP4Z2P, showed the same overall clustering accuracy of 59.4% (Table S6). However, the 5-gene group (NKAIN1-CYP4Z2P-NRIP3-SEC14L2-TFAP2B) had the most significant discriminatory value (MCC = 0.2334, P = 0.0002) correctly clustering 66 of the 92 mutated tumors and 82 of the 157 non mutated tumors. Notably, this 5-gene set contained the two genes that were differently expressed between exon 9- and exon 20-mutated tumors, and thus had the best capacity to distinguish between these two tumor categories (data not shown). The other two gene sets, both comprising genes involved in Wnt signaling (NKAIN1-CYP4Z2P-WNT5A-TMC5 and NKAIN1-CYP4Z2P-WNT5A-MAPT-MSX2-TFAP2B), classified 65 mutated and 83 non mutated tumors correctly (MCC = 0.2286, P = 0.0003).
Discussion
We used a two-step strategy to identify downstream target genes and signaling pathways affected by PIK3CA mutations in breast tumors. We first applied a pangenomic oligonucleotide microarray approach to a series of 43 ERα-positive tumors with and without PIK3CA mutations, and then validated genes of interest by RT-PCR in an independent series of 249 ERα-positive tumors. A robust set of 19 genes differentially expressed in PIK3CA-mutated and wild-type tumors was thus identified.
Over-expression of several genes involved in Wnt signaling (WNT5A, TCF7L2, MSX2 and TNFRSF11B), regulation of gene transcription (SEC14L2, MSX2, TFAP2B and NRIP3) and metal ion binding (CYP4Z1, CYP4Z2P, SLC40A1, LTF and LIMCH1) was observed in PIK3CA-mutated tumors. Several of these genes have been linked to breast cancer (MAPT, HMGCS2, NR2F2, TFAP2B, NTN4, SEC14L2 and LTF).
The human Wnt signaling network is important for regulation of proliferation, differentiation, growth and survival from the embryo stage [20], [21]. Crosstalk of complex pathways belonging to Wnt signaling has been observed leading to, when altered, disparate effects in different tumor types [22]-[24]. We observed over-expression of four major genes involved in the Wnt pathway, namely WNT5A, TCF7L2, MSX2 and TNFRSF11B. WNT5A encodes a major Wnt ligand affecting tumor cell motility and metastasis, but its role in breast cancer is controversial [23]. The emerging view is that, in breast cancer, WNT5A has a suppressive effect, inhibiting migration and invasion of breast cancer cell lines [24]. Moreover, WNT5A over-expression observed in invasive breast tumors has been associated with a favorable outcome [24]. PIK3CA mutations have also been associated with favorable outcome of breast cancer patients [1], [3], [4], [25]. We can thus suggest a link between gain-of-function mutation in PIK3CA, up-regulation of WNT5A and favorable outcome in breast cancer. We also observed over-expression of TCF7L2, which encodes one of the four major transcription factors involved in the Wnt signaling pathway [20], [26], as well as two other genes (MSX2 and TNFRSF11B) known to be downstream targets of the Wnt signaling pathway [27]–[29]. Wnt signaling has a major role in cancer stem cell self-renewal and tumor maintenance [20], [30] and contributes to tumor invasion, metastasis and angiogenesis [31]. Recent studies have identified a role of Wnt pathway in epidermal-mesenchymal transition during breast cancer development [32], [33]. Thus, Wnt pathway activation appears to be an important consequence of PIK3CA mutations in breast tumors, in keeping with recently observed crosstalk between the PI3K/Akt and Wnt pathways in both physiological (myeloid progenitor cells) [21] and pathological conditions (medulloblastoma) [34].
Better understanding of the biological functions of the Wnt and PI3K/Akt pathways and their interplay could have therapeutic implications for breast cancer. Drugs targeting the PI3K/Akt pathway have given promising preliminary results in human malignancies [35], [36]. However, as the PI3K pathway is crucial for metabolic processes, PI3K inhibitors might also have side effects, especially affecting insulin signaling and cardiac functions [36], [37]. In contrast, targeting of downstream Wnt signaling events might have fewer adverse effects, considering their crucial importance in embryonic development [23], [38].
Genes encoding metal ion-binding proteins were also over-expressed in PIK3CA-mutated tumors. Such metal ion-binding proteins have regulatory roles in central cellular processes such as gene expression, proliferation, differentiation and survival. Increased expression of these proteins in ERα-positive breast tumors has also been reported by Abba et al. [39]. We observed over-expression of LIMCH1, a gene encoding zinc-binding protein, and also four genes encoding iron-binding proteins (LTF, SLC40A1, CYP4Z1 and CYP4Z2P) previously linked to breast cancer. LTF encodes lactoferrin, a protein involved in non specific immunity and that may inhibit carcinogenesis and tumor growth [40]. CYP4Z1 and its pseudogene CYP4Z2P are two members of cytochrome P450 family 4 which have been found to be over-expressed in about 50% of breast cancers relative to normal breast tissue from the same patients [41]. Here, we confirm that the pseudogene CYP4Z2P is expressed in both PIK3CA-mutated and -non mutated ERα-positive breast tumors, by using specific primers unambiguously distinguishing CYP4Z2P from CYP4Z1. Thus, CYP4Z2P is transcriptionally active, but its translation remains to be studied. CYP4Z2P is located in a head-to-head orientation close to CYP4Z1 in chromosome region 1p33 [41], raising the possibility that expression of these two genes is co-regulated in PIK3CA-mutated breast tumors.
We identified several genes previously implicated in breast cancer development or outcome. The proteins encoded by TFAP2B, NTN4 and SEC14L2 have been linked to tumors with less aggressive features and better outcome [12], [13], [42]. MAPT has been proposed as a predictive marker of taxane responsiveness in breast cancer [43]. NR2F2 has been also detected up-regulated in breast cancer, but its involvement in tumor development remains elusive because of its ability to affect both pro-oncogenic and anti-oncogenic proteins [44], [45]. HMGCS2 was recently shown to be regulated in response to hormonal stimulation [18].
NRIP3, TMC5, REEP1 and NKAIN1, whose expression had not previously been described in breast cancer, were also deregulated in the PIK3CA-mutated breast tumors. NRIP3, TMC5 and REEP1 are differentially expressed in various other tumor types [46]–[48]. Interestingly, NKAIN1 was the only gene under-expressed in PIK3CA-mutated tumors and was also the most discriminatory gene for these tumors. The role of these genes in breast cancer development remains to be evaluated in following studies.
Recently, Loi et al. identified a 278 gene-expression signature associated specifically with PIK3CA exon 20-mutated ER-positive/ERBB2-negative tumors [25]. These authors observed an unexpected significant down-expression of some Akt-regulated genes such as RPS6KB1 in their PIK3CA-mutated tumor series, but a normal level of AKT1 and mTOR transcripts. They also showed that phosphor-Akt expression was not significantly up-regulated at the protein level. In the present study, we did not identify RPS6KB1, AKT1 and mTOR in our final 19-gene set nor in the list of 6124 genes differentially expressed in PIK3CA-mutated tumors (Table S2 and S3). Interestingly, among the 168 significantly up-regulated genes detected by Loi et al., WNT5A and MSX2, as well as HMGCS2 and LTF, were identified in agreement with our results. The data of Loi et al. [25] confirm thus the positive association between PIK3CA mutation and Wnt signaling pathway activation reported in the present manuscript.
In conclusion, this gene expression profiling study suggests that over-expression of genes belonging to the Wnt signaling pathway may play a pivotal role in PIK3CA-mutated breast tumors, in particular WNT5A. Further studies of biological mechanisms affected by PIK3CA mutations may have therapeutic implications.
Supporting Information
Table S1.
Molecular, pathological and clinical characteristics of patients in relation to metastasis free survival (MFS) in the 43 ERα-positive and 249 ERα-positive patient series.
https://doi.org/10.1371/journal.pone.0015647.s001
(PDF)
Table S2.
2538 probes up-regulated in tumors with PIK3CA mutations (Mutated) compared to tumors without PIK3CA mutation (Normal) with a P value <0.05 identified by a parametric t test using BRB ArrayTools. These genes were ranked according to fold change (FC) calculated between expression intensities of tumors with PIK3CA mutations and those of tumors without PIK3CA mutation. The 216 probes with a FC≥2 are put in bold. In this list, the probes belonging to Wnt signaling pathway are shaded in light grey, and PIK3R1, HMGCS2 and MAPT are shaded in dark grey.
https://doi.org/10.1371/journal.pone.0015647.s002
(PDF)
Table S3.
3586 probes down-regulated in tumors with PIK3CA mutations (Mutated) compared to tumors without PIK3CA mutation (Normal) with a P value <0.05 identified by a parametric t test using BRB ArrayTools. These genes were ranked according to fold change (FC) calculated between expression intensities of tumors with PIK3CA mutations and those of tumors without PIK3CA mutation. The 28 probes with a FC≥2 are put in bold.
https://doi.org/10.1371/journal.pone.0015647.s003
(PDF)
Table S4.
List of 56 probes (39 unique genes) deregulated in REα-positive breast tumors with PIK3CA mutations compared to those without PIK3CA mutation identified by PAM. These genes are presented according to the rank in PAM output. For each gene, we report the fold change (FC) calculated between expression intensities of tumors with and without PIK3CA mutations using BRB Arrays Tools. The genes with a FC≥3 are put in bold.
https://doi.org/10.1371/journal.pone.0015647.s004
(PDF)
Table S5.
Relative mRNA expression levels of the 29 genes in the 44 ERα-positive breast tumors with exon 9 PIK3CA mutations compared to the 47 tumors with exon 20 PIK3CA mutations. The tumor with PIK3CA mutations in both exon 9 and exon 20 was excluded from the analysis. For each gene, we report the median (range) of the mRNA levels of each gene relative to normal breast tissue samples, the fold change (FC) between tumors with exon 9-mutated and exon 20-mutated PIK3CA and the P value associated to Mann-Whitney U test.
https://doi.org/10.1371/journal.pone.0015647.s005
(PDF)
Table S6.
Supervised hierarchical clustering analysis of the 249 ERα-positive breast tumors. Classification performance of discriminating gene sets identified from the 19 significantly deregulated genes in tumors with PIK3CA mutations (18 up-regulated genes + NKAIN1). Each gene set separates the 157 tumors without PIK3CA mutation (N) and the 92 tumors with PIK3CA mutations (M) in two main clusters (cluster 1 and 2).
https://doi.org/10.1371/journal.pone.0015647.s006
(PDF)
Acknowledgments
We thank Aurélien de Reyniès and Laetitia Marisa of the Ligue Nationale Contre le Cancer. We are extremely thankful for the contributions made by the CIT platforms (Affymetrix, IGBMC: Christelle Thibault, Philippe Kastner; BAC arrays, Institut Curie: Olivier Delattre; RNA analysis, Saint-Louis Hospital: Mira Ayadi).
Author Contributions
Conceived and designed the experiments: IB RL. Performed the experiments: SV AS CA. Analyzed the data: GCC. Contributed reagents/materials/analysis tools: GCC IB RL. Wrote the paper: MC GCC IB MS GCC.
References
- 1. Zhao L, Vogt PK (2008) Class I PI3K in oncogenic cellular transformation. Oncogene 27: 5486–5496.
- 2. Dillon RL, White DE, Muller WJ (2007) The phosphatidyl inositol 3-kinase signaling network: implications for human breast cancer. Oncogene 26: 1338–1345.
- 3. Samuels Y, Wang Z, Bardelli A, Silliman N, Ptak J, et al. (2004) High frequency of mutations of the PIK3CA gene in human cancers. Science 304: 554.
- 4. Karakas B, Bachman KE, Park BH (2006) Mutation of the PIK3CA oncogene in human cancers. Br J Cancer 94: 455–459.
- 5. Zhao JJ, Liu Z, Wang L, Shin E, Loda MF, et al. (2005) The oncogenic properties of mutant p110alpha and p110beta phosphatidylinositol 3-kinases in human mammary epithelial cells. Proc Natl Acad Sci U S A 102: 18443–18448.
- 6. Bader AG, Kang S, Vogt PK (2006) Cancer-specific mutations in PIK3CA are oncogenic in vivo. Proc Natl Acad Sci U S A 103: 1475–1479.
- 7. Saal LH, Holm K, Maurer M, Memeo L, Su T, et al. (2005) PIK3CA mutations correlate with hormone receptors, node metastasis, and ERBB2, and are mutually exclusive with PTEN loss in human breast carcinoma. Cancer Res 65: 2554–2559.
- 8. Stemke-Hale K, Gonzalez-Angulo AM, Lluch A, Neve RM, Kuo WL, et al. (2008) An integrative genomic and proteomic analysis of PIK3CA, PTEN, and AKT mutations in breast cancer. Cancer Res 68: 6084–6091.
- 9. Tibshirani R, Hastie T, Narasimhan B, Chu G (2002) Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci U S A 99: 6567–6572.
- 10. Bieche I, Onody P, Laurendeau I, Olivi M, Vidaud D, et al. (1999) Real-time reverse transcription-PCR assay for future management of ERBB2-based clinical applications. Clin Chem 45: 1148–1156.
- 11. Matthews BW (1975) Comparison of the predicted and observed secondary structure of T4 phage lysozyme. Biochim Biophys Acta 405: 442–451.
- 12. Gee JM, Robertson JF, Ellis IO, Nicholson RI, Hurst HC (1999) Immunohistochemical analysis reveals a tumour suppressor-like role for the transcription factor AP-2 in invasive breast cancer. J Pathol 189: 514–520.
- 13. Johnykutty S, Tang P, Zhao H, Hicks DG, Yeh S, et al. (2009) Dual expression of alpha-tocopherol-associated protein and estrogen receptor in normal/benign human breast luminal cells and the downregulation of alpha-tocopherol-associated protein in estrogen-receptor-positive breast carcinomas. Mod Pathol 22: 770–775.
- 14. Shan L, Yu M, Qiu C, Snyderwine EG (2003) Id4 regulates mammary epithelial cell growth and differentiation and is overexpressed in rat mammary gland carcinomas. Am J Pathol 163: 2495–2502.
- 15. Tang W, Dodge M, Gundapaneni D, Michnoff C, Roth M, et al. (2008) A genome-wide RNAi screen for Wnt/beta-catenin pathway components identifies unexpected roles for TCF transcription factors in cancer. Proc Natl Acad Sci U S A 105: 9697–9702.
- 16. Satoh K, Ginsburg E, Vonderhaar BK (2004) Msx-1 and Msx-2 in mammary gland development. J Mammary Gland Biol Neoplasia 9: 195–205.
- 17. Cheng SL, Shao JS, Cai J, Sierra OL, Towler DA (2008) Msx2 exerts bone anabolism via canonical Wnt signaling. J Biol Chem 283: 20505–20522.
- 18. Riggins RB, Lan JP-J, Klimach U, Zwart A, Cavalli LR, et al. (2008) ERR Mediates Tamoxifen Resistance in Novel Models of Invasive Lobular Breast Cancer. Cancer Res 68: 8908–8917.
- 19. Andre F, Mazouni C, Liedtke C, Kau SW, Frye D, et al. (2008) HER2 expression and efficacy of preoperative paclitaxel/FAC chemotherapy in breast cancer. Breast Cancer Res Treat 108: 183–190.
- 20. Reya T, Clevers H (2005) Wnt signalling in stem cells and cancer. Nature 434: 843–850.
- 21. Nteliopoulos G, Marley SB, Gordon MY (2009) Influence of PI-3K/Akt pathway on Wnt signalling in regulating myeloid progenitor cell proliferation. Evidence for a role of autocrine/paracrine Wnt regulation. Br J Haematol 146: 637–651.
- 22. Veeman MT, Axelrod JD, Moon RT (2003) A second canon. Functions and mechanisms of beta-catenin-independent Wnt pathway. Dev Cell 5: 367–77.
- 23. Leris AC, Roberts TR, Jiang WG, Newbold RF, Mokbel K (2005) WNT5A expression in human breast cancer. Anticancer Res 25: 731–734.
- 24. Ford CE, Ekström EJ, Howlin J, Andersson T (2009) The WNT-5a derived peptide, Foxy-5, possesses dual properties that impair progression of ERalpha negative breast cancer. Cell Cycle 8: 1838–42.
- 25. Loi S, Haibe-Kains B, Majjaj S, Lallemand F, Durbecq V, et al. (2010) Proc Natl Acad Sci U S A 107: 10208–13.
- 26. Katoh M (2005) WNT/PCP signaling pathway and human cancer. Oncol Rep 14: 1583–1588.
- 27. Willert J, Epping M, Pollack JR, Brown PO, Nusse R (2002) A transcriptional response to Wnt protein in human embryonic carcinoma cells. BMC Dev Biol 2: 8.
- 28. Krishnan V, Bryant HU, Macdougald OA (2006) Regulation of bone mass by Wnt signaling. J Clin Invest 116: 1202–1209.
- 29. Bu G, Lu W, Liu CC, Selander K, Yoneda T, et al. (2008) Breast cancer-derived Dickkopf1 inhibits osteoblast differentiation and osteoprotegerin expression: implication for breast cancer osteolytic bone metastases. Int J Cancer 123: 1034–1042.
- 30. Malanchi I, Huelsken J (2009) Cancer stem cells: never Wnt away from the niche. Curr Opin Oncol 21: 41–46.
- 31. Wang Y (2009) Wnt/Planar cell polarity signaling: A new paradigm for cancer therapy. Mol Cancer Ther 8: 2103–2109.
- 32. Yook JI, Li XY, Ota I, Hu C, Kim HS, et al. (2006) A Wnt-Axin2-GSK3beta cascade regulates Snail1 activity in breast cancer cells. Nat Cell Biol 8: 1398–1406.
- 33. DiMeo TA, Anderson K, Phadke P, Fan C, Perou CM, et al. (2009) A novel lung metastasis signature links Wnt signaling with cancer cell self-renewal and epithelial-mesenchymal transition in basal-like breast cancer. Cancer Res 69: 5364–5373.
- 34. Baryawno N, Sveinbjörnsson B, Eksborg S, Chen CS, Kogner P, et al. (2010) Small-molecule inhibitors of phosphatidylinositol 3-kinase/Akt signaling inhibit Wnt/beta-catenin pathway cross-talk and suppress medulloblastoma growth. Cancer Res 70: 266–276.
- 35. Samuels Y, Diaz LA Jr, Schmidt-Kittler O, Cummins JM, Delong L, et al. (2005) Mutant PIK3CA promotes cell growth and invasion of human cancer cells. Cancer Cell 7: 561–573.
- 36. Carnero A (2009) Novel inhibitors of the PI3K family. Expert Opin Investig Drugs 18: 1265–1277.
- 37. Weng LP, Smith WM, Brown JL, Eng C (2001) PTEN inhibits insulin-stimulated MEK/MAPK activation and cell growth by blocking IRS-1 phosphorylation and IRS-1/Grb-2/Sos complex formation in a breast cancer model. Hum Mol Genet 10: 605–616.
- 38. Barker N, Clevers H (2006) Mining the Wnt pathway for cancer therapeutics. Nat Rev Drug Discov 5: 997–1014.
- 39. Abba MC, Hu Y, Sun H, Drake JA, Gaddis S, et al. (2005) Gene expression signature of estrogen receptor alpha status in breast cancer. BMC Genomics 6: 37.
- 40. Benaïssa M, Peyrat JP, Hornez L, Mariller C, Mazurier J, et al. (2005) Expression and prognostic value of lactoferrin mRNA isoforms in human breast cancer. Int J Cancer 114: 299–306.
- 41. Rieger MA, Ebner R, Bell DR, Kiessling A, Rohayem J, et al. (2004) Identification of a novel mammary-restricted cytochrome P450, CYP4Z1, with overexpression in breast carcinoma. Cancer Res 64: 2357–2364.
- 42. Esseghir S, Kennedy A, Seedhar P, Nerurkar A, Poulsom R, et al. (2007) Identification of NTN4, TRA1, and STC2 as prognostic markers in breast cancer in a screen for signal sequence encoding proteins. Clin Cancer Res 13: 3164–3173.
- 43. Rouzier R, Rajan R, Wagner P, Hess KR, Gold DL, et al. (2005) Microtubule-associated protein tau: a marker of paclitaxel sensitivity in breast cancer. Proc Natl Acad Sci U S A 102: 8315–8320.
- 44. Nakshatri H, Mendonca MS, Bhat-Nakshatri P, Patel NM, Goulet RJ Jr, et al. (2000) The orphan receptor COUP-TFII regulates G2/M progression of breast cancer cells by modulating the expression/activity of p21(WAF1/CIP1), cyclin D1, and cdk2. Biochem Biophys Res Commun 270: 1144–1153.
- 45. Navab R, Gonzalez-Santos JM, Johnston MR, Liu J, Brodt P, et al. (2004) Expression of chicken ovalbumin upstream promoter-transcription factor II enhances invasiveness of human lung carcinoma cells. Cancer Res 64: 5097–5105.
- 46. Maxwell GL, Chandramouli GV, Dainty L, Litzi TJ, Berchuck A, et al. (2005) Microarray analysis of endometrial carcinomas and mixed mullerian tumors reveals distinct gene expression profiles associated with different histologic types of uterine cancer. Clin Cancer Res 11: 4056–4066.
- 47. Stanbrough M, Bubley GJ, Ross K, Golub TR, Rubin MA, et al. (2006) Increased expression of genes converting adrenal androgens to testosterone in androgen-independent prostate cancer. Cancer Res 66: 2815–2825.
- 48. Wu F, Ivanov I, Xu R, Safe S (2009) Role of SP transcription factors in hormone-dependent modulation of genes in MCF-7 breast cancer cells: microarray and RNA interference studies. J Mol Endocrinol 42: 19–33.