J Bagh College Dentistry Vol. 29(1), March 2017 Computer assisted Oral Diagnosis 83 Computer Assisted Immunohistochemical Score Prediction Via Simplified Image Acquisition Technique Salam N. Jawad, B.D.S., M.Sc. (1) Bashar H. Abdullah, B.D.S., M.Sc., PhD. (2) ABSTRACT Background: techniques of image analysis have been used extensively to minimize interobserver variation of immunohistochemical scoring, yet; image acquisition procedures are often demanding, expensive and laborious. This study aims to assess the validity of image analysis to predict human observer’s score with a simplified image acquisition technique. Materials and methods: formalin fixed- paraffin embedded tissue sections for ameloblastomas and basal cell carcinomas were immunohistochemically stained with monoclonal antibodies to MMP-2 and MMP-9. The extent of antibody positivity was quantified using Imagej® based application on low power photomicrographs obtained with a conventional camera. Results of the software were employed to predict human visual scoring results with stepwise multiple regression analysis. Results: the overall prediction of epithelial score depicted as r square value was 0.26 (p<0.001) which was obviously higher than that of stromal score (0.10; p<0.01). Epithelial and stromal MMP-2 score prediction was generally higher than that of MMP-9. Collectively, ameloblastomas had a more efficient score prediction compared to basal cell carcinomas. Conclusion: there is a considerable variability in the prediction capacity of the technique with respect to different antibodies, different tumors and cellular versus stromal score. Keywords: immunohistochemistry, image analysis (J Bagh Coll Dentistry 2017; 29(1):83-88) INTRODUCTION Quantification of staining results is a crucial step within the methodology of thousands of immunohistochemical study designs; where attempts are always made to reach objective, concise and reproducible results. Semiquantitative scoring is the most common approach that involves one or more trained observer to convert what he/she sees into numeric or logical values(1). Advancements in computer hardware capabilities through the last few decades facilitated the use of image processing programs to collect data from captured images(1, 2, 3). One of the commonly employed multiplatform programs is “Imagej”; which is an open source software developed by NHS in 1997 and was involved in many microscopic research designs. The software had many built in or optional “add on” functions in the form of combined analytic steps merged into single click; they are called “plugins” and “macros” (4). Varghese et al. (2014) developed an Imagej® compatible macro that would separate and quantify diaminobenzidine(DAB) and Hematoxylin colors according to pixel intensity within images of immunohistochemically stained (1) PhD. student, Department of Oral Diagnosis, College of Dentistry, University of Baghdad. (2) Professor, Department of Oral Diagnosis, College of Dentistry, University of Baghdad. tissue sections(5). This study aims to evaluate the validity of this specific macro function with respect to human observer’s percentage scoring of two epithelial tumors (Ameloblastoma & Basal cell carcinoma) stained with two monoclonal antibodies (MMP-2, MMP-9) with a simplified image acquisition technique. MATERIALS AND METHODS Formalin fixed, paraffin embedded tissue specimens of ameloblastomas (AB) and basal cell carcinoma (BCC) were retrieved from the archives of oral pathology laboratory at the college of dentistry/Baghdad university and medical city laboratories. Five um thick sections were mounted on positively charged slides and stained immunohistochemically with monoclonal antibodies for MMP-2 (AB: n=32; BCC: n=26; ab2462,1:200) and MMP-9 (AB: n=33; BCC: n=29; ab76003,1:100) using EXPOSE Mouse and Rabbit Specific HRP/DAB Detection IHC kit (ab80436) . Stained tissue sections were evaluated by two observers with at least 5 high power fields (X40) for cellular and stromal positivity as a percentage of positive cells to the total high power field (HPF) then, a mean of the fields was calculated. Visual scoring values of 2 experienced observers were correlated by Pearson correlation coefficient, the values were accepted when R between the observers’ readings exceeds (0.80), J Bagh College Dentistry Vol. 29(1), March 2017 Computer assisted Oral Diagnosis 84 otherwise; the sections were re-evaluated to reach a consensus. At least, 3 photomicrographs for the DAB and hematoxylin stained tissue sections were captured with a 10X objective using an ordinary digital camera fixed to the eyepiece (Samsung-GT- N7100). With Imagej’s region of interest (ROI) selection options; multiple cellular areas were surrounded with selection polygons followed by IHC macro application (Fig. 1). For any given case, a mean value of each IHC profiler parameters result was calculated (Pixel count, percentage of high positive, percentage of positive, percentage of low positive, percentage of negative). The same preceding steps (selection and macro application) were applied to stromal areas only. Stepwise and multiple regression analysis for human score prediction was employed for statistical analysis using SPSS.22 software, adopting the Imagej’s output data as independent variables (predictors). . RESULTS Tumor epithelial and stromal expression was found in tissue sections of ameloblastoma and basal cell carcinoma for MMP-2 and MMP-9 to different extents. Combined visual and computer assisted immunohistochmical scoring characteristics of the two antibodies are summarized in tables 1 and 2. The predictive ability of the software algorithm, reflected by the (R square) varied across different tumors, antibodies and the evaluated compartments (epithelial or stromal). Tumor cellular evaluation with the software showed the highest R square values for the multiple stepwise regression in which the percentage of negative contribution (N%) and the pixel count (PC) were the most common independent variables that predict the human visual score. On the other hand, evidently lower R square values were found upon applying the software to the stromal elements, with more versatile predicting variables between different antibodies and different tumors. Figure 1: (A) Applying multiple polygon selection to regions of interest (ROI) of an image. (B) Choosing “IHC profiler macro” from the Imagej’s plugins drop down menu. (C) Software’s Results window A J Bagh College Dentistry Vol. 29(1), March 2017 Computer assisted Oral Diagnosis 85 DISCUSSION Reproducibility is a major issue in interpretation of immunohistochemical results quantification, both intra-observer and inter- observer agreements are a pre-requisite for any research to be valid (6). Several attempts have been made to neutralize human observer’s possible bias through image based computer assisted analytic approaches with variable success (1, 3, 7, 8). A research group at the Molecular Functional Imaging Laboratory in India developed a macro (a single computer instruction that expands automatically into a set of instructions to perform a particular task) named “IHC profiler” ; compatible with Imagej®, the macro would separate images of DAB stained immunohistochemical sections to 3 fractions of varied brown positivity and one more for blue negativity depending on individual pixel color values and display the results as a percentage contribution of high positive, positive, low positive and negative(5). In this study; the macro was applied to low power micrographs of immunohistochemically stained tissue sections of AB and BCC and the resulting output in addition to the number of pixels analyzed were plotted to predict the corresponding human observer’s results according to standard scoring methods. In order to be more time effective and simple for the non-experienced observer, this study had major design differences from the original work by Varghese et al. (5); the images were acquired with a conventional digital camera rather than a slide scanner in an attempt for method simplification to be accessible for any low budget laboratory, moreover; the study used a low magnification power for photomicrographs for a generalized view that would minimize the number of photographs needed for analysis. Ameloblastoma and basal cell carcinoma were selected in this study as prototypical examples for benign and malignant epithelial tumors that bear a morphological mimicry. Stepwise regression analysis is a special variant of multiple regression that enters the independent variables (macro results in this study) in the regression model one variable at a time with omission of the effect of other variables to specify the most significant variables that predict the value of the dependent variable (human scoring in this study)(9). At first glance, the predictive ability of IHC profiler in the cellular compartment was mostly associated with (N) variable which actually represents all the three macro variables (H, P and L) subtracted from 100%, however; pixel count had an effect in predicting AB human results more than BCC which is most probably resulting from the relatively less cellular islands selected for AB as compared to the BCC that leads to a lower pixel count. MMP-2 showed relatively high predictive results where the model explained a minimum of 40% of human scoring variability, a slightly higher R square value was noted in BCC in contrast to MMP-9 predictions which had notably lower values overall and a higher R square for AB. Regarding stromal expression, though generally less than cellular expression; R square values showed more predictive power in MMP-2 than MMP-9. Again, predictive values of MMP-2 human expression were higher in AB with conversely higher values for MMP-9 in BCC Those results suggest that the model is definitely better predicting epithelial MMP-2 scores than MMP-9 regardless of the tumor type, moreover; it was more predictive concerning MMP-9 in AB and MMP-2 in BCC. However, overall results of the study approach indicated marked variability of the software’s prediction capacity with respect to different antibodies and tumors. Keeping a minimally laborious and complicated image acquisition technique probably mandates the employment of more sophisticated computer analytic procedures to approach human visual perception. REFERENCES: 1- Heindl A, Nawaz S, Yuan Y. Mapping spatial heterogeneity in the tumor microenvironment: a new era for digital pathology. Lab Invest 2015; 95: 377- 384. 2- Kaczmarek E, Górna A, Majewski P. Techniques of image analysis for quantitative immunohistochemistry. RoczAkad Med Bialymst 2004; 49: 155-158. 3- Matos LLD, Stabenow E, Tavares MR, Ferraz AR, Capelozzi VL, Pinhal MS. Immunohistochemistry quantification by a digital computer-assisted method compared to semiquantitative analysis. Clinics 2006; 61: 417-424. 4- Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat methods 2012; 9: 671-675. 5- Varghese F, Bukhari AB, Malhotra R, De A. IHC Profiler: an open source plugin for the quantitative evaluation and automated scoring of immunohistochemistry images of human tissue samples. PloS one 2014; 9: e96801. 6- Zlobec I, Steele R, Michel RP, Compton CC,Lugli A, Jass JR. Scoring of p53, VEGF, Bcl-2 and APAF-1 immunohistochemistry and interobserver reliability in colorectal cancer. Mod pathol 2006; 19: 1236-1242. 7- Lehr HA, Mankoff DA, Corwin D, Santeusanio G, Gown AM. Application of photoshop-based image analysis to quantification of hormone receptor expression in breast cancer. J HistochemCytochem1997; 45: 1559-1565. J Bagh College Dentistry Vol. 29(1), March 2017 Computer assisted Oral Diagnosis 86 8- Laurinavicius A,Laurinaviciene A, Ostapenko V, Dasevicius D, Jarmalaite S, LazutkaJ. Immunohistochemistry profiles of breast ductal carcinoma: factor analysis of digital image analysis data. DiagnPathol 2012; 7: 27. 9- Dejmek A,Hjerpe A. Immunohistochemical reactivity in mesothelioma and adenocarcinoma: a stepwise logistic regression analysis. Apmis 1994; 102: 255- 264. الخالصة في التفسير تباين غرض تجنب الاستعملت تقنيات التحليل الصوري باستخدام الكومبيوتر بشكل مسهب في عدة دراسات سابقة ل الخلفية: نتائج الصبغات المناعية عند تعدد الفاحصين, مع ذلك؛ امتازت الطرق المختلفة بكونها تتطلب عمليات معقدة ومكلفة اللتقاط البصري ل عدة للتحليل. تهدف هذة الدراسة لفحص امكانية توقع نتائج الفحص البصري بناءا على التحليل الصوري لصور شرائح الصور الم مجهرية ملتقطة بطريقة مبسطة. ومي الَ َوَرم لل رة بشمع البارافين: تضمنت الدراسة عينات مثبتة بالفورمالين ومطمومواد وطرق العمل )اميلوبالستوما( ومينائي ال رم بطريقة التصبيغ الكيميائي MMP-9و MMP-2سرطان الخلية القاعدية صبغت شرائحها باجسام مضادة احادية النسل موجهة الى اجريت والتي ®Imagejالمناعي. تم تقييم نتائج االصطباغ المناعي بالطرق البصرية المعتادة و بوساطة التحليل الصوري لبرنامج ملتقطة باستخدام كاميرا اعتيادية من ثم قورنت نتائج التحليل الصوري مع نتائج الفحص البصري على صور مجهرية واطئة التكبير المجهري عن طريق احتساب صحة توقع التحليل الصوري باسلوب االنحدار الخطي المتسلسل. R)2 ;0.10=مقارنة باالجزاء السدوية منهاR =0.26 (p<0.001)2 ظهرت نتائج اعلى للتوقع في اجزاء االورام الظهارية النتائج: p<0.01) وبشكل عام؛ كانت نتائج التوقع لالجزاء الظهارية والسدوية لمؤشر ,MMP-2 اعلى من نتائج مؤشرMMP-9 كذلك؛ , مقارنة بسرطان الخلية القاعدية.)االميلوبالستوما( ورم قيمتها اعلى فيكانت داخل تعددةوالمناطق الم المختلفة المؤشرات واالورام بين بشكل ملموس ليلي للطريقة المستخدمةنتائج التوقع التح تتبايناالستتنتاجات: االورام. J Bagh College Dentistry Vol. 29(1), March 2017 Computer assisted Oral Diagnosis 87 J Bagh College Dentistry Vol. 29(1), March 2017 Computer assisted Oral Diagnosis 88