A model for accurate quantification of CRISPR effects in pooled FACS screens.

CRISPR screens are powerful tools to identify key genes that underlie biological processes. One important type of screen uses fluorescence activated cell sorting (FACS) to sort perturbed cells into bins based on the expression level of marker genes, followed by guide RNA (gRNA) sequencing. Analysis of these data presents several statistical challenges due to multiple factors including the discrete nature of the bins and typically small numbers of replicate experiments. To address these challenges, we developed a robust and powerful Bayesian random effects model and software package called Waterbear. Furthermore, we used Waterbear to explore how various experimental design parameters affect statistical power to establish principled guidelines for future screens. Finally, we experimentally validated our experimental design model findings that, when using Waterbear for analysis, high power is maintained even at low cell coverage and a high multiplicity of infection. We anticipate that Waterbear will be of broad utility for analyzing FACS-based CRISPR screens.

Long noncoding RNA LIRIL2R modulates FOXP3 levels and suppressive function of human CD4+ regulatory T cells by regulating IL2RA.

Regulatory T cells (Tregs) are central in controlling immune responses, and dysregulation of their function can lead to autoimmune disorders or cancer. Despite extensive studies on Tregs, the basis of epigenetic regulation of human Treg development and function is incompletely understood. Long intergenic noncoding RNAs (lincRNA)s are important for shaping and maintaining the epigenetic landscape in different cell types. In this study, we identified a gene on the chromosome 6p25.3 locus, encoding a lincRNA, that was up-regulated during early differentiation of human Tregs. The lincRNA regulated the expression of interleukin-2 receptor alpha (IL2RA), and we named it the lincRNA regulator of IL2RA (LIRIL2R). Through transcriptomics, epigenomics, and proteomics analysis of LIRIL2R-deficient Tregs, coupled with global profiling of LIRIL2R binding sites using chromatin isolation by RNA purification, followed by sequencing, we identified IL2RA as a target of LIRIL2R. This nuclear lincRNA binds upstream of the IL2RA locus and regulates its epigenetic landscape and transcription. CRISPR-mediated deletion of the LIRIL2R-bound region at the IL2RA locus resulted in reduced IL2RA expression. Notably, LIRIL2R deficiency led to reduced expression of Treg-signature genes (e.g., FOXP3, CTLA4, and PDCD1), upregulation of genes associated with effector T cells (e.g., SATB1 and GATA3), and loss of Treg-mediated suppression.

Systematic decoding of cis gene regulation defines context-dependent control of the multi-gene costimulatory receptor locus in human T cells.

Cis-regulatory elements (CREs) interact with trans regulators to orchestrate gene expression, but how transcriptional regulation is coordinated in multi-gene loci has not been experimentally defined. We sought to characterize the CREs controlling dynamic expression of the adjacent costimulatory genes CD28, CTLA4 and ICOS, encoding regulators of T cell-mediated immunity. Tiling CRISPR interference (CRISPRi) screens in primary human T cells, both conventional and regulatory subsets, uncovered gene-, cell subset- and stimulation-specific CREs. Integration with CRISPR knockout screens and assay for transposase-accessible chromatin with sequencing (ATAC-seq) profiling identified trans regulators influencing chromatin states at specific CRISPRi-responsive elements to control costimulatory gene expression. We then discovered a critical CCCTC-binding factor (CTCF) boundary that reinforces CRE interaction with CTLA4 while also preventing promiscuous activation of CD28. By systematically mapping CREs and associated trans regulators directly in primary human T cell subsets, this work overcomes longstanding experimental limitations to decode context-dependent gene regulatory programs in a complex, multi-gene locus critical to immune homeostasis.

A proximal enhancer regulates RORA expression during early human Th17 cell differentiation.

Gene regulatory elements, such as enhancers, greatly influence cell identity by tuning the transcriptional activity of specific cell types. Dynamics of enhancer landscape during early human Th17 cell differentiation remains incompletely understood. Leveraging ATAC-seq-based profiling of chromatin accessibility and comprehensive analysis of key histone marks, we identified a repertoire of enhancers that potentially exert control over the fate specification of Th17 cells. We found 23 SNPs associated with autoimmune diseases within Th17-enhancers that precisely overlapped with the binding sites of transcription factors actively engaged in T-cell functions. Among the Th17-specific enhancers, we identified an enhancer in the intron of RORA and demonstrated that this enhancer positively regulates RORA transcription. Moreover, CRISPR-Cas9-mediated deletion of a transcription factor binding site-rich region within the identified RORA enhancer confirmed its role in regulating RORA transcription. These findings provide insights into the potential mechanism by which the RORA enhancer orchestrates Th17 differentiation.

