Bimodality of intratumor Ki67 expression is an independent prognostic factor of overall survival in patients with invasive breast carcinoma

Advancing Ki67 hotspot detection in breast cancer: a comparative analysis of automated digital image analysis algorithms

AIM

Manual detection and scoring of Ki67 hotspots is difficult and prone to variability, limiting its clinical utility. Automated hotspot detection and scoring by digital image analysis (DIA) could improve the assessment of the Ki67 hotspot proliferation index (PI). This study compared the clinical performance of Ki67 hotspot detection and scoring DIA algorithms based on virtual dual staining (VDS) and deep learning (DL) with manual Ki67 hotspot PI assessment.

METHODS

Tissue sections of 135 consecutive invasive breast carcinomas were immunohistochemically stained for Ki67. Two DIA algorithms, based on VDS and DL, automatically determined the Ki67 hotspot PI. For manual assessment; two independent observers detected hotspots and calculated scores using a validated scoring protocol.

RESULTS

Automated hotspot detection and assessment by VDS and DL could be performed in 73% and 100% of the cases, respectively. Automated hotspot detection by VDS and DL led to higher Ki67 hotspot PIs (mean 39.6% and 38.3%, respectively) compared to manual consensus Ki67 PIs (mean 28.8%). Comparing manual consensus Ki67 PIs with VDS Ki67 PIs revealed substantial correlation (r = 0.90), while manual consensus versus DL Ki67 PIs demonstrated high correlation (r = 0.95).

CONCLUSION

Automated Ki67 hotspot detection and analysis correlated strongly with manual Ki67 assessment and provided higher PIs compared to manual assessment. The DL-based algorithm outperformed the VDS-based algorithm in clinical applicability, because it did not depend on virtual alignment of slides and correlated stronger with manual scores. Use of a DL-based algorithm may allow clearer Ki67 PI cutoff values, thereby improving the clinical usability of Ki67.

Spatial distributions of CD8 and Ki67 cells in the tumor microenvironment independently predict breast cancer-specific survival in patients with ER+HER2- and triple-negative breast carcinoma

INTRODUCTION

Breast cancer (BC) presents diverse malignancies with varying biological and clinical behaviors, driven by an interplay between cancer cells and tumor microenvironment. Deciphering these interactions is crucial for personalized diagnostics and treatment. This study explores the prognostic impact of tumor proliferation and immune response patterns, assessed by computational pathology indicators, on breast cancer-specific survival (BCSS) models in estrogen receptor-positive HER2-negative (ER+HER2-) and triple-negative BC (TNBC) patients.

MATERIALS AND METHODS

Whole-slide images of tumor surgical excision samples from 252 ER+HER2- patients and 63 TNBC patients stained for estrogen and progesterone receptors, Ki67, HER2, and CD8 were analyzed. Digital image analysis (DIA) was performed for tumor tissue segmentation and quantification of immunohistochemistry (IHC) markers; the DIA outputs were subsampled by hexagonal grids to assess the spatial distributions of Ki67-positive tumor cells and CD8-positive (CD8+) cell infiltrates, expressed as Ki67-entropy and CD8-immunogradient indicators, respectively. Prognostic models for BCSS were generated using multivariable Cox regression analysis, integrating clinicopathological and computational IHC indicators.

RESULTS

In the ER+HER2- BC, multivariable Cox regression revealed that high CD8+ density within the tumor interface zone (IZ) (HR: 0.26, p = 0.0056), low immunodrop indicator of CD8+ density (HR: 2.93, p = 0.0051), and low Ki67-entropy (HR: 5.95, p = 0.0.0061) were independent predictors of better BCSS, while lymph node involvement predicted worse BCSS (HR: 3.30, p = 0.0013). In TNBC, increased CD8+ density in the IZ stroma (HR: 0.19, p = 0.0119) and Ki67-entropy (HR: 3.31, p = 0.0250) were independent predictors of worse BCSS. Combining these independent indicators enhanced prognostic stratification in both BC subtypes.

CONCLUSIONS

Computational biomarkers, representing spatial properties of the tumor proliferation and immune cell infiltrates, provided independent prognostic information beyond conventional IHC markers in BC. Integrating Ki67-entropy and CD8-immunogradient indicators into prognostic models can improve patient stratification with regard to BCSS.

Characterization of Breast Cancer Intra-Tumor Heterogeneity Using Artificial Intelligence

Intra-tumor heterogeneity (ITH) is a fundamental characteristic of breast cancer (BC), influencing tumor progression, prognosis, and therapeutic responses. However, the complexity of ITH in BC makes its accurate characterization challenging. This study leverages deep learning (DL) techniques to comprehensively evaluate ITH in early-stage luminal BC and provide a nuanced understanding of its impact on tumor behavior and patient outcomes. A large cohort (n = 2561) of early-stage luminal BC was evaluated using whole slide images (WSIs) of hematoxylin and eosin-stained slides of excision specimens. Morphological features of both the tumor and stromal components were meticulously annotated by a panel of pathologists in a subset of cases. A DL model was applied to develop an algorithm to assess the degree of heterogeneity of various morphological features per individual case utilizing defined patches. The results of extracted features were used to generate an overall heterogeneity score that was correlated with the clinicopathological features and outcome. Overall, 162 features were quantified and a significant positive correlation between these features was identified. Specifically, there was a significant association between a high degree of intra-tumor heterogeneity and larger tumor size, poorly differentiated tumors, highly proliferative tumors, tumors of no special type (NST), and those with low estrogen receptor (ER) expression. When all features are considered in combination, a high overall heterogeneity score was significantly associated with parameters characteristic of aggressive tumor behavior, and it was an independent predictor of poor patient outcome. In conclusion, DL models can be used to accurately decipher the complexity of ITH and provide extra information for outcome prediction.

