Mozziconacci, Merle and Lesne proposed that the developmental gene regulatory code — the mapping from combinations of transcription-factor (TF) concentrations to gene expression levels — is mediated by an adaptor: the gene domain, the genomic segment delimited by insulator/boundary elements (DNA elements bound by proteins such as CTCF that limit chromatin contacts across them) that contains a developmental gene and its enhancers and manifests as a contact domain in 3D genome maps. Whether this code can be explicitly cracked remains open. The unresolved core is quantitative and mechanistic: (i) how the composition and arrangement of TF binding sites at a single enhancer implement a quantitative reading of TF concentrations; (ii) the integration logic by which multiple enhancers within one domain jointly regulate the promoter — additive, competitive, exclusive-OR-like, or other nonlinear rules — and how this logic depends on enhancer order, genomic distance, and domain topology; and (iii) how the 3D folding dynamics and boundaries of the domain modulate enhancer–promoter contact probabilities in single cells while tissue-level expression stays robust. Substantial progress has been made on these threads — live imaging of the Drosophila even-skipped locus, shadow-enhancer analyses (shadow enhancers being seemingly redundant enhancers driving the same expression pattern), CRISPR-perturbation enhancer–promoter maps, synthetic enhancer platforms, and single-cell joint chromatin-structure-plus-transcription assays — but no general decoding rules or settled integration logic have been established, and expression phase diagrams remain characterized for almost no genes beyond eve.

Background

In embryonic development, cells in the same tissue express the same genes differently because each gene responds to the local concentrations of transcription factors (TFs) — proteins that bind short DNA sequence motifs and activate or repress transcription. Lewis Wolpert's positional-information concept (1969) holds that such local TF concentrations act as a 'positional value' that cells decode; the question of how this decoding works quantitatively is a central problem of developmental biology.

Enhancers are short regulatory DNA regions (typically a few hundred base pairs) bearing clusters of TF binding sites, which activate transcription of a target gene's promoter (the gene's transcription start control region) over large genomic distances in an orientation- and distance-tolerant way. Most developmental genes in animals are regulated not by one enhancer but by several, embedded with the gene in a genomic interval demarcated by insulators (also called boundary elements) — DNA elements bound by proteins such as CTCF that limit contacts across them. Genome-wide chromosome-conformation-capture assays (Hi-C, which cross-links, ligates, and sequences physically close DNA fragments to estimate pairwise contact frequency across a cell population) reveal topologically associating domains (TADs): genomic intervals whose internal contact frequency far exceeds contacts across their borders. Mozziconacci, Merle and Lesne (2020) call the insulator-delimited segment containing a developmental gene and its enhancers the gene domain, and propose it as the adaptor of the gene regulatory code: by analogy with the genetic code and its tRNA adaptors, the code maps codewords (local concentrations of specific TFs) onto encoded objects (gene expression levels), and the gene domain is the intermediary physical entity that makes this mapping gratuitous (not fixed by physicochemical necessity) and evolvable. They visualize the code as an expression phase diagram: a map of gene expression level over the space of TF concentration combinations, with 'active phases' (well-expressed regions) and 'inactive phases'.

The paradigmatic example is the Drosophila even-skipped (eve) gene. In the early embryo, maternal gradients establish nonuniform concentrations of gap-gene TFs (the gap genes Krüppel, Hunchback, Giant and Knirps, expressed in broad overlapping bands along the embryo axis), and eve responds with seven expression stripes along the embryo axis; five enhancers in the eve gene domain are each predominantly responsible for one or two stripes. Yet enhancers are not independent modules: Bothma et al. (2015) showed in living embryos that whether two enhancers act additively or non-additively depends on their strengths; Kim et al. (2013) found that rearranging eve enhancers changes expression output, so enhancer order matters; and Lim et al. (2018) showed that deleting one stripe's enhancer also shifts a neighboring stripe. The gene domain is more than the sum of its enhancers. Transcription at these loci occurs in discontinuous bursts (episodes) of mRNA production, measurable in living embryos by MS2 tagging — a technique that inserts MS2 RNA stem loops into a gene so nascent transcripts fluoresce upon binding of a fluorescent coat protein, giving single-cell, minute-resolution transcription histories. Live imaging (Lammers et al. 2020; Berrocal et al. 2020, 2024) shows eve stripes are shaped through modulation of burst frequency and amplitude, but the molecular mechanism implementing this, and how several simultaneously active enhancers combine their effects, remain unknown.

