## Modelling of survival of patients with colon adenocarcinoma based on multivariable analysis of the state of cancer cell nuclear apparatusGrabovoy A.N. ^{1}National Cancer Institute, Kyiv, 03022, Ukraine^{2}ESC “Institute of Biology” Taras Shevchenko National University, Kyiv 03022, Ukraine
Submitted: October 12, 2015.
*Correspondence: E-mail: angrabovoy@gmail.com Abbreviations used: АС — adenocarcinoma; AUC — area under ROC curve; D — average DNA content rank in nucleus up to 1.2 c.u.; D+ — average DNA content rank in nucleus up to 1.2–2.5 c.u.; DM — medium optical density of nucleus; G — tumor grade; LR — likelihood ratio test; NDNA — DNA content in nucleus; NOR — nuclear organizer region; nNOR — number of NORs in nucleus; Т+ — average DNA content rank in nucleus higher than 2.5 c.u.; vNOR — total volume of NOR in nucleus.
Convenient systems for evaluation based on differentiation grade of colon adenocarcinoma (АС) sometimes are not optimal due to the difficulties of objective assessment and prognostic validity required for a choice of correct treatment of AC patients [1, 2]. Therefore, the search for new criteria and the development of prognostic means is of importance for clinicians [3–6]. In the process of neoplastic transformation nuclear apparatus undergoes significant changes caused by genome damage and its altered replication during cell division and realization of genetic information. In particular, structural organization and size of nuclei, state of nucleolar organizer regions, DNA quantity (polyploidy and aneuploidy), immunohistochemical markers of proliferation and apoptosis regulation, and a great variety of other features, have been studied. On large cohorts of AC patients there has been shown a relation of some patterns of the state of tumor cell nuclei with the properties of neoplasms and prognosis. However, when these schemes have been applied for prognosis in an individual patient, they were found to be of low significance with large error range. Therefore, discrete patterns could be considered as such which use for prognosis is of limited accuracy. Some patterns of nuclear state are directly or indirectly interrelated, in particular, mitotic activity, DNA quantity (synthesis, ploidy, aneuploidy), the state of nuclear organizer regions (NOR), apoptosis regulation. Increased and uncontrolled mitotic activity is among the main features of tumor growth. This could be accompanied with malfunction of mitosis mechanisms (upon chromosomal instability) leading to the development of polyploidy or aneuploidy. Along with this, DNA synthesis preceding division, at conditions of high proliferation activity results also in a certain increase of total DNA quantity in tumor cells. From other side, lesion of mechanisms of apoptosis and survival of the cells, especially the cells with increased content of DNA in nuclei, also promotes an elevation of medium DNA content in tumor cell nuclei. The aim of the work: the creation of a mathematical model of survival in patients with colon AC based on multivariable analysis of the state of cancer cell nuclear apparatus. ## Patients and MethodsThe study was performed on 141 samples of biopsy materials or materials obtained during surgical treatment of the patients with colon AC or benign colon neoplasms: conventionally normal colon tissue (C; n = 11); polyps and adenomas (В; n = 16); AC of different grades: G1 (n = 28); G2 (n = 73); G3 (n = 13). The material was fixed in a buffered 10% formaline (pH 7.4, 4 °С, 24 h) and embedded in paraplast with the use of Histos-5 histoprocessor (Milestone, Italy). From paraffin blocks, 5 µm sections were cut using Microm НМ325 (Thermo Scientific, Germany). The sections were stained with hematoxylin and eosin for common tumor assessment, or with Einarson’s gallocyanin chrome alum stain (pH 1.62, 37 °С, 24 h) for determination of nucleic acid content [7]. In each case, a part of sections was treated with RNAase (Macherey-Nagel GmbH & Co. KG, Germany) for RNA destruction [8]. NOR were detected using silver nitrate impregnation. Immunohistochemical reactions were performed with the use of mouse monoclonal antibodies against antigen Ki‑67 (MIB-1, Dako, Denmark), Bcl-2 oncoprotein (Dako, Denmark), and р53 protein (Dako, Denmark) according to the instructions of the manufacturer and visualized using detection system EnVision The preparations were examined and photographed using microscope Nikon Eclipse 80i supplied with camera DS-5SMc/L2 at standardized conditions, images (× 400, 1280 × 960 pixels RGB) were processed using analysis system ImageJ 1.46. In images of the preparations stained with gallocyanin chrome alums, in each 30 tumor cells there have been determined cross-sectional area of cell nucleus, medium (DM) and integral optical density of nucleus, and the quantity of nuclear DNA (NDNA). For calculation of the latest index, DNA content in lymphocyte nuclei was accepted as a unit (c.u.) [9]. The tumors were distributed into three ranks by average DNA content: D — average DNA content up to 1.2 c.u., D+ — average DNA content from 1.2 to 2.5 c.u., T+ — average DNA content higher than 2.5 c.u. [10]. In silver-impregnated sections, in each 50 tumor cells there was determined cross-sectional area of nucleus, quantity of NOR in nucleus (nNOR), and diameter of each NOR with following calculation of their total volume (vNOR). In each tumor, the percent of cells expressing Ki-67, Bcl-2, and р53 has been determined. Statistical analysis was performed with the use of IBM SPSS Statistics 22. For assessment of normality, one sample Kolmogorov — Smirnov test was used. Correlation analysis was performed using Spearman rank correlation method. The differences between the groups were evaluated by Mann — Whitney U test. Survival analysis was performed by Kaplan — Meier estimator using log-rank test for evaluation of the differences between the groups. For creation of survival models, Cox proportional hazards model (Cox regression model) was used. Standard form of regression equation is:
where: x — the studied characteristics of tumor cell nuclear apparatus: X_{i}_{1} — 0 if AC differentiation grade is G2, 1 — if AC differentiation grade is G3; X_{2} — ratio of nNOR/Ki-67 values; X_{3} — product of nNOR × NDNA values; X_{4} — derivative of nNOR × NDNA/Ki-67; X_{5} — DNA content rank, may be of three values (D, D+ or T+); X_{6} — 0 corresponds to rank D+, 1 — to rank T+.To build the models of prognosis on the basis of Cox regression, different number of nuclear apparatus features were used. Model I includes just a single characteristic X Model II includes X Model III includes all 4 characteristics of models I and II. Model IV includes 5 characteristics: X Model V also includes 5 characteristics: X The prognostic validity of the developed models was compared with that of basic Cox model that doesn’t include any characteristics of tumor process with the use of likelihood ratio test (LR). As an index of prognosis model quality (survival of patients) there has been used ROC-analysis where area under the curve (AUC) reflects model quality as follows: 0.5–0.6 — bad, 0.6–0.7 — fair, 0.7–0.8 — good, 0.8–0.9 — very good, 0.9–1.0 — perfect. The differences were considered significant if p < 0.05. ## Results and DiscussionThe performed study has shown that average DNA content in cell nuclei of tumors of different grade varied from 2.06 (G3) to 2.38 (G1) (Table 1). While the differences between groups C and B were significant (р < 0.05), the differences between various G were insignificant (р > 0.05). In the studied sample among AC of different grades (В, G1, G2, G3) there were the tumors with different DNA content ranks (D, D+, T+). If one takes into account the variability of this index, its use as a criterion of tumor grading and independent prognostic factor seems to be improper [9, 11]. Also, the decrease of DNA content ranks in nuclei of tumor cells along with G increase in combination with restriction of cellular heterogeneity via reduction on the number of hyperploid cells [9], supports the realization of clonal evolution during AC progression [8, 12]. Table 1. Quantitative values of features of nuclear apparatus of tumor cells of colon AC and benign colon neoplasms Table 1. Quantitative values of features of nuclear apparatus of tumor cells of colon AC and benign colon neoplasms
Note: NDNA — content of DNA in nucleus (cu); D, D+, T+ – DNA content ranks; nNОR — number of nucleolar organizers; vNОR — total volume of nucleolar organizers; Ki-67, Bcl-2, p53 — expression level of corresponding marker in tumor; p — significance: *p < 0.05, **p < 0.01 (comparing with previous group).Number of NОR in tumor cell nuclei tended to decrease along with increase of grade while their total volume (vNОR) has been increasing; however, the differences between these indexes in approximal tumor groups were insignificant. Therein, the direct dependence between nNОR/vNОR was not found. An absence of proportional increase of nNОR/vNОR in tumor cell nuclei along with elevation of DNA content is evidence of DNA dysfunction or inactivity [9]. This phenomenon could be considered as a manifestation of decreased viability of these cells leading to ageing, apoptosis, and elimination of such morphofunctional types from a tumor [13, 14]. All studied tumors have shown Ki-67 expression index higher than zero, which ranged from 2.17% to practically 100% with a median of 37%. In general, there was recorded a sharp rise of tumor cell Ki-67 expression index in group B The performed correlation analysis has shown that there are a number of statistically significant weak relations between the parameters. G correlated with nNОR (r = −0.340; p < 0.01) and their average volume (r = −0.348; p < 0.01), but had no correlation with vNОR. Apart from this, there was recorded a relation between increase of Ki-67 labeling index and tumor grade (r = 0.193; p < 0.05), and also with nNОR (r = −0.202; p < 0.05). Positive correlation was determined also for p53 and G (r = 0.227; p < 0.05). Survival analysis dependent on discrete patterns has shown that it significantly differed only in the groups distributed only by grade (G) (p < 0.01). No statistically significant differences between Bcl-2(+/−) and p53(+/−) groups has been revealed (Table 2). However, in the tumors with different DNA content ranks there has been observed a tendency for survival change. It was the highest in rank D and the lowest in rank T+. This fact could be explained by the different mechanisms of development of these neoplasms [15]. Tumors from rank D may be refer to ones which develop via epigenetic mechanisms and microsatellite, while ranks D+ and T+ — chromosome instability. The data obtained in correlation analysis and survival analysis in the groups allowed to build a survival model using Cox regression (proportional Table 2. Comparing of patients survival with colon AC in different groups using log-rank test (null hypothesis — there is no difference between groups)
Note: G — tumor grade; D, D+, T+ — DNA content ranks; Bcl-2 — Bcl-2 expression in tumor; p53 — expression in tumor.Table 3. Prognostic significance of discrete features of nuclear apparatus of tumor cells (Cox model)
Note: β — coefficient; SE — standard error of coefficient; p — significance; nNOR — number of nucleolar organizers; vNOR — total volume of nucleolar organizers; Ki-67, Bcl-2, p53 — expression level of corresponding marker in tumor; NDNA — quantity of DNA cu in nucleus; DM — medium optical density of nucleus; G — tumor grade; D, D+, T+ — DNA content ranks.As far as tumor grading (G) is the only one significant prognostic criterion, it has been used in the following analysis as a base for model building. Apart from this, to elevate model accuracy groups B and G1 were excluded because in these groups there has been recorded 100% 5-year survival rate, what preconditions their uninformativeness. For assessment of prognostic validity of a combination of the factors, we have performed a number of procedures of direct and inverse inclusion of the variables in regression model. There have been used either discrete factors or their combinations and their interrelation has been analyzed (Table 4). Inclusion of two tumor grading values — G2 and G3 into regression equation produced a significant model (Model І), but with low AUC (0.339) (Figure). However, positive value β = 1.87 indicates that the higher G increases the likelihood of lethal outcome. The search for significant combination of factors allowed to include such elements as NDNA, nNOR and Ki-67 which formed a base for Model ІІ. Its analysis has shown that simultaneous increase of NDNA and nNOR increases the risk of death. However, increase of nNOR in parallel with the decrease of Ki-67 expression index elevates the risk. The quality of this model is practically equal from the first one what indicates their similarly low prognostic validity (see Table 4). Joining of Models I and II into a new Model ІІІ has been found to be rational because it resulted in a decrease of LR and an increase of AUC up to 0.605 (see Figure), what however couldn’t be considered statistically significant at p = 0.18. Better result was obtained after inclusion of feature of DNA content ranks (D, D+, T+) into regression equation. This procedure led to the model of better quality (Model IV): AUC = 0.663 (see Table 4). Further, the best model was obtained using D+ and T+ ranks (Model V), AUC = 0.722 — fair quality. Thus, Model V based on 5 features: X Table 4. Regression parameters and features included to the model upon regression analysis of survival (Cox model) of AC patients
Note: n — number of cases included into analysis; LR — likelihood ratio (comparing with null model); AUC — area under ROC curve; β — coefficient; SE — standard error of coefficient; p () — significance for LR, AUC and β; X_{1} — tumor grade (G2 or G3); X_{2 }= nNOR/Ki-67; X_{3 }= nNOR × NDNA; X_{4 }= nNOR × NDNA/Ki-67; X_{5 }— DNA content ranks (D, D+ or T+); X_{6 }— DNA content ranks (D+ or T+).Figure. ROС-analysis of different survival models in colon AC. I–V — ROС curves for different models
So, the performed studies have shown that discrete patterns of AC cell nuclei state (DNA quantity, number and volume of NOR, expression rates of Ki-67, Bcl-2 and p53) are prognostically invalid in a case of their separate use and could not be considered as direct prognostic criteria. However, the derivatives of a number of patterns could serve as significant criteria for prognosis (in particular, nNOR/Ki-67, nNOR × NDNA, and nNOR × NDNA/Ki-67). Besides, upon their use in survival model mathematically expressed dependence is not linear and is described by more complex functions (in this case, exponential one). Good quality survival model based on Cox proportional Table 5. Estimated survival function for patients with colon AC (Model V) for mean values of features
Note: t — survival time, months; S_{0}(t) — baseline survival at time t; S(t) — survival function at time t.The proposed survival model in the case of AC demonstrates the quality twice superior to that of a model based just on the use of tumor grade (G), which presently is a single common independent histological criterion of prognosis [3, 2]. ## REFERENCES- Compton CC, Fielding PL, Burgart LJ,
. Prognostic factors in colorectal cancer. Arch Pathol Lab Med 2000;*et al***124**: 979–94. - WHO classification of tumors the digestive system. Lyon: IARC, 2010: 131–81.
- Barresi V, Bonetti LR, Ieni A,
. Histological grading in colorectal cancer: new insights and perspectives. Histol Histopathol 2015;*et al***30**: 1059–67. - Bolocan A, Ion D, Ciocan D,
. Prognostic and predictive factors in colorectal cancer. Chirurgia (Bucur) 2012;*et al***107**: 555–63. - Coppedè F, Lopomo A, Spisni R,
. Genetic and epigenetic biomarkers for diagnosis, prognosis and treatment of colorectal cancer. World J Gastroenterol 2014;*et al***20**: 943–56. - Kulendran M, Stebbing J, Marks C,
Predictive and prognostic factors in colorectal cancer: a personalized approach. Cancers (Basel) 2011;*et al.***29**: 1622–38. - Luppa Kh. Basic histochemistry. Moscow: Myr, 1980. 344 p. (in Russian).
- Sottoriva A, Kang H, Ma Z,
. A Big Bang model of human colorectal tumor growth. Nat Genet 2015;*et al***47**: 209–16. - Grabovoy AN, Antoniuk SA, Vorobey EA. Nucleic acids contents in the colon epithelial tumor cells nuclei. Ukr Morfolohichnyi Almanakh 2013;
**11:**73–6 (in Ukrainian). - Avtandylov HH. Diagnostic medical ploidometry. Moscow: Medytsyna, 2006. 192 p. (in Russian).
- Kørner H. Microsatellite instability and DNA ploidy in colorectal cancer. Cancer 2009;
**115**: 271–82. - Ding L, Raphael BJ, Chen F, Wendl MC. Advances for studying clonal evolution in cancer. Cancer Lett 2013;
**340**: 212–9. - Davoli T, de Lange T. The causes and consequences of polyploidy in normal development and cancer. Ann Rev Cell Dev Biol 2011;
**27**: 585–610. - Holland AJ, Cleveland DW. Losing balance: the origin and impact of aneuploidy in cancer. EMBO 2012;
**1****3**: 501–14. - Haigis K. Molecular pathogenesis of colorectal cancer. Spr Science Business Med 2013;
**6**: 320.
## No Comments » Add comments |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

## Leave a comment