Skip to main content
  • Research article
  • Open access
  • Published:

A polymorphism at the 3'-UTR region of the aromatase gene defines a subgroup of postmenopausal breast cancer patients with poor response to neoadjuvant letrozole



Aromatase (CYP19A1) regulates estrogen biosynthesis. Polymorphisms in CYP19A1 have been related to the pathogenesis of breast cancer (BC). Inhibition of aromatase with letrozole constitutes the best option for treating estrogen-dependent BC in postmenopausal women. We evaluate a series of polymorphisms of CYP19A1 and their effect on response to neoadjuvant letrozole in early BC.


We analyzed 95 consecutive postmenopausal women with stage II-III ER/PgR [+] BC treated with neoadjuvant letrozole. Response to treatment was measured by radiology at 4th month by World Health Organization (WHO) criteria. Three polymorphisms of CYP19A1, one in exon 7 (rs700519) and two in the 3'-UTR region (rs10046 and rs4646) were evaluated on DNA obtained from peripheral blood.


Thirty-five women (36.8%) achieved a radiological response to letrozole. The histopathological and immunohistochemical parameters, including hormonal receptor status, were not associated with the response to letrozole. Only the genetic variants (AC/AA) of the rs4646 polymorphism were associated with poor response to letrozole (p = 0.03). Eighteen patients (18.9%) reported a progression of the disease. Those patients carrying the genetic variants (AC/AA) of rs4646 presented a lower progression-free survival than the patients homozygous for the reference variant (p = 0.0686). This effect was especially significant in the group of elderly patients not operated after letrozole induction (p = 0.009).


Our study reveals that the rs4646 polymorphism identifies a subgroup of stage II-III ER/PgR [+] BC patients with poor response to neoadjuvant letrozole and poor prognosis. Testing for the rs4646 polymorphism could be a useful tool in order to orientate the treatment in elderly BC patients.

Peer Review reports


After cessation of ovarian activity at postmenopause, the aromatase [cytochrome P450 19 (CYP19); OMIM:107910] from bones, adipose tissue and muscle, becomes the key enzyme in estradiol and estrone biosynthesis through the aromatization of testosterone and androstenedione respectively [1]. Aromatase is encoded by CYP19A1 which maps at chromosome 15q21.1 and has a complex structure with a region that contains 10 tissue-specific non-coding upstream exons with separate promoters that regulate transcription in different cells and tissues [2]. Elevated levels of aromatase expression have been observed in breast tumors relative to normal breast tissue [3], and alterations in aromatase expression have been implicated in the pathogenesis of estrogen-dependent diseases, including breast cancer (BC) [4, 5].

The major biological characteristic of BC is that two thirds of cases express estrogen (ER) and/or progesterone receptors (PgR) [2, 6, 7], endocrine manipulation being an effective treatment for these patients. Tamoxifen, a selective ER modulator, has for decades been the gold standard therapy in ER/PgR [+] BC patients [8]. However, in postmenopausal patients, aromatase inhibitors like letrozole have consistently demonstrated to be more effective than tamoxifen [9]. Letrozole is highly specific for the aromatase enzyme and inhibits whole-body aromatization of almost 99%, providing a highly and essentially complete withdrawal of estrogen in postmenopausal women [10].

Neoadjuvant hormonal therapy is a safety option for postmenopausal women with large operable or locally advanced hormone sensitive BC, increasing the possibilities of breast-conserving surgery (BCS) [11]. In this context, letrozole has shown to be more effective than tamoxifen, obtaining higher response and BCS rates following 4 months of preoperative therapy [12, 13]. However, less than 60% of patients will respond to this optimal endocrine therapy [14]. To date, only the presence and intensity of hormonal receptors (ER/PgR) are useful tools as predictive markers in clinical practice. Identification of more accurate markers for more efficient patient selection to exclude non-responsive cases remains crucial.

Genetic polymorphisms in the CYP19A1 gene have been associated with altered sex hormone levels in serum and urine [1517], providing an explanation for an elevated risk for BC in relation to estrogen exposure. For example, the tetranucleotide repeat polymorphism in intron 4 (TTTA)n has been associated with BC risk in initial studies [1822]. The polymorphism rs10046 in the 3' untranslated region (3'-UTR) also showed inconsistent associations with BC risk [2326] and has also been related with tumor stage [25], circulating sex hormone levels [15], HER2 status and disease free and overall survival [27, 28]. rs4646 has been reported to be associated with circulating steroid hormone levels [17] and, as in rs10046, with the HER2 status of the tumor [27]. Hence, it is biologically plausible that the polymorphisms in the CYP19A1 gene may be associated with the response to aromatase inhibitors. Indeed, recently the rs4646 CYP19A1 polymorphism has been associated with letrozole efficacy in advance disease [29]; however, no evidence between CYP19A1 polymorphisms and therapeutic efficacy of aromatase inhibitors in early stage BC has yet been established.

In the present study, we perform a genetic analysis of three CYP19A1 polymorphisms in a series of postmenopausal endocrine-sensitive BC patients treated with neoadjuvant letrozole and describe their association both with radiological response at 4 months (RR4M) of treatment and with progression free survival (PFS).


Patient selection and DNA samples

We identified a total of 153 postmenopausal patients treated in our institution with letrozole as primary systemic therapy for stage II - III BC between April 1999 and March 2008. Only those cases which met the following criteria were eligible for the study: 1) histological diagnosis of invasive BC by tru-cut or biopsy; 2) positive hormonal receptors status measured by immunohistochemistry (> 10% ER and/or PgR); 3) preserved tumor block; 4) bidimensionally measurable disease by mammogram and/or breast ultrasound; 5) adequate radiological follow-up at 4 months (± 14 days); 5) treatment with letrozole (Femara: Novartis Pharma AG, Basel, Switzerland) as neoadjuvant therapy by 2.5 mg/day for a minimum of 4 months in the absence of progression or side effects; and 6) patients able to provide a blood sample.

