Identification, validation and biological characterisation of novel glioblastoma tumour microenvironment subtypes: implications for precision immunotherapy

BACKGROUND

New precision medicine therapies are urgently required for glioblastoma (GBM). However, to date, efforts to subtype patients based on molecular profiles have failed to direct treatment strategies. We hypothesised that interrogation of the GBM tumour microenvironment (TME) and identification of novel TME-specific subtypes could inform new precision immunotherapy treatment strategies.

MATERIALS AND METHODS

A refined and validated microenvironment cell population (MCP) counter method was applied to >800 GBM patient tumours (GBM-MCP-counter). Specifically, partition around medoids (PAM) clustering of GBM-MCP-counter scores in the GLIOTRAIN discovery cohort identified three novel patient clusters, uniquely characterised by TME composition, functional orientation markers and immune checkpoint proteins. Validation was carried out in three independent GBM-RNA-seq datasets. Neoantigen, mutational and gene ontology analysis identified mutations and uniquely altered pathways across subtypes. The longitudinal Glioma Longitudinal AnalySiS (GLASS) cohort and three immunotherapy clinical trial cohorts [treatment with neoadjuvant/adjuvant anti-programmed cell death protein 1 (PD-1) or PSVRIPO] were further interrogated to assess subtype alterations between primary and recurrent tumours and to assess the utility of TME classifiers as immunotherapy biomarkers.

RESULTS

TME^High tumours (30%) displayed elevated lymphocyte, myeloid cell immune checkpoint, programmed cell death protein 1 (PD-1) and cytotoxic T-lymphocyte-associated protein 4 transcripts. TME^High/mesenchymal+ patients featured tertiary lymphoid structures. TME^Med (46%) tumours were enriched for endothelial cell gene expression profiles and displayed heterogeneous immune populations. TME^Low (24%) tumours were manifest as an 'immune-desert' group. TME subtype transitions upon recurrence were identified in the longitudinal GLASS cohort. Assessment of GBM immunotherapy trial datasets revealed that TME^High patients receiving neoadjuvant anti-PD-1 had significantly increased overall survival (P = 0.04). Moreover, TME^High patients treated with adjuvant anti-PD-1 or oncolytic virus (PVSRIPO) showed a trend towards improved survival.

CONCLUSIONS

We have established a novel TME-based classification system for application in intracranial malignancies. TME subtypes represent canonical 'termini a quo' (starting points) to support an improved precision immunotherapy treatment approach.

P10.17.B Androgen receptor signalling has independent prognostic value and is associated with resistance to temozolomide and tumour growth in glioblastoma

Background

The spectrum of sex differences in glioblastoma (GBM) is a rising topic. GBM is more frequent (sex ratio M:F=1.6) and seem to have worse prognosis in males (M) compared to females. Androgen receptor (AR) is expressed in GBM but its role in GBM biology and clinics is not yet fully understood.

Material and Methods

We have selected from our tumour tissue bank: (i) GBM samples with clinical, pathological and molecular annotations and (ii) GBM-patient derived cell lines -GBM-PDCL- with transcriptomic profiling and available tissue from parental tumours.AR expression was investigated using multiplex IHC Opal stained-FFPE sections. Sensitivity of PDCL to hormonotherapy (dihydrotestosterone or enzalutamide EZT) alone or combined with temozolomide (TMZ) was tested using WST1 viability assay. Response to drugs were correlated to AR expression (high: log2 [FPKM] >2 or low: log2 [FPKM] <2). PDCLs dihydrotestosterone secretion was examined by ELISA in both supernatants and cell lysates.

Results

In our cohort of 1188 IDH muted GBM patients, the median follow-up was 27.7 (25.8-31.3). The sex ratio M:F was 1.68 and the median age at diagnosis was 57 years [18-69]. Compared to males, females have better prognosis with a median overall survival of 26.1 months vs 20.8 months (p<0.0001). In multivariate analysis, male gender was an independent prognostic factor with a HR of 1.43 CI95%(1.21-1.69) (p<0.0001) adjusted on age, MGMT status and Karnofsky indice . Although interaction between age and sex was not significant (p = 0.854), HR of gender varied through age category 1.53 CI95%( 1.01-2.33 ), 1.51 CI95%( 0.888-2.55 ) and 1.38 CI95%( 0.857-2.21) for males <45 years, between 46-55 years and > 55 years respectively. In GBM tumours, AR staining was positive in 5- 30% tumour cells. After 72h exposure to 10nM dihydrotestosterone, ARhigh GBM-PDCLs exhibit AR nuclear translocation and c-myc expression upregulation, suggesting functional activation of the receptor. AR pharmacological blockade increased temozolomide sensitivity in ARhigh GBM-PDCLs regardless of the patient’s sex. Moreover, ARhigh PDCLs secreted DHT, suggesting autocrine and paracrine AR stimulation in tumour microenvironment.

