Abstract High-grade glioma (HGG) presents significant treatment challenges, particularly in identifying less responsive tumor areas to target with adaptive radiotherapy (RT). This study introduces a patient-specific computational pipeline designed to predict and spatially map HGG response to RT, focusing on regions demonstrating resilience to RT. Our computational pipeline is centered around a family of 3D biophysical models describing tumor cell invasion, proliferation, and response to RT which can be personalized for individual patients using longitudinal multiparametric MRI data. Utilizing longitudinal multiparametric MRI data reporting on cell density (diffusion weighted MRI) and the extent of disease (T2-FLAIR and contrast-enhanced T1-weighted MRIs) from 21 patients, ranging from baseline to three weeks into RT, we calibrated a family of biophysical models to forecast tumor cell density up to the final week of RT and one-month follow-up. By employing pareto tail analysis on empirical cumulative distribution function fits of voxel-wise changes in apparent diffusion coefficient maps (ΔADC) between the future images (end of RT, one-month follow-up) and pre-treatment images, we identified areas of increased tumor cellularity on both the ground truth and model predicted cellularity maps. The analysis differentiated tumor voxels into “non-responding” and “responding” categories, based on cutoffs determined through ROC analysis across various conditions—including assessments of both enhancing and non-enhancing tumor regions at different time points. Optimal ΔADC cutoffs for predicted and measured data were established, leading to high area under curve (AUC; ranging from 0.803 to 0.934) values and significant Youden’s index scores (ranging from 0.526 to 0.765) across all conditions. Our findings demonstrate the feasibility of this novel pipeline in delineating biologically relevant RT targets, offering a potential tool for enhancing the efficacy of adaptive RT strategies in treating HGG.
Read full abstract