Finally, 95 patients fulfilled all criteria, were contacted by their physician (AGZ, ALC, ARS, JG, MAC, VG) and agreed to participate on the study. All patients gave their written consent for the use of their DNA for this specific research proposal. In addition, this study was reviewed and approved by the FIVO Science and Ethics Committees.

Peripheral blood samples (7 ml) were obtained from these patients and processed for DNA extraction at the Laboratory of Molecular Biology.

Response to treatment at four months was done using the WHO criteria: complete response (CR) required the complete disappearance of all disease, partial response (PR) was defined as a reduction of ≥ 50% in tumor volume and stable disease (SD) was any reduction < 50% in tumor volume. The response evaluation was done by an independent radiologist (JCD).

Pre-treatment study included patient characteristics, menopausal status, disease history, histological grade and diagnosis, evaluation of hormone receptor status, HER2 and Ki67 status by immunohistochemical analysis on core-cut biopsies, ECOG performance status and tumor assessment. The clinical data were prospectively reviewed and stored within a specific database.

Histological and immunohistochemical analysis

Core-cut biopsies were taken with a 14-guage needle before treatment and fixed in neutral buffered formalin for histological studies. The breast tumors were graded according to the modified Bloom and Richardson score on H&E-stained slides [30].

Four-micrometer sections from embedded blocks were cut on poly-L-lysine-coated slides, and dewaxed, endogenous peroxidase was inhibited with 3% hydrogen peroxide for 30 min. Immunoreactivity was enhanced with antigen retrieval treatment by heating the slides in a microwave oven for 10 min (700 W) in 10 mM sodium citrate buffer pH 7, followed by cooling for 20 min at room temperature. Sections were blocked with 20% horse serum in phosphate-buffer saline (PBS) and incubated with primary antibody for 1 h at room temperature. The incubation time for the secondary antibody and avidin-biotin complexes was 30 min at room temperature. Sections were extensively washed and the immunoreactions developed using DAB (0.05% 3'3' diamino-benzidine in 0.1% hydrogen peroxide). Negative controls included substitution of the primary antibody by mouse ascites or PBS. Slides were counterstained in Mayer hematoxylin, dehydrated, and mounted.

ER and PgR expression in primary tumor was immunohistochemically evaluated using the anti-ER (clone 6F11) and anti-PgR (clone 1A6) antibodies (Novocastra Laboratories, Newcastle upon Tyne, England) at a dilution of 1:40 and 1:30 respectively. Ki67 expression was also determined using the anti-Ki67 (clone MIB1) antibody (Dako Corp., Carpinteria, CA) at a 1:50 dilution.

HER2 protein expression was immunohistochemically evaluated using the HercepTest Kit (Dako Corp., Carpinteria, CA) according to the manufacturer's protocol. The scoring system was as follows: 0, tumors with no or weak staining in less than 10% of the cells; +, tumors with a faint or barely perceptible membrane staining in more than 10% of cells or with noncircumferential staining; ++, moderate circumferential membrane staining; +++, strong circumferential membrane staining. Scores 0 and + were considered as negative; scores ++ and +++ were considered as positive for HER2 overexpression [31].

DNA isolation

Genomic DNA used for genotyping studies was isolated from peripheral blood mononucleated cells, previously separated by centrifugation with Histopaque®-1077 (SIGMA-ALDRICH, LTD Irvine, UK), using the UltraClean DNA Blood Isolation Kit (MO BIO, Carlsbad, CA, USA). DNA integrity was evaluated by the A260/A280 absorbance ratio with a Nanodrop-1000 (Thermo Scientific) spectrophotometer.

CYP19A1 Genotyping

The single nucleotide polymorphisms (SNPs) rs10046 (A/G) and rs4646 (C/A) located at the 3'-untranslated region (3'-UTR) of CYP19A1, previously associated with letrozole efficacy in advanced BC [29], and the rs700519 (C/T) in the first base of codon 264 (Arg264Cys), related with low levels of aromatase activity [32], were evaluated. The rs700519 was analyzed by PCR amplification and direct sequencing. Briefly, amplification was carried out in a final volume of 25 μl containing 100 ng genomic DNA, 0.3 μM of forward (5'-CATGAAGTGTAGGGTC-TATGTAAT-3') and reverse (5'-GATCTTTACACACCTCTACACAGT-3') primers, 0.2 mM dNTPs (deoxyribonucleoside triphosphates), 2 mM MgCl2, 1× _Buffer II (Applied Biosystems, New Jersey, USA) and 1 unit AmpliTaq Gold (Applied Biosystems, New Jersey, USA). PCR conditions were 94°C for 10 minutes; 35 cycles with denaturation at 94°C for 30 seconds, annealing at 52°C for 30 seconds, and elongation at 72°C for 30 seconds; and a final extension step at 72°C for 10 minutes. PCR products were visualized on 1.5% agarose gels containing ethidium bromide and after UV radiation. Purified PCR products were then sequenced using BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Warrington WA1 4SR, UK) and electrophoresed on a ABI 3130xl Genetic Analyzer (Applied Biosystems). Sequencing analysis was performed with the Sequencing Analysis 5.2 Software (Applied Biosystems).

