Single-cell sequencing maps gene expression to mutational phylogenies in PDGF- and EGF-driven gliomas

Emerging therapeutic strategies in glioblastsoma: drug repurposing, mechanisms of resistance, precision medicine, and technological innovations

Glioblastoma (GBM) is an aggressive Grade IV brain tumor with a poor prognosis. It results from genetic mutations, epigenetic changes, and factors within the tumor microenvironment (TME). Traditional treatments like surgery, radiotherapy, and chemotherapy provide limited survival benefits due to the tumor's heterogeneity and resistance mechanisms. This review examines novel approaches for treating GBM, focusing on repurposing existing medications such as antipsychotics, antidepressants, and statins for their potential anti-GBM effects. Advances in molecular profiling, including next-generation sequencing, artificial intelligence (AI), and nanotechnology-based drug delivery, are transforming GBM diagnosis and treatment. The TME, particularly GBM stem cells and immune evasion, plays a key role in therapeutic resistance. Integrating multi-omics data and applying precision medicine show promise, especially in combination therapies and immunotherapies, to enhance clinical outcomes. Addressing challenges such as drug resistance, targeting GBM stem cells, and crossing the blood-brain barrier is essential for improving treatment efficacy. While current treatments offer limited benefits, emerging strategies such as immunotherapies, precision medicine, and drug repurposing show significant potential. Technologies like liquid biopsies, AI-powered diagnostics, and nanotechnology could help overcome obstacles like the blood-brain barrier and GBM stem cells. Ongoing research into combination therapies, targeted drug delivery, and personalized treatments is crucial. Collaborative efforts and robust clinical trials are necessary to translate these innovations into effective therapies, offering hope for improved survival and quality of life for GBM patients.

Glioblastoma Cortical Organoids Recapitulate Cell-State Heterogeneity and Intercellular Transfer

Glioblastoma (GBM) is characterized by heterogeneous malignant cells that are functionally integrated within the neuroglial microenvironment. In this study, we model this ecosystem by growing GBM into long-term cultured human cortical organoids that contain the major neuroglial cell types found in the cerebral cortex. Single-cell RNA sequencing analysis suggests that, compared with matched gliomasphere models, GBM cortical organoids more faithfully recapitulate the diversity and expression programs of malignant cell states found in patient tumors. Additionally, we observe widespread transfer of GBM transcripts and GFP to nonmalignant cells in the organoids. Mechanistically, this transfer involves extracellular vesicles and is biased toward defined GBM cell states and astroglia cell types. These results extend previous GBM organoid modeling efforts and suggest widespread intercellular transfer in the GBM neuroglial microenvironment. Significance: Models that recapitulate intercellular communications in GBM are limited. In this study, we leverage GBM cortical organoids to characterize widespread mRNA and GFP transfer from malignant to nonmalignant cells in the GBM neuroglial microenvironment. This transfer involves extracellular vesicles, may contribute to reprogramming the microenvironment, and may extend to other cancer types. See related commentary by Shakya et al., p. 261.

Mapping the developmental trajectory of human astrocytes reveals divergence in glioblastoma

Glioblastoma (GBM) is defined by heterogeneous and resilient cell populations that closely reflect neurodevelopmental cell types. Although it is clear that GBM echoes early and immature cell states, identifying the specific developmental programmes disrupted in these tumours has been hindered by a lack of high-resolution trajectories of glial and neuronal lineages. Here we delineate the course of human astrocyte maturation to uncover discrete developmental stages and attributes mirrored by GBM. We generated a transcriptomic and epigenomic map of human astrocyte maturation using cortical organoids maintained in culture for nearly 2 years. Through this approach, we chronicled a multiphase developmental process. Our time course of human astrocyte maturation includes a molecularly distinct intermediate period that serves as a lineage commitment checkpoint upstream of mature quiescence. This intermediate stage acts as a site of developmental deviation separating IDH-wild-type neoplastic astrocyte-lineage cells from quiescent astrocyte populations. Interestingly, IDH1-mutant tumour astrocyte-lineage cells are the exception to this developmental perturbation, where immature properties are suppressed as a result of D-2-hydroxyglutarate oncometabolite exposure. We propose that this defiance is a consequence of IDH1-mutant-associated epigenetic dysregulation, and we identified biased DNA hydroxymethylation (5hmC) in maturation genes as a possible mechanism. Together, this study illustrates a distinct cellular state aberration in GBM astrocyte-lineage cells and presents developmental targets for experimental and therapeutic exploration.

Research progress and the prospect of using single-cell sequencing technology to explore the characteristics of the tumor microenvironment

In precision cancer therapy, addressing intra-tumor heterogeneity poses a significant obstacle. Due to the heterogeneity of each cell subtype and between cells within the tumor, the sensitivity and resistance of different patients to targeted drugs, chemotherapy, etc., are inconsistent. Concerning a specific tumor type, many feasible treatments or combinations can be used by specifically targeting the tumor microenvironment. To solve this problem, it is necessary to further study the tumor microenvironment. Single-cell sequencing techniques can dissect distinct tumor cell populations by isolating cells and using statistical computational methods. This technology may assist in the selection of targeted combination therapy, and the obtained cell subset information is crucial for the rational application of targeted therapy. In this review, we summarized the research and application advances of single-cell sequencing technology in the tumor microenvironment, including the most commonly used single-cell genomic and transcriptomic sequencing, and their future development direction was proposed. The application of single-cell sequencing technology has been expanded to include epigenomics, proteomics, metabolomics, and microbiome analysis. The integration of these different omics approaches has significantly advanced the development of single-cell multiomics sequencing technology. This innovative approach holds immense potential for various fields, such as biological research and medical investigations. Finally, we discussed the advantages and disadvantages of using single-cell sequencing to explore the tumor microenvironment.

Unveiling chemotherapy-induced immune landscape remodeling and metabolic reprogramming in lung adenocarcinoma by scRNA-sequencing

Chemotherapy is widely used to treat lung adenocarcinoma (LUAD) patients comprehensively. Considering the limitations of chemotherapy due to drug resistance and other issues, it is crucial to explore the impact of chemotherapy and immunotherapy on these aspects. In this study, tumor samples from nine LUAD patients, of which four only received surgery and five received neoadjuvant chemotherapy, were subjected to scRNA-seq analysis. In vitro and in vivo assays, including flow cytometry, immunofluorescence, Seahorse assay, and tumor xenograft models, were carried out to validate our findings. A total of 83,622 cells were enrolled for subsequent analyses. The composition of cell types exhibited high heterogeneity across different groups. Functional enrichment analysis revealed that chemotherapy drove significant metabolic reprogramming in tumor cells and macrophages. We identified two subtypes of macrophages: Anti-mac cells (CD45+CD11b+CD86+) and Pro-mac cells (CD45+CD11b+ARG +) and sorted them by flow cytometry. The proportion of Pro-mac cells in LUAD tissues increased significantly after neoadjuvant chemotherapy. Pro-mac cells promote tumor growth and angiogenesis and also suppress tumor immunity. Moreover, by analyzing the remodeling of T and B cells induced by neoadjuvant therapy, we noted that chemotherapy ignited a relatively more robust immune cytotoxic response toward tumor cells. Our study demonstrates that chemotherapy induces metabolic reprogramming within the tumor microenvironment of LUAD, particularly affecting the function and composition of immune cells such as macrophages and T cells. We believe our findings will offer insight into the mechanisms of drug resistance and provide novel therapeutic targets for LUAD in the future.

Impact of Hashimoto's thyroiditis on the tumor microenvironment in papillary thyroid cancer: insights from single-cell analysis