Bidirectional epigenetic editing reveals hierarchies in gene regulation.

CRISPR perturbation methods are limited in their ability to study non-coding elements and genetic interactions. In this study, we developed a system for bidirectional epigenetic editing, called CRISPRai, in which we apply activating (CRISPRa) and repressive (CRISPRi) perturbations to two loci simultaneously in the same cell. We developed CRISPRai Perturb-seq by coupling dual perturbation gRNA detection with single-cell RNA sequencing, enabling study of pooled perturbations in a mixed single-cell population. We applied this platform to study the genetic interaction between two hematopoietic lineage transcription factors, SPI1 and GATA1, and discovered novel characteristics of their co-regulation on downstream target genes, including differences in SPI1 and GATA1 occupancy at genes that are regulated through different modes. We also studied the regulatory landscape of IL2 (interleukin-2) in Jurkat T cells, primary T cells and chimeric antigen receptor (CAR) T cells and elucidated mechanisms of enhancer-mediated IL2 gene regulation. CRISPRai facilitates investigation of context-specific genetic interactions, provides new insights into gene regulation and will enable exploration of non-coding disease-associated variants.

Improving prime editing with an endogenous small RNA-binding protein.

Prime editing enables the precise modification of genomes through reverse transcription of template sequences appended to the 3' ends of CRISPR-Cas guide RNAs1. To identify cellular determinants of prime editing, we developed scalable prime editing reporters and performed genome-scale CRISPR-interference screens. From these screens, a single factor emerged as the strongest mediator of prime editing: the small RNA-binding exonuclease protection factor La. Further investigation revealed that La promotes prime editing across approaches (PE2, PE3, PE4 and PE5), edit types (substitutions, insertions and deletions), endogenous loci and cell types but has no consistent effect on genome-editing approaches that rely on standard, unextended guide RNAs. Previous work has shown that La binds polyuridine tracts at the 3' ends of RNA polymerase III transcripts2. We found that La functionally interacts with the 3' ends of polyuridylated prime editing guide RNAs (pegRNAs). Guided by these results, we developed a prime editor protein (PE7) fused to the RNA-binding, N-terminal domain of La. This editor improved prime editing with expressed pegRNAs and engineered pegRNAs (epegRNAs), as well as with synthetic pegRNAs optimized for La binding. Together, our results provide key insights into how prime editing components interact with the cellular environment and suggest general strategies for stabilizing exogenous small RNAs therein.

CRISPR-Cas9 screen of E3 ubiquitin ligases identifies TRAF2 and UHRF1 as regulators of HIV latency in primary human T cells.

During HIV infection of CD4+ T cells, ubiquitin pathways are essential to viral replication and host innate immune response; however, the role of specific E3 ubiquitin ligases is not well understood. Proteomics analyses identified 116 single-subunit E3 ubiquitin ligases expressed in activated primary human CD4+ T cells. Using a CRISPR-based arrayed spreading infectivity assay, we systematically knocked out 116 E3s from activated primary CD4+ T cells and infected them with NL4-3 GFP reporter HIV-1. We found 10 E3s significantly positively or negatively affected HIV infection in activated primary CD4+ T cells, including UHRF1 (pro-viral) and TRAF2 (anti-viral). Furthermore, deletion of either TRAF2 or UHRF1 in three JLat models of latency spontaneously increased HIV transcription. To verify this effect, we developed a CRISPR-compatible resting primary human CD4+ T cell model of latency. Using this system, we found that deletion of TRAF2 or UHRF1 initiated latency reactivation and increased virus production from primary human resting CD4+ T cells, suggesting these two E3s represent promising targets for future HIV latency reversal strategies. IMPORTANCE: HIV, the virus that causes AIDS, heavily relies on the machinery of human cells to infect and replicate. Our study focuses on the host cell's ubiquitination system which is crucial for numerous cellular processes. Many pathogens, including HIV, exploit this system to enhance their own replication and survival. E3 proteins are part of the ubiquitination pathway that are useful drug targets for host-directed therapies. We interrogated the 116 E3s found in human immune cells known as CD4+ T cells, since these are the target cells infected by HIV. Using CRISPR, a gene-editing tool, we individually removed each of these enzymes and observed the impact on HIV infection in human CD4+ T cells isolated from healthy donors. We discovered that 10 of the E3 enzymes had a significant effect on HIV infection. Two of them, TRAF2 and UHRF1, modulated HIV activity within the cells and triggered an increased release of HIV from previously dormant or "latent" cells in a new primary T cell assay. This finding could guide strategies to perturb hidden HIV reservoirs, a major hurdle to curing HIV. Our study offers insights into HIV-host interactions, identifies new factors that influence HIV infection in immune cells, and introduces a novel methodology for studying HIV infection and latency in human immune cells.