Genotyping of the rs10046 and rs4646 polymorphisms were carried out by allelic discrimination using specific TaqMan SNP Genotyping Assays and following the manufacturer's instructions. In brief, 20 μl PCR reactions were carried out containing 10.0 μl TaqMan Universal PCR Master Mix (Applied Biosystems, New Jersey, USA), 20.0 ng DNA template and 1.0 μL TaqMan SNP Genotyping Assays (Applied Biosystems, Foster City CA 94404). The references of selected assays were: C_8234731_1 (rs10046) and C_8234730_1 (rs4646). All assays were performed in 96-well plates including non-template controls. PCR reactions were read on a 7500 Fast Real Time PCR System in end-point mode using the Allelic Discrimination Sequence Detector software (Applied Biosystems).

Statistical analysis

Statistical analysis was carried out using the SPSS statistical software package (version 15.0.1, SPSS Inc., Chicago, Illinois, USA). A total of 95 patients were included in the analysis. Data on radiological response at the fourth month of induction therapy, progression and survival status, immunohistochemical markers, histological grade, stage and genotype were available for these patients.

For the statistical analysis we used binary variables reflecting the positivity status of the measures (yes or no) as well as the genotypic status of each polymorphism (homozygous reference or heterozygous/homozygous less frequent). Association with histopathological parameters, all categorical, was also assessed using a chi-square test to determine homogeneity or linear trend for ordinal variables. The significance level was set at 5%.

Deviations of genotype frequencies from those expected under Hardy-Weinberg equilibrium (HWE) were assessed by χ2 (p > 0.05) tests [33].

To study the impact of the biological factors on progression-free survival (PFS), the Kaplan-Meier proportional risk test (Log Rank) was used [34, 35]. PFS was defined as the time interval from the start of treatment to local or distant disease progression. Those patients that reported disease progression during the follow-up were considered as censored, being the data of confirmation of this progression considered in the calculation of PFS. Univariate predictors of PFS were entered into a Cox proportional hazards model using stepwise selection to identify the independent predictors of poor outcome, a confidence interval (CI) of 95% was also considered [36].


Patient and tumor characteristics

Ninety-five stage II-III ER/PgR-positive BC in postmenopausal women treated with letrozole as neoadjuvant therapy were analyzed. The median age of the series was 78.3 years (range: 61.7-89.8). The histopathological and clinical characteristics of the series are summarized in Table 1. As shown, 67 patients presented infiltrating ductal carcinoma (IDC), 16 infiltrating lobular carcinoma (ILC) and the remaining 12 cases corresponded to 9 mucinous carcinomas, 2 tubular carcinomas and one medullar carcinoma. The histological grade was available in only 52 of the IDC due to sample limitations. Most of the cases were stage IIa BC (68.4%), the majority of cases having a T2 tumor (80%). All tumors expressed HR (ER or PgR or both). Considering the whole HR expression for each patient we distinguished two groups of cases: those having low level HR (expressing HR in ≤ 40% of tumor cells) and those with a high HR expression (> 40% of cells). Ten cases were classified in the first group (10.5%). HER2 and Ki67 expression was also evaluated, 32 cases were HER2 ++/+++ (33.7%) and 32 out of 73 cases analyzed expressed Ki67 (44.4%).

Table 1 Clinical and pathological characteristics of the patients.

The median treatment period with letrozole as neoadjuvant therapy was 10.23 months (range: 3.53-85). Response to treatment was evaluated by radiological means (mammogram or ultrasound) at four months of letrozole administration using WHO criteria [37]. A total of 34 partial responses (PR) and one complete response (CR) where achieved, and were classified as responders for the analysis. Hence, the overall radiological response rate at 4 months (RR4M) was 36.8%. On the contrary, 55 women remained with stable disease (SD) and 5 developed progression disease (PD), these all being considered as non-responders.

CYP19A1 polymorphism analysis

A total of three single nucleotide polymorphisms (SNPs) were evaluated: one located in exon 7 (rs700519) and two in the 3'-UTR region of CYP19A1 (rs10046 and rs4646). Polymorphism references, genotypes, gene and chromosome locations, and frequency distribution of observed genotypes in the series are shown in Table 2. In addition, values for Hardy-Weinberg (HW) [33] equilibrium were estimated for each polymorphism and are also listed in Table 2. Both rs10046 and rs4646 were in HW equilibrium; however, for rs700519 this calculation was not possible as no genetic variant (CT, TT) was observed in any of the analyzed cases.

Table 2 CYP19A1 polymorphisms and genotype frequencies

Analysis between RR4M and clinical, pathological and genotypic parameters

No statistically significant association was found between RR4M and the clinical and pathological parameters (Table 3). However, an overrepresentation of the genetic variants (AC/AA) of the rs4646 polymorphism was observed in the non-responder group (48.3%) compared with the responders (25.7%) (p = 0.03) (Table 3).

Table 3 Association of clinical, pathological and genotypic parameters with radiological response at 4th month of letrozole administration.

Clinical, pathological, genotypic parameters and progression free survival (PFS)

The median follow-up of the series was of 40.6 months (range: 4.3-113.8 months). A total of 18 patients (18.9%) reported disease progression during this period (Table 1). Among the analyzed parameters, IIIa-IIIb stages and the non response to letrozole induction therapy were significantly associated with shorter PFS (Table 4; Additional file 1: Supplemental Table S1; Fig. 1). Interestingly, those patients carrying the genetic variants (AC/AA) of rs4646 presented a worse PFS than the patients homozygous for the reference variant (Table 4; Additional file 1: Supplemental Table S1; Fig. 2). Multivariate analysis showed that lymph node involvement [Hazard Ratio (HR) = 4.3 (1.7-11); p = 0.002] and RR4M [HR = 11 (1.5-83); p = 0.012] constituted independent prognostic factors for PFS.

