The ISMRM Open Science Initiative for Perfusion Imaging (OSIPI): Results from the OSIPI-Dynamic Contrast-Enhanced challenge.

PURPOSE: K trans $$ {K}^{\mathrm{trans}} $$ has often been proposed as a quantitative imaging biomarker for diagnosis, prognosis, and treatment response assessment for various tumors. None of the many software tools for K trans $$ {K}^{\mathrm{trans}} $$ quantification are standardized. The ISMRM Open Science Initiative for Perfusion Imaging-Dynamic Contrast-Enhanced (OSIPI-DCE) challenge was designed to benchmark methods to better help the efforts to standardize K trans $$ {K}^{\mathrm{trans}} $$ measurement. METHODS: A framework was created to evaluate K trans $$ {K}^{\mathrm{trans}} $$ values produced by DCE-MRI analysis pipelines to enable benchmarking. The perfusion MRI community was invited to apply their pipelines for K trans $$ {K}^{\mathrm{trans}} $$ quantification in glioblastoma from clinical and synthetic patients. Submissions were required to include the entrants' K trans $$ {K}^{\mathrm{trans}} $$ values, the applied software, and a standard operating procedure. These were evaluated using the proposed OSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score defined with accuracy, repeatability, and reproducibility components. RESULTS: Across the 10 received submissions, the OSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score ranged from 28% to 78% with a 59% median. The accuracy, repeatability, and reproducibility scores ranged from 0.54 to 0.92, 0.64 to 0.86, and 0.65 to 1.00, respectively (0-1 = lowest-highest). Manual arterial input function selection markedly affected the reproducibility and showed greater variability in K trans $$ {K}^{\mathrm{trans}} $$ analysis than automated methods. Furthermore, provision of a detailed standard operating procedure was critical for higher reproducibility. CONCLUSIONS: This study reports results from the OSIPI-DCE challenge and highlights the high inter-software variability within K trans $$ {K}^{\mathrm{trans}} $$ estimation, providing a framework for ongoing benchmarking against the scores presented. Through this challenge, the participating teams were ranked based on the performance of their software tools in the particular setting of this challenge. In a real-world clinical setting, many of these tools may perform differently with different benchmarking methodology.

Radio-pathomic maps of glioblastoma identify phenotypes of non-enhancing tumor infiltration associated with bevacizumab treatment response.

BACKGROUND: Autopsy-based radio-pathomic maps of glioma pathology have shown substantial promise inidentifying areas of non-enhancing tumor presence, which may be able to differentiate subsets of patients that respond favorably to treatments such as bevacizumab that have shown mixed efficacy evidence. We tested the hypthesis that phenotypes of non-enhancing tumor fronts can distinguish between glioblastoma patients that will respond favorably to bevacizumab and will visually capture treatment response. METHODS: T1, T1C, FLAIR, and ADC images were used to generate radio-pathomic maps of tumor characteristics for 79 pre-treatment patients with a primary GBM or high-grade IDH1-mutant astrocytoma for this study. Novel phenotyping (hypercellular, hypocellular, hybrid, or well-circumscribed front) of the non-enhancing tumor front was performed on each case. Kaplan Meier analyses were then used to assess differences in survival and bevacizumab efficacy between phenotypes. Phenotype compartment segmentations generated longitudinally for a subset of 26 patients over the course of bevacizumab treatment, where a mixed effect model was used to detect longitudinal changes. RESULTS: Well-Circumscribed patients showed significant/trending increases in survival compared to Hypercellular Front (HR = 2.0, p = 0.05), Hypocellular Front (HR = 2.02, p = 0.03), and Hybrid Front tumors (HR = 1.75, p = 0.09). Only patients with hypocellular or hybrid fronts showed significant survival benefits from bevacizumab treatment (HR = 2.35, p = 0.02; and HR = 2.45, p = 0.03, respectively). Hypocellular volumes decreased by an average 50.52 mm3 per day of bevacizumab treatment (p = 0.002). CONCLUSION: Patients with a hypocellular tumor front identified by radio-pathomic maps showed improved treatment efficacy when treated with bevacizumab, and reducing hypocellular volumes over the course of treatment may indicate treatment response.

Morphologic features of prostate cancer-encased native vessels: An image analysis study.

