Non-coding Regulatory Elements Research Articles

Locus of (IL-9) control: IL9 epigenetic regulation in cellular function and human disease.

Interleukin-9 (IL-9) is a multifunctional cytokine with roles in a broad cross-section of human diseases. Like many cytokines, IL-9 is transcriptionally regulated by a group of noncoding regulatory elements (REs) surrounding the IL9 gene. These REs modulate IL-9 transcription by forming 3D loops that recruit transcriptional machinery. IL-9-promoting transcription factors (TFs) can bind REs to increase locus accessibility and permit chromatin looping, or they can be recruited to already accessible chromatin to promote transcription. Ample mechanistic and genome-wide association studies implicate this interplay between IL-9-modulating TFs and IL9 cis-REs in human physiology, homeostasis, and disease.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs.

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Engineered CRISPR-Cas12a for higher-order combinatorial chromatin perturbations.

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting one to three genomic sites per cell. We engineer an Acidaminococcus Cas12a (AsCas12a) variant, multiplexed transcriptional interference AsCas12a (multiAsCas12a), that incorporates R1226A, a mutation that stabilizes the ribonucleoprotein-DNA complex via DNA nicking. The multiAsCas12a-KRAB fusion improves CRISPRi activity over DNase-dead AsCas12a-KRAB fusions, often rescuing the activities of lentivirally delivered CRISPR RNAs (crRNA) that are inactive when used with the latter. multiAsCas12a-KRAB supports CRISPRi using 6-plex crRNA arrays in high-throughput pooled screens. Using multiAsCas12a-KRAB, we discover enhancer elements and dissect the combinatorial function of cis-regulatory elements in human cells. These results instantiate a group testing framework for efficiently surveying numerous combinations of chromatin perturbations for biological discovery and engineering.

Deep learning the cis-regulatory code for gene expression in selected model plants

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.

Coronary Artery Disease Risk Variant Dampens the Expression of CALCRL by Reducing HSF Binding to Shear Stress Responsive Enhancer in Endothelial Cells In Vitro.

CALCRL (calcitonin receptor-like) protein is an important mediator of the endothelial fluid shear stress response, which is associated with the genetic risk of coronary artery disease. In this study, we functionally characterized the noncoding regulatory elements carrying coronary artery disease that risks single-nucleotide polymorphisms and studied their role in the regulation of CALCRL expression in endothelial cells. To functionally characterize the coronary artery disease single-nucleotide polymorphisms harbored around the gene CALCRL, we applied an integrative approach encompassing statistical, transcriptional (RNA-seq), and epigenetic (ATAC-seq [transposase-accessible chromatin with sequencing], chromatin immunoprecipitation assay-quantitative polymerase chain reaction, and electromobility shift assay) analyses, alongside luciferase reporter assays, and targeted gene and enhancer perturbations (siRNA and clustered regularly interspaced short palindromic repeats/clustered regularly interspaced short palindromic repeat-associated 9) in human aortic endothelial cells. We demonstrate that the regulatory element harboring rs880890 exhibits high enhancer activity and shows significant allelic bias. The A allele was favored over the G allele, particularly under shear stress conditions, mediated through alterations in the HSF1 (heat shock factor 1) motif and binding. CRISPR deletion of rs880890 enhancer resulted in downregulation of CALCRL expression, whereas HSF1 knockdown resulted in a significant decrease in rs880890-enhancer activity and CALCRL expression. A significant decrease in HSF1 binding to the enhancer region in endothelial cells was observed under disturbed flow compared with unidirectional flow. CALCRL knockdown and variant perturbation experiments indicated the role of CALCRL in mediating eNOS (endothelial nitric oxide synthase), APLN (apelin), angiopoietin, prostaglandins, and EDN1 (endothelin-1) signaling pathways leading to a decrease in cell proliferation, tube formation, and NO production. Overall, our results demonstrate the existence of an endothelial-specific HSF (heat shock factor)-regulated transcriptional enhancer that mediates CALCRL expression. A better understanding of CALCRL gene regulation and the role of single-nucleotide polymorphisms in the modulation of CALCRL expression could provide important steps toward understanding the genetic regulation of shear stress signaling responses.

Chromosome 20p11.2 deletions cause congenital hyperinsulinism via the loss of FOXA2 or its regulatory elements.

