Pixel-Wise Prediction of Glioblastoma Recurrence Using Mid-Treatment MRI and Deep Learning
Sara Rezvanjou, Laurentius Oscar Osapoetra, Daniel Moore-Palhares, Liam S. P. Lawrence, Angus Z. Lau, Arjun Sahgal, Gregory J. CzarnotaIntroduction: Glioblastoma (GBM) commonly recurs close to the original gross tumor volume (GTV), suggesting treatment-resistant subregions may be present before radiologically visible progression. This study investigated whether deep learning on apparent diffusion coefficient (ADC) maps could generate pixel-wise probability maps of treatment-resistant regions within the GTV. Methods: The training/validation set included images of 64 in-house patients after radiation fraction 10 (Fx10) and 40 external patients before treatment used as an auxiliary cross-domain training cohort; the held-out test set included 16 in-house Fx10 patients. Resistant regions were defined as the overlap between the GTV from scans acquired at radiation fraction 10 and subsequent recurrence after finishing treatment. Models were trained on 2D GTV-cropped slices with patient-level separation to prevent data leakage. Single- and multimodal MRI inputs were compared, and ADC-based models (including baseline UNet, pretrained ResNet34, BiomedCLIP-assisted UNet, and a BLIP image-text UNet) were later evaluated. Performance was evaluated using conventional metrics, including ROC AUC, PR AUC, F1-score, recall, specificity, and precision. Results: ADC was the most informative single modality, achieving 69% and 80% for Receiver Operating Characteristic (ROC)- and Precision-Recall (PR)-Area Under the Curve (AUC), respectively. In the ADC-based architecture comparison, no model was consistently superior across all metrics. The baseline UNet achieved the highest accuracy, F1-score, and recall, whereas pretrained ResNet34 and BLIP-assisted models showed slightly higher ROC AUC than the baseline UNet (both 70%). Error analysis showed that performance decreased substantially for small or sparse resistant regions, and false-negative errors were enriched near resistant-region and GTV boundaries. Conclusions: ADC-based deep learning showed feasibility for pixel-wise spatial mapping of treatment-resistant glioblastoma regions within the GTV. A UNet trained from scratch provided the strongest sensitivity/F1 trade-off, whereas pretrained and VLM-assisted models provided only metric-specific benefits. Model comparisons showed that BiomedCLIP-assisted and BLIP image-text UNet models did not clearly improve fine spatial segmentation performance over the baseline UNet.