Transcription factor binding site orientation and order are major drivers of gene regulatory activity
Open accessAbstract
The gene regulatory code and grammar remain largely unknown, precluding our ability to link phenotype to genotype in regulatory sequences. Here, using a massively parallel reporter assay (MPRA) of 209,440 sequences, we examine all possible pair and triplet combinations, permutations and orientations of eighteen liver-associated transcription factor binding sites (TFBS). We find that TFBS orientation and order have a major effect on gene regulatory activity. Corroborating these results with genomic analyses, we find clear human promoter TFBS orientation biases and similar TFBS orientation and order transcriptional effects in an MPRA that tested 164,307 liver candidate regulatory elements. Additionally, by adding TFBS orientation to a model that predicts expression from sequence we improve performance by 7.7%. Collectively, our results show that TFBS orientation and order have a significant effect on gene regulatory activity and need to be considered when analyzing the functional effect of variants on the activity of these sequences.
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Integrity signals
Registry
Cites
Statistics
Claims the paper makes
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Claims other papers make about this one
Each is checked against this paper’s own text: does it say what the citing paper says it does?
The bottom-up synthetic approach can provide data of sufficient statistical power to characterize combinatorial action of TFBSs, as demonstrated by recent studies that have profiled tens of thousands to millions of synthetic or even random DNA sequences in cell lines or yeast.
“and it can provide data of sufficient statistical power to characterize combinatorial action of TFBSs, as demonstrated by recent studies that have profiled tens of thousands to millions of synthetic or even random DNA sequences in cell lines or yeast.33,35,36”
From Design principles of cell-state-specific enhancers in hematopoiesis (2025) · cites this paper as [33]
Measuring enhancer activity of synthetic, systematically designed DNA sequences with TFBSs placed in an otherwise inert sequence context is an alternative way to study the cis-regulatory code.
“An alternative way to study the cis-regulatory code is to measure the enhancer activity of synthetic, systematically designed DNA sequences, where TF binding sites (TFBSs) are placed in an otherwise inert sequence context.27,30,31,32,33,34”
From Design principles of cell-state-specific enhancers in hematopoiesis (2025) · cites this paper as [33]
Self-check
Percentages, sums and reported differences re-computed from the paper’s own numbers.
The self-check needs the paper’s full text.
Cites
Systematic Localization of Common Disease-Associated Variation in Regulatory DNA
Maurano, M. T., Humbert, R., Rynes, E., Thurman, R. E., Haugen, E., Wang, H., Reynolds, A. P., Sandstrom, R., Qu, H., Brody, J., et al. (2012) · Science · doi:10.1126/science.1222794 · Cited by 3,716 · see citing works
Base-resolution models of transcription-factor binding reveal soft motif syntax
Avsec, Ž., Weilert, M., Shrikumar, A., Krueger, S., Alexandari, A., Dalal, K., Fropf, R., McAnany, C., Gagneur, J., Kundaje, A., et al. (2021) · Nature Genetics · doi:10.1038/s41588-021-00782-6 · Cited by 610 · cited 1× · see citing works
“In addition, in previous work, we identi fi ed widespread transcriptional strand asymmetries in the orientation of TFBSs across human transcribed regions 10 , a research area in which there have been a few studies to date 11 – 15 .”
Decoding enhancers using massively parallel reporter assays
Inoue, F., & Ahituv, N. (2015) · Genomics · doi:10.1016/j.ygeno.2015.06.005 · Cited by 272 · see citing works
Self-citationA systematic evaluation of the design and context dependencies of massively parallel reporter assays
Klein, J. C., Agarwal, V., Inoue, F., Keith, A., Martin, B., Kircher, M., Ahituv, N., & Shendure, J. (2020) · Nature Methods · doi:10.1038/s41592-020-0965-y · Cited by 226 · see citing works
Grossman, S. R., Zhang, X., Wang, L., Engreitz, J., Melnikov, A., Rogov, P., Tewhey, R., Isakova, A., Deplancke, B., Bernstein, B. E., et al. (2017) · Proceedings of the National Academy of Sciences · doi:10.1073/pnas.1621150114 · Cited by 165 · see citing works
Systematic interrogation of human promoters
Weingarten-Gabbay, S., Nir, R., Lubliner, S., Sharon, E., Kalma, Y., Weinberger, A., & Segal, E. (2019) · Genome Research · doi:10.1101/gr.236075.118 · Cited by 110 · cited 2× · see citing works
“Previous work has examined the effect of genomic element orientation on expression in both single locus experimental studies 7,32 and via high-throughput assays 9,13,18,33 .”
Flagged first, then by citations
Research works printed but not linked: 39
- [1]Chatterjee S Ahituv N Gene regulatory elements, major drivers of human disease Annu. Rev. Genomics Hum. Genet. 2017 18 45 63 10.1146/annurev-genom-091416-035537 28399667 · DOI
- [2]Manolio TA Finding the missing heritability of complex diseases Nature 2009 461 747 753 10.1038/nature08494 19812666 PMC2831613 · DOI · Paper page
- [4]Shendure J Akey JM The origins, determinants, and consequences of human mutations Science 2015 349 1478 1483 10.1126/science.aaa9119 26404824 · DOI
- [5]Matharu N Ahituv N Modulating gene regulation to treat genetic disorders Nat. Rev. Drug Discov. 2020 19 757 775 10.1038/s41573-020-0083-7 33020616 PMC12922764 · DOI
- [6]Spitz F Furlong EEM Transcription factors: from enhancer binding to developmental control Nat. Rev. Genet. 2012 13 613 626 10.1038/nrg3207 22868264 · DOI · Paper page
- [7]Panne D The enhanceosome Curr. Opin. Struct. Biol. 2008 18 236 242 10.1016/j.sbi.2007.12.002 18206362 · DOI
- [8]Kulkarni MM Arnosti DN Information display by transcriptional enhancers Development 2003 130 6569 6575 10.1242/dev.00890 14660545 · DOI · Paper page
- [9]Jolma A DNA-dependent formation of transcription factor pairs alters their binding specificity Nature 2015 527 384 388 10.1038/nature15518 26550823 · DOI
- [10]Georgakopoulos-Soares I Asymmetron: a toolkit for the identification of strand asymmetry patterns in biological sequences Nucleic Acids Res. 2021 49 e4 10.1093/nar/gkaa1052 33211865 PMC7797064 · DOI
- [11]Zhang L Guarente L The yeast activator HAP1–a GAL4 family member–binds DNA in a directly repeated orientation Genes Dev. 1994 8 2110 2119 10.1101/gad.8.17.2110 7958882 · DOI
- [12]King DA Zhang L Guarente L Marmorstein R Structure of a HAP1-DNA complex reveals dramatically asymmetric DNA binding by a homodimeric protein Nat. Struct. Biol. 1999 6 64 71 10.1038/4940 9886294 · DOI
- [14]Tippens ND Transcription imparts architecture, function and logic to enhancer units Nat. Genet. 2020 52 1067 1075 10.1038/s41588-020-0686-2 32958950 PMC7541647 · DOI
- [18]Roberts BS Genome-wide strand asymmetry in massively parallel reporter activity favors genic strands Genome Res. 2021 31 866 876 10.1101/gr.270751.120 33879525 PMC8092006 · DOI
- [19]Smith RP Massively parallel decoding of mammalian regulatory sequences supports a flexible organizational model Nat. Genet. 2013 45 1021 1028 10.1038/ng.2713 23892608 PMC3775494 · DOI · Paper page
- [20]Castro-Mondragon JA JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles Nucleic Acids Res. 2022 50 D165 D173 10.1093/nar/gkab1113 34850907 PMC8728201 · DOI · Paper page
- [21]Kulakovskiy IV HOCOMOCO: towards a complete collection of transcription factor binding models for human and mouse via large-scale ChIP-Seq analysis Nucleic Acids Res. 2018 46 D252 D259 10.1093/nar/gkx1106 29140464 PMC5753240 · DOI
- [22]Inoue F A systematic comparison reveals substantial differences in chromosomal versus episomal encoding of enhancer activity Genome Res. 2017 27 38 52 10.1101/gr.212092.116 27831498 PMC5204343 · DOI · Paper page
- [23]Gordon MG lentiMPRA and MPRAflow for high-throughput functional characterization of gene regulatory elements Nat. Protoc. 2020 15 2387 2412 10.1038/s41596-020-0333-5 32641802 PMC7550205 · DOI
- [24]Chong JA REST: A mammalian silencer protein that restricts sodium channel gene expression to neurons Cell 1995 80 949 957 10.1016/0092-8674(95)90298-8 7697725 · DOI
- [25]Puig RR Boddie P Khan A Castro-Mondragon JA Mathelier A UniBind: maps of high-confidence direct TF-DNA interactions across nine species BMC Genomics 2021 22 482 10.1186/s12864-021-07760-6 34174819 PMC8236138 · DOI
- [26]ENCODE Project Consortium Expanded encyclopaedias of DNA elements in the human and mouse genomes Nature 2020 583 699 710 10.1038/s41586-020-2493-4 32728249 PMC7410828 · DOI · Paper page
- [27]ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature489, 57–74 (2012). PMC3439153 22955616
- [28]Agarwal, V. et al. Massively parallel characterization of transcriptional regulatory elements in three diverse human cell types. bioRxiv 2023.03.05.531189. 10.1101/2023.03.05.531189. · DOI
- [29]Agarwal V Shendure J Predicting mRNA abundance directly from genomic sequence using deep convolutional neural networks Cell Rep. 2020 31 107663 10.1016/j.celrep.2020.107663 32433972 · DOI · Paper page
- [30]Zhou J Deep learning sequence-based ab initio prediction of variant effects on expression and disease risk Nat. Genet. 2018 50 1171 1179 10.1038/s41588-018-0160-6 30013180 PMC6094955 · DOI · Paper page
- [31]Avsec Ž Effective gene expression prediction from sequence by integrating long-range interactions Nat. Methods 2021 18 1196 1203 10.1038/s41592-021-01252-x 34608324 PMC8490152 · DOI · Paper page
- [32]Natesan S Gilman MZ DNA bending and orientation-dependent function of YY1 in the c-fos promoter Genes Dev. 1993 7 2497 2509 10.1101/gad.7.12b.2497 8276234 · DOI
- [34]Chytil, M., Peterson, B. R., Erlanson, D. A. & Verdine, G. L. The orientation of the AP-1 heterodimer on DNA strongly affects transcriptional potency. Proc. Natl Acad. Sci.95 14076–14081 (1998). PMC24329 9826656
- [35]Falvo JV Parekh BS Lin CH Fraenkel E Maniatis T Assembly of a functional beta interferon enhanceosome is dependent on ATF-2–c-jun heterodimer orientation Mol. Cell. Biol. 2000 20 4814 4825 10.1128/MCB.20.13.4814-4825.2000 10848607 PMC85927 · DOI
- [36]Krivan W Wasserman WW A predictive model for regulatory sequences directing liver-specific transcription Genome Res. 2001 11 1559 1566 10.1101/gr.180601 11544200 PMC311083 · DOI
- [37]Inoue F Kreimer A Ashuach T Ahituv N Yosef N Identification and massively parallel characterization of regulatory elements driving neural induction Cell Stem Cell 2019 25 713 727.e10 10.1016/j.stem.2019.09.010 31631012 PMC6850896 · DOI
- [38]Li H Durbin R Fast and accurate short read alignment with Burrows-Wheeler transform Bioinformatics 2009 25 1754 1760 10.1093/bioinformatics/btp324 19451168 PMC2705234 · DOI · Paper page
- [39]Frankish A GENCODE 2021 Nucleic Acids Res. 2021 49 D916 D923 10.1093/nar/gkaa1087 33270111 PMC7778937 · DOI · Paper page
- [40]Roadmap Epigenomics Consortium Integrative analysis of 111 reference human epigenomes Nature 2015 518 317 330 10.1038/nature14248 25693563 PMC4530010 · DOI · Paper page
- [41]Grant, C. E., Bailey, T. L. & Noble, W. S. FIMO: scanning for occurrences of a given motif. Bioinformatics27, 1017–1018 (2011). PMC3065696 21330290
- [42]Jiang M Anderson J Gillespie J Mayne M uShuffle: a useful tool for shuffling biological sequences while preserving the k-let counts BMC Bioinf. 2008 9 192 10.1186/1471-2105-9-192 PMC2375906 18405375 · DOI
- [43]Glover JN Harrison SC Crystal structure of the heterodimeric bZIP transcription factor c-Fos-c-Jun bound to DNA Nature 1995 373 257 261 10.1038/373257a0 7816143 · DOI
- [44]Carlezon WA Jr Duman RS Nestler EJ The many faces of CREB Trends Neurosci. 2005 28 436 445 10.1016/j.tins.2005.06.005 15982754 · DOI
- [45]Fraser JD Martinez V Straney R Briggs MR DNA binding and transcription activation specificity of hepatocyte nuclear factor 4 Nucleic Acids Res. 1998 26 2702 2707 10.1093/nar/26.11.2702 9592157 PMC147588 · DOI
Cited by
Sorted by citations
Design principles of cell-state-specific enhancers in hematopoiesis
Frömel, R., Rühle, J., Bernal Martinez, A., Szu-Tu, C., Pacheco Pastor, F., Martinez-Corral, R., & Velten, L. (2025) · Cell · doi:10.1016/j.cell.2025.04.017 · Cited by 29 · see citing works
Orthrus: toward evolutionary and functional RNA foundation models
Fradkin, P., Shi, R. “., Dalal, T., Isaev, K., Frey, B. J., Lee, L. J., Morris, Q., & Wang, B. (2026) · Nature Methods · doi:10.1038/s41592-026-03064-3 · Cited by 22 · see citing works
DNA-BFN: Bayesian Flow Networks for Generative Modeling of Regulatory DNA Sequences
Yu, H., Du, Q., Wei, L., & Wang, X. (2026) · Bioinformatics · doi:10.1093/bioinformatics/btag785
- Authors
- Georgakopoulos-Soares, I., Deng, C., Agarwal, V., Chan, C. S. Y., Zhao, J., Inoue, F., & Ahituv, N.
- Affiliations
- University of California, San Francisco (US)
- Kyoto University (JP)
- Corresponding author
- Ahituv, Nadav
- ORCID
- Venue
- Nature Communications
- Published
- 2023-04-22
- Volume
- 14
- Issue
- 1
- Pages
- 2333-2333
- ISSN
- 2041-1723 (ISSN-L)
- PMC
- PMC10122648
- OpenAlex
- W4366769693
- Cited by
- 82
- References
- 6
- Citations by year
- 2026: 31 · 2025: 25 · 2024: 24 · 2023: 2
- FWCI
- 7.09
- Type
- journal article
- Language
- en
- Open access
- gold
- doi.org
- Locations
- Nature Communications (publishedVersion, cc-by)
- PubMed Central (publishedVersion, cc-by)
- PubMed (publishedVersion)
- RePEc: Research Papers in Economics (submittedVersion)
- DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
- eScholarship (California Digital Library) (submittedVersion, cc-by)
- Funders
- U.S. Department of Health & Human Services | NIH | National Human Genome Research Institute (1UM1HG011966)
- Topics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Chromatin Dynamics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene Regulatory Network Analysis
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic Associations and Epidemiology
- Keywords
- transcription factor binding sites, gene expression regulation, promoter, massively parallel reporter assay, regulatory variants, transcriptional regulation
- MeSH terms
- Binding Sites, Gene Expression Regulation, Humans, Promoter Regions, Genetic, Protein Binding, Transcription Factors
About these checks
Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.
Statistics and arithmetic are recomputed, not judged.
Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.
What counts as a claim ›How we measure confidence ›How we check this ›
Are you an author, or do you see a problem with this page? Request a review