Citations

Abstract High-grade gliomas (HGGs) nearly always recur after standard initial treatment, and the resulting mixture of recurrent tumor and treatment-induced reactive changes presents major diagnostic challenges. Anatomical imaging, such as MRI, cannot adequately distinguish progressive disease from treatment effect (pseudo-progression). Furthermore, there is marked intra-tumoral heterogeneity, such that some areas of a tumor may demonstrate necrotic treatment effect and others frank recurrence. Due to this difficulty reliably differentiating between these two clinical findings, analytic methods using multiple modalities are necessary to further our understanding of this disease process. To this end, we sought to correlate radiographic, histopathologic and molecular features of surgically sampled post-treatment suspected recurrence to identify markers distinguishing tumor growth from treatment effect. We performed Stimulated Raman Histology (SRH) imaging and highly multiplexed RNA-sequencing (PLATE-seq) on 84 MRI-localized biopsies from 39 patients with clinically suspected recurrent HGG. The SRH images were classified as recurrent tumor or gliotic/reactive tissue using a convolutional neural network trained on an independent cohort including a large set of recurrent HGG, and an automated cell-counting algorithm was used to quantify cellularity from the SRH image of each sample. Differential gene expression analysis of the PLATE-seq data was used to identify gene sets that distinguish recurrent tumor from treatment effect, and single sample gene set variation analysis (GSVA) was used to further assess the molecular and cellular composition of each MRI-localized sample. The histopathologic and molecular features of each sample were also correlated with the MRI features of the corresponding biopsy sites, and this data is currently being used to train machine learning models that predict the distribution of recurrent tumor and treatment-induced reactive changes within a patient’s radiographic lesion. These predictive radiomic models will help to guide neurosurgical sampling, and improve our ability to monitor glioma progression and response to therapy.

Abstract Glioblastomas (GBMs) are known to be complex tumors comprising multiple subpopulations of genetically-distinct cancer cells; it is thought that this genetic variation is a major factor in the lack of observed survival benefit of treatment regimes that target one of these subpopulations. The field of radiogenomics seeks to study correlations between MRI patterns and genetic features of GBM tumors. Spatial radiogenomic maps produced using machine-learning (ML) methods that are trained against information from image-localized patient biopsies identify regions where particular cancer sub-populations are predicted to occur within a GBM, thus non-invasively characterizing the regional genetic variability of these tumors. These tumor subpopulations may also interact with one another, in ways which may be of a competitive or cooperative nature to varying degrees. It is important to ascertain the nature of these interactions, as they may have implications for treatment response to targeted therapies, and characterization of the spatio-temporal dynamics of these co-evolving sub-populations will shed light on why some therapies fail. Here we combine mathematical modeling techniques and spatially-resolved radiogenomic maps to study the nature of these interactions between molecularly-distinct GBM subpopulations. We model the interactions between cell populations using a partial differential equation based formalism. The model is parameterized using radiogenomic ML maps from which we infer the nature of interactions between subpopulations. Furthermore, using maps as inputs, the model turns static maps into dynamic information, thus providing insight into how these subpopulations composing the tumor change over time and the effect this has on observed treatment response for individual patients.

Abstract Sex differences in brain tumor related epilepsy (BTE) as an initial symptom of glioma are rarely investigated. Sex differences are being explored with increased interest more recently. Epidemiological studies have generally indicated overall incidence of epilepsy to be slightly higher in males than females, but a higher incidence of idiopathic generalized epilepsy has been seen in females in some studies. Clinicians regularly face the challenge of anti-epileptic drug prescription determination in patients with newly identified brain mass and would greatly benefit from knowing clinical characteristics in BTE. We analyzed clinical characteristics in patients with pathologically-determined glioma with known seizure status at clinical presentation using non-parametric and Fisher’s Exact tests with the entire cohort and male and female separately. From our clinical data repository, we identified 223 patients (females, n=86 (39%); males, n=137 (61%)) divided by those that presented with seizure (SP, n=96, 43%; F, n= 28; M, n= 68) and those that presented without seizure (nSP, n=127, 57%, F n=58, M n=69). In sex-specific analysis males had a significantly (x2=6.28, p=0.01) higher proportion of SP (49.6%) than females (32.6%). There was a dependence between seizure presentation status and tumor grade in males (x2=5.94, p=0.015). Tumor kinetics in the female cohort showed significantly (p< 0.0001) more diffuse tumors in SP compared to nSP. Right-sided tumors were predominantly nSP, while left-sided tumors were split more evenly between SP and nSP. The predominance of nSP among right-sided tumors seems to be driven by females (74%), male right-sided tumors were split more evenly between SP and nSP. Seizure status is dependent on tumor location for the combined cohort, however not for males (x2=5.23, p=0.16) or females (x2=4.34, p=0.23). Clearly there are sex differences and similarities in clinical characteristics at glioma presentation. Further characterization could lead to improved adherence of national standards for AED prescription.

