Search

Abstract Nervous system cancers exhibit diverse transcriptional cell states influenced by normal development, injury response, and growth. However, the understanding of these states’ regulation and pharmacological relevance remains limited. Here we present “single-cell regulatory-driven clustering” (scregclust), a method that reconstructs cellular regulatory programs from extensive collections of single-cell RNA sequencing (scRNA-seq) data from both tumors and developing tissues. The algorithm efficiently divides target genes into modules, predicting key transcription factors and kinases with minimal computational time. Applying this method to adult and childhood brain cancers, we identify critical regulators and suggest interventions that could improve temozolomide treatment in glioblastoma. Additionally, our integrative analysis reveals a meta-module regulated by SPI1 and IRF8 linked to an immune-mediated mesenchymal-like state. Finally, scregclust’s flexibility is demonstrated across 15 tumor types, uncovering both pan-cancer and specific regulators. The algorithm is provided as an easy-to-use R package that facilitates the exploration of regulatory programs underlying cell plasticity.

Summary: 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.

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

Abstract 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.

Abstract 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.

The identification of patient-specific disease mechanisms and druggable targets is crucial for precision medicine in glioblastoma. Here, the authors show that comparing patients-matched glioma-initiating cells with neural stem cells enables the discovery of patient-specific mechanisms of disease and the identification of effective drugs

Abstract 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.

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.

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/.

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.

We lack effective treatment for half of children with high-risk neuroblastoma. Here, the authors introduce an algorithm that can predict the effect of interventions on gene expression signatures associated with high disease processes and risk, and identify and validate promising drug targets.

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.

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.

Summary: Glioblastoma (GBM) is a malignant brain tumor with few therapeutic options. The disease presents with a complex spectrum of genomic aberrations, but the pharmacological consequences of these aberrations are partly unknown. Here, we report an integrated pharmacogenomic analysis of 100 patient-derived GBM cell cultures from the human glioma cell culture (HGCC) cohort. Exploring 1,544 drugs, we find that GBM has two main pharmacological subgroups, marked by differential response to proteasome inhibitors and mutually exclusive aberrations in TP53 and CDKN2A/B. We confirm this trend in cell and in xenotransplantation models, and identify both Bcl-2 family inhibitors and p53 activators as potentiators of proteasome inhibitors in GBM cells. We can further predict the responses of individual cell cultures to several existing drug classes, presenting opportunities for drug repurposing and design of stratified trials. Our functionally profiled biobank provides a valuable resource for the discovery of new treatments for GBM.

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.

Glioblastoma (GBM) is the most aggressive malignant primary brain tumor in adults, with a median survival of 14.6 months. Recent efforts have focused on identifying clinically relevant subgroups to improve our understanding of pathogenetic mechanisms and patient stratification. Concurrently, the role of immune cells in the tumor microenvironment has received increasing attention, especially T cells and tumor-associated macrophages (TAM). The latter are a mixed population of activated brain-resident microglia and infiltrating monocytes/monocyte-derived macrophages, both of which express ionized calcium-binding adapter molecule 1 (IBA1). This study investigated differences in immune cell subpopulations among distinct transcriptional subtypes of GBM. Human GBM samples were molecularly characterized and assigned to Proneural, Mesenchymal or Classical subtypes as defined by NanoString nCounter Technology. Subsequently, we performed and analyzed automated immunohistochemical stainings for TAM as well as specific T cell populations. The Mesenchymal subtype of GBM showed the highest presence of TAM, CD8+, CD3+ and FOXP3+ T cells, as compared to Proneural and Classical subtypes. High expression levels of the TAM-related gene AIF1, which encodes the TAM-specific protein IBA1, correlated with a worse prognosis in Proneural GBM, but conferred a survival benefit in Mesenchymal tumors. We used our data to construct a mathematical model that could reliably identify Mesenchymal GBM with high sensitivity using a combination of the aforementioned cell-specific IHC markers. In conclusion, we demonstrated that molecularly distinct GBM subtypes are characterized by profound differences in the composition of their immune microenvironment, which could potentially help to identify tumors amenable to immunotherapy.

