
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
Authors' Affiliations: 1 UCD School of Biomolecular and Biomedical Science and 2 UCD School of Medicine and Medical Science, UCD Conway Institute, University College Dublin; 3 Department of Pathology and Laboratory Medicine, St. Vincent's University Hospital, Dublin, Ireland; 4 Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Department of Biostatistics, Harvard School of Public Health, Boston, Massachusetts; 5 Center for Molecular Pathology, Department of Laboratory Medicine, Lund University, Malmö University Hospital, Malmö, Sweden; 6 Department of Epidemiology, University of North Carolina, Chapel Hill, North Carolina; and 7 Tissue Array Research Program, Laboratory of Pathology, Center for Cancer Research, National Cancer Institute, NIH, Bethesda, Maryland
Requests for reprints: William M. Gallagher, UCD School of Biomolecular and Biomedical Science, UCD Conway Institute, University College Dublin, Belfield, Dublin 4, Ireland. Phone: 353-1-7166743; Fax: 353-1-2837211; E-mail: william.gallagher{at}ucd.ie.
| Abstract |
|---|
|
|
|---|
Experimental Design: A breast cancer tissue microarray containing 102 tumors was stained with an anti-survivin antibody. Whole-slide scanning was used to capture high-resolution images. These images were analyzed using automated algorithms to quantify the staining.
Results: Increased nuclear, but not cytoplasmic, survivin was associated with a reduced overall survival (OS; P = 0.038) and disease-specific survival (P = 0.0015). A high cytoplasmic-to-nuclear ratio (CNR) of survivin was associated with improved OS (P = 0.005) and disease-specific survival (P = 0.05). Multivariate analysis revealed that the survivin CNR was an independent predictor of OS (hazard ratio, 0.09; 95% confidence interval, 0.01-0.76; P = 0.027). A survivin CNR of >5 correlated positively with estrogen receptor (P = 0.019) and progesterone receptor (P = 0.033) levels, whereas it was negatively associated with Ki-67 expression (P = 0.04), p53 status (P = 0.005), and c-myc amplification (P = 0.016).
Conclusion: Different prognostic information is supplied by nuclear and cytoplasmic survivin in breast cancer. Nuclear survivin is a poor prognostic marker in breast cancer. Moreover, CNR of survivin, as determined by image analysis, is an independent prognostic factor.
Studies in breast cancer have produced conflicting data about the prognostic influence of survivin. In breast cancer, its prognostic value has been reported to be nonexistent (4–6), associated with improved outcome (7), or associated with adverse outcome (8–10). Despite these conflicting findings, survivin is one of the 16 cancer-related genes represented in the Oncotype DX assay (11). Additionally, following reanalysis of gene expression microarray data published by van 't Veer et al. (12), we identified that survivin gene expression was specifically up-regulated in breast tumors with poor prognosis. Survivin ranked 18th of the top 100 genes associated with poor prognosis in the supervised reanalysis (13).
Immunohistochemical-based studies have failed to reach a consensus about how survivin staining should be interpreted (i.e., does examination of cytoplasmic fraction, nuclear fraction, or both provide more useful information?). In a review of the literature, Li et al. (14) identified 19 publications that measured nuclear survivin in different cancer types and showed that 9 studies concluded that nuclear survivin was associated with an unfavorable prognosis, whereas 5 showed an association with a favorable prognosis. Manual interpretation of survivin is challenging because of the presence of signal within both the cytoplasm and nucleus. This dilemma can be addressed by image analysis, whereby one can accurately identify the nuclear and cytoplasmic fractions of the stain and investigate any link between these fractions and survival.
In this study, we used high-throughput image analysis to fulfill two objectives: (a) the development of an automated quantitative scoring model for survivin and (b) the identification of new prognostic subgroups that were not evident following manual analysis.
| Materials and Methods |
|---|
|
|
|---|
The median age at diagnosis was 60 years (range, 30-87 years). Patients did not receive neoadjuvant treatment and were treated with modified radical mastectomy or wide local excision and axillary lymph node dissection. All patients who had breast-conserving surgery received adjuvant radiotherapy to the remaining breast tissue. Forty-eight percent of patients were lymph node positive. Node-negative patients did not receive adjuvant therapy; however, node-positive premenopausal patients received adjuvant chemotherapy and node-positive postmenopausal patients received adjuvant tamoxifen, concurrent with local recommendations at the time. Eighty-eight percent of tumors were classified as ductal, 6% as lobular, 3% as medullary, and 3% as mucinous. The median tumor size was 2.2 cm (0.8-10 cm).
Following surgical resection, a piece of fresh tumor was immediately snap frozen in liquid nitrogen and stored at –80°C until protein extraction took place. The hormone receptor content was measured by ELISA, with 72% of tumors being estrogen receptor (ER) positive and 62% being progesterone receptor (PR) positive. Diagnostic specimens were all formalin fixed and paraffin embedded in the Department of Pathology in Umeå University Hospital. All tissue blocks were stored in this department before construction of the TMA.
Cell lines. All cell lines (MDA-MB-231, MCF-7, T47D, SKBR3, Hs578T, and BT474) were obtained from the European Collection of Cell Cultures and maintained as previously described (13, 17).
Western blot analysis. Western blot analysis was done as previously described (17). Lysates were separated by SDS-PAGE, and survivin expression was analyzed using a mouse monoclonal anti-human antibody D8 (Santa Cruz Biotechnology, Inc.) at a dilution of 1:50. An anti-β-actin antibody (Abcam) at a dilution of 1:5,000 was used as a loading control.
TMAs and immunohistochemistry. Breast cancer TMAs were constructed as previously described (16) and the blocks were stored in the Department of Pathology in Malmö University Hospital. TMA sections were cut immediately before staining. Sections (4 µm) were dried, deparaffinized, and rehydrated through descending concentrations of ethanol. Heat-mediated antigen retrieval was done using microwave treatment for 2 x 5 min in a citrate buffer before being processed either in the Ventana Benchmark system (Ventana Medical Systems, Inc.) using a prediluted antibody to Her2 (Pathway CB-USA, 760-2694) or in the Dako Techmate 500 system (Dako) for Ki-67 (1:200, M7240; Dako) and survivin (1:50, D8; Santa Cruz Biotechnology).
Manual evaluation of immunohistochemical staining. Survivin expression levels in tumor specimens were evaluated by two pathologists and scored for both nuclear and cytoplasmic staining. All discordant cases were reevaluated and a consensus was reached. Nuclear survivin was scored using the following categories for percent positive cells: 0, 1 to 5 (1+), 6 to 15 (2+), and 16 to 100 (3+). Nuclear and cytoplasmic survivin intensity was divided into four categories: negative (0), mild (1+), moderate (2+), and high staining intensity (3+; Fig. 1 ). The interpretation of Her2, Ki-67, and p53 immunohistochemistry, as well as c-myc amplification, has been previously described (18, 19).
|
Statistical analysis. Spearman's Rho correlation was used to estimate the relationship between duplicate cores from individual tumors, automated and manual analysis, as well as nuclear and cytoplasmic staining intensity. Differences in distribution of clinical data and tumor characteristics between samples with a high and low cytoplasmic-to-nuclear ratio (CNR; described below) were evaluated using the
2 test. Kaplan-Meier analysis and the log-rank test were used to illustrate differences between overall survival (OS), metastases-free survival, and breast cancer–specific survival (BCSS) according to survivin expression. Cox regression proportional hazards models were used to estimate the relationship to OS and BCSS of survivin, Her2, lymph node status, tumor grade, Ki-67, p53, PR, and ER in the patient cohort. Multivariate models included any variable that displayed a significant association with outcome following univariate analysis. All calculations were done using Statistical Package for the Social Sciences version 11.0 (SPSS, Inc.). A P value of <0.05 was considered statistically significant. Random forest clustering (RFC) was performed using R software.8 Hierarchical clustering of random forest data was done in MatLab version 7.4 (The MathWorks, Inc.).
| Results |
|---|
|
|
|---|
Quantitative determination of survivin expression was then ascertained using an image analysis approach, particularly via the use of a commercial positive pixel count algorithm (Aperio Technologies). A pseudocolor "markup" image was generated as an algorithm result, thus allowing the user to confirm that the algorithm was accurately identifying nuclear and cytoplasmic pixels (Fig. 2A ). A full description of the algorithm is available in the Supplementary Data.
|
Associations between nuclear survivin as determined by automated image analysis and survival. We were initially unable to establish any link between survivin nuclear intensity, cytoplasmic intensity or total intensity, and survival, despite analyzing all functions as continuous variables. We, therefore, proceeded to use RFC in an attempt to identify new prognostic subgroups. RFC is an unsupervised strategy that has been used to profile tumors based on TMA data (20, 21). It is attractive for TMA data as it handles highly skewed nonparametric data well (Supplementary Fig. S2).
RFC was done on survivin nuclear intensity and cytoplasmic intensity data, which revealed four distinct clusters (Fig. 2B). The clusters contained the following number of tumors: cluster 1 had 14, cluster 2 contained 32, cluster 3 contained 20, and cluster 4 contained 30. We then determined the mean survivin nuclear intensity and cytoplasmic intensity in all four clusters. As displayed in Fig. 2C, there was a more distinct separation of survivin nuclear intensity in the clusters than survivin cytoplasmic intensity. We performed Kaplan-Meier analysis on the four clusters, which revealed a trend toward worse survival in cluster 4 compared with clusters 1, 2, and 3. When clusters 1 to 3 (n = 66) were combined, they displayed an improved BCSS compared with cluster 4 (n = 30; P = 0.035; Fig. 2D). Given that the significant difference between the clusters was based around nuclear survivin, we proceeded to examine the relationship between nuclear intensity and survival.
One of the issues to consider when examining staining intensity alone is that even a small fraction of staining artifact can significantly alter intensity values. Therefore, to reduce noise secondary to staining artifact, we examined the relationship between nuclear pixel ratio and survival. By using a cutoff of 25th percentile (8%), we also showed a relationship between a high nuclear pixel ratio and decreased BCSS (P = 0.0051; Supplementary Fig. S1E). A survivin nuclear autoscore (SNAS) was devised based on tumors in cluster 4 (n = 30) and staining in >8% of pixels. Twenty-four percent of tumors (n = 23) met these criteria. A positive SNAS (a marker of nuclear expression) was associated with a decreased BCSS (P = 0.0015) and OS (P = 0.0381; Fig. 2E and F). Given the highly significant univariate analysis of SNAS, we proceeded to conduct a Cox multivariate regression analysis of BCSS (Table 1 ) in relation to SNAS. Multivariate analysis revealed that the SNAS was a significant predictor of BCSS in this cohort [hazard ratio (HR), 4.70; 95% confidence interval (95% CI), 1.45-15.26; P = 0.010], along with lymph node status (HR, 3.90; 95% CI, 0.17-2.37; P = 0.027), tumor size (HR, 1.04; 95% CI, 1.01-1.07; P = 0.019), and Her2 status (HR, 3.94; 95% CI, 1.09-13.01; P = 0.025). The relationship between SNAS and other clinicopathologic variables was also examined, with the only significant relationship being between positive SNAS and high grade (P = 0.034; data not shown).
|
|
|
|
| Discussion |
|---|
|
|
|---|
The manual interpretation of immunohistochemical staining remains a time-consuming, subjective process, to which only limited statistical confidence can be assigned due to inherent interobserver and intraobserver variability and the semiquantitative nature of the data (23). The introduction of digital imaging devices and computer-assisted image analysis has provided a major advance toward quantitative description of biological systems. It is anticipated that quantitative scoring systems will, in time, evolve toward integrated environments where the interface between tissue staining, image acquisition, data analysis, and visualization of results will be totally transparent to the end user (24).
This study is the first to use automated quantitative algorithms to analyze survivin immunohistochemical data. Given the fact that survivin expression is associated with the inhibition of apoptosis and the promotion of angiogenesis and proliferation, there is a rationale that increased survivin levels would be associated with an aggressive phenotype (1). Indeed, this has borne out in multiple studies as indicated previously. In breast cancer, however, there have been several conflicting reports about the prognostic value of survivin. Two relatively large studies of 275 patients (9) and 420 patients (8) using cell-free extracts showed that quantitative measurement of survivin expression was a poor prognostic marker associated with high-grade and hormone receptor negativity (8, 9). In general, our data agree with these studies.
By applying image analysis, we showed that increased expression of nuclear survivin, as opposed to cytoplasmic survivin, was associated with a decreased OS and BCSS. Additionally, we applied two novel scoring models (i.e., the SNAS and CNR). We found a significant relationship between grade and SNAS (P = 0.034) and survivin CNR (P < 0.001). Moreover, we showed a significant relationship between low survivin CNR (which is a marker of increased nuclear staining) and ER (P = 0.019) and PR (P = 0.013) negativity. This agrees with the findings of Span et al. (25) who showed that high levels of survivin protein were associated with a poor response to endocrine therapy.
To exclude any confounding variables, our new scoring models were compared with other well-established breast cancer prognostic factors (Table 1). A Cox multivariate regression analysis revealed that SNAS was a significant predictor of BCSS (HR, 4.70; 95% CI, 1.45-15.26; P = 0.010), along with lymph node status (HR, 3.90; 95% CI, 0.17-2.37; P = 0.027), tumor size (HR, 1.04; 95% CI, 1.01-1.07; P = 0.019), and Her2 status (HR, 3.94; 95% CI, 1.09-13.01; P = 0.025).
Likewise, a second Cox multivariate analysis revealed that CNR was a significant predictor of OS (HR, 0.09; 95% CI, 0.01-0.76; P = 0.027), along with tumor size (HR, 1.05; 95% CI, 1.02-1.08; P = 0.002) and lymph node status (HR, 2.71; 95% CI, 1.18-6.24; P = 0.019). These findings would suggest that increased levels of nuclear survivin as determined by image analysis is a predictor of outcome even when controlling for other well-recognized prognostic markers.
As mentioned previously, the prognostic relevance of survivin as measured by immunohistochemistry in breast cancer is a controversial issue and several smaller qualitative studies have produced conflicting results. It is possible that the quantitative measurement of survivin (either by ELISA or image analysis) is necessary for its use as a biomarker in breast cancer. Interestingly, our study and those of Span et al. (25) and Ryan et al. (8) all used quantitative methods to evaluate survivin expression and found similar results. Our study has an added benefit in that it used formalin-fixed, paraffin-embedded materials as opposed to frozen tissue specimens. The use of a TMA as a biomarker discovery and validation platform raises issues of tumor heterogeneity, and definitive utility of a biomarker requires replication on whole patient samples. The fact that we saw an excellent correlation between nuclear intensity (Spearman's Rho = 0.852; P < 0.001) and cytoplasmic intensity (Spearman's Rho = 0.85; P < 0.001) for duplicate cores from individual tumors would suggest that tumor heterogeneity is not a significant contributor to survivin expression. However, this question could only be answered by doing a study of SNAS and survivin CNR on full sections.
There is much debate in the literature about the nuclear-cytoplasmic transport of survivin and its implications for tumorigenesis (26, 27). This study is the first to examine the relationship between the ratio of cytoplasmic to nuclear survivin and outcome. Our results are consistent with the hypothesis that the nuclear and cytoplasmic fractions of survivin have different biological roles (28).
It is well recognized that nuclear-cytoplasmic shuttling of survivin is controlled by an evolutionary conserved Crm1-dependent nuclear export signal (NES; refs. 26, 29, 30). Knauer et al. (26) used small interfering RNA–mediated Crm1 depletion to show that an active NES is necessary to promote the antiapoptotic functions of survivin and inhibition of the NES makes cells more susceptible to chemotherapy- or radiotherapy-induced apoptosis. Colnaghi et al. (30) showed that point mutations in survivin NES abrogated the antiapoptotic effect of survivin. However, they also showed that the mitotic effect of survivin is maintained in NES mutants (30). It could, therefore, be argued that increased levels of nuclear survivin could lead to a proliferative aggressive phenotype.
The role of survivin as a mitotic regulator is well described. Expression of survivin peaks in the G2-M phase of the cell cycle (31) and it is a key member of the chromosomal passenger complex during mitosis, where it regulates levels of the spindle checkpoint protein BubR1 (32–34). One could, therefore, argue that increased levels of nuclear survivin (possibly secondary to a deficient NES) are associated with dysregulation of the mitotic apparatus due to increased levels of BubR1 and a sustained checkpoint response. Our data support this hypothesis, as we saw a strong correlation between low survivin CNR and c-myc amplification, p53 expression, and Ki-67, all of which are associated with a proliferative phenotype (35–37). Furthermore, by applying two novel scoring models, we showed that increased nuclear survivin was associated with poor prognosis. Based on these data, it is not surprising that immunohistochemical-based studies of survivin have produced conflicting results. This may relate to varying specificity of the antibodies used or possibly due to inherent interobserver and intraobserver variability in the manual interpretation of survivin immunohistochemistry. As mentioned previously, it may also be necessary to quantitatively determine survivin levels to fully unearth its prognostic potential in breast cancer. It is in this field that automated image analysis of survivin immunohistochemistry may play an important role and allow for further stratification of patients.
In conclusion, we have presented an alternative method of quantitatively determining the expression of survivin via immunohistochemistry. In more detail, we applied a combination of automated image analysis of survivin immunohistochemistry and RFC to identify new prognostic subgroups. Our data support the hypothesis that the different subcellular pools of survivin have distinct functions and increased levels of nuclear survivin are associated with a proliferative phenotype. Such an approach may be helpful in further dissecting the debate surrounding the role of survivin as a prognostic marker in breast cancer and also as a putative therapeutic target.
| Footnotes |
|---|
Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).
Presented as poster (abstract 49) at the Society of Academic and Research Surgery, 2006, Edinburgh, United Kingdom.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 7/17/07; revised 12/10/07; accepted 1/22/08.
| References |
|---|
|
|
|---|
Ex3, survivin-2B, galectin-3, bag-1, bax-
and MRP-1 mRNAs in breast cancer. Cancer Lett 2003;201:225–36.[CrossRef][Medline]This article has been cited by other articles:
![]() |
M. Krajewska, L. H. Smith, J. Rong, X. Huang, M. L. Hyer, N. Zeps, B. Iacopetta, S. P. Linke, A. H. Olson, J. C. Reed, et al. Image Analysis Algorithms for Immunohistochemical Assessment of Cell Death Events and Fibrosis in Tissue Sections J. Histochem. Cytochem., July 1, 2009; 57(7): 649 - 663. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |