Adhering towards tumorigenicity: altered adhesion mechanisms in glioblastoma cancer stem cells

The STEMRI trial: Magnetic resonance spectroscopy imaging can define tumor areas enriched in glioblastoma stem-like cells

Despite maximally safe resection of the magnetic resonance imaging (MRI)-defined contrast-enhanced (CE) central tumor area and chemoradiotherapy, most patients with glioblastoma (GBM) relapse within a year in peritumoral FLAIR regions. Magnetic resonance spectroscopy imaging (MRSI) can discriminate metabolic tumor areas with higher recurrence potential as CNI+ regions (choline/N-acetyl-aspartate index >2) can predict relapse sites. As relapses are mainly imputed to glioblastoma stem-like cells (GSCs), CNI+ areas might be GSC enriched. In this prospective trial, 16 patients with GBM underwent MRSI/MRI before surgery/chemoradiotherapy to investigate GSC content in CNI-/+ biopsies from CE/FLAIR. Biopsy and derived-GSC characterization revealed a FLAIR/CNI+ sample enrichment in GSC and in gene signatures related to stemness, DNA repair, adhesion/migration, and mitochondrial bioenergetics. FLAIR/CNI+ samples generate GSC-enriched neurospheres faster than FLAIR/CNI-. Parameters assessing biopsy GSC content and time-to-neurosphere formation in FLAIR/CNI+ were associated with worse patient outcome. Preoperative MRI/MRSI would certainly allow better resection and targeting of FLAIR/CNI+ areas, as their GSC enrichment can predict worse outcomes.

Bradykinin B1 Receptor Affects Tumor-Associated Macrophage Activity and Glioblastoma Progression

Bradykinin is a small active peptide and is considered an inflammatory mediator in several pathological conditions. Bradykinin exerts its effects by coupling to its receptors, including bradykinin B1 (B1R) and bradykinin B2. B1R has been implicated in the development of various cancers. Our previous study reported that B1R promoted glioblastoma (GBM) development by supporting the migration and invasion of GBM cells. However, the mechanisms underlying the effects of B1R on tumor-associated macrophages (TAMs) and GBM progression remain unknown. Accordingly, to explore the regulatory effects of B1R overexpression (OE) in GBM on tumor-associated immune cells and tumor progression, we constructed a B1R wild-type plasmid and developed a B1R OE model. The results reveal that B1R OE in GBM promoted the expression of ICAM-1 and VCAM-1-cell adhesion molecules-in GBM. Moreover, B1R OE enhanced GBM cell migration ability and monocyte attachment. B1R also regulated the production of the protumorigenic cytokines and chemokines IL-6, IL-8, CXCL11, and CCL5 in GBM, which contributed to tumor progression. We additionally noted that B1R OE in GBM increased the expression of CD68 in TAMs. Furthermore, B1R OE reduced the level of reactive oxygen species in GBM cells by upregulating heme oxygenase-1, an endogenous antioxidant protein, thereby protecting GBM cells from oxidative stress. Notably, B1R OE upregulated the expression of programmed death-ligand 1 in both GBM cells and macrophages, thus providing resistance against T-cell response. B1R OE in GBM also promoted tumor growth and reduced survival rates in an intracranial xenograft mouse model. These results indicate that B1R expression in GBM promotes TAM activity and modulates GBM progression. Therefore, B1R could be an effective target for therapeutic methods in GBM.

mTOR Inhibition Is Effective against Growth, Survival and Migration, but Not against Microglia Activation in Preclinical Glioma Models

Initially introduced in therapy as immunosuppressants, the selective inhibitors of mTORC1 have been approved for the treatment of solid tumors. Novel non-selective inhibitors of mTOR are currently under preclinical and clinical developments in oncology, attempting to overcome some limitations associated with selective inhibitors, such as the development of tumor resistance. Looking at the possible clinical exploitation in the treatment of glioblastoma multiforme, in this study we used the human glioblastoma cell lines U87MG, T98G and microglia (CHME-5) to compare the effects of a non-selective mTOR inhibitor, sapanisertib, with those of rapamycin in a large array of experimental paradigms, including (i) the expression of factors involved in the mTOR signaling cascade, (ii) cell viability and mortality, (iii) cell migration and autophagy, and (iv) the profile of activation in tumor-associated microglia. We could distinguish between effects of the two compounds that were overlapping or similar, although with differences in potency and or/time-course, and effects that were diverging or even opposite. Among the latter, especially relevant is the difference in the profile of microglia activation, with rapamycin being an overall inhibitor of microglia activation, whereas sapanisertib was found to induce an M2-profile, which is usually associated with poor clinical outcomes.

