A prognostic model based on cell-cycle control predicts outcome of breast cancer patients

Background A prognostic model combining biomarkers of metaphase-anaphase transition of the cell cycle was developed for invasive breast cancer. The prognostic value and clinical applicability of the model was evaluated in comparison with the routine prognosticators of invasive breast carcinoma. Methods The study comprised 1135 breast cancer patients with complete clinical data and up to 22-year follow-up. Regulators of metaphase-anaphase transition were detected immunohistochemically and the biomarkers with the strongest prognostic impacts were combined into a prognostic model. The prognostic value of the model was tested and evaluated in separate patient materials originating from two Finnish breast cancer centers. Results The designed model comprising immunoexpressions of Securin, Separase and Cdk1 identified 8.4-fold increased risk of breast cancer mortality (p < 0.0001). A survival difference exceeding 15 years was observed between the majority (> 75%) of patients resulting with favorable as opposed to unfavorable outcome of the model. Along with nodal status, the model showed independent prognostic impact for all breast carcinomas and for subgroups of luminal, N+ and N- disease. Conclusions The impact of the proposed prognostic model in predicting breast cancer survival was comparable to nodal status. However, the model provided additional information in N- breast carcinoma in identifying patients with aggressive course of disease, potentially in need of adjuvant treatments. Concerning N+, in turn, the model could provide evidence for withholding chemotherapy from patients with favorable outcome.


Background
Cell proliferation, hormone-regulation and HER2 amplification are considered the main biological processes driving breast cancer progression. Although proliferation has been shown a valid prognosticator in all subtypes, particularly triple-negative breast carcinoma (TNBC) has been characterized by high expression of proliferationrelated genes [1,2]. The prognostic value of proliferation is acknowledged in the clinical pathology practice as part of the traditional histological grading as well as in intrinsic classification and in modern personalized signatures retrieved from microarray-based expression-profiling [3][4][5]. However, the impact of deregulated proliferation is still not accurately reflected in the routine clinical parameters and pathological markers applied to treatment decisions of breast cancer patients.
Genetic integrity of the dividing cell is ensured by complex and intricately monitored cellular events at the metaphase-anaphase transition of the cell cycle [6]. Dysfunction of these regulators can lead into missed sister chromatid separation, chromosomal instability and aneuploidy. Premature sister chromatid separation is prevented by the highly controlled sequential activation and inactivation of a cascade of regulatory proteins, particularly Cdc20 (cell division cycle 20), Cohesin, Separase (Extra Spindle Pole Bodies Like protein 1, ESPL1), Securin (Pituitary tumor-transforming gene 1, Pttg1), Pttg1IP (Pituitary tumor-transforming gene 1 interacting protein, also Pituitary tumor-transforming gene 1 binding factor, PBF), Cdk1 (Cyclin-Dependent Kinase protein 1) and CyclinB1 (G2/mitotic-specific cyclin-B1). In more detail, correct segregation of the chromosomes is triggered at the Spindle Assembly Checkpoint (SAC) controlled by Cdc20 activating Anaphase-Promoting Complex / Cyclosome (APC/C) to create the APC/ C Cdc20 complex. Throughout metaphase, contact between the chromatids is maintained by rings of Cohesin. At the initiation of anaphase, Cohesin is removed triggered by Separase and activated by degradation of Securin and/or the Cdk1/CyclinB1 complex, leading to separation of the sister chromatids [7,8].
In our previous research, PTTG1, the gene of human Securin, was detected with the most significant expression difference between human breast cancer and normal breast glandular tissue on basis of a cDNA microarray analysis involving 4000 cancer related genes [9]. In addition to Securin, also several other regulators of metaphase-anaphase transition have been shown with independent prognostic impacts in breast cancer [10][11][12][13][14][15][16][17][18]. In the present study, we introduce on basis of a total of 1135 patients with up to 22-year follow-up, a clinically applicable combination of biomarkers of metaphaseanaphase transition leading to optimal detection of aggressive course of disease and cancer mortality in invasive breast cancer.

