Skip to main content

Epigenetic silencing of TGFBI confers resistance to trastuzumab in human breast cancer



Acquired resistance to trastuzumab is a major clinical problem in the treatment of HER2-positive (HER2+) breast cancer patients. The selection of trastuzumab-resistant patients is a great challenge of precision oncology. The aim of this study was to identify novel epigenetic biomarkers associated to trastuzumab resistance in HER2+ BC patients.


We performed a genome-wide DNA methylation (450K array) and a transcriptomic analysis (RNA-Seq) comparing trastuzumab-sensitive (SK) and trastuzumab-resistant (SKTR) HER2+ human breast cancer cell models. The methylation and expression levels of candidate genes were validated by bisulfite pyrosequencing and qRT-PCR, respectively. Functional assays were conducted in the SK and SKTR models by gene silencing and overexpression. Methylation analysis in 24 HER2+ human BC samples with complete response or non-response to trastuzumab-based treatment was conducted by bisulfite pyrosequencing.


Epigenomic and transcriptomic analysis revealed the consistent hypermethylation and downregulation of TGFBI, CXCL2, and SLC38A1 genes in association with trastuzumab resistance. The DNA methylation and expression levels of these genes were validated in both sensitive and resistant models analyzed. Of the genes, TGFBI presented the highest hypermethylation-associated silencing both at the transcriptional and protein level. Ectopic expression of TGFBI in the SKTR model suggest an increased sensitivity to trastuzumab treatment. In primary tumors, TGFBI hypermethylation was significantly associated with trastuzumab resistance in HER2+ breast cancer patients.


Our results suggest for the first time an association between the epigenetic silencing of TGFBI by DNA methylation and trastuzumab resistance in HER2+ cell models. These results provide the basis for further clinical studies to validate the hypermethylation of TGFBI promoter as a biomarker of trastuzumab resistance in HER2+ breast cancer patients.


Breast cancer (BC) is the most common cancer in women worldwide and the leading cause of cancer deaths for women [1]. Approximately 15–20% of patients with this tumor overexpress the human epidermal growth factor receptor 2 (HER2) protein [2]. HER2 can activate downstream signaling cascades that induce cell proliferation through the Ras-mitogen-activated protein kinase (MAPK) pathway and inhibit cell death through the phosphatidylinositol 3′-kinase (PI3K)/protein kinase B (Akt)/mammalian target of the rapamycin (mTOR) pathway [3, 4]. Characterizing HER2 as a proto-oncogene, a poor prognostic marker, and eventually as a therapeutic target has dramatically changed BC classification, risk assessment, and treatment [5]. The development of powerful targeted therapies directed specifically at HER2 has improved the survival of patients with both early-stage and metastatic BC.

Trastuzumab (Herceptin®; Genentech, Inc., South San Francisco, CA, USA) is a humanized monoclonal antibody that selectively binds with high affinity to the extracellular domain of human HER2 protein and was the first targeted drug approved for treating HER2+ BC [6]. Although there are currently other anti-HER2 agents available, trastuzumab remains the gold standard for treating this subtype of BC. Despite the clinical benefits trastuzumab treatment in HER2+ BC brings, a large percentage of patients display primary or acquired resistance to the drug. In the last few years, several studies have focused on identifying the molecular mechanisms of trastuzumab resistance, such as aberrant activation of downstream signaling pathways [7] or the HER2 carboxyl-terminal fragments (CTF), also known as p95HER2, which are frequently found in HER2-expressing BC cell lines and tumors [8]. Despite these efforts, the complete picture of the molecular mechanisms triggering trastuzumab resistance in BC remains unclear. Therefore, clarifying these mechanisms and identifying new resistance biomarkers is essential in the advance towards precision oncology in BC and the quest for new treatment options for those patients who do not respond to trastuzumab therapy.

DNA methylation is the most well-known epigenetic modification in humans and has been implicated in regulating the expression of a great variety of critical genes in cancer [9]. For this reason, DNA methylation status has emerged as one of the most promising epigenetic biomarkers for several types of cancer, including BC [10]. These epigenetic markers can be useful in detecting tumors earlier or identifying patients with an increased risk of cancer, as well as evaluating disease progression or predicting the response to anticancer drugs [11]. In the last few years, some significant genes that are inactivated by promoter methylation in BC have been identified, including BRCA1 [12] and RASSF1A [13]. However, the analysis of hypermethylated genes in association with trastuzumab resistance is still a largely unexplored field that holds great potential.

The aim of this study was to evaluate the implication DNA methylation has in trastuzumab resistance and to identify epigenetically regulated genes with potential clinical value as biomarkers for trastuzumab resistance in HER2+ BC patients. To this purpose, we employed an integrative approach with genome-wide DNA methylation (450K array) and a transcriptomic analysis (RNA-Seq) in trastuzumab-sensitive and trastuzumab-resistant cell line models (SK and SKTR). In vitro results were validated in a cohort of 24 HER2+ tumor samples from patients experiencing both sensitivity and resistance to trastuzumab-based neoadjuvant therapy settings. Here, we determine epigenetic inactivation of the TGFBI gene by promoter CpG island hypermethylation, which are CpG-rich regions of DNA that are often associated with the transcription start sites of genes, with possible implications for trastuzumab-resistant BC pathways [14]. The hypermethylation of TGFBI also suggest its potential clinical usefulness as a biomarker for trastuzumab resistance in HER2+ BC patients.


Cell culture

SKBr3 (SK) and AU565 (AU) HER2+ breast carcinoma cells were obtained from Eucellbank (University of Barcelona, Spain) [15] and the American Type Culture Collection (ATCC, Rockville, MD, USA). SKBr3 and AU565 cells were routinely grown in McCoy’s (Gibco) and Dulbecco’s modified Eagle’s medium (DMEM; Gibco), respectively, and supplemented with 10% FBS (HyClone Laboratories), 1% l-glutamine (Gibco), 1% sodium pyruvate (Gibco), and 100 U/mL penicillin/streptomycin (HyClone Laboratories). Cell lines were kept at 37 °C and 5% CO2 atmosphere. Long-term trastuzumab-resistant SK cells (SKTR) and AU565 cells (AUTR) had been previously developed in our laboratory [16, 17]. Resistance was confirmed with cell viability assays. The trastuzumab-resistant SKTR and AUTR cells were maintained in 2 μM of trastuzumab, i.e., a concentration in which parental cells were not viable.

Patients and tissue samples

TGFBI promoter methylation levels were retrospectively evaluated in tumor samples from 24 patients diagnosed with HER2+ BC at the Dr. Josep Trueta University Hospital, Girona (Spain) between 2007 and 2015. The patients were selected from the hospital’s pharmacy database. The selection criterion included patients with early or locally advanced HER2+ BC who had received neoadjuvant treatment with trastuzumab and chemotherapy. Twenty patients had no response or partial response and 4 patients had complete response to trastuzumab plus chemotherapy. For all patients, hematoxylin and eosin (H&E)-stained slides from formalin-fixed paraffin-embedded (FFPE) tumor blocks were examined to determine the representative areas of the invasive tumor. Estrogen receptor (ER), progesterone receptor (PR), and HER2 expression had been previously analyzed in the tumors using immunohistochemistry (IHC). For each patient, clinical and histopathological features were obtained: age, stage (TNM classification [18]), histological grade (Bloom-Richardson grading system), menopause status, type of surgery, and relapse.

5-Aza-2′-deoxycytidine treatment

Epigenetic signatures are characterized by a very dynamic nature, where DNA methylation has often been shown as a reversible mechanism of transcriptional control by inhibition of enzymes such as the DNA methyltransferases [19]. We performed reactivation treatments using the demethylating agent, 5-aza-2′-deoxycytidine (5-aza-dC, Sigma-Aldrich, St Louis, MO) at 3 μM and 5 μM of 5-aza-dC for 72 h. The medium was changed every day to promote DNA demethylation.

DNA and RNA isolation procedures

Genomic DNA extraction from cell lines or formalin-fixed paraffin-embedded (FFPE) core biopsies (10 μm) and tissue sections (5 μm) using a QIAamp DNA Mini Kit and Deparaffinization Solution with a QIAamp DNA FFPE Tissue Kit (Qiagen, Hilden, Germany), respectively, were carried out following the manufacturer’s instructions. For the RT-PCR experiments, cells were washed with PBS and then suspended in 1 mL of Qiazol (Qiagen Hilden, Germany) was added. Total RNA was isolated using a GeneJET RNA Purification Kit (Thermo Fisher Scientific) following the instructions provided by the manufacturer. All DNA and RNA samples were quantified using a NanoDrop 2000 Spectrophotometer (Thermo Fisher Scientific).

