BACH2 regulates diversification of regulatory and proinflammatory chromatin states in TH17 cells.

Interleukin-17 (IL-17)-producing helper T (TH17) cells are heterogenous and consist of nonpathogenic TH17 (npTH17) cells that contribute to tissue homeostasis and pathogenic TH17 (pTH17) cells that mediate tissue inflammation. Here, we characterize regulatory pathways underlying TH17 heterogeneity and discover substantial differences in the chromatin landscape of npTH17 and pTH17 cells both in vitro and in vivo. Compared to other CD4+ T cell subsets, npTH17 cells share accessible chromatin configurations with regulatory T cells, whereas pTH17 cells exhibit features of both npTH17 cells and type 1 helper T (TH1) cells. Integrating single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) and single-cell RNA sequencing (scRNA-seq), we infer self-reinforcing and mutually exclusive regulatory networks controlling different cell states and predicted transcription factors regulating TH17 cell pathogenicity. We validate that BACH2 promotes immunomodulatory npTH17 programs and restrains proinflammatory TH1-like programs in TH17 cells in vitro and in vivo. Furthermore, human genetics implicate BACH2 in multiple sclerosis. Overall, our work identifies regulators of TH17 heterogeneity as potential targets to mitigate autoimmunity.

Systematically characterizing the roles of E3-ligase family members in inflammatory responses with massively parallel Perturb-seq.

E3 ligases regulate key processes, but many of their roles remain unknown. Using Perturb-seq, we interrogated the function of 1,130 E3 ligases, partners and substrates in the inflammatory response in primary dendritic cells (DCs). Dozens impacted the balance of DC1, DC2, migratory DC and macrophage states and a gradient of DC maturation. Family members grouped into co-functional modules that were enriched for physical interactions and impacted specific programs through substrate transcription factors. E3s and their adaptors co-regulated the same processes, but partnered with different substrate recognition adaptors to impact distinct aspects of the DC life cycle. Genetic interactions were more prevalent within than between modules, and a deep learning model, comßVAE, predicts the outcome of new combinations by leveraging modularity. The E3 regulatory network was associated with heritable variation and aberrant gene expression in immune cells in human inflammatory diseases. Our study provides a general approach to dissect gene function.

Mostly natural sequencing-by-synthesis for scRNA-seq using Ultima sequencing.

Here we introduce a mostly natural sequencing-by-synthesis (mnSBS) method for single-cell RNA sequencing (scRNA-seq), adapted to the Ultima genomics platform, and systematically benchmark it against current scRNA-seq technology. mnSBS uses mostly natural, unmodified nucleotides and only a low fraction of fluorescently labeled nucleotides, which allows for high polymerase processivity and lower costs. We demonstrate successful application in four scRNA-seq case studies of different technical and biological types, including 5' and 3' scRNA-seq, human peripheral blood mononuclear cells from a single individual and in multiplex, as well as Perturb-Seq. Benchmarking shows that results from mnSBS-based scRNA-seq are very similar to those using Illumina sequencing, with minor differences in results related to the position of reads relative to annotated gene boundaries, owing to single-end reads of Ultima being closer to gene ends than reads from Illumina. The method is thus compatible with state-of-the-art scRNA-seq libraries independent of the sequencing technology. We expect mnSBS to be of particular utility for cost-effective large-scale scRNA-seq projects.

Disruption of the IL-33-ST2-AKT signaling axis impairs neurodevelopment by inhibiting microglial metabolic adaptation and phagocytic function.

To accommodate the changing needs of the developing brain, microglia must undergo substantial morphological, phenotypic, and functional reprogramming. Here, we examined whether cellular metabolism regulates microglial function during neurodevelopment. Microglial mitochondria bioenergetics correlated with and were functionally coupled to phagocytic activity in the developing brain. Transcriptional profiling of microglia with diverse metabolic profiles revealed an activation signature wherein the interleukin (IL)-33 signaling axis is associated with phagocytic activity. Genetic perturbation of IL-33 or its receptor ST2 led to microglial dystrophy, impaired synaptic function, and behavioral abnormalities. Conditional deletion of Il33 from astrocytes or Il1rl1, encoding ST2, in microglia increased susceptibility to seizures. Mechanistically, IL-33 promoted mitochondrial activity and phagocytosis in an AKT-dependent manner. Mitochondrial metabolism and AKT activity were temporally regulated in vivo. Thus, a microglia-astrocyte circuit mediated by the IL-33-ST2-AKT signaling axis supports microglial metabolic adaptation and phagocytic function during early development, with implications for neurodevelopmental and neuropsychiatric disorders.

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.

Metabolic modeling of single Th17 cells reveals regulators of autoimmunity.

Metabolism is a major regulator of immune cell function, but it remains difficult to study the metabolic status of individual cells. Here, we present Compass, an algorithm to characterize cellular metabolic states based on single-cell RNA sequencing and flux balance analysis. We applied Compass to associate metabolic states with T helper 17 (Th17) functional variability (pathogenic potential) and recovered a metabolic switch between glycolysis and fatty acid oxidation, akin to known Th17/regulatory T cell (Treg) differences, which we validated by metabolic assays. Compass also predicted that Th17 pathogenicity was associated with arginine and downstream polyamine metabolism. Indeed, polyamine-related enzyme expression was enhanced in pathogenic Th17 and suppressed in Treg cells. Chemical and genetic perturbation of polyamine metabolism inhibited Th17 cytokines, promoted Foxp3 expression, and remodeled the transcriptome and epigenome of Th17 cells toward a Treg-like state. In vivo perturbations of the polyamine pathway altered the phenotype of encephalitogenic T cells and attenuated tissue inflammation in CNS autoimmunity.

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.