Ultra-high efficiency T cell reprogramming at multiple loci with SEED-Selection.

Multiplexed reprogramming of T cell specificity and function can generate powerful next-generation cellular therapies. However, current manufacturing methods produce heterogenous mixtures of partially engineered cells. Here, we develop a one-step process to enrich for unlabeled cells with knock-ins at multiple target loci using a family of repair templates named Synthetic Exon/Expression Disruptors (SEEDs). SEED engineering associates transgene integration with the disruption of a paired endogenous surface protein, allowing non-modified and partially edited cells to be immunomagnetically depleted (SEED-Selection). We design SEEDs to fully reprogram three critical loci encoding T cell specificity, co-receptor expression, and MHC expression, with up to 98% purity after selection for individual modifications and up to 90% purity for six simultaneous edits (three knock-ins and three knockouts). These methods are simple, compatible with existing clinical manufacturing workflows, and can be readily adapted to other loci to facilitate production of complex gene-edited cell therapies.

Multiplexed CRISPR gene editing in primary human islet cells with Cas9 ribonucleoprotein.

Successful genome editing in primary human islets could reveal features of the genetic regulatory landscape underlying ß cell function and diabetes risk. Here, we describe a CRISPR-based strategy to interrogate functions of predicted regulatory DNA elements using electroporation of a complex of Cas9 ribonucleoprotein (Cas9 RNP) and guide RNAs into primary human islet cells. We successfully targeted coding regions including the PDX1 exon 1, and non-coding DNA linked to diabetes susceptibility. CRISPR-Cas9 RNP approaches revealed genetic targets of regulation by DNA elements containing candidate diabetes risk SNPs, including an in vivo enhancer of the MPHOSPH9 gene. CRISPR-Cas9 RNP multiplexed targeting of two cis-regulatory elements linked to diabetes risk in PCSK1, which encodes an endoprotease crucial for Insulin processing, also demonstrated efficient simultaneous editing of PCSK1 regulatory elements, resulting in impaired ß cell PCSK1 regulation and Insulin secretion. Multiplex CRISPR-Cas9 RNP provides powerful approaches to investigate and elucidate human islet cell gene regulation in health and diabetes.

Base-editing mutagenesis maps alleles to tune human T cell functions.

CRISPR-enabled screening is a powerful tool for the discovery of genes that control T cell function and has nominated candidate targets for immunotherapies1-6. However, new approaches are required to probe specific nucleotide sequences within key genes. Systematic mutagenesis in primary human T cells could reveal alleles that tune specific phenotypes. DNA base editors are powerful tools for introducing targeted mutations with high efficiency7,8. Here we develop a large-scale base-editing mutagenesis platform with the goal of pinpointing nucleotides that encode amino acid residues that tune primary human T cell activation responses. We generated a library of around 117,000 single guide RNA molecules targeting base editors to protein-coding sites across 385 genes implicated in T cell function and systematically identified protein domains and specific amino acid residues that regulate T cell activation and cytokine production. We found a broad spectrum of alleles with variants encoding critical residues in proteins including PIK3CD, VAV1, LCP2, PLCG1 and DGKZ, including both gain-of-function and loss-of-function mutations. We validated the functional effects of many alleles and further demonstrated that base-editing hits could positively and negatively tune T cell cytotoxic function. Finally, higher-resolution screening using a base editor with relaxed protospacer-adjacent motif requirements9 (NG versus NGG) revealed specific structural domains and protein-protein interaction sites that can be targeted to tune T cell functions. Base-editing screens in primary immune cells thus provide biochemical insights with the potential to accelerate immunotherapy design.

PIM kinases regulate early human Th17 cell differentiation.

