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

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.

BackgroundDespite extensive molecular characterization, we lack a comprehensive understanding of lineage identity, differentiation, and proliferation in high-grade gliomas (HGGs).MethodsWe 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.ResultsWhile 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.ConclusionsThese 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.Electronic supplementary materialThe online version of this article (10.1186/s13073-018-0567-9) contains supplementary material, which is available to authorized users.

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

No abstract

Common approaches to gene signature discovery in single‐cell RNA ‐sequencing (sc RNA ‐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 (sc HPF ), 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 sc RNA ‐seq. sc HPF does not require prior normalization and captures statistical properties of single‐cell data better than other methods in benchmark datasets. Applied to sc RNA ‐seq of the core and margin of a high‐grade glioma, sc HPF uncovers marked differences in the abundance of glioma subpopulations across tumor regions and regionally associated expression biases within glioma subpopulations. sc HFP revealed an expression signature that was spatially biased toward the glioma‐infiltrated margins and associated with inferior survival in glioblastoma.

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.

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.