"Stripe" transcription factors provide accessibility to co-binding partners in mammalian genomes

Iterative deep learning design of human enhancers exploits condensed sequence grammar to achieve cell-type specificity

An important and largely unsolved problem in synthetic biology is how to target gene expression to specific cell types. Here, we apply iterative deep learning to design synthetic enhancers with strong differential activity between two human cell lines. We initially train models on published datasets of enhancer activity and chromatin accessibility and use them to guide the design of synthetic enhancers that maximize predicted specificity. We experimentally validate these sequences, use the measurements to re-optimize the model, and design a second generation of enhancers with improved specificity. Our design methods embed relevant transcription factor binding site (TFBS) motifs with higher frequency than comparable endogenous enhancers while using a more selective motif vocabulary, and we show that enhancer activity is correlated with transcription factor expression at the single-cell level. Finally, we characterize causal features of top enhancers via perturbation experiments and show that enhancers as short as 50 bp can maintain specificity. A record of this paper's transparent peer review process is included in the supplemental information.

Krüppel-like factor 5 remodels lipid metabolism in exercised skeletal muscle

Regular physical activity induces a variety of health benefits, preventing and counteracting diseases caused by a sedentary lifestyle. However, the molecular underpinnings of skeletal muscle plasticity in exercise remain poorly understood. We identified a role of the Krüppel-Like Factor 5 (Klf5) in this process, in particular in the regulation of lipid homeostasis. Surprisingly, gain- and loss-of-function studies in muscle in vivo revealed seemingly opposite functions of Klf5 in the response to an acute exercise bout and chronic training, modulating lipid oxidation and synthesis, respectively. Thus, even though only transiently induced, the function of Klf5 is complex and fundamental for a normal long-term training response. These findings highlight the importance of this mediator of external stress response to adaptive remodeling of skeletal muscle tissue.

CellWalker2: Multi-omic discovery using hierarchical cell type relationships

Tissues are composed of cells with a wide range of similarities to each other, yet existing methods for single-cell genomics treat cell types as discrete labels. To address this gap, we developed CellWalker2, a graph diffusion-based model for the annotation and mapping of multi-modal data. With our open-source software package, hierarchically related cell types can be probabilistically matched across contexts and used to annotate cells, genomic regions, or gene sets. Additional features include estimating statistical significance and enabling gene expression and chromatin accessibility to be jointly modeled. Through simulation studies, we show that CellWalker2 performs better than existing methods in cell-type annotation and mapping. We then use multi-omics data from the brain and immune system to demonstrate CellWalker2's ability to assign high-resolution cell-type labels to regulatory elements and TFs and to quantify both conserved and divergent cell-type relationships between species.

Binding of NF-Y to transposable elements in mouse and human cells

BACKGROUND

Transposable Elements (TEs) represent a sizeable amount of mammalian genomes, providing regulatory sequences involved in shaping gene expression patterns. NF-Y is a Transcription factor -TF- trimer that binds to the CCAAT box, belonging to a selected group implicated in determining initiation of coding and noncoding RNAs.

RESULTS

We focus on NF-Y TE locations in 8 human and 8 mouse cells. Binding is exclusive for retroviral LTR12, MLT1 and MER in human and RLTR10 and IAPLTR in mouse cells. Cobinding and analysis of the DNA matrices signal enrichment of distinct TFs neighboring CCAAT in the three TE classes: MAFK/F/G in LTR12 and USF1/2 in MLT1 with precise alignment of sites, PKNOX1, MEIS2, PBX2/3 TALE TFs in MER57. The presence of "epigenetic" marks in human cells indicate prevalent co-association with open chromatin in MER, closed in LTR12 and mixed in MLT1. Based on chromatin features, these locations are mostly marked as enhancers, as confirmed by analysis of loci predicted to generate eRNAs.

CONCLUSIONS

These results are discussed in the context of functional data, suggesting a complex -positive and potentially-negative role of NF-Y on distinct classes of repetitive sequences.

TF Profiler: a transcription factor inference method that broadly measures transcription factor activity and identifies mechanistically distinct networks

TF Profiler is a method of inferring transcription factor (TF) regulatory activity, i.e., when a TF is present and actively participating in the regulation of transcription, directly from nascent sequencing assays such as PRO-seq and GRO-seq. While ChIP assays have measured DNA localization, they fall short of identifying when and where the effector domain of a transcription factor is active. Our method uses RNA polymerase activity to infer TF effector domain activity across hundreds of data sets and transcription factors. TF Profiler is broadly applicable, providing regulatory insights on any PRO-seq sample for any transcription factor with a known binding motif.

Stimulus-Dependent Expression of Bdnf Is Mediated by ATF2, MYT1L, and EGR1 Transcription Factors

Neurotrophins like BDNF have a key role in the proper functioning of the central nervous system, influencing numerous processes like memory formation and behavior. An imbalance in BDNF levels can lead to a wide range of diseases, including depression and neurodevelopmental disorders. While the potential therapeutic effects of BDNF are well-recognized, there is a knowledge gap in understanding the mechanisms governing BDNF expression levels. Here, we focused on the regulation of Bdnf gene expression in response to different stimuli, specifically studying the effects of neuronal activity and BDNF-TrkB signaling on Bdnf transcription in cultured neurons from rats of either sex. We used in vitro DNA pulldown combined with mass spectrometry to determine transcription factors that interact with the Bdnf promoters upon different stimuli and validated numerous known regulators, such as USF and AP1 family, and novel candidate regulators using reporter assays. We show that the USF family of transcription factors is specifically recruited after membrane depolarization, whereas the AP1 family participates in Bdnf regulation only after BDNF-TrkB signaling. We further describe ATF2, MYT1L, and EGR family as novel regulators of Bdnf expression by demonstrating their direct binding to Bdnf promoters using chromatin immunoprecipitation assays both in vitro and in vivo, showing their functional role in Bdnf gene expression and ultimately identifying their regulatory cis-elements in Bdnf promoters. Furthermore, our results show competition between ATF2, CREB, and AP1 family in regulating Bdnf levels. Collectively, our results provide insight into the regulation of Bdnf expression upon different stimuli.

