Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts

BACKGROUND

A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials.

METHODS

Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry, and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors.

RESULTS

We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components.

CONCLUSIONS

Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.

Early-life influenza A (H1N1) infection independently programs brain connectivity, HPA AXIS and tissue-specific gene expression profiles

Early-life adversity covers a range of physical, social and environmental stressors. Acute viral infections in early life are a major source of such adversity and have been associated with a broad spectrum of later-life effects outside the immune system or "off-target". These include an altered hypothalamus-pituitary-adrenal (HPA) axis and metabolic reactions. Here, we used a murine post-natal day 14 (PND 14) Influenza A (H1N1) infection model and applied a semi-holistic approach including phenotypic measurements, gene expression arrays and diffusion neuroimaging techniques to investigate HPA axis dysregulation, energy metabolism and brain connectivity. By PND 56 the H1N1 infection had been resolved, and there was no residual gene expression signature of immune cell infiltration into the liver, adrenal gland or brain tissues examined nor of immune-related signalling. A resolved early-life H1N1 infection had sex-specific effects. We observed retarded growth of males and altered pre-stress (baseline) blood glucose and corticosterone levels at PND42 after the infection was resolved. Cerebral MRI scans identified reduced connectivity in the cortex, midbrain and cerebellum that were accompanied by tissue-specific gene expression signatures. Gene set enrichment analysis confirmed that these were tissue-specific changes with few common pathways. Early-life infection independently affected each of the systems and this was independent of HPA axis or immune perturbations.

Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts

Background

A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials.

Methods

Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA-sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors.

Results

We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components.

Conclusions

Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.

Cross-Talk between miRNAs from the Dlk1-Dio3 Locus and Histone Methylation to Protect Male Cerebellum from Methyl Donor Deficiency

SCOPE

Disruption of the one carbon metabolism during development, i.e., following a gestational vitamin B9 and B12 deficiencies, is involved in birth defects and brain development delay. Using a rat nutritional model, consisting of pups born to dams fed a vitamin B9 and B12 deficient diet (MDD), the study previously reports molecular and cellular alterations in the brain, in a sex dependent manner, with females being more affected than males. The study hypothesizes that epigenetic modifications could participate in the sex differences is observed.

METHODS AND RESULTS

The study investigates lysine methylation of histones and expression of microRNAs in the cerebellum of MDD male and female pups. The study reports a differential regulation of H3K36Me2 and H4K20Me3 between males and females, in response to MDD. Moreover, distinct regulation of Kmt5b and Kdm2a expression by miR-134-5p and miR-369-5p from the Dlk1-Dio3 locus, contributes to the maintenance of expression of genes involved in synaptic plasticity.

CONCLUSION

These results could explain the neuroprotection to MDD that male pups display. The work will contribute to the understanding of the consequences of vitamin starvation on brain development, as well as how the epigenome is affected by one carbon metabolism disruption.

DeSUMOylation of chromatin-bound proteins limits the rapid transcriptional reprogramming induced by daunorubicin in acute myeloid leukemias

Genotoxicants have been used for decades as front-line therapies against cancer on the basis of their DNA-damaging actions. However, some of their non-DNA-damaging effects are also instrumental for killing dividing cells. We report here that the anthracycline Daunorubicin (DNR), one of the main drugs used to treat Acute Myeloid Leukemia (AML), induces rapid (3 h) and broad transcriptional changes in AML cells. The regulated genes are particularly enriched in genes controlling cell proliferation and death, as well as inflammation and immunity. These transcriptional changes are preceded by DNR-dependent deSUMOylation of chromatin proteins, in particular at active promoters and enhancers. Surprisingly, inhibition of SUMOylation with ML-792 (SUMO E1 inhibitor), dampens DNR-induced transcriptional reprogramming. Quantitative proteomics shows that the proteins deSUMOylated in response to DNR are mostly transcription factors, transcriptional co-regulators and chromatin organizers. Among them, the CCCTC-binding factor CTCF is highly enriched at SUMO-binding sites found in cis-regulatory regions. This is notably the case at the promoter of the DNR-induced NFKB2 gene. DNR leads to a reconfiguration of chromatin loops engaging CTCF- and SUMO-bound NFKB2 promoter with a distal cis-regulatory region and inhibition of SUMOylation with ML-792 prevents these changes.

Impact of IDH Mutations, the 1p/19q Co-Deletion and the G-CIMP Status on Alternative Splicing in Diffuse Gliomas

By generating protein diversity, alternative splicing provides an important oncogenic pathway. Isocitrate dehydrogenase (IDH) 1 and 2 mutations and 1p/19q co-deletion have become crucial for the novel molecular classification of diffuse gliomas, which also incorporates DNA methylation profiling. In this study, we have carried out a bioinformatics analysis to examine the impact of the IDH mutation, as well as the 1p/19q co-deletion and the glioma CpG island methylator phenotype (G-CIMP) status on alternative splicing in a cohort of 662 diffuse gliomas from The Cancer Genome Atlas (TCGA). We identify the biological processes and molecular functions affected by alternative splicing in the various glioma subgroups and provide evidence supporting the important contribution of alternative splicing in modulating epigenetic regulation in diffuse gliomas. Targeting the genes and pathways affected by alternative splicing might provide novel therapeutic opportunities against gliomas.

Actin cytoskeleton depolymerization increases matrix metalloproteinase gene expression in breast cancer cells by promoting translocation of cysteine-rich protein 2 to the nucleus

The actin cytoskeleton plays a critical role in cancer cell invasion and metastasis; however, the coordination of its multiple functions remains unclear. Actin dynamics in the cytoplasm control the formation of invadopodia, which are membrane protrusions that facilitate cancer cell invasion by focusing the secretion of extracellular matrix-degrading enzymes, including matrix metalloproteinases (MMPs). In this study, we investigated the nuclear role of cysteine-rich protein 2 (CRP2), a two LIM domain-containing F-actin-binding protein that we previously identified as a cytoskeletal component of invadopodia, in breast cancer cells. We found that F-actin depolymerization stimulates the translocation of CRP2 into the nucleus, resulting in an increase in the transcript levels of pro-invasive and pro-metastatic genes, including several members of the MMP gene family. We demonstrate that in the nucleus, CRP2 interacts with the transcription factor serum response factor (SRF), which is crucial for the expression of MMP-9 and MMP-13. Our data suggest that CRP2 and SRF cooperate to modulate of MMP expression levels. Furthermore, Kaplan-Meier analysis revealed a significant association between high-level expression of SRF and shorter overall survival and distant metastasis-free survival in breast cancer patients with a high CRP2 expression profile. Our findings suggest a model in which CRP2 mediates the coordination of cytoplasmic and nuclear processes driven by actin dynamics, ultimately resulting in the induction of invasive and metastatic behavior in breast cancer cells.