The identity of the glioblastoma (GBM) cell of origin and its contributions to disease progression and treatment response remain largely unknown. We have analyzed how the phenotypic state of the initially transformed cell affects mouse GBM development and essential GBM cell (GC) properties. We find that GBM induced in neural stem-cell-like glial fibrillary acidic protein (GFAP)-expressing cells in the subventricular zone of adult mice shows accelerated tumor development and produces more malignant GCs (mGC1GFAP) that are less resistant to cancer drugs, compared with those originating from more differentiated nestin- (mGC2NES) or 2,′3′-cyclic nucleotide 3′-phosphodiesterase (mGC3CNP)-expressing cells. Transcriptome analysis of mouse GCs identified a 196 mouse cell origin (MCO) gene signature that was used to partition 61 patient-derived GC lines. Human GC lines that clustered with the mGC1GFAP cells were also significantly more self-renewing, tumorigenic, and sensitive to cancer drugs compared with those that clustered with mouse GCs of more differentiated origin.

Glioblastomas are characterized by transcriptionally distinct subtypes, but despite possible clinical relevance, their regulation remains poorly understood. The commonly used molecular classification systems for GBM all identify a subtype with high expression of mesenchymal marker transcripts, strongly associated with invasive growth. We used a comprehensive data-driven network modeling technique (augmented sparse inverse covariance selection, aSICS) to define separate genomic, epigenetic, and transcriptional regulators of glioblastoma subtypes. Our model identified Annexin A2 (ANXA2) as a novel methylation-controlled positive regulator of the mesenchymal subtype. Subsequent evaluation in two independent cohorts established ANXA2 expression as a prognostic factor that is dependent on ANXA2 promoter methylation. ANXA2 knockdown in primary glioblastoma stem cell-like cultures suppressed known mesenchymal master regulators, and abrogated cell proliferation and invasion. Our results place ANXA2 at the apex of a regulatory cascade that determines glioblastoma mesenchymal transformation and validate aSICS as a general methodology to uncover regulators of cancer subtypes.

Previous genome-wide association studies have identified loci associated with the risk of multiple myeloma. Here, the authors present a meta-analysis of six genome wide association studies of the disease and identify eight new loci; functional studies identify genes as candidates for the basis of these associations.

The transcriptional classification of glioblastoma has proven to be a complex issue. In the absence of strong correlations between underlying genomic lesions and transcriptional subtype, there is a need to systematically understand the origins of the glioblastoma subtypes. A recent integrated analysis of data from both mouse models and patient-derived cells supports that the glioblastoma's cell of origin is important in shaping transcriptional diversity and tumor cell malignancy.

Significant advances have been made in methods to analyze genomes and transcriptomes of single cells, but to fully define cell states, proteins must also be accessed as central actors defining a cell’s phenotype. Methods currently used to analyze endogenous protein expression in single cells are limited in specificity, throughput, or multiplex capability. Here, we present an approach to simultaneously and specifically interrogate large sets of protein and RNA targets in lysates from individual cells, enabling investigations of cell functions and responses. We applied our method to investigate the effects of BMP4, an experimental therapeutic agent, on early-passage glioblastoma cell cultures. We uncovered significant heterogeneity in responses to treatment at levels of RNA and protein, with a subset of cells reacting in a distinct manner to BMP4. Moreover, we found overall poor correlation between protein and RNA at the level of single cells, with proteins more accurately defining responses to treatment.

Glioblastoma (GBM) is the most frequent and malignant form of primary brain tumor. GBM is essentially incurable and its resistance to therapy is attributed to a subpopulation of cells called glioma stem cells (GSCs). To meet the present shortage of relevant GBM cell (GC) lines we developed a library of annotated and validated cell lines derived from surgical samples of GBM patients, maintained under conditions to preserve GSC characteristics. This collection, which we call the Human Glioblastoma Cell Culture (HGCC) resource, consists of a biobank of 48 GC lines and an associated database containing high-resolution molecular data. We demonstrate that the HGCC lines are tumorigenic, harbor genomic lesions characteristic of GBMs, and represent all four transcriptional subtypes. The HGCC panel provides an open resource for in vitro and in vivo modeling of a large part of GBM diversity useful to both basic and translational GBM research.

Abstract DNA copy number aberrations (CNAs) are a hallmark of cancer genomes. However, little is known about how such changes affect global gene expression. We develop a modeling framework, EPoC (Endogenous Perturbation analysis of Cancer), to (1) detect disease‐driving CNAs and their effect on target mRNA expression, and to (2) stratify cancer patients into long‐ and short‐term survivors. Our method constructs causal network models of gene expression by combining genome‐wide DNA‐ and RNA‐level data. Prognostic scores are obtained from a singular value decomposition of the networks. By applying EPoC to glioblastoma data from The Cancer Genome Atlas consortium, we demonstrate that the resulting network models contain known disease‐relevant hub genes, reveal interesting candidate hubs, and uncover predictors of patient survival. Targeted validations in four glioblastoma cell lines support selected predictions, and implicate the p53‐interacting protein Necdin in suppressing glioblastoma cell growth. We conclude that large‐scale network modeling of the effects of CNAs on gene expression may provide insights into the biology of human cancer. Free software in MATLAB and R is provided.