Patient materials
The study comprises patients diagnosed and treated with unilateral invasive breast carcinoma in two different institutions ( For all patients, the biomarkers of the metaphaseanaphase transition were immunohistochemically (IHC) detected, and the routine clinico-pathological prognostic features of breast cancer were collected. For both cohorts of breast cancer patients, a prognostic model was assembled based on the optimal combination of biomarkers of the metaphase-anaphase transition. The prognostic value of the models was evaluated in comparison with the clinically applied routine prognosticators of breast cancer. All patients were treated with surgical resection or mastectomy with axillary evacuation, radiation and/or adjuvant treatment with anti-estrogenic or cytostatic drugs depending on the patients' age, hormone receptor and lymph node (N) status according to the international guidelines for breast cancer treatment at the time of diagnosis [19]. No pre-operative adjuvant treatments were administered. Complete clinical data was collected from pathology reports and patient files and registered applying the criteria presented by WHO [20] and St. Gallen International Expert Consensus [21]. Intrinsic subtypes were approximated by immunohistochemistry according to international guidelines [22]. Causes of death were obtained from autopsy reports, death certificates and from the national cancer registry (Statistics Finland, Helsinki, Finland).

Tissue materials
Tissue materials were prepared according to standard histology practice, i.e. fixed in buffered formalin (pH 7.0) and embedded into paraffin blocks. Tissue micro arrays (TMAs) were prepared using the representative tumor area of each patient. The TMAs included two tissue cores with diameters 0.6 mm (cohort I) or 1 mm (cohort II) from each tumor.

Interpretation of IHC
Immunoexpressions for Securin, Separase, Cdc20, Pttg1IP, SA2, Cdk1 and CyclinB1 were observed as combinations of nuclear and cytoplasmic staining and registered as average fractions (%) of positively staining cancer cells [13,14,16,23]. In each case, the number of immunopositive cells was calculated in sets of one hundred cancer cells (minimum 100 and maximum 3 × 100 cancer cells evaluated) and registered as an average fraction (%) of immunopositivity for each patient. Interpretations for IHC of Ki-67, ER, PR, and IHC and ISH for HER2 followed previous literature and generally accepted international guidelines [22,24]. All IHC interpretations were performed by experienced histopathologists (HR, PK).

Statistical analysis
The cutpoints for immunoexpressions of the studied biomarkers were set based on previous literature, histopathological observations and statistical analyses involving the mean, median and univariate prognostic values of each parameter [9,13,14,16,17]. In prognostic analyses, Kaplan-Meier estimates were performed to demonstrate the cumulative percentages of breast cancer specific mortality. Cox's proportional hazard models were used to test associations between disease outcome, biomarker expressions and clinical prognostic features, i.e. tumor size, axillary lymph node status, histological and intrinsic classifications and histological grade. The risk of breast cancer death associated with the studied proteins and the routine prognosticators was quantitated as hazard ratio (HR) with 95% confidence interval (CI). P-values < 0.05 were considered statistically significant. The computations were performed with SAS for Windows, Version 9.3 (I-III) and 9.4 (IV) (SAS Institute, Cary, NC, USA). Kaplan-Meier survival plots were generated using R 2.15.0.

Results
In IHC of the invasive breast carcinomas, Securin was detected as predominantly nuclear but occasionally, showed both nuclear and/or cytoplasmic  immunoreaction. Separase showed two distinct and apparently mutually excluding expression patterns observed in the nucleus or in the cytoplasm of cancer cells. Nuclear immunoexpression was observed for Cdc20, SA2, Cdk1 and CyclinB1 whereas Pttg1IP showed cytoplasmic expression only. Table 2 summarizes the fractions of immunopositive breast carcinomas among all breast cancer subtypes (cohort I) and in subgroups divided according to tumor size, nodal status, histological grade and survival. The prognostic impacts of the studied regulators of metaphase-anaphase transition were first analyzed for all breast cancer subtypes (cohort I, n = 781) with Cox's proportional hazard model. Among all studied biomarkers, statistically significant prognostic value in univariate analyses was observed for Securin (HR 2.1, p < 0.0001, CI 1.6-2.8), nuclear Separase (HR 2.0, p < 0.0004, CI 1.4-3.0), Cdk1 (HR 2.5, p < 0.0001, CI 1.7-3.6) and CyclinB1 (HR 2.0, p = 0.04, CI 1.0-2.0).
In the next phase, these biomarkers were further tested in combinations in order to assemble a prognostic model producing the most significant prognostic impact among all breast cancer subtypes. The optimal model for detecting favorable outcome of disease was determined as the combination of low expressions for Securin (< 10% of cancer cells), Separase (< 1% of cancer cells) and Cdk1 (< 10% of cancer cells). This model was a significant indicator of survival of disease while the opposite detected patients in risk of breast cancer death (HR 8.4, p < 0.0001). The Kaplan-Meier curves demonstrate the survival difference among all patients and in subgroups with N+ and N-disease (Fig. 1). Concluding from the survival analyses, favorable outcome of the model indicated that the majority (> 75%) of patients were alive 18.4 years after primary diagnosis while unfavorable outcome of the model suggested that one quarter (25%) of the patients were already dead of breast cancer after 2.5 years of diagnosis. Among the subgroup of N-patients, no cancer-related deaths were observed among patients exhibiting favorable outcome of the model. Instead, the unfavorable outcome suggested cancer mortality for every fourth patient within 5.3 years from diagnosis. Correspondingly, the majority of N+ patients with favorable and unfavorable outcome of the model were alive after 17.6 and 2.0 years from the primary diagnosis, respectively.
Finally, in multivariable analyses of all breast cancer subtypes (Table 3), the designed model was compared with the established clinical prognosticators of breast cancer, i.e. tumor size, nodal status, intrinsic classification and histological grade for superior prognostic power in predicting the risk of breast cancer mortality. In the whole material, significant prognostic impact was observed for axillary lymph node status along with the designed model. Independent prognostic value was detected for the model also among N+ and N-patients. Tumor size (diameter < 2 cm vs ≥2 cm), intrinsic classification or histological grade (1-2 vs 3) were not detected with independent prognostic value in any of the performed analyses.
As the cohort I contained a relatively old patient material, a more recent patient material (cohort II) was collected to evaluate the prognostic impact of the regulators of metaphase-anaphase transition in luminal (n = 208) and triple-negative breast carcinomas (n = 146). As the result of luminal breast carcinomas, Securin (HR 1.1, p = 0.02, CI 1.0-1.2) and nuclear Separase (HR 5.7, p = 0.002, CI 1.9-17.2)but not Cdk1 remained the most powerful predictors of cancer mortality, along with nodal status (HR 4.9, p = 0.003, CI 1.7-13.7). Also as combined into a model, Securin and Separase showed significant prognostic impact (p = 0.0006). Similar trends for the prognostic model were observed also in separate analyses of N+ (n = 42) and N-(n = 166) patients although the associations were not statistically significant in the small patient groups.
In TNBC (n = 146), no statistically significant prognostic value could be detected for any of the immunohistochemically studied biomarkers or clinicopathological features. The hazard ratio of breast cancer death ns not statistically significant

Discussion
In the present study, prognostic models involving regulators of the metaphase-anaphase transition of the cell cycle were introduced for predicting survival of breast cancer patients. The model was assembled based on data from a total of 1135 breast cancer patients with complete clinical information and up to 22-year followup.
The results show that combining high immunoexpressions for Securin, Separase and Cdk1 comprises a promising prognostic model indicating 8.4-fold increased risk of breast cancer death (p < 0.0001). In luminal breast carcinomas, the combination of Securin and Separase resulted in independent prognostic impact. In all analyses, the prognostic value of the combination of Securin and Separase with or without Cdk1 outperformed the impact of tumor size and histological grade whereas axillary lymph node status remained a strong and independent prognosticator in all analyses. At highest, the introduced prognostic model predicted 19.5-fold increased risk of breast cancer death among N-(p = 0.006) and 6.5-fold increased risk of mortality among N+ breast carcinomas (p = 0.0003). This suggests that the model may provide additional prognostic information to nodal status in treatment decisions of breast cancer patients. Among Ndisease, the model may provide information to identify patients with aggressive course of disease, potentially in need of adjuvant treatments. In N+, the model could provide evidence for withholding chemotherapy from patients with favorable outcome.
TNBC is an aggressive subtype of breast cancer resulting in high mortality. It is also a therapeutically challenging subgroup as the lack of estrogen and progesterone receptors limits the treatment options available. The ongoing intense research has not yet revealed promising novel prognostic biomarkers or treatment targets for TNBC [2]. The present study did not reveal a prognostic impact for the routine clinicopathological parameters or for the studied biomarkers in TNBC. This may be due to the small cohort size as well as the heterogeneous nature of TNBC comprising several molecular subtypes [2].
The value of the introduced model lies in the pivotal role of the studied biomarkers in cancer progression. Loss of control of the cell cycle is a hallmark of malignancy resulting in aneuploidy and genomic instability [25]. Cell cycle checkpoints, including the mitotic checkpoint SAC, are the major cell cycle control mechanism and specifically deregulated in cancer cells. The presently studied regulators of metaphase-anaphase transition are involved in ensuring the fidelity of chromosome segregation. In previous literature, the studied biomarkers have been shown with prognostic impact in breast carcinoma as well as in other malignancies [14,16,17,[26][27][28][29][30][31][32][33]. In addition, numerous strategies have been proposed for the design of cell cycle-selective therapies in cancer, including targeting the metaphaseanaphase transition [8]. In all, the introduced model appears to identify biological drivers that could add to the conventional histopathological evaluation on the proliferative behavior of the tumor.
The reliability of the introduced prognostic model is increased by repeated statistical analysis of the studied biomarkers. The first approach applied for all breast cancer subtypes (cohort I) was based on testing combinations of biomarkers for their independent prognostic impact in comparison to each other, the clinical prognostic features and disease survival. In the second approach applied for analysis luminal and triple-negative breast carcinomas (cohort II), the optimal combination of biomarkers was extracted together with the clinical prognostic markers in sequential multivariate analyses of breast cancer-specific survival. Each approach was independently applied on separate patient materials originating from two Finnish breast cancer centers. As a conclusion, irrespective of the patient material or statistical procedure, the combination of Securin and Separase showed the most significant prognostic impact among all breast carcinomas and among luminal breast cancer.
From the practical point of view, immunohistochemical detection of Securin and Separase comprises a biology driven, cost-effective and reliable prognostic method (Fig. 2). By comparison, Ki-67 -the established proliferation marker of cancer diagnostics -has been criticized for high variability across individual pathologist and institutions, as well as for poor prognostic value [34,35]. In the literature, numerous prognostic models have been introduced for breast cancer but only a few of them have shown impact exceeding that of the routine clinical prognosticators [36][37][38].

Conclusions
In our results from a total of 1135 breast cancer patients with complete clinical information and up to 22-year follow-up, high immunoexpression for the combination of Securin and Separase comprises a powerful prognostic tool to identify patients in risk of breast cancer death. In our scenario, the proposed model may facilitate personalized clinical decision suggesting less aggressive therapy for patients associated with low risk of mortality and indicating benefits from adjuvant therapy for a subgroup of patients with aggressive disease. Despite the accumulating data on multiparametric prognostic models, the clinical judgement remains the key determinant on selecting between the different treatment schemes.
Additional file 1. Details of IHC for detecting Securin, Separase, Cdc20, Pttg1IP, SA2, Cdk1 and CyclinB1 and routine clinical biomarkers of breast cancer ER, PR, Ki-67 and HER2. The datasheet contains the details of the immunohistochemical methods used in this study. Authors' contributions HR participated in the design, histological and immunohistochemical analyses and writing of the manuscript, EL performed the statistical analyses, SK and LK selected and provided the patient specimen from Auria biobank, TK selected and provided the patient specimen from the Biobank of Central Finland, KT participated in analyzing the data and writing the manuscript, PK supervised the study and prepration of the manuscript. All authors read and approved the final manuscript.

Funding
This study has been financially supported by the Cancer Society of Finland, Turku University Hospital and Finska Läkaresällskapet foundation. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials
The data that support the findings of this study are available from Auria Biobank (www.auriabiobank.fi), Turku University Hospital, Turku, and Biobank of Central Finland, Jyväskylä, Finland but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Biobanks stated above.

Ethics approval and consent to participate
The research was performed in accordance with the ethical standards of institutional and national research committees approved by the Regional Ethical Review Boards of Turku University Hospital and Auria Biobank, Turku, Finland, Central Hospital of Central Finland, Jyväskylä, Finland, and Finnish Cancer Registry, Cancer Society of Finland, Helsinki, Finland (permit numbers 6/2002, AB15-9859 and TK-53-716-16). All research procedures were performed in accordance with the 1964 Helsinki declaration and its amendments (www.wma.net/policies-post/wma-declaration-of-helsinki-ethicalprinciples-for-medical-research-involving-human-subjects) and the Finnish Biobank Act (688/2012, Ministry of Social Affairs and Health, Helsinki, Finland). The patients have been included in the study according to the Finnish Biobank Act 688/2012, Ministry of Social Affiars and Health (https://www.finlex.fi/ en/laki/kaannokset/2012/20120688). According to the Act, the national biobank organization, after ethical pre-evaluation, allows collection, storage and processing of patient specimen and information for medical research.