OS05.5.A Glioblastoma-instructed microglia transit to heterogeneous phenotypic states with dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts

Background

A major contributing factor to Glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive tumor microenvironment (TME). GBM-associated myeloid cells, including resident microglia, macrophages and other peripheral immune cells are generally geared towards tumor-supportive roles. It is however unclear whether such recruited myeloid cells are phenotypically and functionally identical. Here, we aim to understand the heterogeneity of the GBM TME, using an unbiased, marker-free approach to systematically characterize cell type identities at the molecular and functional levels.

Material and Methods

We applied single-cell RNA-sequencing, multicolor flow cytometry, immunohistochemical analyses and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM patient-derived orthotopic xenografts (PDOXs) representing different tumor phenotypes were compared to glioma mouse GL261 model and patient tumors.

Results

We show that PDOX models recapitulate major components of the TME found in human GBM. Human GBM cells reciprocally interact with mouse cells to create a GBM-specific TME. The most prominent transcriptomic adaptations are found in tumor-associated macrophages (TAMs), which are largely of microglial origin. We reveal inter-patient heterogeneity of TAMs and identify key signatures of distinct phenotypic states within the microglia-derived TAMs across distinct GBM landscapes. GBM-educated microglia adapt expression of genes involved in immunosuppression, migration, phagocytosis and antigen presentation, indicating functional cross-talk with GBM cells. We identify novel phenotypic states with astrocytic and endothelial-like features. Identified gene signatures and phenotypic states are confirmed in GBM patient tumor tissue. Finally we show that temozolomide treatment leads to transcriptomic adaptation of not only the GBM tumor cells but also adjacent TME components.

Conclusion

Our data provide insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumor. We confirm a crucial role of microglia in supporting the immunosuppressive TME and show that PDOXs allow to monitor the highly plastic GBM ecosystem and its phenotypic adaptations upon treatment. This work further confirms the clinical relevance of PDOX avatars for testing novel therapeutics including modalities designed to target the myeloid compartment.

Neurodegeneration and neuroinflammation are linked, but independent of alpha‐synuclein inclusions, in a seeding/spreading mouse model of Parkinson's disease

A key pathological process in Parkinson's disease (PD) is the transneuronal spreading of α‐synuclein. Alpha‐synuclein (α‐syn) is a presynaptic protein that, in PD, forms pathological inclusions. Other hallmarks of PD include neurodegeneration and microgliosis in susceptible brain regions. Whether it is primarily transneuronal spreading of α‐syn particles, inclusion formation, or other mechanisms, such as inflammation, that cause neurodegeneration in PD is unclear. We used a model of spreading of α‐syn induced by striatal injection of α‐syn preformed fibrils into the mouse striatum to address this question. We performed quantitative analysis for α‐syn inclusions, neurodegeneration, and microgliosis in different brain regions, and generated gene expression profiles of the ventral midbrain, at two different timepoints after disease induction. We observed significant neurodegeneration and microgliosis in brain regions not only with, but also without α‐syn inclusions. We also observed prominent microgliosis in injured brain regions that did not correlate with neurodegeneration nor with inclusion load. Using longitudinal gene expression profiling, we observed early gene expression changes, linked to neuroinflammation, that preceded neurodegeneration, indicating an active role of microglia in this process. Altered gene pathways overlapped with those typical of PD. Our observations indicate that α‐syn inclusion formation is not the major driver in the early phases of PD‐like neurodegeneration, but that microglia, activated by diffusible, oligomeric α‐syn, may play a key role in this process. Our findings uncover new features of α‐syn induced pathologies, in particular microgliosis, and point to the necessity for a broader view of the process of α‐syn spreading.

Mitochondrial inhibitors circumvent adaptive resistance to venetoclax and cytarabine combination therapy in acute myeloid leukemia

Therapy resistance represents a major clinical challenge in acute myeloid leukemia (AML). Here we define a 'MitoScore' signature, which identifies high mitochondrial oxidative phosphorylation in vivo and in patients with AML. Primary AML cells with cytarabine (AraC) resistance and a high MitoScore relied on mitochondrial Bcl2 and were highly sensitive to venetoclax (VEN) + AraC (but not to VEN + azacytidine). Single-cell transcriptomics of VEN + AraC-residual cell populations revealed adaptive resistance associated with changes in oxidative phosphorylation, electron transport chain complex and the TP53 pathway. Accordingly, treatment of VEN + AraC-resistant AML cells with electron transport chain complex inhibitors, pyruvate dehydrogenase inhibitors or mitochondrial ClpP protease agonists substantially delayed relapse following VEN + AraC. These findings highlight the central role of mitochondrial adaptation during AML therapy and provide a scientific rationale for alternating VEN + azacytidine with VEN + AraC in patients with a high MitoScore and to target mitochondrial metabolism to enhance the sensitivity of AML cells to currently approved therapies.

Mitochondrial metabolism supports resistance to IDH mutant inhibitors in acute myeloid leukemia

Mutations in IDH induce epigenetic and transcriptional reprogramming, differentiation bias, and susceptibility to mitochondrial inhibitors in cancer cells. Here, we first show that cell lines, PDXs, and patients with acute myeloid leukemia (AML) harboring an IDH mutation displayed an enhanced mitochondrial oxidative metabolism. Along with an increase in TCA cycle intermediates, this AML-specific metabolic behavior mechanistically occurred through the increase in electron transport chain complex I activity, mitochondrial respiration, and methylation-driven CEBPα-induced fatty acid β-oxidation of IDH1 mutant cells. While IDH1 mutant inhibitor reduced 2-HG oncometabolite and CEBPα methylation, it failed to reverse FAO and OxPHOS. These mitochondrial activities were maintained through the inhibition of Akt and enhanced activation of peroxisome proliferator-activated receptor-γ coactivator-1 PGC1α upon IDH1 mutant inhibitor. Accordingly, OxPHOS inhibitors improved anti-AML efficacy of IDH mutant inhibitors in vivo. This work provides a scientific rationale for combinatory mitochondrial-targeted therapies to treat IDH mutant AML patients, especially those unresponsive to or relapsing from IDH mutant inhibitors.

Fra-1 regulates its target genes via binding to remote enhancers without exerting major control on chromatin architecture in triple negative breast cancers