This study investigates the impact of Hashimoto's thyroiditis (HT), an autoimmune disorder, on the papillary thyroid cancer (PTC) microenvironment using a dataset of 140,456 cells from 11 patients. By comparing PTC cases with and without HT, we identify HT-specific cell populations (HASCs) and their role in creating a TSH-suppressive environment via mTE3, nTE0, and nTE2 thyroid cells. These cells facilitate intricate immune-stromal communication through the MIF-(CD74+CXCR4) axis, emphasizing immune regulation in the TSH context. In the realm of personalized medicine, our HASC-focused analysis within the TCGA-THCA dataset validates the utility of HASC profiling for guiding tailored therapies. Moreover, we introduce a novel, objective method to determine K-means clustering coefficients in copy number variation inference from bulk RNA-seq data, mitigating the arbitrariness in conventional coefficient selection. Collectively, our research presents a detailed single-cell atlas illustrating HT-PTC interactions, deepening our understanding of HT's modulatory effects on PTC microenvironments. It contributes to our understanding of autoimmunity-carcinogenesis dynamics and charts a course for discovering new therapeutic targets in PTC, advancing cancer genomics and immunotherapy research.

SAA suppresses α-PD-1 induced anti-tumor immunity by driving TH2 polarization in lung adenocarcinoma

Cancer stem cells (CSCs) are believed to be crucial in the initiation, progression, and recurrence of cancer. CSCs are also known to be more resistant to cancer treatments. However, the interaction between CSCs and the immune microenvironment is complex and not fully understood. In current study we used single cell RNA sequence (scRNA-Seq, public dataset) technology to identify the characteristic of CSCs. We found that the lung adenocarcinoma cancer stem population is highly inflammatory and remodels the tumor microenvironment by secreting inflammatory factors, specifically the acute phase protein serum amyloid A (SAA). Next, we developed an ex-vivo autologous patient-derived organoids (PDOs) and peripheral blood mononuclear cells (PBMCs) co-culture model to evaluate the immune biological impact of SAA. We found that SAA not only promotes chemoresistance by inducing cancer stem transformation, but also restricts anti-tumor immunity and promotes tumor fibrosis by driving type 2 immunity, and α-SAA neutralization antibody could restrict treatment resistant and tumor fibrosis. Mechanically, we found that the malignant phenotype induced by SAA is dependent on P2X7 receptor. Our data indicate that cancer stem cells secreted SAA have significant biological impact to promote treatment resistant and tumor fibrosis by driving cancer stemness transformation and type 2 immunity polarization via P2X7 receptor. Notably, α-SAA neutralization antibody shows therapeutic potential by restricting these malignant phenotypes.

Loss of SYNCRIP unleashes APOBEC-driven mutagenesis, tumor heterogeneity, and AR-targeted therapy resistance in prostate cancer

Tumor mutational burden and heterogeneity has been suggested to fuel resistance to many targeted therapies. The cytosine deaminase APOBEC proteins have been implicated in the mutational signatures of more than 70% of human cancers. However, the mechanism underlying how cancer cells hijack the APOBEC mediated mutagenesis machinery to promote tumor heterogeneity, and thereby foster therapy resistance remains unclear. We identify SYNCRIP as an endogenous molecular brake which suppresses APOBEC-driven mutagenesis in prostate cancer (PCa). Overactivated APOBEC3B, in SYNCRIP-deficient PCa cells, is a key mutator, representing the molecular source of driver mutations in some frequently mutated genes in PCa, including FOXA1, EP300. Functional screening identifies eight crucial drivers for androgen receptor (AR)-targeted therapy resistance in PCa that are mutated by APOBEC3B: BRD7, CBX8, EP300, FOXA1, HDAC5, HSF4, STAT3, and AR. These results uncover a cell-intrinsic mechanism that unleashes APOBEC-driven mutagenesis, which plays a significant role in conferring AR-targeted therapy resistance in PCa.

Dissecting Intra-tumor Heterogeneity in the Glioblastoma Microenvironment Using Fluorescence-Guided Multiple Sampling

The treatment of the most aggressive primary brain tumor in adults, glioblastoma (GBM), is challenging due to its heterogeneous nature, invasive potential, and poor response to chemo- and radiotherapy. As a result, GBM inevitably recurs and only a few patients survive 5 years post-diagnosis. GBM is characterized by extensive phenotypic and genetic heterogeneity, creating a diversified genetic landscape and a network of biological interactions between subclones, ultimately promoting tumor growth and therapeutic resistance. This includes spatial and temporal changes in the tumor microenvironment, which influence cellular and molecular programs in GBM and therapeutic responses. However, dissecting phenotypic and genetic heterogeneity at spatial and temporal levels is extremely challenging, and the dynamics of the GBM microenvironment cannot be captured by analysis of a single tumor sample. In this review, we discuss the current research on GBM heterogeneity, in particular, the utility and potential applications of fluorescence-guided multiple sampling to dissect phenotypic and genetic intra-tumor heterogeneity in the GBM microenvironment, identify tumor and non-tumor cell interactions and novel therapeutic targets in areas that are key for tumor growth and recurrence, and improve the molecular classification of GBM.

Glioblastoma heterogeneity at single cell resolution

Glioblastoma (GBM) is one of the deadliest types of cancer and highly refractory to chemoradiation and immunotherapy. One of the main reasons for this resistance to therapy lies within the heterogeneity of the tumor and its associated microenvironment. The vast diversity of cell states, composition of cells, and phenotypical characteristics makes it difficult to accurately classify GBM into distinct subtypes and find effective therapies. The advancement of sequencing technologies in recent years has further corroborated the heterogeneity of GBM at the single cell level. Recent studies have only begun to elucidate the different cell states present in GBM and how they correlate with sensitivity to therapy. Furthermore, it has become clear that GBM heterogeneity not only depends on intrinsic factors but also strongly differs between new and recurrent GBM, and treatment naïve and experienced patients. Understanding and connecting the complex cellular network that underlies GBM heterogeneity will be indispensable in finding new ways to tackle this deadly disease. Here, we present an overview of the multiple layers of GBM heterogeneity and discuss novel findings in the age of single cell technologies.

Applying single cell multi-omic analyses to understand treatment resistance in pediatric high grade glioma

Despite improvements in cancer patient outcomes seen in the past decade, tumor resistance to therapy remains a major impediment to achieving durable clinical responses. Intratumoral heterogeneity related to genetic, epigenetic, transcriptomic, proteomic, and metabolic differences between individual cancer cells has emerged as a driver of therapeutic resistance. This cell to cell heterogeneity can be assessed using single cell profiling technologies that enable the identification of tumor cell clones that exhibit similar defining features like specific mutations or patterns of DNA methylation. Single cell profiling of tumors before and after treatment can generate new insights into the cancer cell characteristics that confer therapeutic resistance by identifying intrinsically resistant sub-populations that survive treatment and by describing new cellular features that emerge post-treatment due to tumor cell evolution. Integrative, single cell analytical approaches have already proven advantageous in studies characterizing treatment-resistant clones in cancers where pre- and post-treatment patient samples are readily available, such as leukemia. In contrast, little is known about other cancer subtypes like pediatric high grade glioma, a class of heterogeneous, malignant brain tumors in children that rapidly develop resistance to multiple therapeutic modalities, including chemotherapy, immunotherapy, and radiation. Leveraging single cell multi-omic technologies to analyze naïve and therapy-resistant glioma may lead to the discovery of novel strategies to overcome treatment resistance in brain tumors with dismal clinical outcomes. In this review, we explore the potential for single cell multi-omic analyses to reveal mechanisms of glioma resistance to therapy and discuss opportunities to apply these approaches to improve long-term therapeutic response in pediatric high grade glioma and other brain tumors with limited treatment options.

The Heterogeneous Cellular States of Glioblastoma Stem Cells Revealed by Single Cell Analysis

