Originally published as JHC exPRESS on June 23, 2008. doi:10.1369/jhc.2008.950345
Volume 56 (10): 873-880, 2008 Copyright ©The Histochemical Society, Inc. Synergistic Tissue Counterstaining and Image Segmentation Techniques for Accurate, Quantitative Immunohistochemistry
Small Animal Imaging Laboratory, McConnell Brain Imaging Centre, Montreal Neurological Institute, McGill University, Montreal, Quebec, Canada Correspondence to: Barry J. Bedell, McGill University, 3801 University Street, BT-209, Montreal, Quebec, Canada H3A 2B4. E-mail: bbedell{at}bic.mni.mcgill.ca
Quantitative analysis of digitized IHC-stained tissue sections is increasingly used in research studies and clinical practice. Accurate quantification of IHC staining, however, is often complicated by conventional tissue counterstains caused by the color convolution of the IHC chromogen and the counterstain. To overcome this issue, we implemented a new counterstain, Acid Blue 129, which provides homogeneous tissue background staining. Furthermore, we combined this counterstaining technique with a simple, robust, fully automated image segmentation algorithm, which takes advantage of the high degree of color separation between the 3-amino-9-ethyl-carbazole (AEC) chromogen and the Acid Blue 129 counterstain. Rigorous validation of the automated technique against manual segmentation data, using Ki-67 IHC sections from rat C6 glioma and β-amyloid IHC sections from transgenic mice with amyloid precursor protein (APP) mutations, has shown the automated method to produce highly accurate results compared with ground truth estimates based on the manually segmented images. The synergistic combination of the novel tissue counterstaining and image segmentation techniques described in this study will allow for accurate, reproducible, and efficient quantitative IHC studies for a wide range of antibodies and tissues. (J Histochem Cytochem 56:873–880, 2008)
Key Words: immunohistochemistry quantitative analysis image segmentation tissue counterstaining
IMAGE ANALYSIS of digitized IHC-stained tissue sections provides a powerful tool for quantification of protein expression. The use of quantitative IHC (qIHC) has increased in recent years for numerous applications, including diagnostic and prognostic determinations in the clinical setting, and correlation with complementary quantitative measures, such as real-time PCR, in research laboratories. Although many investigators perform qIHC analysis on chromogen-stained tissue sections without the use of a counterstain, tissue counterstaining presents a number of practical advantages. Among these advantages are orientation of tissue on the slide, assessment of tissue morphology, specific determination of positively and negatively staining structures and cells, and intensification of chromogen staining for improved visualization. For these reasons, IHC sections are routinely counterstained in most anatomical pathology laboratories. Furthermore, the recent advent of commercially available, ultra-high-resolution, digital slide scanners has heightened interest in whole slide and whole tissue qIHC measures for both clinical and research applications. Accurate three-dimensional reconstruction of qIHC sections also requires counterstained tissue to allow for proper alignment of the two-dimensional sections comprising the three-dimensional volume (Brey et al. 2002
Counterstaining of IHC-stained tissue, however, complicates segmentation of the digital IHC image into chromogen-positive and chromogen-negative pixels caused by the color convolution of the chromogen and counterstain. A number of methods have been developed to allow for semi- or fully automated segmentation of IHC images. These techniques have been thoroughly reviewed and rigorously compared by Brey et al. (2003)
One strategy for improving the accuracy of image segmentation and classification is to increase the contrast between the component-of-interest and the background tissue and use a tailored segmentation algorithm based on the resultant image contrast. Such a strategy has been successfully used in other areas of medical imaging, including analysis of multiple sclerosis (MS) lesions on MRI scans (Bedell et al. 1997 In this study, we describe the development of a new tissue counterstain and associated slide preparation techniques that results in homogeneous background staining of tissue, as well as a straightforward, robust, automated image segmentation and classification algorithm that takes advantage of the high degree of separation between the light blue color of the counterstain and the dark-red color of the AEC chromogen. The accuracy of this technique has been evaluated by comparison of the results of the automated algorithm with manual classification of positively stained nuclei on rat C6 glioma tissue IHC stained for the nuclear cell cycle marker Ki-67, as well as on brain tissue from transgenic mice with amyloid precursor protein (APP) mutations, which has been IHC stained for β-amyloid.
Animal Models and Tissue Preparation Rat C6 glioma cells were purchased from American Type Culture Collection (Rockville, MD) and grown in DMEM supplemented with 10% FBS, 125 U/ml penicillin G, 125 µg/ml streptomycin sulfate, and 2.2 µg/ml amphotericin B (Fungizone). All culture reagents were obtained from Gibco BRL (Invitrogen; Burlington, Ontario, Canada). Cultures were grown in monolayers and maintained at 37C in a humidified atmosphere of 5% CO2. On reaching confluency, spheroids were prepared using the hanging drop method previously described by Del Duca et al. (2004) All animal experiments were conducted in accordance with the guidelines of the Canadian Council on Animal Care and the Montreal Neurological Institute and McGill University Institutional Animal Care and Use Committees. Five male Sprague-Dawley rats (250–300 g; Charles River Canada, St. Constant, Quebec, Canada) were anesthetized with 50 mg/kg ketamine and10 mg/kg xylazine. The right cortical surface in the parietal-occipital region was exposed by craniectomy using a high-powered drill (DREMEL; Racine, WI), and the underlying dura and its vessels were carefully removed under a surgical microscope. A piece of the cortex was removed to expose the underlying white matter, and a single spheroid was placed into the surgical defect. The craniectomy was covered with bone wax (Ethicon; Peterborough, Canada), and the overlying skin was sutured. After recovery from anesthesia, the animals were fed and had access to water ad libitum. After 16–18 days, animals were sacrificed, and brains were removed and immediately immersed in 10% neutral-buffered formalin.
Tissue with β-amyloid deposition was obtained from the brains of four 18-month-old transgenic APPSw/Ind line J20 mice (Mucke et al. 2000
IHC
Tissue Counterstaining and Slide Preparation After slide digitization (described in the following section), the coverslips were removed by soaking slides in distilled water for 5–10 min. The Acid Blue 129 counterstain was readily removed by soaking slides in 2 mM aqueous ammonia for 5–10 min. This technique completely removed the Acid Blue 129 while preserving the AEC chromogen. The slides were subsequently counterstained with hematoxylin for 1 min and mounted in permanent aqueous mounting media (Aquatex; EMD Chemicals, Gibbstown, NJ). These hematoxylin-counterstained sections were digitized for comparison to the Acid Blue 129 counterstaining.
Imaging Regions of interest (ROIs) were selected from each of the slides for comparison of the automated segmentation algorithm with manual segmentation. Specifically, cellular regions and regions adjacent to necrotic areas were identified on each Ki-67 IHC slide, whereas regions containing both parenchymal β-amyloid plaques and vascular β-amyloid deposits were identified on each β-amyloid IHC slide. ROIs were captured at x400 effective magnification and stored in the bitmap (BMP) image file format.
Manual Segmentation
Automated Image Segmentation Algorithm
Analysis
The MetricThe overlap metric was used to assess the level of agreement between the competing technique and the ground truth estimate, where
Typically,
Percentage of True Positives
Percentage of False Positives
The quality of the ground truth estimates was first evaluated to assess the level of agreement between the two manual raters. The TP and FP were calculated by alternating each of the manual rater's labels as the ground truth estimate. These values were examined using a multivariate ANOVA (MANOVA) with the raters' values as the main effect and covariates of IHC stain. Statistical differences were determined with two tailed t-tests.
The automated segmentation technique was assessed against the ground truth estimates provided by the manual raters. The results from the manual raters (i.e., single
Automated Segmentation of Ki-67 and β-amyloid IHC-stained Tissue Sections The automated algorithm was found to detect AEC chromogen-positive pixels with a high level of sensitivity for both the Ki-67 and β-amyloid IHC-stained tissue sections, which had been counterstained with Acid Blue 129. Representative examples of the results generated by the automated algorithm using Acid Blue 129–counterstained, Ki-67 and β-amyloid IHC-stained tissue sections are shown in Figure 2 . Blood was not found to be a confounding factor for the automated segmentation, because the erythrocytes stained dark blue (Figure 2A, arrow), whereas the necrotic regions were pale blue (Figure 2A, arrowhead). As such, the blood had a high blue-to-red ratio and was well separated from the dark-red AEC chromogen.
The Acid Blue 129 counterstain provided a greater contrast-to-noise ratio (CNR) than the hematoxylin counterstain, thereby allowing for the use of higher blue-to-red ratio and green intensity level thresholds in the segmentation algorithm to maximize TP while minimizing FP. Furthermore, the homogeneous nuclear and cytoplasmic staining provided by the Acid Blue 129 eliminated the false negatives resulting from the color convolution of the chromogen and dark nuclear hematoxylin staining in nuclei with condensed chromatin (e.g., vascular endothelial cells). Figure 3 shows a side-by-side comparison of the same IHC ROIs counterstained with Acid Blue 129 and hematoxylin. Note that the use of the Acid Blue 129-counterstain reduced the number of false negatives on both the segmented Ki-67 and β-amyloid IHC images.
Comparison of Manual and Automated Segmentation Data The TP and FP results for the manual segmentations of the data (summarized in Table 1 ) were analyzed to determine the level of agreement between the manual raters. The MANOVA analysis showed no effect of IHC stain (F = 1.0651, df = 1, p<0.3174), but a significant effect of rater (F = 3.0141, df = 1, p<0.0157). The results pooled for each of the stains showed significant inter-rater differences in the β-amyloid stain (t = 3.17, df = 5, p<0.033) but not for the Ki-67 stain (t = –1.17, df = 9, p<0.865). These results indicated high inter-rater variability with respect to the definition of positive β-amyloid staining, but a high level of agreement in the definition of positive Ki-67 staining. FP analysis also showed no effect of IHC (F = 1.0651, df = 1, p<0.3174); however, significant differences between raters were observed (F = 5.865, df = 1, p<0.0353). When these results were pooled across each stain, significant differences between β-amyloid segmentations were also identified (t = 3.17, df = 5, p<0.033), whereas no significant differences for the Ki-67 segmentations were observed (t = –1.09, df = 9, p<0.750). This analysis further confirmed the inter-rater variability in the segmentations of β-amyloid plaques, thereby underscoring the intrinsic limitations of manual segmentation techniques.
The overlap metric was used to evaluate the accuracy of the automated technique compared with the ground truth estimate. The results of the ANOVA analysis for the metric are shown in Table 2
. Although the MANOVA analysis did not show any significant effect of the IHC stain (F = 2.25, df = 1, p<0.1531), it did show a significant effect of the segmentation technique (F = 4.79, df = 2, p<0.025). The Tukey-Kramer HSD post hoc analysis showed that the value of the automated segmentation technique had a high level of agreement compared with Rater 1 (Group A) but a lower level of agreement with Rater 2 (Group B), indicating the variability associated with the manual segmentation technique.
The nature of the differences observed between the manual and automated segmentation techniques was further explored. The results of the ANOVA for the TP tests are shown in Table 3 . The MANOVA analysis did not show any significant effect of the IHC stain (F = 1.0651, df = 1, p<0.3174). However, this analysis indicated a significant effect of the segmentation technique (F = 5.29, df = 3, p<0.0353). The Tukey-Kramer HSD post hoc analysis showed that the TP value of the automated segmentation technique had a higher level of agreement compared with Rater 1 or Rater 2 (Group A) than the inter-rater comparison.
The results of the ANOVA for the FP tests are shown in Table 4 . In this case, significant differences were observed for both the segmentation technique (MANOVA results: F = 4.75, df = 3, p<0.0049) and the IHC stain (MANOVA results: F = 12.11, df = 1, p<0.0009). Further analysis of the individual IHC stains by ANOVA indicated that there was a significant effect of segmentation technique for the β-amyloid stain (F = 8.40, df = 3, p<0.0004) but not for the Ki-67 stain (F = 1.9734, df = 3, p<0.145). The Tukey-Kramer HSD post hoc analysis for the β-amyloid stain showed a high level of variability between the segmentation techniques.
We successfully developed a robust, fully automated method for the quantification of AEC-stained IHC sections. The Acid Blue 129 stain provided homogeneous tissue counterstaining and a high degree of contrast between tissue and chromogen. The use of this new counterstaining system allowed for the implementation of a simple segmentation algorithm, which takes advantage of the high level of color separation between the dark-red AEC chromogen and the uniform, light-blue Acid Blue 129 counterstain. Once the optimal threshold criteria were determined based on trial images, and fixed for the algorithm, this automated technique provided results comparable to the ground truth manual segmentation results.
When the automated segmentation technique was compared with both raters, as assessed by the We determined that the use of the Acid Blue 129 counterstain minimized false positives and false negatives compared with the hematoxylin counterstain. The higher CNR provided by the Acid Blue 129 allowed for a greater detection of true chromogen staining while minimizing false positives in the background tissue. Furthermore, the homogenous staining produced by the Acid Blue 129 eliminated the convolution of dark-blue nuclear staining with chromogen staining, which is often produced with hematoxylin counterstaining and can result in false negatives. The superior performance of the algorithm using Acid Blue 129–counterstained sections compared with hematoxylin-counterstained tissue may be advantageous for IHC studies in which the tissue morphology is dominated by cells with highly condensed nuclear chromatin (e.g., neuroendocrine tumors and lymphomas).
The IHC quantification method used in this study was based on a binary classifier (i.e., chromogen positive or chromogen negative). Significant interest exists, however, in quantification of the absolute or relative level of protein expression (Shi et al. 2005
While Acid Blue 129 is superior to hematoxylin as a counterstain for quantification of IHC sections, it does not provide the same level of morphological detail as hematoxylin. In cases in which both quantification and morphological assessment are desired, Acid Blue 129 can be simply removed by washing with 0.2% aqueous ammonia and the sections subsequently re-counterstained with hematoxylin, as described in the Materials and Methods section. Given that coverslips can be removed from slides prepared with our specially formulated mounting media by soaking in distilled water for
The recent advent of ultra-high-resolution, whole slide scanners has revolutionized IHC data archiving and quantitative analysis. The fact that all of our slides were digitized using the Zeiss MIRAX Scan system allowed for digital archiving, thereby effectively eliminating any concerns regarding the long-term stability of the Acid Blue 129 counterstain. By digitizing the entire slides, qIHC can be performed over the entire tissue section, rather than being limited to selected ROIs. The large size of these images (e.g., 1010 pixels for a rat brain section) necessitates the use of fully automated qIHC methods. We have recently applied the method described in this study for whole brain qIHC analysis in transgenic mouse models of Alzheimer's disease (Chakravarty et al. 2007 Given the robustness of the algorithm for the two very different staining patterns produced by the Ki-67 IHC of rat C6 glioma and β-amyloid IHC of transgenic APP mouse brains in this particular study, we anticipate a broad applicability of this technique. Studies are currently underway to rigorously validate this technique using a wide range antibodies, tissues, and chromogens. We believe that the straightforward, inexpensive methods described for the preparation of high-contrast IHC sections and the fully automated segmentation algorithm presented in this study will allow for simple, robust qIHC studies in both research and clinical settings.
This work was supported by funds from the Montreal Neurological Institute. We thank Dr. Edith Hamel for kindly supplying the transgenic APP mouse brains and Kurt Hemmings for assisting in the preparation of the tissue sections and the manual segmentation data.
Received for publication November 23, 2007; accepted June 3, 2008
Abramoff MD, Magelhaes PJ, Ram SJ (2004) Image processing with ImageJ. Biophotonics Int 11:36–42 Bedell BJ, Narayana PA, Wolinsky JS (1997) A dual approach for minimizing false lesion classifications on magnetic resonance images. Magn Reson Med 37:94–102[Medline] Brey EM, King TW, Johnston C, McIntire LV, Reece GP, Patrick CW (2002) A technique for quantitative three-dimensional analysis of microvascular structure. Microvasc Res 63:279–294[CrossRef][Medline] Brey EM, Lalani Z, Johnston C, Wong M, McIntire LV, Duke PJ, Patrick CW (2003) Automated selection of DAB-labeled tissue for immunohistochemical quantification. J Histochem Cytochem 51:575–584 Chakravarty MM, Collins DL, Zehntner SP, Hemmings K, Chan C, Zijdenbos AP, Hamel E, et al. (2007) Ex vivo molecular imaging in transgenic mouse models of Alzheimer's disease. Alzheimer's Association International Conference on Prevention of Dementia, Washington, DC, June 9–13, 2007 Chakravarty MM, Sadikot AF, Germann J, Bertrand G, Collins DL (2005) Anatomical and electrophysiological validation of an atlas for neurosurgical planning. In Duncan J, Gerig G, eds. Eighth International Conference on Medical Image Computing and Computer Assisted Intervention MICCAI 2005, Volume 2 of Lecture Notes in Computer Science. Palm Springs, FL, Springer, 394–401 Del Duca D, Werbowetski T, Del Maestro RF (2004) Spheroid preparation from hanging drops: characterization of a model of brain tumor invasion. J Neurooncol 67:295–303[CrossRef][Medline] Kitamoto T, Ogomori K, Tateishi J, Prusiner SB (1987) Formic acid pretreatment enhances immunostaining of cerebral and systemic amyloids. Lab Invest 57:230–236[Medline] Mucke L, Masliah E, Yu GQ, Mallory M, Rockenstein EM, Tatsuno G, Hu K, et al. (2000) High-level neuronal expression of Aβ1–42 in wild-type human amyloid protein precursor Tg mice: synaptotoxicity without plaque formation. J Neurosci 20:4050–4058 Qi ZL, Huo X, Xu XL, Zhang B, Du MG, Yang HW, Zheng LK, et al. (2006) Relationship between HPV16/18 E6 and p53, p21WAF1, MDM2, Ki67 and cyclin D1 expression in esophageal squamous cell carcinoma: comparative study by using tissue microarray technology. Exp Oncol 28:235–240[Medline] Ruifrok AC, Johnston DA (2001) Quantification of histochemical staining by color deconvolution. Anal Quant Cytol Histol 23:291–299[Medline] Shi S-R, Liu C, Perez J, Taylor CR (2005) Protein-embedding technique: a potential approach to standardization of immunohistochemistry for formalin-fixed, paraffin-embedded tissue sections. J Histochem Cytochem 53:1167–1170
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||