Intratumoral heterogeneity of Ki67 proliferation index outperforms conventional immunohistochemistry prognostic factors in estrogen receptor-positive HER2-negative breast cancer

In breast cancer (BC), pathologists visually score ER, PR, HER2, and Ki67 biomarkers to assess tumor properties and predict patient outcomes. This does not systematically account for intratumoral heterogeneity (ITH) which has been reported to provide prognostic value. This study utilized digital image analysis (DIA) and computational pathology methods to investigate the prognostic value of ITH indicators in ER-positive (ER+) HER2-negative (HER2-) BC patients. Whole slide images (WSIs) of surgically excised specimens stained for ER, PR, Ki67, and HER2 from 254 patients were used. DIA with tumor tissue segmentation and detection of biomarker-positive cells was performed. The DIA-generated data were subsampled by a hexagonal grid to compute Haralick's texture indicators for ER, PR, and Ki67. Cox regression analyses were performed to assess the prognostic significance of the immunohistochemistry (IHC) and ITH indicators in the context of clinicopathologic variables. In multivariable analysis, the ITH of Ki67-positive cells, measured by Haralick's texture entropy, emerged as an independent predictor of worse BC-specific survival (BCSS) (hazard ratio (HR) = 2.64, p-value = 0.0049), along with lymph node involvement (HR = 2.26, p-value = 0.0195). Remarkably, the entropy representing the spatial disarrangement of tumor proliferation outperformed the proliferation rate per se established either by pathology reports or DIA. We conclude that the Ki67 entropy indicator enables a more comprehensive risk assessment with regard to BCSS, especially in cases with borderline Ki67 proliferation rates. The study further demonstrates the benefits of high-capacity DIA-generated data for quantifying the essentially subvisual ITH properties.

Quantitative expression of oestrogen receptor in breast cancer: Clinical and molecular significance

BACKGROUND

Oestrogen receptor (ER) positive breast cancer (BC) patients are eligible for endocrine therapy (ET), regardless of ER immunohistochemical expression level. There is a wide spectrum of ER expression and the response to ET is not uniform. This study aimed to assess the clinical and molecular consequences of ER heterogeneity with respect to ET-response.

METHODS

ER expression, categorised by percentage and staining intensity in a large BC cohort (n = 7559) was correlated with clinicopathological parameters and patient ET response. The Cancer Genome Atlas Data BC cohort (n = 1047) was stratified by ER expression and transcriptomic analysis completed to better understand the molecular basis of ER heterogeneity.

RESULTS

The quantitative proportional increase in ER expression was positively associated with favourable prognostic parameters. Tumours with 1-9% ER expression were characteristically similar to ER-negative (<1%) tumours. Maximum ET-response was observed in tumours with 100% ER expression, with responses significantly different to tumours exhibiting ER at < 100% and significantly decreased survival rates were observed in tumours with 50% and 10% of ER expression. The Histochemical-score (H-score), which considers both staining intensity and percentage, added significant prognostic value over ER percentage alone with significant outcome differences observed at H-scores of 30, 100 and 200. There was a positive correlation between ER expression and ESR1 mRNA expression and expression of ER-regulated genes. Pathway analysis identified differential expression in key cancer-related pathways in different ER-positive groups.

CONCLUSION

ET-response is statistically proportionally related to ER expression with significant differences observed at 10%, 50% and 100%. The H-score adds prognostic and predictive information.

Digital Image Analysis of Ki67 Heterogeneity Improves the Diagnosis and Prognosis of Gastroenteropancreatic Neuroendocrine Neoplasms

Ki67 is a reliable grading and prognostic biomarker of gastroenteropancreatic neuroendocrine neoplasms (GEP-NENs). The intratumor heterogeneity of Ki67, correlated with tumor progression, is a valuable factor that requires image analysis. The application of digital image analysis (DIA) enables new approaches for the assessment of Ki67 heterogeneity distribution. We investigated the diagnostic utility of Ki67 heterogeneity parameters in the classification and grading of GEP-NENs and explored their clinical values with regard to their prognostic relevance. The DIA algorithm was performed on whole-slide images of 102 resection samples with Ki67 staining. Good agreement was observed between the manual and DIA methods in the hotspot evaluation (R² = 0.94, P < .01). Using the grid-based region of interest approach, score-based heat maps provided a distinctive overview of the intratumoral distribution of Ki67 between neuroendocrine carcinomas and neuroendocrine tumors. The computation of heterogeneity parameters related to DIA-determined Ki67 showed that the coefficient of variation and Morisita-Horn index were directly related to the classification and grading of GEP-NENs and provided insights into distinguishing high-grade neuroendocrine neoplasms (grade 3 neuroendocrine tumor vs neuroendocrine carcinoma, P < .01). Our study showed that a high Morisita-Horn index correlated with poor disease-free survival (multivariate analysis: hazard ratio, 56.69), which was found to be the only independent predictor of disease-free survival in patients with GEP-NEN. These spatial biomarkers have an impact on the classification and grading of tumors and highlight the prognostic associations of tumor heterogeneity. Digitization of Ki67 variations provides a direct and objective measurement of tumor heterogeneity and better predicts the biological behavior of GEP-NENs.