Integrative analysis of cell adhesion molecules in glioblastoma identified prostaglandin F2 receptor inhibitor (PTGFRN) as an essential gene

Background

Glioblastoma (GBM) is the most common primary malignant brain tumor in adults exhibiting infiltration into surrounding tissues, recurrence, and resistance to therapy. GBM infiltration is accomplished by many deregulated factors such as cell adhesion molecules (CAMs), which are membrane proteins that participate in cell-cell and cell-ECM interactions to regulate survival, proliferation, migration, and stemness.

Methods

A comprehensive bioinformatics analysis of CAMs (n = 518) in multiple available datasets revealed genetic and epigenetic alterations among CAMs in GBM. Univariate Cox regression analysis using TCGA dataset identified 127 CAMs to be significantly correlated with survival. The poor prognostic indicator PTGFRN was chosen to study its role in glioma. Silencing of PTGFRN in glioma cell lines was achieved by the stable expression of short hairpin RNA (shRNA) against the PTGFRN gene. PTGFRN was silenced and performed cell growth, migration, invasion, cell cycle, and apoptosis assays. Neurosphere and limiting dilution assays were also performed after silencing of PTGFRN in GSCs.

Results

Among the differentially regulated CAMs (n = 181, 34.9%), major proportion of them were found to be regulated by miRNAs (n = 95, 49.7%) followed by DNA methylation (n = 32, 16.7%), and gene copy number variation (n = 12, 6.2%). We found that PTGFRN to be upregulated in GBM tumor samples and cell lines with a significant poor prognostic correlation with patient survival. Silencing PTGFRN diminished cell growth, colony formation, anchorage-independent growth, migration, and invasion and led to cell cycle arrest and induction of apoptosis. At the mechanistic level, silencing of PTGFRN reduced pro-proliferative and promigratory signaling pathways such as ERK, AKT, and mTOR. PTGFRN upregulation was found to be due to the loss of its promoter methylation and downregulation of miR-137 in GBM. PTGFRN was also found to be higher in glioma stem-like cells (GSCs) than the matched differentiated glioma cells (DGCs) and is required for GSC growth and survival. Silencing of PTGFRN in GSCs reduced transcript levels of reprogramming factors (Olig2, Pou3f2, Sall2, and Sox2).

Conclusion

In this study, we provide a comprehensive overview of the differential regulation of CAMs and the probable causes for their deregulation in GBM. We also establish an oncogenic role of PTGFRN and its regulation by miR-137 in GBM, thus signifying it as a potential therapeutic target.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09682-2.

Systems-Level Proteomics Evaluation of Microglia Response to Tumor-Supportive Anti-Inflammatory Cytokines

Background

Microglia safeguard the CNS against injuries and pathogens, and in the presence of certain harmful stimuli are capable of inducing a disease-dependent inflammatory response. When exposed to anti-inflammatory cytokines, however, these cells possess the ability to switch from an inflammatory to an immunosuppressive phenotype. Cancer cells exploit this property to evade the immune system, and elicit an anti-inflammatory microenvironment that facilitates tumor attachment and growth.

Objective

The tumor-supportive biological processes that are activated in microglia cells in response to anti-inflammatory cytokines released from cancer cells were explored with mass spectrometry and proteomic technologies.

Methods

Serum-depleted and non-depleted human microglia cells (HMC3) were treated with a cocktail of IL-4, IL-13, IL-10, TGFβ, and CCL2. The cellular protein extracts were analyzed by LC-MS/MS. Using functional annotation clustering tools, statistically significant proteins that displayed a change in abundance between cytokine-treated and non-treated cells were mapped to their biological networks and pathways.

Results

