Singleton, Kyle

Abstract Glioblastoma (GBM) is a devastating primary brain tumor known for its heterogeneity with a median survival of 15 months. Clinical imaging remains the primary modality to assess brain tumor response, but it is nearly impossible to distinguish between tumor growth and treatment response. Ki67 is a marker of active cell proliferation that shows inter- and intra-patient heterogeneity and should change under many therapies. In this work, we assessed the utility of a semi-supervised deep learning approach for regionally predicting high-vs-low Ki67 in GBM patients based on MRI. We used both labeled and unlabeled datasets to train the model. Labeled data included 114 MRI-localized biopsies from 43 unique GBM patients with available immunohistochemistry Ki67 labels. Unlabeled data included nine repeat routine pretreatment paired scans of newly-diagnosed GBM patients acquired within three days. Data augmentation techniques were utilized to enhance the size of our data and increase generalizability. Data was split between training, validation, and testing sets using 65-15-20 percent ratios. Model inputs were 16x16x3 patches around biopsies on T1Gd and T2 MRIs for labeled data, and around randomly selected patches inside the T2 abnormal region for unlabeled data. The network was a 4-conv layered VGG-inspired architecture. Training objective was accurate prediction of Ki67 in labeled patches and consistency in predictions across repeat unlabeled patches. We measured final model accuracy on held-out test samples. Our promising preliminary results suggest potential for deep learning in deconvolving the spatial heterogeneity of proliferative GBM subpopulations. If successful, this model can provide a non-invasive readout of cell proliferation and reveal the effectiveness of a given cytotoxic therapy dynamically during the patient's routine follow up. Further, the spatial resolution of our approach provides insights into the intra-tumoral heterogeneity of response which can be related to heterogeneity in localization of therapies (e.g. radiation therapy, drug dose heterogeneity).

Abstract Choosing effective chemotherapies for intravenous delivery to brain tumors is challenging, especially given the protective nature of the blood brain barrier (BBB). Connecting drug distribution to non-invasive, pre-surgical magnetic resonance imaging (MRI) could allow for predictive insight into drug distribution. In a previous study, we found that T2Gd images were predictive of a low BBB penetrant drug (Cefazolin), and FLAIR images were predictive of a high BBB penetrant drug (Levetiracetam). While these results are promising, we further seek to explore how advanced MRI sequences might inform image-based models of drug distribution. Prior to surgery, we acquired advanced dynamic contrast enhanced (DCE) and diffusion weighted imaging (DWI) MRI sequences for eight brain tumor patients (7 gliomas and 1 metastatic adenocarcinoma) in addition to the anatomic MRIs. All resulting quantitative maps and acquired images were co-registered. Prior to incision, patients received injections of cefazolin and levetiracetam. Next, multiple blood samples and biopsies were collected during surgery. Biopsies and plasma samples were analyzed for drug concentration using liquid chromatography mass spectrometry (LCMS), and biopsy drug levels were reported as Brain-Plasma Ratio (BPR). Mean image intensity was extracted from a 15x15 voxel window surrounding the biopsy location. We performed linear regression analyses to determine which combination of images were predictive of BPR. We found that considering quantitative imaging improved our initial ability to predict BPR for both drugs. For cefazolin, the third diffusion tensor eigenvalue (L3) map was significantly correlated with BPR (p< 0.001, R2= 0.36). For levetiracetam, the best model consisted of a combination of images and maps with the L3 map and the isotropic diffusion map (P) being the most influential (p= 0.001, R2= 0.63). Advanced MRI-based modeling is a promising tool for forecasting drug distribution in brain tumors and could be of great importance for understanding efficacy and selecting therapeutic strategies.

