Quantification of spatial tumor heterogeneity in immunohistochemistry staining images.
Ontology highlight
ABSTRACT: Motivation
Quantitative immunofluorescence is often used for immunohistochemistry quantification of proteins that serve as cancer biomarkers. Advanced image analysis systems for pathology allow capturing expression levels in each individual cell or subcellular compartment. However, only the mean signal intensity within the cancer tissue region of interest is usually considered as biomarker completely ignoring the issue of tumor heterogeneity.Results
We propose using immunohistochemistry image-derived information on the spatial distribution of cellular signal intensity (CSI) of protein expression within the cancer cell population to quantify both mean expression level and tumor heterogeneity of CSI levels. We view CSI levels as marks in a marked point process of cancer cells in the tissue and define spatial indices based on conditional mean and conditional variance of the marked point process. The proposed methodology provides objective metrics of cell-to-cell heterogeneity in protein expressions that allow discriminating between different patterns of heterogeneity. The prognostic utility of new spatial indices is investigated and compared to the standard mean signal intensity biomarkers using the protein expressions in tissue microarrays incorporating tumor tissues from 1000+ breast cancer patients.Availability and implementation: the r code for computing the proposed spatial indices is included as supplementary material
.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Chervoneva I
PROVIDER: S-EPMC8208754 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The quantification of the expression of different molecules is a key question in both basic and applied sciences. While protein quantification through molecular techniques leads to the loss of spatial information and resolution, immunohistochemistry is usually associated with time-consuming image analysis and human bias. In addition, the scarce automatic software analysis is often proprietary and expensive and relies on a fixed threshold binarization. Here we describe and share a set of macros ready for automated fluorescence analysis of large batches of fixed tissue samples using FIJI/ImageJ. The quantification of the molecules of interest are based on an automatic threshold analysis of immunofluorescence images to automatically identify the top brightest structures of each image. These macros measure several parameters commonly quantified in basic neuroscience research, such as neuropil density and fluorescence intensity of synaptic puncta, perisomatic innervation and col-localization of different molecules and analysis of the neurochemical phenotype of neuronal subpopulations. In addition, these same macro functions can be easily modified to improve similar analysis of fluorescent probes in human biopsies for diagnostic purposes based on the expression patterns of several molecules.
| S-EPMC6039854 | biostudies-other
Project description:With immunohistochemical (IHC) staining increasingly being used to guide clinical decisions, variability in staining quality and reproducibility are becoming essential factors in generating diagnoses using IHC tissue preparations. The current study tested a method to track and quantify the interrun, intrarun, and intersite variability of IHC staining intensity. Our hypothesis was that staining precision between laboratory sites, staining runs, and individual slides may be verified quantitatively, efficiently and effectively utilizing algorithm-based, automated image analysis. To investigate this premise, we tested the consistency of IHC staining in 40 routinely processed (formalin-fixed, paraffin-embedded) human tissues using 10 common antibiomarker antibodies on 2 Dako Omnis instruments at 2 locations (Carpinteria, CA: 30 m above sea level and Longmont, CO: 1500 m above sea level) programmed with identical, default settings and sample pretreatments. Digital images of IHC-labeled sections produced by a whole slide scanner were analyzed by a simple commercially available algorithm and compared with a board-certified veterinary pathologist's semiquantitative scoring of staining intensity. The image analysis output correlated well with pathology scores but had increased sensitivity for discriminating subtle variations and providing reproducible digital quantification across sites as well as within and among staining runs at the same site. Taken together, our data indicate that digital image analysis offers an objective and quantifiable means of verifying IHC staining parameters as a part of laboratory quality assurance systems.
| S-EPMC7368846 | biostudies-literature
Project description:Tumor-associated macrophages (TAMs), one of the most abundant immune components in gastric cancer (GC), are difficult to characterize due to their heterogeneity. Multiple approaches have been used to elucidate the issue, however, due to the tissue-destructive nature of most of these methods, the spatial distribution of TAMs in situ remains unclear. Here we probe the relationship between tumor context and TAM heterogeneity by multiplex immunohistochemistry of 56 human GC cases. Using distinct expression marker profiles on TAMs, we report seven predominant populations distributed between tumor and non-tumor tissue. TAM population-associated gene signatures reflect their heterogeneity and polarization in situ. Increased density of CD163+ (CD206-) TAMs with concurrent high CD68 expression is associated with upregulated immune-signaling and improved patient survival by univariate, but not multivariate analysis. CD68-only and CD206+ TAMs are correlated with high PDL1 expression.
| S-EPMC6718690 | biostudies-literature
Project description:Cancer marker discovery is an emerging topic in high-throughput quantitative proteomics. However, the omics technology usually generates a long list of marker candidates that requires a labor-intensive filtering process in order to screen for potentially useful markers. Specifically, various parameters, such as the level of overexpression of the marker in the cancer type of interest, which is related to sensitivity, and the specificity of the marker among cancer groups, are the most critical considerations. Protein expression profiling on the basis of immunohistochemistry (IHC) staining images is a technique commonly used during such filtering procedures. To systematically investigate the protein expression in different cancer versus normal tissues and cell types, the Human Protein Atlas is a most comprehensive resource because it includes millions of high-resolution IHC images with expert-curated annotations. To facilitate the filtering of potential biomarker candidates from large-scale omics datasets, in this study we have proposed a scoring approach for quantifying IHC annotation of paired cancerous/normal tissues and cancerous/normal cell types. We have comprehensively calculated the scores of all the 17219 tested antibodies deposited in the Human Protein Atlas based on their accumulated IHC images and obtained 457110 scores covering 20 different types of cancers. Statistical tests demonstrate the ability of the proposed scoring approach to prioritize cancer-specific proteins. Top 100 potential marker candidates were prioritized for the 20 cancer types with statistical significance. In addition, a model study was carried out of 1482 membrane proteins identified from a quantitative comparison of paired cancerous and adjacent normal tissues from patients with colorectal cancer (CRC). The proposed scoring approach demonstrated successful prioritization and identified four CRC markers, including two of the most widely used, namely CEACAM5 and CEACAM6. These results demonstrate the potential of this scoring approach in terms of cancer marker discovery and development. All the calculated scores are available at http://bal.ym.edu.tw/hpa/.
| S-EPMC3841220 | biostudies-literature
Project description: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.
| S-EPMC7606448 | biostudies-literature
Project description:The spatial organization of different types of cells in tumor tissues reveals important information about the tumor microenvironment (TME). To facilitate the study of cellular spatial organization and interactions, we developed Histology-based Digital-Staining, a deep learning-based computation model, to segment the nuclei of tumor, stroma, lymphocyte, macrophage, karyorrhexis, and red blood cells from standard hematoxylin and eosin-stained pathology images in lung adenocarcinoma. Using this tool, we identified and classified cell nuclei and extracted 48 cell spatial organization-related features that characterize the TME. Using these features, we developed a prognostic model from the National Lung Screening Trial dataset, and independently validated the model in The Cancer Genome Atlas lung adenocarcinoma dataset, in which the predicted high-risk group showed significantly worse survival than the low-risk group (P = 0.001), with a HR of 2.23 (1.37-3.65) after adjusting for clinical variables. Furthermore, the image-derived TME features significantly correlated with the gene expression of biological pathways. For example, transcriptional activation of both the T-cell receptor and programmed cell death protein 1 pathways positively correlated with the density of detected lymphocytes in tumor tissues, while expression of the extracellular matrix organization pathway positively correlated with the density of stromal cells. In summary, we demonstrate that the spatial organization of different cell types is predictive of patient survival and associated with the gene expression of biological pathways. SIGNIFICANCE: These findings present a deep learning-based analysis tool to study the TME in pathology images and demonstrate that the cell spatial organization is predictive of patient survival and is associated with gene expression.See related commentary by Rodriguez-Antolin, p. 1912.
| S-EPMC7919065 | biostudies-literature
Project description:BackgroundTumour heterogeneity is considered an important mechanism of treatment failure. Imaging-based assessment of tumour heterogeneity is showing promise but the relationship between these mathematically derived measures and accepted 'gold standards' of tumour biology such as immunohistochemical measures is not established.MethodsA total of 20 women with primary breast cancer underwent a research dynamic contrast-enhanced computed tomography prior to treatment with data being available for 15 of these. Texture analysis was performed of the primary tumours to extract 13 locoregional and global parameters. Immunohistochemical analysis associations were assessed by the Spearman rank correlation.ResultsHypoxia-inducible factor-1α was correlated with first-order kurtosis (r = -0.533, P = .041) and higher order neighbourhood grey-tone difference matrix coarseness (r = 0.54, P = .038). Vascular maturity-related smooth muscle actin was correlated with higher order grey-level run-length long-run emphasis (r = -0.52, P = .047), fractal dimension (r = 0.613, P = .015), and lacunarity (r = -0.634, P = .011). Micro-vessel density, reflecting angiogenesis, was also associated with lacunarity (r = 0.547, P = .035).ConclusionsThe associations suggest a biological basis for these image-based heterogeneity features and support the use of imaging, already part of standard care, for assessing intratumoural heterogeneity.
| S-EPMC6552350 | biostudies-literature
Project description:ImportanceLissamine green (LG) staining of the conjunctiva is a key biomarker in evaluating ocular surface disease. The disease currently is assessed using relatively coarse subjective scales. Objective assessment would standardize comparisons over time and between clinicians.ObjectiveTo develop a semiautomated, quantitative system to assess lissamine green staining of the bulbar conjunctiva on digital images.Design, setting, and participantsUsing a standard photography protocol, 35 digital images of the conjunctiva of 11 patients with a diagnosis of dry eye disease based on characteristic signs and symptoms were obtained after topical administration of preservative-free LG, 1%, solution. Images were scored independently by 2 masked ophthalmologists in an academic medical center using the van Bijsterveld and National Eye Institute (NEI) scales. The region of interest was identified by manually marking 7 anatomic landmarks on the images. An objective measure was developed by segmenting the images, forming a vector of key attributes, and then performing a random forest regression. Subjective scores were correlated with the output from a computer algorithm using a cross-validation technique. The ranking of images from least to most staining was compared between the algorithm and the ophthalmologists. The study was conducted from April 26, 2012, through June 2, 2016.Main outcomes and measuresCorrelation and level of agreement among computerized algorithm scores, van Bijsterveld scale clinical scores, and NEI scale clinical scores.ResultsThe scores from the automated algorithm correlated well with the mean scores obtained from the gradings of 2 ophthalmologists for the 35 images using the van Bijsterveld scale (Spearman correlation coefficient, rs = 0.79), and moderately with the NEI scale (rs = 0.61) scores. For qualitative ranking of staining, the correlation between the automated algorithm and the 2 ophthalmologists was rs = 0.78 and rs = 0.83.Conclusions and relevanceThe algorithm performed well when evaluating LG staining of the conjunctiva, as evidenced by good correlation with subjective gradings using 2 different grading scales. Future longitudinal studies are needed to assess the responsiveness of the algorithm to change of conjunctival staining over time.
| S-EPMC5710489 | biostudies-literature
Project description:BackgroundCurrent intra-tumoral heterogeneous feature extraction in radiology is limited to the use of a single slice or the region of interest within a few context-associated slices, and the decoding of intra-tumoral spatial heterogeneity using whole tumor samples is rare. We aim to propose a mathematical model of space-filling curve-based spatial correspondence mapping to interpret intra-tumoral spatial locality and heterogeneity.MethodsA Hilbert curve-based approach was employed to decode and visualize intra-tumoral spatial heterogeneity by expanding the tumor volume to a two-dimensional (2D) matrix in voxels while preserving the spatial locality of the neighboring voxels. The proposed method was validated using three-dimensional (3D) volumes constructed from lung nodules from the LIDC-IDRI dataset, regular axial plane images, and 3D blocks.ResultsDimensionality reduction of the Hilbert volume with a single regular axial plane image showed a sparse and scattered pixel distribution on the corresponding 2D matrix. However, for 3D blocks and lung tumor inside the volume, the dimensionality reduction to the 2D matrix indicated regular and concentrated squares and rectangles. For classification into benign and malignant masses using lung nodules from the LIDC-IDRI dataset, the Inception-V4 indicated that the Hilbert matrix images improved accuracy (85.54% vs. 73.22%, p < 0.001) compared to the original CT images of the test dataset.ConclusionsOur study indicates that Hilbert curve-based spatial correspondence mapping is promising for decoding intra-tumoral spatial heterogeneity of partial or whole tumor samples on radiological images. This spatial-locality-preserving approach for voxel expansion enables existing radiomics and convolution neural networks to filter structured and spatially correlated high-dimensional intra-tumoral heterogeneity.
| S-EPMC8557226 | biostudies-literature
Project description:A new workflow for protein-based tumor heterogeneity probing in tissues is here presented. Tumor heterogeneity is believed to be key for therapy failure and differences in prognosis in cancer patients. Comprehending tumor heterogeneity, especially at the protein level, is critical for tracking tumor evolution, and showing the presence of different phenotypical variants and their location with respect to tissue architecture. Although a variety of techniques is available for quantifying protein expression, the heterogeneity observed in the tissue is rarely addressed. The proposed method is validated in breast cancer fresh-frozen tissues derived from five patients. Protein expression is quantified on the tissue regions of interest (ROI) with a resolution of up to 100 μm in diameter. High heterogeneity values across the analyzed patients in proteins such as cytokeratin 7, β-actin and epidermal growth factor receptor (EGFR) using a Shannon entropy analysis are observed. Additionally, ROIs are clustered according to their expression levels, showing their location in the tissue section, and highlighting that similar phenotypical variants are not always located in neighboring regions. Interestingly, a patient with a phenotype related to increased aggressiveness of the tumor presents a unique protein expression pattern. In summary, a workflow for the localized extraction and protein analysis of regions of interest from frozen tissues, enabling the evaluation of tumor heterogeneity at the protein level is presented.
| S-EPMC8604290 | biostudies-literature