DNA was bisulfite-modified using an EZ DNA Methylation-Gold Kit (Zymo Research) in accordance with the manufacturer’s recommendations.

DNA methylation array

Genome-wide DNA methylation analysis was performed using an Illumina 450K DNA methylation microarray (Infinium HumanMethylation450 BeadChip) as previously described [20]. This technique is an epigenomic approach that allows analyzing the methylation profile of the human genome. Bisulfite-converted DNA from SK and SKTR models were used to hybridized on an Illumina Infinium HumanMethylation450 (450K) BeadChip array, following the Illumina Infinium HD Methylation protocol. Sex chromosomes were excluded from the analysis because they usually represent a high source of variation in DNA methylation levels [21]. Data was deposited into the NCBI Gene Omnibus, accession number: GSE123754.

RNA sequencing

DNA methylation is an epigenetic mechanism that usually affects CpG islands and promoters regulating the expression of the genes. For this reason, a transcriptomic analysis by RNA-Seq was performed. Library construction was performed using an Illumina TruSeq Stranded mRNA Sample Preparation Kit (CAT. No. RS-122-2101, RS-122-2102). Following the transfer of the flow cell to an Illumina HiSeq instrument, a 101-cycle paired-end sequence run was performed using the HiSeq SBS Kit v4 sequencing reagents (FC-401-4002). A detailed description of the analysis is provided in the Additional file 1. Sequencing data have been posted in the Gene Expression Omnibus (GEO) database under accession number GSE114575.

Bisulfite pyrosequencing

To validate the results obtained from the arrays, bisulfite pyrosequencing analyses were performed. The primers for the PCR amplification and sequencing were designed with the PyroMark assay design software version DNA-bs was amplified by PCR under standard conditions with biotinylated primers (see Additional file 1: Table S1) to convert the PCR product into single-stranded DNA templates. We used a Vacuum Prep Tool (Biotage, Sweden) to prepare single-stranded PCR products following the manufacturer’s instructions. The PCR products were observed in 2% agarose gels before pyrosequencing. Reactions were performed in a PyroMark Q96 System version 2.0.6 (Qiagen, Hilden, Germany) using appropriate reagents and protocols. The methylation value was obtained from the average of the CpG dinucleotides included in the sequences analyzed using a Pyro Q-CpG 1.0.9 (Qiagen, Hilden, Germany).

Methylation-specific PCR (MSP)

DNA methylation was also validated using the methylation-specific PCR (MSP) using a set of primers designed by the Methyl Primer Express program for each gene analyzed. For each sample, two PCR reactions were produced using the same DNA template: one reaction with methylated primers (M: methylated) and the other with unmethylated primers (U: unmethylated). Commercial methylated human male genomic DNA (CpGenome Universal Methylated DNA, Millipore) was used as a methylated positive control (IVD), and DNA from normal lymphocytes (NL) as a positive control for unmethylated alleles. Primer sequences can be found in Additional file 1: Table S2.

Gene expression analysis (qRT-PCR)

To confirm the differential regulation of TGFBI, CXCL2 and SLC38A1 observed in SK and SKTR models used gold standard methodologies for gene expression such as quantitative real-time PCR. Total RNA was reverse-transcribed into cDNA using a High Capacity cDNA Archive Kit (Applied Biosystems). Gene expression levels of selected genes were assessed using a LightCycler 480 Real-time PCR System (Roche) with a LightCycler 480 SYBR Green I Master (Roche). Primers are depicted in Additional file 1: Table S3. RT-PCR analyses were performed at least four times, and each gene was run in triplicate. GAPDH was used as an endogenous control to enable normalization.

Short hairpin interference and ectopic expression assays

For the TGFBI long-term knockdown, two different ShRNAs were specifically designed against TGFBI mRNA (NM_000358) in two different loci, by considering a 19-base target sequence and inserted in the pLVX-shRNA2 vector (Clontech). For stable overexpression, the cDNA sequence of TGFBI was amplified from SKWT cDNA and cloned into the pLVX-IRES-tdTomato vector (Clontech). This construction (wild type) was then used as a template to create its mutated version following a PCR-based strategy. Mutations were performed in the RGD domain and the NKDIL (amino acids 354–358), YH18 (amino acids 563–580), and EPDIM (amino acids 617–621) motif in the second and fourth FAS-1 domains [22]. Lentiviral production was performed using the protocol described in Additional file 1. After lentiviral transduction, ZsGreen1-positive or tdTomato-positive cells were sorted by flow cytometry.

Cell viability assays

Cell viability was determined using a 3-(4, 5-dimethyl-2-thiazolyl)-2, 5-diphenyl-2 h-tetrazolium bromide (MTT) assay. Briefly, cells were plated in their growth medium in 96-well plates at a cell density of 1.5 × 103 cells per well. After 24 h, the growth medium was removed, and, 100 μL of fresh medium containing the corresponding concentration of trastuzumab was added to each well for 5 days. Cell viability was measured using the standard colorimetric MTT assay as previously described [16]. Using the multi-well plate reader Benchmark Plus (Bio-Rad), absorbance was determined to be 570 nm. Data presented are from three separate wells per assay, and the assay was performed at least three times.

Western blot analysis

We previously demonstrated that different activation patterns of some HER receptors and their downstream signaling proteins are associated with the molecular mechanisms of acquired trastuzumab resistance (SK and SKTR) [16]. For this reason, a panel of the different proteins were analyzed in our models using Western blot analysis. The parental (SK and AU) and resistant (SKTR and AUTR) models were synchronized by starvation in serum-deprived medium (0.5% FBS) for 24 h. Cells were lysed in ice-cold lysis buffer (Cell Signaling Technology, Inc.) with 100 μg/mL phenylmethylsulfonylfluoride (PMSF). Protein concentration was determined with Lowry (DC Protein Assay, Bio-Rad). Equal amounts of protein were heated in LDS Sample Buffer with Sample Reducing Agent (Invitrogen) for 10 min at 70 °C, separated on SDS-PAGE and transferred to nitrocellulose membranes. Protein was detected using primary antibodies (Additional file 1: Table S5). Secondary antibodies, α-tubulin and β-actin, conjugated to horseradish peroxidase were used. The immune complexes were detected using a chemiluminescent HRP substrate [SuperSignal™ West Femto Maximum Sensitivity Substrate (Thermo Scientific™ Inc.) or Immobilon Western HRP Substrate (Millipore Sigma)] and in a Bio-Rad ChemiDoc™ MP Imaging System. Western blot analyses were repeated at least three times and representative results are shown.

Statistical analysis

Data were analyzed with Student’s t test when two groups were being compared or with one-way analysis of variance (ANOVA) followed by Tukey’s honest significant difference (HSD) post hoc test or Tamhane’s T2 post hoc test for multiple comparisons. The non-parametric Kruskal-Wallis test was used when data did not follow normal distribution. When two groups were being compared, the data were analyzed with the Mann-Whitney U tests for non-normally independent variables; otherwise, the Kruskal-Wallis test was used for more than two groups. In the patient cohort, we analyzed the TGFBI promoter methylation status and its association with the clinical-histopathological characteristics. Patient data were summarized as median (first quartile-third quartile) for continuous variables, and frequencies and percentages for categorical variables. The potential association between clinical-histopathological characteristics and levels of TGFBI methylation (low ≤ 20%, high ≥ 20%) or differences in TGFBI methylation before and after treatment were analyzed using the Mann-Whitney U and Kruskal-Wallis tests for continuous variables and the chi-square or Fisher exact tests, as appropriate, for categorical variables. The correlation between variables was observed using Spearman’s rho coefficient. Characteristic (ROC) curves were used to assess the predictive capacity of the TGFBI marker. ROC curve analysis was performed using MedCalc software Version 18.11.3. The Hosmer-Lemeshow test and calibration plot was assessed with MedCalc and STATA software. Levels of significance were set at p < 0.05 and are represented by asterisks, as follows: p < 0.05 (denoted as *), p < 0.01 (denoted as **), and p < 0.001 (denoted as ***). The statistical analysis was performed using the IBM SPSS software (Version 21.0; SPSS Inc.)