BACKGROUND: Vasculature plays a crucial role in the progression of prostate cancer (PC). Changes to the prostatic native vessels have not been studied since 2000 when Garcia et al. demonstrated marked media hypercellularity and increased artery thickness in prostatic native arteries within PC. We aim to further evaluate and characterize prostatic native vessels with a more accurate method with the use of virtual slides and digital analysis. DESIGN: Pathologist-annotated whole-mount digital slides from 96 entirely submitted prostatectomies were annotated for PC (color-coded by Gleason) using Omero platform. A subset of 44 cases met criteria for further analysis of media thickness, cellularity, and wall thickness to lumen ratio. Cases were included based on containing ≥5 native arteries (≥100 µm diameter) encased on at least 3 sides by PC, with vessels (≥100 µm diameter) designated as controls if they were ≥ 1000 µm away from PC. Annotated vessels were segmented and processed using Matlab 2023b. Mean media thickness (corrected for oblique sections), media: lumen ratio (based on numbers of pixels), and media cellularity (nuclei count) were studied by analysis with SPSS by linear mixed model with nested random effects for subject and slide to account for repeated measures. RESULTS: Vessels encased by PC showed greater media thickness (p=0.02), cellularity (p=0.02) and wall thickness/lumen ratio (p= <0.001) compared to vessels away from PC. These values showed an increasing trend according to stage in cellularity (p=0.14), media thickness (p=0.12) and wall thickness/ lumen ratio (p= 0.33) with higher stage (pT3). A Gleason group comparison showed a borderline-significant gradewise trend when analyzing wall thickness/lumen ratio (p=0.06). Grade 5 emerged as significantly different (p=0.02) from grades 3 or 4 non-cribriform. CONCLUSIONS: Similar to the 2000 study, increased media thickness and hypercellularity of vessels encased by PC were evident compared to controls. Borderline grade-dependent increased vessel cellularity changes were seen, suggesting a possible role in PC progression; the predictive value of these changes for outcome is uncertain. Whether the etiology of changes reflects locally increased intravascular pressure of vessels within tumor should be investigated.

Radio-pathomic maps of glioblastoma identify phenotypes of non-enhancing tumor infiltration associated with bevacizumab treatment response.

Background: Autopsy-based radio-pathomic maps of glioma pathology have shown substantial promise inidentifying areas of non-enhancing tumor presence, which may be able to differentiate subsets of patients that respond favorably to treatments such as bevacizumab that have shown mixed efficacy evidence. We tested the hypthesis that phenotypes of non-enhancing tumor fronts can distinguish between glioblastoma patients that will respond favorably to bevacizumab and will visually capture treatment response. Methods: T1, T1C, FLAIR, and ADC images were used to generate radio-pathomic maps of tumor characteristics for 79 pre-treatment patients with a primary GBM or high-grade IDH1-mutant astrocytoma for this study. Novel phenotyping (hypercellular, hypocellular, hybrid, or well-circumscribed front) of the non-enhancing tumor front was performed on each case. Kaplan Meier analyses were then used to assess differences in survival and bevacizumab efficacy between phenotypes. Phenotype compartment segmentations generated longitudinally for a subset of 26 patients over the course of bevacizumab treatment, where a mixed effect model was used to detect longitudinal changes. Results: Well-Circumscribed patients showed significant/trending increases in survival compared to Hypercellular Front (HR = 2.0, p = 0.05), Hypocellular Front (HR = 2.02, p = 0.03), and Hybrid Front tumors (HR = 1.75, p = 0.09). Only patients with hypocellular or hybrid fronts showed significant survival benefits from bevacizumab treatment (HR = 2.35, p = 0.02; and HR = 2.45, p = 0.03, respectively). Hypocellular volumes decreased by an average 50.52 mm3 per day of bevacizumab treatment (p = 0.002). Conclusion: Patients with a hypocellular tumor front identified by radio-pathomic maps showed improved treatment efficacy when treated with bevacizumab, and reducing hypocellular volumes over the course of treatment may indicate treatment response.

Noninvasive Autopsy-Validated Tumor Probability Maps Identify Glioma Invasion Beyond Contrast Enhancement.

BACKGROUND AND OBJECTIVES: This study identified a clinically significant subset of patients with glioma with tumor outside of contrast enhancement present at autopsy and subsequently developed a method for detecting nonenhancing tumor using radio-pathomic mapping. We tested the hypothesis that autopsy-based radio-pathomic tumor probability maps would be able to noninvasively identify areas of infiltrative tumor beyond traditional imaging signatures. METHODS: A total of 159 tissue samples from 65 subjects were aligned to MRI acquired nearest to death for this retrospective study. Demographic and survival characteristics for patients with and without tumor beyond the contrast-enhancing margin were computed. An ensemble algorithm was used to predict pixelwise tumor presence from pathological annotations using segmented cellularity (Cell), extracellular fluid, and cytoplasm density as input (6 train/3 test subjects). A second level of ensemble algorithms was used to predict voxelwise Cell, extracellular fluid, and cytoplasm on the full data set (43 train/22 test subjects) using 5-by-5 voxel tiles from T1, T1 + C, fluid-attenuated inversion recovery, and apparent diffusion coefficient as input. The models were then combined to generate noninvasive whole brain maps of tumor probability. RESULTS: Tumor outside of contrast was identified in 41.5% of patients, who showed worse survival outcomes (hazard ratio = 3.90, P < .001). Tumor probability maps reliably tracked nonenhancing tumor on a range of local and external unseen data, identifying tumor outside of contrast in 69% of presurgical cases that also showed reduced survival outcomes (hazard ratio = 1.67, P = .027). CONCLUSION: This study developed a multistage model for mapping gliomas using autopsy tissue samples as ground truth, which was able to identify regions of tumor beyond traditional imaging signatures.