Citations

Abstract Background : Sex is recognized as a significant determinant of outcome among glioblastoma patients, but the relative prognostic importance of glioblastoma features has not been thoroughly explored for sex differences. Methods : Combining multi-modal MR images, biomathematical models, and patient clinical information, this investigation assesses which pretreatment variables have a sex-specific impact on the survival of glioblastoma patients (299 males and 195 females). Results : Among males, tumor (T1Gd) radius was a predictor of overall survival (HR=1.027, p=0.044). Among females, higher tumor cell net invasion rate was a significant detriment to overall survival (HR=1.011, p<0.001). Female extreme survivors had significantly smaller tumors (T1Gd) (p=0.010 t-test), but tumor size was not correlated with female overall survival (p=0.955 CPH). Both male and female extreme survivors had significantly lower tumor cell net proliferation rates than other patients (M p=0.004, F p=0.001, t-test). Conclusion : Despite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes.

Abstract Background : Sex is recognized as a significant determinant of outcome among glioblastoma patients, but the relative prognostic importance of glioblastoma features has not been thoroughly explored for sex differences. Methods : Combining multi-modal MR images, biomathematical models, and patient clinical information, this investigation assesses which pretreatment variables have a sex-specific impact on the survival of glioblastoma patients (299 males and 195 females). Results : Among males, tumor (T1Gd) radius was a predictor of overall survival (HR=1.027, p=0.044). Among females, higher tumor cell net invasion rate was a significant detriment to overall survival (HR=1.011, p<0.001). Female extreme survivors had significantly smaller tumors (T1Gd) (p=0.010 t-test), but tumor size was not correlated with female overall survival (p=0.955 CPH). Both male and female extreme survivors had significantly lower tumor cell net proliferation rates than other patients (M p=0.004, F p=0.001, t-test). Conclusion : Despite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes.

Abstract Background : Sex is recognized as a significant determinant of outcome among glioblastoma patients, but the relative prognostic importance of glioblastoma features has not been thoroughly explored for sex differences. Methods : Combining multi-modal MR images, biomathematical models, and patient clinical information, this investigation assesses which pretreatment variables have a sex-specific impact on the survival of glioblastoma patients (299 males and 195 females). Results : Among males, tumor (T1Gd) radius was a predictor of overall survival (HR=1.027, p=0.044). Among females, higher tumor cell net invasion rate was a significant detriment to overall survival (HR=1.011, p<0.001). Female extreme survivors had significantly smaller tumors (T1Gd) (p=0.010 t-test), but tumor size was not correlated with female overall survival (p=0.955 CPH). Both male and female extreme survivors had significantly lower tumor cell net proliferation rates than other patients (M p=0.004, F p=0.001, t-test). Conclusion : Despite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes.

Abstract INTRODUCTION Neuroimaging dogma for glioblastoma asserts that hyperintensity on T1Gd MRI reveals the bulk of the tumor, while T2/FLAIR signal indicates edema. However, it is unclear whether this edema results from immune response or increased tumor cells. Further, one significant driver of the known sex differences in glioblastoma may be differences in immune response, due to the X-linkage of many immune genes. Based on this, we hypothesized that assumptions regarding tumor cellularity in T2/FLAIR images should be tailored to the biological sex of the patient. METHODS Using a retrospective cohort of 18 primary glioblastoma patients receiving multiple image-localized biopsies (82 total) and standard MRI, we assessed: distance of biopsy from T1Gd and T2 areas; a pathologist’s score of percent tumor cell density; and an imaging-based invasion metric, D/ρ. This metric is derived from the biomathematical Proliferation-Invasion model of glioma growth, which features two parameters, net growth rate (ρ) and net invasion rate (D). Their ratio D/ρ is related to degree of invasion, and is estimated from volumetric measurements of MRI abnormalities. Additionally, 25 patient-derived xenograft models implanted in females were grown until moribund, at which point brains were excised and stained for DAPI (to show all cells) and Lamin (to highlight tumor cells). Image processing of lamin-stained sections defines contours of intensity correlating with cell density. RESULTS Outside both the T1Gd and T2 region, male patient biopsies had higher tumor cell densities than females. Males also tended to have higher invasion metrics. Although each set derived from different patients, preclinical metrics of invasion were positively correlated with clinical invasion in females but negatively correlated in males. CONCLUSION Our preliminary finding that cell distribution patterns correlate with imaging metrics differently between the sexes supports the hypothesis that the degree of tumor cell density represented on certain MRI sequences may be sex-specific.

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 PURPOSE Brain tumor related epilepsy (BTE) is a major co-morbidity in patients with glioma. It is difficult to determine whether the use of anti-epileptic drugs is necessary. We attempted to build a machine-learning model to predict the probability of seizure presentation (SP) with glioma. METHODS We trained a random forest classifier using the following variables: volumetric data of pre-treatment MR images (T1Gd and T2-FLAIR sequences), patient demographics (age; sex), and measurements of tumor proliferation (log(ρ)), invasiveness (log(D)) and their relative ratio (log(ρ/D)). Our cohort consisted of 221 patients total. Using an 80-20 ratio, we used 176 patients (76 SP, 100 nSP) for training and the remaining 45 patients (19 SP, 26 nSP) were used for testing. We also trained on male-only and female-only cohorts to evaluate any sex differences in prediction. For training, 108 males (53 SP, 55 nSP) were used and 28 for testing (14 SP, 14 nSP). We used 72 females (21 SP, 49 nSP) for training and 15 (7 SP, 8 nSP) for testing. We corrected for class imbalance in the female cohort before training. Using 10-fold cross-validation and a separate testing set, we measured performance by ROC curve (AUC), accuracy, sensitivity, and specificity of predictions (average of folds in cross validation). RESULTS The female model achieved the highest AUC (0.853) followed by the mixed model (0.726) and the male model (0.651). In the validation set, the accuracy/sensitivity/specificity of the three cohorts were as follows: mixed (0.726/0.696/0.750), female (0.853/0.830/0.875), and male (0.651/0.577/0.722). The performance of the testing set, in terms of accuracy/sensitivity/specificity were: mixed (0.733/0.74/0.73), female (0.8/0.57/1), and male (0.714/0.64/0.79). CONCLUSION We found a negative correlation between seizure probability and size and invasiveness of tumors. Our model shows promising performance on testing set data. Further cohort studies and training is warranted.

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.