WDR5 represents a therapeutically exploitable target for cancer stem cells in glioblastoma.

Glioblastomas (GBMs) are heterogeneous, treatment-resistant tumors driven by populations of cancer stem cells (CSCs). However, few molecular mechanisms critical for CSC population maintenance have been exploited for therapeutic development. We developed a spatially resolved loss-of-function screen in GBM patient-derived organoids to identify essential epigenetic regulators in the SOX2-enriched, therapy-resistant niche and identified WDR5 as indispensable for this population. WDR5 is a component of the WRAD complex, which promotes SET1 family-mediated Lys4 methylation of histone H3 (H3K4me), associated with positive regulation of transcription. In GBM CSCs, WDR5 inhibitors blocked WRAD complex assembly and reduced H3K4 trimethylation and expression of genes involved in CSC-relevant oncogenic pathways. H3K4me3 peaks lost with WDR5 inhibitor treatment occurred disproportionally on POU transcription factor motifs, including the POU5F1(OCT4)::SOX2 motif. Use of a SOX2/OCT4 reporter demonstrated that WDR5 inhibitor treatment diminished cells with high reporter activity. Furthermore, WDR5 inhibitor treatment and WDR5 knockdown altered the stem cell state, disrupting CSC in vitro growth and self-renewal, as well as in vivo tumor growth. These findings highlight the role of WDR5 and the WRAD complex in maintaining the CSC state and provide a rationale for therapeutic development of WDR5 inhibitors for GBM and other advanced cancers.

Glioblastoma stem cells generate vascular pericytes to support vessel function and tumor growth.

Glioblastomas (GBMs) are highly vascular and lethal brain tumors that display cellular hierarchies containing self-renewing tumorigenic glioma stem cells (GSCs). Because GSCs often reside in perivascular niches and may undergo mesenchymal differentiation, we interrogated GSC potential to generate vascular pericytes. Here, we show that GSCs give rise to pericytes to support vessel function and tumor growth. In vivo cell lineage tracing with constitutive and lineage-specific fluorescent reporters demonstrated that GSCs generate the majority of vascular pericytes. Selective elimination of GSC-derived pericytes disrupts the neovasculature and potently inhibits tumor growth. Analysis of human GBM specimens showed that most pericytes are derived from neoplastic cells. GSCs are recruited toward endothelial cells via the SDF-1/CXCR4 axis and are induced to become pericytes predominantly by transforming growth factor ß. Thus, GSCs contribute to vascular pericytes that may actively remodel perivascular niches. Therapeutic targeting of GSC-derived pericytes may effectively block tumor progression and improve antiangiogenic therapy.

A TOX-ic axis of epigenetic stem cell maintenance and chemoresistance in colon cancer.

Cancer stem cells drive tumor growth and survival via self-renewal and therapeutic resistance, but the upstream mechanisms are not well defined. In this issue of PLOS Biology, a study in colon cancer reveals a new signalling network that links epigenetic regulation to these phenotypes.

Glioma stem cell proliferation and tumor growth are promoted by nitric oxide synthase-2.

Malignant gliomas are aggressive brain tumors with limited therapeutic options, and improvements in treatment require a deeper molecular understanding of this disease. As in other cancers, recent studies have identified highly tumorigenic subpopulations within malignant gliomas, known generally as cancer stem cells. Here, we demonstrate that glioma stem cells (GSCs) produce nitric oxide via elevated nitric oxide synthase-2 (NOS2) expression. GSCs depend on NOS2 activity for growth and tumorigenicity, distinguishing them from non-GSCs and normal neural progenitors. Gene expression profiling identified many NOS2-regulated genes, including the cell-cycle inhibitor cell division autoantigen-1 (CDA1). Further, high NOS2 expression correlates with decreased survival in human glioma patients, and NOS2 inhibition slows glioma growth in a murine intracranial model. These data provide insight into how GSCs are mechanistically distinct from their less tumorigenic counterparts and suggest that NOS2 inhibition may be an efficacious approach to treating this devastating disease.

Metaorganismal choline metabolism shapes olfactory perception.

