Lineage specification in glioblastoma is regulated by METTL7B.

Glioblastomas are the most common malignant brain tumors in adults; they are highly aggressive and heterogeneous and show a high degree of plasticity. Here, we show that methyltransferase-like 7B (METTL7B) is an essential regulator of lineage specification in glioblastoma, with an impact on both tumor size and invasiveness. Single-cell transcriptomic analysis of these tumors and of cerebral organoids derived from expanded potential stem cells overexpressing METTL7B reveal a regulatory role for the gene in the neural stem cell-to-astrocyte differentiation trajectory. Mechanistically, METTL7B downregulates the expression of key neuronal differentiation players, including SALL2, via post-translational modifications of histone marks.

Correction: Single-cell RNAseq and longitudinal proteomic analysis of a novel semi-spontaneous urothelial cancer model reveals tumor cell heterogeneity and pretumoral urine protein alterations.

[This corrects the article DOI: 10.1371/journal.pone.0253178.].

Live Detection of Neural Progenitors and Glioblastoma Cells by an Oligothiophene Derivative.

There is an urgent need for simple and non-invasive identification of live neural stem/progenitor cells (NSPCs) in the developing and adult brain as well as in disease, such as in brain tumors, due to the potential clinical importance in prognosis, diagnosis, and treatment of diseases of the nervous system. Here, we report a luminescent conjugated oligothiophene (LCO), named p-HTMI, for non-invasive and non-amplified real-time detection of live human patient-derived glioblastoma (GBM) stem cell-like cells and NSPCs. While p-HTMI stained only a small fraction of other cell types investigated, the mere addition of p-HTMI to the cell culture resulted in efficient detection of NSPCs or GBM cells from rodents and humans within minutes. p-HTMI is functionalized with a methylated imidazole moiety resembling the side chain of histidine/histamine, and non-methylated analogues were not functional. Cell sorting experiments of human GBM cells demonstrated that p-HTMI labeled the same cell population as CD271, a proposed marker for stem cell-like cells and rapidly migrating cells in glioblastoma. Our results suggest that the LCO p-HTMI is a versatile tool for immediate and selective detection of neural and glioma stem and progenitor cells.

Inference of glioblastoma migration and proliferation rates using single time-point images.

Cancer cell migration is a driving mechanism of invasion in solid malignant tumors. Anti-migratory treatments provide an alternative approach for managing disease progression. However, we currently lack scalable screening methods for identifying novel anti-migratory drugs. To this end, we develop a method that can estimate cell motility from single end-point images in vitro by estimating differences in the spatial distribution of cells and inferring proliferation and diffusion parameters using agent-based modeling and approximate Bayesian computation. To test the power of our method, we use it to investigate drug responses in a collection of 41 patient-derived glioblastoma cell cultures, identifying migration-associated pathways and drugs with potent anti-migratory effects. We validate our method and result in both in silico and in vitro using time-lapse imaging. Our proposed method applies to standard drug screen experiments, with no change needed, and emerges as a scalable approach to screen for anti-migratory drugs.

Global hypo-methylation in a proportion of glioblastoma enriched for an astrocytic signature is associated with increased invasion and altered immune landscape.

We describe a subset of glioblastoma, the most prevalent malignant adult brain tumour, harbouring a bias towards hypomethylation at defined differentially methylated regions. This epigenetic signature correlates with an enrichment for an astrocytic gene signature, which together with the identification of enriched predicted binding sites of transcription factors known to cause demethylation and to be involved in astrocytic/glial lineage specification, point to a shared ontogeny between these glioblastomas and astroglial progenitors. At functional level, increased invasiveness, at least in part mediated by SRPX2, and macrophage infiltration characterise this subset of glioblastoma.

Knowledge of the perturbation design is essential for accurate gene regulatory network inference.

The gene regulatory network (GRN) of a cell executes genetic programs in response to environmental and internal cues. Two distinct classes of methods are used to infer regulatory interactions from gene expression: those that only use observed changes in gene expression, and those that use both the observed changes and the perturbation design, i.e. the targets used to cause the changes in gene expression. Considering that the GRN by definition converts input cues to changes in gene expression, it may be conjectured that the latter methods would yield more accurate inferences but this has not previously been investigated. To address this question, we evaluated a number of popular GRN inference methods that either use the perturbation design or not. For the evaluation we used targeted perturbation knockdown gene expression datasets with varying noise levels generated by two different packages, GeneNetWeaver and GeneSpider. The accuracy was evaluated on each dataset using a variety of measures. The results show that on all datasets, methods using the perturbation design matrix consistently and significantly outperform methods not using it. This was also found to be the case on a smaller experimental dataset from E. coli. Targeted gene perturbations combined with inference methods that use the perturbation design are indispensable for accurate GRN inference.

