Probing the interactions of NK cell receptors with ligand expressed in trans and cis.

Certain receptors on natural killer (NK) cells, which are specific for MHC class I (MHC-I) molecules, do not only interact with ligand expressed on opposing cell membranes (in trans) but also interact with those on the same cell membrane (in cis). Cis interactions have been demonstrated for only a small number of cell surface receptors. However, this has not been tested systematically, raising the possibility that additional receptors may be able to bind ligand expressed in cis. Here we describe a number of approaches to evaluate trans and cis binding of the Ly49A NK cell receptor to its H-2D(d) ligand. These procedures should facilitate the investigation of cis/trans interactions of other receptor-ligand pairs and simplify the analysis of NK cell receptor variants.

The loss of GLUT2 expression in the pancreatic beta-cells of diabetic db/db mice is associated with an impaired DNA-binding activity of islet-specific trans-acting factors.

GLUT2 expression is reduced in the pancreatic beta-cells of several diabetic animals. The transcriptional control of the gene in beta-cells involves at least two islet-specific DNA-binding proteins, GTIIa and PDX-1, which also transactivates the insulin, somatostatin and glucokinase genes. In this report, we assessed the DNA-binding activities of GTIIa and PDX-1 to their respective cis-elements of the GLUT2 promoter using nuclear extracts prepared from pancreatic islets of 12 week old db/db diabetic mice. We show that the decreased GLUT2 mRNA expression correlates with a decrease of the GTIIa DNA-binding activity, whereas the PDX-1 binding activity is increased. In these diabetic animals, insulin mRNA expression remains normal. The adjunction of dexamethasone to isolated pancreatic islets, a treatment previously shown to decrease PDX-1 expression in the insulin-secreting HIT-T15 cells, has no effect on the GTIIa and PDX-1 DNA-binding activities. These data suggest that the decreased activity of GTIIa, in contrast to PDX-1, may be a major initial step in the development of the beta-cell dysfunction in this model of diabetes.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells.

BACKGROUND: Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian. RESULTS: We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types. CONCLUSIONS: The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Posttranscriptional repression of GacS/GacA-controlled genes by the RNA-binding protein RsmE acting together with RsmA in the biocontrol strain Pseudomonas fluorescens CHA0.

In the plant-beneficial soil bacterium Pseudomonas fluorescens CHA0, the production of biocontrol factors (antifungal secondary metabolites and exoenzymes) is controlled at a posttranscriptional level by the GacS/GacA signal transduction pathway involving RNA-binding protein RsmA as a key regulatory element. This protein is assumed to bind to the ribosome-binding site of target mRNAs and to block their translation. RsmA-mediated repression is relieved at the end of exponential growth by two GacS/GacA-controlled regulatory RNAs RsmY and RsmZ, which bind and sequester the RsmA protein. A gene (rsmE) encoding a 64-amino-acid RsmA homolog was identified and characterized in strain CHA0. Overexpression of rsmE strongly reduced the expression of target genes (hcnA, for a hydrogen cyanide synthase subunit; aprA, for the main exoprotease; and phlA, for a component of 2,4-diacetylphloroglucinol biosynthesis). Single null mutations in either rsmA or rsmE resulted in a slight increase in the expression of hcnA, aprA, and phlA. By contrast, an rsmA rsmE double mutation led to strongly increased and advanced expression of these target genes and completely suppressed a gacS mutation. Both the RsmE and RsmA levels increased with increasing cell population densities in strain CHA0; however, the amount of RsmA showed less variability during growth. Expression of rsmE was controlled positively by GacA and negatively by RsmA and RsmE. Mobility shift assays demonstrated specific binding of RsmE to RsmY and RsmZ RNAs. The transcription and stability of both regulatory RNAs were strongly reduced in the rsmA rsmE double mutant. In conclusion, RsmA and RsmE together account for maximal repression in the GacS/GacA cascade of strain CHA0.

Regulatory RNA as mediator in GacA/RsmA-dependent global control of exoproduct formation in Pseudomonas fluorescens CHA0.

In Pseudomonas fluorescens CHA0, an antagonist of root-pathogenic fungi, the GacS/GacA two-component system tightly controls the expression of antifungal secondary metabolites and exoenzymes at a posttranscriptional level, involving the RNA-binding protein and global regulator of secondary metabolism RsmA. This protein was purified from P. fluorescens, and RNA bound to it was converted to cDNA, which served as a probe to isolate the corresponding chromosomal locus, rsmZ. This gene encoded a regulatory RNA of 127 nucleotides and a truncated form lacking 35 nucleotides at the 3' end. Expression of rsmZ depended on GacA, increased with increasing population density, and was stimulated by the addition of a solvent-extractable extracellular signal produced by strain CHA0 at the end of exponential growth. This signal appeared to be unrelated to N-acyl-homoserine lactones. A conserved upstream element in the rsmZ promoter, but not the stress sigma factor RpoS, was involved in rsmZ expression. Overexpression of rsmZ effectively suppressed the negative effect of gacS and gacA mutations on target genes, i.e., hcnA (for hydrogen cyanide synthase) and aprA (for the major exoprotease). Mutational inactivation of rsmZ resulted in reduced expression of these target genes in the presence of added signal. Overexpression of rsmA had a similar, albeit stronger negative effect. These results support a model in which GacA upregulates the expression of regulatory RNAs, such as RsmZ of strain CHA0, in response to a bacterial signal. By a titration effect, RsmZ may then alleviate the repressing activity of RsmA on the expression of target mRNAs.

Genetic Variants of Serum Uric Acid and Gout: An Analysis of > 170,000 Individuals

Background/Purpose: Gout is a common and excruciatingly painful inflammatory arthritis caused by hyperuricemia. In addition to various lifestyle risk factors, a substantial genetic predisposition to gout has long been recognized. The Global Urate Genetics Consortium (GUGC) has aimed to comprehensively investigate the genetics of serum uric acid and gout using data from _ 140,000 individuals of European-ancestry, 8,340 individuals of Indian ancestry, 5,820 African-Americans, and 15,286 Japanese. Methods: We performed discovery GWAS meta-analyses of serum urate levels (n_110,347 individuals) followed by replication analyses (n_32,813 different individuals). Our gout analysis involved 3,151 cases and 68,350 controls, including 1,036 incident gout cases that met the American College of Rheumatology Criteria. We also examined the association of gout with fractional excretion of uric acid (n_6,799). A weighted genetic urate score was constructed based on the number of risk alleles across urate-associated loci, and their association with the risk of gout was evaluated. Furthermore, we examined implicated transcript expression in cis (expression quantitative trait loci databases) for potential insights into the gene underlying the association signal. Finally, in order to further identify urate-associated genomic regions, we performed functional network analyses that incorporated prior knowledge on molecular interactions in which the gene products of implicated genes operate. Results: We identified and replicated 28 genome-wide significant loci in association with serum urate (P 5_10_8), including all previously-reported loci as well as 18 novel genetic loci. Unlike the majority of previouslyidentified loci, none of the novel loci appeared to be obvious candidates for urate transport. Rather, they were mapped to genes that encode for purine production, transcription, or growth factors with broad downstream responses. Besides SLC2A9 and ABCG2, no additional regions contained SNPs that differed significantly (P _ 5_10_8) between sexes. Urateincreasing alleles were associated with an increased risk of gout for all loci. The urate genetic risk score (ranging from 10 to 45) was significantly associated with an increased odds of prevalent gout (OR per unit increase, 1.11; 95% CI, 1.09-1.14) and incident gout (OR, 1.10; 95% CI, 1.08-1.13). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. Detailed characterization of the loci revealed associations with transcript expression and the fractional excretion of urate. Network analyses implicated the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. Conclusion: The novel genetic candidates identified in this urate/gout consortium study, the largest to date, highlight the importance of metabolic control of urate production and urate excretion. The modulation by signaling processes that influence metabolic pathways such as glycolysis and the pentose phosphate pathway appear to be central mechanisms underpinned by the novel GWAS candidates. These findings may have implications for further research into urate-lowering drugs to treat and prevent gout.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Identification of Novel Craniofacial Regulatory Domains Located far Upstream of SOX9 and Disrupted in Pierre Robin Sequence.

Mutations in the coding sequence of SOX9 cause campomelic dysplasia (CD), a disorder of skeletal development associated with 46,XY disorders of sex development (DSDs). Translocations, deletions, and duplications within a ∼2 Mb region upstream of SOX9 can recapitulate the CD-DSD phenotype fully or partially, suggesting the existence of an unusually large cis-regulatory control region. Pierre Robin sequence (PRS) is a craniofacial disorder that is frequently an endophenotype of CD and a locus for isolated PRS at ∼1.2-1.5 Mb upstream of SOX9 has been previously reported. The craniofacial regulatory potential within this locus, and within the greater genomic domain surrounding SOX9, remains poorly defined. We report two novel deletions upstream of SOX9 in families with PRS, allowing refinement of the regions harboring candidate craniofacial regulatory elements. In parallel, ChIP-Seq for p300 binding sites in mouse craniofacial tissue led to the identification of several novel craniofacial enhancers at the SOX9 locus, which were validated in transgenic reporter mice and zebrafish. Notably, some of the functionally validated elements fall within the PRS deletions. These studies suggest that multiple noncoding elements contribute to the craniofacial regulation of SOX9 expression, and that their disruption results in PRS.

