Next Article in Journal
Macrophage-like Cells Are Increased in Retinal Vein Occlusion and Correlate with More Intravitreal Injections and Worse Visual Acuity Outcomes
Previous Article in Journal
Assessment of Two Commercial Comprehensive Gene Panels for Personalized Cancer Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Web-Based Dynamic Nomograms for Predicting Overall Survival and Cancer-Specific Survival in Breast Cancer Patients with Lung Metastases

1
Department of General Surgery, Xiangya Hospital, Central South University, Changsha 410008, China
2
National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha 410008, China
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2023, 13(1), 43; https://doi.org/10.3390/jpm13010043
Submission received: 20 October 2022 / Revised: 12 December 2022 / Accepted: 21 December 2022 / Published: 26 December 2022

Abstract

:
Background: 60–70% of patients who die from breast cancer have lung metastases. However, there is a lack of readily available tools for accurate risk stratification in patients with breast cancer lung metastases (BCLM). Therefore, a web-based dynamic nomogram was developed for BCLM to quickly, accurately, and intuitively assess overall and cancer-specific survival rates. Methods: Patients diagnosed with BCLM between 2004 and 2016 were extracted from the Surveillance, Epidemiology, and Final Results (SEER) database. After excluding incomplete data, all patients were randomly assigned to training and validation cohorts (2:1). Patients’ basic clinical information, detailed pathological staging and treatment information, and sociological information were included in further analysis. Nomograms were constructed following the evaluations of the Cox regression model and verified using the concordance index (C-index), calibration curves, time-dependent receiver operating characteristic (ROC) curves, and decision curve analysis (DCA). Web-based dynamic nomograms were published online. Results: 3916 breast cancer patients with lung metastases were identified from the SEER database. Based on multivariate Cox regression analysis, overall survival (OS) and cancer-specific survival (CSS) are significantly correlated with 13 variables: age, marital status, race, grade, T stage, surgery, chemotherapy, bone metastatic, brain metastatic, liver metastatic, estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor-2 (HER2). These are included in the construction of the nomogram of OS and CSS. The time-dependent receiver operating characteristic curve, decision curve analysis, consistency index, and calibration curve prove the distinct advantages of the nomogram. Conclusions: Our web-based dynamic nomogram effectively integrates patient molecular subtype and sociodemographic characteristics with clinical characteristics and guidance and can be easily used. ER-Negative should receive attention in diagnosing and treating BCLM.

1. Background

Breast cancer is the most diagnosed cancer and the second leading cause of cancer-related deaths among women in developed countries [1]. More than 90% of deaths caused by breast cancer are due to metastasis-related complications [2]. However, metastasis is a little-known process that begins when tumor cells are separated from the primary tumor and injected into the bloodstream [3]. These circulating tumor cells (CTCs) eventually stay in the capillary bed of distant organs and extravasate into the parenchyma through the blood vessel wall, resulting in metastatic colonies in secondary sites [4]. The bone, brain, liver, and lungs are more likely to metastasize by breast cancer, knowns as organ orientation [5]. For breast cancer patients with metastases, 30–60% of patients have lesions in the bones, 4–10% in the brain, 15–32% in the liver, 21–32% in the lungs, and in many reported cases—multiple simultaneous transfers [6].
Breast cancer lung metastases (BCLM) often occurs within five years of the initial breast cancer diagnosis and significantly impacts the morbidity and mortality of patients [7]. Physiologically, these metastases can disrupt normal lung function, leading to coughing, breathing difficulties, hemoptysis, and death. The prognosis for patients with metastases confined to the lungs is inferior, with a median survival of only 25 months [8]. It is estimated that 60–70% of patients who die of breast cancer have lung metastases [9]. Therefore, accurate risk stratification for identifying BCLM is essential for treatment selection and prognostic evaluation. In addition, many clinicopathological factors, especially the molecular subtypes of breast cancer, could affect the occurrence and prognosis of BC lung metastasis [10]. ER-positive and PR-positive patients are more likely to develop bone metastases, while patients with more aggressive subtypes such as HER-2-positive and triple-negative tumors are more likely to develop visceral and brain metastases [11]. Moreover, the prognostic value of molecular subtypes in breast cancer patients with lung metastases remains controversial [12]. In addition, the influence of sociodemographic factors and clinicopathological characteristics on the incidence and prognosis of newly diagnosed breast cancer lung metastasis is still unclear.
Accurate risk stratification of BCLM patients is essential for treatment selection and prognosis assessment. We intend to explore the risk factors for BCLM through big data, combined with sociodemographic factors and molecular subtypes of breast cancer, to determine the high-risk factors for BCLM patients at the diagnosis. As a visual representation of mathematical models, nomograms can integrate certain features to estimate specific endpoints and provide practical and comprehensive predictions for clinical practice. Therefore, we aim to create a SEER-based prognostic nomogram for BCLM patients and accurately and conveniently assess overall survival (OS) and cancer-specific survival (CSS).

2. Materials and Methods

2.1. Study Design and Patient Selection

We used the SEER*Stat 8.3.8 program to retrieve and download patient information from the SEER 18 database. The SEER*Stat program is a public country registration database containing cancer incidence data in 18 regions of the United States, accounting for approximately 34.6% of the total population, and detailed patient data.
The trial population included adult female breast cancer patients diagnosed with lung metastases from 2010 to 2016. We included the patients’ race, age, marital status, detailed tumor node metastasis (TNM) staging, and whether they had undergone primary site radical resection, chemotherapy, ER, PR, HER2 receptor status, bone metastasis, brain metastasis, or liver metastasis. We excluded patients who had no pathological results and no complete information. Finally, 3916 patients were included in our study, and the selection flow chart is shown in Figure 1. This type of retrospective research does not require an ethics committee review.