Microbes living in the intestine can regulate key signaling processes in the central nervous system that directly impact brain health. This gut-brain signaling axis is partially mediated by microbe-host-dependent immune regulation, gut-innervating neuronal communication, and endocrine-like small molecule metabolites that originate from bacteria to ultimately cross the blood-brain barrier. Given the mounting evidence of gut-brain crosstalk, a new therapeutic approach of "psychobiotics" has emerged, whereby strategies designed to primarily modify the gut microbiome have been shown to improve mental health or slow neurodegenerative diseases. Diet is one of the most powerful determinants of gut microbiome community structure, and dietary habits are associated with brain health and disease. Recently, the metaorganismal (i.e., diet-microbe-host) trimethylamine N-oxide (TMAO) pathway has been linked to the development of several brain diseases including Alzheimer's, Parkinson's, and ischemic stroke. However, it is poorly understood how metaorganismal TMAO production influences brain function under normal physiological conditions. To address this, here we have reduced TMAO levels by inhibiting gut microbe-driven choline conversion to trimethylamine (TMA), and then performed comprehensive behavioral phenotyping in mice. Unexpectedly, we find that TMAO is particularly enriched in the murine olfactory bulb, and when TMAO production is blunted at the level of bacterial choline TMA lyase (CutC/D), olfactory perception is altered. Taken together, our studies demonstrate a previously underappreciated role for the TMAO pathway in olfactory-related behaviors.

Cell surface CD55 traffics to the nucleus leading to cisplatin resistance and stemness by inducing PRC2 and H3K27 trimethylation on chromatin in ovarian cancer.

BACKGROUND: Platinum resistance is the primary cause of poor survival in ovarian cancer (OC) patients. Targeted therapies and biomarkers of chemoresistance are critical for the treatment of OC patients. Our previous studies identified cell surface CD55, a member of the complement regulatory proteins, drives chemoresistance and maintenance of cancer stem cells (CSCs). CSCs are implicated in tumor recurrence and metastasis in multiple cancers. METHODS: Protein localization assays including immunofluorescence and subcellular fractionation were used to identify CD55 at the cell surface and nucleus of cancer cells. Protein half-life determinations were used to compare cell surface and nuclear CD55 stability. CD55 deletion mutants were generated and introduced into cancer cells to identify the nuclear trafficking code, cisplatin sensitivity, and stem cell frequency that were assayed using in vitro and in vivo models. Detection of CD55 binding proteins was analyzed by immunoprecipitation followed by mass spectrometry. Target pathways activated by CD55 were identified by RNA sequencing. RESULTS: CD55 localizes to the nucleus of a subset of OC specimens, ascites from chemoresistant patients, and enriched in chemoresistant OC cells. We determined that nuclear CD55 is glycosylated and derived from the cell surface pool of CD55. Nuclear localization is driven by a trafficking code containing the serine/threonine (S/T) domain of CD55. Nuclear CD55 is necessary for cisplatin resistance, stemness, and cell proliferation in OC cells. CD55 S/T domain is necessary for nuclear entry and inducing chemoresistance to cisplatin in both in vitro and in vivo models. Deletion of the CD55 S/T domain is sufficient to sensitize chemoresistant OC cells to cisplatin. In the nucleus, CD55 binds and attenuates the epigenetic regulator and tumor suppressor ZMYND8 with a parallel increase in H3K27 trimethylation and members of the Polycomb Repressive Complex 2. CONCLUSIONS: For the first time, we show CD55 localizes to the nucleus in OC and promotes CSC and chemoresistance. Our studies identify a therapeutic mechanism for treating platinum resistant ovarian cancer by blocking CD55 nuclear entry.

Multimodal single-cell/nucleus RNA sequencing data analysis uncovers molecular networks between disease-associated microglia and astrocytes with implications for drug repurposing in Alzheimer's disease.

