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.

Propensity of undulatory swimmers, such as worms, to go against the flow.

The ability to orient oneself in response to environmental cues is crucial to the survival and function of diverse organisms. One such orientation behavior is the alignment of aquatic organisms with (negative rheotaxis) or against (positive rheotaxis) fluid current. The questions of whether low-Reynolds-number, undulatory swimmers, such as worms, rheotax and whether rheotaxis is a deliberate or an involuntary response to mechanical forces have been the subject of conflicting reports. To address these questions, we use Caenorhabditis elegans as a model undulatory swimmer and examine, in experiment and theory, the orientation of C. elegans in the presence of flow. We find that when close to a stationary surface the animal aligns itself against the direction of the flow. We elucidate for the first time to our knowledge the mechanisms of rheotaxis in worms and show that rheotaxis can be explained solely by mechanical forces and does not require sensory input or deliberate action. The interaction between the flow field induced by the swimmer and a nearby surface causes the swimmer to tilt toward the surface and the velocity gradient associated with the flow rotates the animal to face upstream. Fluid mechanical computer simulations faithfully mimic the behavior observed in experiments, supporting the notion that rheotaxis behavior can be fully explained by hydrodynamics. Our study highlights the important role of hydrodynamics in the behavior of small undulating swimmers and may assist in developing control strategies to affect the animals' life cycles.

Gait synchronization in Caenorhabditis elegans.

Collective motion is observed in swarms of swimmers of various sizes, ranging from self-propelled nanoparticles to fish. The mechanisms that govern interactions among individuals are debated, and vary from one species to another. Although the interactions among relatively large animals, such as fish, are controlled by their nervous systems, the interactions among microorganisms, which lack nervous systems, are controlled through physical and chemical pathways. Little is known, however, regarding the mechanism of collective movements in microscopic organisms with nervous systems. To attempt to remedy this, we studied collective swimming behavior in the nematode Caenorhabditis elegans, a microorganism with a compact nervous system. We evaluated the contributions of hydrodynamic forces, contact forces, and mechanosensory input to the interactions among individuals. We devised an experiment to examine pair interactions as a function of the distance between the animals and observed that gait synchronization occurred only when the animals were in close proximity, independent of genes required for mechanosensation. Our measurements and simulations indicate that steric hindrance is the dominant factor responsible for motion synchronization in C. elegans, and that hydrodynamic interactions and genotype do not play a significant role. We infer that a similar mechanism may apply to other microscopic swimming organisms and self-propelled particles.

Orienting actin filaments for directional motility of processive myosin motors.

To utilize molecular motors in manmade systems, it is necessary to control the motors' motion. We describe a technique to orient actin filaments so that their barbed ends point in the same direction, enabling same-type motors to travel unidirectionally. Myosin-V and myosin-VI were observed to travel, respectively, toward and away from the filaments' barbed ends. When both motors were present, they occasionally passed each other while "walking" in opposite directions along single actin filaments.

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.

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.

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.

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.

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.

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.

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.

Impact of exercise on seniors' motor control response to external dynamics.

Research on postural stability, motor control, and fall occurrence in seniors is common, but few studies address the influence of exercise and external dynamics on elderly balance. Using pre- and post-training tests, the effects of a Fitball exercise program on performance in eight subjects was documented. The exercise program focused on improving dynamic balance and postural stability of seniors. To evaluate progress-related changes, pre- and post-tests in a dynamic environment were applied. Center of gravity (COG) excursion, catch success rate, and balance success rate were quantified, and synchronized data collection of 3D motion capture (VICON v8i) and ground reaction force (2 KISTLER platforms) was analyzed. During pre- and post-tests, participants stood in a walklike stance and were asked to catch a weighted ball, which dropped unexpectedly. Results showed no significant changes in balance success rate. Significant improvements were found, however, in both COG control and catch success rate following training (p<0.05).

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 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.

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.