Molecular signature of primate astrocytes reveals pathways and regulatory changes contributing to human brain evolution

Astrocytes contribute to the development and regulation of the higher-level functions of the brain, the critical targets of evolution. However, how astrocytes evolve in primates is unsettled. Here, we obtain human, chimpanzee, and macaque induced pluripotent stem-cell-derived astrocytes (iAstrocytes). Human iAstrocytes are bigger and more complex than the non-human primate iAstrocytes. We identify new loci contributing to the increased human astrocyte. We show that genes and pathways implicated in long-range intercellular signaling are activated in the human iAstrocytes and partake in controlling iAstrocyte complexity. Genes downregulated in human iAstrocytes frequently relate to neurological disorders and were decreased in adult brain samples. Through regulome analysis and machine learning, we uncover that functional activation of enhancers coincides with a previously unappreciated, pervasive gain of "stripe" transcription factor binding sites. Altogether, we reveal the transcriptomic signature of primate astrocyte evolution and a mechanism driving the acquisition of the regulatory potential of enhancers.

Stripe out for brain evolution: Astrocytes go human

Evolution of the human brain involved neuron and glial cell changes. In this issue of Cell Stem Cell, Ciuba et al.1 identify an epigenetic mechanism central in the evolution of human astrocytes, where small enhancer modifications increased transcription factor binding, orchestrating transcriptomic changes that influence astrocyte structure and function.

Massively parallel characterization of transcriptional regulatory elements

The human genome contains millions of candidate cis-regulatory elements (cCREs) with cell-type-specific activities that shape both health and many disease states1. However, we lack a functional understanding of the sequence features that control the activity and cell-type-specific features of these cCREs. Here we used lentivirus-based massively parallel reporter assays (lentiMPRAs) to test the regulatory activity of more than 680,000 sequences, representing an extensive set of annotated cCREs among three cell types (HepG2, K562 and WTC11), and found that 41.7% of these sequences were active. By testing sequences in both orientations, we find promoters to have strand-orientation biases and their 200-nucleotide cores to function as non-cell-type-specific 'on switches' that provide similar expression levels to their associated gene. By contrast, enhancers have weaker orientation biases, but increased tissue-specific characteristics. Utilizing our lentiMPRA data, we develop sequence-based models to predict cCRE function and variant effects with high accuracy, delineate regulatory motifs and model their combinatorial effects. Testing a lentiMPRA library encompassing 60,000 cCREs in all three cell types further identified factors that determine cell-type specificity. Collectively, our work provides an extensive catalogue of functional CREs in three widely used cell lines and showcases how large-scale functional measurements can be used to dissect regulatory grammar.

TFIIH kinase CDK7 drives cell proliferation through a common core transcription factor network

How cyclin-dependent kinase 7 (CDK7) coordinately regulates the cell cycle and RNA polymerase II transcription remains unclear. Here, high-resolution cryo-electron microscopy revealed how two clinically relevant inhibitors block CDK7 function. In cells, CDK7 inhibition rapidly suppressed transcription, but constitutively active genes were disproportionately affected versus stimulus-responsive. Distinct transcription factors (TFs) regulate constitutive versus stimulus-responsive genes. Accordingly, stimulus-responsive TFs were refractory to CDK7 inhibition whereas constitutively active "core" TFs were repressed. Core TFs (n = 78) are predominantly promoter associated and control cell cycle and proliferative gene expression programs across cell types. Mechanistically, rapid suppression of core TF function can occur through CDK7-dependent phosphorylation changes in core TFs and RB1. Moreover, CDK7 inhibition depleted core TF protein levels within hours, consistent with durable target gene suppression. Thus, a major but unappreciated biological function for CDK7 is regulation of a TF cohort that drives proliferation, revealing an apparent universal mechanism by which CDK7 coordinates RNAPII transcription with cell cycle CDK regulation.

Histone Deacetylase Inhibitor Largazole Deactivates A Subset of Superenchancers and Causes Mitotic Chromosome Mis-alignment by Suppressing SP1 and BRD4

Histone deacetylase inhibitors have been investigated as potential therapeutic agents for cancer and other diseases. HDIs are known to promote histone acetylation, resulting in an open chromatin conformation and generally increased gene expression. In previous work, we reported that a subset of genes, particularly those regulated by superenhancers, can be suppressed by the HDAC inhibitor largazole. To elucidate the molecular mechanisms underlying gene repression by largazole, we conducted transposase-accessible chromatin sequencing, ChIP-seq, and RNA-seq studies. Our findings revealed that while largazole treatment generally enhances chromatin accessibility, it selectively decreases the accessibility of a subset of superenhancer regions. These genomic regions, showing the most significant changes in the presence of largazole, were enriched with transcription factor binding motifs for SP1, BRD4, CTCF, and YY1. ChIP-seq analysis confirmed reduced binding of BRD4 and SP1 at their respective sites on chromatin, particularly at superenhancers regulating genes such as ID1, c-Myc and MCMs. Largazole exerts its effects by inhibiting DNA replication, RNA processing, and cell cycle progression, partially through the suppression of SP1 expression. Depletion of SP1 by shRNA mimics several key biological effects of largazole and increases cellular sensitivity to the drug. Specific to cell cycle regulation, we demonstrated that largazole disrupts G/M transition by interfering with chromosome alignment during metaphase, a phenotype also observed with SP1 depletion. Our results suggest that largazole exerts its growth-inhibitory effect by suppressing BRD4 and SP1 at super-enhancers, leading to cytostatic responses and mitotic dysfunction.

Transcription factor specificity protein (SP) family in renal physiology and diseases