The ArgR regulatory protein, a helper to the anaerobic regulator ANR during transcriptional activation of the arcD promoter in Pseudomonas aeruginosa.

Pseudomonas aeruginosa, when deprived of oxygen, generates ATP from arginine catabolism by enzymes of the arginine deiminase pathway, encoded by the arcDABC operon. Under conditions of low oxygen tension, the transcriptional activator ANR binds to a site centered 41.5 bp upstream of the arcD transcriptional start. ANR-mediated anaerobic induction was enhanced two- to threefold by extracellular arginine. This arginine effect depended, in trans, on the transcriptional regulator ArgR and, in cis, on an ArgR binding site centered at -73.5 bp in the arcD promoter. Binding of purified ArgR protein to this site was demonstrated by electrophoretic mobility shift assays and DNase I footprinting. This ArgR recognition site contained a sequence, 5'-TGACGC-3', which deviated in only 1 base from the common sequence motif 5'-TGTCGC-3' found in other ArgR binding sites of P. aeruginosa. Furthermore, an alignment of all known ArgR binding sites confirmed that they consist of two directly repeated half-sites. In the absence of ANR, arginine did not induce the arc operon, suggesting that ArgR alone does not activate the arcD promoter. According to a model proposed, ArgR makes physical contact with ANR and thereby facilitates initiation of arc transcription.

Clonal acquisition of the Ly49A NK cell receptor is dependent on the trans-acting factor TCF-1.

Families of clonally expressed major histocompatibility complex (MHC) class I-specific receptors provide specificity to and regulate the function of natural killer (NK) cells. One of these receptors, mouse Ly49A, is expressed by 20% of NK cells and inhibits the killing of H-2D(d) but not D(b)-expressing target cells. Here, we show that the trans-acting factor TCF-1 binds to two sites in the Ly49A promoter and regulates its activity. Moreover, we find that TCF-1 determines the size of the Ly49A NK cell subset in vivo in a dosage-dependent manner. We propose that clonal Ly49A acquisition during NK cell development is regulated by TCF-1.

Structural variation and its effect on expression.

Structural variation, whether it is caused by copy number variants or present in a balanced form, such as reciprocal translocations and inversions, can have a profound and dramatic effect on the expression of genes mapping within and close to the rearrangement, as well as affecting others genome wide. These effects can be caused by altering the copy number of one or more genes or regulatory elements (dosage effect) or from physical disruption of links between regulatory elements and their associated gene or genes, resulting in perturbation of expression. Similarly, large-scale structural variants can result in genome-wide expression changes by altering the positions that chromosomes occupy within the nucleus, potentially disrupting not only local cis interactions, but also trans interactions that occur throughout the genome. Structural variation is, therefore, a significant factor in the study of gene expression and is discussed here in more detail.

Spatio-temporal sequence of cross-regulatory events in root meristem growth.

A central question in developmental biology is how multicellular organisms coordinate cell division and differentiation to determine organ size. In Arabidopsis roots, this balance is controlled by cytokinin-induced expression of SHORT HYPOCOTYL 2 (SHY2) in the so-called transition zone of the meristem, where SHY2 negatively regulates auxin response factors (ARFs) by protein-protein interaction. The resulting down-regulation of PIN-FORMED (PIN) auxin efflux carriers is considered the key event in promoting differentiation of meristematic cells. Here we show that this regulation involves additional, intermediary factors and is spatio-temporally constrained. We found that the described cytokinin-auxin crosstalk antagonizes BREVIS RADIX (BRX) activity in the developing protophloem. BRX is an auxin-responsive target of the prototypical ARF MONOPTEROS (MP), a key promoter of vascular development, and transiently enhances PIN3 expression to promote meristem growth in young roots. At later stages, cytokinin induction of SHY2 in the vascular transition zone restricts BRX expression to down-regulate PIN3 and thus limit meristem growth. Interestingly, proper SHY2 expression requires BRX, which could reflect feedback on the auxin responsiveness of SHY2 because BRX protein can directly interact with MP, likely acting as a cofactor. Thus, cross-regulatory antagonism between BRX and SHY2 could determine ARF activity in the protophloem. Our data suggest a model in which the regulatory interactions favor BRX expression in the early proximal meristem and SHY2 prevails because of supplementary cytokinin induction in the later distal meristem. The complex equilibrium of this regulatory module might represent a universal switch in the transition toward differentiation in various developmental contexts.