Genome-wide DNA methylation analysis in trastuzumab-sensitive (SK) and trastuzumab-resistant (SKTR) breast cancer models

The Infinium HumanMethylation 450 BeadChip (450k array) is a validated tool for carrying out epigenomic projects because it allows the methylation status of approximately 450,000 CpGs located throughout the human genome to be detected [20]. Our group recently employed this approach in previous studies aimed at identifying genes regulated by DNA methylation for clinical applications such as biomarkers [23, 24]. Therefore, we took advantage of the 450K array methodology to characterize the DNA methylation profile associated with trastuzumab resistance in BC, and compared a trastuzumab-sensitive (SK) and trastuzumab-resistant (SKTR) BC models (Fig. 1a). The global analysis of the DNA methylation (β values) corresponding to the CpG sites with p value < 0.01 (469,927 CpGs) showed differences in the scatter plot (Fig. 1b) between SK and SKTR models (r2 = 0.93) and revealed 27,314 differentially methylated CpGs (Δβ ≥ 0.20) between both models. In particular, the number of CpGs that gained (red triangle in scatter plot) and lost (green triangle in scatter plot) a methylation level ≥ 0.20 in SKTR with respect to SK was 14,845 and 12,469, respectively.

Fig. 1
figure 1

Analysis of the DNA methylation profile associated to trastuzumab resistance in breast cancer cell lines. a Schematic flow chart used to identify differentially methylated genes associated to trastuzumab resistance after comparing the sensitive (SK) and resistant (SKTR) to trastuzumab human breast cancer cell lines. b Scatter plot representing DNA methylation normalized levels (β values) of SK and SKTR cell lines. Red and green triangles indicate the CpGs that gained and lost, respectively, a methylation level ≥ 0.20 in SKTR with respect to SK cells. c Supervised hierarchical clustering of the most variable CpGs (Δβ ≥ 0.20) from island and promoter regions between the SK and SKTR cell lines. d Summary of gene ontology (GO) analysis of the biological process categories for the 152 differentially methylated genes at CpG island and promoter levels between the SK and SKTR cell lines. e Venn diagram showing the differentially methylated and differentially expressed genes (obtained by RNA-Seq) between the SK and SKTR cell lines. The name of the 31 hypermethylated and downregulated genes in SKTR is indicated

Therefore, we focused our study on analyzing the methylation levels of the CpG sites located at the regulatory regions corresponding to 89,168 CpGs and 13,508 genes. The supervised hierarchical clustering of the most variable CpGs from promoters and islands (Δβ ≥ 0.20) showed a methylation pattern that clearly discriminated between the SK and SKTR models (Fig. 1c). Next, we used more stringent criteria to determine the most differentially methylated CpGs in the promoters and islands, considering in SKTR the CpGs with a methylation level (β) in SK < 0.20 and in SKTR > 0.60 as hypermethylated, and in SKTR the CpGs with a methylation level (β) in SKTR < 0.20 and in SK > 0.60 as hypomethylated. This analysis revealed 184 differentially methylated CpGs corresponding to 152 genes which, according to a Gene Ontology analysis (GO), were significantly associated (FDR < 0.05) with several biological processes related to cancer, such as cell adhesion pathways (GO:0007155), regulation of transcription (GO:0006355), development (GO:0007275), and control of apoptosis (GO:0043065) (Fig. 1d). To verify whether the methylation observed in these genes had an impact on their transcriptional expression or not, we performed a transcriptomic analysis with RNA-Seq comparing the SK and SKTR models and obtained 1995 overexpressed and 3191 downregulated genes in SKTR displaying a fold change ≥ 1.5 over the SK model. Then, we correlated the genes differentially methylated at the promoter and island with those differentially expressed between SK and SKTR (Fig. 1e). We identified 31 hypermethylated and downregulated genes and no gene hypomethylated and overexpressed in SKTR with respect to the SK model (Table 1). From this list of 31 genes, we selected three genes: TGFBI (transforming growth factor beta induced), CXCL2 (C-X-C motif chemokine ligand 2), and SLC38A1 (Solute carrier Family 38 Member 1) for further analysis. TGFBI, CXCL2, and SLC38A1 were selected because of the high number of differentially methylated CpGs (≥ 3 CpG sites) they presented and their previous implications in BC [25,26,27].

Table 1 Thirty-one differentially methylated and differentially expressed genes between the SK and SKTR models

Epigenetic silencing of TGFBI, CXCL2, and SLC38A1 in trastuzumab-resistant cells (SKTR)

First, we analyzed the DNA methylation of the selected genes by bisulfite pyrosequencing, which is a technique that allows single genes to be analyzed on a single CpG level [28]. We observed significantly higher methylation levels for the three genes in SKTR relative to the SK (Fig. 2a; p < 0.05 for TGFBI and SLC38A1; p < 0.01 for CXCL2). Like the previous transcriptomic analysis with RNA-Seq, the SKTR model also showed a significant reduction in the transcript expression levels of TGFBI, CXCL2, and SLC38A1 (p < 0.05) when compared to the SK cells analyzed by qRT-PCR (Fig. 2b). Interestingly, TGFBI was the gene that showed the greatest decrease (fold change = 6.7) in transcriptional levels in the SKTR model. Similar results were obtained by methylation-specific polymerase chain reaction (MSP) analysis (a very specific and sensitive method [29]) that showed higher methylation in SKTR with respect to the SK cells for all the genes evaluated (Fig. 2d). These results confirmed that hypermethylation of the three selected genes at the promoter and CpG island in SKTR cells is associated with reduce transcriptional levels.

Fig. 2
figure 2

DNA methylation-associated silencing of selected genes comparing the trastuzumab-resistant and trastuzumab-sensitive cell line. a DNA methylation levels of TGFBI, CXCL2, and SLC38A1 in SK and SKTR cell lines by bisulfite pyrosequencing and b methylation-specific polymerase chain reaction (MSP) analysis. c Expression levels of TGFBI, CXCL2, and SLC38A1 in the unmethylated (SK) and methylated (SKTR) cell lines determined by qRT-PCR. d Restored expression of selected genes (TGFBI, CXCL2, and SLC38A1) after DNA demethylating agent 5-aza-2′-deoxycytidine (5-aza-dC) in the SKTR methylated cell line by qRT-PCR. e Protein expression of all selected genes (TGFBI, CXCL2, and SLC38A1) in SK and SKTR cells before and after 5-aza-dC treatment by Western blot. f Analysis in AU and AUTR cell lines of TGFBI methylation by MSP (top) and (middle) its transcriptional (qRT-PCR) and (bottom) protein levels (Western blot) before and after 5-aza-dC treatment. In MSP, the presence of visible polymerase chain reaction products in lanes marked U indicates unmethylated sequences; the presence of products in lanes marked M indicates methylated sequences. In vitro methylated DNA (IVD) was used as a positive control for methylated sequences. DNA from normal lymphocytes (NL) was used as a negative control for methylated sequences. Results shown are representative of those obtained from three independent experiments, and b-actin was used as a control. Values from pyrosequencing and qRT-PCR were determined from triplicates and are expressed as the mean ± SEM. Significance of Mann-Whitney U test, **p < 0.01; *p < 0.05

In this sense, after treating the SKTR model with the demethylating agent 5-aza-2′-deoxycytidine (5-aza-dC) at 3 μM and 5 μM (Fig. 2c), the transcriptional levels of TGFBI, CXCL2, and SLC38A1 were significantly restored (p < 0.05), indicating that methylation is an epigenetic mechanism that has a functional role in the transcriptional control of the three genes in question. Importantly, this epigenetic silencing by DNA methylation observed was also confirmed at the protein level in two (TGFBI and CXCL2) of the three selected genes (Fig. 2e), showing a reduction of the protein levels in SKTR with respect to SK cells that was significantly restored in the hypermethylated SKTR model after the treatment with 5-aza-dC 3 μM and 5 μM. The change observed in the protein levels after the 5-aza-dC treatment was especially drastic in TGFBI gene, where the absence of protein expression induced by the promoter hypermethylation in SKTR was completely recovered after the in vitro demethylation.

Taken together, these results suggest that TGFBI, an extracellular matrix (ECM) protein whose secretion is induced by transforming growth factor-β (TGF-β), is the most epigenetically regulated gene of the three selected.