Clinical implication of low estrogen receptor (ER-low) expression in breast cancer

Breast cancer is a heterogeneous disease, and the estrogen receptor (ER) remains the most important biomarker in breast oncology. Most guidelines set a positive expression threshold of 1% staining in immunohistochemistry (IHC) to define ER positivity. However, different expression levels may be associated with diverse degrees of sensitivity to endocrine therapy as ER expression may impact breast cancer molecular biology as a continuous variable. ER-lo tumors, defined as those with 1-10% ER expression, represent a relatively small subgroup of breast cancer patients, with an estimated prevalence of 2-7%. These tumors are similar to ERneg disease in their molecular landscape, clinicopathological characteristics, prognosis, and response to therapy. Nevertheless, a proportion may retain some degree of ER signaling dependency, and the possibility of responding to some degree to endocrine therapy cannot be completely ruled out. This review article discusses the most important considerations regarding the definition of ER positivity, pathology assessment, prognosis, and therapeutic implication of ERlo breast cancer from the medical oncology perspective.

Clinical validity and clinical utility of Ki67 in early breast cancer

Ki67 represents an immunohistochemical nuclear localized marker that is widely used in surgical pathology. Nuclear immunoreactivity for Ki67 indicates that cells are cycling and are in G1- to S-phase. The percentage of Ki67-positive tumor cells (Ki67 index) therefore provides an estimate of the growth fraction in tumor specimens. In breast cancer (BC), tumor cell proliferation rate is one of the most relevant prognostic markers and Ki67 is consequently helpful in prognostication similar to histological grading and mRNA profiling-based BC risk stratification. In BCs treated with short-term preoperative endocrine therapy, Ki67 dynamics enable distinguishing between endocrine sensitive and resistant tumors. Despite its nearly universal use in pathology laboratories worldwide, no internationally accepted consensus has yet been achieved for some methodological details related to Ki67 immunohistochemistry (IHC). Controversial issues refer to choice of IHC antibody clones, scoring methods, inter-laboratory reproducibility, and the potential value of computer-assisted imaging analysis and/or artificial intelligence for Ki67 assessment. Prospective clinical trials focusing on BC treatment have proven that Ki67, as determined by standardized central pathology assessment, is of clinical validity. Clinical utility has been demonstrated in huge observational studies.

The state of the art for artificial intelligence in lung digital pathology

Lung diseases carry a significant burden of morbidity and mortality worldwide. The advent of digital pathology (DP) and an increase in computational power have led to the development of artificial intelligence (AI)-based tools that can assist pathologists and pulmonologists in improving clinical workflow and patient management. While previous works have explored the advances in computational approaches for breast, prostate, and head and neck cancers, there has been a growing interest in applying these technologies to lung diseases as well. The application of AI tools on radiology images for better characterization of indeterminate lung nodules, fibrotic lung disease, and lung cancer risk stratification has been well documented. In this article, we discuss methodologies used to build AI tools in lung DP, describing the various hand-crafted and deep learning-based unsupervised feature approaches. Next, we review AI tools across a wide spectrum of lung diseases including cancer, tuberculosis, idiopathic pulmonary fibrosis, and COVID-19. We discuss the utility of novel imaging biomarkers for different types of clinical problems including quantification of biomarkers like PD-L1, lung disease diagnosis, risk stratification, and prediction of response to treatments such as immune checkpoint inhibitors. We also look briefly at some emerging applications of AI tools in lung DP such as multimodal data analysis, 3D pathology, and transplant rejection. Lastly, we discuss the future of DP-based AI tools, describing the challenges with regulatory approval, developing reimbursement models, planning clinical deployment, and addressing AI biases. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Intratumoral Heterogeneity and Immune Response Indicators to Predict Overall Survival in a Retrospective Study of HER2-Borderline (IHC 2+) Breast Cancer Patients

Breast cancer (BC) categorized as human epidermal growth factor receptor 2 (HER2) borderline [2+ by immunohistochemistry (IHC 2+)] presents challenges for the testing, frequently obscured by intratumoral heterogeneity (ITH). This leads to difficulties in therapy decisions. We aimed to establish prognostic models of overall survival (OS) of these patients, which take into account spatial aspects of ITH and tumor microenvironment by using hexagonal tiling analytics of digital image analysis (DIA). In particular, we assessed the prognostic value of Immunogradient indicators at the tumor–stroma interface zone (IZ) as a feature of antitumor immune response. Surgical excision samples stained for estrogen receptor (ER), progesterone receptor (PR), Ki67, HER2, and CD8 from 275 patients with HER2 IHC 2+ invasive ductal BC were used in the study. DIA outputs were subsampled by HexT for ITH quantification and tumor microenvironment extraction for Immunogradient indicators. Multiple Cox regression revealed HER2 membrane completeness (HER2 MC) (HR: 0.18, p = 0.0007), its spatial entropy (HR: 0.37, p = 0.0341), and ER contrast (HR: 0.21, p = 0.0449) as independent predictors of better OS, with worse OS predicted by pT status (HR: 6.04, p = 0.0014) in the HER2 non-amplified patients. In the HER2-amplified patients, HER2 MC contrast (HR: 0.35, p = 0.0367) and CEP17 copy number (HR: 0.19, p = 0.0035) were independent predictors of better OS along with worse OS predicted by pN status (HR: 4.75, p = 0.0018). In the non-amplified tumors, three Immunogradient indicators provided the independent prognostic value: CD8 density in the tumor aspect of the IZ and CD8 center of mass were associated with better OS (HR: 0.23, p = 0.0079 and 0.14, p = 0.0014, respectively), and CD8 density variance along the tumor edge predicted worse OS (HR: 9.45, p = 0.0002). Combining these three computational indicators of the CD8 cell spatial distribution within the tumor microenvironment augmented prognostic stratification of the patients. In the HER2-amplified group, CD8 cell density in the tumor aspect of the IZ was the only independent immune response feature to predict better OS (HR: 0.22, p = 0.0047). In conclusion, we present novel prognostic models, based on computational ITH and Immunogradient indicators of the IHC biomarkers, in HER2 IHC 2+ BC patients.