The ubiquitous family of dimeric transcription factors AP-1 is made up of Fos and Jun family proteins. It has long been thought to operate principally at gene promoters and how it controls transcription is still ill-understood. The Fos family protein Fra-1 is overexpressed in triple negative breast cancers (TNBCs) where it contributes to tumor aggressiveness. To address its transcriptional actions in TNBCs, we combined transcriptomics, ChIP-seqs, machine learning and NG Capture-C. Additionally, we studied its Fos family kin Fra-2 also expressed in TNBCs, albeit much less. Consistently with their pleiotropic effects, Fra-1 and Fra-2 up- and downregulate individually, together or redundantly many genes associated with a wide range of biological processes. Target gene regulation is principally due to binding of Fra-1 and Fra-2 at regulatory elements located distantly from cognate promoters where Fra-1 modulates the recruitment of the transcriptional co-regulator p300/CBP and where differences in AP-1 variant motif recognition can underlie preferential Fra-1- or Fra-2 bindings. Our work also shows no major role for Fra-1 in chromatin architecture control at target gene loci, but suggests collaboration between Fra-1-bound and -unbound enhancers within chromatin hubs sometimes including promoters for other Fra-1-regulated genes. Our work impacts our view of AP-1.

A New Synuclein-Transgenic Mouse Model for Early Parkinson's Reveals Molecular Features of Preclinical Disease

Understanding Parkinson's disease (PD), in particular in its earliest phases, is important for diagnosis and treatment. However, human brain samples are collected post-mortem, reflecting mainly end-stage disease. Because brain samples of mouse models can be collected at any stage of the disease process, they are useful in investigating PD progression. Here, we compare ventral midbrain transcriptomics profiles from α-synuclein transgenic mice with a progressive, early PD-like striatal neurodegeneration across different ages using pathway, gene set, and network analysis methods. Our study uncovers statistically significant altered genes across ages and between genotypes with known, suspected, or unknown function in PD pathogenesis and key pathways associated with disease progression. Among those are genotype-dependent alterations associated with synaptic plasticity and neurotransmission, as well as mitochondria-related genes and dysregulation of lipid metabolism. Age-dependent changes were among others observed in neuronal and synaptic activity, calcium homeostasis, and membrane receptor signaling pathways, many of which linked to G-protein coupled receptors. Most importantly, most changes occurred before neurodegeneration was detected in this model, which points to a sequence of gene expression events that may be relevant for disease initiation and progression. It is tempting to speculate that molecular changes similar to those changes observed in our model happen in midbrain dopaminergic neurons before they start to degenerate. In other words, we believe we have uncovered molecular changes that accompany the progression from preclinical to early PD.

Patient-derived organoids and orthotopic xenografts of primary and recurrent gliomas represent relevant patient avatars for precision oncology

Patient-based cancer models are essential tools for studying tumor biology and for the assessment of drug responses in a translational context. We report the establishment a large cohort of unique organoids and patient-derived orthotopic xenografts (PDOX) of various glioma subtypes, including gliomas with mutations in IDH1, and paired longitudinal PDOX from primary and recurrent tumors of the same patient. We show that glioma PDOXs enable long-term propagation of patient tumors and represent clinically relevant patient avatars that retain histopathological, genetic, epigenetic, and transcriptomic features of parental tumors. We find no evidence of mouse-specific clonal evolution in glioma PDOXs. Our cohort captures individual molecular genotypes for precision medicine including mutations in IDH1, ATRX, TP53, MDM2/4, amplification of EGFR, PDGFRA, MET, CDK4/6, MDM2/4, and deletion of CDKN2A/B, PTCH, and PTEN. Matched longitudinal PDOX recapitulate the limited genetic evolution of gliomas observed in patients following treatment. At the histological level, we observe increased vascularization in the rat host as compared to mice. PDOX-derived standardized glioma organoids are amenable to high-throughput drug screens that can be validated in mice. We show clinically relevant responses to temozolomide (TMZ) and to targeted treatments, such as EGFR and CDK4/6 inhibitors in (epi)genetically defined subgroups, according to MGMT promoter and EGFR/CDK status, respectively. Dianhydrogalactitol (VAL-083), a promising bifunctional alkylating agent in the current clinical trial, displayed high therapeutic efficacy, and was able to overcome TMZ resistance in glioblastoma. Our work underscores the clinical relevance of glioma organoids and PDOX models for translational research and personalized treatment studies and represents a unique publicly available resource for precision oncology.

Reference-free deconvolution, visualization and interpretation of complex DNA methylation data using DecompPipeline, MeDeCom and FactorViz

DNA methylation profiling offers unique insights into human development and diseases. Often the analysis of complex tissues and cell mixtures is the only feasible option to study methylation changes across large patient cohorts. Since DNA methylomes are highly cell type specific, deconvolution methods can be used to recover cell type-specific information in the form of latent methylation components (LMCs) from such 'bulk' samples. Reference-free deconvolution methods retrieve these components without the need for DNA methylation profiles of purified cell types. Currently no integrated and guided procedure is available for data preparation and subsequent interpretation of deconvolution results. Here, we describe a three-stage protocol for reference-free deconvolution of DNA methylation data comprising: (i) data preprocessing, confounder adjustment using independent component analysis (ICA) and feature selection using DecompPipeline, (ii) deconvolution with multiple parameters using MeDeCom, RefFreeCellMix or EDec and (iii) guided biological inference and validation of deconvolution results with the R/Shiny graphical user interface FactorViz. Our protocol simplifies the analysis and guides the initial interpretation of DNA methylation data derived from complex samples. The harmonized approach is particularly useful to dissect and evaluate cell heterogeneity in complex systems such as tumors. We apply the protocol to lung cancer methylomes from The Cancer Genome Atlas (TCGA) and show that our approach identifies the proportions of stromal cells and tumor-infiltrating immune cells, as well as associations of the detected components with clinical parameters. The protocol takes slightly >3 d to complete and requires basic R skills.

Extracellular ATP and CD39 Activate cAMP-Mediated Mitochondrial Stress Response to Promote Cytarabine Resistance in Acute Myeloid Leukemia

Relapses driven by chemoresistant leukemic cell populations are the main cause of mortality for patients with acute myeloid leukemia (AML). Here, we show that the ectonucleotidase CD39 (ENTPD1) is upregulated in cytarabine-resistant leukemic cells from both AML cell lines and patient samples in vivo and in vitro. CD39 cell-surface expression and activity is increased in patients with AML upon chemotherapy compared with diagnosis, and enrichment in CD39-expressing blasts is a marker of adverse prognosis in the clinics. High CD39 activity promotes cytarabine resistance by enhancing mitochondrial activity and biogenesis through activation of a cAMP-mediated adaptive mitochondrial stress response. Finally, genetic and pharmacologic inhibition of CD39 ecto-ATPase activity blocks the mitochondrial reprogramming triggered by cytarabine treatment and markedly enhances its cytotoxicity in AML cells in vitro and in vivo. Together, these results reveal CD39 as a new residual disease marker and a promising therapeutic target to improve chemotherapy response in AML. SIGNIFICANCE: Extracellular ATP and CD39-P2RY13-cAMP-OxPHOS axis are key regulators of cytarabine resistance, offering a new promising therapeutic strategy in AML.This article is highlighted in the In This Issue feature, p. 1426.