Conclusion

Our study highlights the independent prognostic significance of sex in GBM patients, where males exhibit a worse prognosis. In GBM-PDCL, we show that androgen receptor is (i) functional in vitro and (ii) associated with resistance to temozolomide. Further studies are needed to investigate the role of AR in tumour tissue microenvironment.

Intestinal UDP-Glucuronosyltransferase 1A1 and Protection against Irinotecan-Induced Toxicity in a Novel UDP-Glucuronosyltransferase 1A1 Tissue-Specific Humanized Mouse Model

The human UDP-glucuronosyltransferases (UGTs) represent an important family of drug-metabolizing enzymes, with UGT1A1 targeting the conjugation and detoxification of many exogenous substances, including pharmaceutical drugs. In this study we generated humanized UGT1A1 mice expressing the human UGT1A1 gene in either liver (hUGT1A1^HEP ) or intestine (hUGT1A1^GI ), enabling experiments to examine tissue-specific properties of UGT1A1-specific glucuronidation. Hepatic and intestinal tissue-specific expression and function of UGT1A1 were demonstrated. Although the liver is considered a major organ for detoxification, intestinal UGT1A1 is an important contributor for drug clearance. Mice were challenged with irinotecan (CPT-11), a prodrug hydrolyzed by carboxylesterases to form the active metabolite 7-ethyl-10-hydroxycamptothecin (SN-38) and detoxified by UGT1A1. Humanized UGT1A1^HEP mice that have no intestinal UGT1A1 displayed a greater lethality rate when exposed to CPT-11 than hUGT1A1^GI mice. When exposed to a low dose of CPT-11 (10 mg/kg), hUGT1A1^HEP mice displayed greater intestinal inflammatory (IL-1β and IL-6) insult in addition to p53-triggered apoptotic responses. In vitro studies with intestinal crypt organoids exposed to CPT-11 confirmed the results observed in vivo and indicated that CPT-11 impacts stemness, apoptosis, and endoplasmic reticulum (ER) stress in organoids deficient in UGT1A1. When we examined the induction of ER stress in organoids with thapsigargin, an inhibitor of sarco/endoplasmic reticulum Ca²⁺ ATPase, apoptosis and the caspase surge that occurred in hUGT1A1^HEP mice were blocked in hUGT1A1^GI organoids. This study reveals the importance of intestinal UGT1A1 in preventing inflammation, apoptosis, and loss of stemness capacity upon systemic challenge with an important chemotherapeutic agent. SIGNIFICANCE STATEMENT: Hepatic and intestinal UGT1A1 play a key role in the metabolism and detoxification of endogenous and exogenous compounds. The use of tissue-specific humanized models expressing UGT1A1 in liver or intestine has confirmed the relevance of the intestinal tract in the detoxification of irinotecan. Mechanistic studies using intestinal organoids highlighted the importance of UGT1A1 in reducing inflammation, apoptosis, and loss of stemness. These new models provide valuable tools for studying tissue-specific glucuronidation of substances that are metabolized by human UGT1A1.

New hints towards a precision medicine strategy for IDH wild-type glioblastoma

Glioblastoma represents the most common primary malignancy of the central nervous system in adults and remains a largely incurable disease. The elucidation of disease subtypes based on mutational profiling, gene expression and DNA methylation has so far failed to translate into improved clinical outcomes. However, new knowledge emerging from the subtyping effort in the IDH-wild-type setting may provide directions for future precision therapies. Here, we review recent learnings in the field, and further consider how tumour microenvironment differences across subtypes may reveal novel contexts of vulnerability. We discuss recent treatment approaches and ongoing trials in the IDH-wild-type glioblastoma setting, and propose an integrated discovery stratagem incorporating multi-omics, single-cell technologies and computational approaches.

