Perturb-Seq: Dissecting Molecular Circuits with Scalable Single-Cell RNA Profiling of Pooled Genetic Screens.

Genetic screens help infer gene function in mammalian cells, but it has remained difficult to assay complex phenotypes-such as transcriptional profiles-at scale. Here, we develop Perturb-seq, combining single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-based perturbations to perform many such assays in a pool. We demonstrate Perturb-seq by analyzing 200,000 cells in immune cells and cell lines, focusing on transcription factors regulating the response of dendritic cells to lipopolysaccharide (LPS). Perturb-seq accurately identifies individual gene targets, gene signatures, and cell states affected by individual perturbations and their genetic interactions. We posit new functions for regulators of differentiation, the anti-viral response, and mitochondrial function during immune activation. By decomposing many high content measurements into the effects of perturbations, their interactions, and diverse cell metadata, Perturb-seq dramatically increases the scope of pooled genomic assays.

A Distinct Gene Module for Dysfunction Uncoupled from Activation in Tumor-Infiltrating T Cells.

Reversing the dysfunctional T cell state that arises in cancer and chronic viral infections is the focus of therapeutic interventions; however, current therapies are effective in only some patients and some tumor types. To gain a deeper molecular understanding of the dysfunctional T cell state, we analyzed population and single-cell RNA profiles of CD8(+) tumor-infiltrating lymphocytes (TILs) and used genetic perturbations to identify a distinct gene module for T cell dysfunction that can be uncoupled from T cell activation. This distinct dysfunction module is downstream of intracellular metallothioneins that regulate zinc metabolism and can be identified at single-cell resolution. We further identify Gata-3, a zinc-finger transcription factor in the dysfunctional module, as a regulator of dysfunction, and we use CRISPR-Cas9 genome editing to show that it drives a dysfunctional phenotype in CD8(+) TILs. Our results open novel avenues for targeting dysfunctional T cell states while leaving activation programs intact.

Poised chromatin at the ZEB1 promoter enables breast cancer cell plasticity and enhances tumorigenicity.

The recent discovery that normal and neoplastic epithelial cells re-enter the stem cell state raised the intriguing possibility that the aggressiveness of carcinomas derives not from their existing content of cancer stem cells (CSCs) but from their proclivity to generate new CSCs from non-CSC populations. Here, we demonstrate that non-CSCs of human basal breast cancers are plastic cell populations that readily switch from a non-CSC to CSC state. The observed cell plasticity is dependent on ZEB1, a key regulator of the epithelial-mesenchymal transition. We find that plastic non-CSCs maintain the ZEB1 promoter in a bivalent chromatin configuration, enabling them to respond readily to microenvironmental signals, such as TGFß. In response, the ZEB1 promoter converts from a bivalent to active chromatin configuration, ZEB1 transcription increases, and non-CSCs subsequently enter the CSC state. Our findings support a dynamic model in which interconversions between low and high tumorigenic states occur frequently, thereby increasing tumorigenic and malignant potential.

Induction and transcriptional regulation of the co-inhibitory gene module in T cells.

The expression of co-inhibitory receptors, such as CTLA-4 and PD-1, on effector T cells is a key mechanism for ensuring immune homeostasis. Dysregulated expression of co-inhibitory receptors on CD4+ T cells promotes autoimmunity, whereas sustained overexpression on CD8+ T cells promotes T cell dysfunction or exhaustion, leading to impaired ability to clear chronic viral infections and diseases such as cancer1,2. Here, using RNA and protein expression profiling at single-cell resolution in mouse cells, we identify a module of co-inhibitory receptors that includes not only several known co-inhibitory receptors (PD-1, TIM-3, LAG-3 and TIGIT) but also many new surface receptors. We functionally validated two new co-inhibitory receptors, activated protein C receptor (PROCR) and podoplanin (PDPN). The module of co-inhibitory receptors is co-expressed in both CD4+ and CD8+ T cells and is part of a larger co-inhibitory gene program that is shared by non-responsive T cells in several physiological contexts and is driven by the immunoregulatory cytokine IL-27. Computational analysis identified the transcription factors PRDM1 and c-MAF as cooperative regulators of the co-inhibitory module, and this was validated experimentally. This molecular circuit underlies the co-expression of co-inhibitory receptors in T cells and identifies regulators of T cell function with the potential to control autoimmunity and tumour immunity.