Abstract OBJECTIVE Glioblastoma (GBM) is the most common malignant primary brain tumor in adults, with males more commonly affected than females(1.6:1). Despite advancements in treatments, prognosis is dismal with a median overall survival of 15 months. Our aim was to investigate sex as a variable in GBM patient survival after receiving incremental levels of standard-of-care treatment regimens – different extents of surgical resection and different numbers of cycles of adjuvant temozolomide chemotherapy. METHODS Drawing from our extensive multi-institutional brain tumor repository, we investigated GBM subjects with overall survival (OS), extent of resection (EOR), number of temozolomide (TMZ) cycles, and sex data (n=620, males: n=387, females: n=233). Cox proportional hazard ratios were computed to investigate the multivariable predictive value of the patient variables with OS. Patients were then divided into groups based on their sex, EOR (either biopsy, subtotal resection (STR) or gross total resection (GTR)), and TMZ cycles (I: < 6 cycles, II: 7-11 cycles and III: >12 cycles). RESULTS We observed that STR was beneficial for females (HR=0.52; CI=0.33-0.83; p-value=0.013), while for males the benefit was not detected (HR=0.73; CI=0.46-1.15; p-value=0.173) for STR but was detectable for GTR (HR=0.58, CI=0.37-0.90; p-value=0.014). Females receiving 7-11 cycles of TMZ showed a survival benefit (HR=0.52; CI=0.12-0.53; p-value=0.048) while males in the same group did not (HR=0.74; CI=0.46-1.19; p-value=0.21), in comparison to those in group I of TMZ cycles. No sex differences were identified in patients receiving < =6 cycles or >=12 cycles. CONCLUSION Together, our results contribute to the growing literature that sex differences exist in GBM patients, even in response to standard-of-care therapies. This should be accounted for when designing clinical trials for GBM so that we may advance our pursuit to deliver personalized medicine.

Abstract Background Accurate assessments of patient response to therapy are a critical component of personalized medicine. In glioblastoma (GBM), the most aggressive form of brain cancer, tumor growth dynamics are heterogenous across patients, complicating assessment of treatment response. This study aimed to analyze days gained (DG), a burgeoning model-based dynamic metric, for response assessment in patients with recurrent GBM who received bevacizumab-based therapies. Methods DG response scores were calculated using volumetric tumor segmentations for patients receiving bevacizumab with and without concurrent cytotoxic therapy (N = 62). Kaplan–Meier and Cox proportional hazards analyses were implemented to examine DG prognostic relationship to overall (OS) and progression-free survival (PFS) from the onset of treatment for recurrent GBM. Results In patients receiving concurrent bevacizumab and cytotoxic therapy, Kaplan–Meier analysis showed significant differences in OS and PFS at DG cutoffs consistent with previously identified values from newly diagnosed GBM using T1-weighted gadolinium-enhanced magnetic resonance imaging (T1Gd). DG scores for bevacizumab monotherapy patients only approached significance for PFS. Cox regression showed that increases of 25 DG on T1Gd imaging were significantly associated with a 12.5% reduction in OS hazard for concurrent therapy patients and a 4.4% reduction in PFS hazard for bevacizumab monotherapy patients. Conclusion DG has significant meaning in recurrent therapy as a metric of treatment response, even in the context of anti-angiogenic therapies. This provides further evidence supporting the use of DG as an adjunct response metric that quantitatively connects treatment response and clinical outcomes.