Optimal Sparsity Selection Based on an Information Criterion for Accurate Gene Regulatory Network Inference.

Accurate inference of gene regulatory networks (GRNs) is important to unravel unknown regulatory mechanisms and processes, which can lead to the identification of treatment targets for genetic diseases. A variety of GRN inference methods have been proposed that, under suitable data conditions, perform well in benchmarks that consider the entire spectrum of false-positives and -negatives. However, it is very challenging to predict which single network sparsity gives the most accurate GRN. Lacking criteria for sparsity selection, a simplistic solution is to pick the GRN that has a certain number of links per gene, which is guessed to be reasonable. However, this does not guarantee finding the GRN that has the correct sparsity or is the most accurate one. In this study, we provide a general approach for identifying the most accurate and sparsity-wise relevant GRN within the entire space of possible GRNs. The algorithm, called SPA, applies a "GRN information criterion" (GRNIC) that is inspired by two commonly used model selection criteria, Akaike and Bayesian Information Criterion (AIC and BIC) but adapted to GRN inference. The results show that the approach can, in most cases, find the GRN whose sparsity is close to the true sparsity and close to as accurate as possible with the given GRN inference method and data. The datasets and source code can be found at https://bitbucket.org/sonnhammergrni/spa/.

Autocrine signaling can explain the emergence of Allee effects in cancer cell populations.

In many human cancers, the rate of cell growth depends crucially on the size of the tumor cell population. Low, zero, or negative growth at low population densities is known as the Allee effect; this effect has been studied extensively in ecology, but so far lacks a good explanation in the cancer setting. Here, we formulate and analyze an individual-based model of cancer, in which cell division rates are increased by the local concentration of an autocrine growth factor produced by the cancer cells themselves. We show, analytically and by simulation, that autocrine signaling suffices to cause both strong and weak Allee effects. Whether low cell densities lead to negative (strong effect) or reduced (weak effect) growth rate depends directly on the ratio of cell death to proliferation, and indirectly on cellular dispersal. Our model is consistent with experimental observations from three patient-derived brain tumor cell lines grown at different densities. We propose that further studying and quantifying population-wide feedback, impacting cell growth, will be central for advancing our understanding of cancer dynamics and treatment, potentially exploiting Allee effects for therapy.

Fast and accurate gene regulatory network inference by normalized least squares regression.

MOTIVATION: Inferring an accurate gene regulatory network (GRN) has long been a key goal in the field of systems biology. To do this, it is important to find a suitable balance between the maximum number of true positive and the minimum number of false-positive interactions. Another key feature is that the inference method can handle the large size of modern experimental data, meaning the method needs to be both fast and accurate. The Least Squares Cut-Off (LSCO) method can fulfill both these criteria, however as it is based on least squares it is vulnerable to known issues of amplifying extreme values, small or large. In GRN this manifests itself with genes that are erroneously hyper-connected to a large fraction of all genes due to extremely low value fold changes. RESULTS: We developed a GRN inference method called Least Squares Cut-Off with Normalization (LSCON) that tackles this problem. LSCON extends the LSCO algorithm by regularization to avoid hyper-connected genes and thereby reduce false positives. The regularization used is based on normalization, which removes effects of extreme values on the fit. We benchmarked LSCON and compared it to Genie3, LASSO, LSCO and Ridge regression, in terms of accuracy, speed and tendency to predict hyper-connected genes. The results show that LSCON achieves better or equal accuracy compared to LASSO, the best existing method, especially for data with extreme values. Thanks to the speed of least squares regression, LSCON does this an order of magnitude faster than LASSO. AVAILABILITY AND IMPLEMENTATION: Data: https://bitbucket.org/sonnhammergrni/lscon; Code: https://bitbucket.org/sonnhammergrni/genespider. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MYCMI-7: A Small MYC-Binding Compound that Inhibits MYC: MAX Interaction and Tumor Growth in a MYC-Dependent Manner.