The serine/threonine-specific Moloney murine leukemia virus (PIM) kinase family (i.e., PIM1, PIM2, and PIM3) has been extensively studied in tumorigenesis. PIM kinases are downstream of several cytokine signaling pathways that drive immune-mediated diseases. Uncontrolled T helper 17 (Th17) cell activation has been associated with the pathogenesis of autoimmunity. However, the detailed molecular function of PIMs in human Th17 cell regulation has yet to be studied. In the present study, we comprehensively investigated how the three PIMs simultaneously alter transcriptional gene regulation during early human Th17 cell differentiation. By combining PIM triple knockdown with bulk and scRNA-seq approaches, we found that PIM deficiency promotes the early expression of key Th17-related genes while suppressing Th1-lineage genes. Further, PIMs modulate Th cell signaling, potentially via STAT1 and STAT3. Overall, our study highlights the inhibitory role of PIMs in human Th17 cell differentiation, thereby suggesting their association with autoimmune phenotypes.

The lncRNA Malat1 inhibits miR-15/16 to enhance cytotoxic T cell activation and memory cell formation.

Proper activation of cytotoxic T cells via the T cell receptor and the costimulatory receptor CD28 is essential for adaptive immunity against viruses, intracellular bacteria, and cancers. Through biochemical analysis of RNA:protein interactions, we uncovered a non-coding RNA circuit regulating activation and differentiation of cytotoxic T cells composed of the long non-coding RNA Malat1 (Metastasis Associated Lung Adenocarcinoma Transcript 1) and the microRNA family miR-15/16. miR-15/16 is a widely and highly expressed tumor suppressor miRNA family important for cell proliferation and survival. miR-15/16 play important roles in T cell responses to viral infection, including the regulation of antigen-specific T cell expansion and memory. Comparative Argonaute-2 high-throughput sequencing of crosslinking immunoprecipitation (AHC) combined with gene expression profiling in normal and miR-15/16-deficient mouse T cells revealed a large network of hundreds of direct miR-15/16 target mRNAs, many with functional relevance for T cell activation, survival and memory formation. Among these targets, Malat1 contained the largest absolute magnitude miR-15/16-dependent AHC peak. This binding site was among the strongest lncRNA:miRNA interactions detected in the T cell transcriptome. We used CRISPR targeting with homology directed repair to generate mice with a 5-nucleotide mutation in the miR-15/16-binding site in Malat1. This mutation interrupted Malat1:miR-15/16 interaction, and enhanced the repression of other miR-15/16 target genes, including CD28. Interrupting Malat1 interaction with miR-15/16 decreased cytotoxic T cell activation, including the expression of interleukin 2 (IL-2) and a broader CD28-responsive gene program. Accordingly, Malat1 mutation diminished memory cell persistence in mice following LCMV Armstrong and Listeria monocytogenes infection. This study marks a significant advance in the study of long non-coding RNAs in the immune system by ascribing cell-intrinsic, sequence-specific in vivo function to Malat1. These findings have implications for T cell-mediated autoimmune diseases, antiviral and anti-tumor immunity, as well as lung adenocarcinoma and other malignancies where Malat1 is overexpressed.

HIC1 interacts with FOXP3 multi protein complex: Novel pleiotropic mechanisms to regulate human regulatory T cell differentiation and function.

Transcriptional repressor, hypermethylated in cancer 1 (HIC1) participates in a range of important biological processes, such as tumor repression, immune suppression, embryonic development and epigenetic gene regulation. Further to these, we previously demonstrated that HIC1 provides a significant contribution to the function and development of regulatory T (Treg) cells. However, the mechanism by which it regulates these processes was not apparent. To address this question, we used affinity-purification mass spectrometry to characterize the HIC1 interactome in human Treg cells. Altogether 61 high-confidence interactors were identified, including IKZF3, which is a key transcription factor in the development of Treg cells. The biological processes associated with these interacting proteins include protein transport, mRNA processing, non-coding (ncRNA) transcription and RNA metabolism. The results revealed that HIC1 is part of a FOXP3-RUNX1-CBFB protein complex that regulates Treg signature genes thus improving our understanding of HIC1 function during early Treg cell differentiation.

Modular pooled discovery of synthetic knockin sequences to program durable cell therapies.