Therefore, to verify that the results obtained for these genes were not specific to the SK and SKTR models, we extended their epigenetic analyses to other trastuzumab-sensitive (AU) and trastuzumab-resistant (AUTR) human BC cell models (Fig. 2f and Additional file 2: Figure S1). Likewise, only TGFBI gene showed a high correlation between methylation, expression, and protein levels in the AUTR model. With MSP, we observed that hypermethylation of TGFBI in AUTR (compared to AU) was associated with a significant decrease (p < 0.05) in the transcriptional and protein levels of the gene. Importantly, the transcriptional and protein levels of TGFBI gene were both restored in the hypermethylated AUTR model after the treatment with 5-aza-dC 3 μM and 5 μM. These results confirm that the epigenetic silencing of TGFBI by DNA methylation is associated to trastuzumab resistance in human BC cells.

The role of TGFBI expression in trastuzumab resistance cells (SKTR)

We have shown that TGFBI promoter hypermethylation is highly associated with the downregulation of its transcript and protein levels in two trastuzumab-resistant models when compared to the corresponding sensitive models. Therefore, we next examined the functional contribution of epigenetic inactivation of the TGFBI gene to trastuzumab resistance. We depleted the endogenous TGFBI gene expression in the SK model by stable transfection with two different short hairpin RNAs (shTGFBI A and shTGFBI B). We also rescued the TGFBI expression in the TGFBI-hypermethylated and the downregulated SKTR models by stable transfection with a plasmid encoding the full-length TGFBI cDNA (TGFBI).

TGFBI is an extracellular matrix (ECM) protein and plays a role in mediating cell adhesion to the ECM, cell proliferation, adhesion, migration, and differentiation through interacting with collagen, fibronectin, laminin, and several integrins [30,31,32]. These integrin-binding properties of TGFBI have been related to different integrin-binding motifs located in fasciclin-1 domains, including NKDIL, YH18, and EPDIM as well as an Arg-Gly-Asp (RGD) domain [22]. For this reason, we became interested in determining whether this interaction between TGFBI with the ECM and integrins could be involved in the trastuzumab resistance in our models. With this objective in mind, we also expressed a mutated form of TGFBI with four different altered integrin binding motifs (NKDIL, YH18, EPDIM, and RGD), affecting cellular adhesion through the integrin interactions [22]. The transfection and validation protocols were the same as for the overexpression vector.

Despite its known function in cell adhesion and integrin-mediated signaling, we did not observe morphological changes after depletion, overexpression, or mutagenesis of TGFBI (Fig. 3a). The efficiency of the transfection was assessed by measuring the TGFBI gene expression using quantitative RT-PCR and Western blotting (Fig. 3b). A significant change in its levels was detected both after depletion and overexpression. Upon TGFBI transfection, we analyzed the resistance or sensitivity to trastuzumab in comparison to the parental and empty vector-transfected cells through an MTT assay over 5 days (Fig. 3c; Additional file 3: Figure S2). In the SK model, we observed that TGFBI depletion did not affect the trastuzumab response (Fig. 3c, upper panel). In contrast, overexpression of TGFBI in SKTR led to a significantly higher sensitivity at a trastuzumab concentration of 10 μM compared to the resistant empty vector cell line (Mock: p = 0.020) and the SKTR model (p = 0.004) (Fig. 3c, middle panel). The TGFBI mutated form (TGFBImut) showed the same response to trastuzumab treatment as the mock vector did, suggesting that TGFBI-mediated trastuzumab sensitivity requires the NKDIL, YH, EPDIM, and RGD binding motifs (Fig. 3c, lower panel).

Fig. 3
figure 3

Impact of TGFBI expression in trastuzumab-resistant cells. a Representative bright-field microscopy images of TGFBI depletion (Scramble, ShTGFBI A, and ShTGFBI B), overexpression (Mock and TGFBI), and mutagenesis (TGFBImut) in SK and SKTR cells. b Expression analysis by qRT-PCR and Western blot showing the in vitro stable depletion of TGFBI in SK cells (left) and overexpression or mutagenesis of TGFBI in SKTR cells (right). Values of qRT-PCR were determined from triplicates and are expressed as the mean ± SEM. c Cell viability determined by MTT assays upon the use of increasing concentrations of trastuzumab (10−6 to 10 μM) for 5 days. (Upper) The TGFBI depletion in SK did not affect cell viability upon trastuzumab treatment. (Middle) The TGFBI overexpression in SKTR cells (TGFBI) give rise to a major sensitivity to trastuzumab at 10 μM. (Bottom) The TGFBI mutagenesis in SKTR cells (TGFBImut) did not affect cell viability after trastuzumab treatment. Results are expressed as percentage of surviving cells after drug treatment (mean ± SEM). One-way ANOVA using a Tukey HSD post hoc test, **p < 0.01; *p < 0.05 indicate levels of statistical significance. d HER family receptors and their downstream proteins related to PI3K/AKT and MAPK/ERK1/2 pathway characterization in TGFBI-depletion SK cells and TGFBI-overexpression SKTR cells. Western blot showing an enhanced of phosphorylation levels of HER1, HER2, and AKT upon overexpression and mutagenesis of TGFBI in SKTR cells. TGFBI depletion did not produce any change in the HER receptors and their related downstream proteins. Results shown are representative of those obtained from 3 independent experiments and β-actin was used as a control

Therefore, we examined HER family protein receptors and their downstream proteins related to PI3K/AKT and MAPK/ERK1/2 pathways for changes after TGFBI depletion, overexpression, and mutation (Fig. 3d). Consistent with the results of the trastuzumab cell viability, TGFBI depletion in SK cells showed no apparent changes in either total protein or activation levels for any of the proteins examined. Unlike TGFBI depletion, overexpression of TGFBI and TGFBImut in the SKTR model produced changes in the activation levels of some of the HER receptors and downstream signaling proteins in comparison with the mock control vector. In particular, overexpression of either wild type TGFBI or its mutated form resulted in a significant increase in the activation levels of HER1 (pHER1), HER2 (pHER2), and AKT (pAKT) compared to the SKTR model (Mock) that contains hypermethylated TGFBI without change in total levels of the respective proteins. Hence, TGFBI overexpression and its mutated form allow the trastuzumab-resistant model to adopt similar activation levels for HER1, HER2, and AKT as the trastuzumab-sensitive cells do.

These data suggest that selective overexpression of TGFBI in the SKTR model (hypermethylated for TGFBI) induces an increased sensitivity to trastuzumab and the activation of HER1 and HER2 receptors and the AKT downstream protein. Furthermore, the mutated domains from the second and fourth FAS-1 regions of TGFBI are involved in its response to trastuzumab treatment, but not in its activation or interaction with HER2 downstream proteins. In summary, the SKTR model with an overexpression of TGFBI presents a behavior similar to the SK model.

TGFBI hypermethylation in trastuzumab-resistant HER2+ breast cancer patients

Given the previous in vitro results obtained, we wanted to determine whether the presence of the TGFBI promoter CpG island hypermethylation-associated to trastuzumab resistance also occurring in BC patients. Therefore, we evaluated TGFBI gene methylation in human tumors from 24 HER2+ early BC patients as a preliminary study. As with our sensitive and resistant models, we analyzed TGFBI methylation by bisulfite pyrosequencing in pre-treatment samples from patients who were complete responders and non-responders to trastuzumab-based therapy (Fig. 4a). With the non-responder patients, we also evaluated TGFBI methylation after treatment (post-treatment samples), i.e., as we did in our SKTR model. For the non-response patients, we obtained paired samples (pre-treatment and post-treatment samples) from 10 patients, pre-treatment samples only from 3 patients, and post-treatment samples only from 7 patients. For the patients with pathological complete response, only pre-treatment samples were evaluated. The limited number of samples obtained was, in part, due to the exhaustion of the sample during diagnostic procedures and the poor preservation of the DNA in the paraffin tissue (FFPE) blocks.

Fig. 4
figure 4