Glioblastoma stem cells (GSCs) contributed to the progression, treatment resistance, and relapse of glioblastoma (GBM). However, current researches on GSCs were performed usually outside the human tumor microenvironment, ignoring the importance of the cellular states of primary GSCs. In this study, we leveraged single-cell transcriptome sequencing data of 6 independent GBM cohorts from public databases, and combined lineage and stemness features to identify primary GSCs. We dissected the cell states of GSCs and correlated them with the clinical outcomes of patients. As a result, we constructed a cellular hierarchy where GSCs resided at the center. In addition, we identified and characterized 2 different and recurrent GSCs subpopulations: proliferative GSCs (pGSCs) and quiescent GSCs (qGSCs). The pGSCs showed high cell cycle activity, indicating rapid cell division, while qGSCs showed a quiescent state. Then we traced the processes of tumor development by pseudo-time analysis and tumor phylogeny, and found that GSCs accumulated throughout the whole tumor development period. During the process, pGSCs mainly contributed to the early stage and qGSCs were enriched in the later stage. Finally, we constructed an 8-gene prognostic signature reflecting pGSCs activity and found that patients whose tumors were enriched for the pGSC signature had poor clinical outcomes. Our study highlights the primary GSCs heterogeneity and its correlation to tumor development and clinical outcomes, providing the potential targets for GBM treatment.

Single-cell transcriptomic analysis deciphers key transitional signatures associated with oncogenic evolution in human intramucosal oesophageal squamous cell carcinoma

BACKGROUND AND AIMS

The early diagnosis and intervention of oesophageal squamous cell carcinoma (ESCC) are particularly important because of the lack of effective therapies and poor prognosis. Comprehensive research on early ESCC at the single-cell level is rare due to the need for fresh and high-quality specimens obtained from ESD. This study aims to systematically describe the cellular atlas of human intramucosal ESCC.

METHODS

Five paired samples of intramucosal ESCC, para-ESCC oesophageal tissues from endoscopically resected specimens and peripheral blood mononuclear cells were adopted for scRNA-seq analysis. Computational pipeline scMetabolism was applied to quantify the metabolic diversity of single cells.

RESULTS

A total of 164 715 cells were profiled. Epithelial cells exhibited high intra-tumoural heterogeneity and two evolutionary trajectories during ESCC tumorigenesis initiated from proliferative cells, and then through an intermediate state, to two different terminal states of normally differentiated epithelial cells or malignant cells, respectively. The abundance of CD8+ TEX s, Tregs and PD1+ CD4+ T cells suggested an exhausted and suppressive immune microenvironment. Several genes in immune cells, such as CXCL13, CXCR5 and PADI4, were identified as new biomarkers for poor prognosis. A new subcluster of malignant cells associated with metastasis and angiogenesis that appeared at an early stage compared with progressive ESCC was also identified in this study. Intercellular interaction analysis based on ligand-receptor pairs revealed the subcluster of malignant cells interacting with CAFs via the MDK-NCL pathway, which was verified by cell proliferation assay and IHC. This indicates that the interaction may be an important hallmark in the early change of tumour microenvironment and serves as a sign of CAF activation to stimulate downstream pathways for facilitating tumour invasion.

CONCLUSION

This study demonstrates the changes of cell subsets and transcriptional levels in human intramucosal ESCC, which may provide unique insights into the development of novel biomarkers and potential intervention strategies.

Single Cell Analysis of Human Thyroid Reveals the Transcriptional Signatures of Aging

The thyroid gland plays a critical role in the maintenance of whole-body metabolism. However, aging frequently impairs homeostatic maintenance by thyroid hormones due to increased prevalence of subclinical hypothyroidism associated with mitochondrial dysfunction, inflammation, and fibrosis. To understand the specific aging-related changes of endocrine function in thyroid epithelial cells, we performed single-cell RNA sequencing (RNA-seq) of 54 726 cells derived from pathologically normal thyroid tissues from 7 patients who underwent thyroidectomy. Thyroid endocrine epithelial cells were clustered into 5 distinct subpopulations, and a subset of cells was found to be particularly vulnerable with aging, showing functional deterioration associated with the expression of metallothionein (MT) and major histocompatibility complex class II genes. We further validated that increased expression of MT family genes are highly correlated with thyroid gland aging in bulk RNAseq datasets. This study provides evidence that aging induces specific transcriptomic changes across multiple cell populations in the human thyroid gland.

Complex Analysis of Single-Cell RNA Sequencing Data

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.

A single-cell atlas of glioblastoma evolution under therapy reveals cell-intrinsic and cell-extrinsic therapeutic targets

Recent longitudinal studies of glioblastoma (GBM) have demonstrated a lack of apparent selection pressure for specific DNA mutations in recurrent disease. Single-cell lineage tracing has shown that GBM cells possess a high degree of plasticity. Together this suggests that phenotype switching, as opposed to genetic evolution, may be the escape mechanism that explains the failure of precision therapies to date. We profiled 86 primary-recurrent patient-matched paired GBM specimens with single-nucleus RNA, single-cell open-chromatin, DNA and spatial transcriptomic/proteomic assays. We found that recurrent GBMs are characterized by a shift to a mesenchymal phenotype. We show that the mesenchymal state is mediated by activator protein 1. Increased T-cell abundance at recurrence was prognostic and correlated with hypermutation status. We identified tumor-supportive networks of paracrine and autocrine signals between GBM cells, nonmalignant neuroglia and immune cells. We present cell-intrinsic and cell-extrinsic targets and a single-cell multiomics atlas of GBM under therapy.

Single-cell transcriptome reveals cellular hierarchies and guides p-EMT-targeted trial in skull base chordoma

Skull base chordoma (SBC) is a bone cancer with a high recurrence rate, high radioresistance rate, and poorly understood mechanism. Here, we profiled the transcriptomes of 90,691 single cells, revealed the SBC cellular hierarchies, and explored novel treatment targets. We identified a cluster of stem-like SBC cells that tended to be distributed in the inferior part of the tumor. Combining radiated UM-Chor1 RNA-seq data and in vitro validation, we further found that this stem-like cell cluster is marked by cathepsin L (CTSL), a gene involved in the packaging of telomere ends, and may be responsible for radioresistance. Moreover, signatures related to partial epithelial-mesenchymal transition (p-EMT) were found to be significant in malignant cells and were related to the invasion and poor prognosis of SBC. Furthermore, YL-13027, a p-EMT inhibitor that acts through the TGF-β signaling pathway, demonstrated remarkable potency in inhibiting the invasiveness of SBC in preclinical models and was subsequently applied in a phase I clinical trial that enrolled three SBC patients. Encouragingly, YL-13027 attenuated the growth of SBC and achieved stable disease with no serious adverse events, underscoring the clinical potential for the precision treatment of SBC with this therapy. In summary, we conducted the first single-cell RNA sequencing of SBC and identified several targets that could be translated to the treatment of SBC.

Relationship Between CNVs and Immune Cells Infiltration in Gastric Tumor Microenvironment

Gastric cancer (GC) is a highly fatal and common malignancy of the digestive system. Recent therapeutic advancements have significantly improved the clinical outcomes in GC, but due to the unavailability of suitable molecular targets, a large number of patients do not respond to the immune checkpoint inhibitors (ICI) therapy. To identify and validate potential therapeutic and prognostic targets of gastric cancer, we used the “inferCNV” R package for analyzing single-cell sequencing data (GSE112302) of GC and normal epithelial cells. First, by using LASSO, we screened genes that were highly correlated with copy number variations (CNVs). Therefrom, five gene signature (CPVL, DDC, GRTP1, ONECUT2, and PRSS21) was selected by cross-validating the prognosis and risk management with the GC RNA-seq data obtained from GEO and TCGA. Moreover, the correlation analyses between CNVs of these genes and immune cell infiltration in gastric cancer identified CPVL as a potential prognostic marker. Finally, CPVL showed high expression in gastric cancer samples and cell lines, then siRNA-mediated silencing of CPVL expression in gastric cancer cells showed significant proliferation arrest in MGC803 cells. Here, we conclude that CNVs are key regulators of the immune cells infiltration in gastric TME as well as cancer development, and CPVL could potentially be used as a prognostic and therapeutic marker in gastric cancer.

