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.

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.

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.

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

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.

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.

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.

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.

Highly specific epigenome editing by CRISPR-Cas9 repressors for silencing of distal regulatory elements.

Epigenome editing with the CRISPR (clustered, regularly interspaced, short palindromic repeats)-Cas9 platform is a promising technology for modulating gene expression to direct cell phenotype and to dissect the causal epigenetic mechanisms of gene regulation. Fusions of nuclease-inactive dCas9 to the Krüppel-associated box (KRAB) repressor (dCas9-KRAB) can silence target gene expression, but the genome-wide specificity and the extent of heterochromatin formation catalyzed by dCas9-KRAB are not known. We targeted dCas9-KRAB to the HS2 enhancer, a distal regulatory element that orchestrates the expression of multiple globin genes, and observed highly specific induction of H3K9 trimethylation (H3K9me3) at the enhancer and decreased chromatin accessibility of both the enhancer and its promoter targets. Targeted epigenetic modification of HS2 silenced the expression of multiple globin genes, with minimal off-target changes in global gene expression. These results demonstrate that repression mediated by dCas9-KRAB is sufficiently specific to disrupt the activity of individual enhancers via local modification of the epigenome.

RNA-guided gene activation by CRISPR-Cas9-based transcription factors.

Technologies for engineering synthetic transcription factors have enabled many advances in medical and scientific research. In contrast to existing methods based on engineering of DNA-binding proteins, we created a Cas9-based transactivator that is targeted to DNA sequences by guide RNA molecules. Coexpression of this transactivator and combinations of guide RNAs in human cells induced specific expression of endogenous target genes, demonstrating a simple and versatile approach for RNA-guided gene activation.

Correction of dystrophin expression in cells from Duchenne muscular dystrophy patients through genomic excision of exon 51 by zinc finger nucleases.

Duchenne muscular dystrophy (DMD) is caused by genetic mutations that result in the absence of dystrophin protein expression. Oligonucleotide-induced exon skipping can restore the dystrophin reading frame and protein production. However, this requires continuous drug administration and may not generate complete skipping of the targeted exon. In this study, we apply genome editing with zinc finger nucleases (ZFNs) to permanently remove essential splicing sequences in exon 51 of the dystrophin gene and thereby exclude exon 51 from the resulting dystrophin transcript. This approach can restore the dystrophin reading frame in ~13% of DMD patient mutations. Transfection of two ZFNs targeted to sites flanking the exon 51 splice acceptor into DMD patient myoblasts led to deletion of this genomic sequence. A clonal population was isolated with this deletion and following differentiation we confirmed loss of exon 51 from the dystrophin mRNA transcript and restoration of dystrophin protein expression. Furthermore, transplantation of corrected cells into immunodeficient mice resulted in human dystrophin expression localized to the sarcolemmal membrane. Finally, we quantified ZFN toxicity in human cells and mutagenesis at predicted off-target sites. This study demonstrates a powerful method to restore the dystrophin reading frame and protein expression by permanently deleting exons.

Reading frame correction by targeted genome editing restores dystrophin expression in cells from Duchenne muscular dystrophy patients.

Genome editing with engineered nucleases has recently emerged as an approach to correct genetic mutations by enhancing homologous recombination with a DNA repair template. However, many genetic diseases, such as Duchenne muscular dystrophy (DMD), can be treated simply by correcting a disrupted reading frame. We show that genome editing with transcription activator-like effector nucleases (TALENs), without a repair template, can efficiently correct the reading frame and restore the expression of a functional dystrophin protein that is mutated in DMD. TALENs were engineered to mediate highly efficient gene editing at exon 51 of the dystrophin gene. This led to restoration of dystrophin protein expression in cells from Duchenne patients, including skeletal myoblasts and dermal fibroblasts that were reprogrammed to the myogenic lineage by MyoD. Finally, exome sequencing of cells with targeted modifications of the dystrophin locus showed no TALEN-mediated off-target changes to the protein-coding regions of the genome, as predicted by in silico target site analysis. This strategy integrates the rapid and robust assembly of active TALENs with an efficient gene-editing method for the correction of genetic diseases caused by mutations in non-essential coding regions that cause frameshifts or premature stop codons.

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.

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.

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.

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.

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.

Environmental regulation of notochordal gene expression in nucleus pulposus cells.

Cells of the nucleus pulposus (NP) in the intervertebral disc are derived directly from the embryonic notochord. In humans, a shift in NP cell population coincides with the beginning of age-related changes in the extracellular matrix that can lead to spinal disorders. To begin identifying the bases of these changes, the manner by which relevant environmental factors impact cell function must be understood. This study investigated the roles of biochemical, nutritional, and physical factors in regulating immature NP cells. Specifically, we examined cell morphology, attachment, proliferation, and expression of genes associated with the notochord and immature NP (Sox9, CD24, and type IIA procollagen). Primary cells isolated from rat caudal discs were exposed to different media formulations and physical culture configurations either in 21% (ambient) or 2% (hypoxic) O2. As expected, cells in alginate beads retained a vacuolated morphology similar to chordocytes, with little change in gene expression. Interestingly, NP tissues not enzymatically digested were more profoundly influenced by oxygen. In monolayer, alpha-MEM preserved vacuolated morphology, produced the highest efficiency of attachment, and best maintained gene expression. DMEM and Opti-MEM cultures resulted in high levels of proliferation, but these appeared to involve small non-vacuolated cells. Gene expression patterns for cells in DMEM monolayer cultures were consistent with chondrocyte de-differentiation, with the response being delayed by hypoxia. Overall, results indicate that certain environmental conditions induce cellular changes that compromise the notochordal phenotype in immature NP. These results form the foundation on which the mechanisms of such changes can be elucidated.