Abstract

Aims. The available research on the association between estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), ER−/PR+ status, and the occurrence of lung cancer subsequent to breast cancer in patients (referred to as BC-LuC) had been limited. Consequently, there is a need to examine whether ER, PR, HER2, and ER−/PR+ have independent correlations with the risk and outcomes of BC-LuC, while appropriately adjusting for other potential covariates. Methods. The present study employed a cohort design and utilized data from the Surveillance, Epidemiology, and End Results (SEER) program spanning from 2010 to 2015. The study population consisted of 683,336 individuals who were diagnosed with breast cancer (referred to as BC). Various covariates were assessed at baseline, including age, sex, race, marital status, CS tumor size, laterality, radiation, chemotherapy, months from diagnosis to treatment, breast subtype, AJCC 7th edition (2010–2015), and combined summary stage (2004+). The primary objective of this study was to investigate the association between ER, PR, HER2, ER−/PR+ status, and the risk of developing BC-LuC. Logistic regression analysis was employed to assess this association. Furthermore, multivariable Cox regression analyses were conducted to calculate adjusted hazard ratios (HRs) along with their respective 95% confidence intervals (CIs). Kaplan–Meier plots and log-rank tests were utilized to estimate the outcomes, specifically overall survival (OS), disease-specific survival (DSS), and metastasis. Results. The average age of 198,972 selected participants was 59.8 ± 13.1 years, and about 99.3% of them were female. Result of fully adjusted binary logistic regression showed PR+ and HER2+ were positively associated with lower risk BC-LuC after adjusting confounders (ORs = 0.84, 95% CI: 0.73–0.96,  = 0.011 and ORs = 0.83, 95% CI: 0.72–0.96,  = 0.012, respectively). ER+ and ER−/PR+ were detected no significant relationship with BC-LuC (ORs = 1.03, 95% CI: 0.87–1.22,  = 0.718 and ORs = 1.02, 95% CI: 0.61–1.72,  = 0.936, respectively). In subgroups analyses, the results remain stable. Multivariable Cox regression showed that BC-LuC patients with ER and PR were significantly associated with OS and DSS. However, ER, PR, HER2, and ER−/PR+ were significantly associated with OS and DSS in breast cancer patients. The relationship between ER, PR, HER2, and ER−/PR+ and metastasis in breast cancer patients was different. Conclusion. The results of this study indicated a potential correlation between PR- and HER2- status and a risk of developing BC-LuC. Furthermore, it appears that the prognosis of BC-LuC may be influenced by the presence of ER+ and PR+. Therefore, additional research is warranted to fully investigate and validate this association.

1. Introduction

Breast cancer has the first highest cancer incidence and mortality rates among women, posing a significant public health challenge and burden in both the United States and China [1, 2]. However, advancements in screening modalities, such as molybdenum target X-ray and color Doppler ultrasound, coupled with improvements in diagnosis, have contributed to an increasing number of breast cancer survivors. Numerous studies had documented a heightened risk of developing a second primary lung cancer among individuals who have overcome an initial breast cancer diagnosis [3, 4]. Notably, approximately ten percent of breast cancer patients who surpass the ten-year survival threshold are afflicted by a subsequent primary carcinoma, with lung cancer being a predominant manifestation. The emergence of a second primary lung cancer can complicate clinical assessments and necessitate aggressive interventions, further amplifying the burdens faced by breast cancer survivors.