Transcriptional regulatory patterns of the myelin basic protein and malic enzyme genes by the thyroid hormone receptors alpha1 and beta1.

While there is evidence that the two ubiquitously expressed thyroid hormone (T3) receptors, TRalpha1 and TRbeta1, have distinct functional specificities, the mechanism by which they discriminate potential target genes remains largely unexplained. In this study, we demonstrate that the thyroid hormone response elements (TRE) from the malic enzyme and myelin basic protein genes (METRE and MBPTRE) respectively, are not functionally equivalent. The METRE, which is a direct repeat motif with a 4-base pair gap between the two half-site hexamers binds thyroid hormone receptor as a heterodimer with 9-cis-retinoic acid receptor (RXR) and mediates a high T3-dependent activation in response to TRalpha1 or TRbeta1 in NIH3T3 cells. In contrast, the MBPTRE, which consists of an inverted palindrome formed by two hexamers spaced by 6 base pairs, confers an efficient transactivation by TRbeta1 but a poor transactivation by TRalpha1. While both receptors form heterodimers with RXR on MBPTRE, the poor transactivation by TRalpha1 correlates also with its ability to bind efficiently as a monomer. This monomer, which is only observed with TRalpha1 bound to MBPTRE, interacts neither with N-CoR nor with SRC-1, explaining its functional inefficacy. However, in Xenopus oocytes, in which RXR proteins are not detectable, the transactivation mediated by TRalpha1 and TRbeta1 is equivalent and independent of a RXR supply, raising the question of the identity of the thyroid hormone receptor partner in these cells. Thus, in mammalian cells, the binding characteristics of TRalpha1 to MBPTRE (i.e. high monomer binding efficiency and low transactivation activity) might explain the particular pattern of T3 responsiveness of MBP gene expression during central nervous system development.

Control of gene expression in melanocytes and RPE by conserved distal regulatory elements

