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. Keywords: single-cell RNA-seq; pooled screen; CRISPR; epistasis; genetic interaction

Cancer immunogenomic approach to neoantigen discovery in a checkpoint blockade responsive murine model of oral cavity squamous cell carcinoma

Head and neck squamous cell carcinomas (HNSCC) are an ideal immunotherapy target due to their high mutation burden and frequent infiltration with lymphocytes. Preclinical models to investigate targeted and combination therapies as well as defining biomarkers to guide treatment represent an important need in the field. Immunogenomics approaches have illuminated the role of mutation-derived tumor neoantigens as potential biomarkers of response to checkpoint blockade as well as representing therapeutic vaccines. Here, we aimed to define a platform for checkpoint and other immunotherapy studies using syngeneic HNSCC cell line models (MOC2 and MOC22), and evaluated the association between mutation burden, predicted neoantigen landscape, infiltrating T cell populations and responsiveness of tumors to anti-PD1 therapy. We defined dramatic hematopoietic cell transcriptomic alterations in the MOC22 anti-PD1 responsive model in both tumor and draining lymph nodes. Using a cancer immunogenomics pipeline and validation with ELISPOT and tetramer analysis, we identified the H-2Kb-restricted ICAM1P315L (mICAM1) as a neoantigen in MOC22. Finally, we demonstrated that mICAM1 vaccination was able to protect against MOC22 tumor development defining mICAM1 as a bona fide neoantigen. Together these data define a pre-clinical HNSCC model system that provides a foundation for future investigations into combination and novel therapeutics

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

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 co-inhibitory receptor expression on CD4+ T cells promotes autoimmunity while sustained overexpression on CD8+ T cells promotes T cell dysfunction or exhaustion, leading to impaired ability to clear chronic viral infections and cancer1,2. Here, we used RNA and protein expression profiling at single-cell resolution to 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 a number of novel surface receptors. We functionally validated two novel 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 multiple 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, which we validated experimentally. This molecular circuit underlies the co-expression of co-inhibitory receptors in T cells and identifies novel regulators of T cell function with the potential to regulate autoimmunity and tumor immunity

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.Advanced Medical Research Foundation (Brookline, Mass.)Breast Cancer Research FoundationNational Institutes of Health (U.S.) (NIH grant HG002668)National Institutes of Health (U.S.) (NIH grant CA146445

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. Maus, Regev, and colleagues performed a deep transcriptional analysis of CAR T cells bearing different signaling domains, at rest and after activation. Findings include CAR-specific signatures and a distinct program activated in CARs with 4-1BB domains that was characterized by expression of central memory markers, MHC class II, and IL-21, among others

Phenotypes Determined by Cluster Analysis and Their Survival in the Prospective European Scleroderma Trials and Research Cohort of Patients With Systemic Sclerosis

Objective Systemic sclerosis (SSc) is a heterogeneous connective tissue disease that is typically subdivided into limited cutaneous SSc (lcSSc) and diffuse cutaneous SSc (dcSSc) depending on the extent of skin involvement. This subclassification may not capture the entire variability of clinical phenotypes. The European Scleroderma Trials and Research (EUSTAR) database includes data on a prospective cohort of SSc patients from 122 European referral centers. This study was undertaken to perform a cluster analysis of EUSTAR data to distinguish and characterize homogeneous phenotypes without any a priori assumptions, and to examine survival among the clusters obtained. Methods A total of 11,318 patients were registered in the EUSTAR database, and 6,927 were included in the study. Twenty-four clinical and serologic variables were used for clustering. Results Clustering analyses provided a first delineation of 2 clusters showing moderate stability. In an exploratory attempt, we further characterized 6 homogeneous groups that differed with regard to their clinical features, autoantibody profile, and mortality. Some groups resembled usual dcSSc or lcSSc prototypes, but others exhibited unique features, such as a majority of lcSSc patients with a high rate of visceral damage and antitopoisomerase antibodies. Prognosis varied among groups and the presence of organ damage markedly impacted survival regardless of cutaneous involvement. Conclusion Our findings suggest that restricting subsets of SSc patients to only those based on cutaneous involvement may not capture the complete heterogeneity of the disease. Organ damage and antibody profile should be taken into consideration when individuating homogeneous groups of patients with a distinct prognosis