Dysregulated specificity proteins (SPs), members of the C2H2 zinc-finger family, are crucial transcription factors (TFs) with implications for renal physiology and diseases. This comprehensive review focuses on the role of SP family members, particularly SP1 and SP3, in renal physiology and pathology. A detailed analysis of their expression and cellular localization in the healthy human kidney is presented, highlighting their involvement in fatty acid metabolism, electrolyte regulation, and the synthesis of important molecules. The review also delves into the diverse roles of SPs in various renal diseases, including renal ischemia/reperfusion injury, diabetic nephropathy, renal interstitial fibrosis, and lupus nephritis, elucidating their molecular mechanisms and potential as therapeutic targets. The review further discusses pharmacological modulation of SPs and its implications for treatment. Our findings provide a comprehensive understanding of SPs in renal health and disease, offering new avenues for targeted therapeutic interventions and precision medicine in nephrology.

A new pipeline SPICE identifies novel JUN-IKZF1 composite elements

Transcription factor partners can cooperatively bind to DNA composite elements to augment gene transcription. Here, we report a novel protein-DNA binding screening pipeline, termed Spacing Preference Identification of Composite Elements (SPICE), that can systematically predict protein binding partners and DNA motif spacing preferences. Using SPICE, we successfully identified known composite elements, such as AP1-IRF composite elements (AICEs) and STAT5 tetramers, and also uncovered several novel binding partners, including JUN-IKZF1 composite elements. One such novel interaction was identified at CNS9, an upstream conserved noncoding region in the human IL10 gene, which harbors a non-canonical IKZF1 binding site. We confirmed the cooperative binding of JUN and IKZF1 and showed that the activity of an IL10-luciferase reporter construct in primary B and T cells depended on both this site and the AP1 binding site within this composite element. Overall, our findings reveal an unappreciated global association of IKZF1 and AP1 and establish SPICE as a valuable new pipeline for predicting novel transcription binding complexes.

ChromatinHD connects single-cell DNA accessibility and conformation to gene expression through scale-adaptive machine learning

Gene regulation is inherently multiscale, but scale-adaptive machine learning methods that fully exploit this property in single-nucleus accessibility data are still lacking. Here, we develop ChromatinHD, a pair of scale-adaptive models that uses the raw accessibility data, without peak-calling or windows, to link regions to gene expression and determine differentially accessible chromatin. We show how ChromatinHD consistently outperforms existing peak and window-based approaches and find that this is due to a large number of uniquely captured, functional accessibility changes within and outside of putative cis-regulatory regions. Furthermore, ChromatinHD can delineate collaborating regulatory regions, including their preferential genomic conformations, that drive gene expression. Finally, our models also use changes in ATAC-seq fragment lengths to identify dense binding of transcription factors, a feature not captured by footprinting methods. Altogether, ChromatinHD, available at https://chromatinhd.org , is a suite of computational tools that enables a data-driven understanding of chromatin accessibility at various scales and how it relates to gene expression.

Rise and SINE: roles of transcription factors and retrotransposons in zygotic genome activation

In sexually reproducing organisms, life begins with the fusion of transcriptionally silent gametes, the oocyte and sperm. Although initiation of transcription in the embryo, known as zygotic genome activation (ZGA), is universally required for development, the transcription factors regulating this process are poorly conserved. In this Perspective, we discuss recent insights into the mechanisms of ZGA in totipotent mammalian embryos, namely ZGA regulation by several transcription factors, including by orphan nuclear receptors (OrphNRs) such as the pioneer transcription factor NR5A2, and by factors of the DUX, TPRX and OBOX families. We performed a meta-analysis and compiled a list of pan-ZGA genes, and found that most of these genes are indeed targets of the above transcription factors. Remarkably, more than a third of these ZGA genes appear to be regulated both by OrphNRs such as NR5A2 and by OBOX proteins, whose motifs co-occur in SINE B1 retrotransposable elements, which are enriched near ZGA genes. We propose that ZGA in mice is activated by recruitment of multiple transcription factors to SINE B1 elements that function as enhancers, and discuss a potential relevance of this mechanism to Alu retrotransposable elements in human ZGA.

Footprint-C reveals transcription factor modes in local clusters and long-range chromatin interactions

The proximity ligation-based Hi-C and derivative methods are the mainstream tools to study genome-wide chromatin interactions. These methods often fragment the genome using enzymes functionally irrelevant to the interactions per se, restraining the efficiency in identifying structural features and the underlying regulatory elements. Here we present Footprint-C, which yields high-resolution chromatin contact maps built upon intact and genuine footprints protected by transcription factor (TF) binding. When analyzed at one-dimensional level, the billions of chromatin contacts from Footprint-C enable genome-wide analysis at single footprint resolution, and reveal preferential modes of local TF co-occupancy. At pairwise contact level, Footprint-C exhibits higher efficiency in identifying chromatin structural features when compared with other Hi-C methods, segregates chromatin interactions emanating from adjacent TF footprints, and uncovers multiway interactions involving different TFs. Altogether, Footprint-C results suggest that rich regulatory modes of TF may underlie both local residence and distal chromatin interactions, in terms of TF identity, valency, and conformational configuration.

A new pipeline SPICE identifies novel JUN-IKZF1 composite elements

Transcription factor partners can cooperatively bind to DNA composite elements to augment gene transcription. Here, we report a novel protein-DNA binding screening pipeline, termed Spacing Preference Identification of Composite Elements (SPICE), that can systematically predict protein binding partners and DNA motif spacing preferences. Using SPICE, we successfully identified known composite elements, such as AP1-IRF composite elements (AICEs) and STAT5 tetramers, and also uncovered several novel binding partners, including JUN-IKZF1 composite elements. One such novel interaction was identified at CNS9, an upstream conserved noncoding region in the human IL10 gene, which harbors a non-canonical IKZF1 binding site. We confirmed cooperative binding of JUN and IKZF1 and showed that the activity of an IL10 -luciferase reporter construct in primary B and T cells depended on both this site and the AP1 binding site within this composite element. Overall, our findings reveal an unappreciated global association of IKZF1 and AP1 and establish SPICE as a valuable new pipeline for predicting novel transcription binding complexes.