The proteomic analysis of HMC3 cells enabled the identification of ~10,000 proteins. Stimulation with anti-inflammatory cytokines resulted in the activation of distinct, yet integrated clusters of proteins that trigger downstream a number of tumor-promoting biological processes. The observed changes could be classified into four major categories, i.e., mitochondrial gene expression, ECM remodeling, immune response, and impaired cell cycle progression. Intracellular immune activation was mediated mainly by the transducers of MAPK, STAT, TGFβ, NFKB, and integrin signaling pathways. Abundant collagen formation along with the expression of additional receptors, matrix components, growth factors, proteases and protease inhibitors, was indicative of ECM remodeling processes supportive of cell-cell and cell-matrix adhesion. Overexpression of integrins and their modulators was reflective of signaling processes that link ECM reorganization with cytoskeletal re-arrangements supportive of cell migration. Antigen processing/presentation was represented by HLA class I histocompatibility antigens, and correlated with upregulated proteasomal subunits, vesicular/viral transport, and secretory processes. Immunosuppressive and proangiogenic chemokines, as well as anti-angiogenic factors, were detectable in low abundance. Pronounced pro-inflammatory, chemotactic or phagocytic trends were not observed, however, the expression of certain receptors, signaling and ECM proteins indicated the presence of such capabilities.

Conclusions

Comprehensive proteomic profiling of HMC3 cells stimulated with anti-inflammatory cytokines revealed a spectrum of microglia phenotypes supportive of cancer development in the brain via microenvironment-dependent biological mechanisms.

CHI3L2 Is a Novel Prognostic Biomarker and Correlated With Immune Infiltrates in Gliomas

CHI3L2 (Chitinase-3-Like Protein 2) is a member of chitinase-like proteins (CLPs), which belong to the glycoside hydrolase 18 family. Its homologous gene, CHI3L1, has been extensively studied in various tumors and has been shown to be related to immune infiltration in breast cancer and glioblastoma. High CHI3L2 expression was reported to be associated with poor prognosis in breast cancer and renal cell carcinoma. However, the prognostic significance of CHI3L2 in glioma and its correlation between immune infiltration remains unclear. In this study, we examined 288 glioma samples by immunohistochemistry to find that CHI3L2 is expressed in tumor cells and macrophages in glioma tissues and highly expressed in glioblastoma and IDH wild-type gliomas. Relationships between CHI3L2 expression and clinical features (grade, age, Ki67 index, P53, PHH3 (mitotic figures), ATRX, TERTp, MGMTp, IDH, and 1p/19q co-deleted status) were evaluated. Kaplan-Meier survival was conducted to show high CHI3L2 expression in tumor cells (TC) and macrophage cells (MC) indicated poor prognosis in diffusely infiltrating glioma (DIG), lower-grade glioma (LGG), and IDH wild-type gliomas (IDH-wt). The overall survival time was higher in patients with dual-low CHI3L2 expression in TC and MC compared to those in patients with non-dual CHI3L2 expression and dual high expression in DIG and IDH wild-type gliomas. By univariate and multivariate analysis, we found that high CHI3L2 expression in tumor cells was an independent unfavorable prognostic factor in glioma patients. Moreover, we used two datasets (TCGA and CGGA) to verify the results of our study and explore the potential functional role of CHI3L2 by GO and KEGG analyses in gliomas. TIMER platform analysis indicated CHI3L2 expression was closely related to diverse marker genes of tumor immune infiltrating cells, including monocytes, TAMs, M1 macrophages, M2 macrophages, TGFβ1+ Treg and T cell exhaustion in GBM and LGG. Western Blot validated CHI3L2 is expressed in glioma cells and microglia cells. The results of flow cytometry showed that CHI3L2 induces the apoptosis of CD8+ T cells. In conclusion, these results demonstrate CHI3L2 is related to poor prognosis and immune infiltrates in gliomas, suggesting it may serve as a promising prognostic biomarker and represent a new target for glioma patients.

A retrospective analysis of GSE84010: Cell adhesion molecules might contribute to bevacizumab resistance in glioblastoma