In parallel, genome-scale perturbation work connected 3D folding to output: Fulco et al. (2019) derived the activity-by-contact (ABC) rule — an enhancer's effect on a gene scales with the product of its regulatory activity and its Hi-C contact frequency with the promoter — an accurate but phenomenological, population-averaged predictor, not a mechanistic decoding rule. Whether enhancer–promoter contact even is the causal event, and how 3D folding mechanisms mediate enhancer selectivity for promoters, is still debated (Yang and Hansen 2024); for the eve domain boundary itself, Bing et al. (2024) concluded that domain borders in Drosophila are set by physical pairing between boundary elements rather than by the loop-extrusion mechanism (a model in which the cohesin protein complex progressively enlarges chromatin loops until halted at a boundary) prevailing in vertebrates. New single-cell methods that measure chromatin architecture and transcription in the same cell (imaging: Mateo et al. 2019; sequencing: Qu et al. 2023) now make the joint folding-plus-expression measurements foreseeable that the source paper called for. Against this background, the source paper's open question stands: can the adaptor-mediated code be made explicit — quantitatively and mechanistically — for developmental gene domains in general, beyond the exemplified eve locus?

Problem Statement

For developmental gene domains — insulator-delimited genomic segments containing a developmental gene and its enhancers, exemplified by the Drosophila even-skipped (eve) locus — determine the quantitative decoding rules of the proposed gene-domain adaptor of the regulatory code. A complete answer must specify: (1) how the composition and arrangement of transcription-factor (TF) binding sites at a single enhancer implement a quantitative reading of TF concentrations, i.e., the local input function from the TF concentration codeword to enhancer activity; (2) the integration logic by which the activities of multiple enhancers within the same domain are combined at the shared promoter — whether additive, competitive, exclusive-OR-like, or another nonlinear rule — including how this logic depends on enhancer order, genomic distance, and domain topology, and what the transcriptional output is when two or more enhancers are jointly active; and (3) how the 3D folding of the domain and its boundary elements modulate enhancer-promoter contact probabilities and transcriptional output at the single-cell level, and how cell-to-cell variability in domain conformation is compatible with robust average expression at tissue level. The common objective is an explicit, predictive input-output code mapping TF-concentration combinations to gene expression (an experimentally determined expression phase diagram, or a quantitative model that predicts expression levels over the relevant TF-concentration combinations) for developmental genes beyond the exemplified eve locus: general decoding rules instantiated per gene domain, in the sense of the source's program of characterizing expression phase diagrams for the many genes involved in development, rather than a single universal mapping fixed at once for all genes; a complete answer requires these rules to be instantiated for at least one developmental gene domain beyond eve, while extending them across the full set of developmental genes and to other developmental contexts and species remains a longer-term goal.