Artificial intelligence applied to breast pathology

The convergence of digital pathology and computer vision is increasingly enabling computers to perform tasks performed by humans. As a result, artificial intelligence (AI) is having an astoundingly positive effect on the field of pathology, including breast pathology. Research using machine learning and the development of algorithms that learn patterns from labeled digital data based on "deep learning" neural networks and feature-engineered approaches to analyze histology images have recently provided promising results. Thus far, image analysis and more complex AI-based tools have demonstrated excellent success performing tasks such as the quantification of breast biomarkers and Ki67, mitosis detection, lymph node metastasis recognition, tissue segmentation for diagnosing breast carcinoma, prognostication, computational assessment of tumor-infiltrating lymphocytes, and prediction of molecular expression as well as treatment response and benefit of therapy from routine H&E images. This review critically examines the literature regarding these applications of AI in the area of breast pathology.

Machine-Learning-Based Evaluation of Intratumoral Heterogeneity and Tumor-Stroma Interface for Clinical Guidance

Assessment of intratumoral heterogeneity and tumor-host interaction within the tumor microenvironment is becoming increasingly important for innovative cancer therapy decisions because of the unique information it can generate about the state of the disease. However, its assessment and quantification are limited by ambiguous definitions of the tumor-host interface and by human cognitive capacity in current pathology practice. Advances in machine learning and artificial intelligence have opened the field of digital pathology to novel tissue image analytics and feature extraction for generation of high-capacity computational disease management models. A particular benefit is expected from machine-learning applications that can perform extraction and quantification of subvisual features of both intratumoral heterogeneity and tumor microenvironment aspects. These methods generate information about cancer cell subpopulation heterogeneity, potential tumor-host interactions, and tissue microarchitecture, derived from morphologically resolved content using both explicit and implicit features. Several studies have achieved promising diagnostic, prognostic, and predictive artificial intelligence models that often outperform current clinical and pathology criteria. However, further effort is needed for clinical adoption of such methods through development of standardizable high-capacity workflows and proper validation studies.

Assessment of proliferation in breast cancer: cell cycle or mitosis? An observational study

BACKGROUND AND AIMS

Proliferation is an important indicator of breast cancer (BC) prognosis, but is assessed using different approaches. Not all cells in the cell cycle are committed to division. This study aimed to characterise quantitative differences between BC cells in the cell cycle and those in mitosis and assess their relationship with other pathological parameters.

METHODS AND RESULTS

A cohort of BC sections (n = 621) was stained with haematoxylin and eosin and immunohistochemistry for Ki-67. The proportion of mitotic cells and Ki-67-positive cells was assessed in the same areas. The Cancer Genome Atlas (TCGA) BC cohort was used to assess MKI-67 transcriptome level and its association with the mitotic counts. The mean proportion of BC cells in the cell cycle was 24% (range = 1-90%), while the mean proportion of BC cells in mitosis was 5% (range = 0-73%). A low proportion of mitoses to whole cycling cells was associated with low histological grade tumours and the luminal A molecular subtype, while tumours with a high proportion of mitoses to the overall cycling cells were associated with triple-negative subtype, larger tumour size, grade 3 tumours and lymph node metastasis. The high mitosis/low Ki-67-positive cells tumours showed a significant association with variables of poor prognosis, including high-grade and triple-negative subtypes.

CONCLUSION

The proportion of BC cells in the cell cycle and mitosis is variable. We show that not only the number of cells in the cell cycle or mitosis, but also the difference between them, provides valuable information on tumour aggressiveness.

How the variability between computer-assisted analysis procedures evaluating immune markers can influence patients' outcome prediction

Differences between computer-assisted image analysis (CAI) algorithms may cause discrepancies in the identification of immunohistochemically stained immune biomarkers in biopsies of breast cancer patients. These discrepancies have implications for their association with disease outcome. This study aims to compare three CAI procedures (A, B and C) to measure positive marker areas in post-neoadjuvant chemotherapy biopsies of patients with triple-negative breast cancer (TNBC) and to explore the differences in their performance in determining the potential association with relapse in these patients. A total of 3304 digital images of biopsy tissue obtained from 118 TNBC patients were stained for seven immune markers using immunohistochemistry (CD4, CD8, FOXP3, CD21, CD1a, CD83, HLA-DR) and were analyzed with procedures A, B and C. The three methods measure the positive pixel markers in the total tissue areas. The extent of agreement between paired CAI procedures, a principal component analysis (PCA) and Cox multivariate analysis was assessed. Comparisons of paired procedures showed close agreement for most of the immune markers at low concentration. The probability of differences between the paired procedures B/C and B/A was generally higher than those observed in C/A. The principal component analysis, largely based on data from CD8, CD1a and HLA-DR, identified two groups of patients with a significantly lower probability of relapse than the others. The multivariate regression models showed similarities in the factors associated with relapse for procedures A and C, as opposed to those obtained with procedure B. General agreement among the results of CAI procedures would not guarantee that the same predictive breast cancer markers were consistently identified. These results highlight the importance of developing additional strategies to improve the sensitivity of CAI procedures.

