Development of a nomogram based on serum cytokine-related riskscore in breast cancer

Background: Cytokines are involved in many inflammatory diseases and thus play an important role in tumor immune regulation. In recent years, researchers have found that breast cancer is not only related to genetic and environmental factors, but also to the chronic inflammation and immunity. However, the correlation between serum cytokines and blood tests indicators remain unclear.

Methods: A total of 84 serum samples and clinicopathological data of breast cancer patients from Tianjin Cancer Institute & Hospital, Tianjin Medical University, Tianjin, P. R. China were collected. The expression levels of the 12 cytokines were detected by immunofluorescence method. Blood tests results were obtained from medical records. By stepwise Cox regression analysis, a cytokine-related gene signature was generated. Univariate and multivariate Cox regression were used to analyze the influence on the prognosis of patients. A nomogram was constructed to illustrate the cytokine-related riskscore predicting 5-year OS, which was further evaluated and validated by C-index and ROC curve. The correlation between the expression of cytokines in serum and other blood indicators was studied by using Spearman’s test.

Results: The riskscore was calculated as IL-4×0.99069 + TNF-α×0.03683. Patients were divided into high and low risk groups according to the median riskscore, with the high-risk group has a shorter survival time by log-rank test (training set, P=0.017; validation set, P=0.013). Combined with the clinical characteristics, the riskscore was found to be an independent factor for predicting the OS of breast cancer patients in both training cohort (HR=1.2, P<0.01) and validation cohort (HR=1.6, P=0.023). The 5-year C-index and AUC of the nomogram were 0.78 and 0.68, respectively. IL-4 was further found to be negatively correlated with ALB.

Conclusion: In summary, we have developed a nomogram based on two cytokines including IL-4 and TNF-α to predict OS of breast cancer and investigated their correlation with blood test indicators.

Introduction

According to data released by the International Agency for Research on Cancer (IARC) in 2020, breast cancer has surpassed lung cancer for the first time, becoming the world’s most common cancer, accounting for approximately 11.7 per cent of new cancer cases. Breast cancer is not only the most common cancer in women, but also the leading cause of cancer-related death, which seriously threatens women’s life and health. Breast cancer is a disease possessed different morphology, molecular characteristics, behaviors, and therapeutic responses (1). There are many traditional clinical pathological variables and molecular predictive markers to reflect the prognosis. The feasibility of serum tumor markers in the diagnosis, prognosis and clinical treatment of malignant tumors has been extensively studied nowadays.

Cytokines are small molecular glycoproteins that are synthesized and secreted by different cells (mainly monocytes and activated T lymphocytes) in response to stimulation. Cancer cells can also produce cytokines. When cytokines affect the growth and function of immunocompetent cells, they can activate or modulate the anti-tumor response and act to promote or inhibit tumor growth (2).The serum levels of IL-6, IL-8 and Il-10 may be used as potential markers to predict the metastasis and prognosis of breast cancer in a Chinese population study (3). But the effect of cytokine expression on the prognosis of breast cancer patients is still controversial. IL-6 is one of the most studied cytokines in breast cancer. A meta-analysis of 3224 breast cancer patients showed that IL-6 expression was not associated with DFS (4). Though Ahmad N and his colleagues demonstrated that high IL-6/IL-10 expression was associated with better disease-free survival (5). And another study found that circulating IL-6 is associated with poor survival in patients with metastatic breast cancer (6). A study on TNBC indicated that simultaneous inhibition of IL-6 and IL-8 expression significantly inhibit tumor growth in vitro and in vivo. Multivariate Cox analysis showed that IL-6 and IL-8 expression can predict the survival time of patients (7). Increased serum IL-8 expression in patients with advanced breast cancer and multivariate analysis showed that serum IL-8 was an independent significant factor affecting the post relapse survival (8). Furthermore, in a study of 207 patients with invasive breast cancer, high IL-17 status was found to be a significant prognostic factor for DFS (9).

