Samuel A Bobholz1, Alisha Hoefs1, Jordyn Hamburger1, Allison K Lowman1, Savannah R Duenweg2, Aleksandra Winiarz2, Margaret Stebbins2, Fitzgerald Kyereme1, Jennifer Connelly3, Dylan Coss4, Wade M Mueller5, Mohit Agarwal1, Anjishnu Banerjee6, and Peter S LaViolette1,7
1Radiology, Medical College of Wisconsin, Milwaukee, WI, United States, 2Biophysics, Medical College of Wisconsin, Milwaukee, WI, United States, 3Neurology, Medical College of Wisconsin, Milwaukee, WI, United States, 4Pathology, Medical College of Wisconsin, Milwaukee, WI, United States, 5Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, United States, 6Biostatistics, Medical College of Wisconsin, Milwaukee, WI, United States, 7Biomedical Engineering, Medical College of Wisconsin, Milwaukee, WI, United States
Synopsis
Keywords: Tumors, Tumor
We tested the hypothesis that autopsy-based
radio-pathomic maps of glioblastoma pathology reveal
distinct phenotypes (hypercellular, hypocellular, hybrid, and well-circumscribed fronts) that differ in patient survival and bevacizumab treatment response. Patients with tumor invasion beyond contrast showed worse survival outcomes compared to patients with well-circumscribed tumors. Additionally, patients with hypocellular components of the non-enhancing front selectively benefit from bevacizumab treatment, with an observable reduction in the hypocellular volume over the course of bevacizumab use.
Introduction
Multi-parametric
MRI is used to guide treatment delivery and monitor disease progression in
glioblastoma patients. However traditional imaging signatures such as T1-weighted
gadolinium contrast enhancement and T2-weighted fluid attenuated inversion
recovery (FLAIR) hyperintensity are known to be confounded by the presence of treatments
such as bevacizumab and other antiangiogenic agents. Our previously developed autopsy-based
non-invasive maps of tumor pathology have been used to identify areas of tumor
beyond traditional imaging signatures in the post-treated state, which may be
able to more precisely monitor treatment response later in a patient’s clinical
history. This study used characteristics identified by autopsy tissue-based
radio-pathomic maps of tumor pathology to examine heterogenous characteristics
of bevacizumab response in glioma patients.Methods
Radio-pathomic
models developed in a recently published manuscript were the primary contrast
of interest in this study1,2 (see Figure 1 for an overview of the
radio-pathomic mapping process). Briefly, pre- and post-contrast T1-weighted
images (T1, T1C), T2-weighted FLAIR images, and apparent diffusion coefficient
(ADC) images were used as input to predict cellularity, extracellular fluid
density, cytoplasm density, and tumor probability using tissue samples aligned
to the last clinical imaging prior to death as ground truth. A training dataset
of 43 patients was used to train a bagging regression ensemble using 5 by 5
voxel tiles from the MRI as input and voxelwise pathological charact-eristics
as labels, with a held-out test set of 22 subjects used for model validation. These
models demonstrated good quantitative performance and successfully identified
areas of tumor outside the contrast-enhancing region. A system of phenotypes
was then developed based on the visual appearance of the non-enhancing tumor
front on radio-pathomic maps: Hypercellular Front, where portions of
high cellularity extent beyond the contrast-enhancing margin; Hypocellular Front,
where portions of low cellularity and high extracellular fluid density extend
beyond the contrast-enhancing margin and are thought to indicate areas of edema
or hypoxia; Hybrid Front, where both hypocellular and hypercellular
regions spread beyond the primary enhancing area, and Well-Circumscribed
tumors, where all abnormal pathological findings occur within the contrast-enhancing
area (see Figure 2 for examples of each phenotype). Radio-pathomic maps
were then generated for an independent dataset of 80 GBM cases with imaging prior
to first surgery and graded using the described phenotyping system. Kaplan-Meier
analysis was then used to compare differences in survival between phenotypes,
as well as patients who had and had not received bevacizumab treatment within
each phenotype. Additionally, masks for hyper/hypocellular presence were drawn
for a subset of 26 GBM patients imaged over the course of bevacizumab
treatment, and linear mixed effect models were used to assess for longitudinal
changes in hyper/hypocellular volumes associated with treatment.Results
Figure
3 shows the
Kaplan-Meier curves for overall survival within each phenotype. Well-circumscribed
patients showed significant/trending increases in survival compared to Hyercellular
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). Figure 4 shows survival curves
within patients who have and have not received bevacizumab treatment, divided
by phenotype. 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). Figure 5 shows the hypocellular volume
per subject with respect to days since initiating bevacizumab treatment. Hypocellular
volumes decreased by an average 50.52 mm3 per day of bevacizumab
treatment (p=0.002).Discussion
This study
investigated radio-pathomic phenotypes of glioblastoma patients to identify
differences in patient survival and bevacizumab treatment response. We found
that patients with active tumor presence beyond the contrast enhancing region
survived less long than patients with tumors well-circumscribed by contrast
enhancement, concurring with prior findings observed in our autopsy-based
dataset. This suggests that occult glioma invasion may be a hallmark of more aggressive
disease, requiring data-driven approaches with rich pathological ground truth
to visualize non-invasively. This study also found that patients with
hypocellular presence (Hypocellular front and Hybrid front) demonstrate survival
benefit from bevacizumab treatment not seen in patients without hypocellular
presence. This indicates that these patients may present with pathological
features that allow for more efficacious treatment with antiangiogenic agents,
and that our non-invasive maps of tumor characteristics could inform clinical
decision-making regarding bevacizumab use. Hypocellular regions were also seen
to decrease over the course of bevacizumab treatment, further suggesting that
these maps may be useful in monitoring antiangiogenic treatment response and
underscoring their utility as a response assessment tool for non-enhancing
portions of tumor. Future research using this technique may be able to identify
subsets of patients within retrospective clinical trials that selectively
respond to failed treatments.Acknowledgements
No acknowledgement found.References
1. Bobholz, S. A., Lowman, A.
K., Connelly, J. M., Duenweg, S. R., Winiarz, A., Brehler, M., … LaViolette, P.
S. (2022). Non-invasive tumor probability maps developed using autopsy tissue
identify novel areas of tumor beyond the imaging-defined margin. MedRxiv,
2022.08.17.22278910. https://doi.org/10.1101/2022.08.17.22278910
2. Bobholz, S. A., Lowman, A.
K., Brehler, M., Kyereme, F., Duenweg, S. R., Sherman, J., … LaViolette, P. S.
(2022). Radio-Pathomic Maps of Cell Density Identify Brain Tumor Invasion
beyond Traditional MRI-Defined Margins. AJNR. American Journal of
Neuroradiology, 43(5), 682–688. https://doi.org/10.3174/ajnr.A7477