Because disease-associated microglia (DAM) and disease-associated astrocytes (DAA) are involved in the pathophysiology of Alzheimer's disease (AD), we systematically identified molecular networks between DAM and DAA to uncover novel therapeutic targets for AD. Specifically, we develop a network-based methodology that leverages single-cell/nucleus RNA sequencing data from both transgenic mouse models and AD patient brains, as well as drug-target network, metabolite-enzyme associations, the human protein-protein interactome, and large-scale longitudinal patient data. Through this approach, we find both common and unique gene network regulators between DAM (i.e., PAK1, MAPK14, and CSF1R) and DAA (i.e., NFKB1, FOS, and JUN) that are significantly enriched by neuro-inflammatory pathways and well-known genetic variants (i.e., BIN1). We identify shared immune pathways between DAM and DAA, including Th17 cell differentiation and chemokine signaling. Last, integrative metabolite-enzyme network analyses suggest that fatty acids and amino acids may trigger molecular alterations in DAM and DAA. Combining network-based prediction and retrospective case-control observations with 7.2 million individuals, we identify that usage of fluticasone (an approved glucocorticoid receptor agonist) is significantly associated with a reduced incidence of AD (hazard ratio [HR] = 0.86, 95% confidence interval [CI] 0.83-0.89, P < 1.0 × 10-8). Propensity score-stratified cohort studies reveal that usage of mometasone (a stronger glucocorticoid receptor agonist) is significantly associated with a decreased risk of AD (HR = 0.74, 95% CI 0.68-0.81, P < 1.0 × 10-8) compared to fluticasone after adjusting age, gender, and disease comorbidities. In summary, we present a network-based, multimodal methodology for single-cell/nucleus genomics-informed drug discovery and have identified fluticasone and mometasone as potential treatments in AD.

Iron inhibits glioblastoma cell migration and polarization.

Glioblastoma is one of the deadliest malignancies facing modern oncology today. The ability of glioblastoma cells to diffusely spread into neighboring healthy brain makes complete surgical resection nearly impossible and contributes to the recurrent disease faced by most patients. Although research into the impact of iron on glioblastoma has addressed proliferation, there has been little investigation into how cellular iron impacts the ability of glioblastoma cells to migrate-a key question, especially in the context of the diffuse spread observed in these tumors. Herein, we show that increasing cellular iron content results in decreased migratory capacity of human glioblastoma cells. The decrease in migratory capacity was accompanied by a decrease in cellular polarization in the direction of movement. Expression of CDC42, a Rho GTPase that is essential for both cellular migration and establishment of polarity in the direction of cell movement, was reduced upon iron treatment. We then analyzed a single-cell RNA-seq dataset of human glioblastoma samples and found that cells at the tumor periphery had a gene signature that is consistent with having lower levels of cellular iron. Altogether, our results suggest that cellular iron content is impacting glioblastoma cell migratory capacity and that cells with higher iron levels exhibit reduced motility.

Microglial immunometabolism endophenotypes contribute to sex difference in Alzheimer's disease.

INTRODUCTION: The molecular mechanisms that contribute to sex differences, in particular female predominance, in Alzheimer's disease (AD) prevalence, symptomology, and pathology, are incompletely understood. METHODS: To address this problem, we investigated cellular metabolism and immune responses ("immunometabolism endophenotype") across AD individuals as a function of sex with diverse clinical diagnosis of cognitive status at death (cogdx), Braak staging, and Consortium to Establish a Registry for AD (CERAD) scores using human cortex metabolomics and transcriptomics data from the Religious Orders Study / Memory and Aging Project (ROSMAP) cohort. RESULTS: We identified sex-specific metabolites, immune and metabolic genes, and pathways associated with the AD diagnosis and progression. We identified female-specific elevation in glycerophosphorylcholine and N-acetylglutamate, which are AD inflammatory metabolites involved in interleukin (IL)-17 signaling, C-type lectin receptor, interferon signaling, and Toll-like receptor pathways. We pinpointed distinct microglia-specific immunometabolism endophenotypes (i.e., lipid- and amino acid-specific IL-10 and IL-17 signaling pathways) between female and male AD subjects. In addition, female AD subjects showed evidence of diminished excitatory neuron and microglia communications via glutamate-mediated immunometabolism. DISCUSSION: Our results point to new understanding of the molecular basis for female predominance in AD, and warrant future independent validations with ethnically diverse patient cohorts to establish a likely causal relationship of microglial immunometabolism in the sex differences in AD. HIGHLIGHTS: Sex-specific immune metabolites, gene networks and pathways, are associated with Alzheimer's disease pathogenesis and disease progression. Female AD subjects exhibit microglial immunometabolism endophenotypes characterized by decreased glutamate metabolism and elevated interleukin-10 pathway activity. Female AD subjects showed a shift in glutamate-mediated cell-cell communications between excitatory neurons to microglia and astrocyte.

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 stem cells in glioblastoma.