Deregulated expression of MYC family oncogenes occurs frequently in human cancer and is often associated with aggressive disease and poor prognosis. While MYC is a highly warranted target, it has been considered "undruggable," and no specific anti-MYC drugs are available in the clinic. We recently identified molecules named MYCMIs that inhibit the interaction between MYC and its essential partner MAX. Here we show that one of these molecules, MYCMI-7, efficiently and selectively inhibits MYC:MAX and MYCN:MAX interactions in cells, binds directly to recombinant MYC, and reduces MYC-driven transcription. In addition, MYCMI-7 induces degradation of MYC and MYCN proteins. MYCMI-7 potently induces growth arrest/apoptosis in tumor cells in a MYC/MYCN-dependent manner and downregulates the MYC pathway on a global level as determined by RNA sequencing. Sensitivity to MYCMI-7 correlates with MYC expression in a panel of 60 tumor cell lines and MYCMI-7 shows high efficacy toward a collection of patient-derived primary glioblastoma and acute myeloid leukemia (AML) ex vivo cultures. Importantly, a variety of normal cells become G1 arrested without signs of apoptosis upon MYCMI-7 treatment. Finally, in mouse tumor models of MYC-driven AML, breast cancer, and MYCN-amplified neuroblastoma, treatment with MYCMI-7 downregulates MYC/MYCN, inhibits tumor growth, and prolongs survival through apoptosis with few side effects. In conclusion, MYCMI-7 is a potent and selective MYC inhibitor that is highly relevant for the development into clinically useful drugs for the treatment of MYC-driven cancer. Significance: Our findings demonstrate that the small-molecule MYCMI-7 binds MYC and inhibits interaction between MYC and MAX, thereby hampering MYC-driven tumor cell growth in culture and in vivo while sparing normal cells.

Real-time evaluation of glioblastoma growth in patient-specific zebrafish xenografts.

BACKGROUND: Patient-derived xenograft (PDX) models of glioblastoma (GBM) are a central tool for neuro-oncology research and drug development, enabling the detection of patient-specific differences in growth, and in vivo drug response. However, existing PDX models are not well suited for large-scale or automated studies. Thus, here, we investigate if a fast zebrafish-based PDX model, supported by longitudinal, AI-driven image analysis, can recapitulate key aspects of glioblastoma growth and enable case-comparative drug testing. METHODS: We engrafted 11 GFP-tagged patient-derived GBM IDH wild-type cell cultures (PDCs) into 1-day-old zebrafish embryos, and monitored fish with 96-well live microscopy and convolutional neural network analysis. Using light-sheet imaging of whole embryos, we analyzed further the invasive growth of tumor cells. RESULTS: Our pipeline enables automatic and robust longitudinal observation of tumor growth and survival of individual fish. The 11 PDCs expressed growth, invasion and survival heterogeneity, and tumor initiation correlated strongly with matched mouse PDX counterparts (Spearman R = 0.89, p < 0.001). Three PDCs showed a high degree of association between grafted tumor cells and host blood vessels, suggesting a perivascular invasion phenotype. In vivo evaluation of the drug marizomib, currently in clinical trials for GBM, showed an effect on fish survival corresponding to PDC in vitro and in vivo marizomib sensitivity. CONCLUSIONS: Zebrafish xenografts of GBM, monitored by AI methods in an automated process, present a scalable alternative to mouse xenograft models for the study of glioblastoma tumor initiation, growth, and invasion, applicable to patient-specific drug evaluation.

p53-Mediated Radiosensitization of 177Lu-DOTATATE in Neuroblastoma Tumor Spheroids.

p53 is involved in DNA damage response and is an exciting target for radiosensitization in cancer. Targeted radionuclide therapy against somatostatin receptors with 177Lu-DOTATATE is currently being explored as a treatment for neuroblastoma. The aim of this study was to investigate the novel p53-stabilizing peptide VIP116 in neuroblastoma, both as monotherapy and together with 177Lu-DOTATATE. Five neuroblastoma cell lines, including two patient-derived xenograft (PDX) lines, were characterized in monolayer cultures. Four out of five were positive for 177Lu-DOTATATE uptake. IC50 values after VIP116 treatments correlated with p53 status, ranging between 2.8-238.2 µM. IMR-32 and PDX lines LU-NB-1 and LU-NB-2 were then cultured as multicellular tumor spheroids and treated with 177Lu-DOTATATE and/or VIP116. Spheroid growth was inhibited in all spheroid models for all treatment modalities. The most pronounced effects were observed for combination treatments, mediating synergistic effects in the IMR-32 model. VIP116 and combination treatment increased p53 levels with subsequent induction of p21, Bax and cleaved caspase 3. Combination treatment resulted in a 14-fold and 1.6-fold induction of MDM2 in LU-NB-2 and IMR-32 spheroids, respectively. This, together with differential MYCN signaling, may explain the varying degree of synergy. In conclusion, VIP116 inhibited neuroblastoma cell growth, potentiated 177Lu-DOTATATE treatment and could, therefore, be a feasible treatment option for neuroblastoma.

Comparative epigenetic analysis of tumour initiating cells and syngeneic EPSC-derived neural stem cells in glioblastoma.