TGFBI promoter hypermethylation in HER2+ breast cancer patients with sensitivity and resistance to trastuzumab. a Schematic representation of selected patient samples. TGFBI methylation levels were evaluated in tumor samples of 24 HER2+ breast cancer samples. From this cohort, after trastuzumab plus chemotherapy in neoadjuvant regimen, 20 patients developed partial or no response and 4 patients presented complete treatment response. Of the 20 patients with non-response, 10 patients had pre-treatment and post-treatment samples, 3 patients with pre-treatment only samples, and 7 with post-treatment only samples (Top). The 4 patients with complete response to treatment only had pre-treatment samples (Bottom). b TGFBI methylation of 3 consecutive CpG sites in the 5′-end promoter CpG island in HER2+ breast cancer treated with trastuzumab analyzed by bisulfite pyrosequencing. The central solid line indicates the median and the limits of the vertical lines show the upper and lower percentiles. c TGFBI methylation of 3 consecutive CpG sites in the 5′-end promoter CpG island in resistant patients with paired pre- and post-treatment samples analyzed by bisulfite pyrosequencing. d Diagnostic accuracy of TGFBI hypermethylation for resistant HER2+ breast cancer samples. ROC analysis was applied to the TGFBI methylation levels analyzed by pyrosequencing for trastuzumab-resistant samples (pre- and post-treatment samples). Area under the curve (AUC) was 0.9502 (p < 0.0001). TGFBI showed great potential for monitoring trastuzumab response in HER2+ breast cancer patients. Significance of Mann-Whitney U test is indicated as ***p < 0.001; **p < 0.01

Our results showed similar pre-treatment TGFBI methylation levels (p = 0.651) in the patients with complete response to trastuzumab (6.45% ± 1.90) and in the non-responders (6.08% ± 1.51; Fig. 4b). These results indicate that the TGFBI methylation levels of the pre-treatment samples are not associated with the absence of response to trastuzumab. In contrast, the non-responsive patients showed significantly higher (p = 0.001) methylation levels of TGFBI in tumors following treatment with trastuzumab (30.26% ± 3.52) than before starting the therapy (6.08% ± 1.51), suggesting that acquired resistance to trastuzumab is associated with increased methylation levels of TGFBI. In particular, when considering the non-responsive patients with pre- and post-treatment paired samples, we observed a significant (p < 0.001) TGFBI promoter hypermethylation after trastuzumab (8 out of 10 = 80%) compared to the pre-treatment samples, indicating that the increase in TGFBI methylation levels is associated with trastuzumab resistance (Fig. 4c). Importantly, the ROC curve analysis with an AUC of 0.95 (p < 0.0001; 95% CI 0.803 to 0.996) allowed us to clearly differentiate the methylation levels of TGFBI between the pre- and post-treatment samples of the non-responsive patients. This result suggests that TGFBI could be a potential biomarker for monitoring trastuzumab response during HER2+ BC patients’ treatment (Fig. 4d). In addition, the Hosmer-Lemeshow test (chi-square = 89,606; p = 0.3456) suggests that the methylation of TGFBI has a certain ability to predict the resistance to trastuzumab when compared with an ideal model (Additional file 4: Figure S3). On the other hand, no significant association between TGFBI methylation levels before and after treatment and the clinical-histopathological characteristics was identified (Table 2).

Table 2 Clinical and pathological characteristics in HER2-positive early breast cancer according to TGFBI promoter methylation levels


HER2 is overexpressed in 20–30% of BC and is associated with aggressive phenotype and poor prognosis. Despite the initial good response to the treatment, a large portion of patients present de novo or acquired treatment resistance. Currently, HER2 detection is the only validated biomarker for predicting the benefit of anti-HER2 therapies [33]. While different cellular and molecular mechanisms involved in trastuzumab resistance have been described, none of them are used to detect, predict, or monitor BC treatment in a clinical routine [7, 8, 34]. Thus, there is great interest in finding potential biomarkers able to detect, predict, or monitor trastuzumab response which, in turn, would allow us to stratify patients and save costs and toxicities for those who do not respond to treatment.

Epigenetic biomarker studies are focused on analyzing the methylation changes in the promoter regions of candidate genes in specific tumor types for the implication in gene silencing [9, 34]. In general, this gene repression results in an adaptive advantage for the cells, allowing a more aggressive and invasive phenotype to be adopted. After analyzing the promoter methylation profile of long-term trastuzumab-sensitive and trastuzumab-resistant models developed in our laboratory [16, 17], three different methylated genes (TGFBI, CXCL2, and SLC38A1) implicated in cancer were identified and subsequently validated using different methylation and expression approaches. Additionally, the selected genes were validated in two other long-term trastuzumab-sensitive (AU) and trastuzumab-resistant (AUTR) human HER2+ BC cell lines developed by our group [17]. Although both cell lines analyzed derive from the same patient and have similar genetic and transcriptomic profiles, they have some differences in transcriptional levels [35]. Moreover, SKBr3 and AU565 have been described as having distinct responses to drug treatment depending on culture conditions [36]. From all the genes selected, TGFBI (also known as Big-H3 or keratoepithelin) was the only one with a consistent methylation, expression, and protein pattern in both resistant models. This epigenetic silencing of the TGFBI gene observed was in line with previous studies of other cancers [37, 38].

As explained before, TGFBI encodes a 68-kDa secretory protein induced by transforming growth factor-β (TGF-β) in human adenocarcinoma cells as well as other human cell types [39]. According to the bibliography, TGFBI protein is composed by 683 amino acids containing a secretory signal (SP) in the N-terminal cysteine-rich domain (CRD), and four fasciclin-1 domains (FAS1-1, FAS1-2, FAS1-3, and FAS1-4), which contain several known integrin-binding motifs including NKDIL, YH18, and EPDIM as well as an Arg-Gly-Asp (RGD) domain [22]. TGFBI has been reported to function as an ECM protein to mediate cell adhesion and migration through interacting with collagen, fibronectin, and laminin and several integrins including α1β1, α3β1, αvβ3, and αvβ5 [30,31,32]. While TGFBI has been reported to be involved in tumorigenesis, its role is not clear. In this paper, we have linked the promoter DNA methylation-associated silencing of TGFBI with a possible tumor suppressor function in trastuzumab-resistant models. Several previous reports have indicated that TGFBI plays the role of a tumor suppressor gene in various cancers such as lung, breast, and ovarian [25, 40,41,42]. However, in other cancers, such as colon or pancreas, TGFBI has been described as having a tumor-promoting function [43,44,45,46]. These opposing effects of TGFBI suggest that its expression and function are dependent on cell type [22].

TGFBI promotes cell adhesion by interacting with different integrins, which, in turn, has been shown to play an important role in tumor progression in humans. TGFBI’s inhibition of invasive capacity correlates with a previous analysis by our group, which revealed a greater invasive capacity for SKTR compared to the SK [16]. To determine if TGFBI has a potential role in trastuzumab resistance, different functional studies were carried out. Our results show that TGFBI overexpression in SKTR significantly sensitizes the cells to the treatment and activates some HER2 downstream protein pathways (pAKT, pHER1, and pHER2), adopting a similar pattern to SK. In contrast, TGFBI knockdown in TGFBI endogenous expressed in the SK cells did not produce any effect. A direct interaction between the HER2 receptor and different integrins has been described in multiple reports [47, 48]. Moreover, different integrins have been related to trastuzumab resistance such as β1 or α6β4 [49, 50]. Therefore, TGFBI mutated vector was performed to investigate the role of the RGD motif and the three integrin-binding motifs NKDIL, EPDIM, and YH18 in trastuzumab response [22]. Like TGFBI overexpression, the TGFBI mutated form induced the activation of different HER2 downstream proteins, but there were no differences in trastuzumab treatment response. Therefore, it is likely that the four mutated domains of TGFBI are required for its function in trastuzumab response, which is consistent with other studies that described the involvement of FAS1 or RGD motifs in tumor angiogenesis and tumor growth inhibition, as well as in promoting apoptosis [31, 51,52,53]. Moreover, it has also been described that TGFBI is involved in mesothelioma progression through the AKT/mTOR pathway which could be related to the HER2 pathway changes observed in TGFBI overexpression and mutagenesis [54]. In summary, this study suggests that TGFBI expression may promote the effectiveness of trastuzumab treatment.

The methylation of TGFBI showed certain capacity to predict the resistance to trastuzumab according to the Hosmer-Lemeshow test; however, these results should be interpreted with caution because the number of grouped patients was small. Although more in vitro functional analysis are necessary to elucidate the TGFBI role in trastuzumab resistance, we suggest some hypothetical situations taking into account all TGFBI action mechanisms. It has been described that TGFBI could inhibit cell adhesion to various ECM proteins inhibiting cell proliferation and invasion in neuroblastoma [55]. Moreover, cell adhesion to ECM has been consistently reported as one of the mechanisms used by tumor cells to resist chemotherapy [56]. These observations are in correlation with previous studies by our group, where it is exposed that our SKTR model presented a high significant adhesion capacity to bind to extracellular matrix proteins like fibronectin, collagen I, collagen IV, or laminin I compared to SK [16]. Therefore, we hypothesized that in the absence of TGFBI, the integrins could bind to the ECM providing attachment sites for cells inducing more invasion, or acting as a physical drug barrier, restricting drug transport and limiting their efficacy. The intratumoral diffusion and physical masking of HER receptors by ECM proteins have also been described as a mechanism that can affect the therapeutic efficacy of some drugs including trastuzumab [57].

