ORIGINAL RESEARCH article
The Prognostic Impact of Intratumoral Aryl Hydrocarbon Receptor in Primary Breast Cancer Depends on the Type of Endocrine Therapy: A Population-Based Cohort Study
- 1Division of Oncology, Department of Clinical Sciences, Lund, Lund University and Skåne University Hospital, Lund, Sweden
- 2Department of Molecular & Cellular Biology and Dan L. Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, TX, United States
- 3Division of Oncology and Therapeutic Pathology, Department of Clinical Sciences, Lund, Lund University and Skåne University Hospital, Lund, Sweden
- 4Division of Surgery, Department of Clinical Sciences, Lund, Lund University, Lund, Sweden
- 5Department of Surgery, Kristianstad Hospital, Kristianstad, Sweden
- 6Department of Oncology, Aarhus University and Aarhus University Hospital, Aarhus, Denmark
The aryl hydrocarbon receptor (AhR) is a master regulator of multiple pathways involved in breast cancer, and influences the estrogen receptor alpha (ER) and aromatase/CYP19A1. The purpose of this study was to elucidate the interplay between intratumoral levels of AhR and aromatase, patient characteristics (including AhR and CYP19A1 genotypes), clinicopathological features, and prognosis in breast cancer patients receiving adjuvant treatments. A prospective cohort of 1116 patients with primary breast cancer in Sweden, included 2002–2012, was followed until June 30th 2019 (median 8.7 years). Tumor‐specific AhR (n=920) and aromatase levels (n=816) were evaluated on tissue microarrays using immunohistochemistry. Associations between cytoplasmatic (AhRcyt) and nuclear (AhRnuc) AhR levels, intratumoral aromatase, clinicopathological features, and prognosis in different treatment groups were analyzed. Low AhRcyt levels (n=183) and positive intratumoral aromatase (n=69) were associated with estrogen receptor (ER)– status and more aggressive tumors. Genotypes were not associated with their respective protein levels. The functional AhRArg554Lys GG genotype was associated with recurrence-free survival in switch-therapy (sequential tamoxifen/aromatase inhibitors (AI) or AI/tamoxifen) treated patients (HRadj 0.42; 95% CI 0.22–0.83). High AhRcyt levels were associated with longer recurrence-free survival during the first 10 years of follow-up among tamoxifen-only treated patients (HRadj 0.40; 95% CI 0.23–0.71) compared to low AhRcyt levels, whereas an almost inverse association was seen in patients with switch-therapy (Pinteraction=0.023). Intratumoral aromatase had little prognostic impact. These findings warrant confirmation in an independent cohort, preferably in a randomized clinical trial comparing different endocrine regimens. They might also guide the selection of breast cancer patients for clinical trials with selective AhR modulators.
Breast cancer remains an important cause of disease burden and death in women despite novel therapeutic options (1). Along with new treatments, important predictive markers have emerged, guiding the selection of targeted breast cancer therapies (2). However, novel prognostic and predictive tumor markers that can be targeted with new or repurposed treatment choices are urgently needed to minimize both over- and undertreatment. Most breast cancer patients receive adjuvant endocrine therapy, and some up to 10 years following surgery (2). Two potential intratumoral markers that merit further investigation in the adjuvant setting are the master regulator aryl hydrocarbon receptor (AhR) and aromatase, the key enzyme in androgen to estrogen conversion.
We have previously shown that the AhRArg554Lys polymorphism modifies the relationship between aromatase inhibitor (AI) response and CYP1A2 (3), which is an important enzyme in estrogen metabolism. AhR is a ligand-activated transcription factor with a wide range of endogenous and exogenous ligands (e.g., toxins, such as dioxin, and several drugs including raloxifene, 4-hydroxy-tamoxifen, in addition to tobacco smoke, and cruciferous vegetables) (4–7). After binding to ligands, the cytoplasmic AhR (AhRcyt) translocates to the nucleus and dimerizes with its transcriptional partner AhR nuclear translocator (ARNT), whereby the complex recognizes the dioxin response elements in the promoter of downstream genes that include several cytochrome P450 (CYP) genes (8). These genes include CYP19A1, which encodes aromatase, in addition to CYP1A1, CYP1A2, and CYP1B1 (9, 10). Therefore, the subcellular AhR localization may be of importance.
The functional AhRArg554Lys polymorphism influences the mRNA expression level of AhR (11). We have previously reported an association between the same AhRArg554Lys polymorphism and ER status in breast cancer, and that this polymorphism impacted response to endocrine switch-therapy with sequential tamoxifen/AI or AI/tamoxifen (12). Some groups have reported that various CYP19A1 genotypes might impact response to AIs (13–16), but this association has been confirmed neither by our group nor others (3, 17).
AhR regulates multiple pathways that might influence all the major stages of carcinogenesis (18). High expression of AhR in breast cancer has been associated with signaling pathways related to metabolism and insulin-like growth factor (IGF) signaling (19). AhR acts as an immunomodulator (8) and is a suggested link between inflammation and breast cancer (20). Several studies have investigated the prognostic impact of AhR expression in breast cancer but with inconsistent results (5, 19, 21–23). A couple of these studies suggested a positive association between AhR levels and good prognosis (5, 22), whereas others indicated a negative prognostic impact (21), that possibly differed between subgroups of patients and according to the intracellular localization of AhR (23).
AhR can be effectively modulated by its ligands, e.g., selective AhR modulators (24), resulting in either agonistic or antagonistic effects on many of the hallmarks of cancer (6). A clinical phase I trial investigating an AhR inhibitor in patients with advanced solid tumors is ongoing (NCT04069026).
The cross-talk between the AhR and the estrogen receptor alpha (ER) plays a major role in signaling processes in female reproductive organs, and it has been shown in in vivo models that ligand-activated AhR confers anti-estrogenic effects partly due to lower ER levels in ductal epithelial cells (25). Furthermore, studies have shown AhR-mediated degradation of ER through activation of the proteasome pathway (26, 27). Selective estrogen receptor modulators (SERMs) can act as AhR agonists in some cases (28), such as the active metabolite of tamoxifen, 4-hydroxy-TAM, that modulates the transcriptional activity of AhR (29). Another SERM, raloxifene, induced apoptosis in ER– breast cancer cells (5), which implies that AhR also plays a role in the hormone-independent setting. A small study reported that AhR might induce intratumoral aromatase and thereby stimulate estrogen-dependent breast cancer progression (30). AhR is thus a potential target for new drugs in breast cancer (28, 31).
We hypothesized that high tumor levels of AhR would be associated with good prognosis in breast cancer, but that the prognostic impact might depend on the subcellular AhR localization, treatments, as well as body constitution. We hypothesized that intratumoral AhR and aromatase levels would be associated with each other but not with polymorphisms in their respective genes. The aim was to study associations between intratumoral levels of AhR and aromatase, patient characteristics, including AhR and CYP19A1 genotypes, and clinicopathological features and prognosis in different treatment groups of primary breast cancer patients.
Materials and Methods
Primary breast cancer patients in Lund, Sweden were invited to participate in an ongoing prospective study, the Breast Cancer (BC) blood study. Patients with any previous breast cancer or another cancer diagnosis within 10 years were not eligible for participation. Between October 2002 and June 2012, 1116 patients were enrolled. After excluding 51 patients who had received preoperative treatment, 39 patients with in situ carcinoma and eight patients with early recurrences within 0.3 years, 1018 patients with invasive breast cancer remained. A flowchart of the selection criteria is presented in Figure 1.
Figure 1 Inclusion flowchart showing the number of included and excluded patients and representative images of AhR and aromatase staining intensities (40×). The bar represents 20 µm.
Preoperatively, the patients completed a questionnaire as previously described (32). Blood samples were obtained, and body measurements were taken by a trained research nurse. Body mass index (BMI) was calculated, and the cut-off for overweight was set to ≥25 kg/m2 according to the WHO classification (33). The waist-to-hip ratio (WHR) was calculated, and the cut-off for central obesity was set to >0.85 since WHO recommends that a woman’s WHR is ≤0.85 (34). Breast size was measured with plastic cups, as described previously (35). Clinical information, including medication use, was retrieved from medical records, and combined with information from questionnaires.
Information on clinical tumor markers, such as ER and progesterone receptor (PR) expression (cut-off at >10% positively stained nuclei according to current clinical routine in Sweden), was collected from pathology reports. HER2 status (amplified/non-amplified) was incorporated into a clinical routine in November 2005 and was available for most tumors after that date. Clinical HER2 status was supplemented with retrospective HER2 analysis by gene protein assay on tissue microarrays (TMA) (36). Information on histological type and grade, invasive tumor size, and axillary node involvement (ALNI) was retrieved from the patient charts and pathology reports.
Patients were followed until June 30th, 2019. Information on survival and breast cancer events was obtained from the Swedish Population Registry, the Regional Tumor Registry, pathology reports, and patient charts. Local or regional recurrences, contralateral cancers, or distant metastasis and death due to any cause were considered as endpoints in recurrence-free survival (RFS) analyses. For analyses of distant metastasis-free survival (DMFS), both distant metastasis and death due to any cause were used as an endpoint. Patients who emigrated were followed until the first breast cancer event, if occurring before emigration, and otherwise until the last follow-up before emigration. For all patients living in Sweden, information on subsequent death was obtained from the Swedish Population Registry. Patients without events were censored at the time of the last follow-up before emigration. Other patients were censored at the time of the last follow-up by June 30th, 2019.
Adjuvant breast cancer treatment was administered according to clinical recommendations and was only considered before the first breast cancer event. In patients without any breast cancer events, treatments were recorded until the last follow-up or death before July 1st, 2019. Written informed consents were obtained from all participating patients, and the study was approved by the local ethics committee at Lund University (Dnr 75-02, Dnr 37-08, Dnr 658-09, Dnr 58-12, Dnr 379-12, Dnr 227-13, Dnr 277-15, and Dnr 458-15). The study adhered to Reporting Recommendations for Tumor Markers Prognostic Studies (REMARK) criteria (37).
Tissue Microarray Construction and Immunohistochemistry
Duplicate 1-mm cores from representative tumor regions of formalin-fixed paraffin-embedded tissue blocks were collected from surgical specimens and assembled in a TMA block using a semi-automated tissue array device (Beecher Instruments Inc., Sun Prairie, WI). The TMA blocks were stored at room temperature before sectioning, and the 4-μm thick TMA sections were kept at -20°C until immunohistochemical staining. The sections were automatically deparaffinized before pretreatment using the PT Link system (DAKO, Glostrup, Denmark).
Immunohistochemistry was performed using the Autostainer Plus from DAKO with the EnVision FLEX high-pH kit, according to the manufacturer’s instructions (DAKO, Glostrup, Denmark). A rabbit polyclonal antibody against the AhR, diluted 1:1000 (BML-SA550, Enzo Life Sciences, Farmingdale, NY) was used. The characteristics of the AhR antibody have been reported elsewhere (22, 30). Sections were also stained with a mouse monoclonal antibody to aromatase (clone 677, provided by professor D. Edwards) at 1.9 mg/ml, diluted 1:250, for 1 h at room temperature followed by EnVision FLEX high-pH kit, in an Autostainer Plus, according to the manufacturer’s instructions. The antibody has previously been thoroughly validated for the detection of aromatase by immunohistochemistry (38, 39). Human placenta was used as a positive control for aromatase.
AhR staining was evaluated by two independent observers (SBj and AB), without knowledge of tumor characteristics and patient information. In case of discrepancy, a re-examination was done until consensus was reached or a senior pathologist (KJ) was consulted. Scoring included cytoplasmic staining intensity score (AhRcyt); negative, weak, moderate, strong, and nuclear staining (AhRnuc); negative, positive. AhRcyt was negative in only 12 tumors (1.3%), and this group was therefore combined with the group with weak staining (n=171) to form a category denoted as low (n=183), and the groups with moderate (n=413) and strong (n=324) staining were combined into a group with high intensity (n=737), thus creating a dichotomized variable for the statistical analyses. In the case of several different AhRcyt intensities within the same tumor cores, the highest intensity with a fraction of >20% of the invasive tumor cells was selected.
Similarly, intratumoral aromatase staining was evaluated by two independent observers (ES and HT), and KJ was consulted in case of discrepancy. Only cytoplasmic staining of aromatase in invasive cells was evaluated. Each tumor was assigned an intensity score (negative, weak, moderate, strong) and a percentage of stained cells. Since intratumoral aromatase was positive in at least 1% of the cells in only 69 cases (8.5%), a dichotomized variable (negative, positive) was used for all analyses.
For both markers, duplicate cores were evaluated jointly. In the case of bilateral tumors (n=17), 12 had an evaluable AhR staining. AhRcyt staining differed on the contralateral side in two cases, but both remained in the same category. AhRnuc staining also differed in two cases, and these changed category. Intratumoral aromatase did not differ in any of the bilateral cases. Data from the tumor on the side with the highest AhR level were used, and all tumor characteristics were taken from the corresponding side. Since only two cases differed, no sensitivity analysis was performed.
Genomic DNA was obtained for buffy coats from the patients’ leukocyte portion of frozen peripheral blood using the Wizard Genomic DNA Purification Kit (Promega, Madison, USA). Genotyping was performed at the Region Skåne Competence Centre (RSKC Malmö), Skåne University Hospital, Malmö, Sweden, according to the manufacturer’s instructions with reagents included in the iPLEX™ genotyping kit (Sequenom, Inc., San Diego, CA, USA) and the software and equipment in the MassARRAY® platform (Sequenom, Inc., San Diego, CA, USA). The functional AhR SNPArg554Lys (rs2066853) and the CYP19A1 SNPs (rs10046), (rs4646), Aro1 (rs4775936), and Aro2 (rs10459592) were analyzed. In case of missing values for CYP19A1, most could be imputed based on the other SNPs as previously described (3). Haplotypes of CYP19A1 were constructed by cross-tabulation of the genotypes of the CYP19A1 SNPs, which resulted in nine haplotypes assembled into four diplotypes and a combined group of rare diplotypes (<10%) as previously described (3). Genotyping was performed in 2008 and was only available for patients included between 2002 and 2008 (n=576).
All statistical analyses were performed using SPSS software version 26 (IBM Corp, Armonk, NY, USA). Variables were dichotomized as follows: age (≥50 years), where age <50 years was used as a proxy for premenopausal status, BMI (≥25 kg/m2), waist circumference (≥80 cm), WHR (>0.85), breast volume (≥850 ml), nulliparous, ever use of oral contraceptives (OCs), ever use of menopausal hormone therapy (MHT), coffee intake (≥2 cups/day), current smoker before surgery, alcohol abstainer, adjuvant treatments (chemotherapy, radiotherapy, tamoxifen, AI, and/or herceptin) before any breast cancer event, last follow-up, or death.
Tumor characteristics included invasive pathologic tumor size (≤20 mm, >20 mm, or skin or muscular involvement), any ALNI, histological grade (grade I–III), hormone receptor status (ER, PR), HER2 amplification, and/or triple negativity.
Patient and tumor characteristics were analyzed in relation to AhR levels (AhRcyt; high versus low, AhRnuc; positive versus negative) and intratumoral aromatase expression (positive versus negative). Chi-Square tests were used for dichotomous variables, and test for trends was assessed by linear-by-linear association tests. Correlations between staining intensity and time between surgery and staining (years) were calculated with Spearman’s Rho (Rs). Staining intensity of AhRcyt and AhRnuc was negatively correlated with the time between surgery and staining (P<0.0001). Therefore, logistic regression models with patient and tumor characteristics, as well as other analyses, were performed with adjustments for the time between surgery and staining for AhRcyt and AhRnuc. Staining of intratumoral aromatase was not correlated with the time between surgery and staining (P=0.76).
The impact of AhRcyt, AhRnuc, and intratumoral aromatase on RFS and DMFS was estimated with Kaplan-Meier curves and assessed with the LogRank test. Groups were formed according to the expression of AhR in different compartments as follows: group 1, high AhRcyt and negative AhRnuc; group 2, high AhRcyt and positive AhRnuc; group 3, low AhRcyt and positive AhRnuc; and group 4, low AhRcyt and negative AhRnuc.
Cox regression models were used for multivariable analyses providing hazard ratios (HRs) with 95% confidence intervals (CI). Adjustments were performed in two different models: Model 1, the time between surgery and staining (years; continuous), age (years; continuous), BMI ≥25 kg/m2, and tumor characteristics (tumor size >20 mm or skin/muscular involvement irrespective of size, histological grade III, ER status); and Model 2: model 1 with the addition of preoperative smoking, alcohol abstention, and adjuvant treatments (chemotherapy, radiotherapy, tamoxifen, and AI). When comparing four groups of AhR levels in different cellular compartments, as previously described, group 1 (high AhRcyt and negative AhRnuc) was used as a reference since it was the largest group.
To examine whether there were any effect modifications by age ≥50 years at inclusion, preoperative BMI ≥25 kg/m2, WHR >0.85, tumor characteristics (tumor size >20 mm or skin/muscular involvement, any ALNI, histological grade III, or ER status), preoperative smoking, alcohol abstention, adjuvant treatments (chemotherapy, radiotherapy, tamoxifen and AI) on the associations between AhRcyt and AhRnuc and prognosis, multiplicative interaction variables between these factors and the categories of AhRcyt and AhRnuc were calculated. Interaction analyses were adjusted according to model 1.
Power calculations including 900 patients, of which 20% had low AhRcyt levels and 33% had positive AhRnuc status, with an accrual interval of 10 years and additional follow-up time of seven years, 80% power, and α of 0.05, showed that with a mean survival time of nine years it was possible to detect true HRs of ≤0.75 or ≥1.38 for AhRcyt and of ≤0.78 or ≥1.31 for AhRnuc. Power calculations, including 804 patients, of which 8.3% had intratumoral aromatase, showed that it was possible to detect true HRs of ≤0.65 or ≥1.66. The power calculations were performed with the PS Power and Sample Size Calculation Program, version 3.1.2 (40). All P-values were two-sided, and each P-value should be interpreted as the level of evidence against each null hypothesis. Since this is an exploratory study, nominal P-values are presented without adjustments for multiple testing (41).
Patient characteristics are presented in Table 1. The median age at inclusion was 61 years (range, 24–99 years). All P-values were adjusted for the time between surgery and staining in AhR analyses. Patients with high AhRcyt levels (n=737/920) had somewhat larger WHR (Padj=0.046). Positive AhRnuc status (n=292/920) was associated with alcohol abstention (Padj=0.046). Otherwise, patient characteristics were similar between patients with different levels of AhRcyt and AhRnuc. Patients with positive intratumoral aromatase (n=69/816) were younger (P=0.001) and more likely to be nulliparous (P=0.003) than those with negative status.
Tumor characteristics and clinical data are presented in Table 2. AhRcyt levels and AhRnuc status were positively correlated (rs=0.22, P<0.0001; Figure 1). High AhRcyt levels were associated with several favorable tumor characteristics, such as lower histological grade, ER+ and PR+ status, and lower frequency of triple-negative tumors (all adjusted Ps ≤ 0.0004). In contrast, positive AhRnuc status was not associated with tumor characteristics apart from the lower frequency of PR+ tumors (Padj=0.036). The highest frequency of ER+ tumors (92.7%) was found in patients with high AhRcyt and negative AhRnuc status, while the lowest frequency of ER+ tumors (75.0%) was observed in tumors with low AhRcyt and positive AhRnuc status (P<0.0001).
Table 2 Tumor characteristics, treatments, and events in relation to AhRcyt, AhRnuc, and intratumoral aromatase levels.
A lower proportion of patients with high compared with low AhRcyt levels received adjuvant chemotherapy (Padj=0.0004). Adjuvant treatment with tamoxifen was associated with low AhRcyt levels (P=0.010) and positive AhRnuc status (P=0.0032) while radiotherapy or AI treatment showed no association with either AhRcyt or AhRnuc status.
Intratumoral aromatase status was not correlated with either AhRcyt or AhRnuc status (see Figure 1). Positive intratumoral aromatase status was strongly associated with ER– and PR– status, triple negativity, and higher histological grade (all Ps ≤ 0.0002). In line with this, positive intratumoral aromatase status was associated with adjuvant chemotherapy (P=0.001) but not with adjuvant endocrine therapy or radiotherapy.
Tumors With Non-Evaluable Staining
Patients whose tumors could not be evaluated for AhR (n=98) were somewhat younger than included patients and had smaller tumors and less ALNI than patients with evaluable AhR. Also, they had more HER2+ tumors, no difference in ER and PR status, but received less adjuvant endocrine therapy. Patients whose tumors could not be evaluated for intratumoral aromatase expression (n=202) had less ALNI and lower histological grade, no difference in hormone receptors, but received less endocrine therapy than patients with evaluable aromatase expression.
AhR and CYP19A1 Genotypes and AhR and Intratumoral Aromatase Status
There was no correlation between the AhRArg554Lys and any of the CYP19A1 genotypes. The functional AhRArg554Lys genotype (n=576) was not clearly associated with either AhRcyt or AhRnuc levels (both Ps≥0.088). None of the CYP19A1 genotypes or diplotypes (n=575) were associated with intratumoral aromatase. However, the CYP19A1 diplotypes were associated with AhRnuc status (P=0.033, 4 degrees of freedom (d.f.)). Rare CYP19A1 diplotypes were more common in patients with positive AhRnuc tumors.
Follow-up and Events
Patients received follow-up questionnaires for up to 15 years from inclusion. For the 920 patients with invasive cancer included in the AhR survival analyses, 660 patients were still at risk by the end of follow-up, and for these patients, the median follow-up time was 8.7 years (IQR 6.9–11.0). During the follow-up time, breast cancer events occurred in 178 patients, of whom 113 had distant metastases. Also, 176 patients died, and 94 of these had a prior breast cancer event.
Intratumoral AhR Levels in Relation to Prognosis
High AhRcyt levels were associated with a better 10-year prognosis in terms of both RFS and DMFS in univariable analyses. The curves merged by the 15-year follow-up (Figures 2A, B). In multivariable analyses, high AhRcyt levels were associated with approximately half the risk of an event by the 5-year follow-up, HRadj 0.55 (95% CI 0.36–0.84), and 5-year RFS of patients with high AhRcyt levels was 89.6% compared with 79.6% for patients with low levels. The association became weaker by the 10-year follow-up, HRadj 0.72 (95% CI 0.52–1.00). The results were mainly the same for DMFS, with 5-year DMFS of 92.5% in patients with high AhRcyt levels compared with 82.9% in patients with low levels. In contrast, AhRnuc was not associated with RFS or DMFS (Figures 2C, D).
Figure 2 High AhRcyt levels compared to low AhRcyt in relation to 5-year, 10-year, and the entire follow-up: (A) recurrence-free survival and (B) distant metastasis-free survival. Positive AhRnuc levels compared to negative AhRnuc status in relation to 5-year, 10-year, and the entire follow-up (C) recurrence-free survival and (D) distant metastasis-free survival. Associations (3 d.f.) between four groups of AhRcyt levels and AhRnuc status (high AhRcyt and positive AhRnuc; high AhRcyt and negative AhRnuc; low AhRcyt and positive AhRnuc; low AhRcyt and negative AhRnuc) and 5-year, 10-year and the entire follow-up (E) recurrence-free survival and (F) distant metastasis-free survival.
Patients were divided into four groups depending on the subcellular localization of AhR (Table 3). By the 10-year follow-up, patients with low AhRcyt and negative AhRnuc had the poorest prognosis, both in terms of RFS, HRadj 1.65 (95% CI 1.16–2.33), and DMFS, HRadj 1.69 (95% CI 1.12–2.55). Beyond that, no clear associations were found (Figures 2E, F). When only considering patients with more than five years of follow-up, those with high AhRcyt and positive AhRnuc had a tendency towards more late recurrences.
Table 3 Multivariable models with crude and adjusted HR (95% CIs) for combined AhRcyt and AhRnuc in relation to RFS and DMFS for 5-yr, 10-yr, and the entire follow-up.
Effect Modifications Between AhRcyt Levels and Endocrine Therapy on Prognosis
Formal interaction analyses were conducted to study effect modifications between patient and tumor characteristics, treatments, and the AhRcyt and AhRnuc on prognosis. No effect modifications were found between AhRcyt levels and patients or tumor characteristics. However, an interaction was found between AhRcyt and any AI treatment (Pinteraction=0.030). No interactions were found between AhRcyt and chemotherapy, radiotherapy, or tamoxifen. Several of the patients who received any AI treatment also had received tamoxifen (switch-therapy). Therefore, three new variables were formed; tamoxifen only, AI only, and switch-therapy.
Interaction analyses with patients with ER+ tumors were conducted, including these variables. In this model, the effect modification seemed driven by switch-therapy rather than by AI only (Pinteraction=0.82). The prognostic impact of AhRcyt in endocrine treatment groups is shown in Figures 3A, F, G, J. High AhRcyt levels were strongly associated with good 10-year prognosis in tamoxifen-treated patients, HRadj 0.40 (95% CI 0.23–0.71), 10-year RFS was 78.2% in patients with high AhRcyt levels compared with 52.9% in patients with low levels, while almost an inverse association was seen in patients who received switch-therapy (Pinteraction=0.023).
Figure 3 High AhRcyt compared to low AhRcyt levels in patients with ER+ tumors in relation to 5-year, 10-year and the entire follow-up recurrence-free survival in patients with (A) no endocrine therapy, (D) tamoxifen only, (G) aromatase inhibitors only, and (J) switch-therapy. Positive AhRnuc levels compared to negative AhRnuc status in relation to 5-year, 10-year, and the entire follow-up recurrence-free survival in patients with (B) no endocrine therapy, (E) tamoxifen only, (H) aromatase inhibitors only, and (K) switch-therapy. AhRArg554Lys GG compared any A genotype in relation to 5-year, 10-year and the entire follow-up recurrence-free survival in patients with (C) no endocrine therapy, (F) tamoxifen only, (I) aromatase inhibitors only, and (L) switch-therapy.
Effect Modifications Between AhRnuc Status and Clinicopathological Factors on Prognosis
An interaction was found between AhRnuc and age ≥50 years on RFS where younger patients with positive AhRnuc status had a better prognosis compared to patients with negative AhRnuc status, HRadj=0.40 (95% CI 0.17–0.99), an association not seen in older patients (Pinteraction=0.011). Another interaction between AhRnuc and any use of chemotherapy was found where chemonaïve patients with positive AhRnuc status had a poorer prognosis than patients with negative AhRnuc status, HRadj=1.54 (95% CI 1.06–2.25), while no difference was seen in chemotherapy-treated patients (Pinteraction=0.008). No interaction was found between AhRnuc and any adjuvant use of radiotherapy, tamoxifen, or AI.
Similarly, we conducted an additional interaction analysis in patients with ER+ tumors with tamoxifen only, AI only, and switch-therapy. In this analysis, an interaction was found between AhRnuc and the use of tamoxifen only (Pinteraction=0.041) on RFS, but not with AI only or switch-therapy (Figures 3B, E, H, K).
Effect Modifications Between AhRArg554Lys and Endocrine Therapy on Prognosis
In line with our previously reported findings, switch-treated patients with AhRArg554Lys (GG) genotype had a better prognosis than any A genotypes, HRadj 0.42 (95% CI 0.22–0.83). This association was strongest during a 10-year follow-up (Figure 3L). No clear associations were found in other treatment groups (Figures 3C, F, I). However, no significant effect modifications were found between the AhRArg554Lys genotype and endocrine therapy.
AhR and Intratumoral Aromatase Status in Relation to Prognosis
Overall, intratumoral aromatase was not associated with RFS (LogRank P=0.43) or DMFS (LogRank P=0.17). Since there were few patients with tumors with positive aromatase status, no further analyses were conducted.
The main findings in this study were that high AhRcyt levels were associated with favorable tumor characteristics and clinical outcomes in primary breast cancer patients, while AhRnuc status was not associated with overall prognosis. The prognostic impact of AhRcyt was substantially modified by the type of endocrine therapy. The largest difference was seen between tamoxifen only and switch-therapy. Likewise, the prognostic impact of AhRArg554Lys genotype also differed between endocrine treatment groups. AhRnuc status was prognostic in chemonaïve but not in chemotherapy-treated patients. Intratumoral aromatase was associated with several aggressive tumor characteristics but not with AhR levels or prognosis. Neither AhR levels nor intratumoral aromatase status was associated with their respective genotypes.
Associations between favorable tumor characteristics and AhR levels seen in the present study have previously been reported in some (22) but not all (23, 42) studies. In line with our findings, two previous studies reported good prognosis with higher AhR levels (5, 22) but another study reported no association (19). In contrast, a couple of studies suggested associations between high AhR expression and poor prognosis (21, 23). To our knowledge, only one previous study evaluated the impact of AhRcyt and AhRnuc levels separately and found a negative prognostic impact of higher AhRnuc levels in lymph node-negative patients, who are less likely to receive chemotherapy, while AhRcyt levels were not associated with prognosis (23). However, no formal interaction analysis was presented in the paper. Our results do not support their finding since no effect modifications between ALNI and AhR levels on prognosis were found.
Patients with high AhRcyt and negative AhRnuc status had the best prognosis during the entire follow-up period in the present study. Therefore, we hypothesize that it is predominantly the inactivated AhRcyt that confers a good prognosis in breast cancer. Although patients with high AhRcyt and positive AhRnuc status had good initial outcomes during the first five years, when endocrine therapy is given to most patients with ER+ tumors, late recurrences were more common in this group. These patients might benefit from extended endocrine therapy. It is possible that a ligand-activated, translocated AhRnuc and the downstream transcription of multiple pathways confers a negative prognostic impact that is only apparent after completion of endocrine therapy. The worst initial outcomes were associated with low AhRcyt and negative AhRnuc status. Since only a handful of patients had positive AhRnuc status and low AhRcyt levels, we could not draw any conclusions regarding this combination.
In vitro studies suggest that AhR activation induces intratumoral aromatase, which in turn increases intratumoral estrogen synthesis and proliferation (43, 44) while at the same time inhibiting the ER pathway in breast cancer cells (27, 45, 46). Although no overall correlation was found between intratumoral aromatase status and AhR levels in the present study, an exploratory analysis revealed a marginally positive correlation between intratumoral aromatase and AhRcyt levels in patients with ER+ but not ER– tumors (data not shown). However, the number of patients was small, and the level of evidence was low.
Interesting interactions were found between AhRcyt and AhRnuc status and endocrine therapy in relation to RFS. To our knowledge, this is the first study suggesting an interplay between AhR levels and the type of endocrine therapy in breast cancer. Patients with high AhRcyt levels who only received endocrine therapy with tamoxifen had an excellent prognosis during the first 10 years of follow-up, whereas an almost inverse relationship was seen for patients treated with switch-therapy between tamoxifen and AIs. Further, switch-treated patients with the AhRArg554Lys GG genotype had favorable outcomes compared to patients with any A genotype. No prognostic impact of the AhRArg554Lys genotype was seen in other endocrine treatment groups. The homozygous AA genotype has been associated with lower mRNA expression of AhR (11), but in the present study, no correlation between genotypes and protein levels were found.
Positive intratumoral aromatase status was only found in 8.5% of the patients in the cohort. This is in contrast with previous findings in a small study using the same monoclonal antibody, where intratumoral aromatase staining was denoted in 17 of 28 tumors (47). However, aromatase can be heterogeneously expressed in the tumor (30, 47) and the use of TMA instead of whole tumor sections in the present study might have yielded an underestimation of the intratumoral aromatase status (48). Intratumoral aromatase status was associated with hormone receptor negativity and other aggressive tumor characteristics but not AhR in the current study. Thus, we could not confirm the previous in vitro findings that AhR induces intratumoral aromatase in breast cancer (30), but due to the limited number of aromatase-positive tumors, these findings should be interpreted with caution.
In contrast to our hypothesis that anthropometric factors would influence tumor levels of AhR and aromatase, only a marginal association between WHR and AhR was found. Also, a marginal inverse association with alcohol abstention (49) but not smoking (7) was observed. Younger and nulliparous patients were substantially more likely to have tumors positive for intratumoral aromatase, both factors associated with ER– status (50). ER– status was strongly associated with intratumoral aromatase status in the present study, which is in contrast to a previous report (48). Additional staining for classical AhR targets such as CYP1A1 and CYP1B1 to confirm the activity status of AhR would have been interesting but was outside the scope of this study.
During the inclusion period, more than half of the primary breast cancer patients undergoing breast surgery in Lund were included in the cohort (32). The main reason for non-inclusion was a limited number of research nurses. The included patients did not differ substantially from all breast cancer patients operated in Lund during this period (32). Therefore, the findings can be considered generalizable for breast cancer patients treated at this clinic. However, evaluable TMA cores were slightly more often missing for patients with less advanced tumors than for patients with more advanced tumors, but this differed somewhat between the two markers. Missing evaluable tumor tissue could potentially cause selection bias in the analyses. Prognostic associations for AhR became stronger when HER2 was included in the multivariable models (data not shown). However, HER2 status was more often missing from small tumors with longer follow-up times (36) and inclusion of this marker could bias the analyses. Therefore, HER2 status was not used in the multivariable analyses presented. MHT was not associated with either marker. Since we have previously reported that MHT had no overall prognostic impact in this cohort (51), MHT was not included in the adjustment models. In this observational study, we took commonly used prognosticators into account in the multivariable models to evaluate the independent prognostic value of each marker, but residual confounding remains possible. We also adjusted for the time between surgery and staining since AhR levels were lower the longer the tumors had been stored. This time variable also captures other changes occurring during the period, such as changes in treatments. Cox proportional hazards models were used. However, since not all hazards were proportional, we also divided the follow-up period into 5-year and 10-year follow-ups. Although statistical power was decent in the entire cohort, the number of patients and events in some subgroups were smaller, which led to lower power in these analyses.
In conclusion, our results suggest that AhRcyt levels are associated with ER+ status and other favorable tumor characteristics and prognosis. The prognostic impact of AhR was substantially modified by the type of endocrine therapy where high AhRcyt levels were associated with significantly longer recurrence-free survival during the first 10 years of follow-up among patients who received tamoxifen as only endocrine therapy compared to patients with low AhRcyt levels, whereas an almost inverse relationship was seen in patients treated with switch-therapy. Our findings suggest that both AhRcyt levels and AhRArg554Lys genotypes merit further study as to whether they can be used to guide the selection of endocrine therapy for patients with ER+ breast cancer in the clinical setting. These findings warrant confirmation in an independent cohort, preferably in a randomized clinical trial comparing different regimens of tamoxifen and AIs. They might also guide the selection of breast cancer patients for inclusion in trials with selective AhR modulators.
Data Availability Statement
The data analyzed in this study is subject to the following licenses/restrictions: Data is not publicly available due to privacy laws. Questions regarding data should be directed to email@example.com.
The studies involving human participants were reviewed and approved by Lund University Ethics committee. The patients/participants provided their written informed consent to participate in this study.
HT and HJ: study design, manuscript preparation, statistical analysis, data analysis, and interpretation. HJ, KI, and SBo: study supervision. HT, ES, SBj, AB, MR, SK, BN, KJ, and HJ: data collection. DE: provision of the anti-aromatase 677 antibody. HT, ES, SBj, AB, MR, SK, DE, BN, KJ, KI, SBo, and HJ: contributed to the manuscript review, critical revision for important intellectual content, and read and approved the final draft for submission. All authors are responsible for the manuscript content.
This work was supported by grants from The Swedish Cancer Society (CAN2017/368), the Medical Faculty at Lund University, the Mrs. Berta Kamprad Foundation, the South Swedish Health Care Region (Region Skåne ALF 40620) Region Skåne ST-ALF, and the Skåne University Hospital Fund. Open access funding provided by Lund University. The funding agencies played no role in the study’s design, the collection, analysis, and interpretation of the data, the writing of the manuscript, or in the decision to submit the manuscript for publication.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank our research nurses Linda Ågren, Helén Thell, Jessica Åkesson, Anette Ahlin Gullers, Monika Eberhard Mészaros, Maj-Britt Hedenblad, Karin Henriksson, and Anette Möller. We also thank Erika Bågeman, Maria Henningson, Maria Hjertberg, and Andrea Markkula for data entry. We thank Professor H. Sasano for insightful comments regarding staining with the aromatase antibody. We thank Kurt Christensen and the Protein and Monoclonal Antibody Production Core at Baylor College of Medicine for the production and purification of the aromatase monoclonal antibody.
AhR, aryl hydrocarbon receptor; AhRcyt, cytoplasmic aryl hydrocarbon receptor; AhRnuc, nuclear aryl hydrocarbon receptor; AI, aromatase inhibitor; ALNI, axillary lymph node involvement; ARNT, aryl hydrocarbon receptor nuclear translocator; BC, breast cancer; BMI, body mass index; CI, confidence interval; CYP, cytochrome P450; DFS, disease-free survival; DMFS, disease metastasis-free survival; DNA, deoxyribonucleic acid; ER, estrogen receptor; HER2, human-epidermal growth factor receptor 2; HR, hazard ratio; IGF, insulin-like growth factor; MHT, menopausal hormone therapy; OC, oral contraceptives; PR, progesterone receptor; REMARK, recommendations for Tumor Marker Prognostic Studies; SERM, selective estrogen receptor modulator; TMA, tumor tissue microarray; WHR, waist-to-hip ratio.
1. Li N, Deng Y, Zhou L, Tian T, Yang S, Wu Y, et al. Global Burden of Breast Cancer and Attributable Risk Factors in 195 Countries and Territories, From 1990 to 2017: Results From the Global Burden of Disease Study 2017. J Hematol Oncol (2019) 12(1):140. doi: 10.1186/s13045-019-0828-0
2. Burstein HJ, Curigliano G, Loibl S, Dubsky P, Gnant M, Poortmans P, et al. Estimating the Benefits of Therapy for Early-Stage Breast Cancer: The St. Gallen International Consensus Guidelines for the Primary Therapy of Early Breast Cancer 2019. Ann Oncol (2019) 30(10):1541–57. doi: 10.1093/annonc/mdz235
3. Simonsson M, Veerla S, Markkula A, Rose C, Ingvar C, Jernstrom H. CYP1A2–a Novel Genetic Marker for Early Aromatase Inhibitor Response in the Treatment of Breast Cancer Patients. BMC Cancer (2016) 16:256. doi: 10.1186/s12885-016-2284-3
5. O’Donnell EF, Koch DC, Bisson WH, Jang HS, Kolluri SK. The Aryl Hydrocarbon Receptor Mediates Raloxifene-Induced Apoptosis in Estrogen Receptor-Negative Hepatoma and Breast Cancer Cells. Cell Death Dis (2014) 5:e1038. doi: 10.1038/cddis.2013.549
7. Wang GZ, Zhang L, Zhao XC, Gao SH, Qu LW, Yu H, et al. The Aryl Hydrocarbon Receptor Mediates Tobacco-Induced PD-L1 Expression and is Associated With Response to Immunotherapy. Nat Commun (2019) 10(1):1125. doi: 10.1038/s41467-019-08887-7
8. Rothhammer V, Quintana FJ. The Aryl Hydrocarbon Receptor: An Environmental Sensor Integrating Immune Responses in Health and Disease. Nat Rev Immunol (2019) 19(3):184–97. doi: 10.1038/s41577-019-0125-8
9. Beischlag TV, Luis Morales J, Hollingshead BD, Perdew GH. The Aryl Hydrocarbon Receptor Complex and the Control of Gene Expression. Crit Rev Eukaryot Gene Expr (2008) 18(3):207–50. doi: 10.1615/critreveukargeneexpr.v18.i3.20
10. Brauze D, Zawierucha P, Kiwerska K, Bednarek K, Oleszak M, Rydzanicz M, et al. Induction of Expression of Aryl Hydrocarbon Receptor-Dependent Genes in Human HepaRG Cell Line Modified by shRNA and Treated With Beta-Naphthoflavone. Mol Cell Biochem (2017) 425(1-2):59–75. doi: 10.1007/s11010-016-2862-3
11. Helmig S, Seelinger JU, Dohrel J, Schneider J. RNA Expressions of AHR, ARNT and CYP1B1 are Influenced by AHR Arg554Lys Polymorphism. Mol Genet Metab (2011) 104(1-2):180–4. doi: 10.1016/j.ymgme.2011.06.009
12. Simonsson M, Markkula A, Rose C, Ingvar C, Jernström HC. Abstract LB-194: Ahr Arg554Lys Impacts Tumor ER Expression and Endocrine Treatment Response in Breast Cancer. Cancer Res (2016) 76(14 Supplement):LB–194-LB-194. doi: 10.1158/1538-7445.Am2016-lb-194
13. Colomer R, Monzo M, Tusquets I, Rifa J, Baena JM, Barnadas A, 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–6. doi: 10.1158/1078-0432.CCR-07-1923
14. Park IH, Lee YS, Lee KS, Kim SY, Hong SH, Jeong J, et al. Single Nucleotide Polymorphisms of CYP19A1 Predict Clinical Outcomes and Adverse Events Associated With Letrozole in Patients With Metastatic Breast Cancer. Cancer Chemother Pharmacol (2011) 68(5):1263–71. doi: 10.1007/s00280-011-1615-y
15. Liu L, Bai YX, Zhou JH, Sun XW, Sui H, Zhang WJ, et al. A Polymorphism At the 3’-UTR Region of the Aromatase Gene is Associated With the Efficacy of the Aromatase Inhibitor, Anastrozole, in Metastatic Breast Carcinoma. Int J Mol Sci (2013) 14(9):18973–88. doi: 10.3390/ijms140918973
16. Gervasini G, Jara C, Olier C, Romero N, Martinez R, Carrillo JA. Polymorphisms in ABCB1 and CYP19A1 Genes Affect Anastrozole Plasma Concentrations and Clinical Outcomes in Postmenopausal Breast Cancer Patients. Br J Clin Pharmacol (2017) 83(3):562–71. doi: 10.1111/bcp.13130
17. Leyland-Jones B, Gray KP, Abramovitz M, Bouzyk M, Young B, Long B, et al. CYP19A1 Polymorphisms and Clinical Outcomes in Postmenopausal Women With Hormone Receptor-Positive Breast Cancer in the BIG 1-98 Trial. Breast Cancer Res Treat (2015) 151(2):373–84. doi: 10.1007/s10549-015-3378-3
19. Vacher S, Castagnet P, Chemlali W, Lallemand F, Meseure D, Pocard M, et al. High AHR Expression in Breast Tumors Correlates With Expression of Genes From Several Signaling Pathways Namely Inflammation and Endogenous Tryptophan Metabolism. PloS One (2018) 13(1):e0190619. doi: 10.1371/journal.pone.0190619
22. Saito R, Miki Y, Hata S, Takagi K, Iida S, Oba Y, et al. Aryl Hydrocarbon Receptor in Breast Cancer-a Newly Defined Prognostic Marker. Horm Cancer (2014) 5(1):11–21. doi: 10.1007/s12672-013-0160-z
23. Jeschke U, Zhang X, Kuhn C, Jalaguier S, Colinge J, Pfender K, et al. The Prognostic Impact of the Aryl Hydrocarbon Receptor (AhR) in Primary Breast Cancer Depends on the Lymph Node Status. Int J Mol Sci (2019) 20(5):1016. doi: 10.3390/ijms20051016
25. Helle J, Bader MI, Keiler AM, Zierau O, Vollmer G, Chittur SV, et al. Cross-Talk in the Female Rat Mammary Gland: Influence of Aryl Hydrocarbon Receptor on Estrogen Receptor Signaling. Environ Health Perspect (2016) 124(5):601–10. doi: 10.1289/ehp.1509680
27. Wormke M, Stoner M, Saville B, Walker K, Abdelrahim M, Burghardt R, et al. The Aryl Hydrocarbon Receptor Mediates Degradation of Estrogen Receptor Alpha Through Activation of Proteasomes. Mol Cell Biol (2003) 23(6):1843–55. doi: 10.1128/mcb.23.6.1843-1855.2003
29. DuSell CD, Nelson ER, Wittmann BM, Fretz JA, Kazmin D, Thomas RS, et al. Regulation of Aryl Hydrocarbon Receptor Function by Selective Estrogen Receptor Modulators. Mol Endocrinol (2010) 24(1):33–46. doi: 10.1210/me.2009-0339
30. Saito R, Miki Y, Hata S, Ishida T, Suzuki T, Ohuchi N, et al. Aryl Hydrocarbon Receptor Induced Intratumoral Aromatase in Breast Cancer. Breast Cancer Res Treat (2017) 161(3):399–407. doi: 10.1007/s10549-016-4063-x
32. Persson M, Simonsson M, Markkula A, Rose C, Ingvar C, Jernström H. Impacts of Smoking on Endocrine Treatment Response in a Prospective Breast Cancer Cohort. Br J Cancer (2016) 115(3):382–90. doi: 10.1038/bjc.2016.174
33. World Health Organization. Obesity: preventing and managing the global epidemic of obesity. Geneva, Switzerland (1997). Available at: https://apps.who.int/iris/handle/10665/63854.
34. World Health Organization. Waist Circumference and Waist-Hip Ratio: report of a WHO expert consultation (2008). Available at: https://apps.who.int/iris/bitstream/handle/10665/44583/9789241501491_eng.pdf?ua=1.
35. Markkula A, Bromee A, Henningson M, Hietala M, Ringberg A, Ingvar C, et al. Given Breast Cancer, Does Breast Size Matter? Data From a Prospective Breast Cancer Cohort. Cancer Causes Control (2012) 23(8):1307–16. doi: 10.1007/s10552-012-0008-9
36. Sandén E, Khazaei S, Tryggvadottir H, Borgquist S, Isaksson K, Jirström K, et al. Re-Evaluation of HER2 Status in 606 Breast Cancers-Gene Protein Assay on Tissue Microarrays Versus Routine Pathological Assessment. Virchows Arch (2020) 477(2):317–20. doi: 10.1007/s00428-020-02768-x
37. Altman DG, McShane LM, Sauerbrei W, Taube SE. Reporting Recommendations for Tumor Marker Prognostic Studies (REMARK): Explanation and Elaboration. PloS Med (2012) 9(5):e1001216. doi: 10.1371/journal.pmed.1001216
38. Sasano H, Edwards DP, Anderson TJ, Silverberg SG, Evans DB, Santen RJ, et al. Validation of New Aromatase Monoclonal Antibodies for Immunohistochemistry: Progress Report. J Steroid Biochem Mol Biol (2003) 86(3-5):239–44. doi: 10.1016/s0960-0760(03)00363-7
39. Sasano H, Anderson TJ, Silverberg SG, Santen RJ, Conway M, Edwards DP, et al. The Validation of New Aromatase Monoclonal Antibodies for Immunohistochemistry–a Correlation With Biochemical Activities in 46 Cases of Breast Cancer. J Steroid Biochem Mol Biol (2005) 95(1-5):35–9. doi: 10.1016/j.jsbmb.2005.04.027
41. Victor A, Elsasser A, Hommel G, Blettner M. Judging a Plethora of P-Values: How to Contend With the Problem of Multiple Testing–Part 10 of a Series on Evaluation of Scientific Publications. Dtsch Arztebl Int (2010) 107(4):50–6. doi: 10.3238/arztebl.2010.0050
42. Li ZD, Wang K, Yang XW, Zhuang ZG, Wang JJ, Tong XW. Expression of Aryl Hydrocarbon Receptor in Relation to p53 Status and Clinicopathological Parameters in Breast Cancer. Int J Clin Exp Pathol (2014) 7(11):7931–7.
43. Cheshenko K, Brion F, Le Page Y, Hinfray N, Pakdel F, Kah O, et al. Expression of Zebra Fish Aromatase cyp19a and cyp19b Genes in Response to the Ligands of Estrogen Receptor and Aryl Hydrocarbon Receptor. Toxicol Sci (2007) 96(2):255–67. doi: 10.1093/toxsci/kfm003
44. Chan MY, Huang H, Leung LK. 2,3,7,8-Tetrachlorodibenzo-para-dioxin Increases Aromatase (CYP19) mRNA Stability in MCF-7 Cells. Mol Cell Endocrinol (2010) 317(1-2):8–13. doi: 10.1016/j.mce.2009.11.012
45. Klinge CM, Bowers JL, Kulakosky PC, Kamboj KK, Swanson HI. The Aryl Hydrocarbon Receptor (AHR)/AHR Nuclear Translocator (ARNT) Heterodimer Interacts With Naturally Occurring Estrogen Response Elements. Mol Cell Endocrinol (1999) 157(1-2):105–19. doi: 10.1016/s0303-7207(99)00165-3
46. Ohtake F, Takeyama K, Matsumoto T, Kitagawa H, Yamamoto Y, Nohara K, et al. Modulation of Oestrogen Receptor Signalling by Association With the Activated Dioxin Receptor. Nature (2003) 423(6939):545–50. doi: 10.1038/nature01606
47. Geisler J, Suzuki T, Helle H, Miki Y, Nagasaki S, Duong NK, et al. Breast Cancer Aromatase Expression Evaluated by the Novel Antibody 677: Correlations to Intra-Tumor Estrogen Levels and Hormone Receptor Status. J Steroid Biochem Mol Biol (2010) 118(4-5):237–41. doi: 10.1016/j.jsbmb.2009.10.010
48. Lykkesfeldt AE, Henriksen KL, Rasmussen BB, Sasano H, Evans DB, Moller S, et al. In Situ Aromatase Expression in Primary Tumor is Associated With Estrogen Receptor Expression But is Not Predictive of Response to Endocrine Therapy in Advanced Breast Cancer. BMC Cancer (2009) 9:185. doi: 10.1186/1471-2407-9-185
49. Malik DE, David RM, Gooderham NJ. Ethanol Potentiates the Genotoxicity of the Food-Derived Mammary Carcinogen PhIP in Human Estrogen Receptor-Positive Mammary Cells: Mechanistic Support for Lifestyle Factors (Cooked Red Meat and Ethanol) Associated With Mammary Cancer. Arch Toxicol (2018) 92(4):1639–55. doi: 10.1007/s00204-018-2160-9
50. Fortner RT, Sisti J, Chai B, Collins LC, Rosner B, Hankinson SE, et al. Parity, Breastfeeding, and Breast Cancer Risk by Hormone Receptor Status and Molecular Phenotype: Results From the Nurses’ Health Studies. Breast Cancer Res (2019) 21(1):40. doi: 10.1186/s13058-019-1119-y
51. Godina C, Ottander E, Tryggvadottir H, Borgquist S, Isaksson K, Jernstrom H. Prognostic Impact of Menopausal Hormone Therapy in Breast Cancer Differs According to Tumor Characteristics and Treatment. Front Oncol (2020) 10:80. doi: 10.3389/fonc.2020.00080
Keywords: breast cancer, aryl hydrocarbon receptor, intratumoral aromatase, endocrine therapy, polymorphisms, prognosis
Citation: Tryggvadottir H, Sandén E, Björner S, Bressan A, Ygland Rödström M, Khazaei S, Edwards DP, Nodin B, Jirström K, Isaksson K, Borgquist S and Jernström H (2021) The Prognostic Impact of Intratumoral Aryl Hydrocarbon Receptor in Primary Breast Cancer Depends on the Type of Endocrine Therapy: A Population-Based Cohort Study. Front. Oncol. 11:642768. doi: 10.3389/fonc.2021.642768
Received: 13 January 2021; Accepted: 22 April 2021;
Published: 20 May 2021.
Edited by:Laurence Gluch, The Strathfield Breast Centre, Australia
Reviewed by:Stephen Safe, Texas A&M University, United States
Thomas Haarmann-Stemmann, Leibniz-Institut für Umweltmedizinische Forschung (IUF), Germany
Copyright © 2021 Tryggvadottir, Sandén, Björner, Bressan, Ygland Rödström, Khazaei, Edwards, Nodin, Jirström, Isaksson, Borgquist and Jernström. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Helena Jernström, Helena.Jernstrom@med.lu.se