Single-cell characterization of macrophages in glioblastoma reveals MARCO as a mesenchymal pro-tumor marker.

BACKGROUND: Macrophages are the most common infiltrating immune cells in gliomas and play a wide variety of pro-tumor and anti-tumor roles. However, the different subpopulations of macrophages and their effects on the tumor microenvironment remain poorly understood. METHODS: We combined new and previously published single-cell RNA-seq data from 98,015 single cells from a total of 66 gliomas to profile 19,331 individual macrophages. RESULTS: Unsupervised clustering revealed a pro-tumor subpopulation of bone marrow-derived macrophages characterized by the scavenger receptor MARCO, which is almost exclusively found in IDH1-wild-type glioblastomas. Previous studies have implicated MARCO as an unfavorable marker in melanoma and non-small cell lung cancer; here, we find that bulk MARCO expression is associated with worse prognosis and mesenchymal subtype. Furthermore, MARCO expression is significantly altered over the course of treatment with anti-PD1 checkpoint inhibitors in a response-dependent manner, which we validate with immunofluorescence imaging. CONCLUSIONS: These findings illustrate a novel macrophage subpopulation that drives tumor progression in glioblastomas and suggest potential therapeutic targets to prevent their recruitment.

Pathway-based classification of glioblastoma uncovers a mitochondrial subtype with therapeutic vulnerabilities.

The transcriptomic classification of glioblastoma (GBM) has failed to predict survival and therapeutic vulnerabilities. A computational approach for unbiased identification of core biological traits of single cells and bulk tumors uncovered four tumor cell states and GBM subtypes distributed along neurodevelopmental and metabolic axes, classified as proliferative/progenitor, neuronal, mitochondrial and glycolytic/plurimetabolic. Each subtype was enriched with biologically coherent multiomic features. Mitochondrial GBM was associated with the most favorable clinical outcome. It relied exclusively on oxidative phosphorylation for energy production, whereas the glycolytic/plurimetabolic subtype was sustained by aerobic glycolysis and amino acid and lipid metabolism. Deletion of the glucose-proton symporter SLC45A1 was the truncal alteration most significantly associated with mitochondrial GBM, and the reintroduction of SLC45A1 in mitochondrial glioma cells induced acidification and loss of fitness. Mitochondrial, but not glycolytic/plurimetabolic, GBM exhibited marked vulnerability to inhibitors of oxidative phosphorylation. The pathway-based classification of GBM informs survival and enables precision targeting of cancer metabolism.

ERK1/2 phosphorylation predicts survival following anti-PD-1 immunotherapy in recurrent glioblastoma.

Only a subset of recurrent glioblastoma (rGBM) responds to anti-PD-1 immunotherapy. Previously, we reported enrichment of BRAF/PTPN11 mutations in 30% of rGBM that responded to PD-1 blockade. Given that BRAF and PTPN11 promote MAPK/ERK signaling, we investigated whether activation of this pathway is associated with response to PD-1 inhibitors in rGBM, including patients that do not harbor BRAF/PTPN11 mutations. Here we show that immunohistochemistry for ERK1/2 phosphorylation (p-ERK), a marker of MAPK/ERK pathway activation, is predictive of overall survival following adjuvant PD-1 blockade in two independent rGBM patient cohorts. Single-cell RNA-sequencing and multiplex immunofluorescence analyses revealed that p-ERK was mainly localized in tumor cells and that high-p-ERK GBMs contained tumor-infiltrating myeloid cells and microglia with elevated expression of MHC class II and associated genes. These findings indicate that ERK1/2 activation in rGBM is predictive of response to PD-1 blockade and is associated with a distinct myeloid cell phenotype.

A map of tumor-host interactions in glioma at single-cell resolution.

BACKGROUND: Single-cell RNA sequencing is the reference technique for characterizing the heterogeneity of the tumor microenvironment. The composition of the various cell types making up the microenvironment can significantly affect the way in which the immune system activates cancer rejection mechanisms. Understanding the cross-talk signals between immune cells and cancer cells is of fundamental importance for the identification of immuno-oncology therapeutic targets. RESULTS: We present a novel method, single-cell Tumor-Host Interaction tool (scTHI), to identify significantly activated ligand-receptor interactions across clusters of cells from single-cell RNA sequencing data. We apply our approach to uncover the ligand-receptor interactions in glioma using 6 publicly available human glioma datasets encompassing 57,060 gene expression profiles from 71 patients. By leveraging this large-scale collection we show that unexpected cross-talk partners are highly conserved across different datasets in the majority of the tumor samples. This suggests that shared cross-talk mechanisms exist in glioma. CONCLUSIONS: Our results provide a complete map of the active tumor-host interaction pairs in glioma that can be therapeutically exploited to reduce the immunosuppressive action of the microenvironment in brain tumor.

