Despite advances in the molecular exploration of paediatric cancers, approximately 50% of children with high-risk neuroblastoma lack effective treatment. To identify therapeutic options for this group of high-risk patients, we combine predictive data mining with experimental evaluation in patient-derived xenograft cells. Our proposed algorithm, TargetTranslator, integrates data from tumour biobanks, pharmacological databases, and cellular networks to predict how targeted interventions affect mRNA signatures associated with high patient risk or disease processes. We find more than 80 targets to be associated with neuroblastoma risk and differentiation signatures. Selected targets are evaluated in cell lines derived from high-risk patients to demonstrate reversal of risk signatures and malignant phenotypes. Using neuroblastoma xenograft models, we establish CNR2 and MAPK8 as promising candidates for the treatment of high-risk neuroblastoma. We expect that our method, available as a public tool (targettranslator.org), will enhance and expedite the discovery of risk-associated targets for paediatric and adult cancers.

Tumor cell heterogeneity is a crucial characteristic of malignant brain tumors and underpins phenomena such as therapy resistance and tumor recurrence. Advances in single‐cell analysis have enabled the delineation of distinct cellular states of brain tumor cells, but the time‐dependent changes in such states remain poorly understood. Here, we construct quantitative models of the time‐dependent transcriptional variation of patient‐derived glioblastoma (GBM) cells. We build the models by sampling and profiling barcoded GBM cells and their progeny over the course of 3 weeks and by fitting a mathematical model to estimate changes in GBM cell states and their growth rates. Our model suggests a hierarchical yet plastic organization of GBM, where the rates and patterns of cell state switching are partly patient‐specific. Therapeutic interventions produce complex dynamic effects, including inhibition of specific states and altered differentiation. Our method provides a general strategy to uncover time‐dependent changes in cancer cells and offers a way to evaluate and predict how therapy affects cell state composition.

Glioblastoma (GBM) is an aggressive type of brain cancer with a poor prognosis for affected patients. The current line of treatment only gives the patients a survival time of on average 15 months. In this work, we use genome-scale metabolic models (GEMs) together with other systems biology tools to examine the global transcriptomics-data of GBM-patients obtained from The Cancer Genome Atlas (TCGA). We reveal the molecular mechanisms underlying GBM and identify potential therapeutic targets for effective treatment of patients. The work presented consists of two main parts. The first part stratifies the patients into two groups, high and low survival, and compares their gene expression. The second part uses GBM and healthy brain tissue GEMs to simulate gene knockout in a GBM cell model to find potential therapeutic targets and predict their side effect in healthy brain tissue. We (1) find that genes upregulated in the patients with low survival are linked to various stages of the glioma invasion process, and (2) identify five essential genes for GBM, whose inhibition is non-toxic to healthy brain tissue, therefore promising to investigate further as therapeutic targets.