Spatial correlation between in vivo imaging and immunohistochemical biomarkers: A methodological study.

In this study, we present a method that enables voxel-by-voxel comparison of in vivo imaging to immunohistochemistry (IHC) biomarkers. As a proof of concept, we investigated the spatial correlation between dynamic contrast enhanced (DCE-)CT parameters and IHC biomarkers Ki-67 (proliferation), HIF-1α (hypoxia), and CD45 (immune cells). 54 whole-mount tumor slices of 15 laryngeal and hypopharyngeal carcinomas were immunohistochemically stained and digitized. Heatmaps of biomarker positivity were created and registered to DCE-CT parameter maps. The adiabatic approximation to the tissue homogeneity model was used to fit the following DCE parameters: Ktrans (transfer constant), Ve (extravascular and extracellular space), and Vi (intravascular space). Both IHC and DCE maps were downsampled to 4 × 4 × 3 mm[3] voxels. The mean values per tumor were used to calculate the between-subject correlations between parameters. For the within-subject (spatial) correlation, values of all voxels within a tumor were compared using the repeated measures correlation (rrm). No between-subject correlations were found between IHC biomarkers and DCE parameters, whereas we found multiple significant within-subject correlations: Ve and Ki-67 (rrm = -0.17, P < .001), Ve and HIF-1α (rrm = -0.12, P < .001), Ktrans and CD45 (rrm = 0.13, P < .001), Vi and CD45 (rrm = 0.16, P < .001), and Vi and Ki-67 (rrm = 0.08, P = .003). The strongest correlation was found between IHC biomarkers Ki-67 and HIF-1α (rrm = 0.35, P < .001). This study shows the technical feasibility of determining the 3 dimensional spatial correlation between histopathological biomarker heatmaps and in vivo imaging. It also shows that between-subject correlations do not reflect within-subject correlations of parameters.

Improved delineation with diffusion weighted imaging for laryngeal and hypopharyngeal tumors validated with pathology.

OBJECTIVE: This study aims to determine the added value of a geometrically accurate diffusion-weighted (DW-) MRI sequence on the accuracy of gross tumor volume (GTV) delineations, using pathological tumor delineations as a ground truth. METHODS: Sixteen patients with laryngeal or hypopharyngeal carcinoma were included. After total laryngectomy, the specimen was cut into slices. Photographs of these slices were stacked to create a 3D digital specimen reconstruction, which was registered to the in vivo imaging. The pathological tumor (tumorHE) was delineated on the specimen reconstruction. Six observers delineated all tumors twice: once with only anatomical MR imaging, and once (a few weeks later) when DW sequences were also provided. The majority voting delineation of session one (GTVMRI) and session two (GTVDW-MRI), as well as the clinical target volumes (CTVs), were compared to the tumorHE. RESULTS: The mean tumorHE volume was 11.1 cm3, compared to a mean GTVMRI volume of 18.5 cm3 and a mean GTVDW-MRI volume of 15.7 cm3. The median sensitivity (tumor coverage) was comparable between sessions: 0.93 (range: 0.61-0.99) for the GTVMRI and 0.91 (range: 0.53-1.00) for the GTVDW-MRI. The CTV volume also decreased when DWI was available, with a mean CTVMR of 47.1 cm3 and a mean CTVDW-MRI of 41.4 cm3. Complete tumor coverage was achieved in 15 and 14 tumors, respectively. CONCLUSION: GTV delineations based on anatomical MR imaging tend to overestimate the tumor volume. The availability of the geometrically accurate DW sequence reduces the GTV overestimation and thereby CTV volumes, while maintaining acceptable tumor coverage.

Using Intratumor Heterogeneity of Immunohistochemistry Biomarkers to Classify Laryngeal and Hypopharyngeal Tumors Based on Histologic Features.