Integrating single-cell RNA-seq and imaging with SCOPE-seq2.

Live cell imaging allows direct observation and monitoring of phenotypes that are difficult to infer from transcriptomics. However, existing methods for linking microscopy and single-cell RNA-seq (scRNA-seq) have limited scalability. Here, we describe an upgraded version of Single Cell Optical Phenotyping and Expression (SCOPE-seq2) for combining single-cell imaging and expression profiling, with substantial improvements in throughput, molecular capture efficiency, linking accuracy, and compatibility with standard microscopy instrumentation. We introduce improved optically decodable mRNA capture beads and implement a more scalable and simplified optical decoding process. We demonstrate the utility of SCOPE-seq2 for fluorescence, morphological, and expression profiling of individual primary cells from a human glioblastoma (GBM) surgical sample, revealing relationships between simple imaging features and cellular identity, particularly among malignantly transformed tumor cells.

Single-Cell Profiling and SCOPE-Seq Reveal Lineage Dynamics of Adult Ventricular-Subventricular Zone Neurogenesis and NOTUM as a Key Regulator.

In the adult ventricular-subventricular zone (V-SVZ), neural stem cells (NSCs) generate new olfactory bulb (OB) neurons and glia throughout life. To map adult neuronal lineage progression, we profiled >56,000 V-SVZ and OB cells by single-cell RNA sequencing (scRNA-seq). Our analyses reveal the molecular diversity of OB neurons, including fate-mapped neurons, lineage progression dynamics, and an NSC intermediate enriched for Notum, which encodes a secreted WNT antagonist. SCOPE-seq technology, which links live-cell imaging with scRNA-seq, uncovers cell-size transitions during NSC differentiation and preferential NOTUM binding to proliferating neuronal precursors. Consistently, application of NOTUM protein in slice cultures and pharmacological inhibition of NOTUM in slice cultures and in vivo demonstrated that NOTUM negatively regulates V-SVZ proliferation. Timely, context-dependent neurogenesis demands adaptive signaling among neighboring progenitors. Our findings highlight a critical regulatory state during NSC activation marked by NOTUM, which attenuates WNT-stimulated proliferation in NSC progeny.

Single-cell transcriptomics of human T cells reveals tissue and activation signatures in health and disease.

Human T cells coordinate adaptive immunity in diverse anatomic compartments through production of cytokines and effector molecules, but it is unclear how tissue site influences T cell persistence and function. Here, we use single cell RNA-sequencing (scRNA-seq) to define the heterogeneity of human T cells isolated from lungs, lymph nodes, bone marrow and blood, and their functional responses following stimulation. Through analysis of >50,000 resting and activated T cells, we reveal tissue T cell signatures in mucosal and lymphoid sites, and lineage-specific activation states across all sites including distinct effector states for CD8+ T cells and an interferon-response state for CD4+ T cells. Comparing scRNA-seq profiles of tumor-associated T cells to our dataset reveals predominant activated CD8+ compared to CD4+ T cell states within multiple tumor types. Our results therefore establish a high dimensional reference map of human T cell activation in health for analyzing T cells in disease.

Author Correction: Immune and genomic correlates of response to anti-PD-1 immunotherapy in glioblastoma.

In the version of this article originally published, the graph in Extended Data Fig. 2c was a duplication of Extended Data Fig. 2b. The correct version of Extended Data Fig. 2c is now available online.

De novo gene signature identification from single-cell RNA-seq with hierarchical Poisson factorization.

Common approaches to gene signature discovery in single-cell RNA-sequencing (scRNA-seq) depend upon predefined structures like clusters or pseudo-temporal order, require prior normalization, or do not account for the sparsity of single-cell data. We present single-cell hierarchical Poisson factorization (scHPF), a Bayesian factorization method that adapts hierarchical Poisson factorization (Gopalan et al, 2015, Proceedings of the 31st Conference on Uncertainty in Artificial Intelligence, 326) for de novo discovery of both continuous and discrete expression patterns from scRNA-seq. scHPF does not require prior normalization and captures statistical properties of single-cell data better than other methods in benchmark datasets. Applied to scRNA-seq of the core and margin of a high-grade glioma, scHPF uncovers marked differences in the abundance of glioma subpopulations across tumor regions and regionally associated expression biases within glioma subpopulations. scHFP revealed an expression signature that was spatially biased toward the glioma-infiltrated margins and associated with inferior survival in glioblastoma.

