Computational modeling demonstrates that glioblastoma cells can survive spatial environmental challenges through exploratory adaptation

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell-state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing nongenetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupt acquired resistance in GBM.

Prioritization of genes for translation: a computational approach

INTRODUCTION

Gene identification for genetic diseases is critical for the development of new diagnostic approaches and personalized treatment options. Prioritization of gene translation is an important consideration in the molecular biology field, allowing researchers to focus on the most promising candidates for further investigation.

AREAS COVERED

In this paper, we discussed different approaches to prioritize genes for translation, including the use of computational tools and machine learning algorithms, as well as experimental techniques such as knockdown and overexpression studies. We also explored the potential biases and limitations of these approaches and proposed strategies to improve the accuracy and reliability of gene prioritization methods. Although numerous computational methods have been developed for this purpose, there is a need for computational methods that incorporate tissue-specific information to enable more accurate prioritization of candidate genes. Such methods should provide tissue-specific predictions, insights into underlying disease mechanisms, and more accurate prioritization of genes.

EXPERT OPINION

Using advanced computational tools and machine learning algorithms to prioritize genes, we can identify potential targets for therapeutic intervention of complex diseases. This represents an up-and-coming method for drug development and personalized medicine.

The impact of MEIS1 TALE homeodomain transcription factor knockdown on glioma stem cell growth

Myeloid ecotropic virus insertion site 1 (MEIS1) is a HOX co-factor necessary for organ development and normal hematopoiesis. Recently, MEIS1 has been linked to the development and progression of various cancers. However, its role in gliomagenesis particularly on glioma stem cells (GSCs) remains unclear. Here, we demonstrate that MEIS1 is highly upregulated in GSCs compared to normal, and glioma cells and to its differentiated counterparts. Inhibition of MEIS1 expression by shRNA significantly reduced GSC growth in both in vitro and in vivo experiments. On the other hand, integrated transcriptomics analyses of glioma datasets revealed that MEIS1 expression is correlated to cell cycle-related genes. Clinical data analysis revealed that MEIS1 expression is elevated in high-grade gliomas, and patients with high MEIS1 levels have poorer overall survival outcomes. The findings suggest that MEIS1 is a prognostic biomarker for glioma patients and a possible target for developing novel therapeutic strategies against GBM.

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing non-genetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupts acquired resistance in GBM.

Establishment and immune phenotyping of patient-derived glioblastoma models in humanized mice

Glioblastoma (GBM) is the most aggressive and common type of malignant brain tumor diagnosed in adults. Preclinical immunocompetent mouse tumor models generated using mouse tumor cells play a pivotal role in testing the therapeutic efficacy of emerging immune-based therapies for GBMs. However, the clinical translatability of such studies is limited as mouse tumor lines do not fully recapitulate GBMs seen in inpatient settings. In this study, we generated three distinct, imageable human-GBM (hGBM) models in humanized mice using patient-derived GBM cells that cover phenotypic and genetic GBM heterogeneity in primary (invasive and nodular) and recurrent tumors. We developed a pipeline to first enrich the tumor-initiating stem-like cells and then successfully established robust patient-derived GBM tumor engraftment and growth in bone marrow-liver-thymus (BLT) humanized mice. Multiplex immunofluorescence of GBM tumor sections revealed distinct phenotypic features of the patient GBM tumors, with myeloid cells dominating the immune landscape. Utilizing flow cytometry and correlative immunofluorescence, we profiled the immune microenvironment within the established human GBM tumors in the BLT mouse models and showed tumor infiltration of variable human immune cells, creating a unique immune landscape compared with lymphoid organs. These findings contribute substantially to our understanding of GBM biology within the context of the human immune system in humanized mice and lay the groundwork for further translational studies aimed at advancing therapeutic strategies for GBM.

The tumour ecology of quiescence: Niches across scales of complexity

