Search

The influence of biological sex differences on human health and disease, while being increasingly recognized, has long been underappreciated and underexplored. While humans of all sexes are more alike than different, there is evidence for sex differences in the most basic aspects of human biology and these differences have consequences for the etiology and pathophysiology of many diseases. In a disease like cancer, these consequences manifest in the sex biases in incidence and outcome of many cancer types. The ability to deliver precise, targeted therapies to complex cancer cases is limited by our current understanding of the underlying sex differences. Gaining a better understanding of the implications and interplay of sex differences in diseases like cancer will thus be informative for clinical practice and biological research. Here we review the evidence for a broad array of biological sex differences in humans and discuss how these differences may relate to observed sex differences in various diseases, including many cancers and specifically glioblastoma. We focus on areas of human biology that play vital roles in healthy and disease states, including metabolism, development, hormones, and the immune system, and emphasize that the intersection of sex differences in these areas should not go overlooked. We further propose that mathematical approaches can be useful for exploring the extent to which sex differences affect disease outcomes and accounting for those in the development of therapeutic strategies.

High-grade glioma (HGG), and particularly Glioblastoma (GBM), can exhibit pronounced intratumoral heterogeneity that confounds clinical diagnosis and management. While conventional contrast-enhanced MRI lacks the capability to resolve this heterogeneity, advanced MRI techniques and PET imaging offer a spectrum of physiologic and biophysical image features to improve the specificity of imaging diagnoses. Published studies have shown how integrating these advanced techniques can help better define histologically distinct targets for surgical and radiation treatment planning, and help evaluate the regional heterogeneity of tumor recurrence and response assessment following standard adjuvant therapy. Application of texture analysis and machine learning (ML) algorithms has also enabled the emerging field of radiogenomics, which can spatially resolve the regional and genetically distinct subpopulations that coexist within a single GBM tumor. This review focuses on the latest advances in neuro-oncologic imaging and their clinical applications for the assessment of intratumoral heterogeneity.

Lacunarity, a quantitative morphological measure of how shapes fill space, and fractal dimension, a morphological measure of the complexity of pixel arrangement, have shown relationships with outcome across a variety of cancers. However, the application of these metrics to glioblastoma (GBM), a very aggressive primary brain tumor, has not been fully explored. In this project, we computed lacunarity and fractal dimension values for GBM-induced abnormalities on clinically standard magnetic resonance imaging (MRI). In our patient cohort (n = 402), we connect these morphological metrics calculated on pretreatment MRI with the survival of patients with GBM. We calculated lacunarity and fractal dimension on necrotic regions (n = 390), all abnormalities present on T1Gd MRI (n = 402), and abnormalities present on T2/FLAIR MRI (n = 257). We also explored the relationship between these metrics and age at diagnosis, as well as abnormality volume. We found statistically significant relationships to outcome for all three imaging regions that we tested, with the shape of T2/FLAIR abnormalities that are typically associated with edema showing the strongest relationship with overall survival. This link between morphological and survival metrics could be driven by underlying biological phenomena, tumor location or microenvironmental factors that should be further explored.

Digital, OBJECTIVE Recent studies have proposed resection of the T2 FLAIR hyperintensity beyond the T1 contrast enhancement (supramarginal resection [SMR]) for IDH–wild-type glioblastoma (GBM) to further improve patients’ overall survival (OS). GBMs have significant variability in tumor cell density, distribution, and infiltration. Advanced mathematical models based on patient-specific radiographic features have provided new insights into GBM growth kinetics on two important parameters of tumor aggressiveness: proliferation rate (ρ) and diffusion rate (D). The aim of this study was to investigate OS of patients with IDH–wild-type GBM who underwent SMR based on a mathematical model of cell distribution and infiltration profile (tumor invasiveness profile). METHODS Volumetric measurements were obtained from the selected regions of interest from pre- and postoperative MRI studies of included patients. The tumor invasiveness profile (proliferation/diffusion [ρ/D] ratio) was calculated using the following formula: ρ/D ratio = (4π/3)2/3 × (6.106/[VT21/1 − VT11/1])2, where VT2 and VT1 are the preoperative FLAIR and contrast-enhancing volumes, respectively. Patients were split into subgroups based on their tumor invasiveness profiles. In this analysis, tumors were classified as nodular, moderately diffuse, or highly diffuse. RESULTS A total of 101 patients were included. Tumors were classified as nodular (n = 34), moderately diffuse (n = 34), and highly diffuse (n = 33). On multivariate analysis, increasing SMR had a significant positive correlation with OS for moderately and highly diffuse tumors (HR 0.99, 95% CI 0.98–0.99; p = 0.02; and HR 0.98, 95% CI 0.96–0.99; p = 0.04, respectively). On threshold analysis, OS benefit was seen with SMR from 10% to 29%, 10% to 59%, and 30% to 90%, for nodular, moderately diffuse, and highly diffuse, respectively. CONCLUSIONS The impact of SMR on OS for patients with IDH–wild-type GBM is influenced by the degree of tumor invasiveness. The authors’ results show that increasing SMR is associated with increased OS in patients with moderate and highly diffuse IDH–wild-type GBMs. When grouping SMR into 10% intervals, this benefit was seen for all tumor subgroups, although for nodular tumors, the maximum beneficial SMR percentage was considerably lower than in moderate and highly diffuse tumors., Ciencias Médicas y de la Salud

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.

Standard-of-care multiparameter magnetic resonance imaging (MRI) scans of the brain were used to objectively subdivide glioblastoma multiforme (GBM) tumors into regions that correspond to variations in blood flow, interstitial edema, and cellular density. We hypothesized that the distribution of these distinct tumor ecological “habitats” at the time of presentation will impact the course of the disease. We retrospectively analyzed initial MRI scans in 2 groups of patients diagnosed with GBM, a long-term survival group comprising subjects who survived &gt, 36 month postdiagnosis, and a short-term survival group comprising subjects who survived ≤19 month postdiagnosis. The single-institution discovery cohort contained 22 subjects in each group, while the multi-institution validation cohort contained 15 subjects per group. MRI voxel intensities were calibrated, and tumor voxels clustered on contrast-enhanced T1-weighted and fluid-attenuated inversion-recovery (FLAIR) images into 6 distinct “habitats” based on low- to medium- to high-contrast enhancement and low–high signal on FLAIR scans. Habitat 6 (high signal on calibrated contrast-enhanced T1-weighted and FLAIR sequences) comprised a significantly higher volume fraction of tumors in the long-term survival group (discovery cohort, 35% ± 6.5%, validation cohort, 34% ± 4.8%) compared with tumors in the short-term survival group (discovery cohort, 17% ± 4.5%, p &lt, 0.03, validation cohort, 16 ± 4.0%, p &lt, 0.007). Of the 6 distinct MRI-defined habitats, the fractional tumor volume of habitat 6 at diagnosis was significantly predictive of long- or short-term survival. We discuss a possible mechanistic basis for this association and implications for habitat-driven adaptive therapy of GBM.

Purpose:Quantify in vivo biomechanical tissue properties in various breast densities and in normal risk and high risk women using Magnetic Resonance Imaging (MRI)/MRE and examine the association between breast biomechanical properties and cancer risk.Methods: Patients with normal risk or high risk of breast cancer underent 3.0 T breast MR imaging and elastography. Breast parenchymal enhancement (BPE), density (from most recent mammogram), stiffness, elasticity, and viscosity were recorded. Within each breast density group (non-dense versus dense), stiffness, elasticity, and viscosity were compared across risk groups (normal versus high). A multivariable logistic regression model was used to evaluate whether the MRE parameters (separately for stiffness, elasticity, and viscosity) predicted risk status after controlling for clinical factors.Results: 50 normal risk and 86 high risk patients were included. Risk groups were similar on age, density, and menopausal status. Among patients with dense breasts, mean stiffness, elasticity, and viscosity were significantly higher in high risk patients (N = 55) compared to normal risk patients (N = 34; all p < 0.001). Stiffness remained a significant predictor of risk status (OR=4.26, 95% CI [1.96, 9.25]) even after controlling for breast density, BPE, age, and menopausal status. Similar results were seen for elasticity and viscosity.Conclusion: A structurally-based, quantitative biomarker of tissue stiffness obtained from MRE is associated with differences in breast cancer risk in dense breasts. Tissue stiffness could provide a novel prognostic marker to help identify high risk women with dense breasts who would benefit from increased surveillance and/or risk reduction measures.

Intra-tumor genetic heterogeneity is an important cause of treatment failure of GBM. Using MRI and image-localized biopsies, it is possible to train machine learning (ML) models to predict regional genetic status. However, biopsy samples are limited, making it difficult to train a robust ML model. We proposed a data-inclusive model called Weakly Supervised Ordinal Support Vector Machines (WSO-SVM) which leverages the vast amount of MRI data outside the sparsely sampled biopsy regions to augment the biopsy samples to improve ML accuracy. Our study included a unique dataset of 104 image-localized biopsies with spatially matched multiparametric MRI from 30 untreated Glioblastoma (GBM) patients. Each biopsy sample went through genetic sequencing analysis and our study focused on two GBM hallmark genes, EGFR and PDGFRA. For each gene, a biopsy sample was labeled as “altered” if the copy number of this gene was amplified or a mutation was found, and “non-altered” otherwise. From the localized region of six MRI contrasts from T1gd, T2w, diffusion and perfusion imaging, over 300 texture features were extracted. To account for biopsy sample location uncertainty, six neighboring regions of the biopsy sample including four neighbors two pixels away from the biopsy location and two neighbors on adjacent slices were also included in model training. WSO-SVM achieved 0.83 accuracy, 0.77 sensitivity, and 0.86 specificity for classifying EGFR; 0.77 accuracy, 0.74 sensitivity, and 0.79 specificity for classifying PDGFRA, based on 10-fold cross validation. Furthermore, using the trained models, we generated regional EGFR and PDGFRA alteration maps for each patient within the enhancing and non-enhancing tumoral areas. On average we found a greater proportion of enhancing tumor with co-alteration than non-enhancing tumor, while this trend was reversed when considering only one altered gene. The ratio of EGFR vs PDGFRA alterations was higher in non-enhancing tumor than enhancing tumor.

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.

Morphological characteristics have been linked to outcomes across a variety of cancers. Lacunarity is a quantitative morphological measure of how shapes fill space while fractal dimension is a morphological measure of the complexity of pixel arrangement. Glioblastoma is the most aggressive primary brain tumor with a short expected survival given the current standard-of-care treatment. Due to the sensitive location of the tumor, there is a heavy reliance on imaging to assess the state of the disease in the clinic. In this project, we computed lacunarity and fractal dimension values for glioblastoma-induced abnormalities on gadolinium-enhanced T1-weighted magnetic resonance imaging (T1Gd MRI) as well as T2-weighted (T2) and fluid-attenuated inversion recovery (FLAIR) MRIs. In our patient cohort (n=402), we aim to connect these morphological metrics calculated on pretreatment MRI with the survival of patients with GBM. We calculated lacunarity and fractal dimension across all MRI slices on necrotic regions (n=390) and abnormalities on T1Gd MRI (n=402), as well as on enhancing abnormalities present on T2/FLAIR MRI (n=257). We also explored the relationship between these metrics and age at diagnosis, as well as abnormality volume. We found statistically significant relationships to outcome across all three imaging subtypes, with the shape of T2/FLAIR abnormalities showing the strongest relationship with overall survival. The link between morphological and survival metrics could be driven by underlying biological phenomena, tumor location or microenvironmental factors that should be further explored.

Gliomas are brain tumors characterized by highly variable growth patterns. Magnetic resonance imaging (MRI) is the cornerstone of glioma diagnosis and management planning. However, glioma features on MRI do not directly correlate with tumor cell distribution. Additionally, there is evidence that glioma tumor characteristics and prognosis are sex-dependent. Magnetic resonance elastography (MRE) is an imaging technique that allows interrogation of tissue stiffness in-vivo and has found utility in the imaging of several cancers. We investigate the relationship between MRI features, MRE features, and growth parameters derived from an established mathematical model of glioma proliferation and invasion. Results suggest that both the relationship between tumor volume and tumor stiffness as well as the relationship between the parameters derived from the mathematical model and tumor stiffness are sex-dependent. These findings lend evidence to a growing body of knowledge about the clinical importance of sex in the context of cancer diagnosis, prognosis and treatment.

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.

Glioblastoma (GBM) is the most aggressive primary brain tumor and can have cystic components, identifiable through magnetic resonance imaging (MRI). Previous studies suggest that cysts occur in 7-23% of GBMs and report mixed results regarding their prognostic impact. Using our retrospective cohort of 493 patients with first-diagnosis GBM, we carried out an exploratory analysis on this potential link between cystic GBM and survival. Using pretreatment MRIs, we manually identified 88 patients with GBM that had a significant cystic component at presentation and 405 patients that did not. Patients with cystic GBM had significantly longer overall survival and were significantly younger at presentation. Within patients who received the current standard of care (SOC) (N=184, 40 cystic), we did not observe a survival benefit of cystic GBM. Unexpectedly, we did not observe a significant survival benefit between this SOC cystic cohort and patients with cystic GBM diagnosed before the standard was established (N=40 with SOC, N=19 without SOC); this significant SOC benefit was clearly observed in patients with noncystic GBM (N=144 with SOC, N=111 without SOC). When stratified by sex, this significant survival benefit was only preserved in male patients (N=303, 47 cystic). We report differences in the absolute and relative sizes of imaging abnormalities on MRI and the prognostic implication of cysts based on sex. We discuss hypotheses for these differences, including the possibility that the presence of a cyst could indicate a less aggressive tumor.

Glioblastomas (GBMs) are the most aggressive primary brain tumours and have no known cure. Each individual tumour comprises multiple sub-populations of genetically-distinct cells that may respond differently to targeted therapies and may contribute to disappointing clinical trial results. Image-localized biopsy techniques allow multiple biopsies to be taken during surgery and provide information that identifies regions where particular sub-populations occur within an individual GBM, thus providing insight into their regional genetic variability. These sub-populations may also interact with one another through a competitive or cooperative nature; it is important to ascertain the nature of these interactions, as they may have implications for responses to targeted therapies. We combine genetic information from biopsies with a mechanistic model of interacting GBM sub-populations to characterise the nature of interactions between two commonly occurring GBM sub-populations, those with EGFR and PDGFRA genes amplified. We study population levels found across image-localized biopsy data from a cohort of 25 patients and compare this to model outputs under competitive, cooperative and neutral interaction assumptions. We explore other factors affecting the observed simulated sub-populations, such as selection advantages and phylogenetic ordering of mutations, which may also contribute to the levels of EGFR and PDGFRA amplified populations observed in biopsy data.

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.

Background. Neuroanatomic locations of gliomas may influence clinical presentations, molecular profiles, and patients’ prognoses. Methods. We investigated our institutional cancer registry to include patients with glioma over a 10-year period. Statistical tests were used to compare demographic, genetic, and clinical characteristics among patients with gliomas in different locations. Survival analysis methods were then used to assess associations between location and overall survival in the full cohort, as well as in relevant subgroups. Results. 182 gliomas were identified. Of the tumours confined to a single lobe, there were 51 frontal (28.0%), 50 temporal (27.5%), 22 parietal (12.1%), and seven occipital tumours (3.8%) identified. Tumours affecting the temporal lobe were associated with reduced overall survival when compared to all other tumours (11 months vs. 13 months, log-rank p = 0.0068). In subgroup analyses, this result was significant for males [HR (95%CI) 2.05 (1.30, 3.24), p = 0.002], but not for females [HR (95%CI) 1.12 (0.65, 1.93), p = 0.691]. Out of 82 cases tested for IDH-1, 10 were mutated (5.5%). IDH-1 mutation was present in six frontal, two temporal, one thalamic, and one multifocal tumour. Out of 21 cases tested for 1p19q deletions, 12 were co-deleted, nine of which were frontal lobe tumours. MGMT methylation was assessed in 45 cases; 7/14 frontal tumours and 6/13 temporal tumours were methylated. Conclusion. Our results support the hypothesis that the anatomical locations of gliomas influence patients’ clinical courses. Temporal lobe tumours were associated with poorer survival, though this association appeared to be driven by these patients’ more aggressive tumour profiles and higher risk baseline demographics. Independently, female patients who had temporal lobe tumours fared better than males. Molecular analysis was limited by the low prevalence of genetic testing in the study sample, highlighting the importance of capturing this information for all gliomas. Importance of this study. The specific neuroanatomic location of tumours in the brain is thought to be predictive of treatment options and overall prognosis. Despite evidence for the clinical significance of this information, there is relatively little information available regarding the incidence and prevalence of tumours in the different anatomical regions of the brain. This study has more fully characterised tumour prevalence in different regions of the brain. Additionally, we have analysed how this information may affect tumours’ molecular characteristics, treatment options offered to patients, and patients’ overall survival. This information will be informative both in the clinical setting and in directing future research.

Glioblastoma (GBM) is the most aggressive primary brain tumor with a short median survival. Tumor recurrence is a clinical expectation of this disease and usually occurs along the resection cavity wall. However, previous clinical observations have suggested that in cases of ischemia following surgery, tumors are more likely to recur distally. Through the use of a previously established mechanistic model of GBM, the Proliferation Invasion Hypoxia Necrosis Angiogenesis (PIHNA) model, we explore the phenotypic drivers of this observed behavior. We have extended the PIHNA model to include a new nutrient-based vascular efficiency term that encodes the ability of local vasculature to provide nutrients to the simulated tumor. The extended model suggests sensitivity to a hypoxic microenvironment and the inherent migration and proliferation rates of the tumor cells are key factors that drive distal recurrence.

BackgroundSex 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.MethodsCombining 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. Pretreatment MR images of 494 glioblastoma patients (299 males and 195 females) were segmented to quantify tumor volumes. Cox proportional hazard (CPH) models and Student’s t-tests were used to assess which variables were associated with survival outcomes.ResultsAmong 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, pConclusionDespite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes, which emphasizes the importance of considering sex as a biological factor when determining patient prognosis and treatment approach.

We analyze the wave speed of the Proliferation Invasion Hypoxia Necrosis Angiogenesis (PIHNA) model that was previously created and applied to simulate the growth and spread of glioblastoma (GBM), a particularly aggressive primary brain tumor. We extend the PIHNA model by allowing for different hypoxic and normoxic cell migration rates and study the impact of these differences on the wave-speed dynamics. Through this analysis, we find key variables that drive the outward growth of the simulated GBM. We find a minimum tumor wave-speed for the model; this depends on the migration and proliferation rates of the normoxic cells and is achieved under certain conditions on the migration rates of the normoxic and hypoxic cells. If the hypoxic cell migration rate is greater than the normoxic cell migration rate above a threshold, the wave speed increases above the predicted minimum. This increase in wave speed is explored through an eigenvalue and eigenvector analysis of the linearized PIHNA model, which yields an expression for this threshold. The PIHNA model suggests that an inherently faster-diffusing hypoxic cell population can drive the outward growth of a GBM as a whole, and that this effect is more prominent for faster-proliferating tumors that recover relatively slowly from a hypoxic phenotype. The findings presented here act as a first step in enabling patient-specific calibration of the PIHNA model.

Non-invasive magnetic resonance imaging (MRI) is the primary imaging modality for visualizing brain tumor growth and treatment response. While standard MRIs are central to clinical decision making, advanced quantitative imaging sequences like diffusion weighted imaging (DWI) are increasingly relied on. Deciding the best way to interpret DWIs, particularly in the context of treatment, is still an area of intense research. With DWI being indicative of tissue structure, it is important to establish the link between DWI and brain tumor mathematical growth models, which could help researchers and clinicians better understand the tumor’s microenvironmental landscape. Our goal was to demonstrate the potential for creating a DWI patient-specific untreated virtual imaging control (UVICs), which represents an individual tumor’s untreated growth and could be compared with actual patient DWIs. We generated a DWI UVIC by combining a patient-specific mathematical model of tumor growth with a multi-compartmental MRI signal equation. GBM growth was mathematically modeled using the Proliferation-Invasion-Hypoxia-Necrosis-Angiogenesis-Edema (PIHNA-E) model, which simulated tumor as being comprised of multiple cellular phenotypes interacting with vasculature, angiogenic factors, and extracellular fluid. The model’s output consisted of spatial volume fraction maps for each microenvironmental species. The volume fraction maps and corresponding T2 and apparent diffusion coefficient (ADC) values from literature were incorporated into a multi-compartmental signal equation to simulate DWI images. Simulated DWIs were created at multiple b-values and then used to calculate ADC maps. We found that the regional ADC values of simulated tumors were comparable to literature values.