A Distinct Transcriptional Program in Human CAR T Cells Bearing the 4-1BB Signaling Domain Revealed by scRNA-Seq.

T cells engineered to express chimeric antigen receptors (CARs) targeting CD19 have produced impressive outcomes for the treatment of B cell malignancies, but different products vary in kinetics, persistence, and toxicity profiles based on the co-stimulatory domains included in the CAR. In this study, we performed transcriptional profiling of bulk CAR T cell populations and single cells to characterize the transcriptional states of human T cells transduced with CD3ζ, 4-1BB-CD3ζ (BBζ), or CD28-CD3ζ (28ζ) co-stimulatory domains at rest and after activation by triggering their CAR or their endogenous T cell receptor (TCR). We identified a transcriptional signature common across CARs with the CD3ζ signaling domain, as well as a distinct program associated with the 4-1BB co-stimulatory domain at rest and after activation. CAR T cells bearing BBζ had increased expression of human leukocyte antigen (HLA) class II genes, ENPP2, and interleukin (IL)-21 axis genes, and decreased PD1 compared to 28ζ CAR T cells. Similar to previous studies, we also found BBζ CAR CD8 T cells to be enriched in a central memory cell phenotype and fatty acid metabolism genes. Our data uncovered transcriptional signatures related to costimulatory domains and demonstrated that signaling domains included in CARs uniquely shape the transcriptional programs of T cells.

Cell plasticity and heterogeneity in cancer.

BACKGROUND: Heterogeneity within a given cancer arises from diverse cell types recruited to the tumor and from genetic and/or epigenetic differences amongst the cancer cells themselves. These factors conspire to create a disease with various phenotypes. There are 2 established models of cancer development and progression to metastatic disease. These are the clonal evolution and cancer stem cell models. CONTENT: The clonal evolution theory suggests that successive mutations accumulating in a given cell generate clonal outgrowths that thrive in response to microenvironmental selection pressures, dictating the phenotype of the tumor. The alternative cancer stem cell (CSC) model suggests that cancer cells with similar genetic backgrounds can be hierarchically organized according to their tumorigenic potential. Accordingly, CSCs reside at the apex of the hierarchy and are thought to possess the majority of a cancer's tumor-initiating and metastatic ability. A defining feature of this model is its apparent unidirectional nature, whereby CSCs undergo symmetric division to replenish the CSC pool and irreversible asymmetric division to generate daughter cells (non-CSCs) with low tumorigenic potential. However, evolving evidence supports a new model of tumorigenicity, in which considerable plasticity exists between the non-CSC and CSC compartments, such that non-CSCs can reacquire a CSC phenotype. These findings suggest that some tumors may adhere to a plastic CSC model, in which bidirectional conversions are common and essential components of tumorigenicity. SUMMARY: Accumulating evidence surrounding the plasticity of cancer cells, in particular, suggests that aggressive CSCs can be created de novo within a tumor. Given the current focus on therapeutic targeting of CSCs, we discuss the implications of non-CSC-to-CSC conversions on the development of future therapies.

Systematic comparison of single-cell and single-nucleus RNA-sequencing methods.

The scale and capabilities of single-cell RNA-sequencing methods have expanded rapidly in recent years, enabling major discoveries and large-scale cell mapping efforts. However, these methods have not been systematically and comprehensively benchmarked. Here, we directly compare seven methods for single-cell and/or single-nucleus profiling-selecting representative methods based on their usage and our expertise and resources to prepare libraries-including two low-throughput and five high-throughput methods. We tested the methods on three types of samples: cell lines, peripheral blood mononuclear cells and brain tissue, generating 36 libraries in six separate experiments in a single center. To directly compare the methods and avoid processing differences introduced by the existing pipelines, we developed scumi, a flexible computational pipeline that can be used with any single-cell RNA-sequencing method. We evaluated the methods for both basic performance, such as the structure and alignment of reads, sensitivity and extent of multiplets, and for their ability to recover known biological information in the samples.

Author Correction: Systematic comparison of single-cell and single-nucleus RNA-sequencing methods.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.