Quiescence is a state of cell cycle arrest, allowing cancer cells to evade anti-proliferative cancer therapies. Quiescent cancer stem cells are thought to be responsible for treatment resistance in glioblastoma, an aggressive brain cancer with poor patient outcomes. However, the regulation of quiescence in glioblastoma cells involves a myriad of intrinsic and extrinsic mechanisms that are not fully understood. In this review, we synthesise the literature on quiescence regulatory mechanisms in the context of glioblastoma and propose an ecological perspective to stemness-like phenotypes anchored to the contemporary concepts of niche theory. From this perspective, the cell cycle regulation is multiscale and multidimensional, where the niche dimensions extend to extrinsic variables in the tumour microenvironment that shape cell fate. Within this conceptual framework and powered by ecological niche modelling, the discovery of microenvironmental variables related to hypoxia and mechanosignalling that modulate proliferative plasticity and intratumor immune activity may open new avenues for therapeutic targeting of emerging biological vulnerabilities in glioblastoma.

Cancer cell heterogeneity & plasticity in glioblastoma and brain tumors

Brain tumors remain one of the most difficult tumors to treat and, depending on the diagnosis, have a poor prognosis. Of brain tumors, glioblastoma (GBM) is the most common malignant glioma and has a dismal prognosis, with only about 5% of patients alive five years after diagnosis. While advances in targeted therapies and immunotherapies are rapidly improving outcomes in a variety of other cancers, the standard of care for GBM has largely remained unaltered since 2005. There are many well-studied challenges that are either unique to brain tumors (i.e., blood-brain barrier and immunosuppressive environment) or amplified within GBM (i.e., tumor heterogeneity at the cellular and molecular levels, plasticity, and cancer stem cells) that make this disease particularly difficult to treat. While we touch on all these concepts, the focus of this review is to discuss the immense inter- and intra-tumoral heterogeneity and advances in our understanding of tumor cell plasticity and epigenetics in GBM. With each improvement in technology, our understanding of the complexity of tumoral heterogeneity and plasticity improves and we gain more clarity on the causes underlying previous therapeutic failures. However, these advances are unlocking new therapeutic opportunities that scientists and physicians are currently exploiting and have the potential for new breakthroughs.

Cancer progression as a learning process

Drug resistance and metastasis—the major complications in cancer—both entail adaptation of cancer cells to stress, whether a drug or a lethal new environment. Intriguingly, these adaptive processes share similar features that cannot be explained by a pure Darwinian scheme, including dormancy, increased heterogeneity, and stress-induced plasticity. Here, we propose that learning theory offers a framework to explain these features and may shed light on these two intricate processes. In this framework, learning is performed at the single-cell level, by stress-driven exploratory trial-and-error. Such a process is not contingent on pre-existing pathways but on a random search for a state that diminishes the stress. We review underlying mechanisms that may support this search, and show by using a learning model that such exploratory learning is feasible in a high-dimensional system as the cell. At the population level, we view the tissue as a network of exploring agents that communicate, restraining cancer formation in health. In this view, disease results from the breakdown of homeostasis between cellular exploratory drive and tissue homeostasis.

Topographic mapping of the glioblastoma proteome reveals a triple-axis model of intra-tumoral heterogeneity

Glioblastoma is an aggressive form of brain cancer with well-established patterns of intra-tumoral heterogeneity implicated in treatment resistance and progression. While regional and single cell transcriptomic variations of glioblastoma have been recently resolved, downstream phenotype-level proteomic programs have yet to be assigned across glioblastoma's hallmark histomorphologic niches. Here, we leverage mass spectrometry to spatially align abundance levels of 4,794 proteins to distinct histologic patterns across 20 patients and propose diverse molecular programs operational within these regional tumor compartments. Using machine learning, we overlay concordant transcriptional information, and define two distinct proteogenomic programs, MYC- and KRAS-axis hereon, that cooperate with hypoxia to produce a tri-dimensional model of intra-tumoral heterogeneity. Moreover, we highlight differential drug sensitivities and relative chemoresistance in glioblastoma cell lines with enhanced KRAS programs. Importantly, these pharmacological differences are less pronounced in transcriptional glioblastoma subgroups suggesting that this model may provide insights for targeting heterogeneity and overcoming therapy resistance.

The novel BET inhibitor UM-002 reduces glioblastoma cell proliferation and invasion

