Digital image analysis (DIA) is becoming central to the quantitative evaluation of tissue biomarkers for discovery, diagnosis and therapeutic selection for the delivery of precision medicine. In this study, automated DIA using a new purpose-built software platform (QuPath) is applied to a cohort of 293 breast cancer patients to score five biomarkers in tissue microarrays (TMAs): ER, PR, HER2, Ki67 and p53. This software is able to measure IHC expression following fully automated tumor recognition in the same immunohistochemical (IHC)-stained tissue section, as part of a rapid workflow to ensure objectivity and accelerate biomarker analysis. The digital scores produced by QuPath were compared with manual scores by a pathologist and shown to have a good level of concordance in all cases (Cohen's κ40.6), and almost perfect agreement for the clinically relevant biomarkers ER, PR and HER2 (κ40.86). To assess prognostic value, cutoff thresholds could be applied to both manual and automated scores using the QuPath software, and survival analysis performed for 5-year overall survival. DIA was shown to be capable of replicating the statistically significant stratification of patients achieved using manual scoring across all biomarkers (Po0.01, log-rank test). Furthermore, the image analysis scores were shown to consistently lead to statistical significance across a wide range of potential cutoff thresholds, indicating the robustness of the method, and identify sub-populations of cases exhibiting different expression patterns within the p53 and Ki67 data sets that warrant further investigation. These findings have demonstrated QuPath's suitability for fast, reproducible, high-throughput TMA analysis across a range of important biomarkers. This was achieved using our tumor recognition algorithms for IHC-stained sections, trained interactively without the need for any additional tumor recognition markers, for example, cytokeratin, to obtain greater insight into the relationship between biomarker expression and clinical outcome applicable to a range of cancer types.

Patients with extreme p53 IHC expression have a worse OS than those with NE expression. Accounting for EN as well as EP expression improves the prognostic impact. Extreme expression positively correlates with nodal stage and histological grade, and negatively with hormone receptor status. Extreme expression may relate to specific mutational status.

Given the morphological overlap and significant association with high-grade CIN and AIS, we believe that SMILE is best regarded as a form of high-grade reserve cell dysplasia, probably arising from reserve or stem cells with the capacity for multidirectional differentiation. We recommend that management of SMILE should be as for AIS.

Breast cancer is a heterogeneous disease, at both an inter- and intra-tumoural level. Appreciating heterogeneity through the application of biomarkers and molecular signatures adds complexity to tumour taxonomy but is key to personalising diagnosis, treatment and prognosis. The extent to which heterogeneity exists, and its interpretation remains a challenge to pathologists. Using HER2 as an exemplar, we have developed a simple reproducible heterogeneity index. Cell-to-cell HER2 heterogeneity was extensive in a proportion of both reported ‘amplified’ and ‘non-amplified’ cases. The highest levels of heterogeneity objectively identified occurred in borderline categories and higher ratio non-amplified cases. A case with particularly striking heterogeneity was analysed further with an array of biomarkers in order to assign a molecular diagnosis. Broad biological complexity was evident. In essence, interpretation, depending on the area of tumour sampled, could have been one of three distinct phenotypes, each of which would infer different therapeutic interventions. Therefore, we recommend that heterogeneity is assessed and taken into account when determining treatment options.