Tissues with defined cellular hierarchies in development and homeostasis give rise to tumors with cellular hierarchies, suggesting that tumors recapitulate specific tissues and mimic their origins. Glioblastoma (GBM) is the most prevalent and malignant primary brain tumor and contains self-renewing, tumorigenic cancer stem cells (CSCs) that contribute to tumor initiation and therapeutic resistance. As normal stem and progenitor cells participate in tissue development and repair, these developmental programs re-emerge in CSCs to support the development and progressive growth of tumors. Elucidation of the molecular mechanisms that govern CSCs has informed the development of novel targeted therapeutics for GBM and other brain cancers. CSCs are not self-autonomous units; rather, they function within an ecological system, both actively remodeling the microenvironment and receiving critical maintenance cues from their niches. To fulfill the future goal of developing novel therapies to collapse CSC dynamics, drawing parallels to other normal and pathological states that are highly interactive with their microenvironments and that use developmental signaling pathways will be beneficial.

Roadmap on plasticity and epigenetics in cancer.

The role of plasticity and epigenetics in shaping cancer evolution and response to therapy has taken center stage with recent technological advances including single cell sequencing. This roadmap article is focused on state-of-the-art mathematical and experimental approaches to interrogate plasticity in cancer, and addresses the following themes and questions: is there a formal overarching framework that encompasses both non-genetic plasticity and mutation-driven somatic evolution? How do we measure and model the role of the microenvironment in influencing/controlling non-genetic plasticity? How can we experimentally study non-genetic plasticity? Which mathematical techniques are required or best suited? What are the clinical and practical applications and implications of these concepts?

Plasma IL13Rα2 as a novel liquid biopsy biomarker for glioblastoma.

PURPOSE: Glioblastoma (GBM) is the most common and deadliest brain tumor with unrelenting and rapid disease progression. The standard of care for GBM is surgical excision followed by radiation with concurrent and adjuvant temozolomide-centered chemotherapy (TMZ). Treatment failure and resistance is the rule and despite advances in imaging technology, early detection of treatment failure or impending resistance remains a challenge. There is a dire, unmet, need in clinical practice for minimally-invasive diagnostic tools to enable timely understanding of disease progression and treatment response. Here, we aim to address this clinical need by leveraging a unique characteristic of GBM: the overexpression of the α2 variant of the IL-13 receptor in over 75% of GBM tumors. METHODS: In this study we examined patients with primary GBM from Penn State and Cleveland Clinic compared to healthy controls. RESULTS: IL13Rα2 was detectable in plasma of GBM patients using ELISA but detection could be optimized by PEG precipitation to enrich for extracellular vesicles (EVs). Patients with GBM had elevated levels of plasma IL13Rα2, which correlated to levels of this receptor in the tumor tissue. Elevated plasma levels of IL13Rα2 predicted longer overall survival (OS) (19.8 vs. 13.2 months). Similarly, detection of IL13Rα2 + cells in tumor tissue also predicted longer OS (22.1 vs. 12.2 months). CONCLUSION: These findings strongly suggest that expression of the IL13Rα2 receptor confer survival advantage in GBM patients, which can be determined through a minimally-invasive liquid biopsy. Detection of plasma IL13Rα2 can also be used to select GBM patients for targeted tumor therapy.

Pharmacokinetic and brain distribution study of an anti-glioblastoma agent in mice by HPLC-MS/MS.