Machine learning-assisted neurotoxicity prediction in human midbrain organoids

INTRODUCTION

Brain organoids are highly complex multi-cellular tissue proxies, which have recently risen as novel tools to study neurodegenerative diseases such as Parkinson's disease (PD). However, with increasing complexity of the system, usage of quantitative tools becomes challenging.

OBJECTIVES

The primary objective of this study was to develop a neurotoxin-induced PD organoid model and to assess the neurotoxic effect on dopaminergic neurons using microscopy-based phenotyping in a high-content fashion.

METHODS

We describe a pipeline for a machine learning-based analytical method, allowing for detailed image-based cell profiling and toxicity prediction in brain organoids treated with the neurotoxic compound 6-hydroxydopamine (6-OHDA).

RESULTS

We quantified features such as dopaminergic neuron count and neuronal complexity and built a machine learning classifier with the data to optimize data processing strategies and to discriminate between different treatment conditions. We validated the approach with high content imaging data from PD patient derived midbrain organoids.

CONCLUSIONS

The here described model is a valuable tool for advanced in vitro PD modeling and to test putative neurotoxic compounds.

GENE-02. ESTABLISHING PERSONALIZED TREATMENT OPTIONS FOR RECURRENT HIGH-GRADE GLIOMAS

BACKGROUND

High grade glioma (HGG) patients develop resistance to standard treatment leading to disease progression and limited life expectancy. Recent advances in the molecular characterisation of treatment-naïve HGGs based on next generation sequencing and DNA methylation analyses have led to a better delineation of HGG-subtypes and identification of distinct genomic abnormalities opening opportunities for personalized treatment strategies.

METHODS

We collected 300 fresh glioma specimen with approximately 100 longitudinal samples of initial and recurrent tumors from 43 matched patients. We succeeded in generating a live-biobank of HGG patient-derived orthotopic xenografts (PDOX) and 3D tumor organoids that neatly recapitulates the mutational spectrum including structural DNA variations and methylation-based subtypes of gliomas. A highlight is the generation of 19 PDOXs of paired initial and relapse HGGs from 9 glioma patients, enabling high-throughput drug screens. We performed comprehensive molecular profiling using arrayCGH, DNA-methylation and targeted DNA sequencing on patient specimen and their derivatives, 3D tumor organoids and PDOXs.

RESULTS

Detailed analysis of the paired longitudinal samples indicated that PDOXs closely recapitulate the evolutionary trajectory of the parental tumors. Furthermore, targeted genomic sequencing of paired HGGs suggests that relapse tumors also accumulate somatic mutations in epigenetic effectors. Based on patient-derived material we carried out drug response screening on 3D tumor organoids using a compound library matching the majority of genes that were assessed with targeted sequencing. Differential drug responses between initial and recurrent tumors were observed and the prevailing primary drug response profiles were essentially recapitulad in the relapse setting.

CONCLUSIONS

Response assessment of treatment-naïve gliomas and their recurrences provides crucial information on the differential sensitivity between initial and relapsed HGGs and offers novel personalized therapeutic options for the relapse setting. Furthermore, in depth correlation of the profiled somatic molecular landscape with drug response will enable pharmacogenomic predictions of potential inhibitors in the clinical setting.

OS12.2 Targeting epigenetic pathways in the treatment of recurrent high-grade glioma

BACKGROUND

High grade glioma (HGG) patients develop resistance to standard treatment leading to disease progression and limited life expectancy. Advances in the molecular characterisation of treatment-naïve HGGs based on next-generation sequencing and DNA methylation analyses have led to a better delineation of HGG subtypes and the identification of distinct genomic abnormalities. Furthermore, using large patient cohorts of longitudinal tumor samples, comprehensive genomic profiling studies emerged to investigate therapy-associated evolution of gliomas. All together, those studies point out the need for personalised treatment strategies, where applied drugs will be adapted to the unique patient-specific genetic abnormalities.

MATERIAL AND METHODS

We collected fresh samples of more than 800 brain tumors containing almost 300 glioma specimen with approximately 100 longitudinal samples of initial and recurrent tumors from 43 matched patients. By now, we have successfully established 34 patient-derived orthotopic xenografts (PDOXs) in mice. We performed comprehensive molecular profiling using array comparative genomic hybridisation, DNA methylation analysis and targeted DNA sequencing on patient specimen and their derivatives such as 3D tumor organoids and PDOXs. The custom-design sequencing panel comprises 234 genes that reflect both established genetic identifiers for individual glioma subtype classification and novel genes encoding mainly epigenetic effector genes. Based on patient-derived material we carried out drug response screening on 3D tumor organoids using a compound library matching the majority of genes that were assessed by targeted sequencing.

RESULTS

We succeeded in generating a live biobank of HGG patient-derived xenografts and 3D organoids that neatly recapitulates the mutational spectrum including structural DNA variation and methylation-based subtypes of gliomas. A highlight is the generation of 19 PDOXs of paired initial and relapse HGGs from a total of 9 glioma patients. A detailed analysis of the paired longitudinal samples indicated that PDOX models closely recapitulate the evolutionary trajectory of the parental tumors. Targeted sequencing of longitudinal HGG PDOXs suggests that relapse tumors accumulate somatic mutations in epigenetic effectors compared with the Initial. Differential drug responses between initial and relapse tumors were observed after screening of in vitro 3D tumor organoids.

CONCLUSION

Response assessment of naïve initial gliomas and recurrences provides crucial information on the differential sensitivity between initial and relapsed HGGs and offers novel personalised therapeutic options in the relapse setting. Furthermore, in depth correlation of the profiled somatic molecular landscape with drug response will enable pharmacogenomic predictions of potential inhibitors in the clinical setting.

Distribution and prognostic impact of microglia/macrophage subpopulations in gliomas