Persistent congenital hyperinsulinism (HI) is a rare genetically heterogeneous condition characterised by dysregulated insulin secretion leading to life-threatening hypoglycaemia. For up to 50% of affected individuals screening of the known HI genes does not identify a disease-causing variant. Large deletions have previously been used to identify novel regulatory regions causing HI. Here, we used genome sequencing to search for novel large (>1 Mb) deletions in 180 probands with HI of unknown cause and replicated our findings in a large cohort of 883 genetically unsolved individuals with HI using off-target copy number variant calling from targeted gene panels. We identified overlapping heterozygous deletions in five individuals (range 3-8 Mb) spanning chromosome 20p11.2. The pancreatic beta-cell transcription factor gene, FOXA2, a known cause of HI was deleted in two of the five individuals. In the remaining three, we found a minimal deleted region of 2.4 Mb adjacent to FOXA2 that encompasses multiple non-coding regulatory elements that are in conformational contact with FOXA2. Our data suggests that the deletions in these three children may cause disease through the dysregulation of FOXA2 expression. These findings provide new insights into the regulation of FOXA2 in the beta-cell and confirm an aetiological role for chromosome 20p11.2 deletions in syndromic HI.

Regional random mutagenesis driven by multiple sgRNAs and diverse on-target genome editing events to identify functionally important elements in non-coding regions.

Functional regions that regulate biological phenomena are interspersed throughout eukaryotic genomes. The most definitive approach for identifying such regions is to confirm the phenotype of cells or organisms in which specific regions have been mutated or removed from the genome. This approach is invaluable for the functional analysis of genes with a defined functional element, the protein-coding sequence. By contrast, no functional analysis platforms have been established for the study of cis-elements or microRNA cluster regions consisting of multiple microRNAs with functional overlap. Whole-genome mutagenesis approaches, such as via N-ethyl-N-nitrosourea and gene trapping, have greatly contributed to elucidating the function of coding genes. These methods almost never induce deletions of genomic regions or multiple mutations within a narrow region. In other words, cis-elements and microRNA clusters cannot be effectively targeted in such a manner. Herein, we established a novel region-specific random mutagenesis method named CRISPR- and transposase-based regional mutagenesis (CTRL-mutagenesis). We demonstrate that CTRL-mutagenesis randomly induces diverse mutations within target regions in murine embryonic stem cells. Comparative analysis of mutants harbouring subtly different mutations within the same region would facilitate the further study of cis-element and microRNA clusters.

The level of protein in the maternal murine diet modulates the facial appearance of the offspring via mTORC1 signaling

The development of craniofacial skeletal structures is fascinatingly complex and elucidation of the underlying mechanisms will not only provide novel scientific insights, but also help develop more effective clinical approaches to the treatment and/or prevention of the numerous congenital craniofacial malformations. To this end, we performed a genome-wide analysis of RNA transcription from non-coding regulatory elements by CAGE-sequencing of the facial mesenchyme of human embryos and cross-checked the active enhancers thus identified against genes, identified by GWAS for the normal range human facial appearance. Among the identified active cis-enhancers, several belonged to the components of the PI3/AKT/mTORC1/autophagy pathway. To assess the functional role of this pathway, we manipulated it both genetically and pharmacologically in mice and zebrafish. These experiments revealed that mTORC1 signaling modulates craniofacial shaping at the stage of skeletal mesenchymal condensations, with subsequent fine-tuning during clonal intercalation. This ability of mTORC1 pathway to modulate facial shaping, along with its evolutionary conservation and ability to sense external stimuli, in particular dietary amino acids, indicate that the mTORC1 pathway may play a role in facial phenotypic plasticity. Indeed, the level of protein in the diet of pregnant female mice influenced the activity of mTORC1 in fetal craniofacial structures and altered the size of skeletogenic clones, thus exerting an impact on the local geometry and craniofacial shaping. Overall, our findings indicate that the mTORC1 signaling pathway is involved in the effect of environmental conditions on the shaping of craniofacial structures.

Mapping the functional impact of non-coding regulatory elements in primary T cells through single-cell CRISPR screens