Abstract OBJECTIVES Dogma suggests that for brain tumors, regions of enhancement on T1-weighted gadolinium contrast enhanced (T1Gd) magnetic resonance imaging (MRI) correlate with intravenously delivered drug distribution as enhancement indicates a compromised blood-brain barrier (BBB). However, poor response to intravenous therapies highlights the importance of the diffuse disease beyond enhancing regions. This study investigated whether imaging features can provide an accurate prediction of drug distribution. METHODS Eight brain tumor patients (7 gliomas and 1 metastatic adenocarcinoma) were included in this Phase 0 trial. Presurgery T1-weighted, T1Gd, T2-weighted gadolinium contrast enhanced (T2Gd), and T2-weighted Fluid Attenuated Inversion Recovery (T2FLAIR) MRIs were acquired. All images underwent bias correction using the N4 algorithm, standardization of intensities, and registration. Prior to incision, patients received both an antibiotic cefazolin (6% BBB penetrance) and levetiracetam (80% BBB penetrance), an anti-seizure drug. Subsequently, multiple blood samples and image-guided biopsies were taken and analyzed for drug concentration using liquid chromatography mass spectrometry. Biopsy drug levels are reported as Brain-Plasma Ratio (BPR), the ratio of biopsy concentration relative to plasma concentration. Mean image intensity was extracted from an 8x8 mm window surrounding each biopsy location. Regression analysis was performed to determine which combination of image types were linearly predictive of BPR for both drugs. Correlations were also analyzed according to the biopsy location radiographic appearance. RESULTS Regression analysis revealed that T2Gd intensity was linearly predictive of cefazolin BPR and FLAIR intensity was linearly predictive of levetiracetam BPR (p=0.009 and 0.041, respectively). Grouping samples according the the radiographic appearance revealed that levetiracetam BPR had a similar pattern of values to that of FLAIR intensities and cefazolin BPR had a similar pattern to T2, further supporting the regression analysis results. CONCLUSIONS Local concentrations of drug may be related to T2-weighted signals (T2Gd and T2FLAIR) rather than the gadolinium distribution on T1Gd images.

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). Additionally, extent of resection, tumor laterality, and IDH1 mutation status were also found to have sex-specific effects on overall survival. 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.

Although glioblastoma (GBM) is a fatal primary brain cancer with short median survival of 15 months, a small number of patients survive >5 years after diagnosis; they are known as extreme survivors (ES). Because of their rarity, very little is known about what differentiates these outliers from other patients with GBM. For the purpose of identifying unknown drivers of extreme survivorship in GBM, the ENDURES consortium (ENvironmental Dynamics Underlying Responsive Extreme Survivors of GBM) was developed. This consortium is a multicenter collaborative network of investigators focused on the integration of multiple types of clinical data and the creation of patient-specific models of tumor growth informed by radiographic and histologic parameters. Leveraging our combined resources, the goals of the ENDURES consortium are 2-fold: (1) to build a curated, searchable, multilayered repository housing clinical and outcome data on a large cohort of ES patients with GBM; and (2) to leverage the ENDURES repository for new insights into tumor behavior and novel targets for prolonging survival for all patients with GBM. In this article, the authors review the available literature and discuss what is already known about ES. The authors then describe the creation of their consortium and some preliminary results.

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). Additionally, extent of resection, tumor laterality, and IDH1 mutation status were also found to have sex-specific effects on overall survival. 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.

Tumor recurrence in glioblastoma multiforme (GBM) is often attributed to acquired resistance to the standard chemotherapeutic agent, temozolomide (TMZ). Promoter methylation of the DNA repair gene MGMT (O6-methylguanine–DNA methyltransferase) has been associated with sensitivity to TMZ, whereas increased expression of MGMT has been associated with TMZ resistance. Clinical studies have observed a downward shift in MGMT methylation percentage from primary to recurrent stage tumors; however, the evolutionary processes that drive this shift and more generally the emergence and growth of TMZ-resistant tumor subpopulations are still poorly understood. Here, we develop a mathematical model, parameterized using clinical and experimental data, to investigate the role of MGMT methylation in TMZ resistance during the standard treatment regimen for GBM—surgery, chemotherapy, and radiation. We first found that the observed downward shift in MGMT promoter methylation status between detection and recurrence cannot be explained solely by evolutionary selection. Next, our model suggests that TMZ has an inhibitory effect on maintenance methylation of MGMT after cell division. Finally, incorporating this inhibitory effect, we study the optimal number of TMZ doses per adjuvant cycle for patients with GBM with high and low levels of MGMT methylation at diagnosis.

Purpose Glioblastomas, lethal primary brain tumors, are known for their heterogeneity and invasiveness. A growing body of literature has been developed demonstrating the clinical relevance of a biomathematical model, the proliferation-invasion model, of glioblastoma growth. Of interest here is the development of a treatment response metric, days gained (DG). This metric is based on individual tumor kinetics estimated through segmented volumes of hyperintense regions on T1-weighted gadolinium-enhanced and T2-weighted magnetic resonance images. This metric was shown to be prognostic of time to progression. Furthermore, it was shown to be more prognostic of outcome than standard response metrics. Although promising, the original article did not account for uncertainty in the calculation of the DG metric, leaving the robustness of this cutoff in question. Methods We harnessed the Bayesian framework to consider the impact of two sources of uncertainty: (1) image acquisition and (2) interobserver error in image segmentation. We first used synthetic data to characterize what nonerror variants are influencing the final uncertainty in the DG metric. We then considered the original patient cohort to investigate clinical patterns of uncertainty and to determine how robust this metric is for predicting time to progression and overall survival. Results Our results indicate that the key clinical variants are the time between pretreatment images and the underlying tumor growth kinetics, matching our observations in the clinical cohort. Finally, we demonstrated that for this cohort, there was a continuous range of cutoffs between 94 and 105 for which the prediction of the time to progression was over 80% reliable. Conclusion Although additional validation must be performed, this work represents a key step in ascertaining the clinical utility of this metric.

Male and female glioblastomas are biologically distinct, and improving outcomes may require sex-specific approaches to treatment.