While the central nervous system is considered an immunoprivileged site and brain tumors display immunosuppressive features, both innate and adaptive immune responses affect glioblastoma (GBM) growth and treatment resistance. However, the impact of the major immune cell population in gliomas, represented by glioma‐associated microglia/macrophages (GAMs), on patients’ clinical course is still unclear. Thus, we aimed at assessing the immunohistochemical expression of selected microglia and macrophage markers in 344 gliomas (including gliomas from WHO grade I–IV). Furthermore, we analyzed a cohort of 241 IDH1R132H‐non‐mutant GBM patients for association of GAM subtypes and patient overall survival. Phenotypical properties of GAMs, isolated from high‐grade astrocytomas by CD11b‐based magnetic cell sorting, were analyzed by immunocytochemistry, mRNA microarray, qRT‐PCR and bioinformatic analyses. A higher amount of CD68‐, CD163‐ and CD206‐positive GAMs in the vital tumor core was associated with beneficial patient survival. The mRNA expression profile of GAMs displayed an upregulation of factors that are considered as pro‐inflammatory M1 (eg, CCL2, CCL3L3, CCL4, PTGS2) and anti‐inflammatory M2 polarization markers (eg, MRC1, LGMN, CD163, IL10, MSR1), the latter rather being associated with phagocytic functions in the GBM microenvironment. In summary, we present evidence that human GBMs contain mixed M1/M2‐like polarized GAMs and that the levels of different GAM subpopulations in the tumor core are positively associated with overall survival of patients with IDH1R132H‐non‐mutant GBMs.

Hub genes in a pan-cancer co-expression network show potential for predicting drug responses

Background: The topological analysis of networks extracted from different types of "omics" data is a useful strategy for characterizing biologically meaningful properties of the complex systems underlying these networks. In particular, the biological significance of highly connected genes in diverse molecular networks has been previously determined using data from several model organisms and phenotypes. Despite such insights, the predictive potential of candidate hubs in gene co-expression networks in the specific context of cancer-related drug experiments remains to be deeply investigated. The examination of such associations may offer opportunities for the accurate prediction of anticancer drug responses.  Methods: Here, we address this problem by: a) analyzing a co-expression network obtained from thousands of cancer cell lines, b) detecting significant network hubs, and c) assessing their capacity to predict drug sensitivity using data from thousands of drug experiments. We investigated the prediction capability of those genes using a multiple linear regression model, independent datasets, comparisons with other models and our own in vitro experiments. Results: These analyses led to the identification of 47 hub genes, which are implicated in a diverse range of cancer-relevant processes and pathways. Overall, encouraging agreements between predicted and observed drug sensitivities were observed in public datasets, as well as in our in vitro validations for four glioblastoma cell lines and four drugs. To facilitate further research, we share our hub-based drug sensitivity prediction model as an online tool. Conclusions: Our research shows that co-expression network hubs are biologically interesting and exhibit potential for predicting drug responses in vitro. These findings motivate further investigations about the relevance and application of our unbiased discovery approach in pre-clinical, translationally-oriented research.

Hub genes in a pan-cancer co-expression network show potential for predicting drug responses

Background: The topological analysis of networks extracted from different types of "omics" data is a useful strategy for characterizing biologically meaningful properties of the complex systems underlying these networks. In particular, the biological significance of highly connected genes in diverse molecular networks has been previously determined using data from several model organisms and phenotypes. Despite such insights, the predictive potential of candidate hubs in gene co-expression networks in the specific context of cancer-related drug experiments remains to be deeply investigated. The examination of such associations may offer opportunities for the accurate prediction of anticancer drug responses.  Methods: Here, we address this problem by: a) analyzing a co-expression network obtained from thousands of cancer cell lines, b) detecting significant network hubs, and c) assessing their capacity to predict drug sensitivity using data from thousands of drug experiments. We investigated the prediction capability of those genes using a multiple linear regression model, independent datasets, comparisons with other models and our own in vitro experiments. Results: These analyses led to the identification of 47 hub genes, which are implicated in a diverse range of cancer-relevant processes and pathways. Overall, encouraging agreements between predicted and observed drug sensitivities were observed in public datasets, as well as in our in vitro validations for four glioblastoma cell lines and four drugs. To facilitate further research, we share our hub-based drug sensitivity prediction model as an online tool. Conclusions: Our research shows that co-expression network hubs are biologically interesting and exhibit potential for predicting drug responses in vitro. These findings motivate further investigations about the relevance and application of our unbiased discovery approach in pre-clinical, translationally-oriented research.

Single-cell transcriptomics reveals distinct inflammation-induced microglia signatures

Microglia are specialized parenchymal‐resident phagocytes of the central nervous system (CNS) that actively support, defend and modulate the neural environment. Dysfunctional microglial responses are thought to worsen CNS diseases; nevertheless, their impact during neuroinflammatory processes remains largely obscure. Here, using a combination of single‐cell RNA sequencing and multicolour flow cytometry, we comprehensively profile microglia in the brain of lipopolysaccharide (LPS)‐injected mice. By excluding the contribution of other immune CNS‐resident and peripheral cells, we show that microglia isolated from LPS‐injected mice display a global downregulation of their homeostatic signature together with an upregulation of inflammatory genes. Notably, we identify distinct microglial activated profiles under inflammatory conditions, which greatly differ from neurodegenerative disease‐associated profiles. These results provide insights into microglial heterogeneity and establish a resource for the identification of specific phenotypes in CNS disorders, such as neuroinflammatory and neurodegenerative diseases.

Predicting clinical outcome of neuroblastoma patients using an integrative network-based approach

BACKGROUND

One of the main current challenges in computational biology is to make sense of the huge amounts of multidimensional experimental data that are being produced. For instance, large cohorts of patients are often screened using different high-throughput technologies, effectively producing multiple patient-specific molecular profiles for hundreds or thousands of patients.

RESULTS

We propose and implement a network-based method that integrates such patient omics data into Patient Similarity Networks. Topological features derived from these networks were then used to predict relevant clinical features. As part of the 2017 CAMDA challenge, we have successfully applied this strategy to a neuroblastoma dataset, consisting of genomic and transcriptomic data. In particular, we observe that models built on our network-based approach perform at least as well as state of the art models. We furthermore explore the effectiveness of various topological features and observe, for instance, that redundant centrality metrics can be combined to build more powerful models.

CONCLUSION

We demonstrate that the networks inferred from omics data contain clinically relevant information and that patient clinical outcomes can be predicted using only network topological data.

REVIEWERS

This article was reviewed by Yang-Yu Liu, Tomislav Smuc and Isabel Nepomuceno.

Inhibition of HIF1α-Dependent Upregulation of Phospho-l-Plastin Resensitizes Multiple Myeloma Cells to Frontline Therapy