Bevacizumab (BEV) is an anti-angiogenesis antibody which has shown favorable therapeutic effects on some solid tumors. However, many clinical trials showed that BEV could only improve PFS instead of OS in glioblastoma (GBM) patients. However, some studies indicate that specific molecular subtypes of GBM could still benefit from combination treatment of BEV and Stupp protocol. Through the subgroup analysis of GSE84010 dataset, we found the neural and proneural subgroup can benefit from the administration of BEV in terms of OS, which is statistically significant. The further KEGG pathway enrichment analysis showed cell adhesion molecules (CAMs) pathway was enriched, and the expression of ITGAM has a predictive value for prognosis. These findings can provide some hints for future administration of BEV in newly diagnosed GBM patients.

Recent Advances in Cancer Plasticity: Cellular Mechanisms, Surveillance Strategies, and Therapeutic Optimization

The processes of recurrence and metastasis, through which cancer relapses locally or spreads to distant sites in the body, accounts for more than 90% of cancer-related deaths. At present there are very few treatment options for patients at this stage of their disease. The main obstacle to successfully treat advanced cancer is the cells' ability to change in ways that make them resistant to treatment. Understanding the cellular mechanisms that mediate this cancer cell plasticity may lead to improved patient survival. Epigenetic reprogramming, together with tumor microenvironment, drives such dynamic mechanisms favoring tumor heterogeneity, and cancer cell plasticity. In addition, the development of new approaches that can report on cancer plasticity in their native environment have profound implications for studying cancer biology and monitoring tumor progression. We herein provide an overview of recent advancements in understanding the mechanisms regulating cell plasticity and current strategies for their monitoring and therapy management.

Cell adhesion in cancer: Beyond the migration of single cells

Homeostasis in healthy tissues strongly relies on cell-to-cell adhesion and cell-to-extracellular matrix interactions. For instance, normal epithelial cells maintain tissue structure by adhering to each other and to the extracellular matrix. The proteins that mediate these distinct interactions are collectively called cell adhesion molecules and are divided into four major groups: cadherins, integrins, selectins, and immunoglobulins. They not only physically anchor cells, but also critically integrate signaling between the extracellular microenvironment and cells. These signals include biochemical cues, as adhesion proteins can both act as ligand-activated receptors and activate mechanotransduction triggered by changes in the physical environment. Molecular mechanisms related to cell adhesion signaling have been extensively studied, especially because mutations and changes in expression of these proteins, particularly cadherins and integrins, are frequently associated with diseases ranging from developmental intellectual disability to cancer. In fact, two major hallmarks of cancer, loss of cell-to-cell adhesion and anchorage-independent growth, are both dependent on cell adhesion molecules. Despite many studies elucidating the relationships between malignant transformation and metastasis and cellular adhesion processes, several areas still await exploration. Here, we highlight recently discovered roles of adhesion molecules in collective cancer cell migration and discuss the utility of three-dimensional models in studying cell-cell adhesion. We also describe recent therapeutic approaches targeting adhesion molecules.

Mutant-allele tumor heterogeneity in malignant glioma effectively predicts neoplastic recurrence

Intra-tumor heterogeneity (ITH) is one of the most important causes of therapy resistance, which eventually leads to the poor outcomes observed in patients with glioma. Mutant-allele tumor heterogeneity (MATH) values are based on whole-exon sequencing and precisely reflect genetic ITH. However, the significance of MATH values in predicting glioma recurrence remains unclear. Information of patients with glioma was obtained from The Cancer Genome Atlas database. The present study calculated the MATH value for each patient, analyzed the distributions of MATH values in different subtypes and investigated the rates of clinical recurrence in patients with different MATH values. Gene enrichment and Cox regression analyses were performed to determine which factors influenced recurrence. A nomogram table was established to predict 1-, 2- and 5-year recurrence probabilities. MATH values were increased in patients with glioma with the wild-type isocitrate dehydrogenase (NADP⁽⁺⁾) (IDH)1/2 (IDH-wt) gene (P=0.001) and glioblastoma (GBM; P=0.001). MATH values were negatively associated with the 2- and 5-year recurrence-free survival (RFS) rates in patients with glioma, particularly in the IDH1/2-wt and GBM cohorts (P=0.001 and P=0.017, respectively). Furthermore, glioma cases with different MATH levels had distinct patterns of gene mutation frequencies and gene expression enrichment. Finally, a nomogram table that contained MATH values could be used to accurately predict the probabilities of the 1-, 2- and 5-year RFS of patients with glioma. In conclusion, the MATH value of a patient may be an independent predictor that influences glioma recurrence. The nomogram model presented in the current study was an appropriate method to predict 1-, 2- and 5-year RFS probabilities in patients with glioma.

