Abstract
Triple-negative breast cancer (TNBC) is a subgroup of breast cancer that is negative for estrogen and progesterone receptor and ERBB2 protein expression. It is characterized by its aggressive behavior and by the lack of targeted therapies. To identify new therapeutic targets in TNBC, we used real-time quantitative RT-PCR to analyze 63 TNBC samples in terms of their mRNA expression of 26 genes coding for the major proteins currently targeted by drugs used to treat other cancers or undergoing clinical trials in breast cancer. Six of the 26 genes tested (VEGFA, SRC, PARP1, PTK2, RAF1, and FGFR3) were significantly upregulated in 13% to 46% of the TNBCs. None of the 6 genes was specifically upregulated in the TNBCs compared with 3 other classical breast tumor subtypes. No association was observed between overexpression of these 6 genes (except for FGFR3) and PIK3CA mutation status. These results confirm the interest of targeting VEGFA and PARP1 in ongoing clinical trials in TNBC patients and also identify new target genes (SRC, PTK2, RAF1, and FGFR3). Clinical trials could be initiated easily with existing drugs. Our results also suggest that these target genes might serve as predictive biomarkers of the TNBC treatment response.
Introduction
Breast cancer is the most common female malignancy in industrialized countries, affecting more than a million women per year worldwide. Breast cancer exhibits heterogeneous behavior, clinical outcomes, and treatment responses. Better understanding of breast cancer biology has led to significant improvements in patient survival. The discovery of steroid hormone dependence resulted in the development of estrogen receptor antagonists and aromatase inhibitors, which are currently the gold standard treatments for patients with hormone-receptor–positive breast tumors. 1 Similarly, the detection of ERBB2 overexpression/amplification led to ERBB2 targeting with trastuzumab, and lapatinib was recently shown to significantly improve the survival of patients with ERBB2-overexpressing tumors. Hormone (estrogen and progesterone) receptor status, ERBB2 overexpression/amplification, and Ki67 expression are now used to predict the prognosis of breast cancers and to guide treatment. 2 However, chemotherapy is the only available systemic therapy for women with so-called triple-negative breast cancer (TNBC), which lacks estrogen receptor (ER) and progesterone receptor (PR) expression and ERBB2 overexpression/gene amplification.
TNBC represents about 15% to 20% of breast cancers. It is characterized by an aggressive clinical course and poor prognosis, owing partly to the lack of targeted therapies.3-5 Most TNBCs have the “basal-like” molecular profile in gene expression arrays, but recent studies have suggested that TNBC is a heterogenous type of cancer; indeed, less common subtypes including “Claudin-low,” “HER2-enriched but without HER2 gene amplification,” and “molecular apocrine” have also been described in TNBC.6-8 Histologically and transcriptionally, TNBCs have many similarities to BRCA1-associated breast cancers, and most BRCA1-associated breast tumors are triple-negative and basal-like. 9 BRCA1 is rarely mutated in sporadic breast cancer, but it has been suggested that BRCA1 (or associated pathways) is inactivated in triple-negative tumors via other molecular mechanisms. Toyama et al. 10 showed that BRCA1 mRNA expression was significantly decreased in TNBCs compared with luminal subtype breast cancers.
The overlap between BRCA1-mutated breast cancers and triple-negative tumors suggests that some triple-negative tumors might respond to therapeutics targeting BRCA1-deficient cells, such as PARP inhibitors. 11
The addition of iniparib, a PARP inhibitor, to chemotherapy improved the clinical benefit and survival of patients with metastatic TNBC without significantly increased toxic effects. 12 On the basis of these results, a phase 3 trial evaluated overall survival and progression-free survival for women with metastatic TNBC but was negative. Given the structural and mechanistic differences between iniparib and other PARP inhibitors, these negative results do not necessarily imply a class effect, and further study of TNBCs with other PARP inhibitors should be encouraged, 13 unless the drugs work in some molecular subtypes of TNBC but not others.
There is increasing evidence that the DNA-repair defects characteristic of BRCA1-related cancers, and especially defective homologous recombination, confer tumor sensitivity to certain systemic agents. Indeed, patients with TNBC have higher pCR (partial Complete Remission) rates than patients with non-TNBC, and neoadjuvant trials have shown higher relapse-free survival in TNBC patients who achieve pCR than in patients with residual disease. 14 Despite this relative chemosensitivity, local and systemic TNBC relapse rates remain higher than in other breast cancer subtypes.15,16
Molecularly directed therapy targets tumor cells and the tumor microenvironment by blocking the effects of tumor-specific molecular changes. Targeted treatments are directed at a specific molecular target that is not present in normal breast cells and that is important for tumor growth and progression.
Targeted treatments tend to have fewer adverse effects, but their use must be guided by biomarker assays. For example, immunohistochemical assays are used to identify the therapeutic target in the breast tumor before prescribing hormone therapy or trastuzumab, and detection of EGFR activating mutations is an obligatory prerequisite to EGFR inhibitor prescription in lung cancer.
To identify new therapeutic targets in TNBC, we applied real-time quantitative RT-PCR assays to 63 triple-negative tumor samples. We quantified the mRNA expression of a panel of 26 genes coding for the major proteins that are currently targeted by drugs used to treat other cancer types or that are undergoing clinical trials in breast cancer.
Results
mRNA expression of the 26 target genes in the 63 triple-negative breast tumors
We used real-time quantitative RT-PCR to analyze mRNA expression of the 26 target genes in a series of 63 TNBCs and 12 normal breast tissues. The mRNA levels of all 26 target genes were high in both the normal and tumorous breast tissues and were thus reliably quantifiable by real-time quantitative RT-PCR based on fluorescence SYBR Green method (Cycle Threshold, Ct < 32). Target gene mRNA levels in the 63 TNBCs were expressed relative to the median mRNA levels observed in the 12 normal breast tissues. For each gene, normalized mRNA values of 3 or more were considered to represent gene overexpression in tumor samples, and values 0.33 or less represented gene underexpression. Medians and ranges of mRNA levels for the 26 target genes are shown in Table 1, along with the percentages of overexpression or underexpression.
mRNA Expression of MKI67 and the 26 Target Genes in TNBC Tissues Relative to Normal Breast Tissues and Percentages of Overexpression and Underexpression
Note: NS = not significant; TNBC = triple-negative breast cancer.
Kruskal-Wallis H test.
Median (range) of mRNA levels relative to normal.
Eighteen (69.2%) of the 26 genes were significantly dysregulated in the TNBCs. Six (23.1%) genes were mainly upregulated (FGFR3, PARP1, PTK2, RAF1, SRC and VEGFA), and 12 (46.1%) were down-regulated.
In the same set of 63 TNBC samples, we also examined the expression of MKI67, which encodes the proliferation-related antigen Ki-67. As expected, MKI67 was upregulated in all tumor samples.
Table 1 shows the mRNA expression levels of the 26 target genes in breast tumors relative to the TBP endogenous control. The same results were obtained when other endogenous RNA controls (RPLP0 or PPIA) were used.
Comparison of mRNA levels of the 6 upregulated genes according to the tumor subtype
We then examined whether the 6 genes upregulated in TNBCs were specific to this tumor subtype by analyzing their mRNA expression in the other 3 major breast tumor subtypes: NNP (ERα-negative, PR-negative, ERBB2-positive tumors), PPN (ERα-positive, PR-positive, ERBB2-negative tumors), and PPP (ERα-positive, PR-positive, ERBB2-positive tumors). The results are shown in Figure 1.