Bromodomain and extraterminal domain (BET) proteins have emerged as therapeutic targets in multiple cancers, including the most common primary adult brain tumor glioblastoma (GBM). Although several BET inhibitors have entered clinical trials, few are brain penetrant. We have generated UM-002, a novel brain penetrant BET inhibitor that reduces GBM cell proliferation in vitro and in a human cerebral brain organoid model. Since UM-002 is more potent than other BET inhibitors, it could potentially be developed for GBM treatment. Furthermore, UM-002 treatment reduces the expression of cell-cycle related genes in vivo and reduces the expression of invasion related genes within the non-proliferative cells present in tumors as measured by single cell RNA-sequencing. These studies suggest that BET inhibition alters the transcriptional landscape of GBM tumors, which has implications for designing combination therapies. Importantly, they also provide an integrated dataset that combines in vitro and ex vivo studies with in vivo single-cell RNA-sequencing to characterize a novel BET inhibitor in GBM.

Spatiotemporal dynamics of a glioma immune interaction model

We report a mathematical model which depicts the spatiotemporal dynamics of glioma cells, macrophages, cytotoxic-T-lymphocytes, immuno-suppressive cytokine TGF-β and immuno-stimulatory cytokine IFN-γ through a system of five coupled reaction-diffusion equations. We performed local stability analysis of the biologically based mathematical model for the growth of glioma cell population and their environment. The presented stability analysis of the model system demonstrates that the temporally stable positive interior steady state remains stable under the small inhomogeneous spatiotemporal perturbations. The irregular spatiotemporal dynamics of gliomas, macrophages and cytotoxic T-lymphocytes are discussed extensively and some numerical simulations are presented. Performed some numerical simulations in both one and two dimensional spaces. The occurrence of heterogeneous pattern formation of the system has both biological and mathematical implications and the concepts of glioma cell progression and invasion are considered. Simulation of the model shows that by increasing the value of time, the glioma cell population, macrophages and cytotoxic-T-lymphocytes spread throughout the domain.

Association of Circadian Clock Gene Expression with Glioma Tumor Microenvironment and Patient Survival

Simple Summary

Gliomas are the most common type of malignant primary brain tumors and are classified according to the cell of origin and genetic features, which can help predict the prognosis and treatment sensitivity. Improving the prognosis remains a challenge; however, chronobiology is a promising field for future works, as circadian clock genes are linked to the tumor biology and outcomes in multiple cancers, including glioma. Here, we examined the relationship of circadian clock genes, IDH mutational status, and prognosis in glioma patients by using unsupervised clustering of the expression of 13 clock genes. We further explored the expression of the clock genes across the tumor regions and cell subpopulations, highlighting the importance of the tumor microenvironment in researching circadian rhythms in cancer. Our research is important for understanding how best to target circadian rhythms to improve patient outcomes in neuro-oncology.

Abstract

Circadian clock genes have been linked to clinical outcomes in cancer, including gliomas. However, these studies have not accounted for established markers that predict the prognosis, including mutations in Isocitrate Dehydrogenase (IDH), which characterize the majority of lower-grade gliomas and secondary high-grade gliomas. To demonstrate the connection between circadian clock genes and glioma outcomes while accounting for the IDH mutational status, we analyzed multiple publicly available gene expression datasets. The unsupervised clustering of 13 clock gene transcriptomic signatures from The Cancer Genome Atlas showed distinct molecular subtypes representing different disease states and showed the differential prognosis of these groups by a Kaplan–Meier analysis. Further analyses of these groups showed that a low period (PER) gene expression was associated with the negative prognosis and enrichment of the immune signaling pathways. These findings prompted the exploration of the relationship between the microenvironment and clock genes in additional datasets. Circadian clock gene expression was found to be differentially expressed across the anatomical tumor location and cell type. Thus, the circadian clock expression is a potential predictive biomarker in glioma, and further mechanistic studies to elucidate the connections between the circadian clock and microenvironment are warranted.

Tumor cell plasticity, heterogeneity, and resistance in crucial microenvironmental niches in glioma