Proteomic Analysis Implicates Vimentin in Glioblastoma Cell Migration

We previously showed lithium chloride (LiCl) and other inhibitors of glycogen synthase kinase-3 (GSK-3) including 6-bromo-indirubin-3-oxime (BIO), can block glioblastoma (GBM) cell migration. To investigate the mechanisms involved we used two-dimensional difference in-gel electrophoresis (2D-DIGE) and mass spectrometry to identify proteins altered after treatment of U251 GBM cells with 20 mM LiCl. Downregulation of the intermediate filament protein vimentin was the most significant change identified. Analysis of patient tumor samples revealed that vimentin is expressed abundantly in GBM, and is prognostic especially in lower grade tumors. Additionally, siRNA-mediated vimentin knockdown impaired GBM migration. Western blotting showed that treatment with LiCl or small molecule GSK-3 inhibitors led to the rapid downregulation of detergent soluble vimentin levels across a panel of GBM-derived cells. Fluorescence reactivation after photobleaching (FRAP) microscopy studies showed a significant reduction in the ability of the vimentin cytoskeleton to recover from photo-bleaching in the presence of LiCl or BIO. Biochemical studies revealed that GSK-3 and vimentin directly interact, and analysis of vimentin revealed a GSK-3 consensus phosphorylation site. We conclude that anti-migratory compounds with the ability to inhibit GSK-3 have effects on vimentin cytoskeletal dynamics, which may play a role in their anti-invasive activity.

The underlying molecular mechanism and potential drugs for treatment in papillary renal cell carcinoma: A study based on TCGA and Cmap datasets

Papillary renal cell carcinoma (PRCC) accounts for 15–20% of all kidney neoplasms and continually attracts attention due to the increase in the incidents in which it occurs. The molecular mechanism of PRCC remains unclear and the efficacy of drugs that treat PRCC lacks sufficient evidence in clinical trials. Therefore, it is necessary to investigate the underlying mechanism in the development of PRCC and identify additional potential anti-PRCC drugs for its treatment. The differently expressed genes (DEGs) of PRCC were identified, followed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses for functional annotation. Then, potential drugs for PRCC treatment were predicted by Connectivity Map (Cmap) based on DEGs. Furthermore, the latent function of query drugs in PRCC was explored by integrating drug-target, drug-pathway and drug-protein interactions. In total, 627 genes were screened as DEGs, and these DEGs were annotated using KEGG pathway analyses and were clearly associated with the complement and coagulation cascades, amongst others. Then, 60 candidate drugs, as predicted based on DEGs, were obtained from the Cmap database. Vorinostat was considered as the most promising drug for detailed discussion. Following protein-protein interaction (PPI) analysis and molecular docking, vorinostat was observed to interact with C3 and ANXN1 proteins, which are the upregulated hub genes and may serve as oncologic therapeutic targets in PRCC. Among the top 20 metabolic pathways, several significant pathways, such as complement and coagulation cascades and cell adhesion molecules, may greatly contribute to the development and progression of PRCC. Following the performance of the PPI network and molecular docking tests, vorinostat exhibited a considerable and promising application in PRCC treatment by targeting C3 and ANXN1.

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.

CD24 expression and stem-associated features define tumor cell heterogeneity and tumorigenic capacities in a model of carcinogenesis

Background

Most carcinomas are composed of heterogeneous populations of tumor cells with distinct and apparently stable phenotypic characteristics.

Methods

Using an in vitro model of carcinogenesis we aimed at experimentally elucidating the significance of heterogeneity in the expression of CD24, a marker frequently overexpressed in various cancers and correlated with poor prognosis.

Results