Immune and genomic correlates of response to anti-PD-1 immunotherapy in glioblastoma.

Immune checkpoint inhibitors have been successful across several tumor types; however, their efficacy has been uncommon and unpredictable in glioblastomas (GBM), where <10% of patients show long-term responses. To understand the molecular determinants of immunotherapeutic response in GBM, we longitudinally profiled 66 patients, including 17 long-term responders, during standard therapy and after treatment with PD-1 inhibitors (nivolumab or pembrolizumab). Genomic and transcriptomic analysis revealed a significant enrichment of PTEN mutations associated with immunosuppressive expression signatures in non-responders, and an enrichment of MAPK pathway alterations (PTPN11, BRAF) in responders. Responsive tumors were also associated with branched patterns of evolution from the elimination of neoepitopes as well as with differences in T cell clonal diversity and tumor microenvironment profiles. Our study shows that clinical response to anti-PD-1 immunotherapy in GBM is associated with specific molecular alterations, immune expression signatures, and immune infiltration that reflect the tumor's clonal evolution during treatment.

Single-Cell Analysis of Regional Differences in Adult V-SVZ Neural Stem Cell Lineages.

The ventricular-subventricular zone (V-SVZ) harbors adult neural stem cells. V-SVZ neural stem cells exhibit features of astrocytes, have a regional identity, and depending on their location in the lateral or septal wall of the lateral ventricle, generate different types of neuronal and glial progeny. We performed large-scale single-cell RNA sequencing to provide a molecular atlas of cells from the lateral and septal adult V-SVZ of male and female mice. This revealed regional and sex differences among adult V-SVZ cells. We uncovered lineage potency bias at the single-cell level among lateral and septal wall astrocytes toward neurogenesis and oligodendrogenesis, respectively. Finally, we identified transcription factor co-expression modules marking key temporal steps in neurogenic and oligodendrocyte lineage progression. Our data suggest functionally important spatial diversity in neurogenesis and oligodendrogenesis in the adult brain and reveal molecular correlates of adult NSC dormancy and lineage specialization.

SCOPE-Seq: a scalable technology for linking live cell imaging and single-cell RNA sequencing.

Optically decodable beads link the identity of a sample to a measurement through an optical barcode, enabling libraries of biomolecules to be captured on beads in solution and decoded by fluorescence. This approach has been foundational to microarray, sequencing, and flow-based expression profiling technologies. We combine microfluidics with optically decodable beads and show that phenotypic analysis of living cells can be linked to single-cell sequencing. As a proof-of-concept, we demonstrate the accuracy and scalability of our tool called Single Cell Optical Phenotyping and Expression sequencing (SCOPE-Seq) to combine live cell imaging with single-cell RNA sequencing.

Single-cell transcriptome analysis of lineage diversity in high-grade glioma.

BACKGROUND: Despite extensive molecular characterization, we lack a comprehensive understanding of lineage identity, differentiation, and proliferation in high-grade gliomas (HGGs). METHODS: We sampled the cellular milieu of HGGs by profiling dissociated human surgical specimens with a high-density microwell system for massively parallel single-cell RNA-Seq. We analyzed the resulting profiles to identify subpopulations of both HGG and microenvironmental cells and applied graph-based methods to infer structural features of the malignantly transformed populations. RESULTS: While HGG cells can resemble glia or even immature neurons and form branched lineage structures, mesenchymal transformation results in unstructured populations. Glioma cells in a subset of mesenchymal tumors lose their neural lineage identity, express inflammatory genes, and co-exist with marked myeloid infiltration, reminiscent of molecular interactions between glioma and immune cells established in animal models. Additionally, we discovered a tight coupling between lineage resemblance and proliferation among malignantly transformed cells. Glioma cells that resemble oligodendrocyte progenitors, which proliferate in the brain, are often found in the cell cycle. Conversely, glioma cells that resemble astrocytes, neuroblasts, and oligodendrocytes, which are non-proliferative in the brain, are generally non-cycling in tumors. CONCLUSIONS: These studies reveal a relationship between cellular identity and proliferation in HGG and distinct population structures that reflects the extent of neural and non-neural lineage resemblance among malignantly transformed cells.