Identification of co-localised transcription factors based on paired motifs analysis

The interaction of transcription factors (TFs) with DNA precisely regulates gene transcription. In mammalian cells, thousands of TFs often interact with DNA cis-regulatory elements in a combinatorial manner rather than act alone. The identification of cooperativity between TFs can help to explore the mechanism of transcriptional regulation. However, little is known about the cooperative patterns of TFs in the genome. To identify which TFs prefer co-localisation, the authors conducted a paired motif analysis in the accessible regions of the human genome based on the Poisson background model. Especially, the authors distinguish the cooperative binding TFs and the competitive binding TFs according to the distance between TF motifs. In the K562 cell line, the authors find that TFs from a same family are always competing the same binding sites, such as FOS_JUN family, whereas KLF family TFs show significant cooperative binding in the adjacency region. Furthermore, the comparative analysis across 16 human cell lines indicates that most TF combination patterns are conserved, but there are still some cell-line-specific patterns. Finally, in human prostate cancer cells (PC-3) and human prostate normal cells (RWPE-2), the authors investigate the specific TF combination patterns in the disease cell and normal cell. The results show that the cooperative binding TF pairs shared by PC-3 and RWPE-2 account for over 90%. Simultaneously, the authors also identify 26 specific TF combination pairs in PC-3 cancer cells.

Single cell variant to enhancer to gene map for coronary artery disease

Although genome wide association studies (GWAS) in large populations have identified hundreds of variants associated with common diseases such as coronary artery disease (CAD), most disease-associated variants lie within non-coding regions of the genome, rendering it difficult to determine the downstream causal gene and cell type. Here, we performed paired single nucleus gene expression and chromatin accessibility profiling from 44 human coronary arteries. To link disease variants to molecular traits, we developed a meta-map of 88 samples and discovered 11,182 single-cell chromatin accessibility quantitative trait loci (caQTLs). Heritability enrichment analysis and disease variant mapping demonstrated that smooth muscle cells (SMCs) harbor the greatest genetic risk for CAD. To capture the continuum of SMC cell states in disease, we used dynamic single cell caQTL modeling for the first time in tissue to uncover QTLs whose effects are modified by cell state and expand our insight into genetic regulation of heterogenous cell populations. Notably, we identified a variant in the COL4A1/COL4A2 CAD GWAS locus which becomes a caQTL as SMCs de-differentiate by changing a transcription factor binding site for EGR1/2. To unbiasedly prioritize functional candidate genes, we built a genome-wide single cell variant to enhancer to gene (scV2E2G) map for human CAD to link disease variants to causal genes in cell types. Using this approach, we found several hundred genes predicted to be linked to disease variants in different cell types. Next, we performed genome-wide Hi-C in 16 human coronary arteries to build tissue specific maps of chromatin conformation and link disease variants to integrated chromatin hubs and distal target genes. Using this approach, we show that rs4887091 within the ADAMTS7 CAD GWAS locus modulates function of a super chromatin interactome through a change in a CTCF binding site. Finally, we used CRISPR interference to validate a distal gene, AMOTL2, liked to a CAD GWAS locus. Collectively we provide a disease-agnostic framework to translate human genetic findings to identify pathologic cell states and genes driving disease, producing a comprehensive scV2E2G map with genetic and tissue level convergence for future mechanistic and therapeutic studies.

Evidence for compensatory evolution within pleiotropic regulatory elements

Pleiotropy, measured as expression breadth across tissues, is one of the best predictors for protein sequence and expression conservation. In this study, we investigated its effect on the evolution of cis-regulatory elements (CREs). To this end, we carefully reanalyzed the Epigenomics Roadmap data for nine fetal tissues, assigning a measure of pleiotropic degree to nearly half a million CREs. To assess the functional conservation of CREs, we generated ATAC-seq and RNA-seq data from humans and macaques. We found that more pleiotropic CREs exhibit greater conservation in accessibility, and the mRNA expression levels of the associated genes are more conserved. This trend of higher conservation for higher degrees of pleiotropy persists when analyzing the transcription factor binding repertoire. In contrast, simple DNA sequence conservation of orthologous sites between species tends to be even lower for pleiotropic CREs than for species-specific CREs. Combining various lines of evidence, we propose that the lack of sequence conservation in functionally conserved pleiotropic CREs is owing to within-element compensatory evolution. In summary, our findings suggest that pleiotropy is also a good predictor for the functional conservation of CREs, even though this is not reflected in the sequence conservation of pleiotropic CREs.

Identifying the "stripe" transcription factors and cooperative binding related to DNA methylation

DNA methylation plays a critical role in gene regulation by modulating the DNA binding of transcription factors (TFs). This study integrates TFs' ChIP-seq profiles with WGBS profiles to investigate how DNA methylation affects protein interactions. Statistical methods and a 5-letter DNA motif calling model have been developed to characterize DNA sequences bound by proteins, while considering the effects of DNA modifications. By employing these methods, 79 significant universal "stripe" TFs and cofactors (USFs), 2360 co-binding protein pairs, and distinct protein modules associated with various DNA methylation states have been identified. The USFs hint a regulatory hierarchy within these protein interactions. Proteins preferentially bind to non-CpG sites in methylated regions, indicating binding affinity is not solely CpG-dependent. Proteins involved in methylation-specific USFs and cobinding pairs play essential roles in promoting and sustaining DNA methylation through interacting with DNMTs or inhibiting TET binding. These findings underscore the interplay between protein binding and methylation, offering insights into epigenetic regulation in cellular biology.

Context transcription factors establish cooperative environments and mediate enhancer communication