We show that CD24^Neg and CD24^Pos cells issued from the same tumorigenic cell line display striking differences in stem-related properties, expression of epithelial-mesenchymal transition/mesenchymal-epithelial transition markers, and tumorigenic capacity. Indeed, while CD24^Neg cells were as tumorigenic as the parental cell line, CD24^Pos cells, although unable to form tumors, were unexpectedly more mesenchymal, displayed enhanced stemness-related properties, and expressed a proinflammatory signature.

Conclusion

Our findings support the view that acquisition of stem-like cell, CD24-associated, attributes like migration, invasion, and plasticity by a tumor subpopulation is not necessarily related to local tumor growth but may be required for escaping the niche and colonizing distant sites.

Glioma Stem Cell Niches in Human Glioblastoma Are Periarteriolar

Survival of primary brain tumor (glioblastoma) patients is seriously hampered by glioma stem cells (GSCs) that are distinct therapy-resistant self-replicating pluripotent cancer cells. GSCs reside in GSC niches, which are specific protective microenvironments in glioblastoma tumors. We have recently found that GSC niches are hypoxic periarteriolar, whereas in most studies, GSC niches are identified as hypoxic perivascular. The aim of this review is to critically evaluate the literature on perivascular GSC niches to establish whether these are periarteriolar, pericapillary, perivenular, and/or perilymphatic. We found six publications showing images of human glioblastoma tissue containing perivascular GSC niches without any specification of the vessel type. However, it is frequently assumed that these vessels are capillaries which are exchange vessels, whereas arterioles and venules are transport vessels. Closer inspection of the figures of these publications showed vessels that were not capillaries. Whether these vessels were arterioles or venules was difficult to determine in one case, but in the other cases, these were clearly arterioles and their perivascular niches were similar to the periarteriolar niches we have found. Therefore, we conclude that in human glioblastoma tumors, GSC niches are hypoxic periarteriolar and are structurally and functionally look-alikes of hematopoietic stem cell niches in the bone marrow.

The hypoxic peri-arteriolar glioma stem cell niche, an integrated concept of five types of niches in human glioblastoma

Glioblastoma is the most lethal primary brain tumor and poor survival of glioblastoma patients is attributed to the presence of glioma stem cells (GSCs). These therapy-resistant, quiescent and pluripotent cells reside in GSC niches, which are specific microenvironments that protect GSCs against radiotherapy and chemotherapy. We previously showed the existence of hypoxic peri-arteriolar GSC niches in glioblastoma tumor samples. However, other studies have described peri-vascular niches, peri-hypoxic niches, peri-immune niches and extracellular matrix niches of GSCs. The aim of this review was to critically evaluate the literature on these five different types of GSC niches. In the present review, we describe that the five niche types are not distinct from one another, but should be considered to be parts of one integral GSC niche model, the hypoxic peri-arteriolar GSC niche. Moreover, hypoxic peri-arteriolar GSC niches are structural and functional look-alikes of hematopoietic stem cell (HSC) niches in the bone marrow. GSCs are maintained in peri-arteriolar niches by the same receptor-ligand interactions as HSCs in bone marrow. Our concept should be rigidly tested in the near future and applied to develop therapies to expel and keep GSCs out of their protective niches to render them more vulnerable to standard therapies.

Functional Virtual Flow Cytometry: A Visual Analytic Approach for Characterizing Single-Cell Gene Expression Patterns

We presented a novel workflow for detecting distribution patterns in cell populations based on single-cell transcriptome study. With the fast adoption of single-cell analysis, a challenge to researchers is how to effectively extract gene features to meaningfully separate the cell population. Considering that coexpressed genes are often functionally or structurally related and the number of coexpressed modules is much smaller than the number of genes, our workflow uses gene coexpression modules as features instead of individual genes. Thus, when the coexpressed modules are summarized into eigengenes, not only can we interactively explore the distribution of cells but also we can promptly interpret the gene features. The interactive visualization is aided by a novel application of spatial statistical analysis to the scatter plots using a clustering index parameter. This parameter helps to highlight interesting 2D patterns in the scatter plot matrix (SPLOM). We demonstrated the effectiveness of the workflow using two large single-cell studies. In the Allen Brain scRNA-seq dataset, the visual analytics suggested a new hypothesis such as the involvement of glutamate metabolism in the separation of the brain cells. In a large glioblastoma study, a sample with a unique cell migration related signature was identified.