Figure 1
figure 1

Kaplan Meier plot showing the progression-free survival with regard the radiological response to letrozole at 4 th month (RR4M) of the induction therapy.

Figure 2
figure 2

Kaplan Meier plot showing the progression-free survival with regard the genetic variants of rs4646 in the global series.

Table 4 Association of clinical, pathological and genotypic parameters with progression free survival (PFS).

Clinical, pathological, genotypic parameters and PFS according to the surgery status

Sixty-five out of the 95 patients (68%) underwent surgery after letrozole induction therapy and no association was observed between the progression status and surgery (Table 1). In fact, 20% of non-surgery cases progressed, while 18.5% of the surgery group experienced progression events (p = 0.861). In addition, no significant differences in terms of PFS were observed between the non-surgery (74.5%) and surgery (60.4%) groups (p = 0.600). However, when we compared clinical, pathological and genetic parameters between the two groups of patients according to surgery status, we observed that 77% of patients of the non-surgery group were over the median age, compared with 37% of the patients from the surgery group (p = 0.00032). Moreover, rs4646 variant (AC/AA) was overrepresented in the non-surgery group (57%) in contrast to the surgery group (32%) (p = 0.024).

Univariate analysis performed for PFS for each group of patients revealed that in the non-surgery group, RR4M and r4646 where associated with a worse prognosis (Table 5; Additional file 2: Supplemental Table S2). Indeed, the 6 patients who progressed were genetically characterized by the variant rs4646 (AC/AA) (Fig. 3A). However, in the multivariate analysis only the polymorphism rs4646 emerged as an independent prognostic factor although it did not reach statistical significance (p = 0.191).

Figure 3
figure 3

Kaplan Meier plot showing the progression-free survival with regard the genetic variants of rs4646 in the no surgery (A) and surgery group of patients (B).

Table 5 Association of clinical and genotypic parameters with progression free survival (PFS) according to surgery status.

In the group of patients operated after letrozole induction, only Ki67 expression and RR4M were associated with worse PFS (Table 5; Additional file 2: Supplemental Table S2), although none of these variables were independent indicators of poor outcome. Contrary to that observed in the non-surgery group, no association of the rs4646 variants with PFS was observed (Fig. 3B).


We describe a relationship between the genetic variants (AC/AA) of the rs4646 polymorphism of CYP19A1 and a poor response to neoadjuvant treatment with letrozole in postmenopausal women with ER/PgR [+] BC. The radiological response was measured at the 4th month of letrozole induction, the observed response rate being 36.8%, similar to the radiological response reported in the P024 trial [12].

At genetic level, mutations in the aromatase gene may lead to a functionally less sensitive aromatase phenotype. Although mutations in the CYP19A1 gene can be generated in vitro, to date no somatic mutation has been found in clinical samples [38]. However, several genetic polymorphisms of CYP19A1 have been reported so far, although the possible functional significance of most of these remains undefined. In the case of BC, population-based studies of common CYP19A1 polymorphisms have generated inconsistent results with regard to their possible association with sex hormone levels, cancer risk, HER2 status or survival [15, 19, 21, 27, 28]. No association between the HER2 or hormonal status of either rs10046 or rs4646 was observed in our series.

Whether aromatase inhibitor agents interact with the different CY19A1 genotypes has not so far been clearly established. Only three previous studies have attempted to demonstrate the effect of certain polymorphisms in CYP19A1 on the efficacy of aromatase inhibitors. In the first, Fasching et al., in a series of patients treated with hormonal therapy for more than 2 years, did not observe significant differences between the rs700519, rs10046 and rs4646 polymorphisms with therapy. Unfortunately, the type of antihormonal therapy was not recorded for this analysis [27]. In the second study, a population-based and in vitro study, Ma et al. revealed reductions in the functional activity of aromatase for four phenotypes resulting from non-synonymous changes [rs700519 (Arg264Cys), rs28757184 (Thr201Met), rs2236722 (Trp39Arg), and rs56658716 (Met364Thr)] [32]. These authors found that levels or aromatase enzyme activity decreased dramatically for the Thr364, and also observed a slight decrease in Cys264 allozyme activity. They also demonstrated that the mechanism by which non-synonymous SNPs affect the enzymatic activity is a consequence of an alteration in the enzyme protein level [32]. In the same study, these authors reported that these variants had no significant differences in their affinity to the aromatase inhibitors letrozole and exemestane [32]. Interestingly, the polymorphism rs700519 has also been reported as a prognostic factor in a Chinese study, mainly in a subgroup of premenopausal women [28], but not in the Caucasian population [27]. In addition, it is important to note the differences between these two populations regarding the minor allele frequency of rs700519, being only 3.2% for the Caucasian population (1,257 patients) [27], in comparison with 15.1% in the Asian Cohort (1,136 patients) [28]. In our series, the genetic variants of this polymorphism were underrepresented and not found in any of the cases.