2.2. Statistical Analysis and Dynamic Nomogram Publication

Patients were randomly divided into 2:1 training and verification groups, 2611 and 1305 cases. Univariate Cox proportional hazards regression analysis was developed to identify independent prognostic factors to construct prognostic factors. Based on univariate analysis results (p-value < 0.1), a multivariate Cox proportional hazard regression analysis was performed to construct a nomogram of significant variables (p-value < 0.05) in the training group. We use 1-year, 3-year, and 5-year OS and CSS for analysis in the nomogram. The consistency index (C-index) is used to evaluate the discrimination ability and accuracy of the nomogram. The identification and calibration are evaluated by bootstrapping 1000 times. Decision curve analysis (DCA) is used to evaluate the advantages and benefits of our new forecasting model compared with another single factor. In addition, time-dependent ROC is used to evaluate the effectiveness of the OS and CSS of the models we built. After validation, the web-based dynamic nomograms for BCLM OS and CSS are published online through the r package “DynNom” and shinyapps.io (https://www.shinyapps.io/) (accessed on 20 December 2022).
All these statistical methods are packaged using r version 4.2.0 (http://www.r-project.org) (accessed on 20 December 2022). In the two-tailed test, statistical significance is set to p < 0.05.

3. Results

3.1. Patient Characteristics

Our study excluded 2144 cases initially (lack of autopsy or death certificate diagnosis (n = 5); survival months are incomplete, or 0 (n = 482); lack of histological grade (n = 725); lacking ER, PR, and HER2 status (n = 358); T0 and Tx according to the 6th edition AJCC staging (n = 574) (Figure 1)), and 3916 cases were finally included (2611 cases were training cohort and 1305 cases were Verification cohort). Table 1 summarizes the patient characteristics. Among all BCLM patients, the median age was 62, 41.2% were married, and 72.3% were white. Pathologically, more than half of the patients were poorly differentiated tumors (55.8%), and 42.2% of patients had T4. Only 30.2% of patients underwent surgery in terms of treatment, but most patients underwent chemotherapy (62.7%). Bone metastasis is the most common in patients with BCLM (53.5%), followed by liver metastasis (26.7%). The probability of HER2 negative is the highest, more than 70.7%. In addition, there are 460 triple-negative breast cancer (TNBC) patients, accounting for 11.7%. The longest follow-up of the entire cohort is 82 months, the average OS is 20.9 months, and the CSS is similar to the OS, which is 21.0 months.

3.2. Establishment of a Prognostic Nomogram

Univariate and multivariate Cox regression analysis is used to calculate the weights of variables in OS and CSS (expressed as OR) and identify independent risk factors. Variables with significant differences in univariate analysis are included in the Cox regression model of multivariate analysis, in which OS and CSS are significantly related to 13 variables, namely age, marital status, race, grade, T stage, surgery, chemotherapy, bone metastatic, brain metastatic, liver metastatic, ER, PR, and HER2. Patient gender, insurance status, lymph node metastasis, and radiation therapy were excluded. (Table 2 and Table 3). All important variables are integrated to build a nomogram of OS and CSS. The nomograms of 1-year, 3-year, and 5-year OS are shown in Figure 2A, and the nomograms of 1-year, 3-year, and 5-year CSS are shown in Figure 3A. After adding up the scores associated with each variable and placing the total score in the bottom tier, the probability of OS and CSS can be estimated at 1, 3, and 5 years.

3.3. Verification of Prognostic Nomogram

The concordance index (C-index), calibration curves, time-dependent receiver operating characteristic (ROC) curves, and decision curve analysis (DCA) methods were used to verify the superiority of our nomogram.
In the C index evaluation, the training group’s OS and CSS C indexes were 0.702 (95% CI, 0.687–0.716) and 0.710 (95% CI, 0.693–0.727). The OS and CSS C indexes of the validation group were, respectively, 0.703 (95% CI, 0.682–0.723), and 0.715 (95% CI, 0.692–0.737) (Table 4).
Furthermore, high reliability was demonstrated by calibration curves in OS (Figure 2B,D) and CSS (Figure 3B,D). Time-dependent ROC curves determine the sensitivity and specificity of predicting the prognosis of BCLM. Figure 2C,F shows the 1-year, 3-year, and 5-year values of the area under the curve of the OS nomogram (training group: 1-year operating system, 75.38%; 3-year OS, 74.37%; and 5-year OS, 74.18%; Verification Group: 1-year OS, 75.81%; 3-years OS, 76.33%; and 5-years OS, 74.98%). The AUC value of the predicted CSS nomogram is shown in Figure 3 (training group: 1-year CSS, 75.85%; 3-year CSS, 75.64%; and 5-year CSS, 76.59%; Verification group: 1-year CSS, 77.10%; 3-year period CSS, 77.78%; and 5-year CSS, 77.97%). In addition, due to the slight deviation from the reference line, Figure 2D,E is for the OS and Figure 3D,E is for CSS.
DCA can compare nomograms and other factors to help clinicians make beneficial decisions. Figure 2D,G shows the DCA curves of OS for nomogram and other clinical factors, and Figure 3D,G represents CSS. Compared with all other clinical factors, the DCA of the nomogram showed a superior net benefit, suggesting an excellent clinical application of the nomogram in this study.

3.4. Risk Stratification and Web-Based Dynamic Nomogram Publication

X-tile software (version 3.6.1; Yale University, New Haven, CT, USA) was used to calculate the cut-off value by adding the scores associated with each variable to obtain the total score for patients with BCLM. The cut-off values for OS are 347 and 464, and the cut-off values for CSS are 394 and 548 (Figure 4A,B). Therefore, BCLM patients are classified as high-risk (score > 464), medium-risk (347 < score ≤ 464), and low-risk (score ≤ 347) for OS. In addition, BCLM patients were classified as high-risk (score > 548), medium-risk (394 < score ≤ 548), and low-risk (score ≤ 394) for CSS (Figure 4C,D).
Based on risk stratification, Kaplan–Meier survival curves were plotted for all BCLM patients, as shown in Figure 5. The 5-year OS rate is the highest in the low-risk group at 34.5%, followed by the intermediate-risk group at 12.9%, and the worst at 4.9% in the high-risk group (Figure 5A). Likewise, the low-risk group had the highest 5-year CSS rate at 40.1%, followed by the intermediate-risk group with 14.5%, and the high-risk group with 4.6% (Figure 5B). There was a statistically significant difference in survival outcomes among the three groups (p < 0.001). We finally published our web-based nomogram for predicting OS in BCLM patients (Figure 5C) and CSS (Figure 5D), which can be accessed at (https://nomogram-xiangyahospital.shinyapps.io/BCLMOSnomogram) (accessed on 20 December 2022) and (https://cssnomogram-xiangyahospital.shinyapps.io/DynNomapp) (accessed on 20 December 2022) to conduct your data calculations and risk analysis. A step-by-step instruction was produced to help researchers use the site (Supplementary File).

4. Discussion

Risk assessment tools and surveillance methods for BCLM remain limited. We constructed a nomogram based on extensive population data from the SEER database and included breast cancer molecular subtypes, treatment modalities, and sociodemographic factors. Although Xie’s study found disease risk factors based on similar data, it did not build a complete risk assessment tool [12]. We developed and constructed a nomogram of an accurate scoring system with clinical value and achieved good classification results for OS and CSS in BCLM patients. We validated our model with multiple methods and found that it works well to provide comprehensive guidance for clinical practice. We published it online, so all readers can easily access and use it.
Our study found that surgical treatment is protective for OS and CSS in BCLM patients. Surgery can prolong the survival time of patients and reduce the risk of death by 29.1%. Several randomized clinical trials and retrospective studies have shown that primary tumor surgery can improve cancer survival by reducing the tumor burden of metastatic breast cancer [13,14,15]. However, selection bias cannot be ignored, such as selecting patients with smaller tumor burden, younger age, or other more favorable characteristics for surgery, which may lead to the final result being more biased towards the surgery group [16,17]. Some basic research suggests that surgical injury may upregulate genes involved in adhesion, invasion, and angiogenesis, leading to breast cancer metastasis in the lung [18]. Although the benefit of primary tumor surgery in patients with newly diagnosed breast cancer and lung metastases is unclear, our results suggest that BCLM patients undergoing surgical treatment have a survival benefit.
Our study found significant differences in the prognosis of BCLM patients of different molecular subtypes. We separately analyzed the three crucial therapeutic targets of breast cancer, ER, PR, and HER2, and the positive of the three reduced the risk of death by 48.9%, 44.2%, and 25.8%. This suggests that BCLM patients pay more attention to evaluating ER molecular targets. There is evidence of an interaction between the ER and HER2 pathways [19]. Zhao’s research also found that molecular typing strongly impacts metastatic breast cancer [20]. The ER pathway can be used as a bypass activation mechanism for downstream signals of the HER2 pathway [21]. Activating the HER2 signaling pathways could further promote the activity of the ER pathway, leading to impaired endocrine therapy response and possibly changing the tumor’s response to HER2-targeted therapy [22]. Some neoadjuvant clinical studies targeting anti-HER2 have shown that, compared with HER2+/ER- patients, HER2+/ER+ patients have a lower pathological complete response rate (pCR) [23]. It also suggests that our breast cancer patients with lung metastasis may be closely related to ER expression in molecular mechanisms.
In sociodemography, we found that marital status is an independent factor for BCLM patients. Marriage could reduce the risk of death by 25%, similar to our previous research on colorectal cancer [24]. However, more than 99% of breast cancer patients are women, and the hormone cycle and stable level are closely related to the occurrence and development of breast cancer [25]. Unmarried patients often lack the care and support of their spouses and are more likely to suffer from chronic psychological distress, such as depression and helplessness and unhealthy living habits, such as smoking and alcoholism [26]. Reducing psychosocial support and psychological stress will affect the immune and endocrine systems and accelerate tumor growth and patient death [27]. Furthermore, with the further development of the industrialization process, the human marriage rate may be further reduced [28]. Therefore, in clinical work, providing psychosocial support from medical workers may help improve the survival of patients.
Our study randomly selected one-third of the patients as the validation group to confirm the superiority of the nomogram in this study. The excellent results of the verification group, including the C index, time-dependent ROC curve, DCA, and calibration curve, ensure the generalizability of the new nomogram. We made our models as web-based dynamic nomograms, making our models available to everyone online. However, there are still some limitations in our research. First, as a retrospective study, the nomogram still needs to be verified by prospective studies in the future. Secondly, much information is missing in our research, such as the lack of metastatic tumor size, location, and treatment response, but with the further improvement of the database, these problems will be further solved. In addition, we need more accurate data to verify the validity of the nomogram. Despite these limitations, the findings demonstrate the clinical value of the sensitivity and specificity of our constructed nomogram.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jpm13010043/s1, File: Web-based dynamic nomograms for predicting overall survival and cancer-specific survival in breast cancer patients with lung metastases.

Author Contributions

Conception and design: K.W., Y.L., D.W. and Z.Z.; administrative support: Z.Z.; provision of study materials or patients: Y.L.; collection and assembly of data: K.W. and Y.L.; data analysis and Interpretation: K.W., Y.L. and D.W.; writing—original draft: K.W., Y.L., D.W. and Z.Z.: writing—review and editing: all authors final approval of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Yuqiang Li was granted for this project, supported by funds from the Natural Science Foundation of Hunan Province (No. 2022JJ40799). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All patient files are available from the SEER database. More details supporting the findings of this study are available from the corresponding author upon request.

Acknowledgments

We would like to thank Yuan Zhou for administrative support.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Redig, A.J.; McAllister, S.S. Breast cancer as a systemic disease: A view of metastasis. J. Intern. Med. 2013, 274, 113–126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Richard, V.; Davey, M.G.; Annuk, H.; Miller, N.; Kerin, M.J. The double agents in liquid biopsy: Promoter and informant biomarkers of early metastases in breast cancer. Mol. Cancer 2022, 21, 95. [Google Scholar] [CrossRef] [PubMed]
  4. Fu, A.; Yao, B.; Dong, T.; Chen, Y.; Yao, J.; Liu, Y.; Li, H.; Bai, H.; Liu, X.; Zhang, Y.; et al. Tumor-resident intracellular microbiota promotes metastatic colonization in breast cancer. Cell 2022, 185, 1356–1372. [Google Scholar] [CrossRef]
  5. Molnar, I.A.; Molnar, B.A.; Vizkeleti, L.; Fekete, K.; Tamas, J.; Deak, P.; Szundi, C.; Szekely, B.; Moldvay, J.; Vari-Kakas, S.; et al. Breast carcinoma subtypes show different patterns of metastatic behavior. Virchows Arch. 2017, 470, 275–283. [Google Scholar] [CrossRef]
  6. Treska, V.; Cerna, M.; Liska, V.; Treskova, I.; Narsanska, A.; Bruha, J. Surgery for breast cancer liver metastases—Factors determining results. Anticancer Res. 2014, 34, 1281–1286. [Google Scholar]
  7. Cardoso, F.; Senkus, E.; Costa, A.; Papadopoulos, E.; Aapro, M.; Andre, F.; Harbeck, N.; Aguilar Lopez, B.; Barrios, C.H.; Bergh, J.; et al. 4th ESO-ESMO International Consensus Guidelines for Advanced Breast Cancer (ABC 4)dagger. Ann. Oncol. 2018, 29, 1634–1657. [Google Scholar] [CrossRef]
  8. Medeiros, B.; Allan, A.L. Molecular Mechanisms of Breast Cancer Metastasis to the Lung: Clinical and Experimental Perspectives. Int. J. Mol. Sci. 2019, 20, 2272. [Google Scholar] [CrossRef] [Green Version]
  9. DeSantis, C.E.; Ma, J.; Gaudet, M.M.; Newman, L.A.; Miller, K.D.; Goding Sauer, A.; Jemal, A.; Siegel, R.L. Breast cancer statistics, 2019. CA Cancer J. Clin. 2019, 69, 438–451. [Google Scholar] [CrossRef] [Green Version]
  10. Hess, K.R.; Varadhachary, G.R.; Taylor, S.H.; Wei, W.; Raber, M.N.; Lenzi, R.; Abbruzzese, J.L. Metastatic patterns in adenocarcinoma. Cancer 2006, 106, 1624–1633. [Google Scholar] [CrossRef]
  11. Mathew, A.; Rajagopal, P.S.; Villgran, V.; Sandhu, G.S.; Jankowitz, R.C.; Jacob, M.; Rosenzweig, M.; Oesterreich, S.; Brufsky, A. Distinct Pattern of Metastases in Patients with Invasive Lobular Carcinoma of the Breast. Geburtshilfe Frauenheilkd 2017, 77, 660–666. [Google Scholar] [CrossRef] [PubMed]
  12. Xiao, W.; Zheng, S.; Liu, P.; Zou, Y.; Xie, X.; Yu, P.; Tang, H.; Xie, X. Risk factors and survival outcomes in patients with breast cancer and lung metastasis: A population-based study. Cancer Med. 2018, 7, 922–930. [Google Scholar] [CrossRef] [PubMed]
  13. Badwe, R.; Hawaldar, R.; Nair, N.; Kaushik, R.; Parmar, V.; Siddique, S.; Budrukkar, A.; Mittra, I.; Gupta, S. Locoregional treatment versus no treatment of the primary tumour in metastatic breast cancer: An open-label randomised controlled trial. Lancet Oncol. 2015, 16, 1380–1388. [Google Scholar] [CrossRef] [PubMed]
  14. Reinhorn, D.; Mutai, R.; Yerushalmi, R.; Moore, A.; Amir, E.; Goldvaser, H. Locoregional therapy in de novo metastatic breast cancer: Systemic review and meta-analysis. Breast 2021, 58, 173–181. [Google Scholar] [CrossRef] [PubMed]
  15. Li, K.; Zhou, C.; Yu, Y.; Niu, L.; Zhang, W.; Wang, B.; He, J.; Ge, G. Metastatic Pattern Discriminates Survival Benefit of Type of Surgery in Patients With De Novo Stage IV Breast Cancer Based on SEER Database. Front. Surg. 2021, 8, 696628. [Google Scholar] [CrossRef] [PubMed]
  16. Khan, S.A. Surgical Management of de novo Stage IV Breast Cancer. Semin. Radiat. Oncol. 2016, 26, 79–86. [Google Scholar] [CrossRef] [PubMed]
  17. Soran, A.; Ozmen, V.; Ozbas, S.; Karanlik, H.; Muslumanoglu, M.; Igci, A.; Canturk, N.Z.; Utkan, Z.; Evrensel, T.; Sezgin, E.; et al. Primary Surgery with Systemic Therapy in Patients with de Novo Stage IV Breast Cancer: 10-year Follow-up; Protocol MF07-01 Randomized Clinical Trial. J. Am. Coll Surg. 2021, 233, 742–751e5. [Google Scholar] [CrossRef]
  18. Al-Sahaf, O.; Wang, J.H.; Browne, T.J.; Cotter, T.G.; Redmond, H.P. Surgical injury enhances the expression of genes that mediate breast cancer metastasis to the lung. Ann. Surg. 2010, 252, 1037–1043. [Google Scholar] [CrossRef]
  19. Yamamura, J.; Kamigaki, S.; Fujita, J.; Osato, H.; Manabe, H.; Tanaka, Y.; Shinzaki, W.; Hashimoto, Y.; Komoike, Y. New insights into patterns of first metastatic sites influencing survival of patients with hormone receptor-positive, HER2-negative breast cancer: A multicenter study of 271 patients. BMC Cancer 2021, 21, 476. [Google Scholar] [CrossRef]
  20. Zhao, W.; Wu, L.; Zhao, A.; Zhang, M.; Tian, Q.; Shen, Y.; Wang, F.; Wang, B.; Wang, L.; Chen, L.; et al. A nomogram for predicting survival in patients with de novo metastatic breast cancer: A population-based study. BMC Cancer 2020, 20, 982. [Google Scholar] [CrossRef]
  21. Deluche, E.; Antoine, A.; Bachelot, T.; Lardy-Cleaud, A.; Dieras, V.; Brain, E.; Debled, M.; Jacot, W.; Mouret-Reynier, M.A.; Goncalves, A.; et al. Contemporary outcomes of metastatic breast cancer among 22,000 women from the multicentre ESME cohort 2008-2016. Eur. J. Cancer 2020, 129, 60–70. [Google Scholar] [CrossRef] [PubMed]
  22. Hoogstraat, M.; Lips, E.H.; Mayayo-Peralta, I.; Mulder, L.; Kristel, P.; van der Heijden, I.; Annunziato, S.; van Seijen, M.; Nederlof, P.M.; Sonke, G.S.; et al. Comprehensive characterization of pre- and post-treatment samples of breast cancer reveal potential mechanisms of chemotherapy resistance. NPJ Breast Cancer 2022, 8, 60. [Google Scholar] [CrossRef]
  23. Bidard, F.C.; Kaklamani, V.G.; Neven, P.; Streich, G.; Montero, A.J.; Forget, F.; Mouret-Reynier, M.A.; Sohn, J.H.; Taylor, D.; Harnden, K.K.; et al. Elacestrant (oral selective estrogen receptor degrader) Versus Standard Endocrine Therapy for Estrogen Receptor-Positive, Human Epidermal Growth Factor Receptor 2-Negative Advanced Breast Cancer: Results From the Randomized Phase III EMERALD Trial. J. Clin. Oncol. 2022, 40, JCO2200338. [Google Scholar] [CrossRef] [PubMed]
  24. Li, Y.; Liu, W.; Zhao, L.; Gungor, C.; Xu, Y.; Song, X.; Wang, D.; Zhou, Z.; Zhou, Y.; Li, C.; et al. Nomograms predicting Overall Survival and Cancer-specific Survival for Synchronous Colorectal Liver-limited Metastasis. J. Cancer 2020, 11, 6213–6225. [Google Scholar] [CrossRef] [PubMed]
  25. Goldzweig, G.; Andritsch, E.; Hubert, A.; Brenner, B.; Walach, N.; Perry, S.; Baider, L. Psychological distress among male patients and male spouses: What do oncologists need to know? Ann. Oncol. 2010, 21, 877–883. [Google Scholar] [CrossRef]
  26. Chen, Z.H.; Yang, K.B.; Zhang, Y.Z.; Wu, C.F.; Wen, D.W.; Lv, J.W.; Zhu, G.L.; Du, X.J.; Chen, L.; Zhou, G.Q.; et al. Assessment of Modifiable Factors for the Association of Marital Status With Cancer-Specific Survival. JAMA Netw. Open 2021, 4, e2111813. [Google Scholar] [CrossRef]
  27. Ellis, L.; Canchola, A.J.; Spiegel, D.; Ladabaum, U.; Haile, R.; Gomez, S.L. Racial and Ethnic Disparities in Cancer Survival: The Contribution of Tumor, Sociodemographic, Institutional, and Neighborhood Characteristics. J. Clin. Oncol. 2018, 36, 25–33. [Google Scholar] [CrossRef]
  28. Hellstrand, J.; Nisen, J.; Myrskyla, M. Less Partnering, Less Children, or Both? Analysis of the Drivers of First Birth Decline in Finland Since 2010. Eur. J. Popul. 2022, 38, 191–221. [Google Scholar] [CrossRef]
Figure 1. The flowchart of patients was identified in the study. A total of 6060 BCML patients were included from the SEER database, and we excluded 2144 patients based on exclusion criteria. The 3916 patients who met the requirements were randomly divided into training and validation cohort by 2:1. The training cohort was used to construct Norman diagrams, and the validation cohort was used to evaluate the effect of the building model. We validated our constructed Norman plots using the Calibration curve, the C-index decision curve analysis (DAC), and time-dependent ROC methods.
Figure 1. The flowchart of patients was identified in the study. A total of 6060 BCML patients were included from the SEER database, and we excluded 2144 patients based on exclusion criteria. The 3916 patients who met the requirements were randomly divided into training and validation cohort by 2:1. The training cohort was used to construct Norman diagrams, and the validation cohort was used to evaluate the effect of the building model. We validated our constructed Norman plots using the Calibration curve, the C-index decision curve analysis (DAC), and time-dependent ROC methods.
Jpm 13 00043 g001
Figure 2. Nomogram and model assessment for predicting 1-, 3-, and 5-year overall survival (OS) in BCLM patients. (A) The OS nomogram for BCLM patients. The calibration plots in the (B) training and (E) validation cohorts for 1-year, 3-year, and 5-year OS. AUC values of ROCs of the nomograms for 1-year, 3-year, and 5-year OS. The values in brackets B and E represent the area under the ROC curves (AUC). (C) for the training cohort and (F) for the validation cohort. Decision curve analyses (DCA) of the nomogram and other factors’ overall survival. (D) for the training cohort and (G) for the validation cohort.
Figure 2. Nomogram and model assessment for predicting 1-, 3-, and 5-year overall survival (OS) in BCLM patients. (A) The OS nomogram for BCLM patients. The calibration plots in the (B) training and (E) validation cohorts for 1-year, 3-year, and 5-year OS. AUC values of ROCs of the nomograms for 1-year, 3-year, and 5-year OS. The values in brackets B and E represent the area under the ROC curves (AUC). (C) for the training cohort and (F) for the validation cohort. Decision curve analyses (DCA) of the nomogram and other factors’ overall survival. (D) for the training cohort and (G) for the validation cohort.
Jpm 13 00043 g002
Figure 3. Nomogram and model assessment for predicting 1-, 3-, and 5-year cancer-specific survival (CSS) in BCLM patients. (A) The CSS nomogram for BCLM patients. The calibration plots in the (B) training and (E) validation cohorts for 1-year, 3-year, and 5-year CSS. AUC values of ROCs of the nomograms for 1-year, 3-year, and 5-year OS. The values in brackets B and E represent the area under the ROC curves (AUC). (C) for the training cohort and (F) for the validation cohort. Decision curve analyses (DCA) of the nomogram and other factors’ overall survival. (D) for the training cohort and (G) for the validation cohort.
Figure 3. Nomogram and model assessment for predicting 1-, 3-, and 5-year cancer-specific survival (CSS) in BCLM patients. (A) The CSS nomogram for BCLM patients. The calibration plots in the (B) training and (E) validation cohorts for 1-year, 3-year, and 5-year CSS. AUC values of ROCs of the nomograms for 1-year, 3-year, and 5-year OS. The values in brackets B and E represent the area under the ROC curves (AUC). (C) for the training cohort and (F) for the validation cohort. Decision curve analyses (DCA) of the nomogram and other factors’ overall survival. (D) for the training cohort and (G) for the validation cohort.
Jpm 13 00043 g003
Figure 4. The cut-off values were calculated using X-tile based on the total scores. (A) The black dot in the figure indicates the OS distinction point selected by X-tile software. (B) The cut-off values of OS are 347 and 464. The blue color is the low-risk group, the gray color is the medium-risk group, and the purple color is the high-risk group. (C) The black dots in the figure represent the CSS distinguishing points selected by X-tile software. (D) The cut-off values of CSS are 394 and 548. OS, overall survival; the blue color is the low-risk group, the gray color is the medium-risk group, and the purple color is the high-risk group. CSS, cancer-specific survival. All patients under our group price adjustment are under the OS, CSS curve.
Figure 4. The cut-off values were calculated using X-tile based on the total scores. (A) The black dot in the figure indicates the OS distinction point selected by X-tile software. (B) The cut-off values of OS are 347 and 464. The blue color is the low-risk group, the gray color is the medium-risk group, and the purple color is the high-risk group. (C) The black dots in the figure represent the CSS distinguishing points selected by X-tile software. (D) The cut-off values of CSS are 394 and 548. OS, overall survival; the blue color is the low-risk group, the gray color is the medium-risk group, and the purple color is the high-risk group. CSS, cancer-specific survival. All patients under our group price adjustment are under the OS, CSS curve.
Jpm 13 00043 g004
Figure 5. Risk stratification by the total score for subgroup survival analysis and web-based dynamic nomogram. (A) The 5-year OS rate is the highest in the low-risk group at 34.5%, followed by the intermediate-risk group at 12.9%, and the worst at 4.9% in the high-risk group. (B) The low-risk group had the highest 5-year CSS rate at 40.1%, followed by the intermediate-risk group with 14.5%, and the high-risk group with 4.6%. (C) Web-based dynamic nomogram of OS of BCLM patients (https://nomogram-xiangyahospital.shinyapps.io/BCLMOSnomogram) (accessed on 20 December 2022). (D) Web-based dynamic nomogram of CSS of BCLM patients (https://cssnomogram-xiangyahospital.shinyapps.io/DynNomapp) (accessed on 20 December 2022).
Figure 5. Risk stratification by the total score for subgroup survival analysis and web-based dynamic nomogram. (A) The 5-year OS rate is the highest in the low-risk group at 34.5%, followed by the intermediate-risk group at 12.9%, and the worst at 4.9% in the high-risk group. (B) The low-risk group had the highest 5-year CSS rate at 40.1%, followed by the intermediate-risk group with 14.5%, and the high-risk group with 4.6%. (C) Web-based dynamic nomogram of OS of BCLM patients (https://nomogram-xiangyahospital.shinyapps.io/BCLMOSnomogram) (accessed on 20 December 2022). (D) Web-based dynamic nomogram of CSS of BCLM patients (https://cssnomogram-xiangyahospital.shinyapps.io/DynNomapp) (accessed on 20 December 2022).
Jpm 13 00043 g005
Table 1. Demographics and clinicopathologic characteristics of the cohort with breast cancer lung metastatic.
Table 1. Demographics and clinicopathologic characteristics of the cohort with breast cancer lung metastatic.
VariablesTotal Cohort (n = 3916)Training Cohort (n = 2611)Validation Cohort (n = 1305)
No%No%No%
Age
≤5078820.151619.827220.8
51–65155639.7104540.051139.2
>65157240.1105040.252240.0
Gender
Female385398.4257398.5128098.1
Male631.6381.5251.9
Marital status
Married161341.2107341.154041.4
Unmarried230358.8153858.976558.6
Insurance
No2285.81465.6826.3
Yes368894.2246594.4122393.7
Race
White283072.3188672.294472.3
Black72918.646918.026020.0
Other3579.12569.81017.7
Grade
I/II172944.2114944.058044.4
III/IV218755.8146256.072555.6
T stage
T143311.129411.313910.7
T2114029.177929.836127.7
T369117.645217.323918.3
T4165242.2108641.656643.3
N stage
N086122.058022.228121.5
N1193149.3127949.065250.0
N251613.233612.918013.8
N360815.541615.919214.7
Surgery
No271969.4180769.291269.9
Yes119730.680430.839330.1
Chemotherapy
No146237.396537.049738.1
Yes245462.7164663.080861.9
Radiation
No168943.175428.993571.6
Yes222756.9185771.137028.4
Bone metastatic
No182946.7124447.658544.8
Yes208753.3136752.472055.2
Brain metastatic
No355990.9237090.8118991.1
Yes3579.12419.21168.9
Liver metastatic
No287273.3192173.695172.9
Yes 104426.769026.435427.1
ER
Negative119430.580330.839130.0
Positive272269.5180869.291470.0
PR
Negative170543.5113243.457343.9
Positive221156.5147956.673256.1
HER2
Negative277070.7182469.994672.5
Positive114629.378730.135927.5
OS (months) 20.9420.920.9
CSS (months) 21.0121.020.9
For marital status, unmarried consists of single, divorced, separated, and widowed. For race, ‘other’ includes American Indian, AK Native, Asian, and Pacific Islander. ER, estrogen receptor. PR, progesterone receptor. HER2, human epidermal growth factor receptor 2.
Table 2. Univariable and multivariable Cox regression model analyses of OS for nomogram.
Table 2. Univariable and multivariable Cox regression model analyses of OS for nomogram.
Characteristics Univariable Analysis Multivariable Analysis
OR95% CIpOR95% CIp
Age <0.001 <0.001
≤50 Reference1 Reference1
51–651.1320.976–1.3120.1011.1430.984–1.3280.081
>651.4391.244–1.665<0.0011.5171.297–1.774<0.001
Gender 0.642
Female Reference1 -
Male0.9020.597–1.362
Marital status <0.001 <0.001
Mariied Reference1 Reference1
Unmarried1.3271.192–1.476<0.0011.2541.124–1.399<0.001
Insurance 0.21
No Reference1 -
Yes1.1440.927–1.411
Race <0.001 <0.001
White Reference1 Reference1
Black1.3171.157–1.499<0.0011.3201.153–1.511<0.001
Other0.8480.700–1.0290.0950.8380.690–1.0180.075
Grade <0.001 0.008
I/II Reference1 Reference1
III/IV1.5261.372–1.696<0.0011.4801.311–1.670<0.001
T stage <0.001 <0.001
T1 Reference1 Reference1
T21.1490.950–1.3900.1531.0940.904–1.3250.355
T31.3431.093–1.649<0.0011.1840.963–1.4570.109
T41.4861.239–1.782<0.0011.3011.083–1.5630.005
N stage 0.82
N0 Reference1 -
N11.008757430.884–1.1520.897
N20.9518465990.794–1.1410.593
N31.0443373760.883–1.2350.612
Surgery <0.001 <0.001
No Reference1 Reference1
Yes0.7110.635–0.796<0.0010.6670.589–0.755<0.001
Chemotherapy <0.001 <0.001
No Reference1 Reference1
Yes0.8110.730–0.901<0.0010.6670.589–0.755<0.001
Radiation 0.356
No Reference1 -
Yes0.9490.849–1.0610.356
Bone metastatic 0.006 <0.001
No Reference1 Reference1
Yes1.1561.042–1.2820.0061.2571.122–1.409<0.001
Brain metastatic <0.001 <0.001
No Reference1 Reference1
Yes1.8451.567–2.171<0.0011.7531.482–2.072<0.001
Liver metastatic <0.001 <0.001
No Reference1 Reference1
Yes 1.7071.527–1.908<0.0011.6941.507–1.905<0.001
ER <0.001 <0.001
Negative Reference1 Reference1
Positive0.5270.473–0.587<0.0010.5270.448–0.619<0.001
PR <0.001 <0.001
Negative Reference1 Reference1
Positive0.5390.535–0.658<0.0010.7560.647–0.882<0.001
HER2 <0.001 <0.001
Negative Reference1 Reference1
Positive0.7290.647–0.821<0.0010.6080.535–0.692<0.001
Table 3. Univariable and multivariable Cox regression model analyses of CSS for nomogram.
Table 3. Univariable and multivariable Cox regression model analyses of CSS for nomogram.
CharacteristicsUnivariable Analysis Multivariable Analysis
OR95% CIpOR95% CIp
Age 0.007 0.002
≤50 Reference1 Reference1
51–651.1140.948–1.3090.1011.113 0.945–1.3120.201
>651.2851.092–1.513<0.0011.349 1.132–1.6080.001
Gender 0.642
Female Reference1 -
Male0.9020.597–1.362
Marital status <0.001 <0.001
Mariied Reference11.270 Reference1
Unmarried1.3411.186–1.517<0.001 1.12–1.441<0.001
Insurance 0.21
No Reference1 -
Yes1.1440.927–1.411
Race <0.001 <0.001
White Reference11.316 Reference1
Black1.3691.183–1.583<0.0010.818 1.129–1.534<0.001
Other0.860.689–1.0730.181 0.654–1.0240.08
Grade <0.001 <0.001
I/II Reference11.534 Reference1
III/IV1.5691.388–1.774<0.001 1.333–1.764<0.001
T stage <0.001 0.002
T1 Reference11.219 Reference1
T21.3191.017–1.7110.0371.340 0.938–1.5830.1390728
T31.6171.231–2.1230.0011.507 1.018–1.7640.0365868
T41.7751.383–2.279<0.001 1.172–1.9380.0013814
N stage 0.82
N0 Reference1 -
N11.008757430.884–1.1520.897
N20.9518465990.794–1.1410.593
N31.0443373760.883–1.2350.612
Surgery <0.001 <0.001
No Reference10.667 Reference1
Yes0.6660.584–0.759<0.001 0.590–0.777<0.001
Chemotherapy <0.001 <0.001
No Reference10.687 Reference1
Yes0.8660.767–0.9790.021 0.594–0.794<0.001
Radiation 0.356
No Reference1 -
Yes0.9490.849–1.0610.356
Bone metastatic 0.002 <0.001
No Reference11.325 Reference1
Yes1.211.072–1.3630.002 1.161–1.511<0.001
Brain metastatic <0.001 <0.001
No Reference11.868 Reference1
Yes2.0391.700–2.446<0.001 1.549–2.253<0.001
Liver metastatic <0.001 <0.001
No Reference11.655 Reference1
Yes 1.7691.559–2.006<0.001 1.449–1.891<0.001
ER <0.001 <0.001
Negative Reference10.550 Reference1
Positive0.5110.452–0.578<0.001 0.458–0.661<0.001
PR <0.001 <0.001
Negative Reference10.694 Reference1
Positive0.5580.496–0.629<0.001 0.580–0.830<0.001
HER2 <0.001 <0.001
Negative Reference10.595 Reference1
Positive0.7420.649–0.849<0.001 0.514–0.688<0.001
Table 4. The C-indices for predictions of OS and CSS.
Table 4. The C-indices for predictions of OS and CSS.
GroupOSCSS
C-Index95% CIC-Index95% CI
Training group0.7020.687–0.7160.7100.693–0.727
Validation group0.7030.682–0.7230.7150.692–0.737
OS, overall survival; CSS, cancer-specific survival; C-index, index of concordance; CI, confidence interval.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, K.; Li, Y.; Wang, D.; Zhou, Z. Web-Based Dynamic Nomograms for Predicting Overall Survival and Cancer-Specific Survival in Breast Cancer Patients with Lung Metastases. J. Pers. Med. 2023, 13, 43. https://doi.org/10.3390/jpm13010043

AMA Style

Wang K, Li Y, Wang D, Zhou Z. Web-Based Dynamic Nomograms for Predicting Overall Survival and Cancer-Specific Survival in Breast Cancer Patients with Lung Metastases. Journal of Personalized Medicine. 2023; 13(1):43. https://doi.org/10.3390/jpm13010043

Chicago/Turabian Style

Wang, Kangtao, Yuqiang Li, Dan Wang, and Zhongyi Zhou. 2023. "Web-Based Dynamic Nomograms for Predicting Overall Survival and Cancer-Specific Survival in Breast Cancer Patients with Lung Metastases" Journal of Personalized Medicine 13, no. 1: 43. https://doi.org/10.3390/jpm13010043

APA Style

Wang, K., Li, Y., Wang, D., & Zhou, Z. (2023). Web-Based Dynamic Nomograms for Predicting Overall Survival and Cancer-Specific Survival in Breast Cancer Patients with Lung Metastases. Journal of Personalized Medicine, 13(1), 43. https://doi.org/10.3390/jpm13010043

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop