Rubin, Joshua B

Glioblastoma (GBM) is an aggressive brain tumor that often progresses despite resection and treatment. Timely and continuous assessment of GBM progression is critical to expedite secondary surgery or enrollment in clinical trials. However, current progression detection requires costly and specialized MRI examinations, which, in the absence of new symptoms or signs, are usually scheduled every 2–3 months. Here, we hypothesized that changes in daily activity are associated with GBM growth and disease progression. We found that wheel-running activity in GBM-bearing mice declined as tumors grew and preceded weight loss and circadian breakdown by over a week. Temozolomide treatment in the morning, but not evening, significantly reduced tumor size and restored daily locomotion in mice. In a pilot study of 6 patients with GBM wearing an actigraphy watch, wrist movement provided a feasible and continuous longitudinal indicator of daily activity with 1-minute resolution. After tumor resection and radiation, daily activity declined in 2 patients 19 and 55 days before detection of progression by MRI but did not change for the 4 patients with stable disease. These results suggest that daily activity tracking using wearable devices may serve as a real-time indicator and potential monitoring tool for GBM progression and treatment efficacy.

Abstract Glioblastoma (GBM) is the most common primary tumor in adults with a median survival of 15 months post-treatment. GBM incidence and progression depends on sex, with reports of 1.6x increased incidence and 6-month shorter overall survival in males, compared to females. Regardless of sex, the standard of care for GBM patients includes surgical resection, followed by chemotherapy and radiotherapy, and often dosing of the synthetic glucocorticoid Dexamethasone (DEX) to reduce brain edema post-treatment. Our lab recently found that DEX can suppress or promote GBM growth depending on the time of day the treatment is administered. Because glucocorticoid secretion depends on circadian time and sex, we tested whether endogenous glucocorticoids regulate sex differences in GBM growth. We transduced a murine chromosomal male or female GBM cell line (Nf1-/-DNp53 XX or XY) with a transcriptional luciferase reporter of the circadian clock gene Per2 (Nf1-/-DNp53-P2L). These cells were next treated with DEX at the daily peak or trough of Per2 expression and growth was measured 48 hours post-treatment. We found that administering DEX at peak Per2 suppresses growth in males, while administering DEX at trough Per2 promotes growth in both males and females. We then knocked-down the glucocorticoid receptor (GR) using an shRNA lentiviral construct in two GBM cell models (GL261 and LN229) and implanted WT or GR-KD cells into the basal ganglia of male and female mice. Loss of GR in GBM cells significantly suppressed GBM growth, especially in male mice, as evidenced by postmortem Ki67 staining. These results suggest that GBM growth in males may be more sensitive to endogenous glucocorticoids, and that sensitivity varies with time of day. This work will inform future personalized circadian interventions to improve GBM patient outcomes. Supported by NIH grant R01NS134885.

Abstract Sex differences in cancer outcome, including that of glioblastoma (GBM), are shaped by biological, hormonal, and immunological factors and impact disease progression, treatment responses, survival, and tumor microenvironment (TME) interactions. Platelets regulate the immune responses and tumor progression in many cancers, but it is not clear how they contribute to these sex-based differences by affecting the dynamics of the TME. Here, we show that GBM patients exhibit heightened platelet reactivity driven by PAR4 signaling. In murine GBM models, both pharmacological inhibition of PAR4 using BMS986120 and genetic deletion of PAR4 significantly prolong survival in females but not males. This survival advantage is estrogen dependent: it is preserved in chromosomal male-hormonal female mice within the four-core genotype model and is rescued in ovariectomized mice treated with estrogen. The survival benefit is TME specific and is mediated by platelet-driven enhancement of CD8⁺ T cell infiltration into the tumor. Inhibition of platelet PAR4 signaling increases calcium signaling through an estrogen-dependent interaction between PAR4 and estrogen receptor β (ERβ)—a receptor interaction not previously described. PAR4-activated platelets within the TME suppress CD8⁺ T cell function, and depletion of CD8⁺ T cells abolishes both the tumor-induced platelet reactivity and the survival benefit conferred by PAR4 inhibition. These findings establish platelet-mediated PAR4 signaling as a critical driver of tumor progression and identify sex-specific immune responses as key to therapeutic efficacy.

High-grade gliomas, like glioblastoma multiforme (GBM), are the most common malignant brain tumors in adults and are treated with the chemotherapy drug temozolomide (TMZ). In humans, a retrospective analysis of patients’ overall survival suggests that morning dosing may confer a benefit over evening dosing. Circadian variation in O6-methylguanine-DNA methyltransferase (MGMT) gene expression and promoter methylation has been implicated in increased tumor cell sensitivity to TMZ in the morning. Although patient compliance with timed oral administration of TMZ was high in a prospective trial, it is not known whether differences in daily sleep patterns of patients impact the biological time of drug administration or overall survival. Using wrist actigraphy collected from 10 high-grade glioma patients, we quantified the moment of oral TMZ delivery in terms of wall clock time and internal biological time during the months after surgical tumor resection. We found that variation of daily rhythms within and between individuals caused dosing times to vary more in their internal biological time than wall clock time so that, for example, some doses taken by patients assigned for the evening (2000 h) were closer to the patient’s internal biological morning. We conclude that wrist actigraphy provides a reliable and non-invasive estimate of personal circadian time that could improve efficacy and precision of TMZ delivery. These findings may inform personalized circadian medicine and optimized times for TMZ delivery in the clinic.