Accurate assessment of the steroid hormone receptor status, particularly estrogen and progesterone receptors (ER and PR), is essential for effective management of invasive breast cancer. The HR status is a critical determinant of prognosis and treatment strategies for breast cancer. Currently, immunohistochemistry (IHC) represents the best method for detecting ER and PR status [5]. Given the clinical, morphological, and molecular heterogeneity of breast cancer, a wide range of variations in its treatment and management may occur [6]. Notably, the St. Gallen surrogated for breast cancer subtypes has identified five primary entities based on the molecular classification through immunohistochemical expression of ER, PR, HER2, and proliferation index Ki-67: luminal-A, luminal-B, HER2-negative, luminal-B HER2-positive, HER2 enriched, and TNBC (triple negative: lack of expression of ER, PR, and no overexpression of HER2) [7, 8]. While some studies had categorized ER−/PR+ breast cancer as technical artifacts in IHC, previous research illustrated that this is a unique subtype with distinct molecular and clinical features [9]. Additionally, Gamrani et al.’s study [10] demonstrated a statistically significant association between the ER−/PR+ group and a higher risk of recurrence and death, placing it midway between the double-negative and double-positive HR subtypes. After endocrine therapy in breast cancer patients, those with ER−/PR+ exhibited higher disease-free survival rates and lower all-cause mortality than a modified group without endocrine treatment [11]. In a meta-analysis of early-stage breast cancer patients treated with endocrine therapy, ER+ breast cancer patients did not appear to benefit from endocrine treatment [12].

However, findings from previous studies regarding the relationship between ER, PR, HER2, and ER−/PR+ and lung cancer subsequent in breast cancer patients was limited. In this study, we aimed to explore the relationship between ER, PR, HER2, and ER−/PR+ and the risk and outcomes of BC-LuC. To this end, we analyzed 198,972 breast cancer patients which had been previously diagnosed in SEER database from 2010 to 1015.

2. Participants and Methods

2.1. Database

A retrospective analysis was conducted utilizing the Surveillance, Epidemiology, and End Results (SEER) database, which was publicly released in November 2021. The study population consisted of over 600,000 breast cancer patients diagnosed between 2010 and 2015, providing a comprehensive representation of real-world data. Hongyu Shao, as one of the authors, obtained appropriate authorization (certification number 15391-Nov2021) to access and analyze the SEER database. Notably, institutional review board approval was not necessary for this particular study.

2.2. Study Population

In this study, cancers were defined based on the SEER variable site recode and associated ICD-O-3 values, which enabled the determination of the sequence and temporal aspects of all cancer diagnosis. The methodology employed in this regard has been described comprehensively in previous literature [13]. Initially, a cohort of 198,972 patients were identified for inclusion in the study based on the following predefined criteria: (1) age greater than or equal to 18 years, (2) diagnosed within the time period spanning from 2010 to 2015, and (3) received pathological confirmation of both breast cancer and lung cancer. Patients with a prior diagnosis of a third type of cancer preceding the occurrence of lung cancer, as well as those with incomplete datasets, were excluded from the analysis.

2.3. Covariates

We used covariates which was based on the previous studies as follows: the year of diagnosis, age (≥18 years old), sex (female and male), race (American Indian/Alaska Native, Asian/Pacific Islander, Black, and White), marital status (married, single, divorced, separated, and widowed), clinical staging(CS) tumor size, laterality (left, right, and bilateral), radiation, chemotherapy, months from diagnosis to treatment and breast subtype (luminal-A, luminal-B, HER2 enriched, and triple-negative). We described breast cancer stage using AJCC, 7th edition (2010–2015) and combined summary stage (2004+).

2.4. Outcomes

We measured three outcomes which were OS and DSS, and metastasis (bone, brain, liver, and lung). OS was the interval in months between breast cancer diagnosis and death from any cause. DSS was the interval in months between breast cancer diagnosis and death from breast cancer. Patients were followed until death or December 31, 2015.

2.5. Statistical Analysis

Data were categorized into continuous and categorical variables. Continuous variables were further divided into two types based on the normality of their distribution. Normally distributed continuous variables were presented as the mean ± standard deviation and compared between groups using the Student’s t-test. Categorical variables were presented as percentages. Variables were compared using the Chi-square tests (categorical variables), one-way ANOVA (normal distribution), and Kruskal–Wallis (skewed distribution) test.