The third study evaluated the efficacy of treatment with letrozole in advanced hormone receptor-positive BC patients with respect to two polymorphisms located at the 3'-UTR (rs10046 and rs4646) and one in intron 2 [rs727479 (G/T)] of the CYP19A1 [29]. The authors found that genetic variants of rs4646 were associated with a greater efficacy of letrozole in terms of time to progression. Hence, patients with the variant genotype (AC/AA) had a three times greater time to progression than the patients with the reference genotype (CC). Furthermore, these authors reported that the frequency of the variant alleles for rs4646 was significantly higher in the responder (61%) than in the non-responder group of patients (40%) [29]. These observations clearly differ from those reported in our series, the genotypic variants of rs4646 being more frequently represented in the non-responder group to letrozole after 4 months of induction therapy (48% vs. 26%). In this regard, there are two main issues that could explain the differences between these two studies. First, and independently of the disease stage, are the criteria for patient selection. Whereas in our series the patients were treated with letrozole as a first line, in the study of Colomer et al. the patients had progressed from a previous treatment with tamoxifen [29]. It is well established that tamoxifen metabolism is influenced by the number of mutant alleles of the gene encoding cytochrome P450 2D6 (CYP2D6), because this enzyme affects the levels of endoxifen, the active tamoxifen metabolite. Thus, for patients who are wild type for CYP2D6, the 5-year disease-free survival outcomes are similar to or perhaps even superior with tamoxifen than with aromatase inhibitors [39]. Therefore, previous treatment with tamoxifen could have genetically selected a population of patients more sensitive to letrozole. Secondly, in our study we used constitutive DNA obtained from peripheral blood for the genetic analysis, whereas Colomer et al. performed the genetic study on DNA obtained from the fixed and paraffin-embedded tumors [29]. The polymorphism studies performed on tumor samples are subject to genetic alterations that can affect the chromosomal region of the gene under study. In this way, CYP19A1 is located in the 21.2 region of the long arm of chromosome 15 (15q21.2) [40], and this region has been reported to be a frequent target of allelic imbalance in advanced breast carcinomas [41], and which could affect the frequency distribution of the allelic variants.

From our series, 65 women underwent surgery after letrozole induction, being age the only difference observed between the operated and non-operated patients. In fact 77% of patients within the non-surgery group were over 78.3 years. We have demonstrated that genetic variants of rs4646 have prognostic value, especially in this group of patients. Indeed, the 6 women who progressed within the group of non-surgery patients were genetically characterized by the variant rs4646 (AC/AA). In this elderly group of patients the decision for undergoing surgery is not always easy and very often both the patient and the clinician opt for a local control of the tumor maintaining the treatment with letrozole. Hence, in this group of patients, in which the tumor is not removed, is where the rs4646 polymorphism identifies the women who progressed, indicating, and despite the small number the cases included in our series, that women with the rs4646 (AC/AA) genotype could benefit from another therapeutic approach.

The regulatory interactions between the ER, growth factor receptors and other kinase signalling pathways could also determine the response to endocrine therapy. In the adjuvant setting, several studies suggest that patients overexpressing HER2 may derive relatively less benefit from endocrine therapy [38]. In the P024 trial, ER [+] tumors that were also EGFR and/or HER2 [+] responded significantly better to letrozole than to tamoxifen (88% vs 21% respectively, P = 0.0004) [13]. The IMPACT study, that confronted tamoxifen with anastrozole, or the preoperative combination of both, observed a similar effect in favour of anastrozole for HER2 [+] tumors (Odds Ratio 58% vs. 22%), although not reaching statistical significance [42]. In our series, as in both randomized trials, no differences in activity were observed relating to HER2 expression among the patients receiving letrozole.

Biological studies of tissues obtained during neoadjuvant therapy have thrown a different perspective on this issue. Short-term estrogen deprivation leads to profound changes in transcriptional profiles, and these changes can be used as predictive tools. Tumor expression of the proliferation antigen Ki67 after 2 weeks of endocrine treatment (tamoxifen or anastrozole) predicts for clinical response and recurrence-free survival [42]. However, among cases on anastrozole, only 7% did not present a reduction in the Ki67 expression at 2 weeks. In many cases, this reduction in the tumor proliferation rate is only modest and could not be sufficient to determine a resistance to the treatment. As expected, early tumor changes are not limited to proliferation markers, but to other genes directly regulated by ER, including the aromatase itself [43]. Other, currently ongoing, approaches have used gene expression profiling techniques able to predict response to endocrine therapy, but the results obtained so far are inconsistent [43]. In our series, expression of proliferating marker Ki67 constituted a factor of poor prognosis in the group of patients undergoing surgery following letrozole induction therapy, suggesting that this marker could identify tumors with a more aggressive behavior.


We have found that genetic variants of the polymorphism rs4646 in the 3'-UTR of CYP19A1 are associated to poor response to letrozole after 4 months induction therapy and to poor outcome of elderly patients without surgery after letrozole administration in postmenopausal women with ER/PgR [+] BC. Although our study would suggest that analysis of rs4646 could improve the clinical management of these patients by facilitating a more individualized therapy, further studies based on larger series are necessary to confirm these findings.



3'-untranslated region

CYP19A1 :

aromatase gene


breast cancer


breast conserving surgery


complete response




estrogen recetor


Hardy-Weinberg equilibrium


infiltrating ductal carcinoma


infiltrating lobular carcinoma


partial response


phosphate-buffer saline


progesterone receptor


progression disease


progression free survival


radiological response at 4 months


single nucleotide polymorphisms


stable disease