Abstract Determining effective therapies for patients with recurrent glioblastoma remains difficult especially in small early phase clinical trials. Currently utilized estimators of treatment response rely on two-dimensional measurements of tumor progression (e.g. RANO, RECIST, MacDonald) that have often fallen short in reliably connecting patients who respond with positive survival outcomes. This has prompted consideration of advanced imaging and volumetric measurements by the RANO working group and others. Our proposed approach, Days Gained (DG), uses volumetric measurements to calculate a personalized response metric based on change in tumor growth characteristics pre/post-treatment. DG has shown promise for connecting imageable treatment response with clinically meaningful (statistically significant) survival benefit when applied in upfront and recurrent therapy settings, even in the context of coincident bevacizumab treatment. In this study, we seek to consider a new cohort of patients receiving bevacizumab in combination with other novel therapies to distinguish any subpopulation of patients with meaningful response. In this work, we identified a subset of 16 patients from bevacizumab clinical trials with data available for DG analysis requiring serial imaging (T1Gd and FLAIR) prior to and after the first cycle of therapy. Each trial evaluated bevacizumab alone versus bevacizumab plus either dasatinib or TRC105 in recurrent disease. Discrimination of time to progression (TTP) was evaluated at DG thresholds identified from our prior analysis of patients treated with bevacizumab. Using previously identified DG thresholds for clinically significant treatment benefit, in our cohort, DG measured by serial FLAIR imaging was significant in predicting longer TTP (120 DG, p=0.014), but T1Gd DG scores did not attain significance (p=0.630). Given the limited insights in connecting response metrics to accurately predicted outcomes, these data further support the use of DG as a burgeoning marker of treatment response (or lack of response) that may be useful for more routine integration into clinical trials.

Abstract BACKGROUND Conventional magnetic resonance imaging (MR) guides patient care in GBM. However, there is mounting awareness that MR enhancement is non-specific reflecting either tumor progression or non-tumoral inflammatory changes. Histological evaluation of GBM is held as the gold standard for disease assessment. However, the invasiveness of this methodology and the sample sparsity limit its usefulness. Methods to infer histological underpinnings of MRI are needed to improve clinical care. METHODS A transfer learning mixed effects model based on T1Gd and FLAIR MR voxel based image features trained and cross-validated to predict FPKM values of CD68, CASP3, CD8A and IL13RA2. Training data included RNAseq from 38 image-localized biopsies from 15 newly-diagnosed GBM patients. These models were then applied to two independent patients, chosen based on therapy and the quality of image registration, at three different time points, just prior to receiving IL13RA2 targeted CAR-T therapy for rGBM, after 2 cycles and after 4 cycles, resulting in a voxel-based prediction of the relative expression levels of these genes. RESULTS Cross-validation of the proposed machine learning models demonstrated a Pearson Correlation coefficient between predicted and observed FPKM values of 0.89 (CD68), 0.92 (CASP3), 0.93 (CD8A), and 0.88 (IL13RA2). When the models were applied to the two CAR-T patients, CD68-modeled expression levels were seen to increase throughout therapy for one patient, while remaining stable for the other. Survival from first CAR-T infusion was, respectively, 412 and 83 days suggesting that the increased CD68 was indicative of therapeutic effectivity. Further, the spatial activity predictions were consistent with expected therapeutic action as the correlation coefficient between CASP3 and the product of IL13R2 and CD8A expression was 0.81, suggesting the T-cells targeting IL13RA2 were inducing cell death. CONCLUSIONS Preliminary results with our model highlights the potential of spatially resolved radiomic maps to provide insight into regional therapeutic effectivity.

Abstract Recurrent glioblastoma is challenging to distinguish from so called “treatment effect” on routine clinical imaging. Further, within tumor heterogeneity reveals that some regions can be histologically dominated by tumor progression whilst others can be dominated by secondary effects of treatment response. Apparent tumor progression on MRI can be very difficult to manage clinically as it is unclear the degree to which the imaging changes are actually tumor progression vs response to treatment (including inflammatory response and necrosis). In this analysis, we study a unique cohort of patients for whom image localized-biopsies reveal heterogeneity in response vs progression. Our dataset included 70 biopsy samples from 32 patients with GBM each histolopatholgically characterized for tumor abundance vs immune infiltrate. Six multiparametric MRI contrasts were available, including T1, T1gd, T2, FLAIR, SWI, and ADC. Images were co-registered. Radiomic (statistical + texture) features were extracted from the region of six image contrasts locally matched with each biopsy sample. Machine learning models were built to predict each biomarker using radiomic features. Leave-one-out cross validation was used to evaluate the prediction accuracy. Radiomic features were found to be informative to the prediction of biomarkers. ANOVA tests show significant improvement of using radiomic features compared with the null model. The prediction accuracy was higher when considering the biomarkers on a binary scale using the median as the cutoff than on a numerical scale. Spatially-informed radiomics models for tumor progression vs treatment effect are possible and can play an instrumental role in navigating confounding imaging changed common during treatment progression.