MYBL2 accelerates epithelial-mesenchymal transition and hepatoblastoma metastasis via the Smad/SNAI1 pathway

Hepatoblastoma (HB) accounts for the majority of hepatic malignancies in children. Although the prognosis of patients with HB has improved in past decades, metastasis is an indicator of poor overall survival. Herein, we applied single-cell RNA sequencing to explore the transcriptomic profiling of 25,264 metastatic cells isolated from the lungs of two patients with HB. The transcriptomes uncovered the heterogeneity of malignant cells after metastatic lung colonization, and these cells had varied expression signatures associated with the cell cycle, epithelial-mesenchymal plasticity, and hepatic differentiation. Single-cell regulatory network inference and clustering (SCENIC) was utilized to identify the co-expressed transcriptional factors which regulated and represented the different cell states. We further screened the key factor by bioinformatics analysis and found that MYBL2 upregulation was significantly associated with metastasis and poor prognosis. The relationship between ectopic MYBL2 and metastasis was subsequently proved by immunohistochemistry (IHC) of HB tissues, and the functions of MYBL2 in promoting proliferation, migration, and epithelial-to-mesenchymal transition (EMT) were verified by in vitro and in vivo assays. Importantly, the levels of Smad2/3 phosphorylation and SNAI1 expression were increased in MYBL2-transfected cells. Consequently, these results indicated that the MYBL2-controlled Smad/SNAI1 pathway induced EMT and promoted HB tumorigenesis and metastasis.

Single-Cell Transcriptomics of Cultured Amniotic Fluid Cells Reveals Complex Gene Expression Alterations in Human Fetuses With Trisomy 18

Trisomy 18, commonly known as Edwards syndrome, is the second most common autosomal trisomy among live born neonates. Multiple tissues including cardiac, abdominal, and nervous systems are affected by an extra chromosome 18. To delineate the complexity of anomalies of trisomy 18, we analyzed cultured amniotic fluid cells from two euploid and three trisomy 18 samples using single-cell transcriptomics. We identified 6 cell groups, which function in development of major tissues such as kidney, vasculature and smooth muscle, and display significant alterations in gene expression as detected by single-cell RNA-sequencing. Moreover, we demonstrated significant gene expression changes in previously proposed trisomy 18 critical regions, and identified three new regions such as 18p11.32, 18q11 and 18q21.32, which are likely associated with trisomy 18 phenotypes. Our results indicate complexity of trisomy 18 at the gene expression level and reveal genetic reasoning of diverse phenotypes in trisomy 18 patients.

Advancement of Single-Cell Sequencing in Medulloblastoma

Single-cell sequencing is a promising attempt to investigate the genomic, transcriptomic, and multiomic level of individual cell in the larger population of cells. The outward evolution of the technique from a manual method to the automation of single-cell sequencing is cogent. Lately, single-cell sequencing is widely used in various fields of science and has applications in neurobiology, immunity, cancer, microbiology, reproduction, and digestion. This chapter introduces the reader to the details of single-cell sequencing, currently used in several small-scale and commercial platforms. The advancement of single-cell sequencing in brain cancer sheds light on questions unanswered so far in the field of oncology.

Unravelling glioblastoma heterogeneity by means of single-cell RNA sequencing

Glioblastoma (GBM) is the most invasive and deadliest brain cancer in adults. Its inherent heterogeneity has been designated as the main cause of treatment failure. Thus, a deeper understanding of the diversity that shapes GBM pathobiology is of utmost importance. Single-cell RNA sequencing (scRNA-seq) technologies have begun to uncover the hidden composition of complex tumor ecosystems. Herein, a semi-systematic search of reference literature databases provided all existing publications using scRNA-seq for the investigation of human GBM. We compared and discussed findings from these works to build a more robust and unified knowledge base. All aspects ranging from inter-patient heterogeneity to intra-tumoral organization, cancer stem cell diversity, clonal mosaicism, and the tumor microenvironment (TME) are comprehensively covered in this report. Tumor composition not only differs across patients but also involves a great extent of heterogeneity within itself. Spatial and cellular heterogeneity can reveal tumor evolution dynamics. In addition, the discovery of distinct cell phenotypes might lead to the development of targeted treatment approaches. In conclusion, scRNA-seq expands our knowledge of GBM heterogeneity and helps to unravel putative therapeutic targets.

ATRX regulates glial identity and the tumor microenvironment in IDH-mutant glioma

BACKGROUND

Recent single-cell transcriptomic studies report that IDH-mutant gliomas share a common hierarchy of cellular phenotypes, independent of genetic subtype. However, the genetic differences between IDH-mutant glioma subtypes are prognostic, predictive of response to chemotherapy, and correlate with distinct tumor microenvironments.

RESULTS

To reconcile these findings, we profile 22 human IDH-mutant gliomas using scATAC-seq and scRNA-seq. We determine the cell-type-specific differences in transcription factor expression and associated regulatory grammars between IDH-mutant glioma subtypes. We find that while IDH-mutant gliomas do share a common distribution of cell types, there are significant differences in the expression and targeting of transcription factors that regulate glial identity and cytokine elaboration. We knock out the chromatin remodeler ATRX, which suffers loss-of-function alterations in most IDH-mutant astrocytomas, in an IDH-mutant immunocompetent intracranial murine model. We find that both human ATRX-mutant gliomas and murine ATRX-knockout gliomas are more heavily infiltrated by immunosuppressive monocytic-lineage cells derived from circulation than ATRX-intact gliomas, in an IDH-mutant background. ATRX knockout in murine glioma recapitulates gene expression and open chromatin signatures that are specific to human ATRX-mutant astrocytomas, including drivers of astrocytic lineage and immune-cell chemotaxis. Through single-cell cleavage under targets and tagmentation assays and meta-analysis of public data, we show that ATRX loss leads to a global depletion in CCCTC-binding factor association with DNA, gene dysregulation along associated chromatin loops, and protection from therapy-induced senescence.

CONCLUSIONS

These studies explain how IDH-mutant gliomas from different subtypes maintain distinct phenotypes and tumor microenvironments despite a common lineage hierarchy.

Single-cell transcriptomic analysis of the tumor ecosystems underlying initiation and progression of papillary thyroid carcinoma

The tumor ecosystem of papillary thyroid carcinoma (PTC) is poorly characterized. Using single-cell RNA sequencing, we profile transcriptomes of 158,577 cells from 11 patients’ paratumors, localized/advanced tumors, initially-treated/recurrent lymph nodes and radioactive iodine (RAI)-refractory distant metastases, covering comprehensive clinical courses of PTC. Our data identifies a “cancer-primed” premalignant thyrocyte population with normal morphology but altered transcriptomes. Along the developmental trajectory, we also discover three phenotypes of malignant thyrocytes (follicular-like, partial-epithelial-mesenchymal-transition-like, dedifferentiation-like), whose composition shapes bulk molecular subtypes, tumor characteristics and RAI responses. Furthermore, we uncover a distinct BRAF-like-B subtype with predominant dedifferentiation-like thyrocytes, enriched cancer-associated fibroblasts, worse prognosis and promising prospect of immunotherapy. Moreover, potential vascular-immune crosstalk in PTC provides theoretical basis for combined anti-angiogenic and immunotherapy. Together, our findings provide insight into the PTC ecosystem that suggests potential prognostic and therapeutic implications.