References

  1. Mozziconacci J, Merle M, Lesne A. The 3D genome shapes the regulatory code of developmental genes. Journal of Molecular Biology 432(3):712-723 (2020). doi:10.1016/j.jmb.2019.10.017, https://doi.org/10.1016/j.jmb.2019.10.017 (accepted manuscript: https://hal.science/hal-02357466)
  2. Bothma JP, Garcia HG, Ng S, Perry MW, Gregor T, Levine M. Enhancer additivity and non-additivity are determined by enhancer strength in the Drosophila embryo. eLife 4:e07956 (2015). doi:10.7554/eLife.07956, https://doi.org/10.7554/eLife.07956
  3. Kim AR, Martinez C, Ionides J, Ramos AF, Ludwig MZ, Ogawa N, Sharp DH, Reinitz J. Rearrangements of 2.5 kilobases of noncoding DNA from the Drosophila even-skipped locus define predictive rules of genomic cis-regulatory logic. PLoS Genetics 9(2):e1003243 (2013). doi:10.1371/journal.pgen.1003243, https://doi.org/10.1371/journal.pgen.1003243
  4. Lim B, Fukaya T, Heist T, Levine M. Temporal dynamics of pair-rule stripes in living Drosophila embryos. Proceedings of the National Academy of Sciences USA 115(33):8376-8381 (2018). doi:10.1073/pnas.1810430115, https://doi.org/10.1073/pnas.1810430115
  5. Mateo LJ, Murphy SE, Hafner A, Cinquini IS, Walker CA, Boettiger AN. Visualizing DNA folding and RNA in embryos at single-cell resolution. Nature 568:49-54 (2019). doi:10.1038/s41586-019-1035-4, https://doi.org/10.1038/s41586-019-1035-4
  6. Fulco CP, Nasser J, Jones TR, Munson G, Bergman DT, Subramanian V, Grossman SR, Anyoha R, Doughty BR, Patwardhan TA, Nguyen TH, Kane M, Perez EM, Durand NC, Lareau CA, Stamenova EK, Aiden EL, Lander ES, Engreitz JM. Activity-by-contact model of enhancer-promoter regulation from thousands of CRISPR perturbations. Nature Genetics 51(12):1664-1669 (2019). doi:10.1038/s41588-019-0538-0, https://doi.org/10.1038/s41588-019-0538-0
  7. Lammers NC, Galstyan V, Reimer A, Medin SA, Wiggins CH, Garcia HG. Multimodal transcriptional control of pattern formation in embryonic development. Proceedings of the National Academy of Sciences USA 117(2):836-847 (2020). doi:10.1073/pnas.1912500117, https://doi.org/10.1073/pnas.1912500117
  8. Berrocal A, Lammers NC, Garcia HG, Eisen MB. Kinetic sculpting of the seven stripes of the Drosophila even-skipped gene. eLife 9:e61635 (2020). doi:10.7554/eLife.61635, https://doi.org/10.7554/eLife.61635
  9. Waymack R, Fletcher A, Enciso G, Wunderlich Z. Shadow enhancers can suppress input transcription factor noise through distinct regulatory logic. eLife 9:e59351 (2020). doi:10.7554/eLife.59351, https://doi.org/10.7554/eLife.59351
  10. Kvon EZ, Waymack R, Gad M, Wunderlich Z. Enhancer redundancy in development and disease. Nature Reviews Genetics 22(5):324-336 (2021). doi:10.1038/s41576-020-00311-x, https://doi.org/10.1038/s41576-020-00311-x
  11. Whitney PH, Shrestha B, Xiong J, Zhang T, Rushlow CA. Shadow enhancers modulate distinct transcriptional parameters that differentially effect downstream patterning events. Development 149(21):dev200940 (2022). doi:10.1242/dev.200940, https://doi.org/10.1242/dev.200940
  12. Qu J, Sun J, Zhao C, Liu X, Zhang X, Jiang S, Wei C, Yu H, Zeng X, Fan L, Ding J. Simultaneous profiling of chromatin architecture and transcription in single cells. Nature Structural & Molecular Biology 30(9):1393-1402 (2023). doi:10.1038/s41594-023-01066-9, https://doi.org/10.1038/s41594-023-01066-9
  13. Berrocal A, Lammers NC, Garcia HG, Eisen MB. Unified bursting strategies in ectopic and endogenous even-skipped expression patterns. eLife 12:RP88671 (2024). doi:10.7554/eLife.88671, https://doi.org/10.7554/eLife.88671
  14. Yang JH, Hansen AS. Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation. Nature Reviews Molecular Cell Biology 25(7):574-591 (2024). doi:10.1038/s41580-024-00710-6, https://doi.org/10.1038/s41580-024-00710-6
  15. Bing X, Ke W, Fujioka M, Kurbidaeva A, Levitt S, Levine M, Schedl P, Jaynes JB. Chromosome structure in Drosophila is determined by boundary pairing not loop extrusion. eLife 13:RP94070 (2024). doi:10.7554/eLife.94070, https://doi.org/10.7554/eLife.94070
  16. Lin S, Lim B. Multifaceted effects on even-skipped transcriptional dynamics upon Krüppel dosage changes. Development 151(5):dev202132 (2024). doi:10.1242/dev.202132, https://doi.org/10.1242/dev.202132
  17. Thomas HF, Feng S, Haslhofer F, Huber M, García Gallardo M, Loubiere V, Vanina D, Pitasi M, Stark A, Buecker C. Enhancer cooperativity can compensate for loss of activity over large genomic distances. Molecular Cell 85(2):362-375 (2025). doi:10.1016/j.molcel.2024.11.008, https://doi.org/10.1016/j.molcel.2024.11.008
  18. Friman ET, Bickmore WA. Enhancer cooperativity in the folded genome. Current Opinion in Genetics & Development 96:102416 (2026). doi:10.1016/j.gde.2025.102416, https://doi.org/10.1016/j.gde.2025.102416