Many enhancers control gene expression by assembling regulatory factor clusters, also referred to as condensates. This process is vital for facilitating enhancer communication and establishing cellular identity. However, how DNA sequence and transcription factor (TF) binding instruct the formation of high regulatory factor environments remains poorly understood. Here we developed a new approach leveraging enhancer-centric chromatin accessibility quantitative trait loci (caQTLs) to nominate regulatory factor clusters genome-wide. By analyzing TF-binding signatures within the context of caQTLs and comparing episomal versus endogenous enhancer activities, we discovered a class of regulators, 'context-only' TFs, that amplify the activity of cell type-specific caQTL-binding TFs, that is, 'context-initiator' TFs. Similar to super-enhancers, enhancers enriched for context-only TF-binding sites display high coactivator binding and sensitivity to bromodomain-inhibiting molecules. We further show that binding sites for context-only and context-initiator TFs underlie enhancer coordination, providing a mechanistic rationale for how a loose TF syntax confers regulatory specificity.

An Atlas of Accessible Chromatin in Advanced Prostate Cancer Reveals the Epigenetic Evolution during Tumor Progression

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal disease that resists therapy targeting androgen signaling, the primary driver of prostate cancer. mCRPC resists androgen receptor (AR) inhibitors by amplifying AR signaling or by evolving into therapy-resistant subtypes that do not depend on AR. Elucidation of the epigenetic underpinnings of these subtypes could provide important insights into the drivers of therapy resistance. In this study, we produced chromatin accessibility maps linked to the binding of lineage-specific transcription factors (TF) by performing assay for transposase-accessible chromatin sequencing on 70 mCRPC tissue biopsies integrated with transcriptome and whole-genome sequencing. mCRPC had a distinct global chromatin accessibility profile linked to AR function. Analysis of TF occupancy across accessible chromatin revealed 203 TFs associated with mCRPC subtypes. Notably, ZNF263 was identified as a putative prostate cancer TF with a significant impact on gene activity in the double-negative subtype (AR- neuroendocrine-), potentially activating MYC targets. Overall, this analysis of chromatin accessibility in mCRPC provides valuable insights into epigenetic changes that occur during progression to mCRPC. Significance: Integration of a large cohort of transcriptome, whole-genome, and ATAC sequencing characterizes the chromatin accessibility changes in advanced prostate cancer and identifies therapy-resistant prostate cancer subtype-specific transcription factors that modulate oncogenic programs.

Genomic binding of NF-Y in mouse and human cells

NF-Y is a Transcription Factor that regulates transcription through binding to the CCAAT-box. To understand its strategy, we analyzed 16 ChIP-seq datasets from human and mouse cells. Shared loci, mostly located in promoters of expressed genes of cell cycle, metabolism and gene expression pathways, are associated with histone marks of active chromatin and specific modules of TFs. Other peaks are in enhancers and Transposable Elements -TE- of retroviral origin in human and mouse. We evaluated the relationship with USF1, a common synergistic partner in promoters and MLT1 TEs, upon NF-YB inactivation: USF1 binding decreases in promoters, modestly in MLT1, suggesting a pioneering role of NF-Y in formers, not in the latters. These data define a common set of NF-Y functional targets across different mammalian cell types, suggesting a pioneering role in promoters with respect to TEs.

Mutations in linker-2 of KLF1 impair expression of membrane transporters and cytoskeletal proteins causing hemolysis

The SP/KLF family of transcription factors harbour three C-terminal C2H2 zinc fingers interspersed by two linkers which confers DNA-binding to a 9-10 bp motif. Mutations in KLF1, the founding member of the family, are common. Missense mutations in linker two result in a mild phenotype. However, when co-inherited with loss-of-function mutations, they result in severe non-spherocytic hemolytic anemia. We generate a mouse model of this disease by crossing Klf1+/- mice with Klf1H350R/+ mice that harbour a missense mutation in linker-2. Klf1H350R/- mice exhibit severe hemolysis without thalassemia. RNA-seq demonstrate loss of expression of genes encoding transmembrane and cytoskeletal proteins, but not globins. ChIP-seq show no change in DNA-binding specificity, but a global reduction in affinity, which is confirmed using recombinant proteins and in vitro binding assays. This study provides new insights into how linker mutations in zinc finger transcription factors result in different phenotypes to those caused by loss-of-function mutations.

Concerted regulation of skeletal muscle metabolism and contractile properties by the orphan nuclear receptor Nr2f6

BACKGROUND

The maintenance of skeletal muscle plasticity upon changes in the environment, nutrient supply, and exercise depends on regulatory mechanisms that couple structural and metabolic adaptations. The mechanisms that interconnect both processes at the transcriptional level remain underexplored. Nr2f6, a nuclear receptor, regulates metabolism and cell differentiation in peripheral tissues. However, its role in the skeletal muscle is still elusive. Here, we aimed to investigate the effects of Nr2f6 modulation on muscle biology in vivo and in vitro.

METHODS

Global RNA-seq was performed in Nr2f6 knockdown C2C12 myocytes (N = 4-5). Molecular and metabolic assays and proliferation experiments were performed using stable Nr2f6 knockdown and Nr2f6 overexpression C2C12 cell lines (N = 3-6). Nr2f6 content was evaluated in lipid overload models in vitro and in vivo (N = 3-6). In vivo experiments included Nr2f6 overexpression in mouse tibialis anterior muscle, followed by gene array transcriptomics and molecular assays (N = 4), ex vivo contractility experiments (N = 5), and histological analysis (N = 7). The conservation of Nr2f6 depletion effects was confirmed in primary skeletal muscle cells of humans and mice.

RESULTS