World Health Organization.


  1. Simpson ER, Mahendroo MS, Means GD, Kilgore MW, Hinshelwood MM, Graham-Lorence S, Amarneh B, Ito Y, Fisher CR, Michael MD, et al: Aromatase cytochrome P450, the enzyme responsible for estrogen biosynthesis. Endocr Rev. 1994, 15 (3): 342-355.

    CAS  PubMed  Google Scholar 

  2. Bulun SE, Lin Z, Imir G, Amin S, Demura M, Yilmaz B, Martin R, Utsunomiya H, Thung S, Gurates B, et al: Regulation of aromatase expression in estrogen-responsive breast and uterine disease: from bench to treatment. Pharmacol Rev. 2005, 57 (3): 359-383. 10.1124/pr.57.3.6.

    Article  CAS  PubMed  Google Scholar 

  3. Irahara N, Miyoshi Y, Taguchi T, Tamaki Y, Noguchi S: Quantitative analysis of aromatase mRNA expression derived from various promoters (I.4, I.3, PII and I.7) and its association with expression of TNF-alpha, IL-6 and COX-2 mRNAs in human breast cancer. Int J Cancer. 2006, 118 (8): 1915-1921. 10.1002/ijc.21562.

    Article  CAS  PubMed  Google Scholar 

  4. O'Neill JS, Elton RA, Miller WR: Aromatase activity in adipose tissue from breast quadrants: a link with tumour site. Br Med J (Clin Res Ed). 1988, 296 (6624): 741-743. 10.1136/bmj.296.6624.741.

    Article  Google Scholar 

  5. Agarwal VR, Bulun SE, Leitch M, Rohrich R, Simpson ER: Use of alternative promoters to express the aromatase cytochrome P450 (CYP19) gene in breast adipose tissues of cancer-free and breast cancer patients. J Clin Endocrinol Metab. 1996, 81 (11): 3843-3849. 10.1210/jc.81.11.3843.

    CAS  PubMed  Google Scholar 

  6. Toniolo PG, Levitz M, Zeleniuch-Jacquotte A, Banerjee S, Koenig KL, Shore RE, Strax P, Pasternack BS: A prospective study of endogenous estrogens and breast cancer in postmenopausal women. J Natl Cancer Inst. 1995, 87 (3): 190-197. 10.1093/jnci/87.3.190.

    Article  CAS  PubMed  Google Scholar 

  7. Paridaens R, Thomas J, Wildiers J, Vermeiren P, Lobelle JP, di Salle E, Ornati G, Zurlo MG, Polli A, Lanzalone S, et al: Safety, activity and estrogen inhibition by exemestane in postmenopausal women with advanced breast cancer: a phase I study. Anticancer Drugs. 1998, 9 (8): 675-683. 10.1097/00001813-199809000-00002.

    Article  CAS  PubMed  Google Scholar 

  8. Early Breast Cancer Trialists' Collaborative Group (EBCTCG).: Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: an overview of the randomised trials. Lancet. 2005, 365 (9472): 1687-1717. 10.1016/S0140-6736(05)66544-0.

    Article  Google Scholar 

  9. Smith IE, Dowsett M: Aromatase inhibitors in breast cancer. N Engl J Med. 2003, 348 (24): 2431-2442. 10.1056/NEJMra023246.

    Article  CAS  PubMed  Google Scholar 

  10. Geisler J, Haynes B, Anker G, Dowsett M, Lonning PE: Influence of letrozole and anastrozole on total body aromatization and plasma estrogen levels in postmenopausal breast cancer patients evaluated in a randomized, cross-over study. J Clin Oncol. 2002, 20 (3): 751-757. 10.1200/JCO.20.3.751.

    Article  CAS  PubMed  Google Scholar 

  11. Semiglazov VF, Semiglazov VV, Dashyan GA, Ziltsova EK, Ivanov VG, Bozhok AA, Melnikova OA, Paltuev RM, Kletzel A, Berstein LM: Phase 2 randomized trial of primary endocrine therapy versus chemotherapy in postmenopausal patients with estrogen receptor-positive breast cancer. Cancer. 2007, 110 (2): 244-254. 10.1002/cncr.22789.

    Article  CAS  PubMed  Google Scholar 

  12. Eiermann W, Paepke S, Appfelstaedt J, Llombart-Cussac A, Eremin J, Vinholes J, Mauriac L, Ellis M, Lassus M, Chaudri-Ross HA, et al: Preoperative treatment of postmenopausal breast cancer patients with letrozole: A randomized double-blind multicenter study. Ann Oncol. 2001, 12 (11): 1527-1532. 10.1023/A:1013128213451.

    Article  CAS  PubMed  Google Scholar 

  13. Ellis MJ, Coop A, Singh B, Mauriac L, Llombert-Cussac A, Janicke F, Miller WR, Evans DB, Dugan M, Brady C, et al: Letrozole is more effective neoadjuvant endocrine therapy than tamoxifen for ErbB-1- and/or ErbB-2-positive, estrogen receptor-positive primary breast cancer: evidence from a phase III randomized trial. J Clin Oncol. 2001, 19 (18): 3808-3816.

    CAS  PubMed  Google Scholar 

  14. Smith I, Dowsett M, trialists oboI: Comparison of anastrozole vs tamoxifen alone and in combination as neoadjuvant treatment of estrogen receptor-positive (ER+) operable breast cancer in posmenopausal women: The IMPACT trial. Breast Cancer Res Treat. 2003, 82 (suppl 2): S9-S15.

    Article  Google Scholar 

  15. Dunning AM, Dowsett M, Healey CS, Tee L, Luben RN, Folkerd E, Novik KL, Kelemen L, Ogata S, Pharoah PD, et al: Polymorphisms associated with circulating sex hormone levels in postmenopausal women. J Natl Cancer Inst. 2004, 96 (12): 936-945.

    Article  CAS  PubMed  Google Scholar 

  16. Tworoger SS, Chubak J, Aiello EJ, Ulrich CM, Atkinson C, Potter JD, Yasui Y, Stapleton PL, Lampe JW, Farin FM, et al: Association of CYP17, CYP19, CYP1B1, and COMT polymorphisms with serum and urinary sex hormone concentrations in postmenopausal women. Cancer Epidemiol Biomarkers Prev. 2004, 13 (1): 94-101. 10.1158/1055-9965.EPI-03-0026.

    Article  CAS  PubMed  Google Scholar 

  17. Haiman CA, Dossus L, Setiawan VW, Stram DO, Dunning AM, Thomas G, Thun MJ, Albanes D, Altshuler D, Ardanaz E, et al: Genetic variation at the CYP19A1 locus predicts circulating estrogen levels but not breast cancer risk in postmenopausal women. Cancer Res. 2007, 67 (5): 1893-1897. 10.1158/0008-5472.CAN-06-4123.

    Article  CAS  PubMed  Google Scholar 

  18. Ribeiro FS, de Amorim LM, de Almeida Simao T, Mendonca GA, de Moura Gallo CV, Pinto LF: CYP19 (TTTA)n polymorphism and breast cancer risk in Brazilian women. Toxicol Lett. 2006, 164 (1): 90-95. 10.1016/j.toxlet.2005.11.014.

    Article  CAS  PubMed  Google Scholar 

  19. Kristensen VN, Andersen TI, Lindblom A, Erikstein B, Magnus P, Borresen-Dale AL: A rare CYP19 (aromatase) variant may increase the risk of breast cancer. Pharmacogenetics. 1998, 8 (1): 43-48. 10.1097/00008571-199802000-00006.

    Article  CAS  PubMed  Google Scholar 

  20. Baxter SW, Choong DY, Eccles DM, Campbell IG: Polymorphic variation in CYP19 and the risk of breast cancer. Carcinogenesis. 2001, 22 (2): 347-349. 10.1093/carcin/22.2.347.

    Article  CAS  PubMed  Google Scholar 

  21. Probst-Hensch NM, Ingles SA, Diep AT, Haile RW, Stanczyk FZ, Kolonel LN, Henderson BE: Aromatase and breast cancer susceptibility. Endocr Relat Cancer. 1999, 6 (2): 165-173. 10.1677/erc.0.0060165.

    Article  CAS  PubMed  Google Scholar 

  22. Thyagarajan B, Brott M, Mink P, Folsom AR, Anderson KE, Oetting WS, Gross M: CYP1B1 and CYP19 gene polymorphisms and breast cancer incidence: no association in the ARIC study. Cancer Lett. 2004, 207 (2): 183-189. 10.1016/j.canlet.2003.12.009.

    Article  CAS  PubMed  Google Scholar 

  23. Ralph DA, Zhao LP, Aston CE, Manjeshwar S, Pugh TW, DeFreese DC, Gramling BA, Shimasaki CD, Jupe ER: Age-specific association of steroid hormone pathway gene polymorphisms with breast cancer risk. Cancer. 2007, 109 (10): 1940-1948. 10.1002/cncr.22634.

    Article  CAS  PubMed  Google Scholar 

  24. Olson JE, Ingle JN, Ma CX, Pelleymounter LL, Schaid DJ, Pankratz VS, Vierkant RA, Fredericksen ZS, Wu Y, Couch FJ, et al: A comprehensive examination of CYP19 variation and risk of breast cancer using two haplotype-tagging approaches. Breast Cancer Res Treat. 2007, 102 (2): 237-247. 10.1007/s10549-006-9324-7.

    Article  PubMed  Google Scholar 

  25. Kristensen VN, Harada N, Yoshimura N, Haraldsen E, Lonning PE, Erikstein B, Karesen R, Kristensen T, Borresen-Dale AL: Genetic variants of CYP19 (aromatase) and breast cancer risk. Oncogene. 2000, 19 (10): 1329-1333. 10.1038/sj.onc.1203425.

    Article  CAS  PubMed  Google Scholar 

  26. Haiman CA, Hankinson SE, Spiegelman D, Brown M, Hunter DJ: No association between a single nucleotide polymorphism in CYP19 and breast cancer risk. Cancer Epidemiol Biomarkers Prev. 2002, 11 (2): 215-216.

    CAS  PubMed  Google Scholar 

  27. Fasching PA, Loehberg CR, Strissel PL, Lux MP, Bani MR, Schrauder M, Geiler S, Ringleff K, Oeser S, Weihbrecht S, et al: Single nucleotide polymorphisms of the aromatase gene (CYP19A1), HER2/neu status, and prognosis in breast cancer patients. Breast Cancer Res Treat. 2008, 112 (1): 89-98. 10.1007/s10549-007-9822-2.

    Article  CAS  PubMed  Google Scholar 

  28. Long JR, Kataoka N, Shu XO, Wen W, Gao YT, Cai Q, Zheng W: Genetic polymorphisms of the CYP19A1 gene and breast cancer survival. Cancer Epidemiol Biomarkers Prev. 2006, 15 (11): 2115-2122. 10.1158/1055-9965.EPI-06-0464.

    Article  CAS  PubMed  Google Scholar 

  29. Colomer R, Monzo M, Tusquets I, Rifa J, Baena JM, Barnadas A, Calvo L, Carabantes F, Crespo C, Munoz M, et al: A single-nucleotide polymorphism in the aromatase gene is associated with the efficacy of the aromatase inhibitor letrozole in advanced breast carcinoma. Clin Cancer Res. 2008, 14 (3): 811-816. 10.1158/1078-0432.CCR-07-1923.

    Article  CAS  PubMed  Google Scholar 

  30. Elston CW, Ellis IO: Pathological prognostic factors in breast cancer. I. The value of histological grade in breast cancer: experience from a large study with long-term follow-up. Histopathology. 2002, 41 (3A): 151-10.1046/j.1365-2559.2002.14691.x. Comment on C. W. Elston & I. O. Ellis. Histopathology 1991; 19; 403-410

    Article  CAS  PubMed  Google Scholar 

  31. Harris LN, Liotcheva V, Broadwater G, Ramirez MJ, Maimonis P, Anderson S, Everett T, Harpole D, Moore MB, Berry DA, et al: Comparison of methods of measuring HER-2 in metastatic breast cancer patients treated with high-dose chemotherapy. J Clin Oncol. 2001, 19 (6): 1698-1706.

    CAS  PubMed  Google Scholar 

  32. Ma CX, Adjei AA, Salavaggione OE, Coronel J, Pelleymounter L, Wang L, Eckloff BW, Schaid D, Wieben ED, Adjei AA, et al: Human aromatase: gene resequencing and functional genomics. Cancer Res. 2005, 65 (23): 11071-11082. 10.1158/0008-5472.CAN-05-1218.

    Article  CAS  PubMed  Google Scholar 

  33. Stern C: The Hardy-Weinberg Law. Science. 1943, 97 (2510): 137-138. 10.1126/science.97.2510.137.

    Article  CAS  PubMed  Google Scholar 

  34. Kaplan E, Meier P: Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958, 53: 457-481. 10.2307/2281868.

    Article  Google Scholar 

  35. Peto R, Peto P: Asymptomatically efficient rank invariant test procedures. J Roy Stat Soc. 1987, 135: 185-206.

    Article  Google Scholar 

  36. Cox D: Regression models and life tables. J R Stat Soc. 1972, 34: 187-220.

    Google Scholar 

  37. Miller AB, Hoogstraten B, Staquet M, Winkler A: Reporting results of cancer treatment. Cancer. 1981, 47 (1): 207-214. 10.1002/1097-0142(19810101)47:1<207::AID-CNCR2820470134>3.0.CO;2-6.

    Article  CAS  PubMed  Google Scholar 

  38. Dowsett M, Martin LA, Smith I, Johnston S: Mechanisms of resistance to aromatase inhibitors. J Steroid Biochem Mol Biol. 2005, 95 (1-5): 167-172. 10.1016/j.jsbmb.2005.04.022.

    Article  CAS  PubMed  Google Scholar 

  39. Punglia RS, Burstein HJ, Winer EP, Weeks JC: Pharmacogenomic variation of CYP2D6 and the choice of optimal adjuvant endocrine therapy for postmenopausal breast cancer: a modeling analysis. J Natl Cancer Inst. 2008, 100 (9): 642-648. 10.1093/jnci/djn100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Bulun SE, Takayama K, Suzuki T, Sasano H, Yilmaz B, Sebastian S: Organization of the human aromatase p450 (CYP19) gene. Seminars in reproductive medicine. 2004, 22 (1): 5-9. 10.1055/s-2004-823022.

    Article  CAS  PubMed  Google Scholar 

  41. Rhiem K, Klein A, Munch M, Kreutzfeld R, Ramser J, Wardelmann E, Schackert G, Von Deimling A, Wiestler OD, Schmutzler RK: Chromosomal region 15q21.1 is a frequent target of allelic imbalance in advanced breast carcinomas. Int J Cancer. 2003, 106 (1): 74-77. 10.1002/ijc.11169.

    Article  CAS  PubMed  Google Scholar 

  42. Dowsett M, Ebbs SR, Dixon JM, Skene A, Griffith C, Boeddinghaus I, Salter J, Detre S, Hills M, Ashley S, et al: Biomarker changes during neoadjuvant anastrozole, tamoxifen, or the combination: influence of hormonal status and HER-2 in breast cancer--a study from the IMPACT trialists. J Clin Oncol. 2005, 23 (11): 2477-2492. 10.1200/JCO.2005.07.559.

    Article  CAS  PubMed  Google Scholar 

  43. Mackay A, Urruticoechea A, Dixon JM, Dexter T, Fenwick K, Ashworth A, Drury S, Larionov A, Young O, White S, et al: Molecular response to aromatase inhibitor treatment in primary breast cancer. Breast Cancer Res. 2007, 9 (3): R37-10.1186/bcr1732.

    Article  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references