Detection of Ki67 Hot-Spots of Invasive Breast Cancer Based on Convolutional Neural Networks Applied to Mutual Information of H&E and Ki67 Whole Slide Images

Ki67 hot-spot detection and its evaluation in invasive breast cancer regions play a significant role in routine medical practice. The quantification of cellular proliferation assessed by Ki67 immunohistochemistry is an established prognostic and predictive biomarker that determines the choice of therapeutic protocols. In this paper, we present three deep learning-based approaches to automatically detect and quantify Ki67 hot-spot areas by means of the Ki67 labeling index. To this end, a dataset composed of 100 whole slide images (WSIs) belonging to 50 breast cancer cases (Ki67 and H&E WSI pairs) was used. Three methods based on CNN classification were proposed and compared to create the tumor proliferation map. The best results were obtained by applying the CNN to the mutual information acquired from the color deconvolution of both the Ki67 marker and the H&E WSIs. The overall accuracy of this approach was 95%. The agreement between the automatic Ki67 scoring and the manual analysis is promising with a Spearman’s ρ correlation of 0.92. The results illustrate the suitability of this CNN-based approach for detecting hot-spots areas of invasive breast cancer in WSI.

Spatial transcriptomics inferred from pathology whole-slide images links tumor heterogeneity to survival in breast and lung cancer

Digital analysis of pathology whole-slide images is fast becoming a game changer in cancer diagnosis and treatment. Specifically, deep learning methods have shown great potential to support pathology analysis, with recent studies identifying molecular traits that were not previously recognized in pathology H&E whole-slide images. Simultaneous to these developments, it is becoming increasingly evident that tumor heterogeneity is an important determinant of cancer prognosis and susceptibility to treatment, and should therefore play a role in the evolving practices of matching treatment protocols to patients. State of the art diagnostic procedures, however, do not provide automated methods for characterizing and/or quantifying tumor heterogeneity, certainly not in a spatial context. Further, existing methods for analyzing pathology whole-slide images from bulk measurements require many training samples and complex pipelines. Our work addresses these two challenges. First, we train deep learning models to spatially resolve bulk mRNA and miRNA expression levels on pathology whole-slide images (WSIs). Our models reach up to 0.95 AUC on held-out test sets from two cancer cohorts using a simple training pipeline and a small number of training samples. Using the inferred gene expression levels, we further develop a method to spatially characterize tumor heterogeneity. Specifically, we produce tumor molecular cartographies and heterogeneity maps of WSIs and formulate a heterogeneity index (HTI) that quantifies the level of heterogeneity within these maps. Applying our methods to breast and lung cancer slides, we show a significant statistical link between heterogeneity and survival. Our methods potentially open a new and accessible approach to investigating tumor heterogeneity and other spatial molecular properties and their link to clinical characteristics, including treatment susceptibility and survival.

Report From the International Society of Urological Pathology (ISUP) Consultation Conference on Molecular Pathology of Urogenital Cancers. I. Molecular Biomarkers in Prostate Cancer

The combined clinical and molecular heterogeneity of prostate cancer necessitates the use of prognostic, predictive, and diagnostic biomarkers to assist the clinician with treatment selection. The pathologist plays a critical role in guiding molecular biomarker testing in prostate cancer and requires a thorough knowledge of the current testing options. In the setting of clinically localized prostate cancer, prognostic biomarkers such as Ki-67 labeling, PTEN loss or mRNA-based genomic signatures can be useful to help determine whether definitive therapy is required. In the setting of advanced disease, predictive biomarkers, such as the presence of DNA repair deficiency mediated by BRCA2 loss or mismatch repair gene defects, may suggest the utility of poly-ADP ribosylase inhibition or immune checkpoint blockade. Finally, androgen receptor-related biomarkers or diagnostic biomarkers indicating the presence of small cell neuroendocrine prostate cancer may help guide the use of androgen receptor signaling inhibitors and chemotherapy. In this review, we examine the current evidence for several prognostic, predictive and diagnostic tissue-based molecular biomarkers in prostate cancer management. For each assay, we summarize a recent survey of the International Society of Urology Pathology (ISUP) members on current testing practices and include recommendations for testing that emerged from the ISUP Working Group on Molecular Pathology of Prostate Cancer and the 2019 Consultation Conference on Molecular Pathology of Urogenital Cancers.

Factors influencing agreement of breast cancer luminal molecular subtype by Ki67 labeling index between core needle biopsy and surgical resection specimens