Nr2f6 knockdown upregulated genes associated with muscle differentiation, metabolism, and contraction, while cell cycle-related genes were downregulated. In human skeletal muscle cells, Nr2f6 knockdown significantly increased the expression of myosin heavy chain genes (two-fold to three-fold) and siRNA-mediated depletion of Nr2f6 increased maximal C2C12 myocyte's lipid oxidative capacity by 75% and protected against lipid-induced cell death. Nr2f6 content decreased by 40% in lipid-overloaded myotubes and by 50% in the skeletal muscle of mice fed a high-fat diet. Nr2f6 overexpression in mice resulted in an atrophic and hypoplastic state, characterized by a significant reduction in muscle mass (15%) and myofibre content (18%), followed by an impairment (50%) in force production. These functional phenotypes were accompanied by the establishment of an inflammation-like molecular signature and a decrease in the expression of genes involved in muscle contractility and oxidative metabolism, which was associated with the repression of the uncoupling protein 3 (20%) and PGC-1α (30%) promoters activity following Nr2f6 overexpression in vitro. Additionally, Nr2f6 regulated core components of the cell division machinery, effectively decoupling muscle cell proliferation from differentiation.

CONCLUSIONS

Our findings reveal a novel role for Nr2f6 as a molecular transducer that plays a crucial role in maintaining the balance between skeletal muscle contractile function and oxidative capacity. These results have significant implications for the development of potential therapeutic strategies for metabolic diseases and myopathies.

Rewiring gene circuitry for plant improvement

Aspirations for high crop growth and yield, nutritional quality and bioproduction of materials are challenged by climate change and limited adoption of new technologies. Here, we review recent advances in approaches to profile and model gene regulatory activity over developmental and response time in specific cells, which have revealed the basis of variation in plant phenotypes: both redeployment of key regulators to new contexts and their repurposing to control different slates of genes. New synthetic biology tools allow tunable, spatiotemporal regulation of transgenes, while recent gene-editing technologies enable manipulation of the regulation of native genes. Ultimately, understanding how gene circuitry is wired to control form and function across varied plant species, combined with advanced technology to rewire that circuitry, will unlock solutions to our greatest challenges in agriculture, energy and the environment.

MSC-derived exosomes deliver ZBTB4 to mediate transcriptional repression of ITIH3 in astrocytes in spinal cord injury

BACKGROUND

BMSC-secreted exosomes (BMSC-Exos) have shown potential for promoting behavioral recovery following spinal cord injury (SCI). However, its role in blocking astrocyte activation remains unclear. Thus, this study aimed to determine whether BMSC-Exos impair the function of astrocytes following SCI in mice and to seek the mechanism.

METHODS

BMSC-Exos were collected by ultracentrifugation and identified. The SCI mice were developed by laminectomy combined with spinal cord shock, followed by BMSC-Exos or nerve growth factor (positive control) treatment. HE staining, Nissl staining, and TUNEL were conducted to analyze the pathological structural damage and neuronal damage in the mouse spinal cord. Bioinformatics was used to screen altered molecules under the BMSC-Exos treatment. Effects of BMSC-Exos and changes in ZBTB4 and ITIH3 expression on neuronal damage induced by activated astrocytes in the co-culture system were analyzed by CCK-8 and flow cytometry.

RESULTS

Nerve growth factor and BMSC-Exos promoted motor function recovery, alleviated nerve injury, and reduced apoptosis in mice with SCI. ZBTB4 was enriched in BMSC-Exos and lowly expressed in SCI. Downregulation of ZBTB4 diminished the therapeutic effects of BMSC-Exos against SCI. ITIH3 was a downstream target of ZBTB4. Neurotoxic activation of astrocytes induced neuronal injury, which was alleviated by BMSC-Exos. However, ZBTB4 knockdown overturned the effects of BMSC-Exos in vitro and combined ITIH3 knockdown alleviated the accentuating effects of ZBTB4 knockdown on neuronal injury.

CONCLUSION

BMSC-Exos protected against astrocyte-induced neuronal injury by delivering ZBTB4 to repress ITIH3, ultimately improving motor function in mice with SCI.

Iterative deep learning-design of human enhancers exploits condensed sequence grammar to achieve cell type-specificity

An important and largely unsolved problem in synthetic biology is how to target gene expression to specific cell types. Here, we apply iterative deep learning to design synthetic enhancers with strong differential activity between two human cell lines. We initially train models on published datasets of enhancer activity and chromatin accessibility and use them to guide the design of synthetic enhancers that maximize predicted specificity. We experimentally validate these sequences, use the measurements to re-optimize the predictor, and design a second generation of enhancers with improved specificity. Our design methods embed relevant transcription factor binding site (TFBS) motifs with higher frequencies than comparable endogenous enhancers while using a more selective motif vocabulary, and we show that enhancer activity is correlated with transcription factor expression at the single cell level. Finally, we characterize causal features of top enhancers via perturbation experiments and show enhancers as short as 50bp can maintain specificity.

Ehmt2 inactivation in pancreatic epithelial cells shapes the transcriptional landscape and inflammation response of the whole pancreas

Introduction: The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that Ehmt2 inactivation in the mouse pancreas alters growth and immune gene expression networks, antagonizing Kras-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of Ehmt2 in maintaining a transcriptional landscape that protects organs from inflammation. Methods: Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from Ehmt2 conditional knockout animals (Ehmt2 fl/fl ) targeted to the exocrine pancreatic epithelial cells (Pdx1-Cre and P48 Cre/+ ), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the Ehmt2 fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels. Results and discussion: The Ehmt2 fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional Ehmt2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific Ehmt2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the Ehmt2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

Massively parallel characterization of regulatory elements in the developing human cortex

Nucleotide changes in gene regulatory elements are important determinants of neuronal development and diseases. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 open chromatin regions, including thousands of sequences with cell type-specific accessibility and variants associated with brain gene regulation. In primary cells, we identified 46,802 active enhancer sequences and 164 variants that alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.

CellWalker2: multi-omic discovery of hierarchical cell type relationships and their associations with genomic annotations