Abstract BACKGROUND An important challenge in radiomics research is reproducibility. Images are collected on different image scanners and protocols, which introduces significant variability even for the same type of image across institutions. In the present proof-of-concept study, we address the reproducibility issue by using domain adaptation – an algorithm that transforms the radiomic features of each new patient to align with the distribution of features formed by the patient samples in a training set. METHOD Our dataset included 18 patients in training with a total of 82 biopsy sample. The pathological tumor cell density was available for each sample. Radiomic (statistical + texture) features were extracted from the region of six image contrasts locally matched with each biopsy sample. A Gaussian Process (GP) classifier was built to predict tumor cell density using radiomic features. Another 6 patients were used to test the training model. These patients had a total of 31 biopsy samples. The images of each test patient were purposely normalized using a different approach, i.e., using the CSF instead of the whole brain as the reference. This was to mimic the practical scenario of image source discrepancy between different patients. Domain adaptation was applied to each test patient. RESULTS Among the 18 training patients, the leave-one-patient-out cross validation accuracy is 0.81 AUC, 0.78 sensitivity, and 0.83 specificity. When the trained model was applied to the 6 test patients (purposely normalized using a different approach than that of the training data), the accuracy dramatically reduced to 0.39 AUC, 0.08 sensitivity, and 0.61 specificity. After using domain adaption, the accuracy improved to 0.68 AUC, 0.62 sensitivity, and 0.72 specificity. CONCLUSION We provide candidate enabling tools to address reproducibility in radiomics models by using domain adaption algorithms to account for discrepancy of the images between different patients.

Abstract INTRODUCTION The quantification of intratumoral heterogeneity – through radiomics-based approaches - can help resolve the regionally distinct genetic drug targets that may co-exist within a single Glioblastoma (GBM) tumor. While this offers potential diagnostic value under the paradigm of individualized oncology, clinical decision-making must also consider the degree of uncertainty associated with each model. In this study, we evaluate the performance of a novel machine-learning (ML) algorithm, called Gaussian Process (GP) modeling, that can quantify the impact of multiple sources of uncertainty in ML model development and prediction accuracy, including variabilities in the copy number measurement, radiomics features, training sample characteristics, and training sample size. METHOD We collected 95 image-localized biopsies from 25 primary GBM patients. We coregistered stereotactic locations with preoperative multi-parametric MRI features (conventional MRI, DSC perfusion, Diffusion Tensor Imaging) to generate spatially matched pairs of MRI and copy number variants (CNV) for for each biopsy. We developed a Gaussian Process (GP) model to predict CNV for Epidermal Growth Factor Receptor (EGFR) based on MRI radiomic features in each patient. We used leave-one-patient-out cross validation to quantify prediction accuracy and model uncertainty. Spatial prediction and uncertainty (p-value) maps were overlaid to help visualize regional genetic variation of EGFR and uncertainty of the radiomic predictions. RESULT: The initial GP radiomics model for EGFR amplification (CNV > 3.5) produced a sensitivity of 0.8 and specificity of 0.8. Samples/regions associated with high uncertainty (p-value >0.05) correlated with either 1) extrapolation of radiomic features from the training set-defined feature space or 2) insufficient training samples in the feature space. CONCLUSION We present a ML-based model that quantifies spatial genetic heterogeneity in GBM, while also estimating model uncertainties that result from multi-source data variabilities. This approach lays the groundwork for prospective clinical integration of modeling-based diagnostic approaches in the paradigm of individualized medicine.