Gut CD4+ T cell phenotypes are a continuum molded by microbes, not by TH archetypes.

CD4+ effector lymphocytes (Teff) are traditionally classified by the cytokines they produce. To determine the states that Teff cells actually adopt in frontline tissues in vivo, we applied single-cell transcriptome and chromatin analyses to colonic Teff cells in germ-free or conventional mice or in mice after challenge with a range of phenotypically biasing microbes. Unexpected subsets were marked by the expression of the interferon (IFN) signature or myeloid-specific transcripts, but transcriptome or chromatin structure could not resolve discrete clusters fitting classic helper T cell (TH) subsets. At baseline or at different times of infection, transcripts encoding cytokines or proteins commonly used as TH markers were distributed in a polarized continuum, which was functionally validated. Clones derived from single progenitors gave rise to both IFN-γ- and interleukin (IL)-17-producing cells. Most of the transcriptional variance was tied to the infecting agent, independent of the cytokines produced, and chromatin variance primarily reflected activities of activator protein (AP)-1 and IFN-regulatory factor (IRF) transcription factor (TF) families, not the canonical subset master regulators T-bet, GATA3 or RORγ.

An IL-27-Driven Transcriptional Network Identifies Regulators of IL-10 Expression across T Helper Cell Subsets.

Interleukin-27 (IL-27) is an immunoregulatory cytokine that suppresses inflammation through multiple mechanisms, including induction of IL-10, but the transcriptional network mediating its diverse functions remains unclear. Combining temporal RNA profiling with computational algorithms, we predict 79 transcription factors induced by IL-27 in T cells. We validate 11 known and discover 5 positive (Cebpb, Fosl2, Tbx21, Hlx, and Atf3) and 2 negative (Irf9 and Irf8) Il10 regulators, generating an experimentally refined regulatory network for Il10. We report two central regulators, Prdm1 and Maf, that cooperatively drive the expression of signature genes induced by IL-27 in type 1 regulatory T cells, mediate IL-10 expression in all T helper cells, and determine the regulatory phenotype of colonic Foxp3+ regulatory T cells. Prdm1/Maf double-knockout mice develop spontaneous colitis, phenocopying ll10-deficient mice. Our work provides insights into IL-27-driven transcriptional networks and identifies two shared Il10 regulators that orchestrate immunoregulatory programs across T helper cell subsets.

Calcitonin Gene-Related Peptide Negatively Regulates Alarmin-Driven Type 2 Innate Lymphoid Cell Responses.

Neuroimmune interactions have emerged as critical modulators of allergic inflammation, and type 2 innate lymphoid cells (ILC2s) are an important cell type for mediating these interactions. Here, we show that ILC2s expressed both the neuropeptide calcitonin gene-related peptide (CGRP) and its receptor. CGRP potently inhibited alarmin-driven type 2 cytokine production and proliferation by lung ILC2s both in vitro and in vivo. CGRP induced marked changes in ILC2 expression programs in vivo and in vitro, attenuating alarmin-driven proliferative and effector responses. A distinct subset of ILCs scored highly for a CGRP-specific gene signature after in vivo alarmin stimulation, suggesting CGRP regulated this response. Finally, we observed increased ILC2 proliferation and type 2 cytokine production as well as exaggerated responses to alarmins in mice lacking the CGRP receptor. Together, these data indicate that endogenous CGRP is a critical negative regulator of ILC2 responses in vivo.

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.

Editing the epigenome: technologies for programmable transcription and epigenetic modulation.

Gene regulation is a complex and tightly controlled process that defines cell identity, health and disease, and response to pharmacologic and environmental signals. Recently developed DNA-targeting platforms, including zinc finger proteins, transcription activator-like effectors (TALEs) and the clustered, regularly interspaced, short palindromic repeats (CRISPR)-Cas9 system, have enabled the recruitment of transcriptional modulators and epigenome-modifying factors to any genomic site, leading to new insights into the function of epigenetic marks in gene expression. Additionally, custom transcriptional and epigenetic regulation is facilitating refined control over cell function and decision making. The unique properties of the CRISPR-Cas9 system have created new opportunities for high-throughput genetic screens and multiplexing targets to manipulate complex gene expression patterns. This Review summarizes recent technological developments in this area and their application to biomedical challenges. We also discuss remaining limitations and necessary future directions for this field.

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.

Design, Assembly, and Characterization of TALE-Based Transcriptional Activators and Repressors.

Transcription activator-like effectors (TALEs) are modular DNA-binding proteins that can be fused to a variety of effector domains to regulate the epigenome. Nucleotide recognition by TALE monomers follows a simple cipher, making this a powerful and versatile method to activate or repress gene expression. Described here are methods to design, assemble, and test TALE transcription factors (TALE-TFs) for control of endogenous gene expression. In this protocol, TALE arrays are constructed by Golden Gate cloning and tested for activity by transfection and quantitative RT-PCR. These methods for engineering TALE-TFs are useful for studies in reverse genetics and genomics, synthetic biology, and gene therapy.