Logistic and multivariable Cox regression analyses were adopted to assess the ER, PR, HER2, and ER−/PR+ on the risk and outcomes in BC-LuC. An extended logistic model approach was used for different covariates adjusted models. Subgroup analyses were stratified by some relevant effect covariates. Survival curves were plotted by Kaplan–Meier and log-rank analyses. All the analyses were performed with the statistical software packages R 3.3.2 (https://www.R-project.org, The R Foundation) and free statistics software versions 1.5. A two-tailed test was performed and  < 0.05 was considered statistically significant.

3. Results

3.1. Baseline Characteristics

A total of 198.972 participants were included in the final analysis, as depicted in Figure 1. Baseline characteristics of these participants were presented in Table 1, stratified by the presence or absence of ER+,PR+,HER+ and ER−/PR+ status. On average, the selected participants had a mean age of 59.8 ± 13.1 years, with age categorized into three groups based on cutoffs of 40 and 60 years. Approximately, 99.3% of the participants were female. The prevalence of BC-LuC was comparably distributed across the four groups. Notably, patients aged below 60 years exhibited higher rates of HER2+ and ER−/PR+ subtypes, while the ER−/PR+ subtype was more prevalent among black individuals compared to other racial groups. The interval from diagnosis to treatment was categorized into three groups: ≤1 month, 1–3 months, and >3 months, with approximately 75.6% of the patients falling into the ≤1 month category, with approximately 75.6% of the patients falling into the ≤1 month category. Furthermore, an association was observed between larger tumor sizes (as per the CS) and increased occurrence of HER2+ and ER−/PR+ subtypes. Radiotherapy usage was comparable across all groups, although chemotherapy utilization was relatively higher among HER2+ and ER−/PR+ patients. Analysis of DSS and OS revealed a higher risk of mortality among individuals exhibiting ER−/PR+ status within the selected participant cohort.

3.2. Relationship between ER, PR, HER2, and ER−/PR+ and the Risk of BC-LuC

The results of our study are summarized in Table 2. In this investigation, we employed multivariate binary logistic regression to construct four models aimed at examining the independent effects of ER,PR,HER2 and ER−/PR+ status on the risk of BC-LuC. The odds ratios (ORs) along with their corresponding 95% confidence intervals (CIs) are presented in Table 2.

In the unadjusted model, the observed effect size can be interpreted as the change in the risk of BC-LuC associated with a one-unit decrease in the presence of ER+, PR+, HER2+, or ER−/PR+. Specifically, for HER2+ patients with BC-LuC, the unadjusted model indicated a 29% reduction in risk (OR = 0.71, 95% CI: 0.61–0.81,  < 0.001, Table 2, Figure 2).

Following adjustment for all covariates, a 19% decreased risk was noted among PR+ patients with BC-LuC (HR = 0.81, 95% CI: 0.71–0.93,  = 0.002), while HER2+ patients exhibited a 17% lower risk (HR = 0.83, 95% CI: 0.72–0.96,  = 0.012, Table 2).

3.3. Subgroup Analysis

In order to observe the trend of effect sizes in various variables, namely, CS tumor size, marital status, months from diagnosis to treatment, PR+ status, and HER2+ status, we employed these variables as stratification factors (Figure 2). We assessed for potential interactions within subgroups but did not find any statistically significant interactions, except for the subgroup of individuals with HER2+ status in relation to marital status (>0.05 for all other subgroups).

3.4. The Results of Outcomes of ER, PR, HER2, and ER−/PR+ in BC-LuC and BC

In this study, we conducted an analysis to investigate the association between ER,PR,HER2, ER−/PR+ status, and various clinical outcomes including DSS,OS, and the occurrence of metastases in bone, brain, liver, and lung. Utilizing multivariable Cox models, we identified significant associations between the hazard ratios (HRs) of ER+ and + in relation to BC-LuC for both DSS and OS (Supplementary Table 1 and Figures 3 and 4 in the attachments). Furthermore, we observed that BC patients in the ER+, PR+, and HER2+ subgroups exhibited longer DSS and OS (Supplementary Table 2 and Supplementary Figures 1 and 2 in the attachments). In contrast, individuals with ER−/PR+ status in BC demonstrated lower disease-specific DSS and OS (Supplementary Table 2 and Supplementary Figures 1 and 2 in the attachments).

Due to the limited number of metastatic patients specifically within the BC-LuC subgroup, we conducted our analysis on the entire study population, which included both individuals with BCand BC-LuC. Through the utilization of multivariate binary logistic regression and subgroup analysis (Supplemetary Table 3 and 4 in the attachments), we explored the relationship between ER,PR,HER2 and ER−/PR+ status, and the occurrence of metastases in the bone, brain, liver, and lung.

Regarding bone metastasis, our findings revealed that HER2+ and ER−/PR+ status were associated with a reduced risk of metastasis, whereas ER+ and PR+ status exhibited an opposite effect. In the case of lung metastasis, ER+ status was found to be associated with a reduced risk, while no significant relationships were observed for the other groups. In terms of liver metastasis, both ER2+ and PR+ status were associated with a reduced risk, whereas HER2+ and ER−/PR+ status manifested an opposite effect. Lastly, for brain metastasis, our analysis demonstrated that ER+ and PR+ status were associated with a reduced risk, whereas no significant associations were observed for the other groups.

4. Discussion

Our findings revealed a significant positive association between PR+ and HER2+ status and a lower risk of BC-LuC when adjusting for other covariates. Furthermore, the results of our study identified stronger associations in CS tumor size, months from diagnosis to treatment, PR+ status, and HER2+ status. The use of subgroup analysis allowed for a more comprehensive understanding of the trend towards a lower risk of BC-LuC in different populations. Moreover, our study demonstrated that ER+ and PR+ status were associated with increased DSS and OS in BC-LuC. Similarly, ER+, PR+, and HER2+ status were associated with improved DSS and OS in BC cases. Conversely, an opposite effect was observed in BC cases with ER−/PR+ status in relation to DSS and OS. Moreover, our study demonstrated that ER+ and PR+ status were associated with increased DSS and OS in BC-LuC. Similarly, ER+, PR+, and HER2+ status were associated with improved DSS and OS in BC cases. Conversely, an opposite effect was observed in BC cases with ER−/PR+ status in relation to DSS and OS. Moreover, our study demonstrated that ER+ and PR+ status were associated with increased DSS and OS in BC-LuC. Similarly, ER+, PR+, and HER2+ status were associated with improved DSS and OS in BC cases. Conversely, an opposite effect was observed in BC cases with ER−/PR+ status in relation to DSS and OS.

As per the meta-analysis conducted by Wang, only one study investigated the relationship between the expression statuses of ER, PR, and HER2 and the risk of developing lung cancer subsequently in female patients with breast cancer [14]. In a study involving 535,941 breast cancer patients [15], the results indicated that a positive estrogen receptor status (RR = 0.93 and  = 0.014) and positive progesterone receptor status (RR = 0.86 and  < 0.001) were associated with a lower risk of developing BC-LuC. Another study conducted by Schonfeld et al. [16] demonstrated that the absence of any negative breast cancer receptor marker increased the risk of BC-LuC. Additionally, Lin’s et al. study [17] reported that HER2 status had no significant effect on the development of BC-LuC. These conclusions align precisely with the findings of our own study. Several potential factors may account for the disparity in results observed between our study and previous research. First, variations in the research population could have contributed to these discrepancies. While some studies specifically focused on female patients, our analysis encompassed all patients, irrespective of gender. Second, unlike our study, the aforementioned research failed to consider the impact of ER, PR, and HER2 status on their interrelationships while adjusting for covariates. Furthermore, subgroup analysis was not conducted in the previous study. It is worth noting that the relationship between ER−/PR+ status and BC-LuC has not been definitively established in previous studies.

Our study sheds light on the potential impacts of ER, PR, and HER2 expression status on the overall and disease-specific survival of patients with BC-LuC. Specifically, we found that ER+ and PR+ statuses were associated with prolonged overall survival but decreased disease-specific survival in patients with BC-LuC. Furthermore, ER+, PR+, HER2+ statuses were also associated with prolonged OS but decreased DSS in BC patients. However, the results pertaining to HER2+ status were not entirely consistent. Previous studies [1822] have reported a similar association between ER+ and PR+ status and improved OS, DSS, and progression-free survival, for both male and lactating breast cancer. This aligns with our findings, although we did not observe any correlation between ER−/PR+ status and clinical outcomes in BC-LuC due to the limited sample size. Our study’s novel exploration of this breast cancer subtype warrants further investigation.

The present study involved a limited number of cases with metastasis in BC-LuC and aimed to analyze the role of ER-, -PR-, and HER2, as well as ER−/PR+ status, in both BC and BC-LuC overall. In line with Salvador’s review [23], which examined breast cancer patients with distant disease, particularly those with the ER+ subtype showed a higher prevalence of bone metastasis. To identify breast cancer patients with brain metastases at diagnosis, a retrospective analysis was conducted using the Surveillance, Epidemiology, and End Results (SEER) database of the National Cancer Institute, including 238,726 adult patients diagnosed with invasive breast cancer between 2010 and 2013, with recorded information on the presence or absence of brain metastases [24]. The findings from this study indicated that HR−/HER2+ and triple-negative subtypes were more prone to brain metastasis. Consistently, Leone’s study [25] also utilizing the SEER database, revealed that HR−/HER2+ subtype had higher odds of liver metastases, while triple-negative subtype had higher odds of brain metastases. These results corroborate our findings. Several studies [26, 27] have demonstrated that ER, PR, and HER2 status may change after breast cancer metastasis without significantly impacting progression-free survival. Given the close association between ER, PR, HER2 and metastasis, identifying indicators that are linked to an increased risk of metastasis is crucial. In this context, elucidating the relationship between ER−/PR+ status as an independent phenotype and metastasis holds particular significance for understanding the survival and burden of breast cancer patients.

The clinical significance of this study is twofold: first, it represents the pioneering observation of an independent association between ER−/PR+ status and the risk and outcomes in BC-LuC, specifically in the context of metastasis. To the best of our knowledge, this is the first study to report such an association. Second, the findings of this study hold considerable potential for guiding future research endeavors aimed at developing diagnostic or predictive models specifically tailored for BC-LuC. By shedding light on the relationship between ER−/PR+ status and BC-LuC, these findings serve as a valuable foundation for further investigations in this field.

This study is subject to several noteworthy limitations. First, the SEER database used in this study provided incomplete and limited information regarding smoking, alcohol consumption, and certain blood biochemical indicators associated with breast cancer prognosis. Consequently, the results of logistic and Cox analyses may have been affected. Second despite the large-scale nature of the patient cohort included in the SEER database, the number of eligible patients with metastasis for inclusion in this study remained relatively low, particularly within certain analytical categories. Therefore, it is imperative to validate our findings by including a larger number of patients from diverse geographical locations.

5. Conclusion

In conclusion, this study provides insights into the risk and outcomes of patients with BC-LuC. Our findings suggest that compared to patients with single breast cancer, the presence of PR+ and HER2+ subtypes may potentially be associated with a reduced risk of developing BC-LuC. Furthermore, we observed that ER+, PR+, HER2+, and ER−/PR+ status may have an impact on the outcomes of BC-LuC and BC patients. These observations highlight the need for further investigation into the association between molecular subtypes and BC-LuC. It is crucial to conduct additional clinical trials to confirm and validate these associations in order to enhance our understanding and clinical management of BC-LuC.

Data Availability

The data used to support the findings of this study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors thank all the staff members in our institution. This work was supported by Self-funded Projects in the Key Research and Development Plan of Xingtai City (grant number No. 2022ZC272).

Supplementary Materials

Supplementary Table 1: multivariable Cox regression analyses HR and 95% CI of the ER, PR, HER2, and ER−/PR+ associated with outcomes in BC-LuC. Supplementary Table 2: multivariable Cox regression analyses HR and 95% CI of the ER, PR, HER2, and ER−/PR+ associated with outcomes in BC. Supplementary Table 3: the association of OR and 95% CI between the ER, PR, HER2, and ER−/PR+ with metastasis in BC and BC-LuC. Supplementary Figure 1: Kaplan–Meier survival curves for DSS of ER, PR, and HER2 and ER−/PR+ in BC. Supplementary Figure 2: Kaplan–Meier survival curves for OS of ER, PR, and HER2 and ER−/PR+ in BC. (Supplementary Materials)