CellWalker2 is a graph diffusion-based method for single-cell genomics data integration. It extends the CellWalker model by incorporating hierarchical relationships between cell types, providing estimates of statistical significance, and adding data structures for analyzing multi-omics data so that gene expression and open chromatin can be jointly modeled. Our open-source software enables users to annotate cells using existing ontologies and to probabilistically match cell types between two or more contexts, including across species. CellWalker2 can also map genomic regions to cell ontologies, enabling precise annotation of elements derived from bulk data, such as enhancers, genetic variants, and sequence motifs. Through simulation studies, we show that CellWalker2 performs better than existing methods in cell type annotation and mapping. We then use data from the brain and immune system to demonstrate CellWalker2's ability to discover cell type-specific regulatory programs and both conserved and divergent cell type relationships in complex tissues.

Type I interferons induce an epigenetically distinct memory B cell subset in chronic viral infection

Memory B cells (MBCs) are key providers of long-lived immunity against infectious disease, yet in chronic viral infection, they do not produce effective protection. How chronic viral infection disrupts MBC development and whether such changes are reversible remain unknown. Through single-cell (sc)ATAC-seq and scRNA-seq during acute versus chronic lymphocytic choriomeningitis viral infection, we identified a memory subset enriched for interferon (IFN)-stimulated genes (ISGs) during chronic infection that was distinct from the T-bet+ subset normally associated with chronic infection. Blockade of IFNAR-1 early in infection transformed the chromatin landscape of chronic MBCs, decreasing accessibility at ISG-inducing transcription factor binding motifs and inducing phenotypic changes in the dominating MBC subset, with a decrease in the ISG subset and an increase in CD11c+CD80+ cells. However, timing was critical, with MBCs resistant to intervention at 4 weeks post-infection. Together, our research identifies a key mechanism to instruct MBC identity during viral infection.

DNA methylation regulates RNA m6A modification through transcription factor SP1 during the development of porcine somatic cell nuclear transfer embryos

Epigenetic modifications play critical roles during somatic cell nuclear transfer (SCNT) embryo development. Whether RNA N6-methyladenosine (m6A) affects the developmental competency of SCNT embryos remains unclear. Here, we showed that porcine bone marrow mesenchymal stem cells (pBMSCs) presented higher RNA m6A levels than those of porcine embryonic fibroblasts (pEFs). SCNT embryos derived from pBMSCs had higher RNA m6A levels, cleavage, and blastocyst rates than those from pEFs. Compared with pEFs, the promoter region of METTL14 presented a hypomethylation status in pBMSCs. Mechanistically, DNA methylation regulated METTL14 expression by affecting the accessibility of transcription factor SP1 binding, highlighting the role of the DNA methylation/SP1/METTL14 pathway in donor cells. Inhibiting the DNA methylation level in donor cells increased the RNA m6A level and improved the development efficiency of SCNT embryos. Overexpression of METTL14 significantly increased the RNA m6A level in donor cells and the development efficiency of SCNT embryos, whereas knockdown of METTL14 suggested the opposite result. Moreover, we revealed that RNA m6A-regulated TOP2B mRNA stability, translation level, and DNA damage during SCNT embryo development. Collectively, our results highlight the crosstalk between RNA m6A and DNA methylation, and the crucial role of RNA m6A during nuclear reprogramming in SCNT embryo development.

EHMT2 Inactivation in Pancreatic Epithelial Cells Shapes the Transcriptional Landscape and Inflammation Response of the Whole Pancreas

The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that EHMT2 inactivation alters growth and immune gene expression networks, antagonizing KRAS-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of EHMT2 in maintaining a transcriptional landscape that protects organs from inflammation. Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from EHMT2 conditional knockout animals ( EHMT2 fl/fl ) targeted to the exocrine pancreatic epithelial cells ( Pdx1-Cre and P48 Cre/+ ), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the EHMT2 fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels. The EHMT2 fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional EHMT2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific EHMT2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the EHMT2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

A transcription factor (TF) inference method that broadly measures TF activity and identifies mechanistically distinct TF networks

TF profiler is a method of inferring transcription factor regulatory activity, i.e. when a TF is present and actively regulating transcription, directly directly from nascent sequencing assays such as PRO-seq and GRO-seq. Transcription factors orchestrate transcription and play a critical role in cellular maintenance, identity and response to external stimuli. While ChIP assays have measured DNA localization, they fall short of identifying when and where transcription factors are actively regulating transcription. Our method, on the other hand, uses RNA polymerase activity to infer TF activity across hundreds of data sets and transcription factors. Based on these classifications we identify three distinct classes of transcription factors: ubiquitous factors that play roles in cellular homeostasis, driving basal gene programs across tissues and cell types, tissue specific factors that act almost exclusively at enhancers and are themselves regulated at transcription, and stimulus responsive TFs which are regulated post-transcriptionally but act predominantly at enhancers. TF profiler is broadly applicable, providing regulatory insights on any PRO-seq sample for any transcription factor with a known binding motif.

Chemical-genetics refines transcription factor regulatory circuits

Transcriptional dysregulation is a key step in oncogenesis, but our understanding of transcriptional control has relied on genetic approaches that are slow and allow for compensation. Chemical-genetic approaches have shortened the time frame for the analysis of transcription factors from days or weeks to minutes. These studies show that while DNA-binding proteins bind to thousands of sites, they are directly required to regulate only a small cadre of genes. Moreover, these transcriptional control networks are far more distinct, with much less overlap and interconnectivity than predicted from DNA binding. The identified direct targets can then be used to dissect the mechanism of action of these factors, which could identify ways to therapeutically manipulate these oncogenic transcriptional control networks.

Emerging insights into enhancer biology and function