Determination of an optimal dosing schedule for combining Irinophore C™ and temozolomide in an orthotopic model of glioblastoma

Our laboratory reported that Irinophore C™ (IrC™; a lipid-based nanoparticulate formulation of irinotecan) is effective against an orthotopic model of glioblastoma (GBM) and that treatment with IrC™ was associated with vascular normalization within the tumor. Here, the therapeutic effects of IrC™ when used in combination with temozolomide (TMZ) in concurrent and sequential treatment schedules were tested. It was anticipated that IrC™ engendered vascular normalization would increase the delivery of TMZ to the tumor and that this would be reflected by improved treatment outcomes. The approach compared equally efficacious doses of irinotecan (IRN; 50 mg/kg) and IrC™ (25 mg/kg) in order to determine if there was a unique advantage achieved when combining TMZ with IrC™. The TMZ sensitive U251MG(O) cell line (null expression of O-6-methylguanine-DNA methyltransferase (MGMT)) modified to express the fluorescent protein mKate2 was inoculated orthotopically into NOD.CB17-SCID mice and treatment was initiated 14 days later. Our results demonstrated that IrC™ and TMZ administered concurrently resulted in optimal treatment outcomes, with 50% long term survivors (>180 days) in comparison to 17% long term survivors in animals treated with IRN and TMZ or TMZ alone. Indeed, the different treatments resulted in a 353%, 222% and 280% increase in median survival time (MST) compared to untreated animals for, respectively, IrC™ combined with TMZ, IRN combined with TMZ, and TMZ alone. When TMZ was administered after completion of IRN or IrC™ dosing, an increase in median survival time of 167-174% was observed compared to untreated animals and of 67% and 74%, respectively, when IRN (50 mg/kg) and IrC™ (25mg/kg) were given as single agents. We confirmed in these studies that after completion of the Q7D×3 dosing of IrC™, but not IRN, the tumor-associated vascular was normalized as compared to untreated tumors. Specifically, reductions in the fraction of collagen IV-free CD31 staining (p<0.05) and reductions in tumor vessel diameter were observed in tumors from IrC™-treated animals when compared to tumors from untreated or IRN treated animals. Analysis by transmission electron microscopy of the ultra-structure of tumors from IrC™-treated and untreated animals revealed that tumor-associated vessels from treated animals were smaller, more organized and exhibited a morphology comparable to normal blood vessels. In conclusion, optimal treatment outcomes were achieved when IrC™ and TMZ were administered concurrently, whereas IrC™ followed by TMZ treatment given sequentially did not confer any therapeutic advantage.

Profiling serum bile acid glucuronides in humans: gender divergences, genetic determinants, and response to fenofibrate

Glucuronidation, catalyzed by UDP-glucuronosyltransferase (UGT) enzymes detoxifies cholestatic bile acids (BAs). We aimed at i) characterizing the circulating BA-glucuronide (-G) pool composition in humans, ii) evaluating how sex and UGT polymorphisms influence this composition, and iii) analyzing the effects of lipid-lowering drug fenofibrate on the circulating BA-G profile in 300 volunteers and 5 cholestatic patients. Eleven BA-Gs were determined in pre- and post-fenofibrate samples. Men exhibited higher BA-G concentrations, and various genotype/BA-G associations were discovered in relevant UGT genes. The chenodeoxycholic acid-3G concentration was associated with the UGT2B7 802C>T polymorphism. Glucuronidation assays confirmed the predominant role of UGT2B7 and UGT1A4 in CDCA-3G formation. Fenofibrate exposure increased the serum levels of 5 BA-G species, including CDCA-3G, and up-regulated expression of UGT1A4, but not UGT2B7, in hepatic cells. This study demonstrates that fenofibrate stimulates BA glucuronidation in humans, and thus reduces bile acid toxicity in the liver.

Development of glioblastoma cell lines expressing red fluorescence for non-invasive live imaging of intracranial tumors

Glioblastoma (GBM) cell lines expressing red fluorescent proteins were evaluated as a tool for non-invasive imaging of orthotopic tumors.

MATERIALS AND METHODS

mKate2- and mCherry-transduced U251MG GBM lines were sorted by flow cytometry. The growth rates and drug sensitivity of the resulting cell lines were compared to those of the parental line. Following orthotopic implantation, mKate2-expressing cells were detected using multispectral imaging.

