- Open Access
Quantitative comparison and reproducibility of pathologist scoring and digital image analysis of estrogen receptor β2 immunohistochemistry in prostate cancer
Diagnostic Pathology volume 11, Article number: 63 (2016)
Digital image analysis offers advantages over traditional pathologist visual scoring of immunohistochemistry, although few studies examining the correlation and reproducibility of these methods have been performed in prostate cancer. We evaluated the correlation between digital image analysis (continuous variable data) and pathologist visual scoring (quasi-continuous variable data), reproducibility of each method, and association of digital image analysis methods with outcomes using prostate cancer tissue microarrays (TMAs) stained for estrogen receptor-β2 (ERβ2).
Prostate cancer TMAs were digitized and evaluated by pathologist visual scoring versus digital image analysis for ERβ2 staining within tumor epithelium. Two independent analysis runs were performed to evaluate reproducibility. Image analysis data were evaluated for associations with recurrence-free survival and disease specific survival following radical prostatectomy.
We observed weak/moderate Spearman correlation between digital image analysis and pathologist visual scores of tumor nuclei (Analysis Run A: 0.42, Analysis Run B: 0.41), and moderate/strong correlation between digital image analysis and pathologist visual scores of tumor cytoplasm (Analysis Run A: 0.70, Analysis Run B: 0.69). For the reproducibility analysis, there was high Spearman correlation between pathologist visual scores generated for individual TMA spots across Analysis Runs A and B (Nuclei: 0.84, Cytoplasm: 0.83), and very high correlation between digital image analysis for individual TMA spots across Analysis Runs A and B (Nuclei: 0.99, Cytoplasm: 0.99). Further, ERβ2 staining was significantly associated with increased risk of prostate cancer-specific mortality (PCSM) when quantified by cytoplasmic digital image analysis (HR 2.16, 95 % CI 1.02–4.57, p = 0.045), nuclear image analysis (HR 2.67, 95 % CI 1.20–5.96, p = 0.016), and total malignant epithelial area analysis (HR 5.10, 95 % CI 1.70–15.34, p = 0.004). After adjusting for clinicopathologic factors, only total malignant epithelial area ERβ2 staining was significantly associated with PCSM (HR 4.08, 95 % CI 1.37–12.15, p = 0.012).
Digital methods of immunohistochemical quantification are more reproducible than pathologist visual scoring in prostate cancer, suggesting that digital methods are preferable and especially warranted for studies involving large sample sizes.
Significant advances in digital imaging have enabled automated technologies to reproduce and often outperform pathologist visual scoring of immunohistochemistry (IHC) assays. Visual scoring has been the traditional gold standard method for quantifying IHC staining, but problems with this method include the limited range of resulting data [1, 2], human error , less than optimal reproducibility , and resulting ordinal or quasi-continuous variable data rather than true continuous variable data. Digital image analysis overcomes many of these limitations. For example, digital methods allow algorithm parameters to be locked yielding more reproducible data especially when staining is weak and most linearly related to antigen concentration [2, 5, 6], and output continuous variable data. Previous studies reveal that IHC cut-points of biomarkers with prognostic relevance may be identified using continuous variable digital imaging data that were either undetected  or not as strongly associated [2, 8–10] using visual scoring data. Furthermore, digital methods offer a feasible way to scale experiments to high-throughput sample sizes (e.g., experiments using tissue microarrays) which can be otherwise time-limiting for pathologists to complete .
Numerous studies have demonstrated a high degree of correlation between digital image analysis and pathologist visual scoring. The majority of this research has been performed in breast cancer tissue on human epidermal growth factor receptor, estrogen receptor, and progesterone receptor [8, 12–22]. Similar strong correlations between software algorithms and pathologist visual scoring have been described in other tissue types including esophageal cancer , colorectal cancer , ovarian cancer , and prostate cancer (PCa) .
Pathologist visual scoring data often use a simple ordinal variable scale (e.g., negative “0”, weak “1 + ”, medium “2 + ”, and strong “3 + ” positive staining). More complex pathologist visual scoring systems have been developed to provide quasi-continuous variable data, such as multiplying an ordinal variable of intensity by an estimate of tissue area comprising that intensity level [26, 27]. Although studies examining the correlation and reproducibility of pathologist visual scoring and digital image analysis have been performed in breast cancer, to date there has been little research validating such tools in PCa. Few prognostic biomarkers are available for routine clinical use in PCa and the use of digital methods for evaluating IHC assays in large PCa studies represents a valuable technique for evaluating protein biomarkers of tumor aggressiveness . Estrogen receptor β2 (ERβ2) plays a metastasis-promoting role in PCa and has been demonstrated to have prognostic value for tumor progression [29, 30]. Here, we have evaluated the correlation between digital image analysis and pathologist visual scoring (using a semi-quantitative scoring technique ), the reproducibility of these two methods, as well as the association of digital image analysis with disease-specific survival using a large set of PCa tissue microarray (TMA) slides stained for ERβ2.
Clinical cohort and TMA construction
PCa patients (ages 35–74 at diagnosis) of European or African ancestry residing in King County, WA were identified from the Seattle-Puget Sound Surveillance, Epidemiology, and End Results (SEER) cancer registry for population-based studies of PCa risk factors after approval from the Fred Hutchinson Cancer Research Center Institutional Review Board [30–32]. All men were diagnosed with histologically confirmed PCa during either 1993–1996 or 2002–2005. Demographic information and medical history for each patient were obtained by structured in-person interviews, and clinicopathologic data were obtained from the cancer registry. Of the 831 interviewed patients who underwent radical prostatectomy, 566 (68 %) consented to release of tumor tissue including for TMA construction and had formalin-fixed paraffin-embedded blocks available for inclusion in this study, which maintains active Fred Hutchinson Cancer Research Center Institutional Review Board approval (IRF #4714) at the time of submitting this paper for publication. Vital status and underlying cause of death was available for these patients through the biannual linkages with the SEER registry and review of death certificates to confirm cancer-specific vs. other cause mortality. PCa recurrence was determined by follow-up surveys sent to patients in 2004–2005 and in 2010–2011, review of medical records, and physician follow-up as needed. Biochemical recurrence was defined as postoperative prostate-specific antigen (PSA) of ≥2.0 ng/mL. Metastatic progression was confirmed by bone scan, magnetic resonance imaging, computerized tomography scan, or biopsy.
Formalin-fixed, paraffin-embedded blocks of tumor tissue obtained at the time of radical prostatectomy were used to make hematoxylin and eosin (H&E) stained slides, which were reviewed by an experienced genitourinary pathologist. Duplicate tumor tissue cores of 1.0 mm diameter were taken from a single tumor focus (≥75 % tumor tissue) of the donor blocks and arrayed into a new recipient paraffin block with a manual tissue arrayer (MTA-1; Beecher Instruments, Sun Prairie, WI).
Immunohistochemistry and pathologist visual scoring
Unstained 5 μm-thick TMA sections were deparaffinized and rehydrated using standard methods. IHC was performed on TMA sections using a custom polyclonal antibody specific for the C-terminus of ERβ2 (482-MKMETLLPEATMEQ-495) as previously reported . ERβ2 stained slides were scanned for pathologist visual scoring using an automated Tissuefax microscope (TissueGnostics, Tarzana, CA) and reviewed via an online web gallery. ERβ2 immunostaining within malignant cells was scored for each TMA spot by a pathologist (X.Z.) blinded to clinical parameters. Cytoplasm and nuclei were evaluated separately. As described previously, immunostaining was assessed using a score calculated by multiplying staining intensity (0 for no staining, 1 for light/weak staining, and 2 for strong/intense staining) by the corresponding percentage of cells staining positive at each intensity (totaling to 100 %) . Tissue spots that were missing, damaged, contained staining artifacts, or had uncertain histology were excluded from the analysis. Raw data for pathologist visual scores are included in Additional file 1.
Slide digitization, annotation, and immunohistochemical quantification
For digital image analysis, TMA whole slide images were obtained at 40x magnification (0.0625 μm2/pixel) with a ScanScope CS (Aperio ePathology, Leica Biosystems Imaging, Vista, CA) and Genie Histology Pattern Recognition software (Aperio) was trained to classify tissues into Image Classes (tumor, stroma, and glass) as previously described . ERβ2 staining in total malignant epithelial areas was quantified using the Color Deconvolution algorithm (Aperio) as the product of staining intensity (average optical density [OD] units) multiplied by the percentage of tumor epithelium with positive staining (denoted as AvgOD*%Pos). Cytoplasmic staining of ERβ2 within tumor epithelium was quantified using the Cytoplasmic algorithm (Aperio) as the product of staining intensity multiplied by the percentage of tumor epithelium with positive cytoplasmic staining (denoted as AvgCytoOD*%PosCyto). Similarly, nuclear staining was summarized as the average staining intensity within nuclei of tumor epithelium multiplied by the percentage of positive nuclei in tumor epithelium (denoted as AvgNuclearOD*%PosNuclei). These metrics have been previously described [11, 28, 33]. The amount of staining present is linearly related to OD .
A blinded reproducibility study (Analysis Run B) was performed by the same pathologist (X.Z.) who rescored the TMAs using the above protocol (time period of 24 months between initial and repeat scoring). Similarly, the scientist who originally quantified the TMAs re-annotated the TMA spots, retrained the Genie Histology Pattern Recognition software, and reanalyzed the TMAs using the Cytoplasmic and Nuclear algorithms (time period of 10 months between initial and repeat scoring). Raw data for the digital image analyses are included in Additional file 1.
ERβ2 IHC staining was evaluated by digital image analysis and pathologist visual scoring for comparison of quantification methods. The average score across duplicate spots was calculated for each case within each of the two Analysis Runs A and B. The association between pathologist scores and digital measures (AvgCytoOD*%PosCyto or AvgNuclearOD*%PosNuclei) was determined within and across Analysis Runs (A and B) using Spearman’s correlation coefficients, and point estimates and 95 % confidence intervals are presented. P-values represent a test of whether the correlation coefficients are statistically significantly different than 0 (no correlation). Associations of ERβ2 (quantified by image analysis confined to tumor cytoplasm, tumor nuclei, or total malignant epithelial areas) and PCa outcomes (recurrence-free survival [RFS] and prostate cancer-specific mortality [PCSM]) were evaluated using Kaplan-Meier analysis and the log-rank test. Image analysis methods were evaluated using Cox regression models adjusted for age at diagnosis (continuous), Gleason score (≤6, 7[3 + 4], 7[4 + 3], and ≥8), pathologic stage (local: pT2, N0/NX, M0; regional: pT3/pT4 or N1-3, M0), and diagnostic preoperative PSA level. Hazard ratios (HRs) and 95 % confidence intervals (CIs) were reported. A two-tailed p-value of <0.05 was considered statistically significant.
ERβ2 was evaluated by IHC on the PCa patient cohort TMAs. In PCa tissue, ERβ2 displayed variable nuclear staining and variable finely granular cytoplasmic staining, both in malignant epithelial cells and in fibromuscular stromal cells (Fig. 1). In normal prostate tissue, ERβ2 displayed cytoplasmic staining in basal and luminal epithelial cells, agreeing with previously reported specificity and localization .
Correlation of digital image analysis and pathologist visual scoring
Figure 1 demonstrates the workflow for digital image analysis (annotation, automated tissue classification, and image analysis). The digital image analysis process required approximately 1 min of a technician’s time for analysis of each TMA spot (under pathologist supervision) while the visual scoring process required approximately 1 min of a pathologist’s time for analysis of each TMA spot.
Patients were represented in duplicate on the TMAs. Therefore we first compared the correlation of patient replicates to each other (within an Analysis Run). As shown in Table 1, the correlations between replicates were similar for both Analysis Runs A and B, and were higher for digital IHC measures (0.84 and 0.84 in tumor cytoplasmic areas in Analysis Runs A and B, respectively; and 0.85 and 0.84 in tumor nuclear areas in Analysis Runs A and B, respectively) compared to pathologist scores (0.72 and 0.71 in tumor cytoplasmic areas in Analysis Runs A and B, respectively; and 0.64 and 0.62 in tumor nuclear areas in Analysis Runs A and B, respectively). Both methods had relatively high correlation, indicating that replicates within a patient were similar. Since IHC evaluation is often performed for linking to outcome data on a per-patient level, these high correlations provided a rationale for averaging patient data together when comparing quantification methods. When quantification methods were directly compared (patient replicates averaged together), there was a weak/moderate correlation between digital IHC measures and pathology scores of tumor nuclei (Analysis Run A: 0.42 (0.34–0.49), p < 0.0001; and Analysis Run B: 0.41 (0.34–0.48), p < 0.0001; Fig. 2), and a moderate/strong correlation between digital IHC measures and pathology scores of tumor cytoplasm (Analysis Run A: 0.70 (0.65–0.74), p < 0.0001; and Analysis Run B: 0.69 (0.64–0.74), p < 0.0001; Fig. 3).
Reproducibility of quantification methods
To assess the reproducibility of these methods, we performed a second independent analysis of the ERβ2 stained TMAs (Analysis Run B). Comparing data between Analysis Runs A and B, there was a high correlation between pathologist visual scores generated for individual TMA spots (0.84 for tumor cytoplasmic areas and 0.83 for tumor nuclear areas for Analysis Runs A and B, respectively), and very high correlation between digital IHC measures generated for individual TMA spots (0.99 for tumor cytoplasmic areas and 0.99 for tumor nuclear areas for Analysis Runs A and B, respectively) as shown in Table 2.
Select characteristics from the PCa TMA patient cohort are described in Table 3. A total of 508 subjects met eligibility criteria and had suitable tissue available for analysis on the TMA. The mean age at the time of radical prostatectomy was 59.0 years and the median pre-operative diagnostic PSA was 5.9 ng/mL (IQR: 4.6, 9.0). Out of 508 patients, 111 men (21.9 %) experienced PCa recurrence, including 14 men (2.76 %) who died of PCa.
Kaplan-Meier analysis demonstrated that ERβ2 quantified by total malignant epithelial area image analysis was borderline associated with time to recurrence in univariate analysis (p = 0.057; Table 4 and Fig. 4). ERβ2 quantified separately by cytoplasmic image analysis and nuclear image analysis were not significantly associated with time to recurrence in univariate or multivariate analysis (adjusted for clinicopathologic features including age at diagnosis, Gleason score, pathologic stage, and diagnostic PSA level; Table 4 and Fig. 4).
In the univariate analysis of survival, ERβ2 staining was significantly associated with increased risk of PCSM when quantified by cytoplasmic image analysis (HR for each tertile increase = 2.16, 95 % CI 1.02–4.57, p = 0.045), nuclear image analysis (HR 2.67, 95 % CI 1.20–5.96, p = 0.016), and total malignant epithelial area analysis (HR 5.10, 95 % CI 1.70–15.34, p = 0.004). After adjusting for age at diagnosis, Gleason score, pathologic stage, and diagnostic PSA level, ERβ2 staining was significantly associated with increased risk of PCSM when quantified by total malignant epithelial area image analysis (HR 4.08, 95 % CI 1.37–12.15, p = 0.012). We did observe a borderline significant association with increased risk of PCSM when ERβ2 staining was quantified by nuclear image analysis (HR 2.32, 95 % CI 0.99–5.41, p = 0.052).
The long-term aims of our studies are to examine the association of PCa biomarkers with outcomes at the patient level. Traditionally, we and others in our field have used pathologist visual scoring, which has been shown to have good to excellent intra- and inter-pathologist reproducibility [4, 6, 35]. However, approximations of positive-staining area has only poor to good pathologist reproducibility , highlighting the need for improved methods. Although we have demonstrated that the digital image analysis process takes similar quantities of time to complete compared to visual scoring for the current analysis, digital image analysis can be completed by a technician supervised by a pathologist, whereas visual scoring requires significantly more pathologist time. In this study, pathologist time for supervision of the digital image analysis required approximately 10-fold less pathologist time versus visual scoring of the TMAs, similar to our previously published findings . Pathologist availability and efficiency limits visual scoring approaches, especially for large studies [36, 37]. Although we have not evaluated inter-observer pathologist reproducibility in this study, this is another important source of error and thus inter-pathologist and inter-scientist reproducibility studies warrant further investigation.
An initial finding in this study was that the correlation between patient replicate TMA spots was higher for digital IHC measures (in both tumor cytoplasmic and nuclear areas) compared to pathologist scores. These data held up across two independent Analysis Runs and raise an interesting discussion point for tumor heterogeneity and its relation to outcome. High correlation between patient replicate TMA spots in our study conforms to prior research indicating that a relatively small number of cores adequately represent the tumor, although this is highly dependent on the antigen being evaluated . In PCa specifically, a 12 biomarker signature has demonstrated high correlation between patient replicate TMA spots of varying tumor grades indicating that expression of the signature in either high or low Gleason grade similarly informed outcome [39, 40]. Other tumor types may have a higher degree of molecular heterogeneity, and this may contribute to outcome. For example, melanoma tumor cells distant from vasculature demonstrate altered expression of numerous hypoxia-related genes, and potentially react to these localized environments in ways that may be critical for disease aggressiveness . Similar experiments in breast and rectal cancer demonstrate distinct expression patterns at the leading/invasive edges of tumor compared with trailing/center portions, a phenomenon referred to as a “prairie fire” antigen distribution [42, 43]. These data demonstrate that genetic and histopathologic spatial heterogeneity may be reflected in the biologic behavior of cells within distinct tumor areas. Additional studies directly exploring the relationship between heterogeneity of tumor biomarker expression and outcome are needed.
One previous study in PCa identifies a high correlation between digital analysis and ordinal pathologist scores of ERG, SLC45A3, and TMPRSS2 IHC . However, continuous data allows the use of statistical methods more suitable to identifying IHC cut-points of biomarkers with prognostic relevance [2, 7–9]. For this reason, pathology studies have developed semi-quantitative scoring methods [26, 27]. However, multiplying ordinal by continuous data does not produce a true continuous variable but rather a quasi-continuous variable. Problems with quasi-continuous scoring systems are exemplified by Rimm et al. who showed bimodal distribution of pathologist visual scoring due to over-calling of very weak staining as “negative” rather than recognizing the weaker staining which often displays the most variability when quantitatively evaluated . Ideally, much research would move in the direction of truly quantitative methods where staining intensity (perhaps detected using fluorescence methods that have wider dynamic ranges than IHC assays) is reflected on a standard curve of controls with biochemically known target antigen quantity .
Schade et al. previously demonstrated that pathologist visual scoring of separate nuclear (intense only) and cytoplasmic (intense only) ERβ2 immunohistochemical staining was associated with a higher risk of PCSM in the same PCa cohort evaluated in the current study . Here, we extended this work by assessing the association between ERβ2 quantified by multiple image analysis methods with PCa outcomes, and identified a significant association with PCSM when ERβ2 was quantified by total malignant epithelial area, and identified a borderline significant association with PCSM when ERβ2 was quantified by nuclear-only staining after adjusting for multiple clinicopathologic factors similar to our previous report . While we observed a significant association with PCSM when ERβ2 was quantified by cytoplasmic image analysis in univariate analysis, this result did not remain significant after adjusting for clinicopathologic factors.
Our current findings build upon the work by Schade et al. showing that ERβ2 is associated with adverse outcomes. However, it is unclear why we found that data obtained by pathologist visual scoring, versus data obtained from the same slides using digital image analysis, yielded slightly different associations with patient outcome metrics (RFS and PCSM). It is possible that low level staining that is present and quantifiable by digital methods, may be interpreted as “negative” by a pathologist relying on visual interpretation of staining intensity , resulting in misclassification. Relatively little is published in this area, with some groups suggesting that digital data result in higher associations than visual scoring with outcome metrics [10, 45] and one instance describing a lower association with digital data . Further work, out of scope for the present study, is required to identify factors that may underlie differences of visual versus digital image analysis data and their correlation with outcome metrics.
Our study, to our knowledge, is the first to assess both reproducibility of pathology visual data and reproducibility of digital methods in the same pathology data set. We demonstrated that digital methods are extremely reproducible across two Analysis Runs which involved re-annotation of tissues, retraining of a pattern recognition algorithm to identify tumor epithelium, and reevaluation and compilation of data. We conclude that computer-aided methods may produce improved datasets and lead to higher quality and more reproducible research, especially in studies involving large sample sizes.
CI, confidence interval; ERβ2, estrogen receptor-β2; H&E, hematoxylin and eosin; HR, hazard ratio; IHC, immunohistochemistry; OD, optical density; PCa, prostate cancer; PCSM, prostate cancer-specific mortality; PSA, prostate-specific antigen; RFS, recurrence-free survival; SEER, surveillance, epidemiology, and end results; TMA, tissue microarray
Rimm DL. What brown cannot do for you. Nat Biotechnol. 2006;24(8):914–6. doi:10.1038/nbt0806-914.
Rimm DL, Giltnane JM, Moeder C, Harigopal M, Chung GG, Camp RL, et al. Bimodal population or pathologist artifact? J Clin Oncol Off J Am Soc Clin Oncol. 2007;25(17):2487–8. doi:10.1200/JCO.2006.07.7537.
Weaver DL, Krag DN, Manna EA, Ashikaga T, Harlow SP, Bauer KD. Comparison of pathologist-detected and automated computer-assisted image analysis detected sentinel lymph node micrometastases in breast cancer. Mod Pathol. 2003;16(11):1159–63. doi:10.1097/01.MP.0000092952.21794.AD.
Jaraj SJ, Camparo P, Boyle H, Germain F, Nilsson B, Petersson F, et al. Intra- and interobserver reproducibility of interpretation of immunohistochemical stains of prostate cancer. Virchows Arch. 2009;455(4):375–81. doi:10.1007/s00428-009-0833-8.
Bloom K, Harrington D. Enhanced accuracy and reliability of HER-2/neu immunohistochemical scoring using digital microscopy. Am J Clin Pathol. 2004;121(5):620–30. doi:10.1309/Y73U-8X72-B68T-MGH5.
Gavrielides MA, Gallas BD, Lenz P, Badano A, Hewitt SM. Observer variability in the interpretation of HER2/neu immunohistochemical expression with unaided and computer-aided digital microscopy. Arch Pathol Lab Med. 2011;135(2):233–42. doi:10.1043/1543-2165-135.2.233.
Harigopal M, Barlow WE, Tedeschi G, Porter PL, Yeh IT, Haskell C, et al. Multiplexed assessment of the Southwest Oncology Group-directed Intergroup Breast Cancer Trial S9313 by AQUA shows that both high and low levels of HER2 are associated with poor outcome. Am J Pathol. 2010;176(4):1639–47. doi:10.2353/ajpath.2010.090711.
Atkinson R, Mollerup J, Laenkholm AV, Verardo M, Hawes D, Commins D, et al. Effects of the change in cutoff values for human epidermal growth factor receptor 2 status by immunohistochemistry and fluorescence in situ hybridization: a study comparing conventional brightfield microscopy, image analysis-assisted microscopy, and interobserver variation. Arch Pathol Lab Med. 2011;135(8):1010–6. doi:10.5858/2010-0462-OAR.
Camp RL, Dolled-Filhart M, King BL, Rimm DL. Quantitative analysis of breast cancer tissue microarrays shows that both high and normal levels of HER2 expression are associated with poor outcome. Cancer Res. 2003;63(7):1445–8.
Feuchtinger A, Stiehler T, Jutting U, Marjanovic G, Luber B, Langer R, et al. Image analysis of immunohistochemistry is superior to visual scoring as shown for patient outcome of esophageal adenocarcinoma. Histochem Cell Biol. 2014. doi:10.1007/s00418-014-1258-2.
Rizzardi AE, Johnson AT, Vogel RI, Pambuccian SE, Henriksen J, Skubitz AP, et al. Quantitative comparison of immunohistochemical staining measured by digital image analysis versus pathologist visual scoring. Diagn Pathol. 2012;7:42. doi:10.1186/1746-1596-7-42.
Bolton KL, Garcia-Closas M, Pfeiffer RM, Duggan MA, Howat WJ, Hewitt SM, et al. Assessment of automated image analysis of breast cancer tissue microarrays for epidemiologic studies. Cancer Epidemiol Biomark Prev. 2010;19(4):992–9. doi:10.1158/1055-9965.EPI-09-1023.
Brugmann A, Eld M, Lelkaitis G, Nielsen S, Grunkin M, Hansen JD, et al. Digital image analysis of membrane connectivity is a robust measure of HER2 immunostains. Breast Cancer Res Treat. 2012;132(1):41–9. doi:10.1007/s10549-011-1514-2.
Chosidow O, Bagot M, Vernant JP, Roujeau JC, Cordonnier C, Kuentz M, et al. Sclerodermatous chronic graft-versus-host disease. Analysis of seven cases. J Am Acad Dermatol. 1992;26(1):49–55.
Faratian D, Kay C, Robson T, Campbell FM, Grant M, Rea D, et al. Automated image analysis for high-throughput quantitative detection of ER and PR expression levels in large-scale clinical studies: the TEAM Trial Experience. Histopathology. 2009;55(5):587–93. doi:10.1111/j.1365-2559.2009.03419.x.
Gokhale S, Rosen D, Sneige N, Diaz LK, Resetkova E, Sahin A, et al. Assessment of two automated imaging systems in evaluating estrogen receptor status in breast carcinoma. Appl Immunohistochem Mol Morphol. 2007;15(4):451–5. doi:10.1097/PAI.0b013e31802ee998.
Krecsak L, Micsik T, Kiszler G, Krenacs T, Szabo D, Jonas V, et al. Technical note on the validation of a semi-automated image analysis software application for estrogen and progesterone receptor detection in breast cancer. Diagn Pathol. 2011;6:6. doi:10.1186/1746-1596-6-6.
Laurinaviciene A, Dasevicius D, Ostapenko V, Jarmalaite S, Lazutka J, Laurinavicius A. Membrane connectivity estimated by digital image analysis of HER2 immunohistochemistry is concordant with visual scoring and fluorescence in situ hybridization results: algorithm evaluation on breast cancer tissue microarrays. Diagn Pathol. 2011;6:87. doi:10.1186/1746-1596-6-87.
Masmoudi H, Hewitt SM, Petrick N, Myers KJ, Gavrielides MA. Automated quantitative assessment of HER-2/neu immunohistochemical expression in breast cancer. IEEE Trans Med Imaging. 2009;28(6):916–25. doi:10.1109/TMI.2009.2012901.
Skaland I, Ovestad I, Janssen EA, Klos J, Kjellevold KH, Helliesen T, et al. Comparing subjective and digital image analysis HER2/neu expression scores with conventional and modified FISH scores in breast cancer. J Clin Pathol. 2008;61(1):68–71. doi:10.1136/jcp.2007.046763.
Turashvili G, Leung S, Turbin D, Montgomery K, Gilks B, West R, et al. Inter-observer reproducibility of HER2 immunohistochemical assessment and concordance with fluorescent in situ hybridization (FISH): pathologist assessment compared to quantitative image analysis. BMC Cancer. 2009;9:165. doi:10.1186/1471-2407-9-165.
Turbin DA, Leung S, Cheang MC, Kennecke HA, Montgomery KD, McKinney S, et al. Automated quantitative analysis of estrogen receptor expression in breast carcinoma does not differ from expert pathologist scoring: a tissue microarray study of 3,484 cases. Breast Cancer Res Treat. 2008;110(3):417–26. doi:10.1007/s10549-007-9736-z.
Alexander BM, Wang XZ, Niemierko A, Weaver DT, Mak RH, Roof KS, et al. DNA repair biomarkers predict response to neoadjuvant chemoradiotherapy in esophageal cancer. Int J Radiat Oncol Biol Phys. 2012;83(1):164–71. doi:10.1016/j.ijrobp.2011.05.033.
Messersmith W, Oppenheimer D, Peralba J, Sebastiani V, Amador M, Jimeno A, et al. Assessment of Epidermal Growth Factor Receptor (EGFR) signaling in paired colorectal cancer and normal colon tissue samples using computer-aided immunohistochemical analysis. Cancer Biol Ther. 2005;4(12):1381–6.
Braun M, Kirsten R, Rupp NJ, Moch H, Fend F, Wernert N, et al. Quantification of protein expression in cells and cellular subcompartments on immunohistochemical sections using a computer supported image analysis system. Histol Histopathol. 2013;28(5):605–10.
Allred DC, Clark GM, Elledge R, Fuqua SA, Brown RW, Chamness GC, et al. Association of p53 protein expression with tumor cell proliferation rate and clinical outcome in node-negative breast cancer. J Natl Cancer Inst. 1993;85(3):200–6.
FitzGerald LM, Zhang X, Kolb S, Kwon EM, Liew YC, Hurtado-Coll A, et al. Investigation of the relationship between prostate cancer and MSMB and NCOA4 genetic variants and protein expression. Hum Mutat. 2013;34(1):149–56. doi:10.1002/humu.22176.
Rizzardi AE, Rosener NK, Koopmeiners JS, Isaksson Vogel R, Metzger GJ, Forster CL, et al. Evaluation of protein biomarkers of prostate cancer aggressiveness. BMC Cancer. 2014;14:244. doi:10.1186/1471-2407-14-244.
Leung YK, Lam HM, Wu S, Song D, Levin L, Cheng L, et al. Estrogen receptor beta2 and beta5 are associated with poor prognosis in prostate cancer, and promote cancer cell migration and invasion. Endocr Relat Cancer. 2010;17(3):675–89. doi:10.1677/ERC-09-0294.
Schade GR, Holt SK, Zhang X, Song D, Wright JL, Zhao S, et al. Prostate cancer expression profiles of cytoplasmic ERbeta1 and nuclear ERbeta2 are associated with poor outcomes following radical prostatectomy. J Urol. 2016. doi:10.1016/j.juro.2015.12.101.
Agalliu I, Salinas CA, Hansten PD, Ostrander EA, Stanford JL. Statin use and risk of prostate cancer: results from a population-based epidemiologic study. Am J Epidemiol. 2008;168(3):250–60. doi:10.1093/aje/kwn141.
Stanford JL, Wicklund KG, McKnight B, Daling JR, Brawer MK. Vasectomy and risk of prostate cancer. Cancer Epidemiol Biomark Prev. 1999;8(10):881–6.
Rizzardi AE, Vogel RI, Koopmeiners JS, Forster CL, Marston LO, Rosener NK, et al. Elevated hyaluronan and hyaluronan-mediated motility receptor are associated with biochemical failure in patients with intermediate-grade prostate tumors. Cancer. 2014;120(12):1800–9. doi:10.1002/cncr.28646.
Krajewska M, Smith LH, Rong J, Huang X, Hyer ML, Zeps N, et al. Image analysis algorithms for immunohistochemical assessment of cell death events and fibrosis in tissue sections. J Histochem Cytochem. 2009;57(7):649–63. doi:10.1369/jhc.2009.952812.
Borlot VF, Biasoli I, Schaffel R, Azambuja D, Milito C, Luiz RR, et al. Evaluation of intra- and interobserver agreement and its clinical significance for scoring bcl-2 immunohistochemical expression in diffuse large B-cell lymphoma. Pathol Int. 2008;58(9):596–600. doi:10.1111/j.1440-1827.2008.02276.x.
Ong CW, Kim LG, Kong HH, Low LY, Wang TT, Supriya S, et al. Computer-assisted pathological immunohistochemistry scoring is more time-effective than conventional scoring, but provides no analytical advantage. Histopathology. 2010;56(4):523–9. doi:10.1111/j.1365-2559.2010.03496.x.
Tzeng CH, Chen PM, Chuang MW, Liu JH, Hsieh RK, Liu CJ, et al. Treatment of severe aplastic anemia: comparison of bone marrow transplantation to immunotherapy. Zhonghua Yi Xue Za Zhi (Taipei). 1989;43(1):21–8.
Moeder CB, Giltnane JM, Harigopal M, Molinaro A, Robinson A, Gelmon K, et al. Quantitative justification of the change from 10 to 30% for human epidermal growth factor receptor 2 scoring in the American Society of Clinical Oncology/College of American Pathologists guidelines: tumor heterogeneity in breast cancer and its implications for tissue microarray based assessment of outcome. J Clin Oncol. 2007;25(34):5418–25. doi:10.1200/JCO.2007.12.8033.
Shipitsin M, Small C, Choudhury S, Giladi E, Friedlander S, Nardone J, et al. Identification of proteomic biomarkers predicting prostate cancer aggressiveness and lethality despite biopsy-sampling error. Br J Cancer. 2014;111(6):1201–12. doi:10.1038/bjc.2014.396.
Shipitsin M, Small C, Giladi E, Siddiqui S, Choudhury S, Hussain S, et al. Automated quantitative multiplex immunofluorescence in situ imaging identifies phospho-S6 and phospho-PRAS40 as predictive protein biomarkers for prostate cancer lethality. Proteome Sci. 2014;12:40. doi:10.1186/1477-5956-12-40.
Widmer DS, Hoek KS, Cheng PF, Eichhoff OM, Biedermann T, Raaijmakers MI, et al. Hypoxia contributes to melanoma heterogeneity by triggering HIF1alpha-dependent phenotype switching. J Invest Dermatol. 2013;133(10):2436–43. doi:10.1038/jid.2013.115.
Avoranta ST, Korkeila EA, Syrjanen KJ, Pyrhonen SO, Sundstrom JT. Lack of CD44 variant 6 expression in rectal cancer invasive front associates with early recurrence. World J Gastroenterol. 2012;18(33):4549–56. doi:10.3748/wjg.v18.i33.4549.
Chung GG, Zerkowski MP, Ghosh S, Camp RL, Rimm DL. Quantitative analysis of estrogen receptor heterogeneity in breast cancer. Lab Invest. 2007;87(7):662–9. doi:10.1038/labinvest.3700543.
McCabe A, Dolled-Filhart M, Camp RL, Rimm DL. Automated quantitative analysis (AQUA) of in situ protein expression, antibody concentration, and prognosis. J Natl Cancer Inst. 2005;97(24):1808–15. doi:10.1093/jnci/dji427.
Desmeules P, Hovington H, Nguile-Makao M, Leger C, Caron A, Lacombe L, et al. Comparison of digital image analysis and visual scoring of KI-67 in prostate cancer prognosis after prostatectomy. Diagn Pathol. 2015;10:67. doi:10.1186/s13000-015-0294-0.
Parker AS, Lohse CM, Leibovich BC, Cheville JC, Sheinin YM, Kwon ED. Comparison of digital image analysis versus visual assessment to assess survivin expression as an independent predictor of survival for patients with clear cell renal cell carcinoma. Hum Pathol. 2008;39(8):1176–84. doi:10.1016/j.humpath.2007.12.009.
This study utilized services from NWBioSpecimen which was supported by National Cancer Institute grant P30 CA015704 (F. Appelbaum), Life Sciences Discovery Fund (LSDF) grant Washington Phenotyped Biospecimen Resource (J. Slattery), LSDF grant Consortium Biospecimen Program (P. Porter), Fred Hutchinson Cancer Research Center, the University of Washington School of Medicine, and the University of Washington Department of Pathology. The work was also supported by grants R01-CA056678, R01-CA092579, R03-CA137799, and P50-CA097186 (J. Stanford) from the National Cancer Institute, with additional support from the Fred Hutchinson Cancer Research Center, the Prostate Cancer Foundation, and the Institute for Prostate Cancer Research.
Availability of data and materials
Raw data for pathologist visual scoring and digital image analysis are available in Additional file 1.
AER, SCS, and JLS designed and coordinated the study. AER and JK performed image and data analysis. XZ performed pathologist visual scoring. RIV and MSG provided statistical analysis. SK and JCH participated in study coordination and database administration. YL and SH coordinated and performed immunohistochemical staining. All authors contributed to and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All patients participating in this study consented to release and use of their biospecimens for research as authorized by active Fred Hutchinson Cancer Research Center Institutional Review Board approval (IRF #4714).
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Cite this article
Rizzardi, A.E., Zhang, X., Vogel, R.I. et al. Quantitative comparison and reproducibility of pathologist scoring and digital image analysis of estrogen receptor β2 immunohistochemistry in prostate cancer. Diagn Pathol 11, 63 (2016). https://doi.org/10.1186/s13000-016-0511-5
- Prostate cancer
- Digital pathology
- Estrogen receptor β2