Abstract Characterizing intra- and inter-tumoral heterogeneity of glioblastoma has historically relied on discrete classifications of malignant cell populations leaving immune and other cell populations, known to exist admixed with the malignant tumor cells, relatively neglected. Manifold learning algorithms can manage deconvolving multiple cell populations and are often used to track cell state transitions in single cell transcriptomics. We applied a manifold learning approach to TCGA microarray data (Nf525) and bulk transcriptomics of 134 image localized biopsies across 30 patients with primary and 9 with recurrent glioblastoma to further elucidate how to organize biopsies across a spectrum of possible tissue states. The algorithm revealed a low-dimensional manifold graph for which each biopsy lives across 3 polarizing tissue states - one that is associated with diffusely invaded brain, one that is enriched in mesenchymal genes, and one that is enriched in classical proliferative tumor signatures. We deconvolved the bulk transcriptomics of the image-localized biopsies to reveal the relative abundance of 18 malignant, immune, and other cell subpopulations in each biopsy. Overlaying the cellular decomposition onto the manifold graph visualizing the tissue state distributions revealed that transitions between states correlate with changes in cellular cohabitation composition. The tumor cellular cohabitation ecologies have the lowest diversity, as inferred by ecological measures such as Shannon entropy and evenness, at the distal poles of the graph when compared to the transitional arms. Further, we found that the relationship between imaging appearance of contrast enhancement on T1-weighted MRI and the biopsy cellular composition varies with sex and primary vs recurrent biopsy status. The limited spectrum of possible tissue states revealed by the manifold learning is suggestive of a limited continuum along which tumor and non-tumoral cell populations can cohabitate. Such a limited low-dimensional biological space may constrain the dynamics of tumor biology in a predictable manner.

Abstract BACKGROUND: Endpoints in clinical trials using novel treatments are evaluated by RANO criteria, which provide an estimate of tumor size from two-dimensional measurements along the most prominent axial slice. However, pediatric low-grade gliomas (pLGG) commonly have variegated shapes with solid and cystic components, potentially resulting in misevaluation of true tumor volume and thus, ultimately, trial outcome. OBJECTIVES: We aim to characterize treatment response through volumetric assessment of progressive/recurrent pLGGs treated with single-agent everolimus on PNOC001 clinical trial. We seek to identify clinically relevant criteria that provide added value to response assessment beyond 2D measurements. METHODS: In a cohort of 44 patients we performed 3D-segmentation of solid, cystic and whole tumor within our PACS-framework and compared results to previously carried-out central imaging review by RANO criteria which had yielded 15 PD, 27 SD, 2 PR and 0 CR. RESULTS: 8 tumors were solid only and 36 had a mixed solid-cystic appearance. When evaluating the entire tumor (i.e. solid and cystic components combined) and using the same RANO cutoff criteria, one case changed from PR to SD, one changed from SD to PR, 3 changed from SD to PD, and 7 changed from PD to SD, resulting in an overall discordance of 27% of cases. CONCLUSION: We propose that incorporation of volumetrics into response assessment provides additional and potentially more accurate information beyond RANO-based measurements. It is crucial to note that the above-reported changes represent numerical discrepancies as opposed to true-to-reality changes in clinical outcome. Determining representative thresholds for the deployment of volumetric measures in clinical trials will be critical. Future work will include data from the PNOC002 clinical trial and evaluate inter-reader agreement and reader discordance. With the availability of PACS-based 3D-tools in neuroradiology practice, well-defined volumetric criteria could be incorporated prospectively into treatment response analysis in clinical trials.

Abstract Background Magnetic Resonance Imaging (MRI) is the mainstay for neurosurgical oncology but not for informing us about glioma biology. An obstacle to developing MR-based glioma biomarkers is the absence of rigorous correlation between MRI features and glioma biology as assessed in multi-regional biopsies, within and across patients. Methods We directly addressed this obstacle by collating a unique cohort of 202 MRI-localized biopsies from 58 patients. We define a low-dimensional transcriptional pseudotime continuum along which heterogeneous high-grade glioma (HGG) samples organize both within and across patients. Results We observe three polarized transcriptional tissue states: infiltrated brain, immune/inflammatory, and proliferative associated with patterns of cohabitation of cellular subpopulations. The states and deconvolved populations show correlation with enhancement status on T1Gd MRI. Moreover, discrete MRI habitats, regions sharing common imaging features, defined as combinations of high or low signal intensity across multiparametric MRI revealed 14 MRI habitats. We order the MRI habitats according to the average pseudotime on the transcriptional continuum. We find that MRI habitats with low pseudotime (associated with early tumor development and diffusely invaded brain tissue) localized at the periphery of the tumor whilst high pseudotime either proliferative or immune/inflammatory states were towards the core of the lesion. We find that composition of MRI habitats are impacted by MGMT status and patient sex. Conclusion This suggests that ongoing aggregation of MRI-localized biopsies may augment our projection of biology onto MRI habitats to support tracking of the evolution of cellular ecologies within and across each patient’s tumor over time. Key Points Trajectory inference of HGG samples reveals a continuum connecting three polarized tissue states reflecting transitions in cellular population ecologies and gene set enrichment Combining multiparametric MRI features into regional MRI habitats reveal clinically-distinct patterns of cell co-habitation and tissue state MGMT methylation status shifts the cellular patterns of cohabitation within MRI habitats Importance of the Study Regional heterogeneity is a defining characteristic of high-grade glioma (HGG). Mapping the transcriptomic biology of multiregionally sampled MRI-localized biopsies from high-grade gliomas onto a low-dimensional space reveals a continuum between three primary tissue states: infiltrated brain, immune/inflammatory, and proliferative. Each state is associated with distinct cell cohabitation ecologies. Transitions between the tissue states coordinate with changes in tumor composition, revealing potential direct connections between tumor biology shifts and their manifestation on clinical imaging.

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

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

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