Reliable determination of Ki67 labeling index (Ki67-LI) on core needle biopsy (CNB) is essential for determining breast cancer molecular subtype for therapy planning. However, studies on agreement between molecular subtype and Ki67-LI between CNB and surgical resection (SR) specimens are conflicting. The present study analyzed the influence of clinicopathological and sampling-associated factors on agreement. Molecular subtype was determined visually by Ki67-LI in 484 pairs of CNB and SR specimens of invasive estrogen receptor (ER)-positive, human epidermal growth factor (HER2)-negative breast cancer. Luminal B disease was defined by Ki67-LI > 20% in SR. Correlation of molecular subtype agreement with age, menopausal status, CNB method, Breast Imaging Reporting and Data System imaging category, time between biopsies, type of surgery, and pathological tumor parameters was analyzed. Recurrence-free survival (RFS) and overall survival (OS) were analyzed using the Kaplan-Meier method. CNB had a sensitivity of 77.95% and a specificity of 80.97% for identifying luminal B tumors in CNB, compared with the final molecular subtype determination after surgery. The correlation of Ki67-LI between CNB and SR was moderate (ROC-AUC 0.8333). Specificity and sensitivity for CNB to correctly define molecular subtype of tumors according to SR were significantly associated with tumor grade, immunohistochemical progesterone receptor (PR) and p53 expression (p < 0.05). Agreement of molecular subtype did not significantly impact RFS and OS (p = 0.22 for both). The identified factors likely mirror intratumoral heterogeneity that might compromise obtaining a representative CNB. Our results challenge the robustness of a single CNB-driven measurement of Ki67-LI to identify luminal B breast cancer of low (G1) or intermediate (G2) grade.

Independent Prognostic Value of Intratumoral Heterogeneity and Immune Response Features by Automated Digital Immunohistochemistry Analysis in Early Hormone Receptor-Positive Breast Carcinoma

Immunohistochemistry (IHC) for ER, PR, HER2, and Ki67 is used to predict outcome and therapy response in breast cancer patients. The current IHC assessment, visual or digital, is based mostly on global biomarker expression levels in the tissue sample. In our study, we explored the prognostic value of digital image analysis of conventional breast cancer IHC biomarkers supplemented with their intratumoral heterogeneity and tissue immune response indicators. Surgically excised tumor samples from 101 female patients with hormone receptor-positive breast cancer (HRBC) were stained for ER, PR, HER2, Ki67, SATB1, CD8, and scanned at 20x. Digital image analysis was performed using the HALO™ platform. Subsequently, hexagonal tiling was used to compute intratumoral heterogeneity indicators for ER, PR and Ki67 expression. Multiple Cox regression analysis revealed three independent predictors of the patient's overall survival: Haralick's texture entropy of PR (HR = 0.19, p = 0.0005), Ki67 Ashman's D bimodality (HR = 3.0, p = 0.01), and CD8+SATB1+ cell density in tumor tissue (HR = 0.32, p = 0.02). Remarkably, the PR and Ki67 intratumoral heterogeneity indicators were prognostically more informative than the rates of their expression. In particular, a distinct non-linear relationship between the rate of PR expression and its intratumoral heterogeneity was observed and revealed a non-linear prognostic effect of PR expression. The independent prognostic significance of CD8+SATB1+ cells infiltrating the tumor could indicate their role in anti-tumor immunity. In conclusion, we suggest that prognostic modeling, based entirely on the computational image-based IHC biomarkers, is possible in HRBC patients. The intratumoral heterogeneity and immune response indicators outperformed both conventional breast cancer IHC and clinicopathological variables while markedly increasing the power of the model.

Practical approaches to automated digital image analysis of Ki-67 labeling index in 997 breast carcinomas and causes of discordance with visual assessment

The Ki-67 labeling index (LI) is an important prognostic factor in breast carcinoma. The Ki-67 LI is traditionally calculated via unaided microscopic estimation; however, inter-observer and intra-observer variability and low reproducibility are problems with this visual assessment (VA) method. For more accurate assessment and better reproducibility with Ki-67 LI, digital image analysis was introduced recently. We used both VA and automated digital image analysis (ADIA) (Ventana Virtuoso image management software) to estimate Ki-67 LI for 997 cases of breast carcinoma, and compared VA and ADIA results. VA and ADIA were highly correlated (intraclass correlation coefficient 0.982, and Spearman’s correlation coefficient 0.966, p<0.05). We retrospectively analyzed cases with a greater than 5% difference between VA and ADIA results. The cause of these differences was: (1) tumor heterogeneity (98 cases, 56.0%), (2) VA interpretation error (32 cases, 18.3%), (3) misidentification of tumor cells (26 cases, 14.9%), (4) poor immunostaining or slide quality (16 cases, 9.1%), and (5) Estimation of non-tumor cells (3 cases, 1.7%). There were more discrepancies between VA and ADIA results in the group with a VA value of 10–20% compared to groups with <10% and ≥20%. Although ADIA is more accurate than VA, there are some limitations. Therefore, ADIA findings require confirmation by a pathologist.

Ancillary Prognostic and Predictive Testing in Breast Cancer: Focus on Discordant, Unusual, and Borderline Results

Ancillary testing in breast cancer has become standard of care to determine what therapies may be most effective for individual patients with breast cancer. Single-marker tests are required on all newly diagnosed and newly metastatic breast cancers. Markers of proliferation are also used, and include both single-marker tests like Ki67 as well as panel-based gene expression tests, which have made more recent contributions to prognostic and predictive testing in breast cancers. This review focuses on pathologist interpretation of these ancillary test results, with a focus on expected versus unexpected results and troubleshooting borderline, unusual, or discordant results.