BackgroundDrug targets with genetic evidence are expected to increase clinical success by at least twofold. Yet, translating disease-associated genetic variants into functional knowledge remains a fundamental challenge of drug discovery. A key issue is that the vast majority of complex disease associations cannot be cleanly mapped to a gene. Immune disease-associated variants are enriched within regulatory elements found in T-cell-specific open chromatin regions.ResultsTo identify genes and molecular programs modulated by these regulatory elements, we develop a CRISPRi-based single-cell functional screening approach in primary human T cells. Our pipeline enables the interrogation of transcriptomic changes induced by the perturbation of regulatory elements at scale. We first optimize an efficient CRISPRi protocol in primary CD4+ T cells via CROPseq vectors. Subsequently, we perform a screen targeting 45 non-coding regulatory elements and 35 transcription start sites and profile approximately 250,000 T -cell single-cell transcriptomes. We develop a bespoke analytical pipeline for element-to-gene (E2G) mapping and demonstrate that our method can identify both previously annotated and novel E2G links. Lastly, we integrate genetic association data for immune-related traits and demonstrate how our platform can aid in the identification of effector genes for GWAS loci.ConclusionsWe describe “primary T cell crisprQTL” — a scalable, single-cell functional genomics approach for mapping regulatory elements to genes in primary human T cells. We show how this framework can facilitate the interrogation of immune disease GWAS hits and propose that the combination of experimental and QTL-based techniques is likely to address the variant-to-function problem.

Optimizing sequence design strategies for perturbation MPRAs: a computational evaluation framework.

The advent of perturbation-based massively parallel reporter assays (MPRAs) technique has facilitated the delineation of the roles of non-coding regulatory elements in orchestrating gene expression. However, computational efforts remain scant to evaluate and establish guidelines for sequence design strategies for perturbation MPRAs. In this study, we propose a framework for evaluating and comparing various perturbation strategies for MPRA experiments. Within this framework, we benchmark three different perturbation approaches from the perspectives of alteration in motif-based profiles, consistency of MPRA outputs, and robustness of models that predict the activities of putative regulatory motifs. While our analyses show very similar results across multiple benchmarking metrics, the predictive modeling for the approach involving random nucleotide shuffling shows significant robustness compared with the other two approaches. Thus, we recommend designing sequences by randomly shuffling the nucleotides of the perturbed site in perturbation-MPRA, followed by a coherence check to prevent the introduction of other variations of the target motifs. In summary, our evaluation framework and the benchmarking findings create a resource of computational pipelines and highlight the potential of perturbation-MPRA in predicting non-coding regulatory activities.

Heritable, multinucleotide deletions in plants using viral delivery of a repair exonuclease and guide RNAs.

CRISPR/Cas9-mediated mutagenesis typically results in short insertion/deletion mutations, which are often too small to disrupt the function of cis-acting regulatory elements. Here, we describe a highly efficient in planta gene editing approach called VirTREX2-HLDel that achieves heritable multinucleotide deletions in both protein-coding genes and noncoding DNA regulatory elements. VirTREX2-HLDel uses RNA viruses to deliver both the 3 prime repair exonuclease 2 (TREX2) and single-guide RNAs. Our method enables recovery of multiplexed heritable deletions and increases the heritable gene editing frequency at poorly edited sites. We identified functional conservation and divergence of MICRORNA164 (miR164) in Nicotiana benthamiana and tomato (Solanum lycopersicum) using VirTREX2-HLDel and observed previously uncharacterized phenotypes in plants with large deletions at this locus. Our viral delivery method reduces the need for tissue culture and will accelerate the understanding of protein-coding and regulatory regions in plants.

Role of Genetic Variation in Transcriptional Regulatory Elements in Heart Rhythm.

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

Prognostic analysis of mutated genes and insight into effects of DNA damage and repair on mutational strand asymmetries in gastric cancer

Gastric cancer (GACA) is a complex and multifaceted disease influenced by a variety of environmental and genetic factors. Somatic mutations play a major role in its development, and their characteristics, including the asymmetry between two DNA strands, are of great interest and appear as a signal of information and guidance, revealing mechanisms of DNA damage and repair. Here, we analyzed the impact of High-frequency mutated genes on patient prognosis and found that the effect of expression levels of tumor protein p53 (TP53) and lysine methyltransferase 2C (KMT2C) genes remained high throughout the development of GACA, with similar expression patterns. After investigating mutation asymmetry across mutagenic processes, we found that transcriptional asymmetry was dominated by T > G mutations under the influence of transcription couples repair and damage. The apolipoprotein B mRNA editing enzyme catalytic polypeptide like (APOBEC) enzyme that induces mutations during DNA replication has been identified here and we identified a replicative asymmetry, which was dominated by C > A mutations in left-replicating. Strand bias in different mutation classes at transcription factor binding sites and enhancer regions were also confirmed, which implies the important role of non-coding regulatory elements in the occurrence of mutations. This work systematically describes mutational strand asymmetries in specific genomic regions, shedding light on the DNA damage and repair mechanisms underlying somatic mutations in cohorts of GACA patients with gastric cancer.

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages.

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