Cell type-specific gene expression is coordinated by DNA-encoded enhancers and the transcription factors (TFs) that bind to them in a sequence-specific manner. As such, these enhancers and TFs are critical mediators of normal development and altered enhancer or TF function is associated with the development of diseases such as cancer. While initially defined by their ability to activate gene transcription in reporter assays, putative enhancer elements are now frequently defined by their unique chromatin features including DNase hypersensitivity and transposase accessibility, bidirectional enhancer RNA (eRNA) transcription, CpG hypomethylation, high H3K27ac and H3K4me1, sequence-specific transcription factor binding, and co-factor recruitment. Identification of these chromatin features through sequencing-based assays has revolutionized our ability to identify enhancer elements on a genome-wide scale, and genome-wide functional assays are now capitalizing on this information to greatly expand our understanding of how enhancers function to provide spatiotemporal coordination of gene expression programs. Here, we highlight recent technological advances that are providing new insights into the molecular mechanisms by which these critical cis-regulatory elements function in gene control. We pay particular attention to advances in our understanding of enhancer transcription, enhancer-promoter syntax, 3D organization and biomolecular condensates, transcription factor and co-factor dependencies, and the development of genome-wide functional enhancer screens.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

Physical interaction between MSL2 and CLAMP assures direct cooperativity and prevents competition at composite binding sites

MSL2, the DNA-binding subunit of the Drosophila dosage compensation complex, cooperates with the ubiquitous protein CLAMP to bind MSL recognition elements (MREs) on the X chromosome. We explore the nature of the cooperative binding to these GA-rich, composite sequence elements in reconstituted naïve embryonic chromatin. We found that the cooperativity requires physical interaction between both proteins. Remarkably, disruption of this interaction does not lead to indirect, nucleosome-mediated cooperativity as expected, but to competition. The protein interaction apparently not only increases the affinity for composite binding sites, but also locks both proteins in a defined dimeric state that prevents competition. High Affinity Sites of MSL2 on the X chromosome contain variable numbers of MREs. We find that the cooperation between MSL2/CLAMP is not influenced by MRE clustering or arrangement, but happens largely at the level of individual MREs. The sites where MSL2/CLAMP bind strongly in vitro locate to all chromosomes and show little overlap to an expanded set of X-chromosomal MSL2 in vivo binding sites generated by CUT&RUN. Apparently, the intrinsic MSL2/CLAMP cooperativity is limited to a small selection of potential sites in vivo. This restriction must be due to components missing in our reconstitution, such as roX2 lncRNA.

Massively parallel characterization of transcriptional regulatory elements in three diverse human cell types

The human genome contains millions of candidate cis-regulatory elements (CREs) with cell-type-specific activities that shape both health and myriad disease states. However, we lack a functional understanding of the sequence features that control the activity and cell-type-specific features of these CREs. Here, we used lentivirus-based massively parallel reporter assays (lentiMPRAs) to test the regulatory activity of over 680,000 sequences, representing a nearly comprehensive set of all annotated CREs among three cell types (HepG2, K562, and WTC11), finding 41.7% to be functional. By testing sequences in both orientations, we find promoters to have significant strand orientation effects. We also observe that their 200 nucleotide cores function as non-cell-type-specific 'on switches' providing similar expression levels to their associated gene. In contrast, enhancers have weaker orientation effects, but increased tissue-specific characteristics. Utilizing our lentiMPRA data, we develop sequence-based models to predict CRE function with high accuracy and delineate regulatory motifs. Testing an additional lentiMPRA library encompassing 60,000 CREs in all three cell types, we further identified factors that determine cell-type specificity. Collectively, our work provides an exhaustive catalog of functional CREs in three widely used cell lines, and showcases how large-scale functional measurements can be used to dissect regulatory grammar.

Massively parallel characterization of psychiatric disorder-associated and cell-type-specific regulatory elements in the developing human cortex

Nucleotide changes in gene regulatory elements are important determinants of neuronal development and disease. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 sequences, including differentially accessible cell-type specific regions in the developing cortex and single-nucleotide variants associated with psychiatric disorders. In primary cells, we identified 46,802 active enhancer sequences and 164 disorder-associated variants that significantly alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning, we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.

EKLF/Klf1 regulates erythroid transcription by its pioneering activity and selective control of RNA Pol II pause-release

EKLF/Klf1 is a zinc-finger transcription activator essential for erythroid lineage commitment and terminal differentiation. Using ChIP-seq, we investigate EKLF DNA binding and transcription activation mechanisms during mouse embryonic erythropoiesis. We utilize the Nan/+ mouse that expresses the EKLF-E339D (Nan) variant mutated in its conserved zinc-finger region and address the mechanism of hypomorphic and neomorphic changes in downstream gene expression. First, we show that Nan-EKLF limits normal EKLF binding to a subset of its sites. Second, we find that ectopic binding of Nan-EKLF occurs largely at enhancers and activates transcription through pioneering activity. Third, we find that for a subset of ectopic targets, gene activation is achieved in Nan/+ only by Nan-EKLF binding to distal enhancers, leading to RNA polymerase II pause-release. These results have general applicability to understanding how a DNA binding variant factor confers dominant disruptive effects on downstream gene expression even in the presence of its normal counterpart.

TFSyntax: a database of transcription factors binding syntax in mammalian genomes

In mammals, transcriptional factors (TFs) drive gene expression by binding to regulatory elements in a cooperative manner. Deciphering the rules of such cooperation is crucial to obtain a full understanding of cellular homeostasis and development. Although this is a long-standing topic, there is no comprehensive database for biologists to access the syntax of TF binding sites. Here we present TFSyntax (https://tfsyntax.zhaopage.com), a database focusing on the arrangement of TF binding sites. TFSyntax maps the binding motif of 1299 human TFs and 890 mouse TFs across 382 cells and tissues, representing the most comprehensive TF binding map to date. In addition to location, TFSyntax defines motif positional preference, density and colocalization within accessible elements. Powered by a series of functional modules based on web interface, users can freely search, browse, analyze, and download data of interest. With comprehensive characterization of TF binding syntax across distinct tissues and cell types, TFSyntax represents a valuable resource and platform for studying the mechanism of transcriptional regulation and exploring how regulatory DNA variants cause disease.