The introduction of novel frontline agents in multiple myeloma (MM), like immunomodulatory drugs and proteasome inhibitors, has improved the overall survival of patients. Yet, MM is still not curable, and drug resistance (DR) remains the main challenge. To improve the understanding of DR in MM, we established a resistant cell line (MOLP8/R). The exploration of DR mechanisms yielded an overexpression of HIF1α, due to impaired proteasome activity of MOLP8/R. We show that MOLP8/R, like other tumor cells, overexpressing HIF1α, have an increased resistance to the immune system. By exploring the main target genes regulated by HIF1α, we could not show an overexpression of these targets in MOLP8/R. We, however, show that MOLP8/R cells display a very high overexpression of LCP1 gene (l-Plastin) controlled by HIF1α, and that this overexpression also exists in MM patient samples. The l-Plastin activity is controlled by its phosphorylation in Ser5. We further show that the inhibition of l-Plastin phosphorylation restores the sensitivity of MOLP8/R to immunomodulatory drugs (IMiDs) and proteasome inhibitors (PIs). Our results reveal a new target gene of DR, controlled by HIF1α.

Cytoplasmic overexpression of RNA-binding protein HuR is a marker of poor prognosis in meningioma, and HuR knockdown decreases meningioma cell growth and resistance to hypoxia

HuR regulates cytoplasmic mRNA stability and translatability, and the HuR expression level has been shown to correlate with poor disease outcome in several cancer types; however, the prognostic value and potential pro-oncogenic properties of HuR in meningioma remain unclear. Thus, in the present study, we analysed 85 meningioma tissue samples to establish the relationship between HuR expression, tumour cell proliferation, and/or patient survival. In addition, we examined the anti-proliferative effects of HuR knockdown in two meningioma cell lines (IOMM-Lee and Ben-Men-1) and conducted transcriptome-wide analyses (IOMM-Lee cells) to elucidate the molecular consequences of HuR knockdown. The results of the present study showed HuR cytoplasmic expression to correlate positively with tumour grade (p = 1.2 × 10^-8 ) and negatively with progression-free and overall survival (p = 0.01) time in human meningioma tissues. In vitro, siHuR-induced HuR knockdown was shown to reduce the growth of both Ben-Men-1 (p = 2 × 10^-8 ) and IOMM-Lee (p = 4 × 10^-9 ) cells. Transcriptome analyses revealed HuR knockdown in IOMM-Lee cells to deregulate the HIF1A signalling pathway (p = 1.5 × 10^-6 ) and to up-regulate the expression of genes essential for the assembly of the cytoplasmic mRNA processing body, global genome nucleotide-excision repair, poly(A)-specific ribonuclease activity, the positive regulation of apoptosis and of cell cycle arrest, and the negative regulation of RNA splicing [p(FDR) < 0.001]. Interestingly, HuR knockdown under hypoxic culture conditions further potentiated the effects of HuR knockdown on cell growth, apoptosis, and HIF1A expression. We thus conclude that cytoplasmic HuR expression is a marker of poor prognosis in meningioma and that HuR is a promising potential therapeutic target for use in tumours refractory to standard therapies. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Chemotherapy-Resistant Human Acute Myeloid Leukemia Cells Are Not Enriched for Leukemic Stem Cells but Require Oxidative Metabolism

Chemotherapy-resistant human acute myeloid leukemia (AML) cells are thought to be enriched in quiescent immature leukemic stem cells (LSC). To validate this hypothesis in vivo, we developed a clinically relevant chemotherapeutic approach treating patient-derived xenografts (PDX) with cytarabine (AraC). AraC residual AML cells are enriched in neither immature, quiescent cells nor LSCs. Strikingly, AraC-resistant preexisting and persisting cells displayed high levels of reactive oxygen species, showed increased mitochondrial mass, and retained active polarized mitochondria, consistent with a high oxidative phosphorylation (OXPHOS) status. AraC residual cells exhibited increased fatty-acid oxidation, upregulated CD36 expression, and a high OXPHOS gene signature predictive for treatment response in PDX and patients with AML. High OXPHOS but not low OXPHOS human AML cell lines were chemoresistant in vivo. Targeting mitochondrial protein synthesis, electron transfer, or fatty-acid oxidation induced an energetic shift toward low OXPHOS and markedly enhanced antileukemic effects of AraC. Together, this study demonstrates that essential mitochondrial functions contribute to AraC resistance in AML and are a robust hallmark of AraC sensitivity and a promising therapeutic avenue to treat AML residual disease.Significance: AraC-resistant AML cells exhibit metabolic features and gene signatures consistent with a high OXPHOS status. In these cells, targeting mitochondrial metabolism through the CD36-FAO-OXPHOS axis induces an energetic shift toward low OXPHOS and strongly enhanced antileukemic effects of AraC, offering a promising avenue to design new therapeutic strategies and fight AraC resistance in AML. Cancer Discov; 7(7); 716-35. ©2017 AACR.See related commentary by Schimmer, p. 670This article is highlighted in the In This Issue feature, p. 653.

RNA sequencing and transcriptome arrays analyses show opposing results for alternative splicing in patient derived samples

Background

RNA sequencing (RNA-seq) and microarrays are two transcriptomics techniques aimed at the quantification of transcribed genes and their isoforms. Here we compare the latest Affymetrix HTA 2.0 microarray with Illumina 2000 RNA-seq for the analysis of patient samples - normal lung epithelium tissue and squamous cell carcinoma lung tumours. Protein coding mRNAs and long non-coding RNAs (lncRNAs) were included in the study.

Results

Both platforms performed equally well for protein-coding RNAs, however the stochastic variability was higher for the sequencing data than for microarrays. This reduced the number of differentially expressed genes and genes with predictive potential for RNA-seq compared to microarray data. Analysis of this variability revealed a lack of reads for short and low abundant genes; lncRNAs, being shorter and less abundant RNAs, were found especially susceptible to this issue. A major difference between the two platforms was uncovered by analysis of alternatively spliced genes. Investigation of differential exon abundance showed insufficient reads for many exons and exon junctions in RNA-seq while the detection on the array platform was more stable. Nevertheless, we identified 207 genes which undergo alternative splicing and were consistently detected by both techniques.

Conclusions

Despite the fact that the results of gene expression analysis were highly consistent between Human Transcriptome Arrays and RNA-seq platforms, the analysis of alternative splicing produced discordant results. We concluded that modern microarrays can still outperform sequencing for standard analysis of gene expression in terms of reproducibility and cost.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3819-y) contains supplementary material, which is available to authorized users.

Isocitrate dehydrogenase 1 mutations prime the all-trans retinoic acid myeloid differentiation pathway in acute myeloid leukemia