Finally, we correlated the previous in vitro results in human samples treated with neoadjuvant anthracycline-taxane-based chemotherapy plus trastuzumab. In our preliminary study, a higher TGFBI methylation level after treatment was observed in patients who developed treatment resistance. These results are in accordance with other studies which demonstrated an association between TGFBI hypermethylation and poor prognosis in prostate and lung cancer [41]. Moreover, decreased TGFBI expression was identified in advanced stages of BC and NSCLC tumors [40, 58]. Due to the small size of the patient cohort and the retrospective design, the role of TGFBI methylation as a biomarker requires further validation in a larger and independent cohort. Nevertheless, the present study brings to light for the first time that TGFBI methylation has a significant discriminative value between pre- and post-treatment samples, as demonstrated by ROC curve analyses. Interestingly, the TGFBI promoter methylation analysis in responder patients showed similar methylation levels with pre-treatment samples of non-responders. This observation suggests a possible role of TGFBI as a monitoring biomarker for trastuzumab response in patients with HER2+ BC.

Although no significant association was found between TGFBI methylation before and after treatment and patients’ clinical-histopathologic characteristics, a high number ER-positive patients were observed. In this sense, previous studies have shown that there is crosstalk between ER and HER2 pathways and it affects the response to treatment [59]. In addition, although the type of surgery is not associated with a higher TGFBI methylation level in post-treatment samples, 73% of our patients were treated with mastectomy versus 27% with lumpectomy, probably due to the high prevalence of stage III. The low histological grade and high clinical stage to treatment were probably due to most of the patients having non-operable cancers including inflammatory tumors. Currently, neoadjuvant treatment is generally used for operable HER2+ BC thanks to an improvement in the efficacy of the drugs used for treatment [60]. At the time our cohort was treated, none of these drugs were yet available. That said, trastuzumab remains the gold standard treatment for HER2+ BC.


The results obtained in this study have provided an overview of the DNA methylation pattern in HER2+-resistant BC. TGFBI promoter hypermethylation could be a potential methylation monitoring biomarker for trastuzumab response. In addition, once the resistance is developed, demethylation of TGFBI with compounds that inhibit DNA methyltransferases could contribute to the sensitization of breast cancer cells to the trastuzumab treatment. However, further studies are required to identify the specific role TGFBI plays in trastuzumab resistance in the neoadjuvant settings to definitively show TGFBI as a clinical relevant biomarker. Future studies based on TGFBI promoter methylation analyses in circulating DNA during neoadjuvant treatment could be a potential strategy with which to confirm our findings.

Availability of data and materials