AP-1 signaling modulates cardiac fibroblast stress responses.

Matrix remodeling outcomes largely dictate patient survival post myocardial infarction. Moreover, human-restricted noncoding regulatory elements have been shown to worsen fibrosis, but their mechanism of action remains elusive. Here, we demonstrate, using induced pluripotent stem cell-derived cardiac fibroblasts (iCFs), that inflammatory ligands abundant in the remodeling heart after infarction activate AP-1 transcription factor signaling pathways resulting in fibrotic responses. This observed signaling induces deposition of fibronectin matrix and is further capable of supporting immune cell adhesion; pathway inhibition blocks iCF matrix production and cell adhesion. Polymorphisms in the noncoding regulatory elements within the 9p21 locus (also referred to as ANRIL) redirect stress programs, and in iCFs, they transcriptionally silence the AP-1 inducible transcription factor GATA5. The presence of these polymorphisms modulate iCF matrix production and assembly and reduce cell-cell signaling. These data suggest that this signaling axis is a critical modulator of cardiac disease models and might be influenced by noncoding regulatory elements.

Cooperative Response to Endocardial Notch Reveals Interaction With Hippo Pathway.

The endocardium is a crucial signaling center for cardiac valve development and maturation. Genetic analysis has identified several human endocardial genes whose inactivation leads to bicuspid aortic valve formation and calcific aortic valve disease, but knowledge is very limited about the role played in valve development and disease by noncoding endocardial regulatory regions and upstream factors. We manipulated Notch signaling in mouse embryonic endocardial cells by short-term and long-term coculture with OP9 stromal cells expressing Notch ligands and inhibition of Notch activity. We examined the transcriptional profile and chromatin accessibility landscape for each condition, integrated transcriptomic, transcription factor occupancy, chromatin accessibility, and proteomic datasets. We generated in vitro and in vivo models with CRISPR-Cas9-edited deletions of various noncoding regulatory elements and validated their regulatory potential. We identified primary and secondary transcriptional responses to Notch ligands in the mouse embryonic endocardium, and a NOTCH-dependent transcriptional signature in valve development and disease. By defining the changes in the chromatin accessibility landscape and integrating with the landscape in developing mouse endocardium and adult human valves, we identify potential noncoding regulatory elements, validated selected candidates, propose interacting cofactors, and define the timeframe of their regulatory activity. Additionally, we found cooperative transcriptional repression with Hippo pathway by inhibiting nuclear Yap (Yes-associated protein) activity in the endocardium during cardiac valve development. Sequential Notch-dependent transcriptional regulation in the embryonic endocardium involves multiple factors. Notch activates certain noncoding elements through these factors and simultaneously suppresses elements that could hinder cardiac valve development and homeostasis. Biorxviv: https://www.biorxiv.org/content/10.1101/2023.03.23.533882v1.full.

Characterization of transcriptional enhancers in the chicken genome using CRISPR-mediated activation

DNA regulatory elements intricately control when, where, and how genes are activated. Therefore, understanding the function of these elements could unveil the complexity of the genetic regulation network. Genome-wide significant variants are predominantly found in non-coding regions of DNA, so comprehending the predicted functional regulatory elements is crucial for understanding the biological context of these genomic markers, which can be incorporated into breeding programs. The emergence of CRISPR technology has provided a powerful tool for studying non-coding regulatory elements in genomes. In this study, we leveraged epigenetic data from the Functional Annotation of Animal Genomes project to identify promoter and putative enhancer regions associated with three genes (HBBA, IRF7, and PPARG) in the chicken genome. To identify the enhancer regions, we designed guide RNAs targeting the promoter and candidate enhancer regions and utilized CRISPR activation (CRISPRa) with dCas9-p300 and dCas9-VPR as transcriptional activators in chicken DF-1 cells. By comparing the expression levels of target genes between the promoter activation and the co-activation of the promoter and putative enhancers, we were able to identify functional enhancers that exhibited augmented upregulation. In conclusion, our findings demonstrate the remarkable efficiency of CRISPRa in precisely manipulating the expression of endogenous genes by targeting regulatory elements in the chicken genome, highlighting its potential for functional validation of non-coding regions.