Programmed death-ligand 1 expression by digital image analysis advances thyroid cancer diagnosis among encapsulated follicular lesions

Recognition of noninvasive follicular thyroid neoplasms with papillary-like nuclear features (NIFTP) that distinguishes them from invasive malignant encapsulated follicular variant of papillary thyroid carcinoma (EFVPTC) can prevent overtreatment of NIFTP patients. We and others have previously reported that programmed death-ligand 1 (PD-L1) is a useful biomarker in thyroid tumors; however, all reports to date have relied on manual scoring that is time consuming as well as subject to individual bias. Consequently, we developed a digital image analysis (DIA) protocol for cytoplasmic and membranous stain quantitation (ThyApp) and evaluated three tumor sampling methods [Systemic Uniform Random Sampling, hotspot nucleus, and hotspot nucleus/3,3'-Diaminobenzidine (DAB)]. A patient cohort of 153 cases consisting of 48 NIFTP, 44 EFVPTC, 26 benign nodules and 35 encapsulated follicular lesions/neoplasms with lymphocytic thyroiditis (LT) was studied. ThyApp quantitation of PD-L1 expression revealed a significant difference between invasive EFVPTC and NIFTP; but none between NIFTP and benign nodules. ThyApp integrated with hotspot nucleus tumor sampling method demonstrated to be most clinically relevant, consumed least processing time, and eliminated interobserver variance. In conclusion, the fully automatic DIA algorithm developed using a histomorphological approach objectively quantitated PD-L1 expression in encapsulated thyroid neoplasms and outperformed manual scoring in reproducibility and higher efficiency.

Digital image analysis of Ki67 proliferation index in breast cancer using virtual dual staining on whole tissue sections: clinical validation and inter-platform agreement

Purpose

The Ki67 proliferation index is a prognostic and predictive marker in breast cancer. Manual scoring is prone to inter- and intra-observer variability. The aims of this study were to clinically validate digital image analysis (DIA) of Ki67 using virtual dual staining (VDS) on whole tissue sections and to assess inter-platform agreement between two independent DIA platforms.

Methods

Serial whole tissue sections of 154 consecutive invasive breast carcinomas were stained for Ki67 and cytokeratin 8/18 with immunohistochemistry in a clinical setting. Ki67 proliferation index was determined using two independent DIA platforms, implementing VDS to identify tumor tissue. Manual Ki67 score was determined using a standardized manual counting protocol. Inter-observer agreement between manual and DIA scores and inter-platform agreement between both DIA platforms were determined and calculated using Spearman’s correlation coefficients. Correlations and agreement were assessed with scatterplots and Bland–Altman plots.

Results

Spearman’s correlation coefficients were 0.94 (p < 0.001) for inter-observer agreement between manual counting and platform A, 0.93 (p < 0.001) between manual counting and platform B, and 0.96 (p < 0.001) for inter-platform agreement. Scatterplots and Bland–Altman plots revealed no skewness within specific data ranges. In the few cases with ≥ 10% difference between manual counting and DIA, results by both platforms were similar.

Conclusions

DIA using VDS is an accurate method to determine the Ki67 proliferation index in breast cancer, as an alternative to manual scoring of whole sections in clinical practice. Inter-platform agreement between two different DIA platforms was excellent, suggesting vendor-independent clinical implementability.

Electronic supplementary material

The online version of this article (10.1007/s10549-018-4669-2) contains supplementary material, which is available to authorized users.

Large-scale DNA organization is a prognostic marker of breast cancer survival

Breast cancer is the leading cause of cancer-related deaths among women worldwide. We investigated whether changes in large-scale DNA organization (LDO) of tumor epithelial nuclei are an indicator of the aggressiveness of the tumor. We tested our algorithm on a set of 172 duplicates TMA cores samples coming from 95 breast cancer patients. Thirty-five patients died of breast cancer, and 60 were still alive 10 years after surgery. Duplicates cores were used to create training and test set. The TMA slides were stained with Feulgen-thionin and imaged using our in-house high-resolution Imaging system. Automated segmentation of cell nuclei followed by manual selection of intact, in-focus nuclei resulted in an average of 50 cell nuclei per sample available for analysis. Using forward stepwise linear discriminant analysis, a combination of six features that combined linearly gave the best discrimination between the two groups of cells: cells collected from 'deceased' patients TMA specimens and cells collected from "survivors" patients TMA specimens. Five of these features measure the spatial organization of DNA chromatin. The resulting canonical score is named cell LDO score. A patient LDO score, percentage of cell nuclei with a cell LDO score higher than a predefined cutoff value, was processed for the specimens in the test set, and a cutoff value was defined to classify patients with a low or a high LDO score. Using this binary test, 82.1% of patients were correctly classified are "deceased" or "survivors," with a specificity of 79% and a sensitivity of 88%. The relative risk of death of an individual with a high LDO score was nine times higher than for a patient with a low LDO score. When testing the combination of LDO score, node status, histological grade, and tumor grade to predict breast cancer survival, LDO was the most significant predictor. LDO classification was also highly associated with survival for only grade 1 and 2 patients as well as for only grade 3 patients. Our result confirms the potential of LDO to measure phenotypic changes associated with more aggressive disease and could be evaluated to identify patients more likely to benefit from adjuvant therapies.

Characterizing the heterogeneity of tumor tissues from spatially resolved molecular measures

Background

Tumor heterogeneity can manifest itself by sub-populations of cells having distinct phenotypic profiles expressed as diverse molecular, morphological and spatial distributions. This inherent heterogeneity poses challenges in terms of diagnosis, prognosis and efficient treatment. Consequently, tools and techniques are being developed to properly characterize and quantify tumor heterogeneity. Multiplexed immunofluorescence (MxIF) is one such technology that offers molecular insight into both inter-individual and intratumor heterogeneity. It enables the quantification of both the concentration and spatial distribution of 60+ proteins across a tissue section. Upon bioimage processing, protein expression data can be generated for each cell from a tissue field of view.

Results

The Multi-Omics Heterogeneity Analysis (MOHA) tool was developed to compute tissue heterogeneity metrics from MxIF spatially resolved tissue imaging data. This technique computes the molecular state of each cell in a sample based on a pathway or gene set. Spatial states are then computed based on the spatial arrangements of the cells as distinguished by their respective molecular states. MOHA computes tissue heterogeneity metrics from the distributions of these molecular and spatially defined states. A colorectal cancer cohort of approximately 700 subjects with MxIF data is presented to demonstrate the MOHA methodology. Within this dataset, statistically significant correlations were found between the intratumor AKT pathway state diversity and cancer stage and histological tumor grade. Furthermore, intratumor spatial diversity metrics were found to correlate with cancer recurrence.

Conclusions

MOHA provides a simple and robust approach to characterize molecular and spatial heterogeneity of tissues. Research projects that generate spatially resolved tissue imaging data can take full advantage of this useful technique. The MOHA algorithm is implemented as a freely available R script (see supplementary information).

Evaluation of the metastatic potential of malignant cells by image processing of digital holographic microscopy data

The cell refractive index has been proposed as a putative cancer biomarker of great potential, being correlated with cell content and morphology, cell division rate and membrane permeability. We used digital holographic microscopy to compare the refractive index and dry mass density of two B16 murine melanoma sublines of different metastatic potential. Using statistical methods, the distribution of phase shifts within the reconstructed quantitative phase images was analyzed by the method of bimodality coefficients. The observed correlation of refractive index, dry mass density and bimodality profile with the metastatic potential of the cells was validated by real time impedance-based assay and clonogenic tests. We suggest that the refractive index and bimodality analysis of quantitative phase image histograms could be developed as optical biomarkers useful in label-free detection and quantitative evaluation of cell metastatic potential.

Automated Image Analysis of HER2 Fluorescence In Situ Hybridization to Refine Definitions of Genetic Heterogeneity in Breast Cancer Tissue

Human epidermal growth factor receptor 2 gene- (HER2-) targeted therapy for breast cancer relies primarily on HER2 overexpression established by immunohistochemistry (IHC) with borderline cases being further tested for amplification by fluorescence in situ hybridization (FISH). Manual interpretation of HER2 FISH is based on a limited number of cells and rather complex definitions of equivocal, polysomic, and genetically heterogeneous (GH) cases. Image analysis (IA) can extract high-capacity data and potentially improve HER2 testing in borderline cases. We investigated statistically derived indicators of HER2 heterogeneity in HER2 FISH data obtained by automated IA of 50 IHC borderline (2+) cases of invasive ductal breast carcinoma. Overall, IA significantly underestimated the conventional HER2, CEP17 counts, and HER2/CEP17 ratio; however, it collected more amplified cells in some cases below the lower limit of GH definition by manual procedure. Indicators for amplification, polysomy, and bimodality were extracted by factor analysis and allowed clustering of the tumors into amplified, nonamplified, and equivocal/polysomy categories. The bimodality indicator provided independent cell diversity characteristics for all clusters. Tumors classified as bimodal only partially coincided with the conventional GH heterogeneity category. We conclude that automated high-capacity nonselective tumor cell assay can generate evidence-based HER2 intratumor heterogeneity indicators to refine GH definitions.

An adaptive positivity thresholding method for automated Ki67 hotspot detection (AKHoD) in breast cancer biopsies

The proliferative activity of breast cancer tissue can be estimated using the Ki67 biomarker. The percentage of positivity of such biomarker is correlated with proliferation and consequently with the prognosis of a breast tumor. Ki67 marked tissue samples are analyzed by an experienced pathologist who identifies the most active areas of tumor cell proliferation called hotspots, and estimates the positivity of each case. A method for the Automated Ki67 Hotspot Detection (AKHoD) is presented in this work. The main objective of the AKHoD method is to automatically and efficiently provide the pathologist with suggestions about Ki67 hotspot areas as a decision support. The input of AKHoD is a digital slide that is divided in tiles. For each tile, AKHoD provides a rough estimate of positivity and cellularity, summarized in very low resolution positivity and cellularity images. In a second step, an adaptive thresholding is applied to such positivity image to identify the most positive connected and convex areas, within cellularity limits set by current guidelines (that is, 500-2000). The method has been preliminarily validated on 50 digital slides for which three expert pathologists provided gold standard hotspots. 82% of the gold standard hotspots have been successfully recognized by the system, spending an average of 54s per slide. While further validation is needed taking into account also patients follow-up, this first experimentation suggests that the proposed method could be adequate for supporting the pathologist in hotspot detection.