The microarray and sequencing dataset generated during the current study is available in NCBI GEO ( under accession GSE123754 and GSE114575, respectively.





Protein kinase B


Breast cancer


C-X-C motif chemokine ligand 2


Extracellular matrix


False discovery rate


Gene ontology


Human epidermal growth factor receptor family


Human epidermal growth factor receptor 2




Ras-mitogen-activated protein kinases


Methylation-specific PCR


Phosphatidylinositol 3′-kinase


Solute carrier Family 38 Member 1


Transforming growth factor beta


Transforming growth factor beta induced


  1. Global Burden of Disease Cancer Collaboration, Fitzmaurice C, Allen C, Barber RM, Barregard L, Bhutta ZA, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 32 cancer groups, 1990 to 2015: a systematic analysis for the global burden of disease study. JAMA Oncol. 2017;3(4):524–48.

    PubMed  PubMed Central  Article  Google Scholar 

  2. Bauer K, Parise C, Caggiano V. Use of ER/PR/HER2 subtypes in conjunction with the 2007 St Gallen consensus statement for early breast cancer. BMC Cancer. 2010;10(1):228.

    PubMed  PubMed Central  Article  Google Scholar 

  3. Hudis CA. Trastuzumab--mechanism of action and use in clinical practice. N Engl J Med. 2007;357(1):39–51.

    CAS  PubMed  Article  Google Scholar 

  4. Luque-Cabal M, García-Teijido P, Fernández-Pérez Y, Sánchez-Lorenzo L, Palacio-Vázquez I. Mechanisms behind the resistance to trastuzumab in HER2-amplified breast cancer and strategies to overcome it. Clin Med Insights Oncol. 2016;10(Suppl 1):21–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  5. Nakashoji A, Hayashida T, Yokoe T, Maeda H, Toyota T, Kikuchi M, et al. The updated network meta-analysis of neoadjuvant therapy for HER2-positive breast cancer. Cancer Treat Rev. 2018;62:9–17.

    CAS  PubMed  Article  Google Scholar 

  6. Carter P, Presta L, Gorman CM, Ridgway JB, Henner D, Wong WL, et al. Humanization of an anti-p185HER2 antibody for human cancer therapy. Proc Natl Acad Sci U S A. 1992;89(10):4285–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. Nahta R, Yu D, Hung M-C, Hortobagyi GN, Esteva FJ. Mechanisms of disease: understanding resistance to HER2-targeted therapy in human breast cancer. Nat Clin Pract Oncol. 2006;3(5):269–80.

    CAS  PubMed  Article  Google Scholar 

  8. Scaltriti M, Rojo F, Ocaña A, Anido J, Guzman M, Cortes J, et al. Expression of p95HER2, a truncated form of the HER2 receptor, and response to anti-HER2 therapies in breast cancer. J Natl Cancer Inst. 2007;99(8):628–38.

    CAS  PubMed  Article  Google Scholar 

  9. Esteller M. Epigenetics in Cancer. N Engl J Med. 2008;358(11):1148–59.

    CAS  PubMed  Article  Google Scholar 

  10. Paluszczak J, Baer-Dubowska W. Epigenetic diagnostics of cancer — the application of DNA methylation markers. J Appl Genet. 2006;47(4):365–75.

    PubMed  Article  Google Scholar 

  11. Duffy MJ, Napieralski R, Martens JWM, Span PN, Spyratos F, Sweep FCGJ, et al. Methylated genes as new cancer biomarkers. Eur J Cancer. 2009;45(3):335–46.

    CAS  PubMed  Article  Google Scholar 

  12. Rice JC, Ozcelik H, Maxeiner P, Andrulis I, Futscher BW. Methylation of the BRCA1 promoter is associated with decreased BRCA1 mRNA levels in clinical breast cancer specimens. Carcinogenesis. 2000;21(9):1761–5.

    CAS  PubMed  Article  Google Scholar 

  13. Dammann R, Schagdarsurengin U, Strunnikova M, Rastetter M, Seidel C, Liu L, et al. Epigenetic inactivation of the Ras-association domain family 1 (RASSF1A) gene and its function in human carcinogenesis. Histol Histopathol. 2003;18(2):665–77.

    CAS  PubMed  Google Scholar 

  14. Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.

    CAS  PubMed  Article  Google Scholar 

  15. Trempe GL. Human breast cancer in culture. Recent Results Cancer Res. 1976;57(57):33–41.

    Google Scholar 

  16. Blancafort A, Giró-Perafita A, Oliveras G, Palomeras S, Turrado C, Campuzano Ò, et al. Dual fatty acid synthase and HER2 signaling blockade shows marked antitumor activity against breast cancer models resistant to anti-HER2 drugs. PLoS One. 2015;10(6):e0131241.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  17. Puig T, Aguilar H, Cufí S, Oliveras G, Turrado C, Ortega-gutiérrez S, et al. A novel inhibitor of fatty acid synthase shows activity against HER2 + breast cancer xenografts and is active in anti-HER2 drug-resistant cell lines. Breast Cancer Res. 2011;13(6):R131.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. Edge SB, Compton CC. The American Joint Committee on Cancer: the 7th edition of the AJCC Cancer staging manual and the future of TNM. Ann Surg Oncol. 2010;17(6):1471–4.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Bender CM, Pao MM, Jones PA. Inhibition of DNA methylation by 5-aza-2′-deoxycytidine suppresses the growth of human tumor cell lines. Cancer Res. 1998;58(1):95–101.

    CAS  PubMed  Google Scholar 

  20. 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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. Naumov VA, Generozov EV, Zaharjevskaya NB, Matushkina DS, Larin AK, Chernyshov SV, et al. Genome-scale analysis of DNA methylation in colorectal cancer using Infinium HumanMethylation450 BeadChips. Epigenetics. 2013;8(9):921–34.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. Ween MP, Oehler MK, Ricciardelli C. Transforming growth factor-beta-induced protein (TGFBI)/(βig-H3): a matrix protein with dual functions in ovarian cancer. Int J Mol Sci. 2012;13(8):10461–77.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 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. 2016;113(47):E7535–44.

    CAS  PubMed  Article  Google Scholar 

  24. Diaz-Lagares A, Mendez-Gonzalez J, Hervas D, Saigi M, Pajares MJ, Garcia D, et al. A novel epigenetic signature for early diagnosis in lung cancer. Clin Cancer Res. 2016;22(13):3361–71.

    CAS  PubMed  Article  Google Scholar 

  25. Wen G, Partridge MA, Li B, Hong M, Liao W, Cheng SK, et al. TGFBI expression reduces in vitro and in vivo metastatic potential of lung and breast tumor cells. Cancer Lett. 2011;308(1):23–32.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. Kronski E, Fiori ME, Barbieri O, Astigiano S, Mirisola V, Killian PH, et al. miR181b is induced by the chemopreventive polyphenol curcumin and inhibits breast cancer metastasis via down-regulation of the inflammatory cytokines CXCL1 and −2. Mol Oncol. 2014;8(3):581–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. Wang K, Cao F, Fang W, Hu Y, Chen Y, Ding H, et al. Activation of SNAT1/SLC38A1 in human breast cancer: correlation with p-Akt overexpression. BMC Cancer. 2013;13:343.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. Dejeux E, El Abdalaoui H, Gut IG, Tost J. Identification and quantification of differentially methylated loci by the PyrosequencingTM technology. Methods Mol Biol. 2009;507:189–205.

    CAS  PubMed  Article  Google Scholar 

  29. Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB. Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands (DNA methylation/tumor suppressor genes/pl6/p15). Proc Natl Acad Sci U S A. 1996;93:9821–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. Hashimoto K, Noshiro M, Ohno S, Kawamoto T, Satakeda H, Akagawa Y, et al. Characterization of a cartilage-derived 66-kDa protein (RGD-CAP/βig-h3) that binds to collagen. Biochim Biophys Acta - Mol Cell Res. 1997;1355(3):303–14.

    CAS  Article  Google Scholar 

  31. Kim J-E, Jeong H-W, Nam J-O, Lee B-H, Choi J-Y, Park R-W, et al. Identification of motifs in the fasciclin domains of the transforming growth factor-β-induced matrix protein βig-h3 that interact with the αvβ5 integrin. J Biol Chem. 2002;277(48):46159–65.

    CAS  PubMed  Article  Google Scholar 

  32. Billings PC, Charles Whitbeck J, Adams CS, Abrams WR, Cohen AJ, Engelsberg BN, et al. The transforming growth factor-β-inducible matrix protein βig-h3 interacts with fibronectin. J Biol Chem. 2002;277(31):28003–9.

    CAS  PubMed  Article  Google Scholar 

  33. Slamon DJ, Leyland-Jones B, Shak S, Fuchs H, Paton V, Bajamonde A, et al. Use of chemotherapy plus a monoclonal antibody against HER2 for metastatic breast cancer that overexpresses HER2. N Engl J Med. 2001;344(11):783–92.

    CAS  PubMed  Article  Google Scholar 

  34. Nahta R, Yuan LXHH, Zhang B, Kobayashi R, Esteva FJ. Insulin-like growth factor-I receptor/human epidermal growth factor receptor 2 heterodimerization contributes to trastuzumab resistance of breast cancer cells. Cancer Res. 2005;65(23):11118–28.

    CAS  PubMed  Article  Google Scholar 

  35. Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, et al. A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell. 2006;10(6):515–27.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. Weigelt B, Lo AT, Park CC, Gray JW, Bissell MJ. HER2 signaling pathway activation and response of breast cancer cells to HER2-targeting agents is dependent strongly on the 3D microenvironment. Breast Cancer Res Treat. 2010;122(1):35–43.

    CAS  PubMed  Article  Google Scholar 

  37. Shao G, Berenguer J, Borczuk AC, Powell CA, Hei TK, Zhao Y. Epigenetic inactivation of Betaig-h3 gene in human cancer cells. Cancer Res. 2006;66(9):4566–73.

    CAS  PubMed  Article  Google Scholar 

  38. Wang N, Zhang H, Yao Q, Wang Y, Dai S, Yang X. TGFBI promoter hypermethylation correlating with paclitaxel chemoresistance in ovarian cancer. J Exp Clin Cancer Res. 2012;31(1):6.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  39. Skonier J, Neubauer M, Madisen L, Bennett K, Plowman GD, Purchio AF. cDNA cloning and sequence analysis of βig-h3, a novel gene induced in a human adenocarcinoma cell line after treatment with transforming growth factor-β. DNA Cell Biol. 1992;11(7):511–22.

    CAS  PubMed  Article  Google Scholar 

  40. Calaf GM, Echiburú-Chau C, Zhao YL, Hei TK. BigH3 protein expression as a marker for breast cancer. Int J Mol Med. 2008;21(5):561–8.

    PubMed  Google Scholar 

  41. Shah JN, Shao G, Hei TK, Zhao Y. Methylation screening of the TGFBI promoter in human lung and prostate cancer by methylation-specific PCR. BMC Cancer. 2008;8(1):284.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  42. Ahmed AA, Mills AD, Ibrahim AEKK, Temple J, Blenkiron C, Vias M, et al. The extracellular matrix protein TGFBI induces microtubule stabilization and sensitizes ovarian cancers to paclitaxel. Cancer Cell. 2007;12(6):514–27.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. Kitahara O, Furukawa Y, Tanaka T, Kihara C, Ono K, Yanagawa R, et al. Alterations of gene expression during colorectal carcinogenesis revealed by cDNA microarrays after laser-capture microdissection of tumor tissues and normal epithelia. Cancer Res. 2001;61(9):3544–9.

    CAS  PubMed  Google Scholar 

  44. Schneider D, Kleeff J, Berberat PO, Zhu Z, Korc M, Friess H, et al. Induction and expression of betaig-h3 in pancreatic cancer cells. Biochim Biophys Acta. 2002;1588(1):1–6.

    CAS  PubMed  Article  Google Scholar 

  45. Guo S-K, Shen M-F, Yao H-W, Liu Y-S. Enhanced expression of TGFBI promotes the proliferation and migration of glioma cells. Cell Physiol Biochem. 2018;49(3):1097–109.

    PubMed  Article  CAS  Google Scholar 

  46. Chen W-Y, Tsai Y-C, Yeh H-L, Suau F, Jiang K-C, Shao A-N, et al. Loss of SPDEF and gain of TGFBI activity after androgen deprivation therapy promote EMT and bone metastasis of prostate cancer. Sci Signal. 2017;10(492):eaam6826.

    PubMed  Article  CAS  Google Scholar 

  47. Desgrosellier JS, Cheresh DA. Integrins in cancer: biological implications and therapeutic opportunities. Nat Rev Cancer. 2010;10(1):9–22.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  48. Soung YH, Clifford JL, Chung J. Crosstalk between integrin and receptor tyrosine kinase signaling in breast carcinoma progression. BMB Rep. 2010;43(5):311–8.

    CAS  PubMed  Article  Google Scholar 

  49. Lesniak D, Xu Y, Deschenes J, Lai R, Thoms J, Murray D, et al. 1-integrin circumvents the antiproliferative effects of trastuzumab in human epidermal growth factor receptor-2-positive breast cancer. Cancer Res. 2009;69(22):8620–8.

    CAS  PubMed  Article  Google Scholar 

  50. Guo W, Pylayeva Y, Pepe A, Yoshioka T, Muller WJ, Inghirami G, et al. β4 integrin amplifies ErbB2 signaling to promote mammary tumorigenesis. Cell. 2006;126(3):489–502.

    CAS  PubMed  Article  Google Scholar 

  51. Nam J-O, Jeong H-W, Lee B-H, Park R-W, Kim I-S. Regulation of tumor angiogenesis by fastatin, the fourth FAS1 domain of βig-h3, via αvβ3 integrin. Cancer Res. 2005;65(10):4153–61.

    CAS  PubMed  Article  Google Scholar 

  52. Zhao YL, Piao CQ, Hei TK. Overexpression of Betaig-h3 gene downregulates integrin alpha5beta1 and suppresses tumorigenicity in radiation-induced tumorigenic human bronchial epithelial cells. Br J Cancer. 2002;86(12):1923–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. Kim J-E, Kim S-J, Jeong H-W, Lee B-H, Choi J-Y, Park R-W, et al. RGD peptides released from βig-h3, a TGF-β-induced cell-adhesive molecule, mediate apoptosis. Oncogene. 2003;22(13):2045–53.

    CAS  PubMed  Article  Google Scholar 

  54. Wen G, Hong M, Li B, Liao W, Cheng SK, Hu B, et al. Transforming growth factor-β-induced protein (TGFBI) suppresses mesothelioma progression through the Akt/mTOR pathway. Int J Oncol. 2011;39(4):1001–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Becker J, Volland S, Noskova I, Schramm A, Schweigerer LL, Wilting J. Keratoepithelin reverts the suppression of tissue factor pathway inhibitor 2 by MYCN in human neuroblastoma: a mechanism to inhibit invasion. Int J Oncol. 2008;32(1):235–40.

    CAS  PubMed  Google Scholar 

  56. Bissey P-A, Law JH, Bruce JP, Shi W, Renoult A, Chua MLK, et al. Dysregulation of the MiR-449b target TGFBI alters the TGFβ pathway to induce cisplatin resistance in nasopharyngeal carcinoma. Oncogenesis. 2018;7(5):40.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  57. Beyer I, Li Z, Persson J, Liu Y, van Rensburg R, Yumul R, et al. Controlled extracellular matrix degradation in breast cancer tumors improves therapy by trastuzumab. Mol Ther. 2011;19(3):479–89.

    CAS  PubMed  Article  Google Scholar 

  58. Pajares MJ, Agorreta J, Salvo E, Behrens C, Wistuba II, Montuenga LM, et al. TGFBI expression is an independent predictor of survival in adjuvant-treated lung squamous cell carcinoma patients. Br J Cancer. 2014;110(6):1545–51.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  59. Schettini F, Buono G, Cardalesi C, Desideri I, De Placido S, Del Mastro L. Hormone receptor/human epidermal growth factor receptor 2-positive breast cancer: where we are now and where we are going. Cancer Treat Rev. 2016;46:20–6.

    CAS  PubMed  Article  Google Scholar 

  60. Gianni L, Pienkowski T, Im Y-H, Tseng L-M, Liu M-C, Lluch A, et al. 5-year analysis of neoadjuvant pertuzumab and trastuzumab in patients with locally advanced, inflammatory, or early-stage HER2-positive breast cancer (NeoSphere): a multicentre, open-label, phase 2 randomised trial. Lancet Oncol. 2016;17(6):791–800.

    CAS  PubMed  Article  Google Scholar 

Download references


The authors thank Aitor Rodriguez-Casanova (Translational Medical Oncology (Oncomet), Health Research Institute of Santiago (IDIS), University Clinical Hospital of Santiago (CHUS), Santiago de Compostela, Spain) for his kind support in the methylation analysis. They also thank Maria Buxó (Girona Biomedical Research Institute (IDIBGI), Girona, Spain) and Marc Saez (CIBER of Epidemiology and Public Health (CIBERESP), Madrid, Spain and Research Group on Statistics, Econometrics and Health (GRECS), University of Girona, Girona, Spain) for assistance with the statistical analysis. We thank the Genomics and Epigenomics Service from the Cancer Epigenetics and Biology Program (PRBC), Bellvitge Biomedical Research Institute (IDIBELL)], the Huntsman Cancer Foundation, and the High Throughput Genomics and Bioinformatic Analysis Shared Resource (HTGBA) at Huntsman Cancer Institute (HCI) for their support. HTGBA is supported by HCI Cancer Center Support Grant (5P30CA042014-24; the content is solely the responsibility of the authors and does not necessarily represent the official views of the U.S. National Cancer Institute or the National Institutes of Health).


