- 1Department of Neurosurgery, The First Affiliated Hospital of Jinan University, Guangzhou, China
- 2Department of Clinical Research, The First Affiliated Hospital of Jinan University, Guangzhou, China
- 3Department of Pathology, Shanxi Provincial People's Hospital, Taiyuan, China
Background: The extent of the relationship between age and the presence of breast cancer synchronous brain metastases (BCSBMs) and mortality has not yet been well-identified or sufficiently quantified. We aimed to examine the association of age with the presence of BCSBMs and all-cause and cancer-specific mortality outcomes using the SEER database.
Methods: Age-associated risk of the presence and survival of BCSBMs were evaluated on a continuous scale (restricted cubic spline, RCS) with logistic or Cox regression models. The main endpoints were the presence of BCSBMs and all-cause mortality or cancer-specific mortality. Cox proportional hazards regression and competing risk models were used in survival analysis.
Results: Among 374,132 adult breast cancer patients, 1,441 (0.38%) had BMs. The presence of BCSBMs displayed a U-shaped relationship with age, with the highest point of the curve occurring at the age of 62. In both the younger (age ≤ 61) and older (age ≥ 62) groups, the observed curve showed a nearly linear relationship between age and the presence of BCSBMs. The relationship between age and all-cause mortality (ASM) and cancer-specific mortality (CSM) was linear. Older age at diagnosis was associated with a higher risk of ASM (HR 1.019, 95% CI: 1.013–1.024, p < 0.001) and CSM (HR 1.016, 95% CI: 1.010–1.023, p < 0.001) in multivariable Cox models. Age (sHR 1.007, 95% CI 1–1.013, p = 0.049) was substantially related to a significantly increased risk of CSM in competing risk models.
Conclusion: Age had a non-linear U-shaped relationship with the presence of BCSBMs and a linear relationship with BCSBMs mortality.
Introduction
The second most frequent solid tumor that can metastasize to the central nervous system is breast cancer (BC) (1, 2). Brain metastases (BMs) are expected to occur in 30–50% of people with metastatic BC (3–5). The brain microenvironment is vastly different from that of extracranial lesions, with its distinct cell types, architectural features, metabolic restrictions, and immunological milieu, which influence the metastatic process and treatment responses (6, 7). According to many important studies, HER2-positive disease, the existence of more than two metastatic sites at BC diagnosis, HR negative, and a more advanced stage of the original tumor were all linked to a greater risk of breast cancer BMs (4, 8, 9).
Despite this compelling evidence, there is less evidence about the associations between age and the presence of breast cancer synchronous brain metastases (BCSBMs). Studies on this topic found different results, with some studies showing that breast cancer BMs occurred more frequently among younger women (10, 11), one study considering that age had no impact on the presence of breast cancer BMs (12), some studies suggesting that patients with older age had greater odds of having breast cancer BMs (5, 9), and some studies regarding advanced age as a risk factor for the presence of BMs (13, 14). There are currently no reliable population-based estimations of the relationship between age and the occurrence of breast cancer BMs.
Prior studies establish that BMs confer a life-threatening prognosis for female BC patients (15–18). Diagnosed with BMs represents an independent risk factor for shorter survival time in a large cohort retrospective study, and it has been estimated that there is a 58% rise in the risk of death from all causes (15). Breast cancer patients previously had a median survival period of 3–6 months from the time their BMs were discovered (17). Another prospective study found that women with brain metastases have a median survival of 26.3 months, compared to 44.6 months for women who do not suffer from brain metastases (19). Furthermore, BMs have been shown to be a reliable predictor of bone metastases in patients with infiltrating duct carcinoma of the breast, which is associated with a worsening prognosis (20).
Age is a significant determinant in cancer incidence and survival, and it is also a considerable factor in BM survival (13, 14, 21). Studies have shown that increasing age has generally been associated with poorer survival of BMs (12–14). There were, however, some inconsistent results (11, 22). Furthermore, there is no credible data on the impact of age on the prognosis of breast cancer patients with BMs.
We aimed to examine in detail the association of age with the presence of synchronous BMs and all-cause and cancer-specific mortality outcomes at diagnosis of BC using a large, multicenter, contemporary, population-based cohort in the United States from the Surveillance, Epidemiology, and End Results (SEER) database.
Methods
Data sources and study population
We conducted a cohort study using data from the SEER program, which contains demographic, illness, and treatment-related information for 34.6 percent of cancer patients in the United States at the time of primary malignancy diagnosis (23). For patients diagnosed between 2010 and 2016, information about the presence or absence of brain metastases at the time of the initial systemic malignancy diagnosis was available. We examined 388,413 individuals aged 18 and above who were diagnosed with primary, invasive breast cancer between January 1, 2011, and December 31, 2016, in the SEER database.
We excluded individuals who were male (3 016 individuals), lacked data for education and household income (64 individuals), were diagnosed with carcinoma in situ (561 individuals), lacked information on bone metastases (9 380 individuals), liver metastases (487 individuals), lung metastases (443 individuals), and brain metastases (214 individuals). We also excluded cases if survival time was unknown (26 individuals). The final analytical cohort for the association of age and the presence of BCSBM included 374 132 participants. Of these, a total of 1,441 individuals were identified as having brain metastases. In addition, we excluded 6 participants whose diagnosis was based on an autopsy or death certificate, as well as one participant of undetermined race, leaving 1,434 participants eligible for survival analysis. Our institutional review board granted an exemption for this study since it is a secondary analysis of existing data (SEER).
Covariates
Demographic variables included patients' age, sex, median household income, high school education percentage, year of diagnosis, race/ethnicity, registry region, marital status, and insurance status were reported by the SEER program. Clinical covariates such as tumor site, histological subtypes, T-stage, lymph nodal positive rate (LNPRate), and metastases in other organs (bone, liver, lung, and brain) were also included in the research. All the above variables were treated as factors, with the exception of age, household income, and education proportion, which were treated as continuous variables and given as mean and standard deviation (SD). The presence or absence of other organ metastases was confirmed before the start of treatment.
Endpoints
In the logistic regression models, the primary endpoint was the presence of BMs at diagnosis. In Cox proportional hazard models, the primary endpoints were all-cause mortality (ASM) and cancer-specific mortality (CSM) based on the International Classification of Diseases, 10th revision (ICD-10) code recorded as the underlying cause of death. In the competing risk models that were analyzed using proportional subdistribution hazards models, the primary endpoint was cancer-specific mortality, while other causes of mortality were the competing risk (24, 25). The months to the event were calculated from the time of diagnosis to the end of follow-up or death.
Statistical analysis
The presence of BMs
Independent factors in demographic variables and clinical covariates were used to determine whether independent factors were associated with the presence of BMs at diagnosis. The associations between age and the presence of BMs at diagnosis were evaluated on a continuous scale with restricted cubic spline curves (RCSs) based on logistic regression models with 4 knots at the 5th, 45th, 65th, and 90th percentiles of age (26). The spline model was adjusted for variables that were found to have significance in univariable logistic analysis (p ≤ 0.05). Then, sensitivity analyses were carried out to see if the findings were reliable.
Based on the cut-off value from the result of RCS, we divided the cohort into two age groups. The mean and SD were calculated for continuous variables, and the proportion was calculated for categorical variables in each age group. The t-test or chi-square test was used to calculate statistical differences for continuous and categorical variables. As the associations between age and the presence of BMs were approximately linear below and above the cut-off value, we additionally used multivariable logistic regression models to calculate the odds ratio (OR) and 95% confidence interval (CI).
The survival of BCSBMs
The hazard ratios of mortality were calculated using univariable and multivariable Cox proportional hazards regression models adjusted for possible confounders (27). We used restricted cubic spline models fitted to Cox proportional hazards models with 4 knots at the 5th, 45th, 65th, and 90th percentiles of age (26). ASM and CSM spline models were further adjusted for significant variables in ASM and CSM univariable Cox regression models, respectively.
A competing mortality risk regression analysis on cumulative incidence functions was conducted using Fine and Gray models to better estimate breast CSM and better account for the high rate of competing events (28, 29). The researchers calculated unadjusted and adjusted subdistribution hazard ratios (sHR) with 95% CI. The Cumulative Incidence Function (CIF) allows for estimating the incidence of CSM while accounting for competing risk.
Statistical analyses were conducted using R language program version 4.0.3 released on 10-10-2020 and STATA software version 14 (StataCorp). The 2-tailed α values of <0.05 were considered statistically significant.
Results
Of data from 374 132 adult patients from the SEER database in the United States, Group A comprised 182,033 (48.6%) patients under the age of 61 (≤ 61), with a mean (SD) age of 50.49 (7.77), whereas Group B (age ≥ 62 years old) contained 192,099 (51.35 %) patients with a mean (SD) age of 72.60 (7.96). Among the entire cohort, 1,441 patients were diagnosed with BMs, accounting for 0.38% of the entire study population. Groups A and B had 775 (53.78%) and 666 (46.22%) BCSBMs, respectively, with the incidence proportion of BMs in Groups A and B being 0.43 percent and 0.35 percent. Tables 1, 3, 4 show the baseline characteristics of the entire cohort and the BMs cohort.
Age and the presence of BMs
Figure 1 depicts the full analytical process. In univariable logistic regression, age (increased by per 1, odds ratio [OR], 0.99, 95%CI: 0.99–1, p < 0.001), non-Hispanic Black (NHB) (vs. non-Hispanic White [NHW], p < 0.001), Hispanic (vs. NHW, p = 0.017), West region (vs. Northeast, p = 0.001), other marital status (vs. married, p < 0.001), other insurance status (vs. insured, p < 0.001), tumor located in upper-inner or lower-inner or upper-outer or lower-outer or axillary tail of breast (vs. central, p < 0.001, p = 0.001, p < 0.001, p = 0.002, and p = 0.002, respectively), T-stage 2 or 3 or 4 (vs. 1, p < 0.001 for each one), LNPRate between 21 and 40% or between 41 and 60% or between 61 and 80% or between 81 and 100% (vs. 0–20%, p < 0.001 for each one), HR+/HER2+ or HR−/HER2+ or triple-negative subtype (vs. HR+/HER2− subtype, p < 0.001 for each one), and metastasized to bone or liver or lung at diagnosis (vs. Not to, p < 0.001 for each one) were related to significantly greater odds of having BMs at diagnosis. Median household income (increased by per $10,000 annual, p < 0.001) was at lower risk of presence of BMs at diagnosed. Restricted cubic spline revealed a U-shaped relationship between age and the presence of BMs after controlling for the aforementioned possible confounders (Figure 2). The risk of having BMs increased rapidly until approximately the age of 62, after which it began to decline rapidly (P for non-linearity < 0.001). The results for sensitivity analyses are shown in Supplementary Figure 1. In multivariable logistic regression models, each year of age conferred a 1% (95% CI: 1–1.02, p = 0.036, p for non-linearity = 0.163) increase in the OR to develop BMs in Group A (age ≤ 61) and a 4% (95% CI: 0.95–0.97, p < 0.001, p for non-linearity = 0.067) decrease in the OR to develop BMs in Group B (age ≥ 62). Detailed data is shown in Table 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1000415/fpubh-10-1000415-HTML-r1/image_m/fpubh-10-1000415-g002.jpg)
Figure 2. Association between age and the presence of BCSBMs using a restricted cubic spline regression model.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1000415/fpubh-10-1000415-HTML-r1/image_m/fpubh-10-1000415-t002.jpg)
Table 2. Univariable and multivariable logistic regression models for the presence of breast cancer synchronous brain metastases.
Age and mortality
A median (interquartile range) of 16 (6–32) months was the time period of mortality ascertainment, corresponding to 1,042 deaths from all causes and 807 deaths caused by cancer. RCS indicated an ascending all-cause mortality risk with increasing age (Figure 3A, p for non-linearity = 0.264), controlled for median household income, race, region, marital status, tumor site, LNPRate, histological subtype, and metastasized to bone, liver, and lung. There was also an escalating cancer-specific mortality risk with increasing age (Figure 3B, p for non-linearity = 0.473), controlled for median household income, race, marital status, insurance status, T-stage, histological subtype, and metastasized to the liver and lung. Additionally, in multivariable Cox proportional-hazard models, the adjusted HRs of age for mortality due to all-cause and cancer-specific were HR = 1.019 (95% CI: 1.013–1.024, p < 0.001) and HR = 1.016 (95% CI: 1.010–1.023, p < 0.001), respectively. The corresponding detailed results of the unadjusted and adjusted Cox models are shown in Tables 3, 4.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1000415/fpubh-10-1000415-HTML-r1/image_m/fpubh-10-1000415-g003.jpg)
Figure 3. Association between age and all-cause mortality (A) and cancer-specific mortality (B) using restricted cubic spline regression models.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1000415/fpubh-10-1000415-HTML-r1/image_m/fpubh-10-1000415-t003.jpg)
Table 3. Distribution of breast cancer synchronous brain metastases all-cause mortality and hazard risk associated with various prognostic factors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1000415/fpubh-10-1000415-HTML-r1/image_m/fpubh-10-1000415-t004.jpg)
Table 4. Distribution of breast cancer synchronous brain metastases cancer-specific mortality and hazard risk associated with various prognostic factors.
In addition, we evaluated CSM using a competing risk model. The CIF curves for the observed risk of CSM are shown in Supplementary Figure 2. Table 5 shows the results of multivariable competing-risk regression analyses predicting the time to CSM. Age (increased by per 1, sHR 1.007, 95% CI 1–1.013, p = 0.049), other marital status (vs. married, sHR 1.191, 95% CI 1.016–1.397, p < 0.001), other insurance status (vs. insured, sHR 1.221, 95% CI 1.034–1.442, p < 0.001), T-stage 3 (vs. 1, sHR 1.438, 95% CI 1.057–1.956, p = 0.021), T-stage 4 (vs. 1, sHR 1.612, 95% CI 1.242–2.091, p < 0.001), triple-negative subtype (vs. HR+/HER2− subtype, sHR 1.693, 95% CI 1.374–2.086, p < 0.001), and metastasized to liver (vs. not to, sHR 1.367, 95% CI 1.163–1.607, p < 0.001) were significantly associated with an increased risk for CSM.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1000415/fpubh-10-1000415-HTML-r1/image_m/fpubh-10-1000415-t005.jpg)
Table 5. Proportional subdistribution hazards models for breast cancer synchronous brain metastases.
Discussion
In a large population-based retrospective cohort study, we examined the complicated association between age and the presence and survival of BCSBMs. As far as we know, this is an epidemiologic study with great interest and novelty based on the SEER program that will generally offer a superior understanding of variation in the onset and prognosis of BCSBMs. The presence of BCSBMs was shown to have a U-shaped relationship with age, with the maximum HR occurring at the age of 62. The association between age and all-cause and cause-specific mortality, on the other hand, resembled linear behavior.
Comparison with other studies
In recent years, there has been an increasing amount of literature on the epidemiology of BMs at the diagnosis of systemic malignancy (12–14, 30). Using data from the SEER database from 2010 to 2013, Martin et al. (13) and Cagney et al. (14) found that 0.41 percent of adult breast cancer patients have synchronous BMs in the United States. In earlier research assessing the 2014–2016 SEER data, we found that 0.42 percent of midlife breast cancer diagnosed in the United States have synchronous BMs. The current study found similar results, and the incidence fraction of BCSBMs remained generally steady from 2011 to 2016.
In reviewing the above-mentioned literature, no solid evidence was found for the association between age and the presence of BMs. As mentioned in the research of Cagney et al. (14) among all patients with cancer, age 41–60 years (vs. age 18–40 y, OR 1.55, 95% CI: 1.39–1.71, p < 0.05) was associated with higher risk of having synchronous BMs, whereas age > 80 years (vs. age 18–40 y, OR 0.60, p < 0.05) had significantly lower odds and age 61–80 years had non-significantly greater odds (vs. age 18–40 y, OR 1.08, p > 0.05). According to another study concentrating on BCSBMs, the relationship between age and the occurrence of synchronous BMs was more elusive and harder to reconcile because OR values of age 41–60 years, 61–80 years, and > 80 years (vs. age 18–40 years) were all non-applicable (13). We reason that these seemingly misleading results might be due to the use of age as a categorical variable in the adjusted logistic regression model. Naturally, we assumed that treating age as a continuous variable would help to ameliorate the situation and lead to a positive conclusion. In this case, however, age (increased by per 1, OR 1, 95% CI: 0.99–1.00, p = 0.49) was not related to the presence of BMs (12). Another report discovered that age ≤ 40 years (vs. age > 40 years, HR 2.10, 95% CI 1.02–4.36) was associated with an increased risk of developing metachronous BMs in HER2-positive breast cancer (31). Additional research showed no information on the effect of age on the occurrence of BCSBMs (11, 32, 33).
A possible explanation for this might be that prior research relied on the essential assumption that the presence of BMs was related linearly to age (12–14, 33). For ordinal or continuous factors, the linearity assumption may be inappropriate, and more elaborate interactions may be necessary (34). Conversely, cubic splines are commonly used because they offer a lot of flexibility when it comes to fitting data, are visually smooth due to their continuous first and second derivatives, and have fewer fit constants than higher degree splines (35–39). In this study, we discovered that using a regression spline to solve such situations is a preferable option. Our findings on age were not in line with the previous findings, whereby we observed a U-shaped association (p for non-linearity <0.001) with the presence of BMs after accounting for potential confounders. Further research revealed that the relationships between age and the occurrence of BCSBMs are approximately linear in both the younger (p for non-linearity = 0.163) and older (p for non-linearity = 0.068) age groups. Consequently, we were able to convert a complex non-linear association into a linear one. Undoubtedly, therefore, age may represent a clinical marker for early identification of a population at high risk for having BCSBMs.
The National Comprehensive Cancer Network (NCCN), American Society of Clinical Oncology (ASCO), and European School of Oncology-Metastatic Breast Cancer guidelines for breast cancer do not recommend routine imaging assessment or continued imaging reassessment of BMs for breast cancer because the overall incidence of BMs is relatively low in the general BC population and there is no proven benefit from non-randomized retrospective studies (40–42). Imaging assessment of the brain is recommended only after the appearance of central nervous system symptoms (43, 44). Nonetheless, timely identification of BMs is critical for BC patients since it may allow better therapeutic responses than delays in diagnosis (45, 46). Research finds that the combination of early detection and advanced therapies (both local and systemic) ultimately leads to more favorable outcomes (47). Although routine screening of all BC patients is not justified, screening people who are at high risk might be beneficial.
For all malignancies, age is a well-validated prognosticator related to survival (21, 48, 49). And previous studies have found that being younger at the time of BCSBMs diagnosis is one of the characteristics that predicts a better outcome (12–14, 50, 51). To determine the association between age and the survival of BCSBMs, many researchers used Cox proportional hazards models (13, 22, 52). To our knowledge, none of them has carried out a comprehensive analysis to determine whether a Cox model would be appropriate in this circumstance. Our current study has significant strengths that compensate for the lack of information. Linear associations between age and ASM (p for non-linearity = 0.264) and CSM (p for non-linearity = 0.473) were discovered in this investigation. Those results suggest that the HR and 95%CI of age obtained by Cox models are reasonable.
However, one study found that in multivariate Cox regression analysis, younger age at first diagnosis of BMs in breast cancer patients was a predictor of shorter OS (52). We believe what the authors stated was a mistake after reviewing the original findings in the Supplementary material, which indicated that the risk of all-cause death rose with age (increased by per 1, HR = 1.02, 95% CI 1.01–1.03, p < 0.001). Another study found that in a univariable Cox model, age did not correlate with breast cancer BMs survival time (22). A possible explanation for this might be the potentially inappropriate classification of age groups. The research of Martin et al., on the other hand, performs better (13). Martin's study has strength in using Fine and Gray's competing risk regression models to analyze BCSBMs CSM (13). Yet, the study did not include CIF curves or sHR values for each variable. Here, we presented the sHR for the occurrence of CSM and gave plots of all cumulative incidences for categorical variables, as indicated by a prior study (24). Mortality risk tends to increase with age in all regression models, and age-related increases in the risk of CSM are substantially overestimated by the standard Cox model (53).
Age is commonly utilized as a covariate in investigations to build prognostic models for predicting the survival of BC patients with BMs, such as a recursive partitioning analysis (RPA) strategy (54) and graded prognostic assessment (GPA) (55, 56). In the era of individualized therapy, the accurate prediction of BC patients with BMs is critical for optimizing care (18, 57). To some extent, the current study is significant in determining the ability of prognostic tools in future research and an improved prognosis for those individuals.
Strengths and limitations
The discovery of a U-shaped association between age and the presence of cancer, as well as approximate linear behaviors between age and ASM and CSM, in a large, nationally representative sample of US general cancer patients from the SEER database with rigorous capture of death events, which was undetected by previous excellent work, is a major strength of our study. At the same time, the results of our competing risk model help to compensate for the lack of survival analyses.
Our research, however, has several limitations. First, in the current research, only the presence or absence of BMs at the diagnosis of the study population was provided. Data on whether metastases develop in the brain in the subsequent course of the disease is not available at this time for the SEER program. Our study did not include the data of some patients who later acquired brain metastases, which might have influenced the accuracy of the results (58). Second, the SEER program is a population-based study being carried out primarily in the United States, which has concerns for generalizability outside of the United States (59, 60). Third, we did not report sHR for the competing event, which may have resulted in a bias toward a better understanding of BCSBMs survival (61, 62).
Conclusion
In conclusion, utilizing a nationally representative database from the United States, this study discovered that age had a non-linear U-shaped relationship with the presence of BCSBMs and a linear relationship with BCSBMs mortality. This article lays the groundwork for future studies. And a better understanding of the complex association might aid in developing age-appropriate public health guidelines.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found at: https://seer.cancer.gov.
Ethics statement
This study was reviewed and approved by the First Affiliated Hospital of Jinan University. Written informed consent for participation was not required for this study in accordance with the national legislation and institutional requirements.
Author contributions
WC had the initial idea, analyzed the data, and wrote the paper. WC, YW, XW, and JL contributed to study design and commenting on drafts and revisions. All authors contributed to the article and approved the submitted version.
Acknowledgments
We heartily thank all of the employees and scientists at each of the SEER registry locations.
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.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2022.1000415/full#supplementary-material
References
1. Tabouret E, Chinot O, Metellus P, Tallet A, Viens P, Goncalves AJAr. Recent trends in epidemiology of brain metastases: an overview. Anticancer Res. (2012) 32:4655–62.
2. Pellerino A, Internò V, Mo F, Franchino F, Soffietti R, Rudà R. Management of brain and leptomeningeal metastases from breast cancer Int J Mol Sci. (2020) 21:8534. doi: 10.3390/ijms21228534
3. Brosnan EM, Anders CKJAotm. Understanding patterns of brain metastasis in breast cancer and designing rational therapeutic strategies. Ann Transl Med. (2018) 6:163. doi: 10.21037/atm.2018.04.35
4. Pasquier D, Darlix A, Louvel G, Fraisse J, Jacot W, Brain E, et al. Treatment and outcomes in patients with central nervous system metastases from breast cancer in the real-life ESME MBC cohort. Eur J Cancer. (2020) 125:22–30. doi: 10.1016/j.ejca.2019.11.001
5. Garcia-Alvarez A, Papakonstantinou A, Oliveira MJC. Brain metastases in HER2-positive breast cancer: current and novel treatment strategies. Cancers. (2021) 13:2927. doi: 10.3390/cancers13122927
6. Boire A, Brastianos PK, Garzia L, Valiente M. Brain metastasis. Nat Rev Cancer. (2020) 20:4–11. doi: 10.1038/s41568-019-0220-y
7. Bailleux C, Eberst L, Bachelot T. Treatment strategies for breast cancer brain metastases. Br J Cancer. (2021) 124:142–55. doi: 10.1038/s41416-020-01175-y
8. Pestalozzi BC, Zahrieh D, Price K, Holmberg S, Lindtner J, Collins J, et al. Identifying breast cancer patients at risk for central nervous system (CNS) metastases in trials of the international breast cancer study group (IBCSG). Ann Oncol. (2006) 17:935–44. doi: 10.1093/annonc/mdl064
9. Sperduto PW, Mesko S, Cagney D, Nesbit E, Chan J, Lee J, et al. Tumor subtype and other prognostic factors in breast cancer patients with brain metastases: the updated graded prognostic assessment (Breast-GPA). J Clin Oncol. (2019) 37(15_Suppl) 1079. doi: 10.1200/JCO.2019.37.15_suppl.1079
10. Yeh RH, Yu JC, Chu CH, Ho CL, Kao HW, Liao GS, et al. Distinct MR imaging features of triple-negative breast cancer with brain metastasis. J Neuroimaging. (2015) 25:474–81. doi: 10.1111/jon.12149
11. Rostami R, Mittal S, Rostami P, Tavassoli F, Jabbari BJ. Brain metastasis in breast cancer: a comprehensive literature review. J Neurooncol. (2016) 127:407–14. doi: 10.1007/s11060-016-2075-3
12. Che W, Wang Y, Wang X, Lyu J. Midlife brain metastases in the United States: is male at risk? Cancer Med. (2022) 11:1202–16. doi: 10.1002/cam4.4499
13. Martin AM, Cagney DN, Catalano PJ, Warren LE, Bellon JR, Punglia RS, et al. Brain metastases in newly diagnosed breast cancer: a population-based study. JAMA Oncol. (2017) 3:1069–77. doi: 10.1001/jamaoncol.2017.0001
14. Cagney DN, Martin AM, Catalano PJ, Redig AJ, Lin NU, Lee EQ, et al. Incidence and prognosis of patients with brain metastases at diagnosis of systemic malignancy: a population-based study. Neurooncology. (2017) 19:1511–21. doi: 10.1093/neuonc/nox077
15. Jung SY, Rosenzweig M, Sereika SM, Linkov F, Brufsky A, Weissfeld JL. Factors associated with mortality after breast cancer metastasis. Cancer Causes Control. (2012) 23:103–12. doi: 10.1007/s10552-011-9859-8
16. Gori S, Rimondini S, De Angelis V, Colozza M, Bisagni G, Moretti G, et al. Central nervous system metastases in HER-2–positive metastatic breast cancer patients treated with trastuzumab: incidence, survival, and risk factors. Oncologist. (2007) 12:766–73. doi: 10.1634/theoncologist.12-7-766
17. Mahmoud-Ahmed AS, Suh JH, Lee S-Y, Crownover RL, Barnett G. Results of whole brain radiotherapy in patients with brain metastases from breast cancer: a retrospective study. Int J Radiat Oncol Biol Phys. (2002) 54:810–7. doi: 10.1016/S0360-3016(02)02967-X
18. Marko NF, Xu Z, Gao T, Kattan MW, Weil RJ. Predicting survival in women with breast cancer and brain metastasis: a nomogram outperforms current survival prediction models. Cancer. (2012) 118:3749–57. doi: 10.1002/cncr.26716
19. Brufsky AM, Mayer M, Rugo HS, Kaufman PA, Tan-Chiu E, Tripathy D, et al. Central nervous system metastases in patients with HER2-positive metastatic breast cancer: incidence, treatment, and survival in patients from registHER. Clin Cancer Res. (2011) 17:4834–43. doi: 10.1158/1078-0432.CCR-10-2962
20. Huang Z, Hu C, Liu K, Yuan L, Li Y, Zhao C, et al. Risk factors, prognostic factors, and nomograms for bone metastasis in patients with newly diagnosed infiltrating duct carcinoma of the breast: a population-based study. BMC Cancer. (2020) 20:1–17. doi: 10.1186/s12885-020-07635-1
21. Wang G-M, Cioffi G, Patil N, Waite KA, Lanese R, Ostrom QT, et al. Importance of the intersection of age and sex to understand variation in incidence and survival for primary malignant gliomas. Neurooncology. (2022) 24:302–10. doi: 10.1093/neuonc/noab199
22. Niu L, Lv H, Zhang M, Zeng H, Wang L, Cui S, et al. Clinical diagnosis and treatment of breast cancer with brain metastases and establishment of a prognostic model: a 10-year, single-center, real-world study of 559 cases. Ann Transl Med. (2021) 9:1331. doi: 10.21037/atm-21-3734
23. Yang J, Li Y, Liu Q, Li L, Feng A, Wang T, et al. Brief introduction of medical database and data mining technology in big data era. J Evid Based Med. (2020) 13:57–69. doi: 10.1111/jebm.12373
24. Latouche A, Allignol A, Beyersmann J, Labopin M. Fine J. A competing risks analysis should report results on all cause-specific hazards and cumulative incidence functions. J Clin Epidemiol. (2013) 66:648–53. doi: 10.1016/j.jclinepi.2012.09.017
25. Austin PC, Lee DS, Fine JPJC. Introduction to the analysis of survival data in the presence of competing risks. Circulation. (2016) 133:601–9. doi: 10.1161/CIRCULATIONAHA.115.017719
26. Frank E. Harrell J. Regression Modeling Strategies: With Applications to Linear Models, Logistic and Ordinal Regression, and Survival Analysis. New York, NY: Springer (2001).
27. Wu W-T, Li Y-J, Feng A-Z, Li L, Huang T, Xu A-D, et al. Data mining in clinical big data: the frequently used databases, steps, and methodological models. Mil Med Res. (2021) 8:1–12. doi: 10.1186/s40779-021-00338-z
28. Scrucca L, Santucci A, Aversa F. Regression modeling of competing risk using R: an in depth guide for clinicians. Bone Marrow Transplant. (2010) 45:1388–95. doi: 10.1038/bmt.2009.359
29. Massa ST, Osazuwa-Peters N, Christopher KM, Arnold LD, Schootman M, Walker RJ, et al. Competing causes of death in the head and neck cancer population. Oral Oncology. (2017) 65:8–15. doi: 10.1016/j.oraloncology.2016.12.006
30. Lamba N, Wen PY, Aizer A. Epidemiology of brain metastases and leptomeningeal disease. Neuro Oncol. (2021) 23:1447–56. doi: 10.1093/neuonc/noab101
31. Maurer C, Tulpin L, Moreau M, Dumitrescu C, de Azambuja E, Paesmans M, et al. Risk factors for the development of brain metastases in patients with HER2-positive breast cancer. ESMO Open. (2018) 3:e000440. doi: 10.1136/esmoopen-2018-000440
32. Laakmann E, Witzel I, Neunhöffer T, Weide R, Schmidt M, Park-Simon T-W, et al. Characteristics and clinical outcome of breast cancer patients with asymptomatic brain metastases. Cancers. (2020) 12:2787. doi: 10.3390/cancers12102787
33. Ascha MS, Ostrom QT, Wright J, Kumthekar P, Bordeaux JS, Sloan AE, et al. Lifetime occurrence of brain metastases arising from lung, breast, and skin cancers in the elderly: a SEER-Medicare study. Cancer Epidemiol Biomarkers Prev. (2019) 28:917–25. doi: 10.1158/1055-9965.EPI-18-1116
34. Durrleman S, Simon RJ. Flexible regression models with cubic splines. Stat Med. (1989) 8:551–61. doi: 10.1002/sim.4780080504
35. Bhaskaran K, dos-Santos-Silva I, Leon DA, Douglas IJ, Smeeth L. Association of BMI with overall and cause-specific mortality: a population-based cohort study of 3·6 million adults in the UK. Lancet Diabetes Endocrinol. (2018) 6:944–53. doi: 10.1016/S2213-8587(18)30288-2
36. Inoue K, Ritz B, Brent GA, Ebrahimi R, Rhee CM. Association of subclinical hypothyroidism and cardiovascular disease with mortality. JAMA Netw Open. (2020) 3:e1920745. doi: 10.1001/jamanetworkopen.2019.20745
37. Ho AS, Luu M, Zalt C, Morris LG, Chen I, Melany M, et al. Mortality risk of nonoperative papillary thyroid carcinoma: a corollary for active surveillance. Thyroid. (2019) 29:1409–17. doi: 10.1089/thy.2019.0060
38. Johannesen CDL, Langsted A, Mortensen MB, Nordestgaard B. Association between low density lipoprotein and all cause and cause specific mortality in Denmark: prospective cohort study. BMJ. (2020) 371:m4266. doi: 10.1136/bmj.m4266
39. Lee DH, Keum N, Hu FB, Orav EJ, Rimm EB, Willett WC, et al. Predicted lean body mass, fat mass, and all cause and cause specific mortality in men: prospective US cohort study. BMJ. (2018) 362:k2575. doi: 10.1136/bmj.k2575
40. Gradishar WJ, Anderson BO, Abraham J, Aft R, Agnese D, Allison KH, et al. Breast cancer, version 3.2020, NCCN clinical practice guidelines in oncology. J Natl Compr Canc Netw. (2020) 18:452–78. doi: 10.6004/jnccn.2020.0016
41. Ramakrishna N, Temin S, Chandarlapaty S, Crews JR, Davidson NE, Esteva FJ, et al. Recommendations on disease management for patients with advanced human epidermal growth factor receptor 2–positive breast cancer and brain metastases: American society of clinical oncology clinical practice guideline. J Clin Oncol. (2014) 32:2100. doi: 10.1200/JCO.2013.54.0955
42. Cardoso F, Costa A, Norton L, Senkus E, Aapro M, Andre F, et al. 2nd international consensus guidelines for advanced breast cancer (ABC2). Breast. (2014) 23:489–502. doi: 10.1016/j.breast.2014.08.009
43. Canney P, Murray E, Dixon-Hughes J, Lewsley L-A, Paul J. A prospective randomised phase III clinical trial testing the role of prophylactic cranial radiotherapy in patients treated with trastuzumab for metastatic breast cancer—Anglo Celtic VII. Clin Oncol. (2015) 27:460–4. doi: 10.1016/j.clon.2015.04.033
44. Cagney DN, Martin AM, Catalano PJ, Brown PD, Alexander BM, Lin NU. Implications of screening for brain metastases in patients with breast cancer and non–small cell lung cancer. JAMA Oncol. (2018) 4:1001–3. doi: 10.1001/jamaoncol.2018.0813
45. Mahner S, Schirrmacher S, Brenner W, Jenicke L, Habermann C, Avril N. Comparison between positron emission tomography using 2-[fluorine-18] fluoro-2-deoxy-D-glucose, conventional imaging and computed tomography for staging of breast cancer. Ann Oncol. (2008) 19:1249–54. doi: 10.1093/annonc/mdn057
46. Hodgson NC, Gulenchyn K. Is there a role for positron emission tomography in breast cancer staging? J Clin Oncol. (2008) 26:712–20. doi: 10.1200/JCO.2007.13.8412
47. Soerjomataram I, Louwman MW, Ribot JG, Roukema JA. An overview of prognostic factors for long-term survivors of breast cancer. Breast Cancer Res Treat. (2008) 107:309–30. doi: 10.1007/s10549-007-9556-1
48. Ostrom QT, Gittleman H, Fulop J, Liu M, Blanda R, Kromer C, et al. CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2008-2012. Neuro-oncology. (2015) 17(Suppl. 4):iv1–62. doi: 10.1093/neuonc/nov189
49. Ostrom QT, Patil N, Cioffi G, Waite K, Kruchko C, Barnholtz-Sloan J. CBTRUS statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2013–2017. Neuro-oncology. (2020) 22(Supplement_1):iv1–96. doi: 10.1093/neuonc/noaa200
50. Lamba N, Kearney RB, Catalano PJ, Hassett MJ, Wen PY, Haas-Kogan DA. Population-based estimates of survival among elderly patients with brain metastases. Neuro-oncology. (2021) 23:661–76. doi: 10.1093/neuonc/noaa233
51. Ho V, Gijtenbeek J, Brandsma D, Beerepoot L, Sonke GJ. Survival of breast cancer patients with synchronous or metachronous central nervous system metastases. Eur J Cancer. (2015) 51:2508–16. doi: 10.1016/j.ejca.2015.07.040
52. Witzel I, Laakmann E, Weide R, Neunhoeffer T, Park-Simon T-J, Schmidt M, et al. Treatment and outcomes of patients in the brain metastases in breast cancer network registry. Eur J Cancer. (2018) 102:1–9. doi: 10.1016/j.ejca.2018.07.004
53. Wolbers M, Koller MT, Witteman JC, Steyerberg EWJE. Prognostic models with competing risks: methods and application to coronary risk prediction. Epidemiology. (2009) 20:555–61. doi: 10.1097/EDE.0b013e3181a39056
54. Gaspar L, Scott C, Rotman M, Asbell S, Phillips T, Wasserman T, et al. Recursive partitioning analysis (RPA) of prognostic factors in three radiation therapy oncology group (RTOG) brain metastases trials. Int J Radiat Oncol Biol Phys. (1997) 37:745–51. doi: 10.1016/S0360-3016(96)00619-0
55. Sperduto PW, Berkey B, Gaspar LE, Mehta M, Curran W. A new prognostic index and comparison to three other indices for patients with brain metastases: an analysis of 1,960 patients in the RTOG database. Int J Radiat Oncol Biol Phys. (2008) 70:510–4. doi: 10.1016/j.ijrobp.2007.06.074
56. Sperduto CM, Watanabe Y, Mullan J, Hood T, Dyste G, Watts C, et al. A validation study of a new prognostic index for patients with brain metastases: the graded prognostic assessment. J Neurosurg. (2008) 109(Supplement):87–9. doi: 10.3171/JNS/2008/109/12/S14
57. Huang Z, Sun B, Wu S, Meng X, Cong Y, Shen G, et al. nomogram for predicting survival in patients with breast cancer brain metastasis. Oncol Lett. (2018) 15:7090–6. doi: 10.3892/ol.2018.8259
58. Pestalozzi BC, Holmes E, de Azambuja E, Metzger-Filho O, Hogge L, Scullion M, et al. CNS relapses in patients with HER2-positive early breast cancer who have and have not received adjuvant trastuzumab: a retrospective substudy of the HERA trial (BIG 1-01). Lancet Oncol. (2013) 14:244–8. doi: 10.1016/S1470-2045(13)70017-2
59. Du XL, Lairson DR, Begley CE, Fang S. Temporal and geographic variation in the use of hematopoietic growth factors in older women receiving breast cancer chemotherapy: findings from a large population-based cohort. J Clin Oncol. (2005) 23:8620. doi: 10.1200/JCO.2005.02.6252
60. Warren JL, Butler EN, Stevens J, Lathan CS, Noone A-M, Ward KC. Receipt of chemotherapy among Medicare patients with cancer by type of supplemental insurance. J Clin Oncol. (2015) 33:312. doi: 10.1200/JCO.2014.55.3107
61. Chen D, Liu Z, Liu W, Fu M, Jiang W, Xu S, et al. Predicting postoperative peritoneal metastasis in gastric cancer with serosal invasion using a collagen nomogram. Nat Commun. (2021) 12:1–11. doi: 10.1038/s41467-020-20429-0
Keywords: breast cancer, brain metastases (BMs), restricted cubic spline (RCS), linearity, non-linearity, prognosis, SEER
Citation: Che W, Wang Y, Wang X and Lyu J (2022) Association between age and the presence and mortality of breast cancer synchronous brain metastases in the United States: A neglected SEER analysis. Front. Public Health 10:1000415. doi: 10.3389/fpubh.2022.1000415
Received: 22 July 2022; Accepted: 24 August 2022;
Published: 23 September 2022.
Edited by:
Liyuan Tao, Peking University Third Hospital, ChinaReviewed by:
Chuan Sun, Zhejiang Hospital, ChinaXiaoran Yin, Second Affiliated Hospital of Xi'an Jiaotong University, China
Copyright © 2022 Che, Wang, Wang and Lyu. 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: Jun Lyu, lyujun2020@jnu.edu.cn; Xiangyu Wang, wang_xy123@126.com