The presence of inflammatory cytokines in the tumor environment affects a series of processes at different stages of tumor progression, including initiation, proliferation, promotion, transformation, angiogenesis, invasion, inhibition of apoptosis, immune surveillance, drug resistance and metastasis (10). There was also a correlation between blood tests indicators and prognosis in breast cancer patients. Liu’s team found that ALB, LDH and total bilirubin were significantly associated with 5-year OS. Patients with higher albumin levels, higher total bilirubin levels, and lower LDH levels had a lower risk of death (11).A meta-analysis revealed a significant effect of higher serum LDH on poorer overall and progression-free survival (12). Therefore, the correlation between cytokines and blood indicators is worth exploring. Ultimately, we examined 12 different cytokines in the serum level and explored their prognostic impact, then analyzed their correlation with routine blood tests and liver kidney function tests. We hope our findings could serve to shed some light on the prognosis of breast cancer patients.

Materials and methods

Samples and clinical data collection

A total of 215 patients with available serum sample were initially included in our study from 2003 to 2020 in Tianjin Cancer Institute & Hospital, Tianjin Medical University, Tianjin, P. R. China.

According to the following inclusion criteria (1): pathological diagnosis of breast cancer (2); complete clinical and pathological records (3); no treatment or exclusion criteria (1): any acute or chronic infection (2); combined with other types of primary malignant tumors; and (3) patients with bilateral primary breast cancer, 131 patients were excluded and 84 samples were included in the study. The Medical Ethics Committee of our institute approved this study (ethics number E2013134).

The results of liver and kidney function tests (ALB, GLO, A/G, ChE, ADA, LDH, GPDA, 5’NT, LAP, β2-MG), routine blood examinations (NEUT, EO, BASO) and coagulation factor D-dimer from 84 blood samples were obtained from medical records. OS was defined as the time from the date of diagnosis to death or censoring. (The follow-up ended on 30 April 2021).

Measurement of cytokines

Venous blood samples were collected with standard test tubes, naturally solidified at room temperature or centrifuged at 2000-4000 rpm for 20 minutes, and about 0.5 mL of isolated serum was taken for examination. A Combined Cytokine Detection Kit (Ceger Biology) was adopted to simultaneously analyze the concentrations of 12 different cytokines in the serum of 84 samples following the manufacturer’s recommendations, in a flow cytometer. This kit is based on immunofluorescence technique. The mixture of capture microspheres contained 12 kinds of capture microspheres with different fluorescence intensities. The surface of capture microspheres was coated with 12 different cytokine-specific antibodies, the capture microspheres were specifically bound to the cytokines in the samples to be tested, and then combined with the PE-labeled fluorescent detection reagent. The capture microspheres, detection reagent, and detection antibody formed a double-frame sandwich complex, by analyzing the fluorescence intensity of the complex, the content of cytokines in the samples to be measured was obtained. The 12 cytokines included:IL⁃1β, IL-2, IL-4, IL-5, IL-6, IL-8, IL⁃10, IL-12p70, IL-17A, TNF-α, IFN-α and IFN⁃ γ.

Construction and validation of cytokine-related signature

Using the “sample” function of R, the cases were divided into training set (n = 58) and validation set (n = 26) by 7:3 random sampling without return. In the training set, with serum expression levels of the 12 cytokines as independent variables and overall survival as the outcome, optimized by Cox regression analysis using the “step” function of R, screening for cytokines that are strongly associated with overall survival in breast cancer. The median riskscore was used as the cut-off value to divide breast cancer patients into high-risk group and low-risk group. Log-rank tests and Kaplan-Meier curves were produced using the “survival” and “survminer” packages. The riskscores of the validation cohort were computed with the unified formula established in the training cohort. In the aftermath, subgroup survival analyses were performed on BC patients with different clinical characteristics (age, tumor stage and molecular subtype) to further assess the predictive value of the riskscore.

Construction of nomogram for breast cancer patients based on cytokine-related riskscore

In both training and validation sets, Cox regression analysis was performed with age, clinical stage, molecular subtype, and cytokine-related riskscore as independent variables, and overall survival as the outcome to construct a prognostic prediction model for breast cancer. Hazard ratio (HR), 95% confidence interval (CI), and p-values were separately determined. The “rms” package was applied to create a nomogram containing cytokine-related riskscore to predict the OS of BC. The C-index and score of the nomogram were calculated using the “survcomp” and “nomogramFormula” packages, respectively. ROC curves were plotted using the “timeRoc” package to further evaluate the predictive performance of the nomogram, riskscore and other clinical indicators.

Statistical analysis

