Frequent loss of endothelin-3 (EDN3) expression due to epigenetic inactivation in human breast cancer
Breast Cancer Research volume 11, Article number: R34 (2009)
Endothelin (EDN) signalling plays a crucial role in cell differentiation, proliferation and migration processes. There is compelling evidence that altered EDN signalling is involved in carcinogenesis by modulating cell survival and promoting invasiveness. To date, most reports have focused on the oncogenic potential of EDN1 and EDN2, both of which are overexpressed in various tumour entities. Here, we aimed at a first comprehensive analysis on EDN3 expression and its implication in human breast cancer.
EDN3 mRNA expression was assessed by Northern blotting in normal human tissues (n = 9) as well as in matched pairs of normal and tumourous tissues from breast specimens (n = 50). EDN3 mRNA expression in breast cancer was further validated by real-time polymerase chain reaction (PCR) (n = 77). A tissue microarray was used to study EDN3 protein expression in breast carcinoma (n = 150) and normal breast epithelium (n = 44). EDN3 promoter methylation was analysed by methylation-specific PCR in breast cell lines (n = 6) before and after demethylating treatment, normal breast tissues (n = 17) and primary breast carcinomas (n = 128). EDN3 expression and methylation data were statistically correlated with clinical patient characteristics and patient outcome.
Loss of EDN3 mRNA expression in breast cancer, as initially detected by array-based expression profiling, could be confirmed by Northern blot analysis (> 2-fold loss in 96%) and real-time PCR (> 2-fold loss in 78%). Attenuated EDN3 expression in breast carcinoma was also evident at the protein level (45%) in association with adverse patient outcome in univariate (P = 0.022) and multivariate (hazard ratio 2.0; P = 0.025) analyses. Hypermethylation of the EDN3 promoter could be identified as the predominant mechanism leading to gene silencing. Reversion of the epigenetic lock by 5-aza-2'-deoxycytidine and trichostatin A resulted in EDN3 mRNA re-expression in vitro. Furthermore, EDN3 promoter hypermethylation was detected in 70% of primary breast carcinomas with significant association to loss of EDN3 mRNA expression (P = 0.005), whilst normal matched breast tissues revealed no EDN3 promoter methylation.
EDN3 is a frequent target of epigenetic inactivation in human breast cancer, potentially contributing to imbalanced EDN signalling commonly found in this disease. The clinical implication supports the view that EDN3, in contrast to EDN1 and EDN2, may act as natural tumour suppressor in the human mammary gland.
Endothelins (EDNs) are widely expressed cytokines in a variety of human tissues, including brain, skeletal muscle, pancreas, small intestine, testis and colon . They constitute a family of small, vasoactive, 21-amino acid peptides referred to as EDN1, EDN2 and EDN3 . EDNs are synthesised as large precursor proteins that are post-translationally cleaved to the biologically active 21-amino acid form . They are involved in fundamental cellular networks like cell proliferation, migration and differentiation processes [4, 5] by interacting with their corresponding cell surface-bound EDN-A (EDNRA) and EDN-B (EDNRB) receptors in an autocrine and also a paracrine manner [6–8]. A balanced regulation of this EDNRA/EDNRB interplay – also referred to as the endothelin axis (ET-axis)- is essential for, for example, homing processes to tissue destinations, where cells differentiate into numerous lineages such as the peripheral nervous system, structural and connective tissue components, cardiac cells or pigment-producing melanocytes .
There is now compelling evidence that imbalanced regulation of the ET-axis is implicated in human carcinogenesis, tumour progression and neo-angiogenesis [8, 10–12]. During malignant cell transformation, the basic tissue architecture, which is maintained by basement membrane delineation, becomes disrupted . This indicates the presence of crucial mediators that trigger the exchange of growth factors between the participating cells at the tumour invasion field. Essentially, such growth factor release is thought to enhance invasiveness, stimulate cell migration and promote neo-vascularisation . Multiple signal transduction pathways are affected downstream from EDNRA/B. In the case of interaction of EDNs with EDNRA, a pertussis toxin-insensitive G protein becomes activated and promotes stimulation of phospholipase C, resulting in the transactivation of the mitogen-activated protein kinase pathway . Second, EDN1 and EDN2 binding to EDNRA can activate p125 focal adhesion kinase and paxillin, both of which have been associated with increased tumour cell invasion. Moreover, EDNs are able to transduce the activation of anti-apoptotic signals through phosphatidylinositol-3-kinase and to stimulate neo-angiogenesis through vascular endothelial growth factor signalling . These multiple ET-axis pathway implications may explain its various impairments of normal cellular integrity in case of an aberrant shift from balanced to imbalanced EDN signalling.
Previously, EDN1 and EDN2 were found to be commonly overexpressed in a broad range of human tumour entities [8, 11, 12]. So far, most reports have focused on the role of EDN1 binding to EDNRA and its effects on tumour growth and neo-angiogenesis [8, 11, 13, 15]. A role similar to that of EDN1 has been described for EDN2 in human breast cancer. Increased expression of EDN1 and EDN2, but not of EDN3, induced chemotaxis of breast cancer cells and increased tumour cell invasion through the basement membrane , although conflicting results have been reported by others . In line with this, previous reports described a compensatory effect of EDN3 by negatively modulating the effects transduced by EDN1  and demonstrated that downregulation of EDN3 is associated with upregulation of EDN1 in human tissues .
However, a comprehensive analysis of EDN3 expression in normal and cancerous breast tissues and its potential implication in human breast cancer has not been published so far. In our study, we investigated for the first time EDN3 mRNA and protein expression in a large number of primary breast tissues and breast cell lines. Furthermore, we identified the molecular mechanism by which EDN3 expression is deregulated in breast carcinomas.
Materials and methods
Cryo-conserved clinical patient material
Cryo-conserved clinical samples were obtained from breast cancer patients treated by primary surgery at the University Hospitals of Aachen, Düsseldorf and Regensburg. Patients receiving neo-adjuvant chemotherapy and patients with recurrent breast cancer were excluded. Resected tissue was snap-frozen in liquid nitrogen immediately after surgery. Only samples containing more than 70% of tumourous cells in haematoxylin/eosin-stained control sections were further processed (n = 128). For 17 samples, macroscopically normal breast tissues containing at least 30% of epithelial cells were available. In all cases, two board-certified pathologists agreed on the diagnosis of breast cancer. Tumour histology was determined according to the criteria of the World Health Organization (2003), whereas disease stage was assessed according to the UICC (Union Internationale contre le Cancer) . Tumours were graded according to Bloom and Richardson, as modified by Elston and Ellis . All patients gave informed consent for retention and analysis of their tissue for research purposes, and the institutional review boards of the participating centres approved the study. For 98 patients, follow-up data were available with a median time of 63 months (range 1 to 124 months). Patient characteristics of this cohort are summarised in Table 1.
Formalin-fixed paraffin-embedded clinical patient material
A tissue microarray (TMA) was created as described previously by Bubendorf and colleagues . The formalin-fixed paraffin-embedded (FFPE) tissue sections were obtained from the archive of the Institute of Pathology, University of Regensburg, Germany. In all cases, two board-certified pathologists agreed on the diagnosis of breast cancer. Patients receiving neo-adjuvant chemotherapy and patients with recurrent breast cancer were excluded. All patients gave informed consent for retention and analysis of their tissues for research purposes, and the institutional review board of the participating centre approved the study. The TMA consisted of 150 primary tumours from malignant breast tissue and 44 normal breast specimens. Follow-up data were available for 146 patients with a median time of 77 months (range 1 to 148 months). Detailed tumour characteristics of this cohort are listed in Table 2.
Breast cell lines
The non-cancerous breast cell lines MCF10A and MCF12A as well as the cancerous breast cell lines MCF7, SKBR3, MDA-MB231 and BT20 were obtained from the American Type Culture Collection (Manassas, VA, USA) and cultured as recommended by the vendor.
Northern blot expression analysis
Expression of EDN3 mRNA in various human tissues was tested using the commercial multiple-tissue Northern (MTN) blots I and II (Clontech, Heidelberg, Germany), containing 2 μg of poly A+ RNA per lane from 16 different human tissues (that is, blot I: heart, whole brain, placenta, lung, liver, skeletal muscle, kidney and pancreas; blot II: spleen, thymus, prostate, testis, ovary, small intestine, colon [no mucosa] and peripheral blood lymphocytes). Hybridisation was performed using 25 ng of an EDN3-specific 722-base pair (bp) polymerase chain reaction (PCR) product derived from (GenBank accession number NM_000114.2) (position: 968 to 1,707), which was verified by sequence analysis. 32P-labelling of the DNA probe was achieved using the Megaprime DNA Labeling System (Amersham Biosciences, now part of GE Healthcare, Little Chalfont, Buckinghamshire, UK), and hybridisation was performed in accordance with the recommendation the manufacturer. The cancer profiling array (CPA) I (Clontech) is a matched tumour/normal expression array consisting of cDNA synthesised from 50 breast carcinomas, 50 normal breast tissues and 3 breast cancer lymph node metastasis specimens. Hybridisation was performed in accordance with the recommendations the manufacturer as described above for the MTN blots. Hybridisation signals of both, MTN blots and the CPA, were evaluated by use of a STORM-860 phosphoimager (Molecular Dynamics, now part of GE Healthcare). Intensity ratios were calculated after normalising signals against the background.
Nucleic acid extraction
Frozen tissue samples and cell line pellets were dissolved in lysis buffer for subsequent DNA isolation using the QIAmp DNA Mini kit (Qiagen, Hilden, Germany) or for total RNA isolation by using TRIzol reagent (Invitrogen Corporation, Carlsbad, CA, USA) in accordance with the protocol supplied by the manufacturers.
Reverse transcription of RNA
Of the extracted total RNA, 1 μg was reverse-transcribed using the Reverse Transcription System (Promega Corporation, Madison, WI, USA) by applying a mix of oligo-dT and pdN(6)-hexamer primers (1:2). The obtained cDNA was diluted (20 ng/μL) and test-amplified using intron-spanning primers for glyceraldehyde-3-phosphate-dehydrogenase (GAPDH). Primer sequences are provided in Table 3. PCRs were initiated as 'Hot Start' PCR at 95°C for 5 minutes and a hold at 80°C before the addition of 1 unit of GoTaq DNA polymerase (Promega Corporation). Cycle conditions were 95°C for 5 minutes, 35 cycles of 95°C for 1 minute, 60°C for 1 minute, 72°C for 1 minute and a final extension at 72°C for 10 minutes. PCR analyses were carried out in a PTC-200 cycler (Bio-Rad Laboratories, Inc., formerly MJ Research, Hercules, CA, USA). Amplificates were evaluated under ultraviolet light after 2% agarose gel electrophoresis containing ethidium bromide. Only samples yielding a specific 510-bp amplificate were further subjected to real-time PCR.
Semi-quantitative real-time polymerase chain reaction
The Roche LightCycler system was used for semi-quantitative light cycler analysis in combination with the LightCycler DNA Master SYBR Green I Kit (Roche, Mannheim, Germany) as previously described . Gene expression was quantified by the comparative cycle threshold (CT) method, normalising CT values to the housekeeping gene GAPDH and calculating relative expression values . A commercially available normal breast cDNA pool (Clontech) was used as a breast reference standard. Primer sequences are listed in Table 3.
Paraffin sections of 2 μm were deparaffinised in xylene followed by rehydration in a decreasing ethanol series. Antigen retrieval was performed by pre-treatment in boiling citrate buffer (pH 6.0) in a microwave oven for 30 minutes (200 W). Immunohistochemistry (IHC) was performed using an NEXES Immuno Stainer (Ventana Medical Systems, Inc., Tucson, AZ, USA) in accordance with the specifications of the manufacturer. A goat polyclonal EDN3-specific antibody (sc-21628; Santa Cruz Biotechnology, Inc., Santa Cruz, CA, USA) was used in a 1:150 dilution by using the ChemMate Envision Kit (DAKO, Hamburg, Germany). Counterstaining was performed by using Mayer's haematoxylin. The incubation with primary antibodies was omitted in negative controls. All analysed samples were stained without the knowledge of histopathological data. Cytoplasmic protein staining was semi-quantitatively scored by an experienced breast pathologist according to the well-established scoring system developed by Remmele and Stegner . To verify staining specificity, the primary antibody was incubated with a 200 molar excess of blocking peptide (sc-21628 P; Santa Cruz Biotechnology, Inc.) for 2 hours prior to its application on test samples.
CpG island prediction
EDN1, EDN2 and EDN3 genomic nucleotide sequences were taken from the Ensembl database and analysed for promoter CpG (cytosine-phosphate-guanine dinucleotide) islands in accordance with the method of Li and Dahiya . A fragment of 2 kb in size, beginning 1 kb 5'-upstream from the annotated transcription start site (TSS) and ending 1 kb 3'-downstream, was analysed by applying the following criteria: island size of greater than 200 bp, guanine/cytosine content of greater than 60% and observed/expected CpG ratio of greater than 0.6.
Bisulphite modification and methylation-specific polymerase chain reaction
Of the genomic DNA, 1 μg was bisulphite-modified using the EZ DNA Methylation Kit (Zymo Research Corporation, Orange, CA, USA) and eluted in 20 μL of Tris buffer (10 mM). Methylation-specific PCR (MSP) was performed in accordance with the method of Herman and colleagues . One microlitre of bisulphite-treated DNA was amplified using MSP primers that specifically recognise either the unmethylated or methylated EDN3 promoter sequence after bisulphite conversion (Table 3). To achieve high accuracy, each primer was designed to cover three CpG sites of template DNA. Commercially available universal poly-methylated DNA and unmethylated DNA (EpiTect Control DNA; Qiagen) were used as positive controls for methylated and unmethylated EDN3 sequences, respectively. MSP products were visualised under ultraviolet light after 3% low-range ultra agarose gel electrophoresis containing ethidium bromide (Bio-Rad Laboratories, Inc.). Promoter methylation status was interpreted in a binary qualitative fashion.
In vitrodemethylating treatment
Cells were seeded at a density of 3 × 104 cells/cm2 in a six-well plate. The demethylation agent 5-aza-2'-deoxycytidine (DAC) (Sigma-Aldrich, Deisenheim, Germany) was added to a final concentration of 1 μM in fresh medium at days 1, 2 and 3 after seeding. Additionally, cells were exposed to 300 nM trichostatin A (TSA) (Sigma-Aldrich) on day 3 for 24 hours. Control cells without DAC/TSA were supplied with fresh medium on days 1, 2 and 3. DNA and RNA were extracted on day 4 as mentioned above.
SPSS version 14.0 (SPSS Inc., Chicago, IL, USA) was used for statistical analyses. All tests were performed two-tailed, and P values of below 0.05 were considered statistically significant. The non-parametric Mann-Whitney U test and the Student t test (paired and unpaired) were used to compare expression results between cancer tissues and normal tissues or between EDN3 mRNA expression and EDN3 methylation status. Contingency table analysis and Fisher exact tests were used to study the statistical association between clinicopathological factors and EDN3 protein expression or EDN3 promoter methylation status. Survival curves comparing patients with or without any of the factors were calculated using the Kaplan-Meier method, with significance evaluated by log-rank statistics. Breast cancer-specific survival (CSS) was measured from the day of surgery until tumour-related death and was censored for patients alive at last contact or in case of death unrelated to the tumour. Disease-free survival (DFS) was measured from surgery until disease relapse and censored for patients alive without evidence of relapse at the last follow-up. For EDN3 protein expression, a multivariate Cox proportional hazard model was employed to assess the relative risks on patient CSS and to test for independent prognostic relevance of clinical/investigational factors. Only patients for whom the status of all selected variables was known were included in the proportional hazard model (n = 121). The limit for reverse-selection procedures was P = 0.2. The proportionality assumption for all variables was assessed with log-negative-log survival distribution functions. For analyses, the following variables were categorised into binary values: small-sized (pT1) versus large-sized (pT2 to pT4), node-negative (pN0) versus node-positive (pN1 to pN3) and low-grade (G1 and G2) versus high-grade (G3) tumours.
EDN3mRNA is differentially expressed in human breast cancer
As previously reported, we performed in silico Northern blot analysis and RNA array-based expression profiling to identify novel candidate genes differentially expressed in human breast cancer [27, 28]. In the latter approach, EDN3 was detected as one of the most frequently downregulated genes, showing substantial expression loss in 63% of breast carcinomas (data not shown). Therefore, we started a detailed analysis on EDN3 expression and its potential implication in human breast cancer.
An initial Northern blot analysis of normal human tissues demonstrates that EDN3 mRNA is abundantly expressed in a variety of non-malignant tissues, including pancreas, spleen, prostate, testis, small intestine and colon, with two major mRNA transcripts of 2.4 and 2.7 kb in size (Figure 1a). A breast cancer dot blot cDNA array was then hybridised with the same EDN3-specific probe. This showed a clear loss of EDN3 mRNA expression, as defined by a fold change (tumour versus normal) of greater than 2 (FC2), in 96% (48 of 50) of the analysed breast carcinoma samples and also in all three corresponding lymph node metastases (P < 0.001) (Figure 1b, c). To confirm this result, EDN3 mRNA expression was also assessed by real-time PCR in a set of 77 breast tumour tissues and 17 corresponding normal breast tissues. In breast carcinomas, downregulation of EDN3 expression by FC2 could be detected in 60 of 77 cases (78%; median: 21-fold) as compared with the normal breast reference standard (Figure 1d). This downregulation was still evident when comparing EDN3 expression among all tumour and normal breast tissues as illustrated by box plot analysis (P < 0.001) (Figure 1e) as well as when comparing fold changes of EDN3 expression among the 17 matched pairs (Figure 1f). In the latter analysis, 14 of 17 pairs (82%) showed EDN3 downregulation by FC2, with a median expression change of 13-fold.
Differential EDN3 protein expression in human breast cancer
To analyse whether loss of EDN3 expression in breast cancer is also evident on the protein level, we used a TMA comprising 150 invasive breast carcinomas and 44 normal breast tissue specimens. The specificity of the antibody applied was determined by the simultaneous use of blocking peptide against EDN3 in normal breast samples, which showed a clear decrease of overall EDN3 staining (Additional data file 1). EDN3 protein was clearly detectable in 75.0% (32 of 44) of normal breast tissue samples analysed (Figure 2a, b), as defined by an immunoreactivity score (IRS) of at least 8. The mean normal expression was determined to be IRS = 8.2 (± 3.8 standard deviations [SDs]). Expression was predominantly localised in luminal and basal epithelial cells of the normal breast and was weakly detectable in normal stromal compartments. In contrast, invasive breast carcinomas showed complete loss or reduced EDN3 expression (IRS < 8) in 45.3% (68 of 150) of cases (Figure 2c, d) and abundant EDN3 expression in 54.7% (82 of 159) (Figure 2e, f). Mean EDN3 expression in invasive breast carcinomas was determined to be IRS = 6.7 (± 4.0 SDs). EDN3 protein was rarely observed in tumour stroma (that is, was detectable only in those stromal cells adjacent to tumour cells with strong EDN3 expression). The difference of EDN3 expression between tumours and normal breast tissues was statistically significant (Mann-Whitney U test: P = 0.037; Student unpaired t test: P = 0.039).
Loss of EDN3 expression is associated with adverse patient outcome in human breast cancer
To investigate a potential clinical relevance of EDN3 in breast cancer, we first analysed whether EDN3 protein expression is associated with clinicopathological parameters. In a bivariate analysis, differential EDN3 expression was not associated with patient age at diagnosis, tumour size, lymph node metastasis, histological grade, histological type, tumour focality or oestrogen or progesterone receptor status (Table 4). Interestingly, in a univariate survival analysis, we found a significant association between low EDN3 expression and unfavourable CSS (P = 0.022) (Figure 2g and Table 5). Patients with low EDN3 expression in the tumour had a mean CSS of 99 months (95% confidence interval [CI] 85 to 113 months) as compared with patients with tumours showing high EDN3 expression, which had a prolonged mean CSS of 116 months (95% CI 106 to 126 months). The visual impression of the Kaplan-Meier curves suggests also an association of EDN3 expression with DFS (Figure 2h), but this was statistically not significant (P = 0.167) (Table 5). Next, we performed a multivariate Cox regression analysis to test for independent significance of EDN3 expression as a prognostic factor in patient CSS (Table 6). Patient age at diagnosis, tumour size, lymph node metastasis, histological grade and EDN3 expression were included in the model. After reverse selection, patient age, lymph node status, grade and EDN3 expression remained significant in the Cox model, with low EDN3 expression displaying a twofold elevated risk of dying from breast cancer (hazard ratio 1.98, 95% CI 1.09 to 3.61; P = 0.025).
Loss of EDN3mRNA expression in breast cancer cell lines
The NH2-terminal peptide structure of EDN3 differs considerably in essential amino acids between residues 2 and residues 4 to 7 from those of EDN1 and EDN2 (Figure 3a). This region bulges out of the basic EDN structure and has been reported to represent a major domain for binding specificity to EDN receptors and thus is critical in EDN signalling activity . Because we found that EDN3 expression is frequently lost in breast cancer tissues, we next compared EDN1, EDN2 and EDN3 mRNA expression in breast cell lines in parallel. EDN1 was much more strongly expressed in all cancerous cell lines (MCF7, SKBR3, MDA-MB231 and BT20) than in non-cancerous MCF12A cells (Figure 3b), showing a mean upregulation of 13.9-fold (± 7.8 SDs). EDN2 was more strongly expressed in SKBR3 cells but less expressed in MCF7, MDA-MB231 and BT20 as compared with MCF12A cells. EDN3 expression, however, was clearly downregulated in all cancerous cell lines as compared with non-cancerous cells, revealing a mean expression of 0.008 (± 0.009 SDs) of that found in MCF12A cells (set to 1). Therefore, upregulation of EDN1 appears to coincide with downregulation of EDN3 in breast cancer cell lines.
Methylation of the EDN3promoter in breast cancer cell lines
Since promoter hypermethylation is responsible for transcriptional silencing of important tumour suppressor genes in various human cancer types , we searched all three EDN genes for the presence of CpG islands in their promoter region. A region of high CpG density in the EDN3 nucleotide sequence was identified as a CpG island (Figure 4a), whereas the CpG density in EDN1 and EDN2 is lower and does not define a CpG island under the selected criteria. To analyse the methylation status of the EDN3 CpG promoter (Figure 4b), we performed MSP with DNA after bisulphite treatment from non-malignant human tissues and three non-malignant and four malignant breast cell lines. The utilised MSP primers were tested on a dilution series of poly-methylated DNA, which revealed a sensitivity of 0.01 (1%) in detecting methylated EDN3 DNA molecules in a background of unmethylated EDN3 DNA molecules (Figure 4c). We observed no methylation of the EDN3 promoter in HMEC cells, human placental tissue, peripheral blood lymphocytes (Figure 4c) or non-malignant MCF10A cells (Figure 4e). A weak methylation signal was detected in non-malignant MCF12A cells. Of the malignant breast cell lines, MDA-MB231 and MCF7 harboured a methylated EDN3 promoter. In BT20 cells, both unmethylated and methylated EDN3 promoter sequences could be detected, whereas EDN3 was unmethylated in SKBR3 cells. Next, we analysed the association between EDN3 expression and promoter methylation in six breast cell lines by in vitro demethylating their DNA and assessing EDN3 mRNA expression after the treatment. A clear conversion of methylation could be observed in cell lines that were originally methylated in the EDN3 promoter region (that is, in MDA-MB231 and MCF7 cells) (Figure 4e), resulting in 47-fold and 28-fold increases, respectively, in EDN3 mRNA expression (Figure 4f). In BT20 cells, however, the demethylating effect was weaker and led to a 3-fold induction of EDN3 transcription. In weakly methylated MCF12A cells, the conversion of methylated alleles induced EDN3 expression by 5-fold, whereas no substantially altered change of EDN3 expression could be detected in unmethylated MCF10A (1.6-fold) or SKBR3 (1.0-fold) cells.
Frequent EDN3promoter methylation in primary breast carcinomas
Next, we analysed EDN3 promoter methylation in primary breast cancer as well. In total, 89 of 128 breast carcinoma samples (69.5%) showed EDN3 promoter methylation (for example, #5 in Figure 5a). The remaining breast carcinoma samples (39 of 128; 30.5%) were not affected by this epigenetic modification (for example, #13 in Figure 5b). None of the 17 normal breast tissues exhibited EDN3 promoter methylation. Cancerous tissues yielded a PCR product with primers specific for the unmethylated EDN3 promoter sequence in all cases, due to non-malignant contaminants (stromal cells and endothelial cells) present in the bulk tumour tissue, as has also been described by Suzuki and colleagues .
We next aimed to analyse whether EDN3 promoter methylation was associated with EDN3 mRNA expression in these tissues. For 71 breast cancer specimens, both EDN3 methylation and EDN3 mRNA expression have been investigated in parallel. Figure 5b shows the distribution of EDN3 mRNA expression among these two groups. While carcinoma samples with unmethylated EDN3 promoter exhibited similar EDN3 mRNA expression as compared with normal breast tissues (Figure 1e), breast carcinomas with EDN3 methylation exhibited a significant downregulation of EDN3 mRNA expression as compared with EDN3 unmethylated samples (P = 0.005, Mann-Whitney U test).
Finally, we asked whether EDN3 promoter methylation may be of clinical relevance in human breast cancer, as we have previously found for EDN3 protein expression. In a univariate analysis, the EDN3 methylation status in breast carcinomas was not associated with patient age at diagnosis, tumour size, lymph node metastasis, histological grade, histological type or oestrogen or progesterone receptor positivity (Table 7). In contrast to EDN3 protein expression, EDN3 promoter methylation was significantly associated neither with patient CSS (P = 0.703) nor with patient DFS (P = 0.632) (data not shown).
The involvement of EDNs in tumourigenesis has been described in several reports [4, 10, 15]. In contrast to the potentially oncogenic role of EDN1 and EDN2, there is still little knowledge about the role of EDN3 in cancer initiation or progression. A recent study demonstrated abundant expression of EDN1 and EDN2 but complete absence of EDN3 expression in a representative set of human breast cancer cell lines . Because we have previously found that EDN3 mRNA expression is downregulated in primary breast carcinomas as compared with normal breast tissues [27, 28], we aimed in this report to provide the first comprehensive analysis of EDN3 expression and its potential implication in human breast cancer.
Initially, we screened various non-malignant epithelial tissues for EDN3 mRNA expression and also analysed its expression using a breast cancer cDNA dot blot array. Besides abundant expression in several human tissues, EDN3 was strongly expressed in normal breast samples, providing evidence for a functional role in epithelial tissues such as the mammary gland. In contrast, most matched breast carcinomas showed diminished EDN3 mRNA expression both on the cDNA dot blot array and by real-time PCR analysis. This finding supports the current evidence that EDN3 may exert a functional role divergent to that of EDN1/EDN2 in the human mammary gland . A further TMA analysis revealed that EDN3 protein is abundantly expressed in normal breast whereas its expression is reduced in a large fraction of breast carcinomas. Frequency differences may arise due to the use of different techniques (real-time PCR versus IHC) on separate tumour cohorts (fresh frozen versus FFPE) and accomplishing different scoring systems. Loss of EDN3 protein expression was not associated with relevant clinicopathological factors. For instance, it occurred with almost equal frequency among all tumour sizes (pT1 to pT4), suggesting that it may be an early event in the development of infiltrating breast carcinoma. Since EDN3 is thought to counterbalance the effects mediated by EDN1 and EDN2 [4, 18], we propose that loss of EDN3 expression could actively enhance overexpression of the ET-axis. Recently, upregulation of ET-axis members was found to be associated with higher histological grade, lymph node metastasis and lymphovascular invasion in breast cancer  and also with advanced tumour progression in ovarian cancer , prostate cancer , Ewing sarcoma and neuroblastoma . A systematic expression analysis on larger breast carcinoma cohorts and metastatic deposits is now required, including all three EDNs and EDNRA/EDNRB. This will unravel in detail the inter-relationship between EDN3 expression loss and upregulation of EDN1/2 and EDNRA/B as well as its association with breast tumour progression. In our study, loss of EDN3 expression was associated with adverse patient outcome. So far, overexpressions of EDN1 and EDNRA were already reported as being associated with impaired survival in breast cancer [5, 34]. Our findings support the view that an imbalanced ET-axis is of pivotal relevance in breast cancer biology and that EDN3, unlike other members of the ET-axis, may represent a novel tumour suppressor gene in the human mammary gland.
Addressing the molecular cause by which EDN3 expression becomes abrogated, we found that the EDN3 gene promoter, unlike EDN1 and EDN2, contains a CpG island as a potential substrate to aberrant hypermethylation and consequently gene inactivation. Indeed, we detected EDN3 promoter methylation in cancerous breast cell lines in functional association with loss of EDN3 mRNA expression. Moreover, a hypermethylated EDN3 promoter was also detected in 70% of breast carcinoma specimens in significant association with loss of EDN3 expression. We therefore conclude that aberrant EDN3 methylation is a tumour-specific event and the predominant mechanism leading to EDN3 expression loss in breast cancer. However, it remains elusive why patient survival was not associated with EDN3 methylation as it was with loss of EDN3 protein expression. In fact, only very few studies detected such outcome association on both molecular levels (for example, for SFRP1 [22, 35] or ITIH5 ), probably due to considerable sensitivity differences of the available detection techniques as well as further genetic or epigenetic alterations contributing to the loss of a gene's expression. Interestingly, EDNRB was previously described to be methylated in numerous tumour entities, such as lung, colon, prostate, bladder, kidney, liver, oesophageal, nasopharyngeal cancer and leukemia [37–45], but to the authors' knowledge, never in gynaecological tumours. So far, there has been no evidence that EDNRB becomes methylated in breast carcinomas since a previous study demonstrated strong EDNRB expression in all invasive ductal carcinoma samples and in all analysed cancerous breast cell lines . Notably, an ET-axis expression pattern similar to that of breast cancer was recently found in cervical cancer; that is, upregulation of EDN1, EDN2, EDNRA and EDNRB expression was accompanied by downregulation of EDN3 expression in cancerous cervix as compared with normal cervical epithelium . This suggests that a decrease of EDN3 expression accompanied by an increase of EDNRB expression may be a particular feature of gynaecological tumour entities.
Since the ET-axis represents crucial decisive elements for the direction of tumour growth, invasion and neo-angiogenesis, it provides a promising intervention point for molecular targeted therapies. EDNR antagonists have been proven as potent and specific ET-axis inhibitors that block cellular pathways implicated in tumour growth. The drug bosentan, which targets both EDNRA and EDNRB, inhibits tumour growth, vascularisation and bone metastasis in breast cancer . Atrasenatan, targeting EDNRA, is capable of inhibiting proliferation and cancer growth-promoting processes [48, 49]. In addition, blockers of EDNRA resensitised cancer cells to paclitaxel-induced apoptosis, as observed in ovarian, prostatic, cervical and nasopharyngeal cancer cell lines [49–51] as well as in primary ovarian and breast cancer [5, 52]. Our study adds a novel aspect to therapeutically targeting the ET-axis in breast cancer. Since the epigenetic lock of the EDN3 gene is potentially reversible by DNA methyltransferase (DNMT) or histone acetyltransferase (HDAC) inhibitors or both, these drug classes may provide a future option in a combined treatment consisting of a decrease in EDN1/2 signalling by blocking EDNRs together with the reactivation of EDN3 expression by DNMT and HDAC inhibitors. Apparently, molecular rebalancing of the ET-axis in cancerous cells may be most efficiently achieved by targeting all deregulated axis molecules.
In summary, our study contributes to the understanding of deregulated EDN signalling commonly observed in human tumours. EDN3 expression, in contrast to abundant EDN1 and EDN2 expression, becomes frequently inactivated by promoter methylation in human breast cancer, potentially causing aberrant activation of the ET-axis, which in turn may promote this disease. We therefore conclude that EDN3 may be an interesting future target for an epigenetic therapy. Forced EDN3 re-expression by DNA demethylation agents in conjunction with inhibitors of EDNRs may rebalance ET-axis-mediated cellular signalling to a more normal status, thus having a therapeutic impact in human breast cancer.
cancer profiling array
breast cancer-specific survival
fold change of greater than 2
methylation-specific polymerase chain reaction
multiple-tissue Northern (blot)
polymerase chain reaction
Grimshaw MJ: Endothelins and hypoxia-inducible factor in cancer. Endocr Relat Cancer. 2007, 14: 233-244. 10.1677/ERC-07-0057.
Levin ER: Endothelins. N Engl J Med. 1995, 333: 356-363. 10.1056/NEJM199508103330607.
Masaki T: The endothelin family: an overview. J Cardiovasc Pharmacol. 2000, 35: S3-S5.
Grimshaw MJ, Hagemann T, Ayhan A, Gillet CE, Binder C, Balkwill FR: A role for endothelin-2 and its receptors in breast tumor cell invasion. Cancer Res. 2004, 64: 2461-2468. 10.1158/0008-5472.CAN-03-1069.
Smollich M, Wulfing P: The endothelin axis: a novel target for pharmacotherapy of female malignancies. Curr Vasc Pharmacol. 2007, 5: 239-248. 10.2174/157016107781024082.
Pollock DM, Keith TL, Highsmith RF: Endothelin receptors and calcium signaling. FASEB J. 1995, 9: 1196-1204.
Bagnato A, Catt KJ: Endothelin as autocrine regulators of tumor cell growth. Trends Endocrinol Metab. 1998, 9: 378-383. 10.1016/S1043-2760(98)00094-0.
Bagnato A, Spinella F, Rosano L: Emerging role of the endothelin axis in ovarian tumour progression. Endocr Relat Cancer. 2005, 12: 761-772. 10.1677/erc.1.01077.
Pla P, Larue L: Involvement of endothelin receptors in normal and pathological development of neural crest cells. Int J Dev Biol. 2003, 47: 315-325.
Alanen K, Deng DX, Chakrabarti S: Augmented expression of endothelin-1, endothelin-3 and the endothelin-B receptor in breast carcinoma. Histopathology. 2000, 36: 161-167. 10.1046/j.1365-2559.2000.00795.x.
Berry P, Burchill S: Endothelins may modulate invasion and proliferation of Ewing's sarcoma and neuroblastoma. Clin Sci (Lond). 2002, 103 (Suppl 48): 322S-326S.
Grimshaw MJ: Endothelins in breast tumour cell invasion. Cancer Lett. 2005, 222: 129-138. 10.1016/j.canlet.2004.08.029.
Bagnato A, Natali PG: Endothelin receptors as novel targets in tumor therapy. J Transl Med. 2004, 2: 16-10.1186/1479-5876-2-16.
Spinella F, Rosano L, Di Castro V, Natali PG, Bagnato A: Endothelin-1 induces vascular endothelial growth factor by increasing hypoxia-inducible factor-1alpha in ovarian carcinoma cells. J Biol Chem. 2002, 277: 27850-27855. 10.1074/jbc.M202421200.
Nelson JB, Chan-Tack K, Hedican SP, Magnuson SR, Opgenorth TJ, Bova GS, Simons JW: Endothelin-1 production and decreased endothelin B receptor expression in advanced prostate cancer. Cancer Res. 1996, 56: 663-668.
Hagemann T, Binder C, Binder L, Pukrop T, Trümper L, Grimshaw MJ: Expression of endothelins and their receptors promotes an invasive phenotype of breast tumor cells but is insufficient to induce invasion in benign cells. DNA Cell Biol. 2005, 24: 766-776. 10.1089/dna.2005.24.766.
Mitaka C, Hirata Y, Ichikawa K, Yokoyama K, Emori T: Effects of TNF-alpha on hemodynamic changes and circulating endothelium-derived vasoactive factors in dogs. Am J Physiol. 1994, 267: H1530-H1536.
Thomson E, Kumarathasan P, Vincent R: Pulmonary expression of preproET-1 and preproET-3 mRNAs is altered re ciprocally in rats after inhalation of air pollutants. Exp Biol Med (Maywood). 2006, 231: 979-984.
TNM Classification of Malignant Tumours. Edited by: Sobin LH, Wittekind C. 2002, New York: Wiley-Liss, 6
Elston EW, Ellis IO: Method for grading breast cancer. J Clin Pathol. 1993, 46: 189-190. 10.1136/jcp.46.2.189-b.
Bubendorf L, Nocito A, Moch H, Sauter G: Tissue microarray (TMA) technology: miniaturized pathology archives for high-throughput in situ studies. J Pathol. 2001, 195: 72-79. 10.1002/path.893.
Veeck J, Niederacher D, An H, Klopocki E, Wiesmann F, Betz B, Galm O, Camara O, Dürst M, Kristiansen G, Huszka C, Knüchel R, Dahl E: Aberrant methylation of the Wnt antagonist SFRP1 in breast cancer is associated with unfavourable prognosis. Oncogene. 2006, 25: 3479-3488. 10.1038/sj.onc.1209386.
Fink L, Seeger W, Ermert L, Hänze J, Stahl U, Grimminger F, Kummer W, Bohle RM: Real-time quantitative RT-PCR after laser-assisted cell picking. Nat Med. 1998, 4: 1329-1333. 10.1038/3327.
Remmele W, Stegner HE: Recommendation for uniform definition of an immunoreactive score (IRS) for immunohistochemical estrogen receptor detection (ER-ICA) in breast cancer tissue. Pathologe. 1987, 8: 138-140.
Li LC, Dahiya R: MethPrimer: designing primers for methylation PCRs. Bioinformatics. 2002, 18: 1427-1431. 10.1093/bioinformatics/18.11.1427.
Herman JG, Graff JR, Myöhänen S, Nelkin BD, Baylin SB: Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci USA. 1996, 93: 9821-9826. 10.1073/pnas.93.18.9821.
Dahl E, Sadr-Nabavi A, Klopocki E, Betz B, Grube S, Kreutzfeld R, Himmelfarb M, Gelling S, Klaman I, Hinzmann B, Kristiansen G, Grützmann R, Kuner R, Petschke B, Rhiem K, Wiechen K, Sers C, Wiestler O, Schneider A, Höfler H, Nährig J, Dietel M, Schäfer R, Rosenthal A, Schmutzler R, Dürst M, Meindl A, Niederacher D: Systematic identification and molecular characterization of genes differentially expressed in breast and ovarian cancer. J Pathol. 2005, 205: 21-28.
Dahl E, Kristiansen G, Gottlob K, Klaman I, Ebner E, Hinzmann B, Hermann K, Pilarsky C, Dürst M, Klinkhammer-Schalke M, Blaszyk H, Knuechel R, Hartmann A, Rosenthal A, Wild PJ: Molecular profiling of laser-microdissected matched tumor and normal breast tissue identifies karyopherin alpha2 as a potential novel prognostic marker in breast cancer. Clin Cancer Res. 2006, 12: 3950-3960. 10.1158/1078-0432.CCR-05-2090.
Wallace BA, Janes RW, Bassolino DA, Krystek SR: A comparison of X-ray and NMR structures for human endothelin-1. Protein Sci. 1995, 4: 75-83.
Esteller M: Epigenetic gene silencing in cancer: the DNA hypermethylome. Hum Mol Genet. 2007, 16 (Spec No 1): R50-59. 10.1093/hmg/ddm018.
Suzuki H, Watkins DN, Jair KW, Schuebel KE, Markowitz SD, Chen WD, Pretlow TP, Yang B, Akiyama Y, Van Engeland M, Toyota M, Tokino T, Hinoda Y, Imai K, Herman JG, Baylin SB: Epigenetic inactivation of SFRP genes allows constitutive WNT signaling in colorectal cancer. Nat Genet. 2004, 36: 417-422. 10.1038/ng1330.
Donninger H, Bonome T, Radonovich M, Pise-Masison CA, Brady J, Shih JH, Barrett JC, Birrer MJ: Whole genome expression profiling of advance stage papillary serous ovarian cancer reveals activated pathways. Oncogene. 2004, 23: 8065-8077. 10.1038/sj.onc.1207959.
Kopetz ES, Nelson JB, Carducci MA: Endothelin-1 as a target for therapeutic intervention. Invest New Drugs. 2002, 20: 173-182. 10.1023/A:1015630513908.
Wulfing P, Diallo R, Kersting C, Wülfing C, Poremba C, Rody A, Greb RR, Böcker W, Kiesel L: Expression of endothelin-1, endothelin-A, and endothelin-B receptor in human breast cancer and correlation with long-term follow-up. Clin Cancer Res. 2003, 9: 4125-4131.
Klopocki E, Kristiansen G, Wild PJ, Klaman I, Castanos-Velez E, Singer G, Stöhr R, Simon R, Sauter G, Leibiger H, Essers L, Weber B, Hermann K, Rosenthal A, Hartmann A, Dahl E: Loss of SFRP1 is associated with breast cancer progression and poor prognosis in early stage tumors. Int J Oncol. 2004, 25: 641-649.
Veeck J, Chorovicer M, Naami A, Breuer E, Zafrakas M, Bektas N, Dürst M, Kristiansen G, Wild PJ, Hartmann A, Knuechel R, Dahl E: The extracellular matrix protein ITIH5 is a novel prognostic marker in invasive node-negative breast cancer and its aberrant expression is caused by promoter hypermethylation. Oncogene. 2008, 27: 865-875. 10.1038/sj.onc.1210669.
Chen SC, Lin CY, Chen YH, Fang HY, Cheng CY, Chang CW, Chen RA, Tai HL, Lee CH, Chou MC, Lin TS, Hsu LS: Aberrant promoter methylation of EDNRB in lung cancer in Taiwan. Oncol Rep. 2006, 15: 167-172.
Pao MM, Tsutsumi M, Liang G, Uzvolgyi E, Gonzales FA, Jones PA: The endothelin receptor B (EDNRB) promoter displays heterogeneous, site specific methylation patterns in normal and tumor cells. Hum Mol Genet. 2001, 10: 903-910. 10.1093/hmg/10.9.903.
De Marzo AM, Isaacs WB, Nelson WG: Hypermethylation of CpG islands in primary and metastatic human prostate cancer. Cancer Res. 2004, 64: 1975-1986. 10.1158/0008-5472.CAN-03-3972.
Friedrich MG, Weisenberger DJ, Cheng JC, Chandrasoma S, Siegmund KD, Gonzalgo ML, Toma MI, Huland H, Yoo C, Tsai YC, Nichols PW, Bochner BH, Jones PA, Liang G: Detection of methylated apoptosis-associated genes in urine sediments of bladder cancer patients. Clin Cancer Res. 2004, 10: 7457-7465. 10.1158/1078-0432.CCR-04-0930.
Pflug BR, Zheng H, Udan MS, D'Antonio JM, Marshall FF, Brooks JD, Nelson JB: Endothelin-1 promotes cell survival in renal cell carcinoma through the ET(A) receptor. Cancer Lett. 2007, 246: 139-148. 10.1016/j.canlet.2006.02.007.
Hsu LS, Lee HC, Chau GY, Yin PH, Chi CW, Lui WY: Aberrant methylation of EDNRB and p16 genes in hepatocellular carcinoma (HCC) in Taiwan. Oncol Rep. 2006, 15: 507-511.
Zhao BJ, Sun DG, Zhang M, Tan SN, Ma X: Identification of aberrant promoter methylation of EDNRB gene in esophageal squamous cell carcinoma. Dis Esophagus. 2009, 22: 55-61. 10.1111/j.1442-2050.2008.00848.x.
Lo KW, Tsang YS, Kwong J, To KF, Teo PM, Huang DP: Promoter hypermethylation of the EDNRB gene in nasopharyngeal carcinoma. Int J Cancer. 2002, 98: 651-655. 10.1002/ijc.10271.
Hsiao PC, Liu MC, Chen LM, Tsai CY, Wang YT, Chen J, Hsu LS: Promoter methylation of p16 and EDNRB gene in leukemia patients in Taiwan. Chin J Physiol. 2008, 51: 27-31.
Sun de J, Liu Y, Lu DC, Kim W, Lee JH, Maynard J, Deisseroth A: Endothelin-3 growth factor levels decreased in cervical cancer compared with normal cervical epithelial cells. Hum Pathol. 2007, 38: 1047-1056. 10.1016/j.humpath.2006.12.015.
Dréau D, Karaa A, Culberson C, Wyan H, McKillop IH, Clemens MG: Bosentan inhibits tumor vascularization and bone metastasis in an immunocompetent skin-fold chamber model of breast carcinoma cell metastasis. Clin Exp Metastasis. 2006, 23: 41-53. 10.1007/s10585-006-9016-z.
Nelson J, Bagnato A, Battistini B, Nisen P: The endothelin axis: emerging role in cancer. Nat Rev Cancer. 2003, 3: 110-116. 10.1038/nrc990.
Rosanò L, Di Castro V, Spinella F, Nicotra MR, Natali PG, Bagnato A: ZD4054, a specific antagonist of the endothelin A receptor, inhibits tumor growth and enhances paclitaxel activity in human ovarian carcinoma in vitro and in vivo. Mol Cancer Ther. 2007, 6: 2003-2011. 10.1158/1535-7163.MCT-07-0151.
Bagnato A, Cirilli A, Salani D, Simeone P, Muller A, Nicotra MR, Natali PG, Venuti A: Growth inhibition of cervix carcinoma cells in vivo by endothelin A receptor blockade. Cancer Res. 2002, 62: 6381-6384.
Akhavan A, McHugh KH, Guruli G, Bies RR, Zamboni WC, Strychor SA, Nelson JB, Pflug BR: Endothelin receptor A blockade enhances taxane effects in prostate cancer. Neoplasia. 2006, 8: 725-732. 10.1593/neo.06388.
Jazaeri AA, Awtrey CS, Chandramouli GV, Chuang YE, Khan J, Sotiriou C, Aprelikova O, Yee CJ, Zorn KK, Birrer MJ, Barrett JC, Boyd J: Gene expression profiles associated with response to chemotherapy in epithelial ovarian cancers. Clin Cancer Res. 2005, 11: 6300-6310. 10.1158/1078-0432.CCR-04-2682.
Ortega Mateo A, De Artiñano AA: Highlights on endothelins: a review. Pharmacol Res. 1997, 36: 339-351. 10.1006/phrs.1997.0246.
The excellent technical assistance from Sonja von Serényi, Sevim Alkaya and Inge Losen is greatly appreciated. We thank Dieter Niederacher (Heinrich-Heine University, Düsseldorf, Germany) and Matthias Dürst (Friedrich-Schiller University, Jena, Germany) for kindly providing patient samples. We thank Monika Klinkhammer-Schalke and Felicitas Horn (Tumour Registry, Regensburg, Germany) for continuous help in obtaining clinical patient information and follow-up data. Immunohistochemical scoring was kindly performed by Nuran Bektas (Institute of Pathology, University Hospital of the RWTH, Aachen, Germany). HMEC DNA was a generous gift by Bernhard Lüscher (Institute of Biochemistry, University Hospital of the RWTH, Aachen, Germany). This work is a research project within the German Human Genome Project and was supported by a grant from the Bundesministerium für Bildung und Forschung (BMBF) to ED (01KW040-1).
The authors declare that they have no competing interests.
FW participated in the conception and design of the study and carried out the gene expression analyses, immunohistochemical studies and methylation experiments. JV performed statistical evaluations, participated in data interpretation and wrote the manuscript. OG provided expertise in DNA methylation analysis and critically revised the manuscript. AH provided clinical samples and clinicopathological data, performed data interpretation and critically revised the manuscript. ME participated in data interpretation and critically revised the manuscript. RK participated in the design and coordination of the study and critically revised the manuscript. ED planned and coordinated the study and critically revised the manuscript. All authors read and approved the final manuscript.
Frank Wiesmann, Jürgen Veeck contributed equally to this work.
Electronic supplementary material
Additional data file 1: A PDF file that demonstrates the specificity of the applied EDN3 antibody by use of competitive blocking peptide in immunohistochemistry. (PDF 1 MB)
About this article
Cite this article
Wiesmann, F., Veeck, J., Galm, O. et al. Frequent loss of endothelin-3 (EDN3) expression due to epigenetic inactivation in human breast cancer. Breast Cancer Res 11, R34 (2009). https://doi.org/10.1186/bcr2319