Boutzen et al. show that the IDH1 mutation and its oncometabolite, (R)-2-hydroxyglutarate, dysregulate downstream target pathways of myeloid-specific TFs, especially CEBPα, priming mutant IDH1-R132H AML blasts to the granulomonocytic lineage.

Acute myeloid leukemia (AML) is characterized by the accumulation of malignant blasts with impaired differentiation programs caused by recurrent mutations, such as the isocitrate dehydrogenase (IDH) mutations found in 15% of AML patients. These mutations result in the production of the oncometabolite (R)-2-hydroxyglutarate (2-HG), leading to a hypermethylation phenotype that dysregulates hematopoietic differentiation. In this study, we identified mutant R132H IDH1-specific gene signatures regulated by key transcription factors, particularly CEBPα, involved in myeloid differentiation and retinoid responsiveness. We show that treatment with all-trans retinoic acid (ATRA) at clinically achievable doses markedly enhanced terminal granulocytic differentiation in AML cell lines, primary patient samples, and a xenograft mouse model carrying mutant IDH1. Moreover, treatment with a cell-permeable form of 2-HG sensitized wild-type IDH1 AML cells to ATRA-induced myeloid differentiation, whereas inhibition of 2-HG production significantly reduced ATRA effects in mutant IDH1 cells. ATRA treatment specifically decreased cell viability and induced apoptosis of mutant IDH1 blasts in vitro. ATRA also reduced tumor burden of mutant IDH1 AML cells xenografted in NOD–Scid–IL2rγ^null mice and markedly increased overall survival, revealing a potent antileukemic effect of ATRA in the presence of IDH1 mutation. This therapeutic strategy holds promise for this AML patient subgroup in future clinical studies.

Assessing cellular and circulating miRNA recovery: the impact of the RNA isolation method and the quantity of input material

MicroRNAs (miRNAs) have emerged as promising cancer biomarkers. However, exploiting their informative potential requires careful optimization of their detection. Here, we compared the efficiency of commonly used RNA extraction kits in miRNA recovery from cells, plasma and urine/plasma-derived exosomes, using single-gene RT-qPCR and miRNA profiling. We used increasing amounts of starting material to investigate the impact of the input material size on miRNA extraction. We showed that miRNA recovery was largely influenced by the isolation method and by the amount of input material. In particular, the miRCURY™ kit provided highly pure RNA. However, its columns poorly recovered miRNAs from limiting amounts of cells and plasma, and rapidly saturated by large RNA species and plasma components, thus impeding miRNA recovery from high input amounts. Overall, the miRNeasy® kit permitted a better miRNA detection despite a less pure extracted RNA. Nevertheless, some miRNAs were preferentially or exclusively isolated by either of the methods. Trizol® LS resulted in very low purity RNA which affected RT-qPCR efficiency. In general, miRCURY™ biofluids kit efficiently extracted miRNAs from plasma. A careful selection of the RNA isolation method and the consideration of the type and size of input material are highly recommended to avoid biased results.

L-plastin Ser5 phosphorylation in breast cancer cells and in vitro is mediated by RSK downstream of the ERK/MAPK pathway

Deregulated cell migration and invasion are hallmarks of metastatic cancer cells. Phosphorylation on residue Ser5 of the actin-bundling protein L-plastin activates L-plastin and has been reported to be crucial for invasion and metastasis. Here, we investigate signal transduction leading to L-plastin Ser5 phosphorylation using 4 human breast cancer cell lines. Whole-genome microarray analysis comparing cell lines with different invasive capacities and corresponding variations in L-plastin Ser5 phosphorylation level revealed that genes of the ERK/MAPK pathway are differentially expressed. It is noteworthy that in vitro kinase assays showed that ERK/MAPK pathway downstream ribosomal protein S6 kinases α-1 (RSK1) and α-3 (RSK2) are able to directly phosphorylate L-plastin on Ser5. Small interfering RNA- or short hairpin RNA-mediated knockdown and activation/inhibition studies followed by immunoblot analysis and computational modeling confirmed that ribosomal S6 kinase (RSK) is an essential activator of L-plastin. Migration and invasion assays showed that RSK knockdown led to a decrease of up to 30% of migration and invasion of MDA-MB-435S cells. Although the presence of L-plastin was not necessary for migration/invasion of these cells, immunofluorescence assays illustrated RSK-dependent recruitment of Ser5-phosphorylated L-plastin to migratory structures. Altogether, we provide evidence that the ERK/MAPK pathway is involved in L-plastin Ser5 phosphorylation in breast cancer cells with RSK1 and RSK2 kinases able to directly phosphorylate L-plastin residue Ser5.

Integrated metabolic modelling reveals cell-type specific epigenetic control points of the macrophage metabolic network

BACKGROUND

The reconstruction of context-specific metabolic models from easily and reliably measurable features such as transcriptomics data will be increasingly important in research and medicine. Current reconstruction methods suffer from high computational effort and arbitrary threshold setting. Moreover, understanding the underlying epigenetic regulation might allow the identification of putative intervention points within metabolic networks. Genes under high regulatory load from multiple enhancers or super-enhancers are known key genes for disease and cell identity. However, their role in regulation of metabolism and their placement within the metabolic networks has not been studied.

METHODS

Here we present FASTCORMICS, a fast and robust workflow for the creation of high-quality metabolic models from transcriptomics data. FASTCORMICS is devoid of arbitrary parameter settings and due to its low computational demand allows cross-validation assays. Applying FASTCORMICS, we have generated models for 63 primary human cell types from microarray data, revealing significant differences in their metabolic networks.

RESULTS

To understand the cell type-specific regulation of the alternative metabolic pathways we built multiple models during differentiation of primary human monocytes to macrophages and performed ChIP-Seq experiments for histone H3 K27 acetylation (H3K27ac) to map the active enhancers in macrophages. Focusing on the metabolic genes under high regulatory load from multiple enhancers or super-enhancers, we found these genes to show the most cell type-restricted and abundant expression profiles within their respective pathways. Importantly, the high regulatory load genes are associated to reactions enriched for transport reactions and other pathway entry points, suggesting that they are critical regulatory control points for cell type-specific metabolism.

CONCLUSIONS

By integrating metabolic modelling and epigenomic analysis we have identified high regulatory load as a common feature of metabolic genes at pathway entry points such as transporters within the macrophage metabolic network. Analysis of these control points through further integration of metabolic and gene regulatory networks in various contexts could be beneficial in multiple fields from identification of disease intervention strategies to cellular reprogramming.

Comprehensive analysis of glycolytic enzymes as therapeutic targets in the treatment of glioblastoma