All analyses were performed using IBM SPSS software 25.0 and R 4.2.0. The Chi-square tests were used to compare the clinicopathological characteristics between the training and validation sets. Youden index was used to determine the optimal cut-off point to divided IL-4 and TNF-α into low and high expression group. Univariate and multivariate Cox proportional hazards models were used to analyze prognostic risk factors. Differences in OS in subgroups were calculated using Kaplan-Meier plots and log-rank test by SPSS. We made a thermogram of the correlation between IL-4, TNF-α and blood test indicators by R (4.2.0).Then we used the Kaplan-Meier plotter database (13) (http://kmplot.com/analysis/ ) to calculate OS according to gene expression levels by using specific probes in different datasets. According to IL-4 expression level (Affymetrix ID, 207539_s_at; dataset, GSE20711), a total of 88 breast cancer patients were included. And TNFA (Affymetrix ID, 207113_s_at; dataset, GSE1456) involved 159 patients. All splitting patients by auto select best cutoff. A P-value < 0.05 was considered statistically significant. Because the data are randomly lost, the missing values (< 5%) of covariates are replaced by average values.

Results

The clinicopathological data

Based on the inclusion and exclusion criteria, a total of 84 blood samples were included in the study. They were divided into training set (n = 58) and validation set (n = 26). 28 (48.3%) patients were older than 53 years old in the training set (median age of patients in the training set), and there were 35(60.3%) patients in the I-II stage and 23(39.7%) in III-IV stage. According to molecular typing, patients including 18(31%) cases of TNBC molecular subtype, 10(17.2%) cases of Luminal a type, 17(29.3%) cases of Luminal B type and 13(22.5%) cases of Her-2 positive. The validation set data is also displayed in the Table 1. And we compared the clinicopathological characteristics between the training and validation cohorts, all the variables were comparable (P>0.05). Complete table of baseline characteristics (Supplementary Table 1) and other laboratory test results (Supplementary Table 2) are listed in Supplementary Data.

TABLE 1

Table 1 Characteristics of Breast Cancer patients in the training cohort and validation cohort.

Development and validation of the cytokine-related signature

Through multivariate Cox regression analysis and “step” function in R, a prediction model based on cytokine-related riskscore is obtained. The results showed that IL4 (HR=2.693, 95%CI 1.63-4.47, P<0.01) and TNF-α (HR=1.038, 95%CI 1.01-1.07, P<0.01) were powerful prognostic factors, and then we constructed the riskscore model. According to their coeficients and expression levels, we calculated riskscore for each patient in the training set as follows: riskscore =IL-4*0.99069 + TNF-α*0.03683. We used the median riskscore to divide the training set into risk-high group and risk-low group for survival analysis. The risk-high group had more deaths and shorter survival time than the risk-low group (Figure 1A). Significant differences in survival were detected between the two groups (P=0.017). The riskscore of the validation set was calculated according to the formula above. Next, the validation set is divided into high and low risk groups based on the median riskscore. The Kaplan–Meier analysis showed that patients in the low-risk group exhibited a significantly better OS than high-risk patients, indicating that the model had a good predictive eficacy (Figure 1B). We further analyzed the effect of serum IL-4 and TNF-α on OS and divided these two cytokines into high and low groups by calculating the cutoff value using Youden index. Our data revealed that both high serum IL-4 and TNF-α were indicative of poorer survival outcomes in comparison to low groups in the training set (P=0.006 and P=0.029, Figure 2).

FIGURE 1

Figure 1 The Overall survival curves of the low-risk and high-risk groups in the Training Cohort (A) and Validation Cohort (B).

FIGURE 2

Figure 2 The Overall survival curves of breast cancer patients with high and low serum IL-4 (A)/TNF-α (B) levels in the training cohort.

Survival analysis in subgroups

To evaluate the prognostic predictive value of the cytokine-related riskscore, a stratified subgroup analysis was performed according to clinical characteristics, as shown in Supplementary Figure 1. Both in young and old populations, the high-risk group had a poorer prognosis (Supplementary Figure 1A, age ≤ 53, P=0.017; age>53, P=0.047). And for patients with different tumor stages, the low-risk group had a better prognosis (Supplementary Figure 1B, stage I/II, P=0.067; stage III/IV, P=0.014), however, for patients in stage I/II, this difference was not significant. In different molecular types, the survival prognosis of high and low risk groups has some differences (Supplementary Figure 1C). In TNBC and HER2-positive patients, OS was significantly worse in the high-risk group than in the low-risk group (TNBC, P=0.011; HER2-positive, P=0.041). But for Luminal A and Luminal B patients, there was essentially no difference in OS between the two groups (Luminal A, P=0.726; Luminal B, P=0.49).

Establishment and validation of a nomogram

Univariate and multivariate Cox regression analyses were used to assess whether cytokine-related riskscore could be used as an independent prognostic factor. According to the result of the univariate Cox hazard model, stage (P=0.001), Luminal B subtype (P=0.01) and riskscore (P<0.001) were significantly related to OS. After considering these factors, multivariate Cox regression analysis showed that cytokine-related riskscore was an independent predictor of poor survival in the training set (HR=1.2, P=0.003 Table 2). Then in the validation set, through the univariate Cox regression analysis, only cytokine-related riskscore was proven to be independent risk factors to predict overall survival (HR=1.6, P=0.023 Table 3). After adding all factors to multivariate Cox regression, riskscore (HR=2.25, P=0.031 Table 3) was still an independent prognostic factor. Subsequently, a nomogram of tumor stage, molecular subtype and riskscore was created to predict the 5-year OS of BC patients in the training cohort. (Figure 3A) The C-index of this nomogram was 0.79. The nomogram scores of the training and validation sets were then calculated using the “formula_rd” function in the “nomogramFormula” package. To compare the predictive value of the nomogram with the riskscore and other available clinical characteristics, ROC curves for 5-year OS were generated in the training set and validation set respectively. As shown in Figure 3 the nomogram model had larger AUCs compared to other metrics. In the training set (Figure 3B), the 5-year AUCs for nomogram, riskscore, stage, and subtype were 0.78, 0.76, 0.61, and 0.44; in the validation set (Figure 3C), were 0.68, 0.68, 0.5, and 0.5.

TABLE 2

Table 2 Univariate and multivariate cox analyses in the training cohort.

TABLE 3

Table 3 Univariate and Multivariate Cox Analyses in the Validation Cohort.

FIGURE 3

Figure 3 Nomogram for predicting 5-year OS in the training set (A); ROC curves for 5-year OS based on the nomogram, riskscore, and other clinical traits in the Training Cohort (B) and Validation Cohort (C).

Correlation analysis between cytokines and blood indicators

As in the heatmap, the liver and kidney function, coagulation factors and routine blood test results of 84 breast cancer blood samples were retrospectively analyzed to explore the correlation between the two cytokines, IL-4 and TNF-α, the statistically significant results are marked with * (Figure 4). IL4 was negatively correlated with ALB (Spearman r = - 0.311, p <0.01) and positively correlated with β2-MG (Spearman r = 0.246, p <0.001), D-Dimer (Spearman r = 0.234, p <0.001) and ADA (Spearman r = 0.177, p <0.01). TNF-α was negatively correlated with LAP (Spearman r = - 0.067, p <0.001) and positively correlated with 5’NT (Spearman r = 0.118, p <0.001), D-Dimer (Spearman r = 0.259, p <0.001) and ADA (Spearman r = 0.14, p <0.01).

FIGURE 4

Figure 4 Correlation between cytokines and blood test results in breast cancer. P-value<0.05 is marked with *, P-value<0.01 is marked with **, P-value<0.001 is marked with ***.

Survival analysis based on Kaplan-Meier plotter database

The Kaplan-Meier plotter public mRNA expression dataset was applied to analyze the relationship between IL-4 and TNFA mRNA expression levels and prognosis. The OS of patients with low IL-4 (Figure 5A, HR=3.07, P=0.039) and TNFA (Figure 5B, HR=2.64, P=0.035) gene expression were significantly better than in patients with high gene expression divided by the median value.

FIGURE 5

Figure 5 The Overall survival curves of breast cancer patients according to IL-4 (A) and TNFA (B) gene expression levels in the public dataset.

Discussion

In this study, based on the detected levels of cytokines in serum, we screened out two cytokines by Cox analysis and “step” function to constructed a riskscore model related to breast cancer prognosis. In both training and validation cohorts, we found that the riskscore for IL-4 and TNF-α had a detrimental effect on overall survival in breast cancer. Immediately after, we established a nomogram containing tumor stage, molecular subtype and riskscore, and this nomogram has better prediction effect compared with cytokine-related riskscore.

IL-4 belongs to the interleukin family and interleukins are small molecular proteins secreted mainly by CD3 + and CD4 + T lymphocytes, which mediate the necessary interactions for the progress of cancer cells (14). Il-4 is called “prototype immunoregulatory cytokine”. Like many cytokines, it can affect a variety of target cells in various ways and plays an important role in regulating antibody production, hematopoiesis, inflammation and the development of efficacy T cells response (15). Cytokines can regulate key interactions between immune and non-immune cells in the tumor microenvironment (16).A study proved that IL-4 acts as a survival factor for cancer cells including breast cancer by increasing the expression of anti-apoptotic molecules, making them resistant to induction of apoptosis (17). A study in 1992 concluded that IL-4 is a non-autocrine inhibitor of breast cancer and can inhibit tumor proliferation (18). However, there are some discrepancies between the available studies on the correlation between serum cytokines and breast cancer prognosis. Our findings differ from these results. We obtained the cutoff value of IL-4 by calculating the Youden index and divided the patients into high and low expression groups. The Kaplan-Meier curve showed that patients in the high serum level group had a shorter survival time. The same results were verified in the Kaplan-Meier plotter database on mRNA level. As well as being a risk factor for breast cancer patients, DeNardo DG found that in murine breast cancer models, IL-4 produced by tumor infiltrating CD4 + T cells stimulate M1 to M2 transition of TAMs (Tumor-associated macrophages), thus supporting lung metastasis of breast cancer cells (19).

Tumor necrosis factor α belongs to tumor necrosis factor/tumor necrosis factor receptor (TNF/TNFR) superfamily cytokines and belongs to type II transmembrane protein (20). TNF-α is one of the most important pro-inflammatory cytokines in breast cancer TME, which is secreted by stromal cells (mainly M1 TAM) and cancer cells themselves. TNF-α plays a dual role in cancer and can enhance or inhibit tumor progression depending on the discovery of its specific cellular environment (21).High serum TNF-α remains an independent prognostic factor for poor overall survival in a case-control study of primary breast cancer patients (22). Which was similar to our study results. In our research, after adding age, stage and molecular type to the Cox analysis, the riskscore composed of IL-4 and TNF-α was still identified as a risk factor for OS in both the training and validation sets.

Cytokines have been reported to play a key role in the immune response to cancer cells. Several studies have been conducted to explore the correlation between serum cytokine levels and breast cancer immunity. At present, some studies have reported that cytokines can be used as a prognostic factor in breast cancer treatment. A study of the response to neoadjuvant chemotherapy in locally advanced breast cancer showed that serum TNF-α levels were significantly lower both in partial and complete response compared to pre-treatment (p < 0.01), there was also a significant difference in TNF-α levels between complete and partial response after treatment (23). In univariate and multivariate analyses, low TGF-β1 and TNF-α and high IFN-γ were associated with significantly lower risk ratios for disease progression, and serum TGF-β1, TNF-α, and IFN-γ levels may be potential prognostic biomarkers for breast cancer patients with bone metastases treated with taxane and zoledronic acid chemotherapy (24). In immunotherapy, TNF-α induced mucin 4 expression in HER2-positive breast cancer can cause resistance to trastuzumab (25). Increased IL4 levels found to be associated with better survival in small cell lung cancer patients receiving ipilimumab immunotherapy (26).

In clinical practice, we can combine cytokines and blood indicator analysis to find the cytokine that is related to the blood indicator, by detecting the expression level of this cytokine to prompt us to predict the outcome of patients. We have collected several blood tests results from patients record and correlated them with the two cytokines. In our study, the results showed that ALB was negatively correlated with IL-4, while β2-MG, D-Dimer and ADA were significantly positively correlated with IL-4. As for TNF-α, it was positively associated with D-Dimmer, ADA and 5’ NT and negatively associated with LAP. Albumin is produced by the liver and has been used to assess severity of disease, disease progression and prognosis. Christopher G Lis et al. found that albumin levels<3.5g/dL had a threefold increased risk for death, and baseline serum albumin levels were a powerful prognostic variable (27).A study showed that serum albumin, lactate dehydrogenase and total bilirubin were significantly associated with 5-year OS in multivariate Cox analysis in breast cancer (11). However, there are few studies on the prognostic value of β2-mg, D-dimmer, 5’ NT, ADA and LAP in breast cancer at present. D-dimer is a biomarker used to indicate clotting and fibrinolytic activity and is more widely used in predicting the risk of thrombosis. There were two studies found significant reductions in serum 5’-NT levels after mastectomy, and the 1995 study also found decreases in serum ADA levels (28, 29). Another research suggested that ADA activity in breast tissue can be an indicator for the differential diagnosis of benign and malignant breast disease. Due to the age of these studies and the lack of more verification, the conclusions need to be more supported in future research. In the next study, we will further explore the effect of cytokines combined with blood markers on the prognosis of breast cancer.

For the nomogram we build, age was not included, because it has not been verified as an independent prognostic factor for OS in BC patients in multivariate Cox analysis. However, in the subgroup analyses we could see that the prognosis of patients with high-risk was poorer and significantly different in both young and old populations, suggesting that there is also a correlation between age and the riskscore. In the subgroup analysis, we can observe that for patients with tumor stage III/IV, the low-risk group has a better prognosis. And among the different molecular subtypes, TNBC and HER2-positive patients with high-risk had a shorter survival. This suggests that cytokine-related riskscore correlates with patient age, tumor stage and molecular subtype, but it has less predictive value for patients with early stage tumors and Luminal type.

One limitation of our study is that we did not include healthy controls and lacked data to compare the differences in cytokine levels between patients and normal individuals. Several studies have found that there were statistical differences in cytokine levels between cancer patients and healthy controls. Significantly elevated levels of IL-4, IL-8 and IL-9 were observed in colorectal cancer (30, 31). And IL-4 was also found higher in breast cancer patients (32). Ultimately, due to the small sample size, our findings need to be validated in a larger cohort and multicenter study.

Conclusion

Taken together, our study explored the role of cytokine-related genes in the prognosis of BC patients. We analyzed and validated a cytokine-related riskscore, and developed a prognostic nomogram combining tumor stage, molecular subtype and riskscore for clinical OS prediction. Then analyzed the correlation between cytokines and blood test indicators. It is hoped that this nomogram will provide individualized prediction of survival for BC patients.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethical Committee of the Tianjin Medical University Cancer Institute and Hospital. The patients/participants provided their written informed consent to participate in this study.

Author contributions

YZ: Conceptualization, data collection, formal analysis, and methodology and writing. YZ and YH: data analysis and revised the manuscript. YH and CC: Data curation and experiments. JZ, XY, and YL: Collection and assembly of data. WZ and Y-ZC: Conceptualization and manuscript review and editing. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National Natural Science Foundation of China (grant number 81472183), Natural Science Foundation of Tianjin (18JCYBJC91600), Shenzhen Chipscreen Biosciences project, Tianjin Key Medical Discipline (Specialty) Construction Project and Tianjin Medical University Cancer Hospital “14th Five-Year” Peak Discipline Support Program Project.

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 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/fonc.2023.1146463/full#supplementary-material

Abbreviations

OS, overall survival; C-indexes, corresponding concordance indexes; DFS, disease-free survival; TNBC, triple negative breast cancer; NEUT, neutrophil granulocyte; EO, eosinophils granulocyte; BASO, basophil granulocyte; ALB, albumin; GLO, globulin; A/G, albumin/globulin ratio; ChE, cholinesterase; ADA, adenosine deaminase; LDH, lactate dehydrogenase; GPDA, glycylyl proline dipeptide aminopeptidase; 5’NT, 5’nucleotidase; LAP, leucineaminopeptidase; β2-MG, beta2-microglobulin; ROC, receiver operating characteristic; BC, breast cancer; AUC, area under the curve.

References

1. Planeix F, Siraj M-A, Bidard F-C, Robin B, Pichon C, Sastre-Garau X, et al. Endothelial follicle-stimulating hormone receptor expression in invasive breast cancer and vascular remodeling at tumor periphery. J Exp Clin Cancer Res CR. (2015) 34:12. doi: 10.1186/s13046-015-0128-7

CrossRef Full Text | Google Scholar

2. Nicolini A, Carpi A, Rossi G. Cytokines in breast cancer. Cytokine Growth Factor Rev (2006) 17:325–37. doi: 10.1016/j.cytogfr.2006.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Wang H, Yang X. Association between serum cytokines and progression of breast cancer in Chinese population. Medicine (2017) 96:e8840. doi: 10.1097/MD.0000000000008840

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Lin S, Gan Z, Han K, Yao Y, Min D. Interleukin-6 as a prognostic marker for breast cancer: a meta-analysis. Tumori (2015) 101:535–41. doi: 10.5301/tj.5000357

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ahmad N, Ammar A, Storr SJ, Green AR, Rakha E, Ellis IO, et al. IL-6 and IL-10 are associated with good prognosis in early stage invasive breast cancer patients. Cancer Immunol Immunother CII (2018) 67:537–49. doi: 10.1007/s00262-017-2106-8

CrossRef Full Text | Google Scholar

6. Salgado R, Junius S, Benoy I, Van Dam P, Vermeulen P, Van Marck E, et al. Circulating interleukin-6 predicts survival in patients with metastatic breast cancer. Int J Cancer (2003) 103:642–6. doi: 10.1002/ijc.10833

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Hartman ZC, Poage GM, den Hollander P, Tsimelzon A, Hill J, Panupinthu N, et al. Growth of triple-negative breast cancer cells relies upon coordinate autocrine expression of the proinflammatory cytokines IL-6 and IL-8. Cancer Res (2013) 73:3470–80. doi: 10.1158/0008-5472.CAN-12-4524-T

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Benoy IH, Salgado R, Van Dam P, Geboers K, Van Marck E, Scharpé S, et al. Increased serum interleukin-8 in patients with early and metastatic breast cancer correlates with early dissemination and survival. Clin Cancer Res an Off J Am Assoc For Cancer Res (2004) 10:7157–62. doi: 10.1158/1078-0432.CCR-04-0812

CrossRef Full Text | Google Scholar

9. Chen W-C, Lai Y-H, Chen H-Y, Guo H-R, Su I-J, Chen HHW. Interleukin-17-producing cell infiltration in the breast cancer tumour microenvironment is a poor prognostic factor. Histopathology (2013) 63:225–33. doi: 10.1111/his.12156

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Grivennikov SI, Greten FR, Karin M. Immunity, inflammation, and cancer. Cell (2010) 140:883–99. doi: 10.1016/j.cell.2010.01.025

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Liu X, Meng QH, Ye Y, Hildebrandt MAT, Gu J, Wu X. Prognostic significance of pretreatment serum levels of albumin, LDH and total bilirubin in patients with non-metastatic breast cancer. Carcinogenesis (2015) 36:243–8. doi: 10.1093/carcin/bgu247

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Liu D, Wang D, Wu C, Zhang L, Mei Q, Hu G, et al. Prognostic significance of serum lactate dehydrogenase in patients with breast cancer: a meta-analysis. Cancer Manage Res (2019) 11:3611–9. doi: 10.2147/CMAR.S199260

CrossRef Full Text | Google Scholar

13. Győrffy B. Survival analysis across the entire transcriptome identifies biomarkers with the highest prognostic power in breast cancer. Comput Struct Biotechnol J (2021) 19:4101–9. doi: 10.1016/j.csbj.2021.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Fasoulakis Z, Kolios G, Papamanolis V, Kontomanolis EN. Interleukins associated with breast cancer. Cureus (2018) 10:e3549. doi: 10.7759/cureus.3549

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Brown MA, Hural J. Functions of IL-4 and control of its expression. Crit Rev In Immunol (1997) 17:1–32. doi: 10.1615/critrevimmunol.v17.i1.10

CrossRef Full Text | Google Scholar

16. Briukhovetska D, Dörr J, Endres S, Libby P, Dinarello CA, Kobold S. Interleukins in cancer: from biology to therapy. Nat Rev Cancer (2021) 21:481–99. doi: 10.1038/s41568-021-00363-z

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Todaro M, Lombardo Y, Francipane MG, Alea MP, Cammareri P, Iovino F, et al. Apoptosis resistance in epithelial tumors is mediated by tumor-cell-derived interleukin-4. Cell Death Differ (2008) 15:762–72. doi: 10.1038/sj.cdd.4402305

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Toi M, Bicknell R, Harris AL. Inhibition of colon and breast carcinoma cell growth by interleukin-4. Cancer Res (1992) 52:275–9.

PubMed Abstract | Google Scholar

19. DeNardo DG, Barreto JB, Andreu P, Vasquez L, Tawfik D, Kolhatkar N, et al. CD4(+) T cells regulate pulmonary metastasis of mammary carcinomas by enhancing protumor properties of macrophages. Cancer Cell (2009) 16:91–102. doi: 10.1016/j.ccr.2009.06.018

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Cruceriu D, Baldasici O, Balacescu O, Berindan-Neagoe I. The dual role of tumor necrosis factor-alpha (TNF-α) in breast cancer: molecular insights and therapeutic approaches. Cell Oncol (Dordrecht) (2020) 43:1–18. doi: 10.1007/s13402-019-00489-1

CrossRef Full Text | Google Scholar

21. Balkwill F. Tumour necrosis factor and cancer. Nat Rev Cancer (2009) 9:361–71. doi: 10.1038/nrc2628

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Papadopoulou E, Tripsianis G, Anagnostopoulos K, Tentes I, Kakolyris S, Galazios G, et al. Significance of serum tumor necrosis factor-alpha and its combination with HER-2 codon 655 polymorphism in the diagnosis and prognosis of breast cancer. Int J Biol Markers (2010) 25:126–35. doi: 10.1177/172460081002500302

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Berberoglu U, Yildirim E, Celen O. Serum levels of tumor necrosis factor alpha correlate with response to neoadjuvant chemotherapy in locally advanced breast cancer. Int J Biol Markers (2004) 19:130–4. doi: 10.1177/172460080401900207

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kim JW, Lee S, Kim HS, Choi YJ, Yoo J, Park KU, et al. Prognostic effects of cytokine levels on patients treated with taxane and zoledronic acid for metastatic breast cancer in bone (BEAT-ZO) (KCSG BR 10-13). Cytokine (2021) 142:155487. doi: 10.1016/j.cyto.2021.155487

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Mercogliano MF, De Martino M, Venturutti L, Rivas MA, Proietti CJ, Inurrigarro G, et al. TNFα-induced mucin 4 expression elicits trastuzumab resistance in HER2-positive breast cancer. Clin Cancer Res an Off J Am Assoc For Cancer Res (2017) 23:636–48. doi: 10.1158/1078-0432.CCR-16-0970

CrossRef Full Text | Google Scholar

26. Hardy-Werbin M, Rocha P, Arpi O, Taus Á, Nonell L, Durán X, et al. Serum cytokine levels as predictive biomarkers of benefit from ipilimumab in small cell lung cancer. Oncoimmunology (2019) 8:e1593810. doi: 10.1080/2162402X.2019.1593810

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lis CG, Grutsch JF, Vashi PG, Lammersfeld CA. Is serum albumin an independent predictor of survival in patients with breast cancer? JPEN J Parenteral Enteral Nutr (2003) 27:10–5. doi: 10.1177/014860710302700110

CrossRef Full Text | Google Scholar

28. Dao TL, Ip C, Patel J. Serum sialyltransferase and 5'-nucleotidase as reliable biomarkers in women with breast cancer. J Natl Cancer Institute (1980) 65:529–34.

Google Scholar

29. Walia M, Mahajan M, Singh K. Serum adenosine deaminase, 5'-nucleotidase & alkaline phosphatase in breast cancer patients. Indian J Med Res (1995) 101:247–9.

PubMed Abstract | Google Scholar

30. Yamaguchi M, Okamura S, Yamaji T, Iwasaki M, Tsugane S, Shetty V, et al. Plasma cytokine levels and the presence of colorectal cancer. PloS One (2019) 14:e0213602. doi: 10.1371/journal.pone.0213602

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Czajka-Francuz P, Francuz T, Cisoń-Jurek S, Czajka A, Fajkis M, Szymczak B, et al. Serum cytokine profile as a potential prognostic tool in colorectal cancer patients - one center study. Rep Pract Oncol Radiotherapy J Greatpoland Cancer Center In Poznan Polish Soc Radiat Oncol (2020) 25:867–75. doi: 10.1016/j.rpor.2020.08.004

CrossRef Full Text | Google Scholar

32. Zhu X, Du L, Feng J, Ling Y, Xu S. Clinicopathological and prognostic significance of serum cytokine levels in breast cancer. Clin Laboratory (2014) 60:1145–51. doi: 10.7754/clin.lab.2013.130738

CrossRef Full Text | Google Scholar