Chronic stimulation can cause T cell dysfunction and limit the efficacy of cellular immunotherapies. Improved methods are required to compare large numbers of synthetic knockin (KI) sequences to reprogram cell functions. Here, we developed modular pooled KI screening (ModPoKI), an adaptable platform for modular construction of DNA KI libraries using barcoded multicistronic adaptors. We built two ModPoKI libraries of 100 transcription factors (TFs) and 129 natural and synthetic surface receptors (SRs). Over 30 ModPoKI screens across human TCR- and CAR-T cells in diverse conditions identified a transcription factor AP4 (TFAP4) construct that enhanced fitness of chronically stimulated CAR-T cells and anti-cancer function in vitro and in vivo. ModPoKI's modularity allowed us to generate an ∼10,000-member library of TF combinations. Non-viral KI of a combined BATF-TFAP4 polycistronic construct enhanced fitness. Overexpressed BATF and TFAP4 co-occupy and regulate key gene targets to reprogram T cell function. ModPoKI facilitates the discovery of complex gene constructs to program cellular functions.

Multiplexed CRISPR gene editing in primary human islet cells with Cas9 ribonucleoprotein.

Successful genome editing in primary human islets could reveal features of the genetic regulatory landscape underlying ß cell function and diabetes risk. Here, we describe a CRISPR-based strategy to interrogate functions of predicted regulatory DNA elements using electroporation of a complex of Cas9 ribonucleoprotein (Cas9 RNP) and guide RNAs into primary human islet cells. We successfully targeted coding regions including the PDX1 exon 1, and non-coding DNA linked to diabetes susceptibility. CRISPR/Cas9 RNP approaches revealed genetic targets of regulation by DNA elements containing candidate diabetes risk SNPs, including an in vivo enhancer of the MPHOSPH9 gene. CRISPR/Cas9 RNP multiplexed targeting of two cis-regulatory elements linked to diabetes risk in PCSK1, which encodes an endoprotease crucial for insulin processing, also demonstrated efficient simultaneous editing of PCSK1 regulatory elements, resulting in impaired ß cell PCSK1 regulation and insulin secretion. Multiplex CRISPR/Cas9 RNP provides powerful approaches to investigate and elucidate human islet cell gene regulation in health and diabetes.

Gene regulatory network inference from CRISPR perturbations in primary CD4+ T cells elucidates the genomic basis of immune disease.

The effects of genetic variation on complex traits act mainly through changes in gene regulation. Although many genetic variants have been linked to target genes in cis, the trans-regulatory cascade mediating their effects remains largely uncharacterized. Mapping trans-regulators based on natural genetic variation, including eQTL mapping, has been challenging due to small effects. Experimental perturbation approaches offer a complementary and powerful approach to mapping trans-regulators. We used CRISPR knockouts of 84 genes in primary CD4+ T cells to perturb an immune cell gene network, targeting both inborn error of immunity (IEI) disease transcription factors (TFs) and background TFs matched in constraint and expression level, but without a known immune disease association. We developed a novel Bayesian structure learning method called Linear Latent Causal Bayes (LLCB) to estimate the gene regulatory network from perturbation data and observed 211 directed edges among the genes which could not be detected in existing CD4+ trans-eQTL data. We used LLCB to characterize the differences between the IEI and background TFs, finding that the gene groups were highly interconnected, but that IEI TFs were much more likely to regulate immune cell specific pathways and immune GWAS genes. We further characterized nine coherent gene programs based on downstream effects of the TFs and linked these modules to regulation of GWAS genes, finding that canonical JAK-STAT family members are regulated by KMT2A, a global epigenetic regulator. These analyses reveal the trans-regulatory cascade from upstream epigenetic regulator to intermediate TFs to downstream effector cytokines and elucidate the logic linking immune GWAS genes to key signaling pathways.

CRISPR screens decode cancer cell pathways that trigger γδ T cell detection.

γδ T cells are potent anticancer effectors with the potential to target tumours broadly, independent of patient-specific neoantigens or human leukocyte antigen background1-5. γδ T cells can sense conserved cell stress signals prevalent in transformed cells2,3, although the mechanisms behind the targeting of stressed target cells remain poorly characterized. Vγ9Vδ2 T cells-the most abundant subset of human γδ T cells4-recognize a protein complex containing butyrophilin 2A1 (BTN2A1) and BTN3A1 (refs. 6-8), a widely expressed cell surface protein that is activated by phosphoantigens abundantly produced by tumour cells. Here we combined genome-wide CRISPR screens in target cancer cells to identify pathways that regulate γδ T cell killing and BTN3A cell surface expression. The screens showed previously unappreciated multilayered regulation of BTN3A abundance on the cell surface and triggering of γδ T cells through transcription, post-translational modifications and membrane trafficking. In addition, diverse genetic perturbations and inhibitors disrupting metabolic pathways in the cancer cells, particularly ATP-producing processes, were found to alter BTN3A levels. This induction of both BTN3A and BTN2A1 during metabolic crises is dependent on AMP-activated protein kinase (AMPK). Finally, small-molecule activation of AMPK in a cell line model and in patient-derived tumour organoids led to increased expression of the BTN2A1-BTN3A complex and increased Vγ9Vδ2 T cell receptor-mediated killing. This AMPK-dependent mechanism of metabolic stress-induced ligand upregulation deepens our understanding of γδ T cell stress surveillance and suggests new avenues available to enhance γδ T cell anticancer activity.