Distribution and frequency of the overexpression of the 6 upregulated genes (VEGFA, PARP1, SRC, PTK2, RAF1, and FGFR3) in TNBCs relative to the 3 other main breast cancer subgroups, NNP (ERα-negative, PR-negative, ERBB2-positive tumors), PPN (ERα-positive, PR-positive, ERBB2-negative tumors), and PPP (ERα-positive, PR-positive, ERBB2-positive tumors).
None of the 6 genes were specifically upregulated in the TNBCs compared with the other 3 tumor subtypes. Nevertheless, it is noteworthy that VEGFA was predominantly upregulated in ERα-negative tumors (TNBC and NNP subtypes) (P = 0.00001 vs. the PPN and PPP subtypes). FGFR3 was significantly more strongly overexpressed in the ERα-positive tumors (PPN and PPP) than in the ERα-negative tumors (TNBC and NNP) (P = 0.00002). SRC overexpression was significantly stronger in the ERBB2-positive tumors (NNP and PPP) than in the ERBB2-negative tumors (PPN and TNBC) (P = 0.006).
Comparison of mRNA levels of the 6 upregulated genes according to PIK3CA status
In the TKR–RAS-PI3K-AKT signaling pathway, PIK3CA is the oncogene that shows the highest frequency of gain-of-function mutations in breast cancer. Indeed, EGFR, BRAF, and KRAS, although frequently mutated in other cancers (colon, lung, etc.), are rarely mutated in breast cancer. PIK3CA status may play an important role in the response to therapies targeting tyrosine kinase receptors in breast cancer.17,18 We therefore examined the overall series of 154 tumor samples (4 subgroups) for PIK3CA mutations in exons 9 and 20. PIK3CA mutations were detected in 6 (9.5%) of 63 TNBCs, 7 (21.8%) of 32 NNP tumors, 16 (41.0%) of 39 PPN tumors, and 7 (35.0%) of 20 PPP tumors. PIK3CA mutations were thus detected in 13 (13.7%) of the 95 ERα-negative tumors and in 23 (39.0%) of the 59 ERα-positive tumors, in agreement with the literature. 19 We then looked for an association between the 6 overexpressed genes and PIK3CA mutation status. A positive association, at the limit of statistical significance (P = 0.045), was only observed between FGFR3 overexpression and PIK3CA mutation status.
Discussion
We used real-time RT-PCR to analyze the expression of 26 selected genes in a large sample of breast tumors. We found that 6 genes (VEGFA, FGFR3, PARP1, SRC, PTK2/FAK, and RAF1) were frequently upregulated in the TNBC subgroup relative to normal breast tissue. However, none of these 6 genes was specifically upregulated in TNBCs compared with the other 3 main breast tumor subgroups (NNP, PPN, and PPP), further highlighting the heterogeneity of TNBC subgroup and the difficulty of finding a specific target.6-8
Several genes were also down-regulated in TNBCs. It is more difficult to restore tumor-suppressor protein expression than to inhibit oncoprotein overexpression. Moreover, the observed down-regulations could be partly explained by differences in proportion of epithelial cells and stromal cells (fibroblasts, adipocytes, endothelial cells, and circulating cells) between normal control breast tissue and tumor samples. Likewise, underexpression of stroma cell–specific genes could be explained by a lower abundance of a particular cell type (fibroblasts, adipocytes, endothelial cells, or circulating cells) in TNBCs relative to normal breast tissue.
By using high-resolution oligonucleotide comparative genomic hybridization arrays, Andre et al. 20 demonstrated that 3 genes including VEGFA were specifically gained in TNBC. Moreover, those investigators established that the occurrence of a DNA gain leads to an unregulated overexpression of mRNA and concluded that such dysregulated genes may represent novel therapeutic targets.
It is noteworthy that PARP inhibitors are effective on tumors that carry a DNA-repair defect, such as BRCA1-deficient tumors.11,21 The most appropriate biomarkers for anti-PARP sensitivity will thus be markers of DNA-repair defects, particularly homologous recombination DNA-repair defects, in non–BRCA1- and non–BRCA2-associated tumors. However, PARP upregulation, which was demonstrated at the protein level by Ossovskaya et al., 22 could be an additional biomarker of PARP inhibitor sensitivity in triple-negative breast tumors. These results and ours for VEGFA and PARP1 support the validity of current clinical trials testing bevacizumab and PARP inhibitors such as BSI-201 4 and AZD2281/Olaparib 21 in TNBC patients. Concerning the other 4 genes identified (RAF1, SRC, FGFR3, and FAK), clinical trials in breast cancer are less advanced.
Recently it has been shown that accumulation of recurrent RAF1 gene amplification is due to EZH2 expression–mediated downregulation of DNA damage repair in breast tumor initiating cells (BTICs), which activates p-ERK-β-catenin signaling to promote BTIC expansion. 23 Overexpression of Polycomb protein EZH2, essential in stem cell self-renewal, has been linked to breast cancer progression. This amplification of RAF1 may be compatible with the RAF1 overexpression that we detected in TBNCs. RAF1 is targeted by sorafenib, a drug currently used in the treatment of kidney cancer and hepatocellular carcinoma. Several clinical trials of chemotherapy with or without sorafenib will start shortly in patients with advanced breast cancer. 24
Several inhibitors of SRC, such as dasatinib and bosutinib, are currently being evaluated in breast cancer patients. 25 It has been suggested that SRC could be a potential target for the treatment of TNBC because SRC expression was more intense in cytoplasm and on the membrane of triple-negative samples than in non–triple-negative samples. 26 However, our results suggest that it is mainly patients with ERBB2-overexpressing tumors who are likely to benefit from SRC inhibitors (Figure 1). Recently, Zhang et al. 27 showed that increased SRC activation conferred considerable trastuzumab resistance to breast cancer cells and correlated with clinical trastuzumab resistance. Thus, SRC targeting might overcome trastuzumab resistance. 27
According to our results at the transcriptional level, FGFR3 expression in invasive breast cancer was not significantly associated with specific clinicopathological/molecular parameters. 28 Yom et al. 29 demonstrated that FAK overexpression by immunohistochemistry and amplification by FISH were significantly associated with the triple-negative subgroup, whereas our transcript results did not associate FAK overexpression with a particular subgroup. No drugs targeting FGFR3 or FAK are currently being tested in breast cancer, although results cited here support such trials. Dovitinib (a specific FGFR3 inhibitor), PD173074 (a highly potent selective pan-FGFR inhibitor), and BIBF-1120 (a multityrosine kinase inhibitor), which target FGFR3 as well as PF-00562271 and PND-1186 (two FAK inhibitors), are currently in phase 1 trials.30-34
It is noteworthy that Turner et al. 35 identified FGFR2 gene amplification (associated with strong overexpression) in 4% (6/165) of TNBCs but not in other subtypes (0/214). Our results show FGFR2 overexpression in 8% of 63 TNBCs but also in 8% of 39 PPNs and 5% of 20 PPPs. None of 32 NNPs overexpressed FGFR2 (data not shown). These results confirm this rare genetic alteration as a potential therapeutic target in breast cancer but not specifically for TNBCs.
The upregulated genes identified here are only upregulated in a fraction (13%-46%) of TNBCs. It will therefore be necessary to test their value as predictive biomarkers of the response to targeted drugs. However, the best such a biomarker is not necessarily the target itself. Indeed, VEGFA (overexpressed in 46% of TNBCs), which codes for a secreted protein measurable in serum, does not seem to predict the response to the anti-VEGF antibody bevacizumab (Avastin®) in some cancers. 36
Concerning the mutation status of the PIK3CA gene, which could play a major role in the response to therapies targeting tyrosine kinase receptors in breast cancer,17,18 we found a positive correlation between PIK3CA mutation and FGFR3 overexpression. This correlation seems to be indirect, being related mainly to ERα-positive tumor status. Indeed, FGFR3 is upregulated in ERα-positive tumors, and PIK3CA mutation is also more frequent in our ERα-positive tumor population (39% vs. 13.7% in the ERα-negative population, in agreement with the literature). 19 This suggests that PIK3CA mutation does not influence the expression of these 6 genes.
In conclusion, our results confirm the validity of VEGFA and PARP1 targeting in ongoing clinical trials in TNBC patients and also identify new potential target genes (SRC, PTK2, RAF1, and FGFR3). Clinical trials of drugs inhibiting the products of these genes could be initiated rapidly with existing drugs. Our results also warrant studies of these target genes as putative predictive biomarkers for selecting the patients most likely to benefit from these drugs.
Materials and Methods
Patients and samples
We analyzed tumor samples from 63 TNBC patients treated at Institut Curie, René Huguenin Hospital, Saint-Cloud, France. Tumor samples containing more than 70% of tumor cells were considered suitable for analysis. Immediately after surgery, the tumor samples were placed in liquid nitrogen until RNA extraction.
The patients met the following criteria: primary unilateral nonmetastatic breast carcinoma; complete clinical, histological, and biological information available; no radiotherapy or chemotherapy before surgery; and full follow-up at Institut Curie, René Huguenin Hospital. Patients underwent physical examinations every 3 months for 2 years, then annually. Mammograms were done annually. Median follow-up was 7.5 years (range 8 months to 29 years). Twenty-six patients relapsed; the distribution of first relapse events was as follows: 22 metastases alone and 4 with both metastases and local relapse. Standard prognostic factors are shown in Table 2.
Characteristics of the 63 TNBCs
Note: NS = not significant; SRB = Scarff Bloom Richardson classification.
Log-rank test.
Data available for 58 patients.
Data available for 62 patients.
ER and PR were routinely analyzed at the time of diagnosis on frozen tumors by ligand binding assay until 1988, by enzyme immunoassay (ER-EIA Monoclonal, PgR-EIA Monoclonal, Abbott Laboratories, Abbott Park, IL) between 1988 and 2000, and then by immunohistochemistry on paraffin sections. HER2 status was routinely analyzed by immunohistochemistry (with confirmation by FISH of the 2+ cases). For this study, the triple negative status of the tumors was confirmed by RT-qPCR on frozen tumors.37,38
To investigate specific dysregulation of candidate genes in the TNBC/NNN subtype, we analyzed a panel of RNAs from the 3 other major breast tumor subtypes: 32 NNP (ERα-negative, PR-negative, ERBB2-positive tumors), 39 PPN (ERα-positive, PR-positive, ERBB2-negative tumors), and 20 PPP (ERα-positive, PR-positive, ERBB2-positive tumors).
Twelve specimens of adjacent normal breast tissue from breast cancer patients or normal breast tissue from women undergoing cosmetic breast surgery were used as sources of normal RNA.
Real-time RT-PCR
The theoretical and practical aspects of real-time quantitative PCR have been described in detail elsewhere. 39
By studying the literature, we selected 26 genes coding for the major proteins currently targeted by drugs used to treat other cancers or for proteins targeted in ongoing breast cancer clinical trials (Suppl. Table S1). The 26 target genes tested in this study are listed in Table 3.
The 26 Selected Genes
The precise amount of total RNA added to each reaction mix (based on optical density) and its quality (i.e., lack of extensive degradation) are both difficult to assess. We therefore also quantified transcripts of 3 endogenous RNA control genes involved in various cellular metabolic pathways, namely TBP 37 (Genbank accession NM_003194), which encodes the TATA box-binding protein (a component of the DNA-binding protein complex TFIID); RPLP0 40 (also known as 36B4; NM_001002), which encodes human acidic ribosomal phosphoprotein P0; and PPIA, 41 which encodes peptidylprolyl isomerase A (also known as cyclophilin A; NM_021130).
Primers for TBP, RPLP0, ERα, PR, ERBB2, MKI67, and the 26 target genes were chosen with the assistance of the Oligo 5.0 computer program (National Biosciences, Plymouth, MN). We searched the dbEST and nr databases to confirm the total gene specificity of the nucleotide sequences chosen as primers and the absence of single nucleotide polymorphisms. In particular, the primer pairs were selected to be unique relative to the sequences of closely related family member genes or of the corresponding retropseudogenes. The nucleotide sequences of the oligonucleotide primers used to amplify MKI67 and the 26 target genes are shown in Suppl. Table S2.
Each sample was normalized on the basis of its TBP (or RPLPO or PPIA) content. Results, expressed as N-fold differences in target gene expression relative to the TBP (or RPLPO) gene and termed “Ntarget,” were determined as Ntarget = 2ΔCtsample, where the ΔCt value of the sample is determined by subtracting the average Ct value of the target gene from the average Ct value of the TBP (or RPLP0 or PPIA) gene.37-39 The Ntarget values of the samples were subsequently normalized such that the median of the 12 normal breast tissue Ntarget values was 1.
RNA extraction, cDNA synthesis, and PCR conditions were as described elsewhere. 39
PIK3CA mutation screening
PIK3CA mutations were screened for in cDNA fragments obtained by RT-PCR amplification of exons 9 and 20 with their flanking exons. Details of the primers and PCR conditions are available on request. The amplified products were sequenced by using the BigDye terminator sequencing kit on an ABI Prism 3130 automatic DNA sequencer (Applied Biosystems, Courtabæuf, France). Sequences were compared with the corresponding cDNA reference sequence (NM_006218).
Statistical analysis
As the mRNA levels did not fit a Gaussian distribution, (a) the mRNA levels in each subgroup of samples were characterized by their median values and ranges rather than their mean values and coefficients of variation, and (b) relationships between the molecular markers and clinical and biological parameters were tested by using the nonparametric Kruskal-Wallis test (links between 1 qualitative parameter and 1 quantitative parameter). Differences between 2 populations were considered significant at confidence levels greater than 95% (P < 0.05).
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by Comité départemental des Hauts-de-Seine de la Ligue Nationale Contre le Cancer, Conseil régional d’Ile-de-France, Cancéropôle Ile-de-France and by the Association d’Aide à la Recherche Cancérologique de Saint-Cloud (ARCS).
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