Haralick texture features are used to quantify the spatial distribution of signal intensities within an image. In this study, the heterogeneity of proliferation (Ki-67 expression) and immune cells (CD45 expression) within tumors was quantified and used to classify histologic characteristics of larynx and hypopharynx carcinomas. Of 21 laryngectomy specimens, 74 whole-mount tumor slides were scored on histologic characteristics. Ki-67 and CD45 immunohistochemistry was performed, and all sections were digitized. The tumor area was annotated in QuPath. Haralick features independent of the diaminobenzidine intensity were extracted from the isolated diaminobenzidine signal to quantify intratumor heterogeneity. Haralick features from both Ki-67 and CD45 were used as input for a principal component analysis. A linear support vector machine was fitted to the first 4 principal components for classification and validated with a leave-one-patient-out cross-validation method. Significant differences in individual Haralick features were found between cohesive and noncohesive tumors for CD45 (angular second motion: P =.03, inverse difference moment: P =.009, and entropy: P =.02) and between the larynx and hypopharynx tumors for both CD45 (angular second motion: P =.03, inverse difference moment: P =.007, and entropy: P =.005) and Ki-67 (correlation: P =.003). Therefore, these features were used for classification. The linear classifier resulted in a classification accuracy of 85% for site of origin and 81% for growth pattern. A leave-one-patient-out cross-validation resulted in an error rate of 0.27 and 0.35 for both classifiers, respectively. In conclusion, we show a method to quantify intratumor heterogeneity of immunohistochemistry biomarkers using Haralick features. This study also shows the feasibility of using these features to classify tumors by histologic characteristics. The classifiers created in this study are a proof of concept because more data are needed to create robust classifiers, but the method shows potential for automated tumor classification.

Oncological outcome of vocal cord-only radiotherapy for cT1-T2 glottic laryngeal squamous cell carcinoma.

PURPOSE: Early-stage glottic cancer can be treated with radiotherapy only. Modern radiotherapy solutions allow for individualized dose distributions, hypofractionation and sparing of organs at risk. The target volume used to be the entire voice box. This series describe the oncological outcome and toxicity of individualized vocal cord-only hypofractionated radiotherapy for early stage (cT1a-T2 N0). METHODS: Retrospective cohort study with patients treated in a single center between 2014 and 2020. RESULTS: A total of 93 patients were included. Local control rate was 100% for cT1a, 97% for cT1b and 77% for cT2. Risk factor for local recurrence was smoking during radiotherapy. Laryngectomy-free survival was 90% at 5 years. Grade III or higher late toxicity was 3.7%. CONCLUSION: Vocal cord-only hypofractionated radiotherapy appears to be oncologically safe in early-stage glottic cancer. Modern, image-guided radiotherapy led to comparable results as historical series with very limited late toxicity.

Validation of automated positive cell and region detection of immunohistochemically stained laryngeal tumor tissue using digital image analysis.

Objectives: This study aimed to validate a digital image analysis (DIA) workflow for automatic positive cell detection and positive region delineation for immunohistochemical hypoxia markers with a nuclear (hypoxia-inducible factor 1α [HIF-1α]) and a cytoplasmic (pimonidazole [PIMO]) staining pattern. Materials and methods: 101 tissue fragments from 44 laryngeal tumor biopsies were immunohistochemically stained for HIF-1α and PIMO. QuPath was used to determine the percentage of positive cells and to delineate positive regions automatically. For HIF-1α, only cells with strong staining were considered positive. Three dedicated head and neck pathologists scored the percentage of positive cells using three categories (0: <1%; 1: 1%-33%; 2: >33%;). The pathologists also delineated the positive regions on 14 corresponding PIMO and HIF-1α-stained fragments. The consensus between observers was used as the reference standard and was compared to the automatic delineation. Results: Agreement between categorical positivity scores was 76.2% and 65.4% for PIMO and HIF-1α, respectively. In all cases of disagreement in HIF-1α fragments, the DIA underestimated the percentage of positive cells. As for the region detection, the DIA correctly detected most positive regions on PIMO fragments (false positive area=3.1%, false negative area=0.7%). In HIF-1α, the DIA missed some positive regions (false positive area=1.3%, false negative area=9.7%). Conclusions: Positive cell and region detection on biopsy material is feasible, but further optimization is needed before unsupervised use. Validation at varying DAB staining intensities is hampered by lack of reliability of the gold standard (i.e., visual human interpretation). Nevertheless, the DIA method has the potential to be used as a tool to assist pathologists in the analysis of IHC staining.

Correlation and colocalization of HIF-1α and pimonidazole staining for hypoxia in laryngeal squamous cell carcinomas: A digital, single-cell-based analysis.