The characterisation of the papillary thyroid carcinoma (PTC) tumour microenvironment remains crucial. Here, the authors perform single-cell RNA sequencing in 11 patients and identify potential opportunities for the use of immunotherapy and its combination with anti-angiogenic therapy in PTC.

A deep convolutional neural network for segmentation of whole-slide pathology images identifies novel tumour cell-perivascular niche interactions that are associated with poor survival in glioblastoma

BACKGROUND

Glioblastoma is the most aggressive type of brain cancer with high-levels of intra- and inter-tumour heterogeneity that contribute to its rapid growth and invasion within the brain. However, a spatial characterisation of gene signatures and the cell types expressing these in different tumour locations is still lacking.

METHODS

We have used a deep convolutional neural network (DCNN) as a semantic segmentation model to segment seven different tumour regions including leading edge (LE), infiltrating tumour (IT), cellular tumour (CT), cellular tumour microvascular proliferation (CTmvp), cellular tumour pseudopalisading region around necrosis (CTpan), cellular tumour perinecrotic zones (CTpnz) and cellular tumour necrosis (CTne) in digitised glioblastoma histopathological slides from The Cancer Genome Atlas (TCGA). Correlation analysis between segmentation results from tumour images together with matched RNA expression data was performed to identify genetic signatures that are specific to different tumour regions.

RESULTS

We found that spatially resolved gene signatures were strongly correlated with survival in patients with defined genetic mutations. Further in silico cell ontology analysis along with single-cell RNA sequencing data from resected glioblastoma tissue samples showed that these tumour regions had different gene signatures, whose expression was driven by different cell types in the regional tumour microenvironment. Our results further pointed to a key role for interactions between microglia/pericytes/monocytes and tumour cells that occur in the IT and CTmvp regions, which may contribute to poor patient survival.

CONCLUSIONS

This work identified key histopathological features that correlate with patient survival and detected spatially associated genetic signatures that contribute to tumour-stroma interactions and which should be investigated as new targets in glioblastoma. The source codes and datasets used are available in GitHub: https://github.com/amin20/GBM_WSSM .

Dissecting the single-cell transcriptome network underlying esophagus non-malignant tissues and esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is among the most prevalent causes of cancer-related death in adults. Tumor microenvironment (TME) has been associated with therapeutic failure and lethal outcomes for patients. However, published reports on the heterogeneity and TME in ESCC are scanty.

METHODS

Five tumor samples and five corresponding non-malignant samples were subjected to scRNA-seq analysis. Bulk RNA sequencing data were retrieved in publicly available databases.

FINDINGS

From the scRNA-seq data, a total of 128,688 cells were enrolled for subsequent analyses. Gene expression and CNV status exhibited high heterogeneity of tumor cells. We further identified a list of tumor-specific genes and four malignant signatures, which are potential new markers for ESCC. Metabolic analysis revealed that energy supply-related pathways are pivotal in cancer metabolic reprogramming. Moreover, significant differences were found in stromal and immune cells between the esophagus normal and tumor tissues, which promoted carcinogenesis at both cellular and molecular levels in ESCC. Immune checkpoints, regarded as potential targets for immunotherapy in ESCC were significantly highly expressed in ESCC, including LAG3 and HAVCR2. Eventually, we constructed a cell-to-cell communication atlas based on cancer cells and immune cells and performed the flow cytometry, qRT-PCR, immunofluorescence, and immunohistochemistry analyses to validate the results.

INTERPRETATION

This study demonstrates a widespread reprogramming across multiple cellular elements within the TME in ESCC, particularly in transcriptional states, cellular functions, and cell-to-cell interactions. The findings offer an insight into the exploration of TME and heterogeneity in the ESCC and provide new therapeutic targets for its clinical management in the future.

FUNDING

The work was supported by the Shanghai Pujiang Program (2020PJD009) and Research Development Fund of Zhongshan Hospital, Fudan University (2019ZSFZ002 and 2019ZSFZ19).

A Drug Screening Pipeline Using 2D and 3D Patient-Derived In Vitro Models for Pre-Clinical Analysis of Therapy Response in Glioblastoma

Glioblastoma is one of the most common and lethal types of primary brain tumor. Despite aggressive treatment with chemotherapy and radiotherapy, tumor recurrence within 6-9 months is common. To overcome this, more effective therapies targeting cancer cell stemness, invasion, metabolism, cell death resistance and the interactions of tumor cells with their surrounding microenvironment are required. In this study, we performed a systematic review of the molecular mechanisms that drive glioblastoma progression, which led to the identification of 65 drugs/inhibitors that we screened for their efficacy to kill patient-derived glioma stem cells in two dimensional (2D) cultures and patient-derived three dimensional (3D) glioblastoma explant organoids (GBOs). From the screening, we found a group of drugs that presented different selectivity on different patient-derived in vitro models. Moreover, we found that Costunolide, a TERT inhibitor, was effective in reducing the cell viability in vitro of both primary tumor models as well as tumor models pre-treated with chemotherapy and radiotherapy. These results present a novel workflow for screening a relatively large groups of drugs, whose results could lead to the identification of more personalized and effective treatment for recurrent glioblastoma.

Applying high-dimensional single-cell technologies to the analysis of cancer immunotherapy

Advances in molecular biology, microfluidics and bioinformatics have empowered the study of thousands or even millions of individual cells from malignant tumours at the single-cell level of resolution. This high-dimensional, multi-faceted characterization of the genomic, transcriptomic, epigenomic and proteomic features of the tumour and/or the associated immune and stromal cells enables the dissection of tumour heterogeneity, the complex interactions between tumour cells and their microenvironment, and the details of the evolutionary trajectory of each tumour. Single-cell transcriptomics, the ability to track individual T cell clones through paired sequencing of the T cell receptor genes and high-dimensional single-cell spatial analysis are all areas of particular relevance to immuno-oncology. Multidimensional biomarker signatures will increasingly be crucial to guiding clinical decision-making in each patient with cancer. High-dimensional single-cell technologies are likely to provide the resolution and richness of data required to generate such clinically relevant signatures in immuno-oncology. In this Perspective, we describe advances made using transformative single-cell analysis technologies, especially in relation to clinical response and resistance to immunotherapy, and discuss the growing utility of single-cell approaches for answering important research questions.

Decoding Cancer Biology One Cell at a Time

Human tumors are composed of diverse malignant and nonmalignant cells, generating a complex ecosystem that governs tumor biology and response to treatments. Recent technological advances have enabled the characterization of tumors at single-cell resolution, providing a compelling strategy to dissect their intricate biology. Here we describe recent developments in single-cell expression profiling and the studies applying them in clinical settings. We highlight some of the powerful insights gleaned from these studies for tumor classification, stem cell programs, tumor microenvironment, metastasis, and response to targeted and immune therapies. SIGNIFICANCE: Intratumor heterogeneity (ITH) has been a major barrier to our understanding of cancer. Single-cell genomics is leading a revolution in our ability to systematically dissect ITH. In this review, we focus on single-cell expression profiling and lessons learned in key aspects of human tumor biology.

3D-printed microplate inserts for long term high-resolution imaging of live brain organoids

BACKGROUND

Organoids are a reliable model used in the study of human brain development and under pathological conditions. However, current methods for brain organoid culture generate tissues that range from 0.5 to 2 mm of size, which need to be constantly agitated to allow proper oxygenation. The culture conditions are, therefore, not suitable for whole-brain organoid live imaging, required to study developmental processes and disease progression within physiologically relevant time frames (i.e. days, weeks, months).

RESULTS