We thank Maria Garcia Flores, Tania Mazcuñan and Laura Martinez for their technical assistance.

This study has been financed by grants GV0A6/256 from the Conselleria d'Empresa Universitat i Ciència of the Generalitat Valenciana (JALG) and SOLTI 06/01 (ALLC).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Jose A Lopez-Guerrero.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

ZGC carried out the molecular genetic studies, participated in the genetic analysis and drafted the manuscript. AGZ reviewed the clinical records of the patients included in the study and drafted the manuscript. AFS performed the genetic analysis. ALlC, ARS, JG, MAC and VG treated most of the patients included in the study and updated the clinical follow-up. AC, SA and ALlB carried out the histopathological and immunohistochemical analysis. JCD performed the radiological evaluation of the study. JC and CVA performed the trucut biopsies and informed to the patients about the study. JALG conceived of the study, performed the statistical analysis and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Zaida Garcia-Casado, Angel Guerrero-Zotano contributed equally to this work.

Electronic supplementary material


Additional file 1: Table S1. Association of clinical, pathological and genotypic parameters with progression free survival (PFS). This table shows the associations of the clinicopathological and genetic parameters with PFS that have not been statistically significant. (DOC 38 KB)


Additional file 2: Table S2. Association of clinical, histopathological and genotypic parameters with progression free survival (PFS) according to surgery status. This table shows the associations of the clinicopathological and genetic parameters with PFS that have not been statistically significant according to surgery status. (DOC 54 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Garcia-Casado, Z., Guerrero-Zotano, A., Llombart-Cussac, A. et al. A polymorphism at the 3'-UTR region of the aromatase gene defines a subgroup of postmenopausal breast cancer patients with poor response to neoadjuvant letrozole. BMC Cancer 10, 36 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: