Computational pathology model to assess acute and chronic transformations of the tubulointerstitial compartment in renal allograft biopsies.

Managing patients with kidney allografts largely depends on biopsy diagnosis which is based on semiquantitative assessments of rejection features and extent of acute and chronic changes within the renal parenchyma. Current methods lack reproducibility while digital image data-driven computational models enable comprehensive and quantitative assays. In this study we aimed to develop a computational method for automated assessment of histopathology transformations within the tubulointerstitial compartment of the renal cortex. Whole slide images of modified Picrosirius red-stained biopsy slides were used for the training (n = 852) and both internal (n = 172) and external (n = 94) tests datasets. The pipeline utilizes deep learning segmentations of renal tubules, interstitium, and peritubular capillaries from which morphometry features were extracted. Seven indicators were selected for exploring the intrinsic spatial interactions within the tubulointerstitial compartment. A principal component analysis revealed two independent factors which can be interpreted as representing chronic and acute tubulointerstitial injury. A K-means clustering classified biopsies according to potential phenotypes of combined acute and chronic transformations of various degrees. We conclude that multivariate analyses of tubulointerstitial morphometry transformations enable extraction of and quantification of acute and chronic components of injury. The method is developed for renal allograft biopsies; however, the principle can be applied more broadly for kidney pathology assessment.

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.

Improving HCC Prognostic Models after Liver Resection by AI-Extracted Tissue Fiber Framework Analytics.

Despite advances in diagnostic and treatment technologies, predicting outcomes of patients with hepatocellular carcinoma (HCC) remains a challenge. Prognostic models are further obscured by the variable impact of the tumor properties and the remaining liver parenchyma, often affected by cirrhosis or non-alcoholic fatty liver disease that tend to precede HCC. This study investigated the prognostic value of reticulin and collagen microarchitecture in liver resection samples. We analyzed 105 scanned tissue sections that were stained using a Gordon and Sweet's silver impregnation protocol combined with Picric Acid-Sirius Red. A convolutional neural network was utilized to segment the red-staining collagen and black linear reticulin strands, generating a detailed map of the fiber structure within the HCC and adjacent liver tissue. Subsequent hexagonal grid subsampling coupled with automated epithelial edge detection and computational fiber morphometry provided the foundation for region-specific tissue analysis. Two penalized Cox regression models using LASSO achieved a concordance index (C-index) greater than 0.7. These models incorporated variables such as patient age, tumor multifocality, and fiber-derived features from the epithelial edge in both the tumor and liver compartments. The prognostic value at the tumor edge was derived from the reticulin structure, while collagen characteristics were significant at the epithelial edge of peritumoral liver. The prognostic performance of these models was superior to models solely reliant on conventional clinicopathologic parameters, highlighting the utility of AI-extracted microarchitectural features for the management of HCC.

Systematically higher Ki67 scores on core biopsy samples compared to corresponding resection specimen in breast cancer: a multi-operator and multi-institutional study.

Ki67 has potential clinical importance in breast cancer but has yet to see broad acceptance due to inter-laboratory variability. Here we tested an open source and calibrated automated digital image analysis (DIA) platform to: (i) investigate the comparability of Ki67 measurement across corresponding core biopsy and resection specimen cases, and (ii) assess section to section differences in Ki67 scoring. Two sets of 60 previously stained slides containing 30 core-cut biopsy and 30 corresponding resection specimens from 30 estrogen receptor-positive breast cancer patients were sent to 17 participating labs for automated assessment of average Ki67 expression. The blocks were centrally cut and immunohistochemically (IHC) stained for Ki67 (MIB-1 antibody). The QuPath platform was used to evaluate tumoral Ki67 expression. Calibration of the DIA method was performed as in published studies. A guideline for building an automated Ki67 scoring algorithm was sent to participating labs. Very high correlation and no systematic error (p = 0.08) was found between consecutive Ki67 IHC sections. Ki67 scores were higher for core biopsy slides compared to paired whole sections from resections (p ≤ 0.001; median difference: 5.31%). The systematic discrepancy between core biopsy and corresponding whole sections was likely due to pre-analytical factors (tissue handling, fixation). Therefore, Ki67 IHC should be tested on core biopsy samples to best reflect the biological status of the tumor.

Machine learning-enabled cancer diagnostics with widefield polarimetric second-harmonic generation microscopy.

The extracellular matrix (ECM) collagen undergoes major remodeling during tumorigenesis. However, alterations to the ECM are not widely considered in cancer diagnostics, due to mostly uniform appearance of collagen fibers in white light images of hematoxylin and eosin-stained (H&E) tissue sections. Polarimetric second-harmonic generation (P-SHG) microscopy enables label-free visualization and ultrastructural investigation of non-centrosymmetric molecules, which, when combined with texture analysis, provides multiparameter characterization of tissue collagen. This paper demonstrates whole slide imaging of breast tissue microarrays using high-throughput widefield P-SHG microscopy. The resulting P-SHG parameters are used in classification to differentiate tumor from normal tissue, resulting in 94.2% for both accuracy and F1-score, and 6.3% false discovery rate. Subsequently, the trained classifier is employed to predict tumor tissue with 91.3% accuracy, 90.7% F1-score, and 13.8% false omission rate. As such, we show that widefield P-SHG microscopy reveals collagen ultrastructure over large tissue regions and can be utilized as a sensitive biomarker for cancer diagnostics and prognostics studies.

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.

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.

Immunogradient Indicators for Antitumor Response Assessment by Automated Tumor-Stroma Interface Zone Detection.

The distribution of tumor-infiltrating lymphocytes (TILs) within the tumor microenvironment provides strong prognostic value, which is increasingly important with the arrival of new immunotherapy modalities. Both visual and image analysis-based assays are developed to assess the immune contexture of the tumors. We propose an automated method based on grid subsampling of microscopy image analysis data to extract the tumor-stroma interface zone (IZ) of controlled width. The IZ is a ranking of tissue areas by their distance to the tumor edge, which is determined by a set of explicit rules. TIL density profiles across the IZ are used to compute a set of novel immunogradient indicators that reflect TIL gradient towards the tumor. We applied this method on CD8 immunohistochemistry images of surgically excised hormone receptor-positive breast and colorectal cancers to predict overall patient survival. In both cohorts, the immunogradient indicators enabled strong and independent prognostic stratification, outperforming clinical and pathologic variables. Patients with breast cancer with low immunogradient levels had a prominent decrease in survival probability 5 years after surgery. Our study provides proof of concept that data-driven, automated, operator-independent IZ sampling enables spatial immune response measurement in the tumor-host interaction frontline for prediction of disease outcomes.

Proliferation Tumour Marker Network (PTM-NET) for the identification of tumour region in Ki67 stained breast cancer whole slide images.

Uncontrolled proliferation is a hallmark of cancer and can be assessed by labelling breast tissue using immunohistochemistry for Ki67, a protein associated with cell proliferation. Accurate measurement of Ki67-positive tumour nuclei is of critical importance, but requires annotation of the tumour regions by a pathologist. This manual annotation process is highly subjective, time-consuming and subject to inter- and intra-annotator experience. To address this challenge, we have developed Proliferation Tumour Marker Network (PTM-NET), a deep learning model that objectively annotates the tumour regions in Ki67-labelled breast cancer digital pathology images using a convolution neural network. Our custom designed deep learning model was trained on 45 immunohistochemical Ki67-labelled whole slide images to classify tumour and non-tumour regions and was validated on 45 whole slide images from two different sources that were stained using different protocols. Our results show a Dice coefficient of 0.74, positive predictive value of 70% and negative predictive value of 88.3% against the manual ground truth annotation for the combined dataset. There were minimal differences between the images from different sources and the model was further tested in oestrogen receptor and progesterone receptor-labelled images. Finally, using an extension of the model, we could identify possible hotspot regions of high proliferation within the tumour. In the future, this approach could be useful in identifying tumour regions in biopsy samples and tissue microarray images.

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.

The Role of Serum Adiponectin for Outcome Prediction in Patients with Dilated Cardiomyopathy and Advanced Heart Failure.

Clinical interpretation of patients' plasma adiponectin (APN) remains challenging; its value as biomarker in dilated cardiomyopathy (DCM) is equivocal. We evaluated whether circulating APN level is an independent predictor of composite outcome: death, left ventricle assist device (LVAD) implantation, and heart transplantation (HT) in patients with nonischemic DCM. 57 patients with nonischemic DCM (average LV diastolic diameter 6.85 cm, LV ejection fraction 26.63%, and pulmonary capillary wedge pressure 22.06 mmHg) were enrolled. Patients underwent echocardiography, right heart catheterization, and endomyocardial biopsy. During a mean follow-up of 33.42 months, 15 (26%) patients died, 12 (21%) patients underwent HT, and 8 (14%) patients were implanted with LVAD. APN level was significantly higher in patients who experienced study endpoints (23.4 versus 10.9 ug/ml, p = 0.01). APN was associated with worse outcome in univariate Cox proportional hazards model (HR 1.04, CI 1.02-1.07, p = 0.001) but lost significance adjusting for other covariates. Average global strain (AGS) is an independent outcome predictor (HR 1.42, CI 1.081-1.866, p = 0.012). Increased circulating APN level was associated with higher mortality and may be an additive prognostic marker in DCM with advanced HF. Combination of serum (APN, BNP, TNF-α) and echocardiographic (AGS) markers may increase the HF predicting power for the nonischemic DCM patients.

Impact of tissue sampling on accuracy of Ki67 immunohistochemistry evaluation in breast cancer.

BACKGROUND: Gene expression studies have identified molecular subtypes of breast cancer with implications to chemotherapy recommendations. For distinction of these types, a combination of immunohistochemistry (IHC) markers, including proliferative activity of tumor cells, estimated by Ki67 labeling index is used. Clinical studies are frequently based on IHC performed on tissue microarrays (TMA) with variable tissue sampling. This raises the need for evidence-based sampling criteria for individual IHC biomarker studies. We present a novel tissue sampling simulation model and demonstrate its application on Ki67 assessment in breast cancer tissue taking intratumoral heterogeneity into account. METHODS: Whole slide images (WSI) of 297 breast cancer sections, immunohistochemically stained for Ki67, were subjected to digital image analysis (DIA). Percentage of tumor cells stained for Ki67 was computed for hexagonal tiles super-imposed on the WSI. From this, intratumoral Ki67 heterogeneity indicators (Haralick's entropy values) were extracted and used to dichotomize the tumors into homogeneous and heterogeneous subsets. Simulations with random selection of hexagons, equivalent to 0.75 mm circular diameter TMA cores, were performed. The tissue sampling requirements were investigated in relation to tumor heterogeneity using linear regression and extended error analysis. RESULTS: The sampling requirements were dependent on the heterogeneity of the biomarker expression. To achieve a coefficient error of 10 %, 5-6 cores were needed for homogeneous cases, 11-12 cores for heterogeneous cases; in mixed tumor population 8 TMA cores were required. Similarly, to achieve the same accuracy, approximately 4,000 nuclei must be counted when the intratumor heterogeneity is mixed/unknown. Tumors of low proliferative activity would require larger sampling (10-12 TMA cores, or 6,250 nuclei) to achieve the same error measurement results as for highly proliferative tumors. CONCLUSIONS: Our data show that optimal tissue sampling for IHC biomarker evaluation is dependent on the heterogeneity of the tissue under study and needs to be determined on a per use basis. We propose a method that can be applied to determine the sampling strategy for specific biomarkers, tissues and study targets. In addition, our findings highlight the benefit of high-capacity computer-based IHC measurement techniques to improve accuracy of the testing.

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

Proliferative activity, assessed by Ki67 immunohistochemistry (IHC), is an established prognostic and predictive biomarker of breast cancer (BC). However, it remains under-utilized due to lack of standardized robust measurement methodologies and significant intratumor heterogeneity of expression. A recently proposed methodology for IHC biomarker assessment in whole slide images (WSI), based on systematic subsampling of tissue information extracted by digital image analysis (DIA) into hexagonal tiling arrays, enables computation of a comprehensive set of Ki67 indicators, including intratumor variability. In this study, the tiling methodology was applied to assess Ki67 expression in WSI of 152 surgically removed Ki67-stained (on full-face sections) BC specimens and to test which, if any, Ki67 indicators can predict overall survival (OS). Visual Ki67 IHC estimates and conventional clinico-pathologic parameters were also included in the study. Analysis revealed linearly independent intrinsic factors of the Ki67 IHC variance: proliferation (level of expression), disordered texture (entropy), tumor size and Nottingham Prognostic Index, bimodality, and correlation. All visual and DIA-generated indicators of the level of Ki67 expression provided significant cutoff values as single predictors of OS. However, only bimodality indicators (Ashman's D, in particular) were independent predictors of OS in the context of hormone receptor and HER2 status. From this, we conclude that spatial heterogeneity of proliferative tumor activity, measured by DIA of Ki67 IHC expression and analyzed by the hexagonal tiling approach, can serve as an independent prognostic indicator of OS in BC patients that outperforms the prognostic power of the level of proliferative activity.