Here we designed 3D-printed microplate inserts adaptable to standard 24 multi-well plates, which allow the growth of multiple organoids in pre-defined and fixed XYZ coordinates. This innovation facilitates high-resolution imaging of whole-cerebral organoids, allowing precise assessment of organoid growth and morphology, as well as cell tracking within the organoids, over long periods. We applied this technology to track neocortex development through neuronal progenitors in brain organoids, as well as the movement of patient-derived glioblastoma stem cells within healthy brain organoids.

CONCLUSIONS

This new bioengineering platform constitutes a significant advance that permits long term detailed analysis of whole-brain organoids using multimodal inverted fluorescence microscopy.

CACTUS: integrating clonal architecture with genomic clustering and transcriptome profiling of single tumor cells

Background

Drawing genotype-to-phenotype maps in tumors is of paramount importance for understanding tumor heterogeneity. Assignment of single cells to their tumor clones of origin can be approached by matching the genotypes of the clones to the mutations found in RNA sequencing of the cells. The confidence of the cell-to-clone mapping can be increased by accounting for additional measurements. Follicular lymphoma, a malignancy of mature B cells that continuously acquire mutations in parallel in the exome and in B cell receptor loci, presents a unique opportunity to join exome-derived mutations with B cell receptor sequences as independent sources of evidence for clonal evolution.

Methods

Here, we propose CACTUS, a probabilistic model that leverages the information from an independent genomic clustering of cells and exploits the scarce single cell RNA sequencing data to map single cells to given imperfect genotypes of tumor clones.

Results

We apply CACTUS to two follicular lymphoma patient samples, integrating three measurements: whole exome, single-cell RNA, and B cell receptor sequencing. CACTUS outperforms a predecessor model by confidently assigning cells and B cell receptor-based clusters to the tumor clones.

Conclusions