Summary We present MetaMarker, a pipeline for discovering metagenomic biomarkers from whole-metagenome sequencing samples. Different from existing methods, MetaMarker is based on a de novo approach that does not require mapping raw reads to a reference database. We applied MetaMarker on whole-metagenome sequencing of colorectal cancer (CRC) stool samples from France to discover CRC specific metagenomic biomarkers. We showed robustness of the discovered biomarkers by validating in independent samples from Hong Kong, Austria, Germany and Denmark. We further demonstrated these biomarkers could be used to build a machine learning classifier for CRC prediction. Availability and implementation MetaMarker is freely available at https://bitbucket.org/mkoohim/metamarker under GPLv3 license. Supplementary information Supplementary data are available at Bioinformatics online.

Introduction Distinguishing between tumor and treatment effect in post-treatment glioma, although crucial for clinical management, is difficult because contrast-enhancing regions are mixtures of recurrent tumor and reactive tissue, and definitive histopathological criteria do not exist. This study disentangles the marked intra-tumoral heterogeneity in the treatment-recurrent setting by developing an unsupervised framework to algorithmically categorize intraoperative MRI-localized biopsies into three clinically-relevant tissue clusters based on joint analysis of RNA sequencing and histopathological data. Methods A retrospective cohort of 84 MRI-localized biopsies from 37 patients with post-treatment recurrent glioblastoma underwent mRNA extraction and quantification via PLATEseq protocol. For 48 of 84 biopsies, a neighboring piece of tissue underwent quantitative histopathology based on labeling index (LI) for SOX2, CD68, NeuN, Ki67, and H&E. Correlation between LIs and gene expression for these 48 samples was performed. Genes significantly correlated (p Results Correlation analysis identified 7779 genes significantly correlated with ≥1 histopathological marker. Clustering revealed three gene sets associated with specific markers: SetA-3688 genes associated with SOX2/Ki67/H&E; SetB-2418 genes associated with CD68; SetC-1673 genes associated with NeuN. ssGSEA using these sets categorized each biopsy into one of three tissue types: 27 biopsies enriched in SetA, 28 in SetB, and 29 in SetC. Conclusions Using MRI-localized biopsies with both RNAseq and histopathological data, this algorithmic approach allows development of three orthogonal tissue-specific gene sets that can be applied to characterize the heterogeneity in post-treatment recurrent glioma: SetA: genes correlated with SOX2/Ki67/H&E, representing recurrent tumor; SetB: genes correlated with CD68 (microglia) representing reactive tissue consistent with treatment effect; SetC: genes correlated with NeuN (neurons), representing infiltrated brain.

The blood–brain barrier (BBB) excludes the vast majority of cancer therapeutics from normal brain. However, the importance of the BBB in limiting drug delivery and efficacy is controversial in high-grade brain tumors, such as glioblastoma (GBM). The accumulation of normally brain impenetrant radiographic contrast material in essentially all GBM has popularized a belief that the BBB is uniformly disrupted in all GBM patients so that consideration of drug distribution across the BBB is not relevant in designing therapies for GBM. However, contrary to this view, overwhelming clinical evidence demonstrates that there is also a clinically significant tumor burden with an intact BBB in all GBM, and there is little doubt that drugs with poor BBB permeability do not provide therapeutically effective drug exposures to this fraction of tumor cells. This review provides an overview of the clinical literature to support a central hypothesis: that all GBM patients have tumor regions with an intact BBB, and cure for GBM will only be possible if these regions of tumor are adequately treated.

Gliomas are the most common of all primary brain tumors. They are characterized by their diffuse infiltration of the brain tissue and are uniformly fatal, with glioblastoma being the most aggressive form of the disease. In recent years, the over-expression of platelet-derived growth factor (PDGF) has been shown to produce tumors in experimental rodent models that closely resemble this human disease, specifically the proneural subtype of glioblastoma. We have previously modeled this system, focusing on the key attribute of these experimental tumors-the "recruitment" of oligodendroglial progenitor cells (OPCs) to participate in tumor formation by PDGF-expressing retrovirally transduced cells-in one dimension, with spherical symmetry. However, it has been observed that these recruitable progenitor cells are not uniformly distributed throughout the brain and that tumor cells migrate at different rates depending on the material properties in different regions of the brain. Here we model the differential diffusion of PDGF-expressing and recruited cell populations via a system of partial differential equations with spatially variable diffusion coefficients and solve the equations in two spatial dimensions on a mouse brain atlas using a flux-differencing numerical approach. Simulations of our in silico model demonstrate qualitative agreement with the observed tumor distribution in the experimental animal system. Additionally, we show that while there are higher concentrations of OPCs in white matter, the level of recruitment of these plays little role in the appearance of "white matter disease," where the tumor shows a preponderance for white matter. Instead, simulations show that this is largely driven by the ratio of the diffusion rate in white matter as compared to gray. However, this ratio has less effect on the speed of tumor growth than does the degree of OPC recruitment in the tumor. It was observed that tumor simulations with greater degrees of recruitment grow faster and develop more nodular tumors than if there is no recruitment at all, similar to our prior results from implementing our model in one dimension. Combined, these results show that recruitment remains an important consideration in understanding and slowing glioma growth.

Malignant large artery stroke is associated with high mortality of 70% to 80% with best medical management. Decompressive craniectomy (DC) is a highly effective tool in reducing mortality. Convincing evidence has accumulated from several randomized trials, in addition to multiple retrospective studies, that demonstrate not only survival benefit but also improved functional outcome with DC in appropriately selected patients. This article explores in detail the evidence for DC, nuances regarding patient selection, and applicability of DC for supratentorial intracerebral hemorrhage and posterior fossa ischemic and hemorrhagic stroke.

INTRODUCTION Glioblastoma (GBM) is an aggressive primary brain tumor, known for its poor prognosis. Due to its diffuse invasiveness into normal-appearing brain, localized treatments such as surgical resection and radiotherapy are typically supplemented with chemotherapy. However, to reach invading tumor cells, such antineoplastic drugs must cross the blood-brain barrier (BBB). That is, while angiogenesis induces BBB breakdown in dense tumor regions, the BBB remains rather intact for invading GBM cells. As a result, it is unclear whether BBB-impermeable drugs are delivered at a sufficient level to be effective. METHODS In order to study heterogeneity in BBB breakdown, experiments were conducted using both flank and intracranial patient-derived xenografts (PDXs) treated with the EGFR-targeted monoclonal antibody drug conjugate, depatuxizumab mafodotin (ABT-414). Time-series bioluminescence imaging (BLI) data was used to develop a differential equation model of tumor growth for three PDX cell lines. Data from untreated PDXs, both in flank and intracranially, were used to parameterize tumor proliferation rates. Flank PDX data were used to parameterize individual sensitivity to ABT-414, whereas intracranial PDX data were used to determine the proportion of drug exposed to the tumor. RESULTS Each PDX line differed in response to the study drug ABT-414. As expected, such heterogeneous responses can primarily be attributed to differences in both drug sensitivity and the proportion of drug that reached the tumor. Notably, the estimated proportion of drug that reached the tumor was highest in the PDX line with the longest survival times, despite also having higher estimates of resistance. This suggests that PDXs with greater overall BBB breakdown may respond better to this agent. CONCLUSIONS Although more cell lines are needed to validate our approach, parameterizing this model for PDXs gives valuable insight into the extent of BBB breakdown in patient GBMs and may aid in determining optimal therapies for individual patients.

Glioblastoma (GBM) is the most aggressive primary brain tumor with a median overall survival of 15 months with standard-of-care treatment. GBM patients sometimes present with a cystic component, which can be identified through magnetic resonance imaging (MRI). Previous studies suggest that cysts occur in 7–22% of GBM patients and have reported mixed results regarding whether cystic GBM have a survival benefit compared to noncystic GBM. Using our large retrospective cohort of 493 first-diagnosis GBM patients, we aim to elucidate this link between cystic GBM and survival. Within this cohort, 88 patients had a significant cystic component at presentation as identified on MRI. Compared to noncystic GBM (n=405), cystic GBM patients had significantly better overall survival (15 vs 22 months median, log-rank, p=0.001) and were significantly younger at the time of presentation (t-test, p=0.002). However, within patients that received current standard-of-care treatment (n=184), cystic GBM (n=40) was not as beneficial for outcome (22 vs 25 months, log-rank, p=0.3). We also did not observe a significant survival benefit when comparing this standard-of-care cystic cohort to cystic GBM patients diagnosed before the standard was established (n=19, 25 vs 23 months, log-rank, p=0.3), but the analogous result for noncystic GBM patients gives a sizeable benefit, as expected (n=144, n=111, respectively, 22 vs 12 months, log-rank p < 0.0001). Together, these results on current standard-of-care may explain later studies that note no significant survival benefit for cystic GBM patients receiving current standard-of-care. We also report differences in the absolute and relative sizes of imaging abnormalities on MRI and in prognostic impact of cysts based on sex. We discuss current hypotheses for these observed differences, including the possibility that the presence of a cyst could be indicative of a less aggressive tumor.

INTRODUCTION Glioblastoma (GBM) is the most common malignant primary brain tumor in adults with a median overall survival (OS) of 15months. Despite advancements in treatments, prognosis is dismal and the prognostic significance of tumor location is not entirely understood. METHODOLOGY: In our study, we investigated sex-specific volumetric, tumor growth kinetics, and outcome differences among GBMs in various brain locations. Primary GBM patients with pretreatment magnetic resonance imaging (MRI) data (N=289, 173 males, 116 females) were selected from our brain tumor repository. Tumor abnormality was segmented on T1-weighted post-gadolinium contrast agent (T1Gd) MRIs. We utilized the Harvard-Oxford brain atlases to determine the location of GBMs. RESULTS Overall, our study found smaller tumors in the left hemisphere. This may be expected as left-hemispheric GBM symptoms could present earlier, leading to earlier diagnosis and treatment. However, when the cohort was split by sex, we found this observation significant for females only in the parietal lobe (p < 0.0001). Further, female GBMs demonstrated smaller necrotic volume in the left hemisphere (p = 0.030). Sex-specific differences in incidence were noted in the temporal and occipital lobes (2M:1F). Comparing tumor growth kinetics in different brain locations and hemispheres, females had significantly lower tumor proliferation rates in the left hemisphere (p = 0.009) and lower tumor proliferation rates in the left frontal lobe (p = 0.031). Controlling for treatment, patients with frontal lobe tumors had significantly longer OS compared to those with GBMs in the temporal lobe (p = 0.046, 312 days). Differences in growth velocities were noted between frontal and parietal lobe with frontal GBMs having lower velocities in comparison to parietal lobe GBMs. CONCLUSION Together, our results demonstrate that tumor growth and proliferation rates may vary based on location and sex. Additional research is needed to further explore the clinical significance of tumor location.

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.

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.

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.

Glioblastoma (GBM) is the most aggressive primary brain tumor with an average survival of 15 months with standard of care treatment. GBM patients typically present with necrosis surrounded by enhancement on T1-weighted post gadolinium magnetic resonance imaging (T1gd MRI), however some patients present with a significant cystic component. Cysts are caused by different underlying biological mechanisms to necrosis and are important to identify for future clinical investigations. These cystic components can be manually identified through MRI but this process can be time consuming for large patient cohorts. Over the last two decades, our lab has collected serial MRI data of brain tumor patients. With over 70,000 images now in the database and that number increasing daily, it is clear that we have a unique resource for clinical investigation and a need to automate this process. To this end, the aim of this work was to develop and assess the performance of a convolution neural network (CNN) model for automatic detection of cystic GBMs. In this retrospective IRB-approved work, we collected pretreatment MRIs of a patient cohort consisting of 85 patients with a significant cystic component at presentation along with 400 non-cystic GBM, both identified manually through MRI. Image slices with a view of the cystic component were used as positive samples for training. Data were randomly split into training, validation, and test sets using a 70:15:15 ratio. The proportion of positive to negative cases was comparable between sets. Prior to training, we used image augmentation techniques to compensate for the class imbalance in our data. Our results showed that deep learning networks can automatically detect cystic GBMs on MRIs with high accuracy and thus illustrates the potential use of this technique in clinically relevant settings.

Glioblastoma (GBM) is the most aggressive primary brain tumor with a median survival of only 15 months with standard of care treatment. Lacunarity, a quantitative morphological measure of how shapes fill space, and fractal dimension, another morphological measure of the complexity of pixel arrangement, of segmented necrotic regions on gadolinium-enhanced T1 weighted (T1gd) MRI have previously been shown to distinguish both overall survival (OS) and progression free survival (PFS) in GBM (n = 95). In our larger patient cohort (n = 389), we sought to validate or refute previously published results connecting morphological metrics and patient survival. We identified pretreatment necrotic regions of our retrospective first-diagnosis GBM patient cohort using segmented T1gd MRI enhancing regions. We calculated lacunarity and fractal dimension across all T1gd MRI slices with enhancing tumor, and used the median lacunarity and fractal dimension values for our analysis. We find that a lacunarity threshold can significantly distinguish OS (14 months vs 19 months median, log-rank p = 0.015, n = 389) and a fractal dimension threshold can significantly distinguish PFS (8 months vs 11 months median, log-rank p = 0.015, n = 123). We believe that morphological metrics such as lacunarity and fractal dimension could play a role in standard-of-care prognostic considerations at tumor presentation. This link between morphological and survival metrics could be driven by underlying biological phenomena or microenvironmental factors that should be further explored.

There is a need for new approaches and endpoints in oncology drug development, particularly with the advent of immunotherapies and the multiple drug combinations under investigation. Tumor dynamics modeling, a key component to oncology “model-informed drug development,” has shown a growing number of applications and a broader adoption by drug developers and regulatory agencies in the past years to support drug development and approval in a variety of ways. Tumor dynamics modeling is also being investigated in personalized cancer therapy approaches. These models and applications are reviewed and discussed, as well as the limitations and issues open for further investigations. A close collaboration between stakeholders like clinical investigators, statisticians, and pharmacometricians is warranted to advance clinical cancer therapeutics.

Seizures are common presenting symptoms of primary brain tumors. Mechanisms of epileptogenesis are still unknown and are believed to be multifactorial. Previous studies have indicated correlation of seizure with tumor location. Recent investigations of our group have shown image-based parameters have sex-specific implications for patient outcome. In this retrospective study, we examined the association of tumor location with the probability and risk of seizure in male and female glioma patients.

Although glioblastoma is a fatal primary brain cancer with a short median survival of 15 months, a small number of patients survive more than 5 years after diagnosis; they are known as extreme survivors (ES). Due to their rarity, very little is known about what differentiates these outliers from other glioblastoma patients. For the purpose of identifying unknown drivers of extreme survivorship in glioblastoma, we developed the ENDURES consortium (ENvironmental Dynamics Underlying Responsive Extreme Survivors of glioblastoma). 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 histological parameters. Leveraging our combined resources, the goals of the ENDURES consortium are two-fold: (1) to build a curated, searchable, multilayered repository housing clinical and outcome data on a large cohort of ES patients with glioblastoma and (2) to leverage the ENDURES repository for new insights on tumor behavior and novel targets for prolonging survival for all glioblastoma patients. In this article, we review the available literature and discuss what is already known about ES. We then describe the creation of our consortium and some of our preliminary results.FundingThis review was financially supported by a grant from the James S. McDonnell FoundationConflicts of InterestThe authors have declared that no conflicts of interest exist.AuthorshipConceptualized consortium: LW, RG, KME, PC, and KRS. Built consortium: SKJ, PK, NR, JS, KME, PC, and KRS. Wrote the manuscript: SKJ, PW, SCM, PK, AP, and KME. Reviewed and edited the manuscript: LFGC, MMM, AHD, PRJ, and LSH. Contributed to writing, provided feedback, and approved of final manuscript: All authors.Link to website for ENDUREShttp://mathematicalneurooncology.org/?page_id=2125

BackgroundBrain tumor related epilepsy (BTE) is a major co-morbidity related to the management of patients with brain cancer. Despite published practice guidelines recommending against anti-epileptic drug (AED) utilization in patients with gliomas, there is heterogeneity in prescription practices of AEDs in these patients. In an attempt to impact BTE management, we statistically analyzed clinically relevant attributes (sex, age, tumor size, tumor growth kinetics, and tumor location) pertaining to seizure at presentation and used them to build a computational machine learning model to predict the probability of a seizure (at presentation).MethodsFrom our clinical data repository, we identified 223 patients (females, n=86; males, n=137) with pathologically-determined glioma and known seizure status at clinical presentation. Non-parametric and Fisher’s Exact tests were used to identify statistical differences in clinical characteristics. We utilized a random forest machine learning method for generating our predictive models by entire cohort and separated by male and female.FindingsPatients were divided into 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). Females presented with seizures significantly less often than males (x2=6·28, p=0·01). SP patients had significantly smaller T1Gd radius compared to nSP (SP 11·30mm, nSP 18.66mm, pInterpretationDespite heterogeneity across our patient cohort, there is strong evidence of a role for patient sex, tumor size, tumor invasion, and patient age in predicting the incidence of seizures at diagnosis. Future studies, with prospectively detailed data collection, may provide clearer insights into the incidence of seizures through a patient’s treatment course.

A typical feature of glioblastoma (GBM) growth is local recurrence after surgery. However, some GBMs recur distally. It has been noted that GBM patients with perioperative ischemia are more likely to have distal recurrence and that GBM cells migrate faster under hypoxic conditions. We apply the Proliferation Invasion Hypoxia Necrosis Angiogenesis (PIHNA) model to examine the effect of faster hypoxic cell migration on simulated GBM growth. Results suggest that a highly migratory hypoxic cell population drives the growth of the whole tumor and leads to distant recurrence, as do higher normoxic tumor cell migration and low cellular proliferation rates.

Glioblastoma (GBM) is a heterogeneous and lethal brain cancer. These tumors are followed using magnetic resonance imaging (MRI), which is unable to precisely identify tumor cell invasion, impairing effective surgery and radiation planning. We present a novel hybrid model, based on multiparametric intensities, which combines machine learning (ML) with a mechanistic model of tumor growth to provide spatially resolved tumor cell density predictions. The ML component is an imaging data-driven graph-based semi-supervised learning model and we use the Proliferation-Invasion (PI) mechanistic tumor growth model. We thus refer to the hybrid model as the ML-PI model. The hybrid model was trained using 82 image-localized biopsies from 18 primary GBM patients with pre-operative MRI using a leave-one-patient-out cross validation framework. A Relief algorithm was developed to quantify relative contributions from the data sources. The ML-PI model statistically significantly outperformed (p

Glioblastomas are aggressive primary brain tumors known for their inter- and intratumor heterogeneity. This disease is uniformly fatal, with intratumor heterogeneity the major reason for treatment failure and recurrence. Just like the nature vs nurture debate, heterogeneity can arise from intrinsic or environmental influences. Whilst it is impossible to clinically separate observed behavior of cells from their environmental context, using a mathematical framework combined with multiscale data gives us insight into the relative roles of variation from different sources. To better understand the implications of intratumor heterogeneity on therapeutic outcomes, we created a hybrid agent-based mathematical model that captures both the overall tumor kinetics and the individual cellular behavior. We track single cells as agents, cell density on a coarser scale, and growth factor diffusion and dynamics on a finer scale over time and space. Our model parameters were fit utilizing serial MRI imaging and cell tracking data from ex vivo tissue slices acquired from a growth-factor driven glioblastoma murine model. When fitting our model to serial imaging only, there was a spectrum of equally-good parameter fits corresponding to a wide range of phenotypic behaviors. When fitting our model using imaging and cell scale data, we determined that environmental heterogeneity alone is insufficient to match the single cell data, and intrinsic heterogeneity is required to fully capture the migration behavior. The wide spectrum of in silico tumors also had a wide variety of responses to an application of an anti-proliferative treatment. Recurrent tumors were generally less proliferative than pre-treatment tumors as measured via the model simulations and validated from human GBM patient histology. Further, we found that all tumors continued to grow with an anti-migratory treatment alone, but the anti-proliferative/anti-migratory combination generally showed improvement over an anti-proliferative treatment alone. Together our results emphasize the need to better understand the underlying phenotypes and tumor heterogeneity present in a tumor when designing therapeutic regimens., Author summary Glioblastoma, the most common primary brain tumor, is an aggressive and difficult to treat cancer. A key reason is that the tumors can be very heterogeneous, consisting of many different mutants driving distinct cell behaviors. We believe that treatments for this disease could be significantly improved by understanding and quantifying the functional impact of heterogeneity within the tumor. From a clinical standpoint, the larger tissue-scale dynamics, like growth rate, can be informed from serial MRI imaging, while the cell-scale heterogeneity, can be informed by analysis of biopsies. In this work, we combined information from both scales using a mathematical framework and multiscale data from an animal model of glioblastoma. We found that a wide range of potential tumor compositions matched imaging data alone, as a result the model predicts a wide variety of responses to treatment. Using both imaging and cell-scale data narrowed the range of possible tumor compositions and better predicted responses to treatment. The mathematical model also predicted that while targeting migration alone did not slow tumor growth (in fact it drove a more proliferative tumor), an anti-proliferative/anti-migratory treatment combination improved treatment response.

Poor clinical trial outcomes for glioblastoma (GBM) can be attributed to multiple possible causes. GBM is heterogeneous, such that there is a chance of treatment–resistant cells coming to predominate the tumor, and due to the blood brain barrier (BBB) it is also possible that therapy was inadequately delivered to the tumor. Mathematically modeling the dynamics of therapeutic response in patient–derived xenografts (PDX) and fitting the mathematical model to bioluminescence imaging flux data, we may be able to assess the degree to which both drug resistance and drug penetrance are driving varied responses to these therapies.

BACKGROUND: MRI-based modeling of tumor cell density (TCD) can significantly improve targeted treatment of Glioblastoma (GBM). Unfortunately, interpatient variability limits the predictive ability of many modeling approaches. We present a Transfer Learning (TL) method that generates individualized patient models, grounded in the wealth of population data, while also detecting and adjusting for interpatient variabilities based on each patient’s own histologic data. METHODS: We recruited primary GBM patients undergoing image-guided biopsies and preoperative imaging including contrast-enhanced MRI (CE-MRI), Dynamic-Susceptibility-Contrast (DSC)-MRI, and Diffusion Tensor Imaging (DTI). We calculated relative cerebral blood volume (rCBV) from DSC-MRI and mean diffusivity (MD) and fractional anisotropy (FA) from DTI. Following image coregistration, we assessed TCD for each biopsy and identified corresponding localized MRI measurements. We then explored a range of univariate and multivariate predictive models of TCD based on MRI measurements in a generalized one-model-fits-all (OMFA) approach. We then implemented both univariate and multivariate individualized TL predictive models, which harness the available population level data but allow for individual variability in their predictions. Finally, we compared Pearson correlation coefficients and mean absolute error between the individualized TL versus generalized OMFA models. RESULTS: TCD significantly correlated with rCBV (r=0.33,p

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.

Purpose Despite the intra- and intertumoral heterogeneity seen in glioblastoma multiforme (GBM), there is little definitive data on the underlying cause of the differences in patient survivals. Serial imaging assessment of tumor growth allows quantification of tumor growth kinetics (TGK) measured in terms of changes in the velocity of radial expansion seen on imaging. Because a systematic study of this entire TGK phenotype—growth before treatment and during each treatment to recurrence —has never been coordinately studied in GBMs, we sought to identify whether patients cluster into discrete groups on the basis of their TGK. Patients and Methods From our multi-institutional database, we identified 48 patients who underwent maximally safe resection followed by radiotherapy with imaging follow-up through the time of recurrence. The patients were then clustered into two groups through a k-means algorithm taking as input only the TGK before and during treatment. Results There was a significant survival difference between the clusters ( P = .003). Paradoxically, patients among the long-lived cluster had significantly larger tumors at diagnosis ( P = .027) and faster growth before treatment ( P = .003) but demonstrated a better response to adjuvant chemotherapy ( P = .048). A predictive model was built to identify which cluster patients would likely fall into on the basis of information that would be available to clinicians immediately after radiotherapy (accuracy, 90.3%). Conclusion Dichotomizing the heterogeneity of GBMs into two populations—one faster growing yet more responsive with increased survival and one slower growing yet less responsive with shorter survival—suggests that many patients who receive standard-of-care treatments may get better benefit from select alternative treatments.

Temozolomide (TMZ) has been the standard-of-care chemotherapy for glioblastoma (GBM) patients for more than a decade. Despite this long time in use, significant questions remain regarding how best to optimize TMZ therapy for individual patients. Understanding the relationship between TMZ response and factors such as number of adjuvant TMZ cycles, patient age, patient sex, and image-based tumor features, might help predict which GBM patients would benefit most from TMZ, particularly for those whose tumors are not MGMT methylated. Using a cohort of 90 newly-diagnosed GBM patients treated according to the Stupp protocol, we examined the relationships between several patient and tumor characteristics and volumetric and survival outcomes during adjuvant chemotherapy. Volumetric changes in MR imaging abnormalities during adjuvant therapy were used to assess TMZ response. T1Gd volumetric response is associated with younger patient age, increased number of TMZ cycles, longer time to nadir volume, and decreased tumor invasiveness. Moreover, increased adjuvant TMZ cycles corresponded with improved volumetric response only among more nodular tumors, and this volumetric response was associated with improved survival outcomes. Finally, in a subcohort of patients with known MGMT methylation status, MGMT methylated tumors were more diffusely invasive than unmethylated tumors, suggesting that the improved response in nodular tumors is not driven by a preponderance of MGMT methylated tumors. Our finding that less diffusely invasive tumors are associated with greater volumetric response to TMZ suggests that patients with these tumors may benefit from additional cycles of adjuvant TMZ, even for those without MGMT methylation.

Glioblastoma (GBM) is the one of the most aggressive and common primary brain malignancies, with a poor median overall survival of less than 15 months. While the immune system is activated and brain-resident microglia and blood-derived macrophages combat the tumor, the tumor can convert some microglia and macrophages to instead exhibit an immune-suppressive phenotype. These co-opted immune cells are thereby termed ‘glioma-associated microglia and macrophages’ (GAMMs), as they allow for continued tumor growth. However, limited clinical data has been collected to understand this phenomenon. As a result, we have collected spatially-distributed image-localized biopsies from a cohort of patients and performed RNA sequencing on each sample. Correlations between normalized RNA counts of genetic markers for macrophages (i.e., CD68, CD163), tumor populations (i.e., SOX2, OLIG2), and key cell functions (i.e., KI67, CASP3) were analyzed. To further investigate the temporal effects that GAMMs have on GBM growth, we proposed the Proliferation-Invasion-Macrophage (PIM) model. This system of partial differential equations incorporates the proliferative and invasive behavior of GBM, as well as populations for both ‘healthy’ and ‘glioma-associated’ macrophages. By exploring the parameter space, we classified the various dynamics of tumor progression and how they relate to the immune response. With further insights of the interactions between GBM and macrophage populations, we can begin to parameterize the model on a patient-specific basis and provide insights to personalized immunotherapies and other novel immune-targeted treatments.

INTRODUCTION: Brain tumors are difficult to segment in MRI scans. Consequently, we have developed a new system for completely automatic brain tumor segmentation by combining a state-of-the-art 3D deep convolutional neural network (CNN) with a large collection of curated segmentations of brain tumors. METHODS: Our brain tumor database holds 74,722 MRI series from 2,742 unique patients. Over the last 15 years our image analysis team has segmented brain tumors in 35,710 of these series. This preliminary experiment identified 741 pre-treatment studies that included a T1GD and FLAIR scan, and at least one adjudicated brain tumor segmentation. These studies were randomly assigned into 600 training, 41 validation, and 100 test cases. CNN training was performed in two stages: 1) 50 epochs on minimally modified MRI volumes; and, 2) 24 epochs to tune the CNN on skull-stripped volumes. Whole-tumor Dice coefficients (1=perfect overlap, 0=no overlap) were calculated by comparing CNN segmentations against adjudicated segmentations from trained measurers. Training was performed in-the-cloud using an Amazon Machine Instance equipped with an NVidia Tesla V100 GPU, 8 Intel Xeon processors, and 64 GB of RAM. RESULTS: Training required 74 hours. Afterwards, our network required 800 seconds to segment 100 studies in the test set (8 seconds/study). The mean whole-tumor Dice coefficient on the test studies was 0.885. DISCUSSION: The best result on the highly cited 2017 BraTS brain tumor segmentation challenge was a whole-tumor Dice of 0.886, achieved by an ensemble of 7 CNNs. BraTS included 274 studies, each with T1GD, FLAIR, T1 and T2 contrasts. The performance of our single CNN may be due to our comparatively large training set. Our goal is to train our CNN on all series in our database. This may provide useful tools to monitor each patients journey, from diagnosis through treatment.

Clinical neuro-oncology relies on the hyperintensity of gadolinium (Gd) contrast agent on magnetic resonance imaging (MRI) in tumor regions to confirm that the blood-brain barrier (BBB) is locally compromised. While the extent of Gd hyperintensity may indicate that systemically administered drug is being distributed to the tumor regions, little is known about how a drug is distributed and how it may relate to the Gd hyperintensity. Additionally, glioblastomas (GBMs) are diffusely invading neoplasms with a significant fraction of the overall tumor cells spread peripheral to the Gd abnormality, which raises uncertainty as to how or if the rest of the diffuse tumor is affected by drug. Given the gap in understanding drug delivery to the brain, we propose a quantitative approach to model drug delivery in GBM based on MRI and matrix-assisted laser desorption/ionization mass spectroscopy imaging (MALDI). T2-weighted (T2) and T1-weighted with Gd contrast (T1Gd) MRI images were acquired for an animal with the GBM12 orthotopic GBM patient derived xenograft (PDX) line dosed with 100mg/kg erlotinib. A T1Gd region of interest (ROI) captured the Gd-associated hyperintensity. MALDI was performed and aligned with MRI images. A Drug ROI was created to represent the increased intensity of erlotinib on MALDI images. Since the Drug ROI encompassed the T1Gd ROI, we subtracted the two ROIs to create a ‘Drug No T1Gd’ ROI, which represented the drug region beyond the edge of the T1Gd ROI. A ‘Brain ROI’ was created to represent the region of the brain outside of the drug’s spread. The MALDI intensities for the three ROIs were all significantly different (p

The hallmark of glioblastoma (GBM) is poor survivorship. However, a small subset (5%) of patients live greater than 5 years (extreme survivors (ES)). The pathological and tumor determinants of ES are unknown. Giant cell (GC) morphology occurs in < 1% of all GBM, is typically IDH1 wildtype, shows a high frequency of p53 mutations, and is reportedly associated with a somewhat better survival than other IDH1-wildtype GBM. However, the clear association between ES and GC GBM has not be established. We aimed to describe the characteristics of ES tumors that presented with GC GBM features and examine sex differences. In our retrospective multi-institutional database, we identified 90 ES patients with GBM. We reviewed neuropathological reports for the diagnosis of GC GBM or pathology description that included the descriptor(s) of GC: monstrocellular or giant cell or bizarre multinucleated and evaluated phenotypical features in order to describe ES with GC GBM. Values are presented as means. Sixteen (17.8%) ES patients were characterized as GC GBM (males (n=9, 56.25%), females (n=7, 43.75%)). Males were significantly younger than females (39.67 vs 56.29, p=0.018). Females presented with significantly smaller tumors than males (F: 9.49 mm, M: 21.14 mm, p=0.0008). Males less often presented with seizures (33.33%) compared to females (42.86%, p=0.7686). Calculations of invasion/proliferation ratios from our mathematical model revealed no statistically significant differences between females and males. Overall survival between males and females was not statistically significant. Across the board, GC GBM occurred most often in frontal lobes (All: 56% (n=9), 66.67% of males (n=6), and 42.86% of females (n=3)). Extreme survivors of GBM exhibit giant cell features at greater than five-fold incidence than the general GBM population. Sex differences warrant significant attention in future explorations of pathological and tumor characteristics in GC GBM.

PURPOSE: Seizures are frequent symptoms of gliomas. Predicting which patients are more likely to seize and in turn require anti-epileptic management has been a challenge. Correctly identifying these patients could help optimize care and minimize side effects. To provide guidance for clinicians, we used machine learning techniques to predict seizure presentation in this population. METHODS: We used 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)). We compared the performance of 5 machine learning models in predicting seizure status, using Artificial Neural Network, Naive Bayes (NB), Linear Discriminant Analysis (LDA), Random Forest, and Support Vector Machine. Correlations between probability of seizure presentation (p(SP)) and continuous variables were also analyzed. RESULTS: Our cohort consisted of 59 seizure-presenting and 77 non-seizure-presenting patients. All models consistently demonstrated significant correlations (p < 0.05) between (p(SP)) and the following variables: T1Gd radius (-0.781 to -0.674), T2-FLAIR (-0.674 to -0.611), and log(/D) (0.169 to 0.294). Age was significant (p < 0.05) in 4 of the 5 models (-0.211 to -0.175). Mean performance measures for the models (and the best performer) were: 0.726 for Area under the ROC curve (0.75 with NB), 0.6202 for sensitivity (0.661 with NB), 0.74 for specificity (0.766 with LDA). The 5 features ranked as most important were: T1Gd, T2-FLAIR, log(/D), age, log(). CONCLUSIONS: We found an association in seizure-presenting patients with smaller, more proliferative tumors and younger age. Machine learning predictive modeling can potentially be informative in the clinical arena. Further validation studies, to determine the degree of data overfitting, model versatility, as well as performance on test data, are warranted.