This work was supported in part by the Spanish Instituto de Salud Carlos III (ISCIII; FIS PI11/00692 and PI14/00329; to’ T. Puig), Fundacion Ramon Areces (to T. Puig), the support of Catalonian Government (2017SGR00385), and La Marató de TV3 (20131530, TPuig), financial support was from the University of Girona (MPCUdG2016/036), and the University of Girona and La Caixa Foundation awarded S. Palomeras with a predoctoral grant. A.D.-L. is funded by a “Juan Rodés” contract (JR17/00016) from ISCIII. Ana B Crujeiras is funded by a research contract “Miguel Servet” (CP17/00088) from the ISCIII.

Author information




TP and ME conceived and designed the study. SP, AD-L, FS, HJF, DHL, ALW, GO, ME, and TP developed the methodology. SP, AD-L, DHL, ALW, GV, AH, TP, and ME contributed to data acquisition (provided facilities etc.). SP, AD-L, ABC, GV, and AH participated in data analysis and interpretation. SP, AD-L, GV, FS, HJF, ABC, DL, ALW, TP, and ME wrote, reviewed, and/or revised the manuscript. SP, GV, AH, DL, and ALW participated in administrative, technical, or material support (i.e., reporting or organizing data, construction database). TP and ME supervised the study. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Manel Esteller or Teresa Puig.

Ethics declarations

Ethics approval and consent to participate

This study was performed in compliance with the clinical research guidelines established by the Dr. Josep Trueta University Hospital Ethics Committee. Signed informed consent was obtained from each patient.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Additional methods. (DOCX 30 kb)

Additional file 2:

Figure S1. DNA methylation-associated silencing of CXCL2 and SLC38A1 comparing the trastuzumab-resistant (AUTR) and -sensitive (AU) AU565 cell model. (A) Expression levels of CXCL2 and SLC38A1 in the unmethylated (AU) and methylated (AUTR) models determined by qRT-PCR. (B) Restored expression of CXCL2 and SLC38A1 after DNA demethylating agent 5-aza-2′-deoxycytidine (5-aza-dC) in the AUTR methylated cell line by qRT-PCR. (C) DNA methylation levels of CXCL2 and SLC38A1 in AU and AUTR cell lines by methylation-specific polymerase chain reaction (MSP) analysis. (D) Protein expression of CXCL2 and SLC38A1 in AU and AUTR cells before and after 5-aza-dC treatment by Western blot. In MSP, the presence of visible polymerase chain reaction products in lanes marked U indicate unmethylated sequences; the presence of products in lanes marked M indicate methylated sequences. In vitro methylated DNA (IVD) was used as a positive control for methylated sequences. DNA from normal lymphocytes (NL) was used as a negative control for methylated sequences. Results shown are representative of those obtained from three independent experiments and b- actin was used as a control. Values from pyrosequencing and qRT-PCR were determined from triplicates and are expressed as the mean ± SEM. Significance of Mann-Whitney U test, **p < 0.01; *p < 0.05. (PDF 235 kb)

Additional file 3:

Figure S2. Cell viability determined by MTT assays upon trastuzumab treatment for 5 days. (A) TGFBI depletion in SK cells, (B) TGFBI overexpression in SKTR cells (TGFBI) and (C) TGFBI-mutagenesis in SKTR cells (TGFBImut). The continuous line represents the fitted dose-response curve while discontinuous lines indicate the simultaneous 95% confidence interval (CIs) for the continuous line with the same color. ANOVA using a Tukey HSD post hoc test, **P < 0.01; *P < 0.05 indicate levels of statistical significance. (PDF 397 kb)

Additional file 4:

Figure S3. Calibration curve for the reliability of TGFBI to predict trastuzumab resistance. Calibration was assessed using the Hosmer-Lemeshow test and graphically using a calibration plot. The observed frequency was calculated for each group. Predicted probability refers to the predicted probabilities generated by the model. Patients were grouped according to deciles of the predicted probabilities. Grouped patients refers to patients grouped at each decile of predicted probabilities generated by the model. A LOWESS regression was fit to all data points (blue line). (PDF 119 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Palomeras, S., Diaz-Lagares, Á., Viñas, G. et al. Epigenetic silencing of TGFBI confers resistance to trastuzumab in human breast cancer. Breast Cancer Res 21, 79 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • HER2+ breast cancer
  • Trastuzumab resistance
  • DNA methylation