Epigenetic mechanisms which play an essential role in normal developmental processes, such as self-renewal and fate specification of neural stem cells (NSC) are also responsible for some of the changes in the glioblastoma (GBM) genome. Here we develop a strategy to compare the epigenetic and transcriptional make-up of primary GBM cells (GIC) with patient-matched expanded potential stem cell (EPSC)-derived NSC (iNSC). Using a comparative analysis of the transcriptome of syngeneic GIC/iNSC pairs, we identify a glycosaminoglycan (GAG)-mediated mechanism of recruitment of regulatory T cells (Tregs) in GBM. Integrated analysis of the transcriptome and DNA methylome of GBM cells identifies druggable target genes and patient-specific prediction of drug response in primary GIC cultures, which is validated in 3D and in vivo models. Taken together, we provide a proof of principle that this experimental pipeline has the potential to identify patient-specific disease mechanisms and druggable targets in GBM.

Modeling glioblastoma heterogeneity as a dynamic network of cell states.

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.

NF1 regulates mesenchymal glioblastoma plasticity and aggressiveness through the AP-1 transcription factor FOSL1.

The molecular basis underlying glioblastoma (GBM) heterogeneity and plasticity is not fully understood. Using transcriptomic data of human patient-derived brain tumor stem cell lines (BTSCs), classified based on GBM-intrinsic signatures, we identify the AP-1 transcription factor FOSL1 as a key regulator of the mesenchymal (MES) subtype. We provide a mechanistic basis to the role of the neurofibromatosis type 1 gene (NF1), a negative regulator of the RAS/MAPK pathway, in GBM mesenchymal transformation through the modulation of FOSL1 expression. Depletion of FOSL1 in NF1-mutant human BTSCs and Kras-mutant mouse neural stem cells results in loss of the mesenchymal gene signature and reduction in stem cell properties and in vivo tumorigenic potential. Our data demonstrate that FOSL1 controls GBM plasticity and aggressiveness in response to NF1 alterations.

Single-cell RNAseq and longitudinal proteomic analysis of a novel semi-spontaneous urothelial cancer model reveals tumor cell heterogeneity and pretumoral urine protein alterations.

Bladder cancer, one of the most prevalent malignancies worldwide, remains hard to classify due to a staggering molecular complexity. Despite a plethora of diagnostic tools and therapies, it is hard to outline the key steps leading up to the transition from high-risk non-muscle-invasive bladder cancer (NMIBC) to muscle-invasive bladder cancer (MIBC). Carcinogen-induced murine models can recapitulate urothelial carcinogenesis and natural anti-tumor immunity. Herein, we have developed and profiled a novel model of progressive NMIBC based on 10 weeks of OH-BBN exposure in hepatocyte growth factor/cyclin dependent kinase 4 (R24C) (Hgf-Cdk4R24C) mice. The profiling of the model was performed by histology grading, single cell transcriptomic and proteomic analysis, while the derivation of a tumorigenic cell line was validated and used to assess in vivo anti-tumor effects in response to immunotherapy. Established NMIBC was present in females at 10 weeks post OH-BBN exposure while neoplasia was not as advanced in male mice, however all mice progressed to MIBC. Single cell RNA sequencing analysis revealed an intratumoral heterogeneity also described in the human disease trajectory. Moreover, although immune activation biomarkers were elevated in urine during carcinogen exposure, anti-programmed cell death protein 1 (anti-PD1) monotherapy did not prevent tumor progression. Furthermore, anti-PD1 immunotherapy did not control the growth of subcutaneous tumors formed by the newly derived urothelial cancer cell line. However, treatment with CpG-oligodeoxynucleotides (ODN) significantly decreased tumor volume, but only in females. In conclusion, the molecular map of this novel preclinical model of bladder cancer provides an opportunity to further investigate pharmacological therapies ahead with regards to both targeted drugs and immunotherapies to improve the strategies of how we should tackle the heterogeneous tumor microenvironment in urothelial bladder cancer to improve responses rates in the clinic.

Inferring the experimental design for accurate gene regulatory network inference.

MOTIVATION: Accurate inference of gene regulatory interactions is of importance for understanding the mechanisms of underlying biological processes. For gene expression data gathered from targeted perturbations, gene regulatory network (GRN) inference methods that use the perturbation design are the top performing methods. However, the connection between the perturbation design and gene expression can be obfuscated due to problems, such as experimental noise or off-target effects, limiting the methods' ability to reconstruct the true GRN. RESULTS: In this study, we propose an algorithm, IDEMAX, to infer the effective perturbation design from gene expression data in order to eliminate the potential risk of fitting a disconnected perturbation design to gene expression. We applied IDEMAX to synthetic data from two different data generation tools, GeneNetWeaver and GeneSPIDER, and assessed its effect on the experiment design matrix as well as the accuracy of the GRN inference, followed by application to a real dataset. The results show that our approach consistently improves the accuracy of GRN inference compared to using the intended perturbation design when much of the signal is hidden by noise, which is often the case for real data. AVAILABILITY AND IMPLEMENTATION: https://bitbucket.org/sonnhammergrni/idemax. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Uncovering cancer gene regulation by accurate regulatory network inference from uninformative data.