Résumé Durant le développement embryonnaire, les cellules pigmentaires des mammifères se développent à partir de deux origines différentes : les melanocytes se développent à partir de la crête neurale alors que les cellules de la rétine pigmentaire (RP) ont une origine neuronale. Un grand nombre de gènes sont impliqués dans la pigmentation dont les gènes de la famille tyrosinase à savoir Tyr, Tyrp1 et Dct. Certaines études ont suggéré que les gènes de la pigmentation sont régulés de manière différentielle dans les mélanocytes et dans la RP. Dans ce travail, les gènes de la famille tyrosinase ont été étudiés comme modèle de la régulation des gènes de la pigmentation par des éléments régulateurs agissant à distance. II a été montré que le promoteur du gène Tyrp1pouvait induire l'expression d'un transgène uniquement dans la RP alors que ce gène est aussi exprimé dans les mélanocytes comme le montre le phénotype des souris mutantes pour Tyrp1. Ce résultat suggère que les éléments régulateurs du promoteur sont suffisants pour l'expression dans la RP mais pas pour l'expression dans les mélanocytes. J'ai donc cherché à identifier la séquence qui régule l'expression dans les mélanocytes. Un chromosome artificiel bactérien (CAB) contenant le gène Tyrp1 s'est avéré suffisant pour induire l'expression dans les mélanocytes, comme démontré par la correction du phénotype mutant. La séquence de ce CAB contient plusieurs régions très conservées qui pourraient représenter de nouveaux éléments régulateurs. Par la suite, j'ai focalisé mon analyse sur une séquence située à -I5 kb qui s'est révélée être un amplificateur spécifique aux mélanocytes comme démontré par des expériences de cultures cellulaire et de transgenèse. De plus, une analyse poussée de cet élément a révélé que le facteur de transcription Sox 10 représentait un transactivateur de cet amplificateur. Comme pour Tyrp1, la régulation du gène tyrosinase est contrôlée par différents éléments régulateurs dans les mélanocytes et la RP. Il a été montré que le promoteur de tyrosinase n'était pas suffisant pour une forte expression dans les mélanocytes et la RP. De plus, l'analyse de la région située en amont a révélé la présence d'un amplificateur nécessaire à l'expression dans les mélanocytes à la position -15 kb. Cet amplificateur n'est toutefois pas actif dans la RP mais agit comme un répresseur dans ces cellules. Ces résultats indiquent que certains éléments nécessaires à l'expression dans les deux types de cellules pigmentaires sont absents de ces constructions. Comme pour Tyrp1, j'ai en premier lieu démontré qu'un CAB était capable de corriger le phénotype albinique, puis ai inséré un gène reporter (lacZ) dans le CAB par recombinaison homologue et ai finalement analysé l'expression du reporter en transgenèse. Ces souris ont montré une expression forte du lacZ dans les mélanocytes et la RP, ce qui indique que le CAB contient les séquences régulatrices nécessaires à l'expression correcte de tyrosinase. Afin de localiser plus précisément les éléments régulateurs, j'ai ensuite généré des délétions dans le CAB et analysé l'expression du lacZ en transgenèse. La comparaison de séquences génomiques provenant de différentes espèces a permis par la suite d'identifier des régions représentant de nouveaux éléments régulateurs potentiels. En utilisant cette approche, j'ai identifié une région qui se comporte comme un amplificateur dans la RP et qui est nécessaire à l'expression de tyrosinase dans ce tissu. De plus, j'ai identifié les facteurs de transcription Mitf et Sox10 comme transactivateurs de l'amplificateur spécifique aux mélanocytes situé à -15 kb. L'identification et la caractérisation des ces éléments régulateurs des gènes tyrosinase et Tyrp1confirme donc que la régulation différentielle des gènes dans les mélanocytes et la RP est liée à des éléments régulateurs séparés. Summary Pigment cells of mammals originate from two different lineages: melanocytes arise from the neural crest, whereas cells of the retinal pigment epithelium (RPE) originate from the optic cup of the developing forebrain. A large set of genes are involved in pigmentation, including the members of the tyrosinase gene family, namely tyrosinase, Tyrp1 and Dct. Previous studies have suggested that pigmentation genes are differentially regulated in melanocytes and RPE. In this work, the tyrosinase gene family was used as a model for studying the involvement of distal regulatory elements in pigment cell-specific gene expression. The promoter of the Tyrp1 gene has been shown to drive detectable transgene expression only to the RPE, even though the gene is also expressed in melanocytes as evident from Tyrp1-mutant mice. This indicates that the regulatory elements responsible for Tyrp1 gene expression in the RPE are not sufficient for expression in melanocytes. I thus searched for a putative melanocyte-specific regulatory sequence and demonstrate that a bacterial artificial chromosome (BAC) containing the Tyrp1 gene and surrounding sequences is able to target transgenic expression to melanocytes and to rescue the Tyrp1 b (brown) phenotype. This BAC contains several highly conserved non-coding sequences that might represent novel regulatory elements. I further focused on a sequence located at -15 kb which I identified as amelanocyte-specific enhancer as shown by cell culture and transgenic mice. In addition, further functional analysis identified the transcription factor Sox10 as being able to bind and transactivate this enhancer. As for Tyrp1, tyrosinase gene regulation is mediated by different cis-regulatory elements in melanocytes and RPE. It was shown that the tyrosinase promoter was not sufficient to confer strong and specific expression in melanocytes and RPE. Moreover, analysis of tyrosinase upstream sequence, revealed the presence of a specific enhancer at position -15 kb which was necessary to confer strong expression in melanocytes. This enhancer element however failed to act as an enhancer in the RPE, but rather repressed expression. This indicates that some regulatory elements required for tyrosinase expression in both RPE and melanocytes are still missing from these constructs. As for Tyrp1, I first demonstrated that a BAC containing the Tyr gene is able to rescue the Tyr c (albino) phenotype in mice, then I inserted a lacZ reporter gene in the BAC by homologous recombination, and finally analysed the pattern of lacZ expression in transgenic mice. These mice showed strong lacZ expression in both RPE and melanocytes, indicating that the BAC contains the regulatory sequences required for proper tyrosinase expression. In order to localize more precisely these regulatory elements, I have then generated several deletions in the BAC and analysed lacZ expression in transgenic mice. Multi-species comparative genomic analysis then allowed identifying conserved sequences that potentially represent novel regulatory elements. Using this experimental approach, I identified a region that behaves as a RPE-specific enhancer and that is required for tyrosinase expression in the retina] pigment epithelium. In addition, I identified the transcription factors Mitf and Sox l0 as being transactivators of the melanocyte-specific enhancer located at -l5 kb. The identification and characterization of these tyrosinase and Tyrp1 distal regulatory element supports the idea that separate regulatory sequences mediate differential gene expression in melanocytes and RPE.

Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.