Beyond genome-wide association studies: Investigating the role of noncoding regulatory elements in primary sclerosing cholangitis.

Genome-wide association studies (GWAS) have identified 30 risk loci for primary sclerosing cholangitis (PSC). Variants within these loci are found predominantly in noncoding regions of DNA making their mechanisms of conferring risk hard to define. Epigenomic studies have shown noncoding variants broadly impact regulatory element activity. The possible association of noncoding PSC variants with regulatory element activity has not been studied. We aimed to (1) determine if the noncoding risk variants in PSC impact regulatory element function and (2) if so, assess the role these regulatory elements have in explaining the genetic risk for PSC. Available epigenomic datasets were integrated to build a comprehensive atlas of cell type-specific regulatory elements, emphasizing PSC-relevant cell types. RNA-seq and ATAC-seq were performed on peripheral CD4+ T cells from 10 PSC patients and 11 healthy controls. Computational techniques were used to (1) study the enrichment of PSC-risk variants within regulatory elements, (2) correlate risk genotype with differences in regulatory element activity, and (3) identify regulatory elements differentially active and genes differentially expressed between PSC patients and controls. Noncoding PSC-risk variants are strongly enriched within immune-specific enhancers, particularly ones involved in T-cell response to antigenic stimulation. In total, 250 genes and >10,000 regulatory elements were identified that are differentially active between patients and controls. Mechanistic effects are proposed for variants at 6 PSC-risk loci where genotype was linked with differential T-cell regulatory element activity. Regulatory elements are shown to play a key role in PSC pathophysiology.

Ribonucleic Acid-Mediated Control of Protein Translation Under Stress.

Significance: The need of cells to constantly respond to endogenous and exogenous stress has necessitated the evolution of pathways to counter the deleterious effects of stress and to restore cellular homeostasis. The inability to activate a timely and adequate response can lead to disease and is a hallmark of aging. Besides protein-coding genes, cells contain a plethora of noncoding regulatory elements that allow cells to respond rapidly and efficiently to external stimuli by activating highly specific and tightly controlled mechanisms. Many of these programs converge on the regulation of translation, one of the most energy-consuming processes in cells. Recent Advances: The noncoding dimension of translational regulation includes short and long noncoding ribonucleic acids (ncRNAs), as well as messenger RNA features, such as the sequence and modification status of the 5' and 3' untranslated regions (UTRs), that do not change the amino acid sequence of the produced protein. Critical Issues: In this review, we discuss the regulatory role of the nonprotein-coding components of translation under stress, particularly oxidative stress. We conclude that the regulation of translation through ncRNAs, UTRs, and nucleotide modifications is emerging as a critical component of the stress response. Future Directions: Further areas of study using long-read sequencing technologies will be discussed. Antioxid. Redox Signal. 39, 374-389.

Quantitative control of Ets1 dosage by a multi-enhancer hub promotes Th1 cell differentiation and protects from allergic inflammation

Multi-enhancer hubs are spatial clusters of enhancers present across numerous developmental programs. Here, we studied the functional relevance of these three-dimensional structures in Tcell biology. Mathematical modeling identified a highly connected multi-enhancer hub at the Ets1 locus, comprising a noncoding regulatory element that was a hotspot for sequence variation associated with allergic disease in humans. Deletion of this regulatory element in mice revealed that the multi-enhancer connectivity was dispensable for Tcell development but required for CD4+ T helper 1 (Th1) differentiation. These mice were protected from Th1-mediated colitis but exhibited overt allergic responses. Mechanistically, the multi-enhancer hub controlled the dosage of Ets1 that was required for CTCF recruitment and assembly of Th1-specific genome topology. Our findings establish a paradigm wherein multi-enhancer hubs control cellular competence to respond to an inductive cue through quantitative control of gene dosage and provide insight into how sequence variation within noncoding elements at the Ets1 locus predisposes individuals to allergic responses.