The lncRNA Malat1 Inhibits miR-15/16 to Enhance Cytotoxic T Cell Activation and Memory Cell Formation.

Proper activation of cytotoxic T cells via the T cell receptor and the costimulatory receptor CD28 is essential for adaptive immunity against viruses, many intracellular bacteria and cancers. Through biochemical analysis of RNA:protein interactions, we uncovered a non-coding RNA circuit regulating activation and differentiation of cytotoxic T cells composed of the long non-coding RNA Malat1 (Metastasis Associated Lung Adenocarcinoma Transcript 1) and the microRNA family miR-15/16. miR-15/16 is a widely and highly expressed tumor suppressor miRNA family important for cell proliferation and survival. miR-15/16 also play important roles in T cell responses to viral infection, including the regulation of antigen-specific T cell expansion and T cell memory. Comparative Argonaute-2 high throughput sequencing of crosslinking immunoprecipitation (Ago2 HITS-CLIP, or AHC) combined with gene expression profiling in normal and miR-15/16-deficient T cells revealed a large network of several hundred direct miR-15/16 target mRNAs, many with functional relevance for T cell activation, survival and memory formation. Among these targets, the long non-coding RNA Malat1 contained the largest absolute magnitude miR-15/16-dependent AHC peak in T cells. This binding site was also among the strongest lncRNA:miRNA interactions detected in the T cell transcriptome. We used CRISPR targeting with homology directed repair to generate mice with a 5-nucleotide mutation in the miR-15/16 binding site in Malat1. This mutation interrupted Malat1:miR-15/16 interaction, and enhanced the repression of other miR-15/16 target genes, including CD28. Interrupting Malat1 interaction with miR-15/16 decreased cytotoxic T cell activation, including the expression of IL-2 and a broader CD28-responsive gene program. Accordingly, Malat1 mutation diminished memory cell persistence following LCMV Armstrong and Listeria monocytogenes infection. This study marks a significant advance in the study of long noncoding RNAs in the immune system by ascribing cell-intrinsic, sequence-specific in vivo function to Malat1. These findings have implications for T cell-mediated autoimmune diseases, antiviral and anti-tumor immunity, as well as lung adenocarcinoma and other malignancies where Malat1 is overexpressed.

Non-coding sequence variation reveals fragility within interleukin 2 feedback circuitry and shapes autoimmune disease risk.

Genetic variants associated with human autoimmune diseases commonly map to non-coding control regions, particularly enhancers that function selectively in immune cells and fine-tune gene expression within a relatively narrow range of values. How such modest, cell-type-selective changes can meaningfully shape organismal disease risk remains unclear. To explore this issue, we experimentally manipulated species-conserved enhancers within the disease-associated IL2RA locus and studied accompanying changes in the progression of autoimmunity. Perturbing distinct enhancers with restricted activity in conventional T cells (Tconvs) or regulatory T cells (Tregs)-two functionally antagonistic T cell subsets-caused only modest, cell-type-selective decreases in IL2ra expression parameters. However, these same perturbations had striking and opposing effects in vivo , completely preventing or severely accelerating disease in a murine model of type 1 diabetes. Quantitative tissue imaging and computational modelling revealed that each enhancer manipulation impinged on distinct IL-2-dependent feedback circuits. These imbalances altered the intracellular signaling and intercellular communication dynamics of activated Tregs and Tconvs, producing opposing spatial domains that amplified or constrained ongoing autoimmune responses. These findings demonstrate how subtle changes in gene regulation stemming from non-coding variation can propagate across biological scales due to non-linearities in intra- and intercellular feedback circuitry, dramatically shaping disease risk at the organismal level.