Massively parallel phenotyping of coding variants in cancer with Perturb-seq.

Genome sequencing studies have identified millions of somatic variants in cancer, but it remains challenging to predict the phenotypic impact of most. Experimental approaches to distinguish impactful variants often use phenotypic assays that report on predefined gene-specific functional effects in bulk cell populations. Here, we develop an approach to functionally assess variant impact in single cells by pooled Perturb-seq. We measured the impact of 200 TP53 and KRAS variants on RNA profiles in over 300,000 single lung cancer cells, and used the profiles to categorize variants into phenotypic subsets to distinguish gain-of-function, loss-of-function and dominant negative variants, which we validated by comparison with orthogonal assays. We discovered that KRAS variants did not merely fit into discrete functional categories, but spanned a continuum of gain-of-function phenotypes, and that their functional impact could not have been predicted solely by their frequency in patient cohorts. Our work provides a scalable, gene-agnostic method for coding variant impact phenotyping, with potential applications in multiple disease settings.

Stem-like intestinal Th17 cells give rise to pathogenic effector T cells during autoimmunity.

While intestinal Th17 cells are critical for maintaining tissue homeostasis, recent studies have implicated their roles in the development of extra-intestinal autoimmune diseases including multiple sclerosis. However, the mechanisms by which tissue Th17 cells mediate these dichotomous functions remain unknown. Here, we characterized the heterogeneity, plasticity, and migratory phenotypes of tissue Th17 cells in vivo by combined fate mapping with profiling of the transcriptomes and TCR clonotypes of over 84,000 Th17 cells at homeostasis and during CNS autoimmune inflammation. Inter- and intra-organ single-cell analyses revealed a homeostatic, stem-like TCF1+ IL-17+ SLAMF6+ population that traffics to the intestine where it is maintained by the microbiota, providing a ready reservoir for the IL-23-driven generation of encephalitogenic GM-CSF+ IFN-γ+ CXCR6+ T cells. Our study defines a direct in vivo relationship between IL-17+ non-pathogenic and GM-CSF+ and IFN-γ+ pathogenic Th17 populations and provides a mechanism by which homeostatic intestinal Th17 cells direct extra-intestinal autoimmune disease.

Cycling cancer persister cells arise from lineages with distinct programs.

Non-genetic mechanisms have recently emerged as important drivers of cancer therapy failure1, where some cancer cells can enter a reversible drug-tolerant persister state in response to treatment2. Although most cancer persisters remain arrested in the presence of the drug, a rare subset can re-enter the cell cycle under constitutive drug treatment. Little is known about the non-genetic mechanisms that enable cancer persisters to maintain proliferative capacity in the presence of drugs. To study this rare, transiently resistant, proliferative persister population, we developed Watermelon, a high-complexity expressed barcode lentiviral library for simultaneous tracing of each cell's clonal origin and proliferative and transcriptional states. Here we show that cycling and non-cycling persisters arise from different cell lineages with distinct transcriptional and metabolic programs. Upregulation of antioxidant gene programs and a metabolic shift to fatty acid oxidation are associated with persister proliferative capacity across multiple cancer types. Impeding oxidative stress or metabolic reprogramming alters the fraction of cycling persisters. In human tumours, programs associated with cycling persisters are induced in minimal residual disease in response to multiple targeted therapies. The Watermelon system enabled the identification of rare persister lineages that are preferentially poised to proliferate under drug pressure, thus exposing new vulnerabilities that can be targeted to delay or even prevent disease recurrence.

Scalable, multimodal profiling of chromatin accessibility, gene expression and protein levels in single cells.

Recent technological advances have enabled massively parallel chromatin profiling with scATAC-seq (single-cell assay for transposase accessible chromatin by sequencing). Here we present ATAC with select antigen profiling by sequencing (ASAP-seq), a tool to simultaneously profile accessible chromatin and protein levels. Our approach pairs sparse scATAC-seq data with robust detection of hundreds of cell surface and intracellular protein markers and optional capture of mitochondrial DNA for clonal tracking, capturing three distinct modalities in single cells. ASAP-seq uses a bridging approach that repurposes antibody:oligonucleotide conjugates designed for existing technologies that pair protein measurements with single-cell RNA sequencing. Together with DOGMA-seq, an adaptation of CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) for measuring gene activity across the central dogma of gene regulation, we demonstrate the utility of systematic multi-omic profiling by revealing coordinated and distinct changes in chromatin, RNA and surface proteins during native hematopoietic differentiation and peripheral blood mononuclear cell stimulation and as a combinatorial decoder and reporter of multiplexed perturbations in primary T cells.

Multimodal pooled Perturb-CITE-seq screens in patient models define mechanisms of cancer immune evasion.