Previously compound I showed great anti-glioblastoma activity without toxicity in a mouse xenograft study. In this study, a sensitive and rapid high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) method was developed and validated to investigate the pharmacokinetics and brain distribution of compound I in mice. The protein precipitation method was applied to extract the compound from mouse plasma and brain homogenates, and it was then separated using a Kinetex C18 column with a mobile phase consisting of acetonitrile-0.1% formic acid water (50:50, v/v). The analytes were detected with multiple reaction monitoring for the quantitative response of the compounds. The inter- and intra-day precisions were <8.29 and 3.85%, respectively, and the accuracy range was within ±7.33%. The method was successfully applied to evaluate the pharmacokinetics of compound I in mouse plasma and brain tissue. The peak concentration in plasma was achieved within 1 h. The apparent elimination half-life was 4.06 h. The peak concentration of compound I in brain tissue was 0.88 µg/g. The results indicated that compound I was rapidly distributed and could cross the blood-brain barrier. The pharmacokinetic profile summarized provides valuable information for the further investigation of compound I as a potential anti-glioblastoma agent.

Independently validated sex-specific nomograms for predicting survival in patients with newly diagnosed glioblastoma: NRG Oncology RTOG 0525 and 0825.

BACKGROUND/PURPOSE: Glioblastoma (GBM) is the most common primary malignant brain tumor. Sex has been shown to be an important prognostic factor for GBM. The purpose of this study was to develop and independently validate sex-specific nomograms for estimation of individualized GBM survival probabilities using data from 2 independent NRG Oncology clinical trials. METHODS: This analysis included information on 752 (NRG/RTOG 0525) and 599 (NRG/RTOG 0825) patients with newly diagnosed GBM. The Cox proportional hazard models by sex were developed using NRG/RTOG 0525 and significant variables were identified using a backward selection procedure. The final selected models by sex were then independently validated using NRG/RTOG 0825. RESULTS: Final nomograms were built by sex. Age at diagnosis, KPS, MGMT promoter methylation and location of tumor were common significant predictors of survival for both sexes. For both sexes, tumors in the frontal lobes had significantly better survival than tumors of multiple sites. Extent of resection, and use of corticosteroids were significant predictors of survival for males. CONCLUSIONS: A sex specific nomogram that assesses individualized survival probabilities (6-, 12- and 24-months) for patients with GBM could be more useful than estimation of overall survival as there are factors that differ between males and females. A user friendly online application can be found here- https://npatilshinyappcalculator.shinyapps.io/SexDifferencesInGBM/ .

The Clinical Impact of Cancer Stem Cells.

Patients with cancer can go though many stages in their disease, including diagnosis, recurrence, metastasis, and treatment failure. Cancer stem cells (CSCs) are a subgroup of cells within tumors that may explain the mechanism by which tumors recur and progress. CSCs can both self-renew and produce progenitor cells of more differentiated cancer cells as well as heterogeneously demonstrate resistance and the abilities to migrate and metastasize. These "stemness" characteristics are often the result of dysregulation of one or more pathways, which can be detected by various biomarkers. Although there has been considerable laboratory research conducted on CSCs, its relevance to the practicing oncologist may seem questionable. We sought to determine the clinical impact of CSCs on patients. A systematic literature search was conducted to identify analyses containing survival information based on the expression of known CSC biomarkers in any cancer. Overall, 234 survival analyses were identified, of which 82% reported that high expression of CSC biomarker(s) resulted in poor overall survival and/or disease-free survival compared with low or no expression of the biomarker. Elevated stemness biomarker levels were also associated with decreased tumor differentiation, altered TNM stage, and increased metastasis. This analysis would suggest that CSCs have a clinical impact on patients and that practicing oncologists need to start considering incorporating CSC-targeting therapies into their patients' treatment regimens. IMPLICATIONS FOR PRACTICE: Cancer stem cells (CSCs) may occur at any stage of cancer and are implicated in the occurrence of resistance, recurrence, and metastasis. A systematic literature analysis has shown that the presence of CSCs, identified via the upregulation of stemness pathway biomarkers, results in reduced survival across all cancers studied. Several CSC-targeting agents are currently approved, and several others are in clinical trials. Future treatment regimens will likely include CSC-targeting agents to enable the elimination of these holdouts to current therapies.

Poly(ADP-Ribose) Polymerase Inhibition Sensitizes Colorectal Cancer-Initiating Cells to Chemotherapy.