OBJECTIVE: Tumor hypoxia results in worse local control and patient survival. We performed a digital, single-cell-based analysis to compare two biomarkers for hypoxia (hypoxia-inducible factor 1-alpha [HIF-1α] and pimonidazole [PIMO]) and their effect on outcome in laryngeal cancer patients treated with accelerated radiotherapy with or without carbogen breathing and nicotinamide (AR versus ARCON). MATERIALS AND METHODS: Immunohistochemical staining was performed for HIF-1α and PIMO in consecutive sections of 44 laryngeal cancer patients randomized between AR and ARCON. HIF-1α expression and PIMO-binding were correlated using digital image analysis in QuPath. High-density areas for each biomarker were automatically annotated and staining overlap was analyzed. Kaplan-Meier survival analyses for local control, regional control and disease-free survival were performed to predict a response benefit of ARCON over AR alone for each biomarker. RESULTS: 106 Tissue fragments of 44 patients were analyzed. A weak, significant positive correlation was observed between HIF-1α and PIMO positivity on fragment level, but not on patient level. A moderate strength correlation (r = 0.705, p < 0.001) was observed between the number of high-density staining areas for both biomarkers. Staining overlap was poor. HIF-1α expression, PIMO-binding or a combination could not predict a response benefit of ARCON over AR. CONCLUSION: Digital image analysis to compare positive cell fractions and staining overlap between two hypoxia biomarkers using open-source software is feasible. Our results highlight that there are distinct differences between HIF-1α and PIMO as hypoxia biomarkers and therefore suggest co-existence of different forms of hypoxia within a single tumor.

Prognostic imaging variables for recurrent laryngeal and hypopharyngeal carcinoma treated with primary chemoradiotherapy: A systematic review and meta-analysis.

BACKGROUND: In this systematic review, we aim to identify prognostic imaging variables of recurrent laryngeal or hypopharyngeal carcinoma after chemoradiotherapy. METHODS: A systematic search was performed in PubMed and EMBASE (1990-2020). The crude data and effect estimates were extracted for each imaging variable. The level of evidence of each variable was assessed and pooled risk ratios (RRs) were calculated. RESULTS: Twenty-two articles were included in this review, 17 on computed tomography (CT) and 5 on magnetic resonance imaging (MRI) variables. We found strong evidence for the prognostic value of tumor volume at various cut-off points (pooled RRs ranging from 2.09 to 3.03). Anterior commissure involvement (pooled RR 2.19), posterior commissure involvement (pooled RR 2.44), subglottic extension (pooled RR 2.25), and arytenoid cartilage extension (pooled RR 2.10) were also strong prognostic factors. CONCLUSION: Pretreatment tumor volume and involvement of several subsites are prognostic factors for recurrent laryngeal or hypopharyngeal carcinoma after chemoradiotherapy.

Effects of selectively assisting impaired subtasks of walking in chronic stroke survivors.

BACKGROUND: Recently developed controllers for robot-assisted gait training allow for the adjustment of assistance for specific subtasks (i.e. specific joints and intervals of the gait cycle that are related to common impairments after stroke). However, not much is known about possible interactions between subtasks and a better understanding of this can help to optimize (manual or automatic) assistance tuning in the future. In this study, we assessed the effect of separately assisting three commonly impaired subtasks after stroke: foot clearance (FC, knee flexion/extension during swing), stability during stance (SS, knee flexion/extension during stance) and weight shift (WS, lateral pelvis movement). For each of the assisted subtasks, we determined the influence on the performance of the respective subtask, and possible effects on other subtasks of walking and spatiotemporal gait parameters. METHODS: The robotic assistance for the FC, SS and WS subtasks was assessed in nine mildly impaired chronic stroke survivors while walking in the LOPES II gait trainer. Seven trials were performed for each participant in a randomized order: six trials in which either 20% or 80% of assistance was provided for each of the selected subtasks, and one baseline trial where the participant did not receive subtask-specific assistance. The influence of the assistance on performances (errors compared to reference trajectories) for the assisted subtasks and other subtasks of walking as well as spatiotemporal parameters (step length, width and height, swing and stance time) was analyzed. RESULTS: Performances for the impaired subtasks (FC, SS and WS) improved significantly when assistance was applied for the respective subtask. Although WS performance improved when assisting this subtask, participants were not shifting their weight well towards the paretic leg. On a group level, not many effects on other subtasks and spatiotemporal parameters were found. Still, performance for the leading limb angle subtask improved significantly resulting in a larger step length when applying FC assistance. CONCLUSION: FC and SS assistance leads to clear improvements in performance for the respective subtask, while our WS assistance needs further improvement. As effects of the assistance were mainly confined to the assisted subtasks, tuning of FC, SS and WS can be done simultaneously. Our findings suggest that there may be no need for specific, time-intensive tuning protocols (e.g. tuning subtasks after each other) in mildly impaired stroke survivors.