- Open Access
A two-gene epigenetic signature for the prediction of response to neoadjuvant chemotherapy in triple-negative breast cancer patients
Clinical Epigenetics volume 11, Article number: 33 (2019)
Pathological complete response (pCR) after neoadjuvant chemotherapy (NAC) in triple-negative breast cancer (TNBC) varies between 30 and 40% approximately. To provide further insight into the prediction of pCR, we evaluated the role of an epigenetic methylation-based signature.
Epigenetic assessment of DNA extracted from biopsy archived samples previous to NAC from TNBC patients was performed. Patients included were categorized according to previous response to NAC in responder (pCR or residual cancer burden, RCB = 0) or non-responder (non-pCR or RCB > 0) patients. A methyloma study was performed in a discovery cohort by the Infinium HumanMethylation450 BeadChip (450K array) from Illumina. The epigenetic silencing of those methylated genes in the discovery cohort were validated by bisulfite pyrosequencing (PyroMark Q96 System version 2.0.6, Qiagen) and qRT-PCR in an independent cohort of TN patients and in TN cell lines.
Twenty-four and 30 patients were included in the discovery and validation cohorts, respectively. In the discovery cohort, nine genes were differentially methylated: six presented higher methylation in non-responder patients (LOC641519, LEF1, HOXA5, EVC2, TLX3, CDKL2) and three greater methylation in responder patients (FERD3L, CHL1, and TRIP10). After validation, a two-gene (FER3L and TRIP10) epigenetic score predicted RCB = 0 with an area under the ROC curve (AUC) = 0.905 (95% CI = 0.805–1.000). Patients with a positive epigenetic two-gene score showed 78.6% RCB = 0 versus only 10.7% RCB = 0 if signature were negative.
These results suggest that pCR in TNBC could be accurately predicted with an epigenetic signature of FERD3L and TRIP10 genes. Further prospective validation of these findings is warranted.
At the present time, chemotherapy (CT) is the only proven therapy for triple-negative breast cancer (TNBC) subtype. Anthracycline and taxane-based CT is still the standard of care for TNBC , with pathological complete response (pCR) rates ranging ~ 30–40% [2,3,4,5,6]. pCR rate in TNBC are associated with better outcomes while residual disease after neoadjuvant chemotherapy (NAC) have a higher relapse risk and poor prognosis . With the advantage of high sequencing technology, several molecular signatures have been developed in the recent years to predict response to neoadjuvant chemotherapy (NAC). Oncotype-Dx , MammaPrint , Blue Print , Endopredict , or Prosigna  are some of them. Recently, an initial 199-gene signature, (E2F4 target gene signature), has shown accurate prediction of response to NAC even when reduced to 33-gene panel and has been validated in 1129 patients across five independent data sets . However, all these predictive panels have shown to do better in ER positive breast cancer than in ER negative. Thus, accurate prediction of response in TNBC still remains a medical need. Other study addressed the issue of the prediction to NAC in 94 patients TNBC treated with paclitaxel and carboplatin according to the Lehman’s TNBC type 4 classification. Basal-like 1 (BL1) subtype showed a pCR rate of 65.2%, while basal-like 2 (BL2) was 47.4%, mesenchymal (M) 36.4%, and luminal androgen receptor (LAR) 21.4% . The I-SPY 2 is an ongoing prospective trial of NAC in breast cancer including a cohort of TNBC patients. A 70-gene panel combined with DNA deficient biomarkers have shown a 75% pCR in a subset of TNBC patients treated with carboplatin and veliparib in a recent publication of this cohort . All these panels have been designed based on arrays of gene-expression techniques. Nevertheless, epigenetic modifications of certain genes can lead to silence or activation of different genes [16, 17].
The epigenetic modifications of the DNA, such as methylation, can modulate gene expression with no DNA sequence modification and contribute to disease development . In this context, epigenetic changes in tumor DNA before CT administration could potentially have a predictive role of response to this therapy . The aim of this study was to identify a predictive epigenetic signature of pCR as defined by the residual cancer burden (RCB) index by Symmans et al. (RCB = 0) in patients with TNBC treated with NAC including anthracyclines and/or taxanes-based regimens.
Clinical characteristics of TNBC patients
Fifty-four patients were included: 24 in the discovery cohort (DC) and 30 in the validation cohort (VC). The clinical characteristics of the patients are summarized in Table 1. After biopsy, all patients were treated with NAC based on a taxane and/or anthracycline regimen and all of them were considered TNBC according to immunohistochemistry for ER, PR, and HER2. Patients were classified in responders (R) if RCB = 0 or non-responders (NR) if RCB > 0.
Analysis of DNA methylome in responder and non-responder patients: discovery cohort
A genome-wide DNA methylation study was performed in the DC (N = 24, 10 R (RCB = 0) and 14 NR (RCB > 0)). Figure 1a summarizes the whole process during this study. The analysis of methylation data showed 133 CpGs sites (71 genes) with differences in methylation levels ≥ 20% (p value < 0.05) that distinguished R patients (treatment sensitive) from NR patients (treatment resistant) (Fig. 1b). According to a Gene Ontology (GO) analysis, some of these genes were involved in biological functions and pathways such as DNA repair, cell adhesion, transcription regulation, or signaling mediated by GTPases (Fig. 1c, Additional file 1: Table S1) that have shown to be implicated in chemoresistance of cancer, including the response to anthracyclines and/or taxanes [20,21,22,23].
Thirty-five CpGs located in promoters, islands, or shores from 23 genes were selected to further validations (Additional file 2: Table S2). Of these, taking into account an intra-group SD ≤0.2, 11 CpGs corresponding to 11 genes showed significant methylation differences (Additional file 3: Table S3), and 9 of these genes showed a consistent methylation profile on consecutive CpGs (Fig. 2). These candidate genes were LOC641519, LEF1, HOXA5, EVC2, TLX3, and CDKL2 with high methylation in NR group and genes FERD3L, CHL1, and TRIP10 with high methylation in R group.
Validation of methylation in candidate genes by pyrosequencing: validation cohort
A pyrosequencing study in the DC and in the VC (N = 30, 9 RCB = 0 and 21 RCB > 0) was performed to validate the candidate genes first technically and secondly in an independent cohort, respectively. Methylation was analyzed for each gene taking into account the differentially methylated CpGs identified in the 450K array and other CpGs located close to them in order to obtain a more consistent result (Additional file 4: Table S4).
In the DC, we replicated by pyrosequencing the methylation data obtained in 450K array for LOC641519/LEF1 gene (p value = 0.02) and HOXA5 gene (p value = 0.0001), where we also observed a methylation level significantly higher in NR patients than R patients, and in FERD3L gene (p value = 0.04), TRIP10 (p value = 0.003), and CHL1 (p value = 0.03), where methylation was also significantly higher in R patients than NR patients. However, in EVC2 gene (p value = 0.07), CDKL2 gene (p value = 0.05), and TLX3 gene (p value = 0.07), replication was not statistically significant but showed a trend towards a higher level of methylation in NR patients (Fig. 3a).
Pyrosequencing in the VC validated the results for FERD3L gene (p value = 0.0087) with high methylation in R group versus NR group (Fig. 3b). Moreover, differences in TRIP10 methylation showed a trend towards significance between R and NR (p value = 0.19). Accordingly, by means of a biological pathway analysis of these genes using the publicly available resource Pathway Commons , we observed that FERD3L and TRIP10 are able to interact with several genes (Additional file 5: Figure S1) which have been previously associated with therapy resistance of breast cancer and other types of tumors [23, 25].
FERD3L methylation and gene expression in TNBC cell lines
In order to evaluate whether DNA methylation has a functional role in the transcriptional control of FERD3L, we evaluated FERD3L gene expression by qRT-PCR and methylation by pyrosequencing in a set of TNBC cell lines (Fig. 4a) with the aim of correlating FERD3L methylation levels with gene expression level and corroborate the results from patients. We observed that FERD3L gene was methylated in all the cell lines studied with levels always higher than 40%. The FERD3L expression inversely correlated with the methylation detected as was expected, showing a low gene expression when methylation was high. Thereby, MDA-MB-231 cell line showed the lowest methylation level and correlated with the highest level of gene expression. Conversely, HCC-1143 cell line that showed the higher methylation was the one with the lowest gene expression level.
FERD3L demethylation assays and gene expression in TNBC cell lines
In order to check if changes in FERD3L methylation status also affect to gene expression level, we performed an assay with the MDA-MB-436 cell line treated with AZA demethylating agent. We observed that the treatment modified FERD3L methylation in MDA-MB-436 cell line inducing a decrease when compared with control cells (MDA-MB-436 cell line not treated with AZA) (p value = 0.05). As we expected, this change in FERD3L methylation was correlated with a significant increase (p value = 0.0022) in FERD3L expression (Fig. 4b).
High FERD3L expression levels correlates with low gene methylation in TNBC patients
The analysis of FERD3L gene expression in the 54 patients showed a significant difference (p value = 0.04) in FERD3L gene expression with high expression in NR patients versus R patients (Fig. 4c). Therefore, it suggests an inverse correlation between methylation and gene expression in NR patients and R patients. The Cancer Genome Atlas dataset (TCGA) analysis for 713 breast cancer patients showed negative correlation between methylation and expression, according with our data. It was detected both when all CpGs in the FERD3L gene promoter were included and when only analyzed the CpG cg10043037 validated for FERD3L gene in the study (Fig. 4d).
Statistical model to predict response to neoadjuvant treatment in TNBC patients
Due to FERD3L and TRIP10 showed the higher level of significance in the VC, we selected both genes for therapy response analysis. It is interesting to note that these two genes did not show statistical differences in methylation between R and NR in terms of age, tumor size, nodule affectation, or ki67 expression (Additional file 6: Table S5) indicating that these genes were not associated with any relevant clinicopathological prognostic factor. Importantly, we were able to create a statistical epigenomic predictive model of pathological response (RCB = 0) with the FERD3L and TRIP10 methylation and using the whole cohort. These two genes were selected as both showed the higher level of significance in the VC. The statistical model for the prediction of therapy response was based on the Akaike information criterion (AIC) and by constructing a receiver operating characteristic (ROC) curve. Based on this model, the following rule was constructed:
A × FERD3L methylation level (%) + B × TRIP10 methylation level (%) > 971 (A and B being constants)
Levels > 971 showed a high likelihood for RCB = 0 (8 out of 11 cases 78.6%), while levels ≤ 971 showed a poor probability for RCB = 0 (3 out of 28 patients 10.7%). Value based of the calculation of the ROC curve with AUC = 0.9056 (95% CI = 0.805–1.000) (Fig. 5a, b).
This study evaluated the methylation profile of triple-negative breast cancer (TNBC) patients treated with neoadjuvant systemic chemotherapy (NAC) depending on the treatment response and identified a two-gene epigenetic signature for discriminating responders (R) from non-responders (NR). The discovery cohort included 24 TNBC patients and was analyzed by using Infinium Human Methylation 450 BeadChip array whereas the validation cohort included 30 TNBC patients and was analyzed by pyrosequencing for the significant and most relevant genes identified in the discovery cohort. For each case, DNA was obtained from core biopsies before the beginning of the treatment. The analysis of the discovery cohort identified nine genes differentially methylated. Six genes with higher methylation in NR patients (LOC641519, LEF1, HOXA5, EVC2, TLX3, CDKL2) and three genes with greater methylation in R patients (FERD3L, CHL1 and TRIP10). After technical and analytical validation by pyrosequencing in both discovery and validation cohorts, we identified a two-gene (FERD3L and TRIP10) signature able to predict response to NAC. The role of promoter methylation in the regulation of FERDL3 gene expression was evaluated in TNBC cell lines and in TNBC tissues, demonstrating an inverse correlation between methylation and expression levels. To our knowledge, the two genes epigenetic model shown in this study is the first epigenetic signature for prediction of response to NAC in TNBC patients. It should be noted that this model almost doubles the predictive potential described for the TNBC subtype by other approaches (~ 30–40% versus 78.6% with our model) [2,3,4,5,6].
Pathological complete response (pCR) is an effective surrogate marker for survival among patients with luminal B/HER2 negative, HER2 positive, and TNBC tumors. The results of a meta-analysis of 6377 patients with operable or advanced non-metastatic disease from six prospective neoadjuvant studies support this . In the TNBC group of 911 patients who received anthracycline and taxane based NAC, 31% achieve pCR. Most recently, a retrospective analysis of 452 TNBC patients showed a pCR of 33% . Overall and disease-free survival were significantly longer among patients achieving pCR versus residual disease . In this context, the development of a more accurate predictive signature in TNBC could have an important clinical impact. To date, available gene panels for prediction of pCR to NAC are based in genomic signatures [8,9,10, 13, 27].
Epigenetics changes affect cellular processes such as gene expression and have clearly been related to the development of diseases such as breast cancer. Few methylation studies have been specifically done on TNBC subtype. Stizarken et al. identified differentially methylated regions that could separate TNBC and non-TNBC patients and classified those according to prognosis. This provided the first evidence that changes in the methylation profile of DNA could be useful to identify and stratify TNBC patients . Recently, Mathe et al. performed a study of gene expression and DNA methylation in the same population, demonstrating that DNA methylation contributes to the deregulation of gene expression . In this context, our results suggest that response to NAC can be predicted accurately with an algorithm of the methylation status of FERD3L and TRIP10 genes in patients with TNBC.
The FERD3L (Fer3-like bHLH transcription factor) gene, also named NATO3 or N-TWIST, is a gene located on chromosome 7 and is a basic helix-loop-helix (bHLH) transcription factor. These factors play an essential role in multiple developmental processes, mainly in neurogenesis, where its regulation is essential for the right development . FERD3L is a member of the TWIST genes family that is implicated in epithelial-mesenchymal transition (EMT) in cancer cells, a process also related to metastasis and may lead to chemo-resistance in TNBC .
Several studies have shown that in metastatic carcinomas including aggressive and metastatic breast cancer, there are an overexpression of TWIST . The inactivation of TWIST by siRNA technology or chemotherapeutic approaches has proved successful [36,37,38], so it is presented as a potential therapeutic target for metastatic breast cancer. There is only one work in the literature that relates the FERD3L gene to cancer, specifically neuroblastoma . Promoter CpG islands in FERD3L gene was found to be highly methylated in neuroblastoma cell lines causing gene silencing and poor prognosis. However, to our knowledge, there is no previous evidence of relation between FERD3L expression and response to CT in the clinical setting or any relation with breast cancer.
The TRIP10 (thyroid hormone receptor interactor 10 gene) is located in chromosome 19 and belongs to the minor histocompatibility antigens family and codifies the Cdc42-interacting protein 4 (CIP4). This protein interacts with the GTPase Cdc42 that is related with actin formation and has been implicated in cytoskeleton organization . Paclitaxel, a type of taxane, promotes microtubule stabilization and polymerization leading to a cell cycle arrest and apoptosis . Actin-microtubule crosstalk is particularly important for cell shape and polarity during cell migration and division . In this context, TRIP10 hypermethylation could increase the efficacy of paclitaxel effects on cytoskeleton. Furthermore, TRIP10 has an important role in the cellular motility and cohesion control since it is implicated in E-cadherin regulation . In fact, previous studies have related TRIP10 with cell invasion in TNBC cells in vitro . Interestingly, another study demonstrated also that TRIP10 gene controls EMT . They described a pro-metastatic role of TRIP10, in concordance with the in vitro and in vivo data from Cerqueira et al. . In addition, TRIP10 expression has shown to be regulated by DNA methylation in mesenchymal stem cell differentiation  and in several types of cancer cell lines and tumors . In particular, in breast cancer, it has been also confirmed that TRIP10 expression can be regulated by epigenetic mechanisms such as DNA methylation [47, 48].
Basing a signature on the methylation status of only two genes could be a limitation of this study. However, the identification of this signature was provided by a very stringent analysis to select the most suitable candidate genes and they were further validated in an independent cohort of patients. In fact, there are several examples of other two-gene signatures, including methylation based signatures, that have previously shown clinical utility for cancer in different types of tumors [49,50,51,52,53,54,55,56]. Additionally, the sample size could be considered small; however, it is important to highlight that the group of patients is very homogeneous in order to eliminate potential confusion factors. In our analysis, we compared methylation levels among TNBC patients who responded and did not respond to neoadjuvant chemotherapy. Both groups were homogeneous also regarding to staging and severity, basing the difference only in the response to chemotherapy. In fact, no statistically significant differences were observed in methylation levels according to staging and severity. Therefore, the strict selection of patients could be a strength that give support and power to these results. The low number of cases analyzed could limit the immediate translational relevance but represents a very good start point for future studies in the field.
Triple-negative breast cancer has a high relapse rate after conventional chemotherapy treatment. To date, no predictors of treatment effectiveness have been identified. In this study, we propose an epigenetic signature based on the methylation levels of the FERD3L and TRIP10 genes. Our algorithm has a complete pathological response prediction potential of 78.6% and increases the predictive potential described by other approaches. This is especially relevant if we consider that it could be a predictive tool in clinical practice that will allow selecting the appropriate treatment as well as better stratification of patients for clinical trials.
Patients treatment and tumor samples
Patients treated with anthracyclines and/or taxanes NAC in the Hospital Clínico of Valencia and diagnosed with an early TNBC between 2005 and 2015 (Table 1) were retrospectively selected for the study according to clinical inclusion/exclusion criteria (Additional file 7: Table S6). Tumor samples were obtained before exposure to any systemic anticancer treatment using ultrasound-guided core needle biopsy. The cores were placed on OCT and stored at − 80 °C, or included in FFPE. Tumor percentage, histology and ER, PR, HER2, and Ki67 expression were determined. ER and PR status were considered negative when nuclear staining is < 10%. For the assessment of HER2, ASCO/CAP recommendations were used . Diagnosis of TNBC was done according to IHC results. The pathological response after NAC was evaluated by the Symmans method (residual cancer burden; RCB) . A value of RCB = 0 implies pCR, whereas values of RCB > 0 indicates that there is still residual tumor. From an initial analysis of the database, 70 patients were identified. Only those that had a tumor percentage > 25% and reached 500 ng after DNA extraction were used. Fifty-four patients were included in the study; of these, a group of 24 patients (10 RCB = 0 and 14 RCB > 0) were selected for the discovery cohort (DC) and 30 patients (9 RCB = 0 and 21 RCB > 0) were included in the validation cohort (VC).
Sample size calculation
Sample size calculation was based on data from a recent series reporting that TNBC patients treated with NAC obtaining a RCB = 0 was around 30% . If our methylation data would be able to predict two groups of responders (R) versus non-responders (NR) and we estimate that the proportion of patients with RCB = 0 could be 45% in R while the proportion of RCB = 0 in the NR group would be 10%, with an alpha error of 5% and a power of 80%, the sample size needed to identify this difference between both proportions is 44 patients (22 patients/group). According to this, the theoretical sample size, including a 15% drop-out, should be 51 patients, (https://select-statistics.co.uk/calculators/sample-size-calculator-two-proportions/).
TNBC cellular lines and treatments
Five TNBC cell lines from American Type Culture Collection (ATCC) were cultured (HCC-1937, HCC-1143, HCC-38, MDA-MB-231, and MDA-MB-436) following standard culture conditions. Treatment with the demethylating agent 5-aza-2′-deoxycytidine (AZA) (Sigma, St. Louis, MO, USA) was performed at 5 uM during 72 h. In this assay, results were performed in triplicates and data were compared with the corresponding non-treated cell line.
DNA extraction and bisulfite conversion
Both OCT or FFPE TNBC tissue were used depending on availability and DNA extraction was performed using DNA purification protocol with NaCl or the kit “DNA PPPE QIAamp Tissue” (Qiagen) respectively. DNA from cells lines was extracted using Trizol reagent. DNA samples were quantified using PicoGreen method (Invitrogen) and quality was evaluated using Nanodrop (Thermo Scientific) and electrophoresis gels. For DNA bisulfite modification, 500 ng of DNA were used and modification was performed with EZ-96 DNA Methylation (Zymo Research Corp.)
DNA methylation analysis by Infinium Human Methylation 450 BeadChip array
Microarray-based DNA methylation analysis was conducted with the Infinium Human Methylation 450 BeadChip (450K array; Illumina, San Diego, CA), that covers > 450,000 CpG sites along the human genome . After bisulfite conversion, hybridization was performed following the Illumina Infinium HD methylation protocol.
Methylation score of each CpG was represented as β value that ranged between 0 (unmethylated) and 1 (completely methylated). Color balance adjustment and normalization were performed using GenomeStudio Illumina software (V2010.3). After filtering, differentially methylated CpG sites (DMCpGs) between R and NR groups were identified following this flowchart (Fig. 1a): for each probe/CpG, the sets of methylation β values belonging to both groups were compared to obtain (1) DMCpGs with a significant p value < 0.05 (1030 CpGs); (2) DMCpGs with average β values between R and NR groups ≥ 0.20 (133 out of 1030 CpGs); (3) DMCpGs localized in island or shore regions of promotors (35 out of 133 CpGs); (4) DMCpGs with intragroup standard deviation (SD) ≤ 0.20 in order to select the most relevant positions for validation (11 out of 35 CpGs); (5) DMCpGs with a consistent methylation profile in consecutive CpGs (10 out of 11). This final filter with ten CpGs (9 genes) yielded the best candidates for validation.
DNA methylation analysis by bisulfite pyrosequencing
Quantitative DNA methylation analysis was performed by bisulfite pyrosequencing of consecutive cytosines located in islands or shores of promoter regions of candidate genes using a Pyro Gold SQA™ Reagent Kit (Qiagen) in a PyroMark Q96 System version 2.0.6 (Qiagen) according to the manufacturer’s instructions. CpG site methylation quantification was obtained using Pyro Q-CpG 1.0.9 (Qiagen). Primer sequences (Additional file 8: Table S7) were designed with PyroMark Assay Design 2.0 (Qiagen).
Gene expression studies and correlation with methylation levels
Total RNA was isolated from OCT/ FFPE samples by mirVana Isolation Kit (Ambion) and from cell lines using Trizol (Invitrogen) according to the manufacturer’s protocol. The RNA (500 ng) were retrotranscribed using the High-Capacity cDNA Reverse Transcription kit (Applied biosystems) according to the manufacturer. Quantitative RT-PCR (qRT-PCR) reactions were performed in triplicate on an Applied Biosystems 7900HT Fast Real-Time PCR system using TaqMan expression assays (FERD3L Hs00541737_s1; TRIP10 Hs00182848_m1; GAPDH Hs03929097_g1 (Applied Biosystems). Gene expression was assayed with GAPDH as endogenous control and using the delta delta Ct method.
Gene ontology and The Cancer Genome Atlas (TCGA) database analysis
A gene ontology (GO) analysis was performed to estimate the enrichment of the DMCpGs identified in particular biological processes . This analysis detects the significant over-representation of GO terms in one of the sets with respect to the other for the entire genome. GO terms with p value < 0.05 were considered significant. DNA methylation and expression data from patients with invasive breast carcinoma were obtained from The Cancer Genome Atlas (TCGA) using the MethHc database (http://methhc.mbc.nctu.edu.tw/php/index.php). Paired DNA methylation and expression data from 713 patients obtained from Infinium 450K array and RNA-Seq, respectively, were used. Based in the methylation results of our study, we considered the methylation data of FERD3L promoter region and the methylation of an individual CpG (cg10043037) located at the promoter region.
In the DC, data were summarized by mean, SD, or median. To identify consistent patterns of differentially methylated CpG sites between responders (R) versus non-responders (NR), a non-parametric Wilcoxon rank sum test was performed in the DC. This test demonstrates quite robust results even for a small number of subjects. Globally, a two-tailed p value of less than 0.05 was considered to indicate statistical significance. All statistical analyses were performed using GraphPad Prism 7 and R software (version 3.2.0).
In the VC, differences in DNA methylation and transcript levels of the identified genes between R and NR and gene expression changes in the cell lines after the demethylating treatment were assessed by the non-parametric Mann-Whitney U test. The correlation between methylation and transcript levels was assessed by Spearman’s rank correlation coefficient.
Receiver operating characteristic (ROC) curves were used to assess the diagnostic predictive capacity of the candidate biomarkers. The area under the curve (AUC) was computed for each ROC curve, and 95% confidence intervals (CI) were also estimated by bootstrapping with 1000 iterations. Sensitivity and specificity were estimated at the optimal cut-off point according to Youden criteria.
Data were summarized by mean, SD, or median. Globally, a two-tailed p value of less than 0.05 was considered to indicate statistical significance. All statistical analyses were performed using GraphPad Prism 7 and R software (version 3.2.0).
Akaike information criterion
American Society of Clinical Oncology/College of American Pathologists
American Type Culture Collection
Area under the curve
Cyclin-dependent kinase like 2
Cell adhesion molecule L1 like
Cytosine and guanine separated by a phosphate
Differentially methylated CpG sites
E2F transcription factor 4
EvC ciliary complex subunit 2
Fer3-like bHLH transcription factor
Erb-b2 receptor tyrosine kinase 2
Luminal androgen receptor
Lymphoid enhancer binding factor 1
A database of DNA Methylation and gene expression in Human Cancer
Non-pathological complete response
Pathological complete response
Quantitative reverse transcriptase polymerase chain reaction
Residual cancer burden
RNA sequencing (RNA: ribonucleic acid)
Receiver operating characteristic
Reverse transcriptase polymerase chain reaction
Small interfering RNA
The cancer genome atlas
T cell leukemia homeobox 3
Triple-negative breast cancer
Thyroid hormone receptor interactor 10
Twist family bHLH transcription factor 1
Sousa B, Cardoso F. Neoadjuvant treatment for HER-2-positive and triple-negative breast cancers. Ann Oncol. 2012;23(Suppl 10):x237–42.
Carey LA, Dees EC, Sawyer L, Gatti L, Moore DT, Collichio F, et al. The triple negative paradox: primary tumor chemosensitivity of breast cancer subtypes. Clin Cancer Res. 2007;13(8):2329–34.
Biswas T, Efird JT, Prasad S, Jindal C, Walker PR. The survival benefit of neoadjuvant chemotherapy and pCR among patients with advanced stage triple negative breast cancer. Oncotarget. 2017;8(68):112712–9.
Lee HJ, Park IA, Song IH, Kim SB, Jung KH, Ahn JH, et al. Comparison of pathologic response evaluation systems after anthracycline with/without taxane-based neoadjuvant chemotherapy among different subtypes of breast cancers. PLoS One. 2015;10(9):e0137885.
Sikov WM, Berry DA, Perou CM, Singh B, Cirrincione CT, Tolaney SM, et al. Impact of the addition of carboplatin and/or bevacizumab to neoadjuvant once-per-week paclitaxel followed by dose-dense doxorubicin and cyclophosphamide on pathologic complete response rates in stage II to III triple-negative breast cancer: CALGB 40603 (Alliance). J Clin Oncol. 2015;33(1):13–21.
Gamucci T, Pizzuti L, Sperduti I, Mentuccia L, Vaccaro A, Moscetti L, et al. Neoadjuvant chemotherapy in triple-negative breast cancer: a multicentric retrospective observational study in real-life setting. J Cell Physiol. 2018;233(3):2313–23.
Penault-Llorca F, Viale G. Pathological and molecular diagnosis of triple-negative breast cancer: a clinical perspective. Ann Oncol. 2012;23(Suppl 6):vi19–22.
Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, et al. A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004;351(27):2817–26.
van ‘t Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415(6871):530–6.
Whitworth P, Stork-Sloots L, de Snoo FA, Richards P, Rotkis M, Beatty J, et al. Chemosensitivity predicted by BluePrint 80-gene functional subtype and MammaPrint in the Prospective Neoadjuvant Breast Registry Symphony Trial (NBRST). Ann Surg Oncol. 2014;21(10):3261–7.
Bertucci F, Finetti P, Viens P, Birnbaum D. EndoPredict predicts for the response to neoadjuvant chemotherapy in ER-positive, HER2-negative breast cancer. Cancer Lett. 2014;355(1):70–5.
Prat A, Galvan P, Jimenez B, Buckingham W, Jeiranian HA, Schaper C, et al. Prediction of response to neoadjuvant chemotherapy using core needle biopsy samples with the Prosigna assay. Clin Cancer Res. 2016;22(3):560–6.
Mark KMK, Varn FS, Ung MH, Qian F, Cheng C. The E2F4 prognostic signature predicts pathological response to neoadjuvant chemotherapy in breast cancer patients. BMC Cancer. 2017;17(1):306.
Echavarria I, Lopez-Tarruella S, Picornell A, Garcia-Saenz JA, Jerez Y, Hoadley K, et al. Pathological response in a triple-negative breast cancer cohort treated with neoadjuvant carboplatin and docetaxel according to Lehmann's refined classification. Clin Cancer Res. 2018;24(8):1845–52.
Wolf DM, Yau C, Sanil A, Glas A, Petricoin E, Wulfkuhle J, et al. DNA repair deficiency biomarkers and the 70-gene ultra-high risk signature as predictors of veliparib/carboplatin response in the I-SPY 2 breast cancer trial. NPJ Breast Cancer. 2017;3:31.
Aznar MA, Labiano S, Diaz-Lagares A, Molina C, Garasa S, Azpilikueta A, et al. CD137 (4-1BB) Costimulation modifies DNA methylation in CD8(+) T cell-relevant genes. Cancer Immunol Res. 2018;6(1):69–78.
Diaz-Lagares A, Crujeiras AB, Lopez-Serra P, Soler M, Setien F, Goyal A, et al. Epigenetic inactivation of the p53-induced long noncoding RNA TP53 target 1 in human cancer. Proc Natl Acad Sci U S A. 2016;113(47):E7535–E44.
Stefansson OA, Esteller M. Epigenetic modifications in breast cancer and their role in personalized medicine. Am J Pathol. 2013;183(4):1052–63.
Symmans WF, Peintinger F, Hatzis C, Rajan R, Kuerer H, Valero V, et al. Measurement of residual breast cancer burden to predict survival after neoadjuvant chemotherapy. J Clin Oncol. 2007;25(28):4414–22.
Chekhun VF, Kulik GI, Yurchenko OV, Tryndyak VP, Todor IN, Luniv LS, et al. Role of DNA hypomethylation in the development of the resistance to doxorubicin in human MCF-7 breast adenocarcinoma cells. Cancer Lett. 2006;231(1):87–93.
Liu X, Liu X, Wu Y, Fang Z, Wu Q, Wu C, et al. MicroRNA-34a attenuates metastasis and chemoresistance of bladder cancer cells by targeting the TCF1/LEF1 axis. Cell Physiol Biochem. 2018;48(1):87–98.
Liu X, Luo X, Wu Y, Xia D, Chen W, Fang Z, et al. MicroRNA-34a attenuates paclitaxel resistance in prostate cancer cells via direct suppression of JAG1/Notch1 axis. Cell Physiol Biochem. 2018;50(1):261–76.
Zhang B, Zhang Y, Dagher MC, Shacter E. Rho GDP dissociation inhibitor protects cancer cells against drug-induced apoptosis. Cancer Res. 2005;65(14):6054–62.
Cerami EG, Gross BE, Demir E, Rodchenkov I, Babur O, Anwar N, et al. Pathway commons, a web resource for biological pathway data. Nucleic Acids Res. 2011;39(Database issue):D685–90.
Yang M, Li Y, Shen X, Ruan Y, Lu Y, Jin X, et al. CLDN6 promotes chemoresistance through GSTP1 in human breast cancer. J Exp Clin Cancer Res. 2017;36(1):157.
von Minckwitz G, Untch M, Blohmer JU, Costa SD, Eidtmann H, Fasching PA, et al. Definition and impact of pathologic complete response on prognosis after neoadjuvant chemotherapy in various intrinsic breast cancer subtypes. J Clin Oncol. 2012;30(15):1796–804.
Ueno T, Masuda N, Yamanaka T, Saji S, Kuroi K, Sato N, et al. Evaluating the 21-gene assay recurrence score (R) as a predictor of clinical response to 24 weeks of neoadjuvant exemestane in estrogen receptor-positive breast cancer. Int J Clin Oncol. 2014;19(4):607–13.
Stirzaker C, Zotenko E, Song JZ, Qu W, Nair SS, Locke WJ, et al. Methylome sequencing in triple-negative breast cancer reveals distinct methylation clusters with prognostic value. Nat Commun. 2015;6:5899.
Mathe A, Wong-Brown M, Locke WJ, Stirzaker C, Braye SG, Forbes JF, et al. DNA methylation profile of triple negative breast cancer-specific genes comparing lymph node positive patients to lymph node negative patients. Sci Rep. 2016;6:33435.
Verzi MP, Anderson JP, Dodou E, Kelly KK, Greene SB, North BJ, et al. N-twist, an evolutionarily conserved bHLH protein expressed in the developing CNS, functions as a transcriptional inhibitor. Dev Biol. 2002;249(1):174–90.
Glackin CA. Targeting the twist and Wnt signaling pathways in metastatic breast cancer. Maturitas. 2014;79(1):48–51.
Li S, Kendall SE, Raices R, Finlay J, Covarrubias M, Liu Z, et al. TWIST1 associates with NF-kappaB subunit RELA via carboxyl-terminal WR domain to promote cell autonomous invasion through IL8 production. BMC Biol. 2012;10:73.
Maestro R, Dei Tos AP, Hamamori Y, Krasnokutsky S, Sartorelli V, Kedes L, et al. Twist is a potential oncogene that inhibits apoptosis. Genes Dev. 1999;13(17):2207–17.
Vesuna F, Lisok A, Kimble B, Raman V. Twist modulates breast cancer stem cells by transcriptional regulation of CD24 expression. Neoplasia. 2009;11(12):1318–28.
Vesuna F, Lisok A, Kimble B, Domek J, Kato Y, van der Groep P, et al. Twist contributes to hormone resistance in breast cancer by downregulating estrogen receptor-alpha. Oncogene. 2012;31(27):3223–34.
Yang J, Mani SA, Donaher JL, Ramaswamy S, Itzykson RA, Come C, et al. Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell. 2004;117(7):927–39.
Wang WS, Yang XS, Xia M, Jiang HY, Hou JQ. Silencing of twist expression by RNA interference suppresses epithelial-mesenchymal transition, invasion, and metastasis of ovarian cancer. Asian Pac J Cancer Prev. 2012;13(9):4435–9.
Yang J, Mani SA, Weinberg RA. Exploring a new twist on tumor metastasis. Cancer Res. 2006;66(9):4549–52.
Abe M, Watanabe N, McDonell N, Takato T, Ohira M, Nakagawara A, et al. Identification of genes targeted by CpG island methylator phenotype in neuroblastomas, and their possible integrative involvement in poor prognosis. Oncology. 2008;74(1–2):50–60.
Gururaj AE, Rayala SK. Kumar R. p21-activated kinase signaling in breast cancer. Breast Cancer Res. 2005;7(1):5–12.
Weaver BA. How Taxol/paclitaxel kills cancer cells. Mol Biol Cell. 2014;25(18):2677–81.
Dogterom M, Koenderink GH. Actin-microtubule crosstalk in cell biology. Nat Rev Mol Cell Biol. 2019;20(1):38–54.
Rolland Y, Marighetti P, Malinverno C, Confalonieri S, Luise C, Ducano N, et al. The CDC42-interacting protein 4 controls epithelial cell cohesion and tumor dissemination. Dev Cell. 2014;30(5):553–68.
Pichot CS, Arvanitis C, Hartig SM, Jensen SA, Bechill J, Marzouk S, et al. Cdc42-interacting protein 4 promotes breast cancer cell invasion and formation of invadopodia through activation of N-WASp. Cancer Res. 2010;70(21):8347–56.
Cerqueira OL, Truesdell P, Baldassarre T, Vilella-Arias SA, Watt K, Meens J, et al. CIP4 promotes metastasis in triple-negative breast cancer and is associated with poor patient prognosis. Oncotarget. 2015;6(11):9397–408.
Hsiao SH, Lee KD, Hsu CC, Tseng MJ, Jin VX, Sun WS, et al. DNA methylation of the Trip10 promoter accelerates mesenchymal stem cell lineage determination. Biochem Biophys Res Commun. 2010;400(3):305–12.
Hsu CC, Leu YW, Tseng MJ, Lee KD, Kuo TY, Yen JY, et al. Functional characterization of Trip10 in cancer cell growth and survival. J Biomed Sci. 2011;18:12.
Leu YW, Yan PS, Fan M, Jin VX, Liu JC, Curran EM, et al. Loss of estrogen receptor signaling triggers epigenetic silencing of downstream targets in breast cancer. Cancer Res. 2004;64(22):8184–92.
Jatkoe TA, Karnes RJ, Freedland SJ, Wang Y, Le A, Baden J. A urine-based methylation signature for risk stratification within low-risk prostate cancer. Br J Cancer. 2015;112(5):802–8.
Meseure D, Vacher S, Alsibai KD, Nicolas A, Chemlali W, Caly M, et al. Expression of ANRIL-polycomb complexes-CDKN2A/B/ARF genes in breast tumors: identification of a two-gene (EZH2/CBX7) signature with independent prognostic value. Mol Cancer Res. 2016;14(7):623–33.
Sun C, Cheng X, Wang C, Wang X, Xia B, Zhang Y. Gene expression profiles analysis identifies a novel two-gene signature to predict overall survival in diffuse large B-cell lymphoma. Biosci Rep. 2019;39(1). https://doi.org/10.1042/BSR20181293.
Sun Y, Hou L, Yang Y, Xie H, Yang Y, Li Z, et al. Two-gene signature improves the discriminatory power of IASLC/ATS/ERS classification to predict the survival of patients with early-stage lung adenocarcinoma. Onco Targets Ther. 2016;9:4583–91.
Wang X, Lu Y, Feng W, Chen Q, Guo H, Sun X, et al. A two kinase-gene signature model using CDK2 and PAK4 expression predicts poor outcome in non-small cell lung cancers. Neoplasma. 2016;63(2):322–9.
Xu Q, Tan C, Ni S, Wang Q, Wu F, Liu F, et al. Identification and validation of a two-gene expression index for subtype classification and prognosis in diffuse large B-cell lymphoma. Sci Rep. 2015;5:10006.
Xu X, Li M, Hu J, Chen Z, Yu J, Dong Y, et al. Expression profile analysis identifies a two-gene signature for prediction of head and neck squamous cell carcinoma patient survival. J Cancer Res Ther. 2018;14(7):1525–34.
Zhao F, Olkhov-Mitsel E, Kamdar S, Jeyapala R, Garcia J, Hurst R, et al. A urine-based DNA methylation assay, ProCUrE, to identify clinically significant prostate cancer. Clin Epigenetics. 2018;10(1):147.
Wolff AC, Hammond ME, Hicks DG, Dowsett M, McShane LM, Allison KH, et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists clinical practice guideline update. Arch Pathol Lab Med. 2014;138(2):241–56.
Sandoval J, Heyn H, Moran S, Serra-Musach J, Pujana MA, Bibikova M, et al. Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics. 2011;6(6):692–702.
Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A. GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 2012;40(Web Server issue):W478–83.
The authors thank sample donors and workers who supported the data collection.
This work was supported by the project GV037/2016, from Consellería D′Educació, Investigació, Cultura i Esport, and to CIBERONC (CB16/12/00481), PI15/01617 and PI18/01219 from ISCIII and FEDER founds. ADL. is funded by a contract “Juan Rodés” (JR17/00016) from ISCIII. ABC. is funded by a contract “Miguel Servet” (CP17/00088) from ISCIII.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study was approved by the local clinical Research Ethics Committee of the Hospital Clínico Universitario de Valencia /INCLIVA with the reference number 2014/178. Patients were requested to sign informed consent for tumor tissue molecular analysis and biobanking.
Consent for publication
Authors BP, ADL, JA PF, ME, ALL, and PE are currently applying for a patent (application number EP18382390.5) relating the content of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of all the biological processes enriched for the 71 differentially methylated genes between responder and non-responder patients according to the Gene Ontology analysis (DOCX 32 kb)
Thirty-five differentially methylated CpGs between responders and non-responders group selected from 450k array (delta value ≥ 0.2) corresponding to 23 genes located in promoter and island/shore (PPT 172 kb)
Eleven differentially methylated CpGs, corresponding to 11 genes, showed significant methylation differences between non-responder and responder patients: 6 genes (LOC641518; LEF1; HOXA5; EVC2; CDKL2; TLX3) presented a methylation increase in non-responders group vs responders, and 5 genes (ZFHX4; LOC100192378; FERD3L; CHL1; TRIP10) decreased methylation level in non-responder patients compared to those who responded to NAC treatment (PPT 225 kb)
CpGs studied by pyrosequencing in the DC and in the VC to validate methylation in the candidate genes identified in the 450k array (Illumina). In bold, CpGs from 450k array. Normal type, consecutive CpGs (PPT 140 kb)
Representation of the pathway interaction network of FERD3L and TRIP10 with other genes using Pathway Commons. FERD3L and TRIP10 are able to interact with different genes that have shown to be implicated in cancer drug resistance (PPT 452 kb)
Mean differences in methylation levels according to clinicopathological prognostic factors in both cohorts (DC+VC). cT, clinical tumor size; cN, clinical nodule affectation (PPTX 48 kb)
Clinical inclusion and exclusion criteria followed to select TNBC patients for the methylation study (PPT 89 kb)
Sequence of primers used by pyrosequencing in the validation assay of candidate genes obtained from 450k array (PPT 143 kb)
About this article
Cite this article
Pineda, B., Diaz-Lagares, A., Pérez-Fidalgo, J.A. et al. A two-gene epigenetic signature for the prediction of response to neoadjuvant chemotherapy in triple-negative breast cancer patients. Clin Epigenet 11, 33 (2019). https://doi.org/10.1186/s13148-019-0626-0
- Triple-negative breast cancer
- Epigenetic signature