Colorectal cancer (CRC) remains a leading killer in the U.S. with resistance to treatment as the largest hurdle to cure. Colorectal cancer-initiating cells (CICs) are a self-renewing tumor population that contribute to tumor relapse. Here, we report that patient-derived CICs display relative chemoresistance compared with differentiated progeny. In contrast, conventional cell lines failed model therapeutic resistance. CICs preferentially repaired chemotherapy-induced DNA breaks, prompting us to interrogate DNA damage pathways against which pharmacologic inhibitors have been developed. We found that CICs critically depended on the key single-strand break repair mediator, poly(ADP-ribose) polymerase (PARP), to survive treatment with standard-of-care chemotherapy. Small molecule PARP inhibitors (PARPi) sensitized CICs to chemotherapy and reduced chemotherapy-treated CIC viability, self-renewal, and DNA damage repair. Although PARPi monotherapy failed to kill CICs, combined PARPi therapy with chemotherapy attenuated tumor growth in vivo. Clinical significance of PARPi for CRC patients was supported by elevated PARP levels in colorectal tumors compared with normal colon, with further increases in metastases. Collectively, our results suggest that PARP inhibition serves as a point of fragility for CICs by augmenting therapeutic efficacy of chemotherapy. Stem Cells 2019;37:42-53.

Platelet-derived growth factor receptors differentially inform intertumoral and intratumoral heterogeneity.

Growth factor-mediated proliferation and self-renewal maintain tissue-specific stem cells and are frequently dysregulated in cancers. Platelet-derived growth factor (PDGF) ligands and receptors (PDGFRs) are commonly overexpressed in gliomas and initiate tumors, as proven in genetically engineered models. While PDGFRα alterations inform intertumoral heterogeneity toward a proneural glioblastoma (GBM) subtype, we interrogated the role of PDGFRs in intratumoral GBM heterogeneity. We found that PDGFRα is expressed only in a subset of GBMs, while PDGFRß is more commonly expressed in tumors but is preferentially expressed by self-renewing tumorigenic GBM stem cells (GSCs). Genetic or pharmacological targeting of PDGFRß (but not PDGFRα) attenuated GSC self-renewal, survival, tumor growth, and invasion. PDGFRß inhibition decreased activation of the cancer stem cell signaling node STAT3, while constitutively active STAT3 rescued the loss of GSC self-renewal caused by PDGFRß targeting. In silico survival analysis demonstrated that PDGFRB informed poor prognosis, while PDGFRA was a positive prognostic factor. Our results may explain mixed clinical responses of anti-PDGFR-based approaches and suggest the need for integration of models of cancer as an organ system into development of cancer therapies.

High-Throughput Automated Single-Cell Imaging Analysis Reveals Dynamics of Glioblastoma Stem Cell Population During State Transition.

Cancer stem cells (CSCs) are a heterogeneous and dynamic self-renewing population that stands at the top of tumor cellular hierarchy and contribute to tumor recurrence and therapeutic resistance. As methods of CSC isolation and functional interrogation advance, there is a need for a reliable and accessible quantitative approach to assess heterogeneity and state transition dynamics in CSCs. We developed a high-throughput automated single cell imaging analysis (HASCIA) approach for the quantitative assessment of protein expression with single-cell resolution and applied the method to investigate spatiotemporal factors that influence CSC state transition using glioblastoma (GBM) CSCs (GSCs) as a model system. We were able to validate the quantitative nature of this approach through comparison of the protein expression levels determined by HASCIA to those determined by immunoblotting. A virtue of HASCIA was exemplified by detection of a subpopulation of SOX2-low cells, which expanded in fraction size during state transition. HASCIA also revealed that GSCs were committed to loose stem cell state at an earlier time point than the average SOX2 level decreased. Functional assessment of stem cell frequency in combination with the quantification of SOX2 expression by HASCIA defined a stable cutoff of SOX2 expression level for stem cell state. We also developed an approach to assess local cell density and found that denser monolayer areas possess higher average levels of SOX2, higher cell diversity, and a presence of a sub-population of slowly proliferating SOX2-low GSCs. HASCIA is an open source software that facilitates understanding the dynamics of heterogeneous cell population such as that of GSCs and their progeny. It is a powerful and easy-to-use image analysis and statistical analysis tool available at https://hascia.lerner.ccf.org. © 2019 International Society for Advancement of Cytometry.