The interactions among the components of a living cell that constitute the gene regulatory network (GRN) can be inferred from perturbation-based gene expression data. Such networks are useful for providing mechanistic insights of a biological system. In order to explore the feasibility and quality of GRN inference at a large scale, we used the L1000 data where ~1000 genes have been perturbed and their expression levels have been quantified in 9 cancer cell lines. We found that these datasets have a very low signal-to-noise ratio (SNR) level causing them to be too uninformative to infer accurate GRNs. We developed a gene reduction pipeline in which we eliminate uninformative genes from the system using a selection criterion based on SNR, until reaching an informative subset. The results show that our pipeline can identify an informative subset in an overall uninformative dataset, allowing inference of accurate subset GRNs. The accurate GRNs were functionally characterized and potential novel cancer-related regulatory interactions were identified.

Monotherapy efficacy of blood-brain barrier permeable small molecule reactivators of protein phosphatase 2A in glioblastoma.

Glioblastoma is a fatal disease in which most targeted therapies have clinically failed. However, pharmacological reactivation of tumour suppressors has not been thoroughly studied as yet as a glioblastoma therapeutic strategy. Tumour suppressor protein phosphatase 2A is inhibited by non-genetic mechanisms in glioblastoma, and thus, it would be potentially amendable for therapeutic reactivation. Here, we demonstrate that small molecule activators of protein phosphatase 2A, NZ-8-061 and DBK-1154, effectively cross the in vitro model of blood-brain barrier, and in vivo partition to mouse brain tissue after oral dosing. In vitro, small molecule activators of protein phosphatase 2A exhibit robust cell-killing activity against five established glioblastoma cell lines, and nine patient-derived primary glioma cell lines. Collectively, these cell lines have heterogeneous genetic background, kinase inhibitor resistance profile and stemness properties; and they represent different clinical glioblastoma subtypes. Moreover, small molecule activators of protein phosphatase 2A were found to be superior to a range of kinase inhibitors in their capacity to kill patient-derived primary glioma cells. Oral dosing of either of the small molecule activators of protein phosphatase 2A significantly reduced growth of infiltrative intracranial glioblastoma tumours. DBK-1154, with both higher degree of brain/blood distribution, and more potent in vitro activity against all tested glioblastoma cell lines, also significantly increased survival of mice bearing orthotopic glioblastoma xenografts. In summary, this report presents a proof-of-principle data for blood-brain barrier-permeable tumour suppressor reactivation therapy for glioblastoma cells of heterogenous molecular background. These results also provide the first indications that protein phosphatase 2A reactivation might be able to challenge the current paradigm in glioblastoma therapies which has been strongly focused on targeting specific genetically altered cancer drivers with highly specific inhibitors. Based on demonstrated role for protein phosphatase 2A inhibition in glioblastoma cell drug resistance, small molecule activators of protein phosphatase 2A may prove to be beneficial in future glioblastoma combination therapies.

HCN Channel Activity Balances Quiescence and Proliferation in Neural Stem Cells and Is a Selective Target for Neuroprotection During Cancer Treatment.

Children suffering from neurologic cancers undergoing chemotherapy and radiotherapy are at high risk of reduced neurocognitive abilities likely via damage to proliferating neural stem cells (NSC). Therefore, strategies to protect NSCs are needed. We argue that induced cell-cycle arrest/quiescence in NSCs during cancer treatment can represent such a strategy. Here, we show that hyperpolarization-activated cyclic nucleotide-gated (HCN) ion channels are dynamically expressed over the cell cycle in NSCs, depolarize the membrane potential, underlie spontaneous calcium oscillations and are required to maintain NSCs in the actively proliferating pool. Hyperpolarizing pharmacologic inhibition of HCN channels during exposure to ionizing radiation protects NSCs cells in neurogenic brain regions of young mice. In contrast, brain tumor-initiating cells, which also express HCN channels, remain proliferative during HCN inhibition. IMPLICATIONS: Our finding that NSCs can be selectively rescued while cancer cells remain sensitive to the treatment, provide a foundation for reduction of cognitive impairment in children with neurologic cancers.