Resistance to immune checkpoint inhibitors (ICIs) is a key challenge in cancer therapy. To elucidate underlying mechanisms, we developed Perturb-CITE-sequencing (Perturb-CITE-seq), enabling pooled clustered regularly interspaced short palindromic repeat (CRISPR)-Cas9 perturbations with single-cell transcriptome and protein readouts. In patient-derived melanoma cells and autologous tumor-infiltrating lymphocyte (TIL) co-cultures, we profiled transcriptomes and 20 proteins in ~218,000 cells under ~750 perturbations associated with cancer cell-intrinsic ICI resistance (ICR). We recover known mechanisms of resistance, including defects in the interferon-γ (IFN-γ)-JAK/STAT and antigen-presentation pathways in RNA, protein and perturbation space, and new ones, including loss/downregulation of CD58. Loss of CD58 conferred immune evasion in multiple co-culture models and was downregulated in tumors of melanoma patients with ICR. CD58 protein expression was not induced by IFN-γ signaling, and CD58 loss conferred immune evasion without compromising major histocompatibility complex (MHC) expression, suggesting that it acts orthogonally to known mechanisms of ICR. This work provides a framework for the deciphering of complex mechanisms by large-scale perturbation screens with multimodal, single-cell readouts, and discovers potentially clinically relevant mechanisms of immune evasion.

RNA-guided transcriptional silencing in vivo with S. aureus CRISPR-Cas9 repressors.

CRISPR-Cas9 transcriptional repressors have emerged as robust tools for disrupting gene regulation in vitro but have not yet been adapted for systemic delivery in adult animal models. Here we describe a Staphylococcus aureus Cas9-based repressor (dSaCas9KRAB) compatible with adeno-associated viral (AAV) delivery. To evaluate dSaCas9KRAB efficacy for gene silencing in vivo, we silenced transcription of Pcsk9, a regulator of cholesterol levels, in the liver of adult mice. Systemic administration of a dual-vector AAV8 system expressing dSaCas9KRAB and a Pcsk9-targeting guide RNA (gRNA) results in significant reductions of serum Pcsk9 and cholesterol levels. Despite a moderate host response to dSaCas9KRAB expression, Pcsk9 repression is maintained for 24 weeks after a single treatment, demonstrating the potential for long-term gene silencing in post-mitotic tissues with dSaCas9KRAB. In vivo programmable gene silencing enables studies that link gene regulation to complex phenotypes and expands the CRISPR-Cas9 perturbation toolbox for basic research and gene therapy applications.

* CRISPR-Based Epigenome Editing of Cytokine Receptors for the Promotion of Cell Survival and Tissue Deposition in Inflammatory Environments.

Musculoskeletal diseases have been associated with inflammatory cytokine action, particularly action by TNF-α and IL-1ß. These inflammatory cytokines promote apoptosis and senescence of cells in diseased tissue and extracellular matrix breakdown. Stem cell-based therapies are being considered for the treatment of musculoskeletal diseases, but the presence of these inflammatory cytokines will have similar deleterious action on therapeutic cells delivered to these environments. Methods that prevent inflammatory-induced apoptosis and proinflammatory signaling, in cell and pathway-specific manners are needed. In this study we demonstrate the use of clustered regularly interspaced short palindromic repeats (CRISPR)-based epigenome editing to alter cell response to inflammatory environments by repressing inflammatory cytokine cell receptors, specifically TNFR1 and IL1R1. We targeted CRISPR/Cas9-based repressors to TNFR1 and IL1R1 gene regulatory elements in human adipose-derived stem cells (hADSCs) and investigated the functional outcomes of repression of these genes. Efficient signaling regulation was demonstrated in engineered hADSCs, as activity of the downstream transcription factor NF-κB was significantly reduced or maintained at baseline levels in the presence of TNF-α or IL-1ß. Pellet culture of undifferentiated hADSCs demonstrated improved survival in engineered hADSCs treated with TNF-α or IL-1ß, while having little effect on their immunomodulatory properties. Furthermore, engineered hADSCs demonstrated improved chondrogenic differentiation capacity in the presence of TNF-α or IL-1ß, as shown by superior production of glycosaminglycans in this inflammatory environment. Overall this work demonstrates a novel method for modulating cell response to inflammatory signaling that has applications in engineering cells delivered to inflammatory environments, and as a direct gene therapy to protect endogenous cells exposed to chronic inflammation, as observed in a broad spectrum of degenerative musculoskeletal pathology.

In vivo genome editing improves muscle function in a mouse model of Duchenne muscular dystrophy.

Duchenne muscular dystrophy (DMD) is a devastating disease affecting about 1 out of 5000 male births and caused by mutations in the dystrophin gene. Genome editing has the potential to restore expression of a modified dystrophin gene from the native locus to modulate disease progression. In this study, adeno-associated virus was used to deliver the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system to the mdx mouse model of DMD to remove the mutated exon 23 from the dystrophin gene. This includes local and systemic delivery to adult mice and systemic delivery to neonatal mice. Exon 23 deletion by CRISPR-Cas9 resulted in expression of the modified dystrophin gene, partial recovery of functional dystrophin protein in skeletal myofibers and cardiac muscle, improvement of muscle biochemistry, and significant enhancement of muscle force. This work establishes CRISPR-Cas9-based genome editing as a potential therapy to treat DMD.