Single-Cell Transcriptomic Analysis of Tumor Heterogeneity.

Intratumoral heterogeneity is among the greatest challenges in precision cancer therapy. However, developments in high-throughput single-cell RNA sequencing (scRNA-seq) may now provide the statistical power to dissect the diverse cellular populations of tumors. In the future these technologies might inform the selection of targeted combination therapies and enrollment criteria for clinical trials.

Terrain following and applications: Caenorhabditis elegans swims along the floor using a bump and undulate strategy.

Nematodes such as Caenorhabditis elegans are heavier than water. When submerged in water, they settle to the bottom surface. Observations reveal that the animals do not lie flat on the bottom surface, but remain substantially suspended above the surface through continuous collisions with the surface, while maintaining their swimming gaits. Consequently, the swimming animals follow the bottom surface topography. When the bottom surface is inclined, the animals swim up or down along the incline. As the magnitude of the gravitational force can be easily estimated, this behaviour provides a convenient means to estimate the animal's propulsive thrust. The animals' tendency to follow the surface topography provides a means to control the swimmers' trajectories and direction of motion, which we demonstrate with a saw tooth-like ratchet that biases the animals to swim in a selected direction. The animals can also serve as surface topography probes since their residence time as a function of position provides information on surface features. Finally, we take advantage of surface following to construct a simple motility-based sorter that can sort animals based on genotype and state of health.

An Automated Microwell Platform for Large-Scale Single Cell RNA-Seq.

Recent developments have enabled rapid, inexpensive RNA sequencing of thousands of individual cells from a single specimen, raising the possibility of unbiased and comprehensive expression profiling from complex tissues. Microwell arrays are a particularly attractive microfluidic platform for single cell analysis due to their scalability, cell capture efficiency, and compatibility with imaging. We report an automated microwell array platform for single cell RNA-Seq with significantly improved performance over previous implementations. We demonstrate cell capture efficiencies of >50%, compatibility with commercially available barcoded mRNA capture beads, and parallel expression profiling from thousands of individual cells. We evaluate the level of cross-contamination in our platform by both tracking fluorescent cell lysate in sealed microwells and with a human-mouse mixed species RNA-Seq experiment. Finally, we apply our system to comprehensively assess heterogeneity in gene expression of patient-derived glioma neurospheres and uncover subpopulations similar to those observed in human glioma tissue.

A hydrodynamic mechanism for attraction of undulatory microswimmers to surfaces (bordertaxis).

Although small nematodes significantly impact human and animal health, agriculture, and ecology, little is known about the role of hydrodynamics in their life cycles. Using the nematode Caenorhabditis elegans as a model undulatory microswimmer, we have observed that animals are attracted to and swim along surfaces. The attraction to surfaces does not require mechanosensory neuron function. In dilute swarms, swimmers aggregate near surfaces. Using resistive force-based theory, symmetry arguments, and direct hydrodynamic simulations, we demonstrate that forces resulting from the interaction between the swimmer-induced flow field and a nearby surface cause a short-range hydrodynamic torque that stirs the swimmers towards the surface. When combined with steric forces, this causes locomotion along the surface. This surface attraction may affect nematode mate and food finding behaviour and, in the case of parasitic nematodes, may facilitate host penetration. Surface attraction must be accounted for when studying animals' responses to various stimuli, and suggests means of controlling undulatory microswimmers.

High-throughput, motility-based sorter for microswimmers such as C. elegans.

Animal motility varies with genotype, disease, aging, and environmental conditions. In many studies, it is desirable to carry out high throughput motility-based sorting to isolate rare animals for, among other things, forward genetic screens to identify genetic pathways that regulate phenotypes of interest. Many commonly used screening processes are labor-intensive, lack sensitivity, and require extensive investigator training. Here, we describe a sensitive, high throughput, automated, motility-based method for sorting nematodes. Our method is implemented in a simple microfluidic device capable of sorting thousands of animals per hour per module, and is amenable to parallelism. The device successfully enriches for known C. elegans motility mutants. Furthermore, using this device, we isolate low-abundance mutants capable of suppressing the somnogenic effects of the flp-13 gene, which regulates C. elegans sleep. By performing genetic complementation tests, we demonstrate that our motility-based sorting device efficiently isolates mutants for the same gene identified by tedious visual inspection of behavior on an agar surface. Therefore, our motility-based sorter is capable of performing high throughput gene discovery approaches to investigate fundamental biological processes.