The integration of independent measurements increases model certainty and is the key to improving model performance in the challenging task of charting the genotype-to-phenotype maps in tumors. CACTUS opens the avenue to study the functional implications of tumor heterogeneity, and origins of resistance to targeted therapies. CACTUS is written in R and source code, along with all supporting files, are available on GitHub (https://github.com/LUMC/CACTUS).

Supplementary Information

The online version contains supplementary material available at (10.1186/s13073-021-00842-w).

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.

The evolution of alternative splicing in glioblastoma under therapy

Background

Alternative splicing is a rich source of tumor-specific neoantigen targets for immunotherapy. This holds promise for glioblastomas (GBMs), the most common primary tumors of the adult brain, which are resistant to standard-of-care therapy. Although most clinical trials enroll patients at recurrence, most preclinical studies have been done with specimens from primary disease. There are limited expression data from GBMs at recurrence and surprisingly little is known about the evolution of splicing patterns under therapy.

Result

We profile 37 primary-recurrent paired human GBM specimens via RNA sequencing. We describe the landscape of alternative splicing in GBM at recurrence and contrast that to primary and non-malignant brain-tissue specimens. By screening single-cell atlases, we identify cell-type-specific splicing patterns and novel splicing events in cell-surface proteins that are suitable targets for engineered T cell therapies. We identify recurrent-specific isoforms of mitogen-activated kinase pathway genes that enhance invasiveness and are preferentially expressed by stem-like cells.

Conclusion

These studies shed light on gene expression in recurrent GBM and identify novel targets for therapeutic development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02259-5.

Genetic Architectures and Cell-of-Origin in Glioblastoma

An aggressive primary brain cancer, glioblastoma (GBM) is the most common cancer of the central nervous system in adults. However, an inability to identify its cell-of-origin has been a fundamental issue hindering further understanding of the nature and pathogenesis of GBM, as well as the development of novel therapeutic targets. Researchers have hypothesized that GBM arises from an accumulation of somatic mutations in neural stem cells (NSCs) and glial precursor cells that confer selective growth advantages, resulting in uncontrolled proliferation. In this review, we outline genomic perspectives on IDH-wildtype and IDH-mutant GBMs pathogenesis and the cell-of-origin harboring GBM driver mutations proposed by various GBM animal models. Additionally, we discuss the distinct neurodevelopmental programs observed in either IDH-wildtype or IDH-mutant GBMs. Further research into the cellular origin and lineage hierarchy of GBM will help with understanding the evolution of GBMs and with developing effective targets for treating GBM cancer cells.

Macropinocytosis: mechanism and targeted therapy in cancers

Macropinocytosis is a form of endocytosis which provides an effective way for non-selective uptakes of extracellular proteins, liquids, and particles. The endocytic process is initiated by the activation of the growth factors signaling pathways. After activation of the biochemical signal, the cell starts internalizing extracellular solutes and nutrients into the irregular endocytic vesicles, known as macropinosomes that deliver them into the lysosomes for degradation. Macropinocytosis plays an important role in the nutritional supply of cancer cells. Due to the rapid expansion of cancer cells and the abnormal vascular microenvironment, cancer cells are usually deprived of oxygen and nutrients. Therefore, they must transform their metabolism to survive and grow in this harsh microenvironment. To satisfy their energy needs, cancer cells enhance the activity of macropinocytosis. Therefore, this metabolic adaptation that is used by cancer cells can be exploited to develop new targeted cancer therapies. In this review, we discuss the molecular mechanism that actuates the process of macropinocytosis in a variety of cancers, and the novel anti-cancer therapeutics in targeting macropinocytosis.

Single-cell transcriptomics in cancer: computational challenges and opportunities

Intratumor heterogeneity is a common characteristic across diverse cancer types and presents challenges to current standards of treatment. Advancements in high-throughput sequencing and imaging technologies provide opportunities to identify and characterize these aspects of heterogeneity. Notably, transcriptomic profiling at a single-cell resolution enables quantitative measurements of the molecular activity that underlies the phenotypic diversity of cells within a tumor. Such high-dimensional data require computational analysis to extract relevant biological insights about the cell types and states that drive cancer development, pathogenesis, and clinical outcomes. In this review, we highlight emerging themes in the computational analysis of single-cell transcriptomics data and their applications to cancer research. We focus on downstream analytical challenges relevant to cancer research, including how to computationally perform unified analysis across many patients and disease states, distinguish neoplastic from nonneoplastic cells, infer communication with the tumor microenvironment, and delineate tumoral and microenvironmental evolution with trajectory and RNA velocity analysis. We include discussions of challenges and opportunities for future computational methodological advancements necessary to realize the translational potential of single-cell transcriptomic profiling in cancer.

Single-cell transcriptomics reveals regulators underlying immune cell diversity and immune subtypes associated with prognosis in nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is an aggressive malignancy with extremely skewed ethnic and geographic distributions. Increasing evidence indicates that targeting the tumor microenvironment (TME) represents a promising therapeutic approach in NPC, highlighting an urgent need to deepen the understanding of the complex NPC TME. Here, we generated single-cell transcriptome profiles for 7581 malignant cells and 40,285 immune cells from fifteen primary NPC tumors and one normal sample. We revealed malignant signatures capturing intratumoral transcriptional heterogeneity and predicting aggressiveness of malignant cells. Diverse immune cell subtypes were identified, including novel subtypes such as CLEC9A⁺ dendritic cells (DCs). We further revealed transcriptional regulators underlying immune cell diversity, and cell–cell interaction analyses highlighted promising immunotherapeutic targets in NPC. Moreover, we established the immune subtype-specific signatures, and demonstrated that the signatures of macrophages, plasmacytoid dendritic cells (pDCs), CLEC9A⁺ DCs, natural killer (NK) cells, and plasma cells were significantly associated with improved survival outcomes in NPC. Taken together, our findings represent a unique resource providing in-depth insights into the cellular heterogeneity of NPC TME and highlight potential biomarkers for anticancer treatment and risk stratification, laying a new foundation for precision therapies in NPC.

Cell-level somatic mutation detection from single-cell RNA sequencing

MOTIVATION

Both single-cell RNA sequencing (scRNA-seq) and DNA sequencing (scDNA-seq) have been applied for cell-level genomic profiling. For mutation profiling, the latter seems more natural. However, the task is highly challenging due to the limited input materials from only two copies of DNA molecules, while whole-genome amplification generates biases and other technical noises. ScRNA-seq starts with a higher input amount, so generally has better data quality. There exists various methods for mutation detection from DNA sequencing, it is not clear whether these methods work for scRNA-seq data.

RESULTS

Mutation detection methods developed for either bulk-cell sequencing data or scDNA-seq data do not work well for the scRNA-seq data, as they produce substantial numbers of false positives. We develop a novel and robust statistical method-called SCmut-to identify specific cells that harbor mutations discovered in bulk-cell data. Statistically SCmut controls the false positives using the 2D local false discovery rate method. We apply SCmut to several scRNA-seq datasets. In scRNA-seq breast cancer datasets SCmut identifies a number of highly confident cell-level mutations that are recurrent in many cells and consistent in different samples. In a scRNA-seq glioblastoma dataset, we discover a recurrent cell-level mutation in the PDGFRA gene that is highly correlated with a well-known in-frame deletion in the gene. To conclude, this study contributes a novel method to discover cell-level mutation information from scRNA-seq that can facilitate investigation of cell-to-cell heterogeneity.

AVAILABILITY AND IMPLEMENTATION

The source codes and bioinformatics pipeline of SCmut are available at https://github.com/nghiavtr/SCmut.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Next Generation Sequencing and Machine Learning Technologies Are Painting the Epigenetic Portrait of Glioblastoma

Even with a rare occurrence of only 1.35% of cancer cases in the United States of America, brain tumors are considered as one of the most lethal malignancies. The most aggressive and invasive type of brain tumor, glioblastoma, accounts for 60–70% of all gliomas and presents with life expectancy of only 12–18 months. Despite trimodal treatment and advances in diagnostic and therapeutic methods, there are no significant changes in patient outcome. Our understanding of glioblastoma was significantly improved with the introduction of next generation sequencing technologies. This led to the identification of different genetic and molecular subtypes, which greatly improve glioblastoma diagnosis. Still, because of the poor life expectancy, novel diagnostic, and treatment methods are broadly explored. Epigenetic modifications like methylation and changes in histone acetylation are such examples. Recently, in addition to genetic and molecular characteristics, epigenetic profiling of glioblastomas is also used for sample classification. Further advancement of next generation sequencing technologies is expected to identify in detail the epigenetic signature of glioblastoma that can open up new therapeutic opportunities for glioblastoma patients. This should be complemented with the use of computational power i.e., machine and deep learning algorithms for objective diagnostics and design of individualized therapies. Using a combination of phenotypic, genotypic, and epigenetic parameters in glioblastoma diagnostics will bring us closer to precision medicine where therapies will be tailored to suit the genetic profile and epigenetic signature of the tumor, which will grant longer life expectancy and better quality of life. Still, a number of obstacles including potential bias, availability of data for minorities in heterogeneous populations, data protection, and validation and independent testing of the learning algorithms have to be overcome on the way.

The Glioma Stem Cell Model in the Era of Single-Cell Genomics

Glioma stem cells (GSCs) are thought to underlie glioma initiation, evolution, and resistance to existing therapies. Although functional evidence for GSCs is abundant, tumor heterogeneity and intrinsic limitations in GSC assays have represented barriers for the field. In this perspective, we revisit the GSC model in light of recent single-cell expression profiling studies. We highlight how classes of glioma differ in their cellular architecture and relate the observed cellular states to established GSC markers. We additionally propose a set of single-cell informed definitions as a framework for our understanding of the cellular architecture of gliomas and a potential therapeutic outlook.

DeepSurvNet: deep survival convolutional network for brain cancer survival rate classification based on histopathological images

Histopathological whole slide images of haematoxylin and eosin (H&E)-stained biopsies contain valuable information with relation to cancer disease and its clinical outcomes. Still, there are no highly accurate automated methods to correlate histolopathological images with brain cancer patients' survival, which can help in scheduling patients therapeutic treatment and allocate time for preclinical studies to guide personalized treatments. We now propose a new classifier, namely, DeepSurvNet powered by deep convolutional neural networks, to accurately classify in 4 classes brain cancer patients' survival rate based on histopathological images (class I, 0-6 months; class II, 6-12 months; class III, 12-24 months; and class IV, >24 months survival after diagnosis). After training and testing of DeepSurvNet model on a public brain cancer dataset, The Cancer Genome Atlas, we have generalized it using independent testing on unseen samples. Using DeepSurvNet, we obtained precisions of 0.99 and 0.8 in the testing phases on the mentioned datasets, respectively, which shows DeepSurvNet is a reliable classifier for brain cancer patients' survival rate classification based on histopathological images. Finally, analysis of the frequency of mutations revealed differences in terms of frequency and type of genes associated to each class, supporting the idea of a different genetic fingerprint associated to patient survival. We conclude that DeepSurvNet constitutes a new artificial intelligence tool to assess the survival rate in brain cancer. Graphical abstract A DCNN model was generated to accurately predict survival rates of brain cancer patients (classified in 4 different classes) accurately. After training the model using images from H&E stained tissue biopsies from The Cancer Genome Atlas database (TCGA, left), the model can predict for each patient, based on a histological image (top right), its survival class accurately (bottom right).

Cardelino: computational integration of somatic clonal substructure and single-cell transcriptomes

Bulk and single-cell DNA sequencing has enabled reconstructing clonal substructures of somatic tissues from frequency and cooccurrence patterns of somatic variants. However, approaches to characterize phenotypic variations between clones are not established. Here we present cardelino (https://github.com/single-cell-genetics/cardelino), a computational method for inferring the clonal tree configuration and the clone of origin of individual cells assayed using single-cell RNA-seq (scRNA-seq). Cardelino flexibly integrates information from imperfect clonal trees inferred based on bulk exome-seq data, and sparse variant alleles expressed in scRNA-seq data. We apply cardelino to a published cancer dataset and to newly generated matched scRNA-seq and exome-seq data from 32 human dermal fibroblast lines, identifying hundreds of differentially expressed genes between cells from different somatic clones. These genes are frequently enriched for cell cycle and proliferation pathways, indicating a role for cell division genes in somatic evolution in healthy skin.

Intratumoral Heterogeneity and Longitudinal Changes in Gene Expression Predict Differential Drug Sensitivity in Newly Diagnosed and Recurrent Glioblastoma

Background: Inevitable recurrence after radiochemotherapy is the major problem in the treatment of glioblastoma, the most prevalent type of adult brain malignancy. Glioblastomas are notorious for a high degree of intratumor heterogeneity manifest through a diversity of cell types and molecular patterns. The current paradigm of understanding glioblastoma recurrence is that cytotoxic therapy fails to target effectively glioma stem cells. Recent advances indicate that therapy-driven molecular evolution is a fundamental trait associated with glioblastoma recurrence. There is a growing body of evidence indicating that intratumor heterogeneity, longitudinal changes in molecular biomarkers and specific impacts of glioma stem cells need to be taken into consideration in order to increase the accuracy of molecular diagnostics still relying on readouts obtained from a single tumor specimen. Methods: This study integrates a multisampling strategy, longitudinal approach and complementary transcriptomic investigations in order to identify transcriptomic traits of recurrent glioblastoma in whole-tissue specimens of glioblastoma or glioblastoma stem cells. In this study, 128 tissue samples of 44 tumors including 23 first diagnosed, 19 recurrent and 2 secondary recurrent glioblastomas were analyzed along with 27 primary cultures of glioblastoma stem cells by RNA sequencing. A novel algorithm was used to quantify longitudinal changes in pathway activities and model efficacy of anti-cancer drugs based on gene expression data. Results: Our study reveals that intratumor heterogeneity of gene expression patterns is a fundamental characteristic of not only newly diagnosed but also recurrent glioblastomas. Evidence is provided that glioblastoma stem cells recapitulate intratumor heterogeneity, longitudinal transcriptomic changes and drug sensitivity patterns associated with the state of recurrence. Conclusions: Our results provide a transcriptional rationale for the lack of significant therapeutic benefit from temozolomide in patients with recurrent glioblastoma. Our findings imply that the spectrum of potentially effective drugs is likely to differ between newly diagnosed and recurrent glioblastomas and underscore the merits of glioblastoma stem cells as prognostic models for identifying alternative drugs and predicting drug response in recurrent glioblastoma. With the majority of recurrent glioblastomas being inoperable, glioblastoma stem cell models provide the means of compensating for the limited availability of recurrent glioblastoma specimens.

Technological advances and computational approaches for alternative splicing analysis in single cells

Alternative splicing of RNAs generates isoform diversity, resulting in different proteins that are necessary for maintaining cellular function and identity. The discovery of alternative splicing has been revolutionized by next-generation transcriptomic sequencing mainly using bulk RNA-sequencing, which has unravelled RNA splicing and mis-splicing of normal cells under steady-state and stress conditions. Single-cell RNA-sequencing studies have focused on gene-level expression analysis and revealed gene expression signatures distinguishable between different cellular types. Single-cell alternative splicing is an emerging area of research with the promise to reveal transcriptomic dynamics invisible to bulk- and gene-level analysis. In this review, we will discuss the technological advances for single-cell alternative splicing analysis, computational strategies for isoform detection and quantitation in single cells, and current applications of single-cell alternative splicing analysis and its potential future contributions to personalized medicine.

Pathway Commons 2019 Update: integration, analysis and exploration of pathway data

Pathway Commons (https://www.pathwaycommons.org) is an integrated resource of publicly available information about biological pathways including biochemical reactions, assembly of biomolecular complexes, transport and catalysis events and physical interactions involving proteins, DNA, RNA, and small molecules (e.g. metabolites and drug compounds). Data is collected from multiple providers in standard formats, including the Biological Pathway Exchange (BioPAX) language and the Proteomics Standards Initiative Molecular Interactions format, and then integrated. Pathway Commons provides biologists with (i) tools to search this comprehensive resource, (ii) a download site offering integrated bulk sets of pathway data (e.g. tables of interactions and gene sets), (iii) reusable software libraries for working with pathway information in several programming languages (Java, R, Python and Javascript) and (iv) a web service for programmatically querying the entire dataset. Visualization of pathways is supported using the Systems Biological Graphical Notation (SBGN). Pathway Commons currently contains data from 22 databases with 4794 detailed human biochemical processes (i.e. pathways) and ∼2.3 million interactions. To enhance the usability of this large resource for end-users, we develop and maintain interactive web applications and training materials that enable pathway exploration and advanced analysis.

Genetic Hallmarks and Heterogeneity of Glioblastoma in the Single-Cell Omics Era

Glioblastoma multiforme is the most common and aggressive malignant primary brain tumor. As implied by its name, the disease displays impressive intrinsic heterogeneity. Among other complications, inter- and intratumoral diversity hamper glioblastoma research and therapy, typically leaving patients with little hope for long-term survival. Extensive genetic analyses, including omics, characterize several recurrent mutations. However, confounding factors mask crucial aspects of the pathology to conventional bulk approaches. In recent years, single-cell omics have made their first appearance in cancer research, and the methodology is about to reach its full potential for glioblastoma too. Here, recent glioblastoma single-cell omics investigations are reviewed, and most promising routes toward less grim prognoses and more efficient therapeutics are discussed.

Capture at the single cell level of metabolic modules distinguishing aggressive and indolent glioblastoma cells

Glioblastoma cell ability to adapt their functioning to microenvironment changes is a source of the extensive intra-tumor heterogeneity characteristic of this devastating malignant brain tumor. A systemic view of the metabolic pathways underlying glioblastoma cell functioning states is lacking. We analyzed public single cell RNA-sequencing data from glioblastoma surgical resections, which offer the closest available view of tumor cell heterogeneity as encountered at the time of patients’ diagnosis. Unsupervised analyses revealed that information dispersed throughout the cell transcript repertoires encoded the identity of each tumor and masked information related to cell functioning states. Data reduction based on an experimentally-defined signature of transcription factors overcame this hurdle. It allowed cell grouping according to their tumorigenic potential, regardless of their tumor of origin. The approach relevance was validated using independent datasets of glioblastoma cell and tissue transcriptomes, patient-derived cell lines and orthotopic xenografts. Overexpression of genes coding for amino acid and lipid metabolism enzymes involved in anti-oxidative, energetic and cell membrane processes characterized cells with high tumorigenic potential. Modeling of their expression network highlighted the very long chain polyunsaturated fatty acid synthesis pathway at the core of the network. Expression of its most downstream enzymatic component, ELOVL2, was associated with worsened patient survival, and required for cell tumorigenic properties in vivo. Our results demonstrate the power of signature-driven analyses of single cell transcriptomes to obtain an integrated view of metabolic pathways at play within the heterogeneous cell landscape of patient tumors.

Glioma through the looking GLASS: molecular evolution of diffuse gliomas and the Glioma Longitudinal Analysis Consortium

Adult diffuse gliomas are a diverse group of brain neoplasms that inflict a high emotional toll on patients and their families. The Cancer Genome Atlas and similar projects have provided a comprehensive understanding of the somatic alterations and molecular subtypes of glioma at diagnosis. However, gliomas undergo significant cellular and molecular evolution during disease progression. We review the current knowledge on the genomic and epigenetic abnormalities in primary tumors and after disease recurrence, highlight the gaps in the literature, and elaborate on the need for a new multi-institutional effort to bridge these knowledge gaps and how the Glioma Longitudinal Analysis Consortium (GLASS) aims to systemically catalog the longitudinal changes in gliomas. The GLASS initiative will provide essential insights into the evolution of glioma toward a lethal phenotype, with the potential to reveal targetable vulnerabilities and, ultimately, improved outcomes for a patient population in need.

CONICS integrates scRNA-seq with DNA sequencing to map gene expression to tumor sub-clones

Motivation

Single-cell RNA-sequencing (scRNA-seq) has enabled studies of tissue composition at unprecedented resolution. However, the application of scRNA-seq to clinical cancer samples has been limited, partly due to a lack of scRNA-seq algorithms that integrate genomic mutation data.

Results

To address this, we present.

CONICS

COpy-Number analysis In single-Cell RNA-Sequencing. CONICS is a software tool for mapping gene expression from scRNA-seq to tumor clones and phylogenies, with routines enabling: the quantitation of copy-number alterations in scRNA-seq, robust separation of neoplastic cells from tumor-infiltrating stroma, inter-clone differential-expression analysis and intra-clone co-expression analysis.

Availability and implementation

CONICS is written in Python and R, and is available from https://github.com/diazlab/CONICS.

Supplementary information

Supplementary data are available at Bioinformatics online.

Single-Cell Applications of Next-Generation Sequencing

The single cell is considered the basic unit of biology, and the pursuit of understanding how heterogeneous populations of cells can functionally coexist in tissues, organisms, microbial ecosystems, and even cancer, makes them the subject of intense study. Next-generation sequencing (NGS) of RNA and DNA has opened a new frontier of (single)-cell biology. Hundreds to millions of cells now can be assayed in parallel, providing the molecular profile of each cell in its milieu inexpensively and in a manner that can be analyzed mathematically. The goal of this article is to provide a high-level overview of single-cell sequencing for the nonexpert and show how its applications are influencing both basic and applied clinical studies in embryology, developmental genetics, and cancer.