Both the perivascular niche (PVN) and the integration into multicellular networks by tumor microtubes (TMs) have been associated with progression and resistance to therapies in glioblastoma, but their specific contribution remained unknown. By long-term tracking of tumor cell fate and dynamics in the live mouse brain, differential therapeutic responses in both niches are determined. Both the PVN, a preferential location of long-term quiescent glioma cells, and network integration facilitate resistance against cytotoxic effects of radiotherapy and chemotherapy—independently of each other, but with additive effects. Perivascular glioblastoma cells are particularly able to actively repair damage to tumor regions. Population of the PVN and resistance in it depend on proficient NOTCH1 expression. In turn, NOTCH1 downregulation induces resistant multicellular networks by TM extension. Our findings identify NOTCH1 as a central switch between the PVN and network niche in glioma, and demonstrate robust cross-compensation when only one niche is targeted.

Whether the perivascular niche (PVN) and the integration into multicellular networks by tumor microtubes (TMs) have a different role in glioblastoma progression and resistance to therapies is currently unclear. Here, the authors, by long-term tracking of individual glioma, demonstrate that both niches can partially compensate for each other and that glioma cells localized in both niches are resistant to radio- and chemotherapy.

A mesoscopic simulator to uncover heterogeneity and evolutionary dynamics in tumors

Increasingly complex in silico modeling approaches offer a way to simultaneously access cancerous processes at different spatio-temporal scales. High-level models, such as those based on partial differential equations, are computationally affordable and allow large tumor sizes and long temporal windows to be studied, but miss the discrete nature of many key underlying cellular processes. Individual-based approaches provide a much more detailed description of tumors, but have difficulties when trying to handle full-sized real cancers. Thus, there exists a trade-off between the integration of macroscopic and microscopic information, now widely available, and the ability to attain clinical tumor sizes. In this paper we put forward a stochastic mesoscopic simulation framework that incorporates key cellular processes during tumor progression while keeping computational costs to a minimum. Our framework captures a physical scale that allows both the incorporation of microscopic information, tracking the spatio-temporal emergence of tumor heterogeneity and the underlying evolutionary dynamics, and the reconstruction of clinically sized tumors from high-resolution medical imaging data, with the additional benefit of low computational cost. We illustrate the functionality of our modeling approach for the case of glioblastoma, a paradigm of tumor heterogeneity that remains extremely challenging in the clinical setting.

Benchmarking algorithms for pathway activity transformation of single-cell RNA-seq data

Biological pathway analysis provides new insights for cell clustering and functional annotation from single-cell RNA sequencing (scRNA-seq) data. Many pathway analysis algorithms have been developed to transform gene-level scRNA-seq data into functional gene sets representing pathways or biological processes. Here, we collected seven widely-used pathway activity transformation algorithms and 32 available datasets based on 16 scRNA-seq techniques. We proposed a comprehensive framework to evaluate their accuracy, stability and scalability. The assessment of scRNA-seq preprocessing showed that cell filtering had the less impact on scRNA-seq pathway analysis, while data normalization of sctransform and scran had a consistent well impact across all tools. We found that Pagoda2 yielded the best overall performance with the highest accuracy, scalability, and stability. Meanwhile, the tool PLAGE exhibited the highest stability, as well as moderate accuracy and scalability.

Tumor Cell Invasion in Glioblastoma

Glioblastoma (GBM) is a particularly devastating tumor with a median survival of about 16 months. Recent research has revealed novel insights into the outstanding heterogeneity of this type of brain cancer. However, all GBM subtypes share the hallmark feature of aggressive invasion into the surrounding tissue. Invasive glioblastoma cells escape surgery and focal therapies and thus represent a major obstacle for curative therapy. This review aims to provide a comprehensive understanding of glioma invasion mechanisms with respect to tumor-cell-intrinsic properties as well as cues provided by the microenvironment. We discuss genetic programs that may influence the dissemination and plasticity of GBM cells as well as their different invasion patterns. We also review how tumor cells shape their microenvironment and how, vice versa, components of the extracellular matrix and factors from non-neoplastic cells influence tumor cell motility. We further discuss different research platforms for modeling invasion. Finally, we highlight the importance of accounting for the complex interplay between tumor cell invasion and treatment resistance in glioblastoma when considering new therapeutic approaches.