RESULTS

Flow cytometry-sorted fluorescent populations exhibiting growth curves that were comparable to those of the parental line were selected. mKate2-expressing cells were inoculated orthotopically and formed tumors which were visualized non-invasively, allowing monitoring of tumor growth over time and the assessment of tumor response to temozolomide drug treatment.

CONCLUSION

The strategy reported here led to the successful development of GBM models expressing mKate2 or mCherry. The fluorescence signal intensity measured in the brain of live animals correlates with tumor size, thus providing a method to assess tumor progression and response to treatment.

Role of the Endothelin System in Secondary Pulmonary Hypertension Related to Air Embolism: Lessons Learned from Testing Four Classes of Endothelin Blockers in a Rat Model

A rat model of acute pulmonary air embolism (APAE) was developed. These animals had a higher right ventricular systolic pressure (RVSP) (+ 69% at 15-minute peak, and 21-34% at 30-180 minutes), as well as a reduced mean arterial blood pressure (10-20% at 60-180 minutes), heart rate (20-26% at 60-180 minutes) and PaO2 (9-11% at 30-180 minutes) compared with control rats. The role of the endothelin (ET) system, known to be involved in pulmonary hypertension of various etiologies, was investigated by evaluating the effect of the four classes of ET blockers: ET-converting enzyme inhibitor (ECEi) (CGS 35066), selective endothelin-A receptor antagonist (ETA-Ra) (Atrasentan, ABT-627), endothelin-B receptor antagonist (ETB-Ra) (A-192621) or mixed endothelin-A/endothelin-B receptor antagonist (ETA/B-Ra) (A-182086) in this animal model. All four were effective, to various degrees, at reducing the APAE-induced rise in RVSP. The relative efficacy of those compounds in reducing the acute elevation (15 minutes) of RVSP was ECEi >or= ETA/B-Ra >> ETA-Ra = ETB-Ra. The sustained elevation (30-180 minutes) of RVSP was totally abolished by ECEi and attenuated by other ET blockers with a relative efficacy of ETA-Ra > ETA/B-Ra >or= ETB-Ra. ET receptor antagonists did not affect right ventricular basal tone (control rats) whereas ECEi reduced it by up to 12% after 2 hours. The APAE reduction in mean arterial blood pressure was unaffected by ETARa, was completely normalized by ETB-Ra, but was further reduced by either ETA/B-Ra or ECEi. The basal mean arterial blood pressure in control rats was unaffected by ETA-Ra, was elevated by ETB-Ra, but was depressed by ETA/B-Ra and ECEi. All ET blockers maintained normal oxygen saturation in APAE. These results support a role for ETs in rat APAE, since ET blockers can attenuate the cardiopulmonary deterioration and blood gas exchange. However, modulation of the central hemodynamic profile is more complex and may limit the usefulness of some ET blockers.

Fluctuations in spatial recognition memory across the menstrual cycle in female rhesus monkeys

Findings are inconsistent regarding whether women's cognitive performance fluctuates across phases of the menstrual cycle, but differences in methodology and the use of reported cycle phase rather than precise hormonal measures may underlie these disparities. Studies in monkeys may help resolve these discrepant findings, since hormonal status can be reliably determined. We tested four young (5-7 years old) female rhesus monkeys daily during one entire menstrual cycle on three cognitive tasks displayed on a computerized touch-screen system: a Matching to Sample task with a 30 s delay (MTS-30s), a Matching to Sample task without delay (MTS-no delay) and the spatial condition of the Delayed Recognition Span Test (spatial-DRST). Blood samples were collected at specific time intervals throughout the cycle and assayed for estradiol and progesterone in order to identify hormonal status. There was a nonsignificant trend for the MTS-30s scores to be better during the follicular and luteal phases, when estradiol levels were low, than during the peri-ovulatory phase, when estradiol levels were at their highest. MTS-no delay performance did not vary as a function of hormonal status. Spatial-DRST scores were significantly better during the follicular and luteal phases than during the peri-ovulatory phase of the cycle. These data in the female rhesus monkey support the hypothesis that spatial memory performance is sensitive to estradiol variations across the menstrual cycle, with better performance associated with low estradiol levels.