Major efforts have been put in anti-angiogenic treatment for glioblastoma (GBM), an aggressive and highly vascularized brain tumor with dismal prognosis. However clinical outcome with anti-angiogenic agents has been disappointing and tumors quickly develop escape mechanisms. In preclinical GBM models we have recently shown that bevacizumab, a blocking antibody against vascular endothelial growth factor, induces hypoxia in treated tumors, which is accompanied by increased glycolytic activity and tumor invasiveness. Genome-wide transcriptomic analysis of patient derived GBM cells including stem cell lines revealed a strong up-regulation of glycolysis-related genes in response to severe hypoxia. We therefore investigated the importance of glycolytic enzymes in GBM adaptation and survival under hypoxia, both in vitro and in vivo. We found that shRNA-mediated attenuation of glycolytic enzyme expression interfered with GBM growth under normoxic and hypoxic conditions in all cellular models. Using intracranial GBM xenografts we identified seven glycolytic genes whose knockdown led to a dramatic survival benefit in mice. The most drastic effect was observed for PFKP (PFK1, +21.8%) and PDK1 (+20.9%), followed by PGAM1 and ENO1 (+14.5% each), HK2 (+11.8%), ALDOA (+10.9%) and ENO2 (+7.2%). The increase in mouse survival after genetic interference was confirmed using chemical inhibition of PFK1 with clotrimazole. We thus provide a comprehensive analysis on the importance of the glycolytic pathway for GBM growth in vivo and propose PFK1 and PDK1 as the most promising therapeutic targets to address the metabolic escape mechanisms of GBM.

An essential role of syntaxin 3 protein for granule exocytosis and secretion of IL-1α, IL-1β, IL-12b, and CCL4 from differentiated HL-60 cells

Besides their roles in the killing of pathogens, neutrophils have the capacity to package a variety of cytokines into cytoplasmic granules for subsequent release upon inflammatory conditions. Because the rapid secretion of cytokines orchestrates the action of other immune cells at the infection site and thus, can contribute to the development and chronicity of inflammatory diseases, we aimed to determine the intracellular SNARE machinery responsible for the regulation of cytokine secretion and degranulation. From a constructed gene-expression network, we first selected relevant cytokines for functional validation by the CBA approach. We established a cytokine-secretion profile for human neutrophils and dHL-60 cells, underlining their similar ability to secrete a broad variety of cytokines within proinflammatory conditions mimicked by LPS stimulation. Secondly, after screening of SNARE genes by microarray experiments, we selected STX3 for further functional studies. With the use of a siRNA strategy, we show that STX3 is clearly required for the maximal release of IL-1α, IL-1β, IL-12b, and CCL4 without alteration of other cytokine secretion in dHL-60 cells. In addition, we demonstrate that STX3 is involved in MMP-9 exocytosis from gelatinase granules, where STX3 is partly localized. Our results suggest that the secretion of IL-1α, IL-1β, IL-12b, and CCL4 occurs during gelatinase degranulation, a process controlled by STX3. In summary, these findings provide first evidence that STX3 has an essential role in trafficking pathways of cytokines in neutrophil granulocytes.

NK cell killer Ig-like receptor repertoire acquisition and maturation are strongly modulated by HLA class I molecules

The interaction between clonally distributed inhibitory receptors and their activating counterparts on NK cells and HLA class I molecules defines NK cell functions, but the role of HLA class I ligands in the acquisition of their receptors during NK development is still unclear. Although some studies demonstrated that HLA-C affects the expression of killer Ig-like receptors (KIR), other studies showed that NK cells acquire their KIR repertoire in a stochastic manner. Only when infected with human CMV is an expansion of self-specific KIR(+) NKG2C(+) NK cells detected. To gain more insight into this question, we compared the coexpression of different KIR molecules, NKG2A, CD8, and CD57, on NK cells in healthy donors and seven patients with deficient HLA class I expression due to mutations in one of the TAP genes. Our results show a correlation between the presence/absence of HLA class I molecules and the coexpression of their receptors. In an HLA class I low-expression context, an increase in KIR molecules' coexpression is detected on the NKG2A(+) CD8(+) subset. In functional assays, hyporesponsiveness was observed for TAP-deficient NK cells derived from four patients. In contrast, NK cells from patient five were functional, whereas CD107a(+) and IFN-γ(+) CD56(dim) NK cells presented a different pattern of HLA class I receptors compared with healthy donors. Taken together, our results provide strong evidence for the role of HLA class I molecules in NK cell maturation and KIR repertoire acquisition.

Abstract 3905: Identification of new tumor-specific splice variants in NSCL cancers by whole genome exon arrays

Although lung cancer is one of the major causes of cancer death in men and women worldwide, the molecular pathogenesis of this disease remains elusive. The contribution of alternative splicing (AS) to cancer pathogenesis and progression is emerging as an area of considerable interest. Cancer-specific splice variants seem to play a key role in disease mechanism and etiology.

In this study, we analyzed 20 matched pairs of tumor specimens of Adenocarcinomas (AC), Squamous Cell Carcinomas (SCC) and adjacent normal tissue using Affymetrix GeneChip Exon microarrays. With this technology, it is possible to profile simultaneously gene expression and alternative splicing patterns for the entire genome using a single array. At the exon level, multiple probes for each exon enable the discrimination among different isoforms of a gene. Total RNA from tumors or normal tissue samples were amplified using the Ambion whole transcript expression kit and were hybridized onto GeneChip Human Exon 1.0 ST Array®. Data generated were normalized using RMA and were subjected to several layers of filtering (DABG filter, Multiple mRNAs filter). The delta splice index (ΔSI) was calculated, and alternative splicing events were detected using MIDAS and Rank Product method. Probesets were ranked based on their SI p-values.

By comparing tumor samples to the corresponding normal subgroups, we generated a list of 1689 splice variant candidates for SCC and 1010 for AC showing significant differential expression between normal and tumor tissues. Top candidates in each tumor subtype were selected based on their ΔSI for manual inspection and further bioinformatics analysis. Study of alternatively spliced candidates revealed that approximately 30 % of the genes were detected in both types of cancer. Gene function analysis combining text mining and knowledge-based approaches indicated that the largest subset (70 %) of these genes was related to cancer. Within this group, 30% of genes were related to lung cancers. Interestingly the major part (85%) of variant candidates was not previously associated with AS events in cancer. We focused our attention on this latter group and, through validation by RT/PCR, confirmed the existence of new specific alternative splicing events occurring in each cancer type. In addition, we identified biological pathways and gene networks that are specifically altered in AC or SCC. The cancer-specific AS variants and pathways identified in our study are promising biomarkers and targets for diagnostic or therapeutic purposes.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 3905. doi:10.1158/1538-7445.AM2011-3905