- Research article
- Open access
- Published:
Survival analysis in gastric cancer: a multi-center study among Iranian patients
BMC Surgery volumeĀ 20, ArticleĀ number:Ā 152 (2020)
Abstract
Background
Gastric cancer (GC) has been considered as the 5th most common type of cancer and the third leading cause of cancer-associated death worldwide. The aim of this historical cohort study was to evaluate the survival predictors for all patients with GC using the Cox proportional hazards, extended Cox, and gamma-frailty models.
Methods
This historical cohort study was performed according to documents of 1695 individuals having GC referred to three medical centers in Iran from 2001 to 2018. First, most significant prognostic risk factors on survival were selected, Cox proportional hazards, extended Cox, gamma-frailty models were applied to evaluate the effects of the risk factors, and then these models were compared with the Akaike information criterion.
Results
The age of patients, body mass index (BMI), tumor size, type of treatment and grade of the tumor increased the hazard rate (HR) of GC patients in both the Cox and frailty models (Pā<ā0.05). Also, the size of the tumor and BMI were considered as time-varying variables in the extended Cox model. Moreover, the frailty model showed that there is at least an unknown factor, genetic or environmental factors, in the model that is not measured (Pā<ā0.05).
Conclusions
Some prognostic factors, including age, tumor size, the grade of the tumor, type of treatment and BMI, were regarded as indispensable predictors in patients of GC. Frailty model revealed that there are unknown or latent factors, genetic and environmental factors, resulting in the biased estimates of the regression coefficients.
Background
GC is the 5th leading cause of cancer-related death in spite of its global decrease in incidence and mortality according to GLOBOCAN 2018 data [1]. Despite advances in diagnostic, therapeutic, and screening methods, the mortality rate has not significantly decreased worldwide. In the Iranian population, the prevalence of this cancer is also increasing as it threatens the health of the population [2]. The survival rate of GC is rather low and often the tumor is not diagnosed until an advanced stage that the cause is related to clinicopathological factors [3]. Indeed, malignancy is frequently identified by a variable but poor overall prognosis, particularly in the late clinical stages. The clinicopathological features of GC patients are very imperative in choosing the right therapeutic strategy that can improve patient survival. In addition, diagnosis and treatment of GC considerably depend on prognostic factors and variations of survival over time [4, 5].
Since identifying the prognosis risk factors for GC patients is extensively important, researchers are interested in the survival time until the occurrence of an event in epidemiological and other data [6]. The most reasonable condition in Cox regression is the proportional hazard (CPH) assumption, which is applied in a short-term follow-up. If the proportional assumption does not hold, the results from a CPH model are misleading, and alternative modelling strategies should be carried out [7]. In a common phenomenon in clinical research, time-varying covariance occurs when a given covariate changes over time during the follow-up period that is called an extended Cox model. Indeed, the main characteristic of data with time-dependent covariates is the survivor function for any individual depending on time and the baseline hazard function [8]. The frailty model, the random component, has been designed to account for variability, and it has been used when there is at least one unaccounted predictor in the model [9]. This model assumes that events (e.g., death) happen earlier for individuals who are more frailty. Some factors, including hereditary, genetic characteristics, growth, and living environment, are effective in the caused differences among the patients. When the assumption of proportionality does not hold, applying the CPH regression leads to a biased estimation and underestimation of variance of the parameters [10].
Other studies have been done on survival analysis such as parametric, artificial neural network, Bayesian and parametric, multi-state, and Cure models in GC [11,12,13]. Viduz et al. surveyed the frailty multi-state model on advanced GC data from the Agamenon National Cancer Registry [14]. Lu et al. performed Cox regression model and log-rank test in patients older than 80āyears who underwent radical gastrectomy for primary GC from 2000 to 2012 [15]. Ghadimi et al. analysed the survival rate of the gastrointestinal patients by parametric models such as log-normal, log-logistic, Weibull, and exponential model using with and without frailty, then Akaike information criterion (AIC) was regarded to evaluate the models [16]. Faradmal et al. applied Cox and Frailty models in Breast cancer data and then compared them with the Concordance index [13].
To the best of our knowledge, this is the first multi-center study that investigates the main prognostic factors of GC in Iranian through applying several survival analysis models. The multi-center study is designed to survey multiple cities in Iran. The sample size (nā=ā1695) is large enough to detect a wide range of associations with adequate statistical power.
We aimed to evaluate the effect of important variables on the survival rate of GC patients who registered at three centers in Iran during 2001ā2018, using Cox regression and two semi-parametric models.
Methods
Patient characteristics
In this study, we included 1695 patients who were diagnosed to have GC and were registered to three separated medical centers in Iran, Rasoul Akram hospital (2013ā2018), Taleghani hospital (2003ā2007), and Fars province in southern Iran (2001ā2006) during 2001ā2018. The project was approved in Ethic Committee of Iran University of Medical Sciences (ethical code: IR.IUMS.REC. 1397.481).
-
The right-censored data at Rasoul Akram hospital of Tehran represented a historical cohort of 346 GC patients, who registered from September 2013 to November 2018. Demographic and clinical characteristics of patients were obtained by checklists of patientās records.
-
The right-censored data of Taleghani hospital revealed a retrospective review of 746 GC patients who enrolled in the study from February 2003 until January 2007. The patientās information was gathered by checklists of patientās records [17].
-
The data of the cancer registry of Fars province center demonstrated a historical cohort study of 603 GC patients from March 2001 to March 2006. Prognostic factors of GC and patientsā vital status were collected in March 2006 [18].
The outcome variable was considered as time (months) elapsed since the cancer diagnosis until death. Some important clinical variables containing tumor size, number of involved lymph nodes, distant metastasis, histology, type of treatment, and demographic variables such as age, gender, marital status, education, BMI, and smoking situation were included in the aforementioned models.
Statistical analysis
Kaplan-Meier, named as the product limit estimator, was used to estimate survival function. In the first step, the Univariate CPH model was performed to find the important factors of GC, chosen variables with Pā<ā0.2 in the Univariate analysis were subjected to multivariable regression analysis with Pā<ā0.05. The Cox model is a very useful approach to survival analysis. On the other hand, when the assumption of proportionality does not satisfy, the outputs might be misleading, and then other different models should be used [7]. In the time-varying variables, the effects of some predictors depend on the time that is called time-dependent variables, an extended version of the CPH model [19]. The proportional hazard (PH) assumption was assessed using Schoenfeld residuals, and then an extended CPH model was fitted to the data [7]. Since several substantial factors such as genetic or environmental factors were not reported in these data, there was evidence of unmeasured heterogeneity among the patients, so a frailty model was applied [20]. In the frailty model, an unobserved multiplicative effect was considered on hazard function by presuming a g(Ī±) distribution with the unit mean and unknown variance of Īø. In the third step, gamma distribution was performed as a frailty component in the Cox model [10]. Finally, the evaluation of models was performed based on the AIC. The significance level for the statistical test was 0.05. The Stata-13 and R-3.2.2 were applied for all statistical analysis.
Results
In this study, medical files of 1695 GC patients were retrospectively reviewed. 949 (56%) were male, and 746 (44%) were female. The mean age of patients at diagnosis was 60.28āĀ±ā12.93, and the follow-up time was 18.79āĀ±ā16.67āmonths. The overall median survival rate of 1695 patients was 13.2āmonths. Eight hundred four (47.4%) of patients were censored, and 891 (52.6%) of patients died at the end of follow-up. The characteristics and pathological features of all the GC patients are presented in TableĀ 1.
Based on KaplanāMeier estimates, the 5-year survival rate was assessed 11.3%, while one-year and three-year survival rates were 63.5 and 37.5%, respectively (Fig.Ā 1).
The test based on Schoenfeld residuals revealed that some variables, including the size of the tumor, the grade of the tumor, and BMI (Pā<ā0.05), were violated from the PH assumption (Global Pā=ā0.04). The results of the CPH were provided in TableĀ 2. The table showed that age, sex, smoking, BMI, type of tumor, involved lymph node, metastasis, type of treatment, tumor size, the grade of the tumor are significant in the univariate Cox model (Pā<ā0.15). Moreover, the output of the multivariate Cox model revealed that age, BMI, type of treatment, tumor size, and grade of the tumor are statistically significant (Pā<ā0.05).
A HR of tumor size categories indicated that T2, T3 (larger sizes) had worse prognoses compare to T1 (small size) (HRā=ā1.5; HRā=ā3.6, Pā<ā0.001). Moreover, the HR of metastasis was 0.3% more than non-metastasis (HRā=ā1.3; Pā<ā0.001). In addition, the HR of death in patients with N2 (N3) was 40% (27%) times more than those with N1. Additionally, the HR of treatment type in surgery plus chemotherapy was 0.21% less than chemotherapy plus radiotherapy that seems to be significant (HRā=ā0.79, Pā<ā0.05). Based on the Multivariate CPH model, the HR of death of patients with BMIā<ā18 (underweight) is 2.45 times more than those with overweight (HRā=ā2.45, Pā<ā0.001). Also, the AIC is āā1620.36 in the CPH model. TableĀ 3 showed the results in evaluating PH to apply the extended Cox model and frailty model.
Since the PH assumption was not met for tumor size, age and BMI, CPH model no longer satisfies the PH assumption and therefore an extended Cox model was performed. Furthermore, the Global test demonstrated CPH did not seem to be suitable because of non-proportional hazards (Pā=ā0.04). The result of the extended Cox model revealed that the size of the tumor and BMI are time-varying variables (Pā<ā0.05 and the AIC is āā1452.07 in this model.
Furthermore, Table 3 indicates the significant variables in the frailty model, including age, tumor size, the grade of the tumor, and BMI (Pā<ā0.05). The results showed that the variance of frailty was significantly greater than zero (Īøā=ā1.67; Pā<ā0.05). According to our finding, there were latent factors which affect the hazard of death. The AIC in the frailty model is āā412.72, which is larger than āā1620.36 andāāā1452.07 in CPH and extended Cox models. These results of comparing AIC in three models indicate that the best models are frailty, extended Cox, and CPH, respectively.
Discussion
Survival analysis have mostly performed using common statistical methods such as CPH that have shortcomings [2, 21,22,23]. However, this is the first multi-center study, comparing different survival models to identify prognostic risk factors in GC patients of developing countries. GC is a worldwide cause of cancer death with a low 5-year survival rate in Iran. A number of factors have been identified as predictive prognosis factors until now [24,25,26]. In the recent study, the CPH, extended Cox, gamma frailty Cox models were fitted to determine the adjusted hazard of GC patients who underwent treatments, surgery, chemotherapy, and radiotherapy. The significant variables in the CPH model were age at diagnosis, BMI, tumor size, type of treatment and grade of the tumor. Additionally, time-dependent variables, the size of the tumor, and BMI with time function g(t)ā=āt were considered in the extended Cox model. Furthermore, the frailty model was chosen as the best model and demonstrated that there are latent factors that affect the hazard of death. The results revealed that age, tumor size, the grade of the tumor, type of treatment and BMI have a direct effect on the hazard of death in GC patients.
The results of CPH and frailty models showed a significant relationship between age over 60āyears at the time of surgery and the 5-year survival of patients with GC. Similar to our results, previous studies have reported the relationship between age and the 5-year survival [24, 26]. The main strength of these investigations was the large sample size that was consistent in our study. Whilst a contradictory result of a study indicated that there is not any relationship between older age and patient survival, which could be because of differences in the sample size [25]. A retrospective study revealed how patients over the age of 70āyears differ from younger patients in postoperative courses with a focus on the frequency of surgical and medical [27].
The findings of our study showed that both sizes of the tumor and tumor grade in patients at the time of diagnosis in three models were significantly associated with patientsā survival. The size of tumor and tumor grade have been identified as a risk factor for the survival of GC patients, which are compatible with our study [24, 26] while the result of Nasseri et al. demonstrated that there is no relationship between the grade of tumor and GC [28].
The main result of our study showed a significant relationship between the survival of GC patients with BMI <ā18, which was consistent with the findings of BMI in the study of Liu et al. [29]. In their study, a 320 cohort study was conducted to survey the effect of BMI and recreational physical activity on GC risk. It can be illustrated due to anorexia, weakness and poor health status of under-weight patients, therefore the variable can be regarded as indispensable prognostic factor. In general, studies of BMI and GC have been restricted. Furthermore, combined surgical treatment and chemotherapy were related to higher survival rate that the result was consistent with several previous studies [30,31,32,33]. A meta-analysis was performed to survey the effect of treatment type on GC [33]. Moreover, combination therapy surgery as well as chemotherapy was significant in some cancers such as Hepatocellular carcinoma [34]. However, an inconsistent study presented the survival analysis of GC patients with incomplete data that treatment type was not significant [4].
The frailty model is performed to explain the random variation of the survival function that may exist due to some unobserved genetic prognostic factors such as genetic and other environmental factors. These results are consistent with the results of many studies in this field [16, 35].
Based on the AIC, the frailty model is the best alternative model for the Cox proportional hazard model. This issue is consistent with most studies conducted on GC patients [10, 15, 28].
Strengths and limitations
The main strength of the present investigation was the multi-center study with a large sample size as well as a lack of missing in the data. The main limitation of this study was the short- term follow-up period. Further, studies with longer follow-up periods may provide more determining evidence regarding the Survival predictors in GC patients.
Conclusion
The results of this study indicate that the ageā>ā60āyears, tumor size, the grade of tumor, type of treatment and BMIā<ā18ākg/m2 are the main prognostic factors in the survival rate of GC patients. In fact, they are reducing the survival rate of GC patients. Also, based on our findings from the frailty model, we might conclude that employing a more intricate statistical model that regards the significant role of latent variables on hazard ratio, including unobserved genetic or environmental factors, would expand the importance of the more analyses.
Availability of data and materials
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- HR:
-
Hazard rate
- BMI:
-
Body mass index
- GC:
-
Gastric cancer
- CPH:
-
Cox proportional hazards model
- PH:
-
Proportional hazards
- AIC:
-
Akaike information criterion
References
Ferlay J, Ervik M, Lam F, Colombet M, Mery L, PiƱeros M, Znaor A, Soerjomataram I, Bray F. Global cancer Observatory: cancer today. Lyon: International agency for research on cancer. Cancer Today; 2018.
Jenabi E, Saatchi M, Khazaei S, Mansori K, Ayubi E, Soheylizad M, Khazaei L, Nematollahi S, Beigi AM. National distribution of stomach cancer incidence in Iran: a population-based study. Adv Hum Biol. 2019;9(1):89.
Zarea K, Beiranvand S, Ghanbari S, Tuvesson H. Incidence of gastrointestinal cancers in Iran: a systematic review. Jundishapur J Chronic Dis Care. 2017;6(1):e37224.
Yaprak G, Tataroglu D, Dogan B, Pekyurek M. Prognostic factors for survival in patients with gastric cancer: single-Centre experience. Northern Clin Istanb. 2020;7(2):146.
Yoshida N, Doyama H, Yano T, Horimatsu T, Uedo N, Yamamoto Y, Kakushima N, Kanzaki H, Hori S, Yao K: Early gastric cancer detection in high-risk patients: a multicentre randomised controlled trial on the effect of second-generation narrow band imaging. Gut. 2020;1ā9.
Kleinbaum DG, Klein M. Survival analysis, vol. 3: Springer; 2010.
Klein JP, Van Houwelingen HC, Ibrahim JG, Scheike TH. Handbook of survival analysis. CRC Press; 2016.
Liu M, Li L, Yu W, Chen J, Xiong W, Chen S, Yu L. Marriage is a dependent risk factor for mortality of colon adenocarcinoma without a time-varying effect. Oncotarget. 2017;8(12):20056.
Moore DF. Applied survival analysis using R: Springer; 2016.
Hanagal DD. Frailty Models in Public Health. In: Handbook of Statistics, vol. 37: Elsevier; 2017. p. 209ā47.
Amoori N, Mahdavi S, Enayatrad M. Epidemiology and trend of stomach cancer mortality in Iran. Int J Epidemiol Res. 2016;3(3):268ā75.
Charati JY, Janbabaei G, Alipour N, Mohammadi S, Gholiabad SG, Fendereski A. Survival prediction of gastric cancer patients by artificial neural network model. Gastroenterol Hepatol Bed Bench. 2018;11(2):110.
Faradmal J, Talebi A, Rezaianzadeh A, Mahjub H. Survival analysis of breast cancer patients using cox and frailty models. J Res Health Sci. 2012;12(2):127ā30.
ViĆŗdez A, Carmona-Bayonas A, Gallego J, Lacalle A, HernĆ”ndez R, Cano J, MacĆas I, Custodio A, de Castro EM, SĆ”nchez A. Optimal duration of first-line chemotherapy for advanced gastric cancer: data from the AGAMENON registry. Clin Transl Oncol. 2019:1ā17.
Lu J, Cao L-L, Zheng C-H, Li P, Xie J-W, Wang J-B, Lin J-X, Chen Q-Y, Lin M, Tu R-H. The preoperative frailty versus inflammation-based prognostic score: which is better as an objective predictor for gastric cancer patients 80 years and older? Ann Surg Oncol. 2017;24(3):754ā62.
Ghadimi M, Mahmoodi M, Mohammad K, Zeraati H, Rasouli M, Sheikhfathollahi M. Family history of the cancer on the survival of the patients with gastrointestinal cancer in northern Iran, using frailty models. BMC Gastroenterol. 2011;11(1):104.
Pourhoseingholi MA, Hajizadeh E, Moghimi Dehkordi B, Safaee A, Abadi A, Zali MR. Comparing cox regression and parametric models for survival of patients with gastric carcinoma. Asian Pac J Cancer Prev. 2007;8(3):412.
Moghimi-Dehkordi B, Safaee A, Pourhoseingholi MA, Fatemi R, Tabeie Z, Zali MR. Statistical comparison of survival models for analysis of cancer data. Asian Pac J Cancer Prev. 2008;9(3):417ā20.
Perperoglou A, Keramopoullos A, van Houwelingen HC. Approaches in modelling long-term survival: an application to breast cancer. Stat Med. 2007;26(13):2666ā85.
Pelagia I. Variable Selection of Fixed Effects and Frailties for Cox Proportional Hazard Frailty Models and Competing Risks Frailty Models: The University of Manchester (United Kingdom); 2016.
Baghestani AR, Pourhoseingholi MA, Ashtari S, Najafimehr H, Busani L, Zali MR. Trend of gastric Cancer after Bayesian correction of misclassification error in neighboring provinces of Iran. Galen Med J. 2019;8:1223.
Ghaffari HR, Yunesian M, Nabizadeh R, Nasseri S, Sadjadi A, Pourfarzi F, Poustchi H, Eshraghian A. Environmental etiology of gastric cancer in Iran: a systematic review focusing on drinking water, soil, food, radiation, and geographical conditions. Environ Sci Pollut Res. 2019;26(11):10487ā95.
Rastaghi S, Jafari-Koshki T, Mahaki B, Bashiri Y, Mehrabani K, Soleimani A. Trends and risk factors of gastric cancer in Iran (2005ā2010). Int J Prev Med. 2019;10.
Katai H, Ishikawa T, Akazawa K, Isobe Y, Miyashiro I, Oda I, Tsujitani S, Ono H, Tanabe S, Fukagawa T. Five-year survival analysis of surgically resected gastric cancer cases in Japan: a retrospective analysis of more than 100,000 patients from the nationwide registry of the Japanese gastric Cancer association (2001ā2007). Gastric Cancer. 2018;21(1):144ā54.
Morgan R, Cassidy M, DeGeus SW, Tseng J, McAneny D, Sachs T. Presentation and survival of gastric Cancer patients at an urban academic safety-net hospital. J Gastrointest Surg. 2019;23(2):239ā46.
van Putten M, Nelen S, Lemmens V, Stoot J, Hartgrink H, Gisbertz SS, Spillenaar Bilgen EJ, Heisterkamp J, Verhoeven R, Nieuwenhuijzen G. Overall survival before and after centralization of gastric cancer surgery in the Netherlands. Br J Surg. 2018;105(13):1807ā15.
Henrik Nienhueser RK, Sisic L, Blank S, Strowitzk MJ, Bruckner T, Jager D, Weichert W, Alexis Ulrich MB, Ott K, Schmidt T. Surgery of Gastric Cancer and Esophageal Cancer: Does Age Matter? J Surg Oncol. 2015;112:387ā95.
Naseri P, Baghestani AR, Momenyan N, Akbari ME. Application of a mixture cure fraction model based on the generalized modified weibull distribution for analyzing survival of patients with breast cancer. Int J Cancer Manage. 2018;11(5).
Liu BZ, Tao L, Chen YZ, Li XZ, Dong YL, Ma YJ, Li SG, Li F, Zhang WJ. Preoperative body mass index, blood albumin and triglycerides predict survival for patients with gastric cancer. PLoS One. 2016;11(6).
Dongyun Yang AH, Lenz C, Togawa K, Lenz F, Lurje G, Pohl A, Winder T, Ning Y, Groshen S, Lenz H-J. Survival of metastatic gastric cancer: significance of age, sex and race/ethnicity. J Gastrointest Oncol. 2011;2(2):77ā84.
Zhang X-F, Huang C-M, Lu H-S, Wu X-Y, Wang C, Guang G-X, Zhang J-Z, Zheng C-H. Surgical treatment and prognosis of gastric cancer in 2613 patients. World J Gastroenterol: WJG. 2004;10(23):3405.
Moghimbeigi A, Tapak L, Roshanaei G, Mahjub H. Survival analysis of gastric cancer patients with incomplete data. J Gastric Cancer. 2014;14(4):259ā65.
Orditura M, Galizia G, Sforza V, Gambardella V, Fabozzi A, Laterza MM, Andreozzi F, Ventriglia J, Savastano B, Mabilia A. Treatment of gastric cancer. World J Gastroenterol: WJG. 2014;20(7):1635.
Sarveazad A, Agah S, Babahajian A, Amini N, Bahardoust M. Predictors of 5 year survival rate in hepatocellular carcinoma patients. J Res Med Sci. 2019;24.
Nikaeen R, Khalilian A, Bahrampour A. Determining the effective factors on gastric cancer using frailty model in South-East and North of Iran. 2017;5(3):35-48.
Acknowledgments
Hereby, the authors would like to express gratitude to the Vice Chancellor for research of Iran University of Medical Sciences.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
Conception and design: AT, MB, ShA. Analysis and interpretation of data: AT, MB. Data collection: MB, MAP, BMD. Authors participate in drafting the article or revising: AT, MB, AM, AA, ShA, HD. I attest to the fact that all authors listed on the title page have read and approved the manuscript, attest to 230 the validity and legitimacy of the data and its interpretation, and agree to its submission to BMC Surgery āfor an evaluation and reviewing for maybe publishing.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was approved by the ethics committee of Iran University of medical science.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisherās Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Talebi, A., Mohammadnejad, A., Akbari, A. et al. Survival analysis in gastric cancer: a multi-center study among Iranian patients. BMC Surg 20, 152 (2020). https://doi.org/10.1186/s12893-020-00816-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12893-020-00816-6