Tumor-initiating cells are a subpopulation in aggressive cancers that exhibit traits shared with stem cells, including the ability to self-renew and differentiate, commonly referred to as stemness. In addition, such cells are resistant to chemo- and radiation therapy posing a therapeutic challenge. To uncover stemness-associated functions in glioma-initiating cells (GICs), transcriptome profiles were compared to neural stem cells (NSCs) and gene ontology analysis identified an enrichment of Ca2+ signaling genes in NSCs and the more stem-like (NSC-proximal) GICs. Functional analysis in a set of different GIC lines regarding sensitivity to disturbed homeostasis using A23187 and Thapsigargin, revealed that NSC-proximal GICs were more sensitive, corroborating the transcriptome data. Furthermore, Ca2+ drug sensitivity was reduced in GICs after differentiation, with most potent effect in the NSC-proximal GIC, supporting a stemness-associated Ca2+ sensitivity. NSCs and the NSC-proximal GIC line expressed a larger number of ion channels permeable to potassium, sodium and Ca2+. Conversely, a higher number of and higher expression levels of Ca2+ binding genes that may buffer Ca2+, were expressed in NSC-distal GICs. In particular, expression of the AMPA glutamate receptor subunit GRIA1, was found to associate with Ca2+ sensitive NSC-proximal GICs, and decreased as GICs differentiated along with reduced Ca2+ drug sensitivity. The correlation between high expression of Ca2+ channels (such as GRIA1) and sensitivity to Ca2+ drugs was confirmed in an additional nine novel GIC lines. Calcium drug sensitivity also correlated with expression of the NSC markers nestin (NES) and FABP7 (BLBP, brain lipid-binding protein) in this extended analysis. In summary, NSC-associated NES+/FABP7+/GRIA1+ GICs were selectively sensitive to disturbances in Ca2+ homeostasis, providing a potential target mechanism for eradication of an immature population of malignant cells.

Functionally interacting perturbations, such as synergistic drugs pairs or synthetic lethal gene pairs, are of key interest in both pharmacology and functional genomics. However, to find such pairs by traditional screening methods is both time consuming and costly. We present a novel computational-experimental framework for efficient identification of synergistic target pairs, applicable for screening of systems with sizes on the order of current drug, small RNA or SGA (Synthetic Genetic Array) libraries (>1000 targets). This framework exploits the fact that the response of a drug pair in a given system, or a pair of genes' propensity to interact functionally, can be partly predicted by computational means from (i) a small set of experimentally determined target pairs, and (ii) pre-existing data (e.g. gene ontology, PPI) on the similarities between targets. Predictions are obtained by a novel matrix algebraic technique, based on cyclical projections onto convex sets. We demonstrate the efficiency of the proposed method using drug-drug interaction data from seven cancer cell lines and gene-gene interaction data from yeast SGA screens. Our protocol increases the rate of synergism discovery significantly over traditional screening, by up to 7-fold. Our method is easy to implement and could be applied to accelerate pair screening for both animal and microbial systems.

The brain tumour glioblastoma is characterised by diffuse and infiltrative growth into surrounding brain tissue. At the macroscopic level, the progression speed of a glioblastoma tumour is determined by two key factors: the cell proliferation rate and the cell migration speed. At the microscopic level, however, proliferation and migration appear to be mutually exclusive phenotypes, as indicated by recent in vivo imaging data. Here, we develop a mathematical model to analyse how the phenotypic switching between proliferative and migratory states of individual cells affects the macroscopic growth of the tumour. For this, we propose an individual-based stochastic model in which glioblastoma cells are either in a proliferative state, where they are stationary and divide, or in motile state in which they are subject to random motion. From the model we derive a continuum approximation in the form of two coupled reaction-diffusion equations, which exhibit travelling wave solutions whose speed of invasion depends on the model parameters. We propose a simple analytical method to predict progression rate from the cell-specific parameters and demonstrate that optimal glioblastoma growth depends on a non-trivial trade-off between the phenotypic switching rates. By linking cellular properties to an in vivo outcome, the model should be applicable to designing relevant cell screens for glioblastoma and cytometry-based patient prognostics.