Shared nucleotide flanks confer transcriptional competency to bZip core motifs

The activity of early-life gene regulatory elements is hijacked in aging through pervasive AP-1-linked chromatin opening

A mechanistic connection between aging and development is largely unexplored. Through profiling age-related chromatin and transcriptional changes across 22 murine cell types, analyzed alongside previous mouse and human organismal maturation datasets, we uncovered a transcription factor binding site (TFBS) signature common to both processes. Early-life candidate cis-regulatory elements (cCREs), progressively losing accessibility during maturation and aging, are enriched for cell-type identity TFBSs. Conversely, cCREs gaining accessibility throughout life have a lower abundance of cell identity TFBSs but elevated activator protein 1 (AP-1) levels. We implicate TF redistribution toward these AP-1 TFBS-rich cCREs, in synergy with mild downregulation of cell identity TFs, as driving early-life cCRE accessibility loss and altering developmental and metabolic gene expression. Such remodeling can be triggered by elevating AP-1 or depleting repressive H3K27me3. We propose that AP-1-linked chromatin opening drives organismal maturation by disrupting cell identity TFBS-rich cCREs, thereby reprogramming transcriptome and cell function, a mechanism hijacked in aging through ongoing chromatin opening.

Lineage-determining transcription factor-driven promoters regulate cell type-specific macrophage gene expression

Mammalian promoters consist of multifarious elements, which make them unique and support the selection of the proper transcript variants required under diverse conditions in distinct cell types. However, their direct DNA-transcription factor (TF) interactions are mostly unidentified. Murine bone marrow-derived macrophages (BMDMs) are a widely used model for studying gene expression regulation. Thus, this model serves as a rich source of various next-generation sequencing data sets, including a large number of TF cistromes. By processing and integrating the available cistromic, epigenomic and transcriptomic data from BMDMs, we characterized the macrophage-specific direct DNA-TF interactions, with a particular emphasis on those specific for promoters. Whilst active promoters are enriched for certain types of typically methylatable elements, more than half of them contain non-methylatable and prototypically promoter-distal elements. In addition, circa 14% of promoters-including that of Csf1r-are composed exclusively of 'distal' elements that provide cell type-specific gene regulation by specialized TFs. Similar to CG-rich promoters, these also contain methylatable CG sites that are demethylated in a significant portion and show high polymerase activity. We conclude that this unusual class of promoters regulates cell type-specific gene expression in macrophages, and such a mechanism might exist in other cell types too.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD)

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD)

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Overlapping binding sites underlie TF genomic occupancy

Sequence-specific DNA binding by transcription factors (TFs) is a crucial step in gene regulation. However, current high-throughput in vitro approaches cannot reliably detect lower affinity TF-DNA interactions, which play key roles in gene regulation. Here, we developed PADIT-seq ( p rotein a ffinity to D NA by in vitro transcription and RNA seq uencing) to assay TF binding preferences to all 10-bp DNA sequences at far greater sensitivity than prior approaches. The expanded catalogs of low affinity DNA binding sites for the human TFs HOXD13 and EGR1 revealed that nucleotides flanking high affinity DNA binding sites create overlapping lower affinity sites that together modulate TF genomic occupancy in vivo . Formation of such extended recognition sequences stems from an inherent property of TF binding sites to interweave each other and expands the genomic sequence space for identifying noncoding variants that directly alter TF binding.

One-Sentence Summary

Overlapping DNA binding sites underlie TF genomic occupancy through their inherent propensity to interweave each other.

MAE-seq refines regulatory elements across the genome

Proper cell fate determination relies on precise spatial and temporal genome-wide cooperation between regulatory elements (REs) and their targeted genes. However, the lengths of REs defined using different methods vary, which indicates that there is sequence redundancy and that the context of the genome may be unintelligible. We developed a method called MAE-seq (Massive Active Enhancers by Sequencing) to experimentally identify functional REs at a 25-bp scale. In this study, MAE-seq was used to identify 626879, 541617 and 554826 25-bp enhancers in mouse embryonic stem cells (mESCs), C2C12 and HEK 293T, respectively. Using ∼1.6 trillion 25 bp DNA fragments and screening 12 billion cells, we identified 626879 as active enhancers in mESCs as an example. Comparative analysis revealed that most of the histone modification datasets were annotated by MAE-Seq loci. Furthermore, 33.85% (212195) of the identified enhancers were identified as de novo ones with no epigenetic modification. Intriguingly, distinct chromatin states dictate the requirement for dissimilar cofactors in governing novel and known enhancers. Validation results show that these 25-bp sequences could act as a functional unit, which shows identical or similar expression patterns as the previously defined larger elements, Enhanced resolution facilitated the identification of numerous cell-specific enhancers and their accurate annotation as super enhancers. Moreover, we characterized novel elements capable of augmenting gene activity. By integrating with high-resolution Hi-C data, over 55.64% of novel elements may have a distal association with different targeted genes. For example, we found that the Cdh1 gene interacts with one novel and two known REs in mESCs. The biological effects of these interactions were investigated using CRISPR-Cas9, revealing their role in coordinating Cdh1 gene expression and mESC proliferation. Our study presents an experimental approach to refine the REs at 25-bp resolution, advancing the precision of genome annotation and unveiling the underlying genome context. This novel approach not only advances our understanding of gene regulation but also opens avenues for comprehensive exploration of the genomic landscape.

The Rax homeoprotein in Müller glial cells is required for homeostasis maintenance of the postnatal mouse retina

Müller glial cells, which are the most predominant glial subtype in the retina, play multiple important roles, including the maintenance of structural integrity, homeostasis, and physiological functions of the retina. We have previously found that the Rax homeoprotein is expressed in postnatal and mature Müller glial cells in the mouse retina. However, the function of Rax in postnatal and mature Müller glial cells remains to be elucidated. In the current study, we first investigated Rax function in retinal development using retroviral lineage analysis and found that Rax controls the specification of late-born retinal cell types, including Müller glial cells in the postnatal retina. We next generated Rax tamoxifen-induced conditional KO (Rax iCKO) mice, where Rax can be depleted in mTFP-labeled Müller glial cells upon tamoxifen treatment, by crossing Raxflox/flox mice with Rlbp1-CreERT2 mice, which we have produced. Immunohistochemical analysis showed a characteristic of reactive gliosis and enhanced gliosis of Müller glial cells in Rax iCKO retinas under normal and stress conditions, respectively. We performed RNA-seq analysis on mTFP-positive cells purified from the Rax iCKO retina and found significantly reduced expression of suppressor of cytokinesignaling-3 (Socs3). Reporter gene assays showed that Rax directly transactivates the Socs3 promoter. We observed decreased expression of Socs3 in Müller glial cells of Rax iCKO retinas by immunostaining. Taken together, the present results suggest that Rax suppresses inflammation in Müller glial cells by transactivating Socs3. This study sheds light on the transcriptional regulatory mechanisms underlying retinal Müller glial cell homeostasis.

Genome-wide census of ATF4 binding sites and functional profiling of trait-associated genetic variants overlapping ATF4 binding motifs

Activating Transcription Factor 4 (ATF4) is an important regulator of gene expression in stress responses and developmental processes in many cell types. Here, we catalogued ATF4 binding sites in the human genome and identified overlaps with trait-associated genetic variants. We probed these genetic variants for allelic regulatory activity using a massively parallel reporter assay (MPRA) in HepG2 hepatoma cells exposed to tunicamycin to induce endoplasmic reticulum stress and ATF4 upregulation. The results revealed that in the majority of cases, the MPRA allelic activity of these SNPs was in agreement with the nucleotide preference seen in the ATF4 binding motif from ChIP-Seq. Luciferase and electrophoretic mobility shift assays in additional cellular models further confirmed ATF4-dependent regulatory effects for the SNPs rs532446 (GADD45A intronic; linked to hematological parameters), rs7011846 (LPL upstream; myocardial infarction), rs2718215 (diastolic blood pressure), rs281758 (psychiatric disorders) and rs6491544 (educational attainment). CRISPR-Cas9 disruption and/or deletion of the regulatory elements harboring rs532446 and rs7011846 led to the downregulation of GADD45A and LPL, respectively. Thus, these SNPs could represent examples of GWAS genetic variants that affect gene expression by altering ATF4-mediated transcriptional activation.

Structural basis for cell type specific DNA binding of C/EBPβ: The case of cell cycle inhibitor p15INK4b promoter

C/EBPβ is a key regulator of numerous cellular processes, but it can also contribute to tumorigenesis and viral diseases. It binds to specific DNA sequences (C/EBP sites) and interacts with other transcription factors to control expression of multiple eukaryotic genes in a tissue and cell-type dependent manner. A body of evidence has established that cell-type-specific regulatory information is contained in the local DNA sequence of the binding motif. In human epithelial cells, C/EBPβ is an essential cofactor for TGFβ signaling in the case of Smad2/3/4 and FoxO-dependent induction of the cell cycle inhibitor, p15INK4b. In the TGFβ-responsive region 2 of the p15INK4b promoter, the Smad binding site is flanked by a C/EBP site, CTTAA•GAAAG, which differs from the canonical, palindromic ATTGC•GCAAT motif. The X-ray crystal structure of C/EBPβ bound to the p15INK4b promoter fragment shows how GCGC-to-AAGA substitution generates changes in the intermolecular interactions in the protein-DNA interface that enhances C/EBPβ binding specificity, limits possible epigenetic regulation of the promoter, and generates a DNA element with a unique pattern of methyl groups in the major groove. Significantly, CT/GA dinucleotides located at the 5'ends of the double stranded element maintain local narrowing of the DNA minor groove width that is necessary for DNA recognition. Our results suggest that C/EBPβ would accept all forms of modified cytosine in the context of the CpT site. This contrasts with the effect on the consensus motif, where C/EBPβ binding is modestly increased by cytosine methylation, but substantially decreased by hydroxymethylation.

Flexibility of flanking DNA is a key determinant of transcription factor affinity for the core motif

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code.

Sequence-specific dynamics of DNA response elements and their flanking sites regulate the recognition by AP-1 transcription factors

Activator proteins 1 (AP-1) comprise one of the largest families of eukaryotic basic leucine zipper transcription factors. Despite advances in the characterization of AP-1 DNA-binding sites, our ability to predict new binding sites and explain how the proteins achieve different gene expression levels remains limited. Here we address the role of sequence-specific DNA flexibility for stability and specific binding of AP-1 factors, using microsecond-long molecular dynamics simulations. As a model system, we employ yeast AP-1 factor Yap1 binding to three different response elements from two genetic environments. Our data show that Yap1 actively exploits the sequence-specific flexibility of DNA within the response element to form stable protein–DNA complexes. The stability also depends on the four to six flanking nucleotides, adjacent to the response elements. The flanking sequences modulate the conformational adaptability of the response element, making it more shape-efficient to form specific contacts with the protein. Bioinformatics analysis of differential expression of the studied genes supports our conclusions: the stability of Yap1–DNA complexes, modulated by the flanking environment, influences the gene expression levels. Our results provide new insights into mechanisms of protein–DNA recognition and the biological regulation of gene expression levels in eukaryotes.

Motif grammar: The basis of the language of gene expression

Collaboration of transcription factors (TFs) and their recognition motifs in DNA is the result of coevolution and forms the basis of gene regulation. However, the way how these short genomic sequences contribute to setting the level of gene products is not understood in sufficient detail. The biological problem to be solved by the cell is complex, because each gene requires a unique regulatory network in each cellular condition using the same genome. Thus far, only some components of these networks have been uncovered. In this review, we compiled the features and principles of the motif grammar, which dictates the characteristics and thus the likelihood of the interactions of the binding TFs and their coregulators. We present how sequence features provide specificity using, as examples, two major TF superfamilies, the bZIP proteins and nuclear receptors. We also discuss the phenomenon of “weak” (low affinity) binding sites, which appear to be components of several important genomic regulatory regions, but paradoxically are barely detectable by the currently used approaches. Assembling the complete set of regulatory regions composed of both weak and strong binding sites will allow one to get more comprehensive lists of factors playing roles in gene regulation, thus making possible the deeper understanding of regulatory networks.

TFregulomeR reveals transcription factors' context-specific features and functions

Transcription factors (TFs) are sequence-specific DNA binding proteins, fine-tuning spatiotemporal gene expression. Since genomic occupancy of a TF is highly dynamic, it is crucial to study TF binding sites (TFBSs) in a cell-specific context. To date, thousands of ChIP-seq datasets have portrayed the genomic binding landscapes of numerous TFs in different cell types. Although these datasets can be browsed via several platforms, tools that can operate on that data flow are still lacking. Here, we introduce TFregulomeR (https://github.com/benoukraflab/TFregulomeR), an R-library linked to an up-to-date compendium of cistrome and methylome datasets, implemented with functionalities that facilitate integrative analyses. In particular, TFregulomeR enables the characterization of TF binding partners and cell-specific TFBSs, along with the study of TF’s functions in the context of different partnerships and DNA methylation levels. We demonstrated that TFs’ target gene ontologies can differ notably depending on their partners and, by re-analyzing well characterized TFs, we brought to light that numerous leucine zipper TFBSs derived from ChIP-seq experiments documented in current databases were inadequately characterized, due to the fact that their position weight matrices were assembled using a mixture of homodimer and heterodimer binding sites. Altogether, analyses of context-specific transcription regulation with TFregulomeR foster our understanding of regulatory network-dependent TF functions.

A novel upstream transcription factor 1 target gene N4bp2l1 that regulates adipogenesis

N4BP2l1, which is highly expressed in oral squamous cell carcinoma, is associated with poor prognosis. However, N4bp2l1's role in adipogenesis remains unknown. We aimed to clarify the expression profile and transcriptional regulation of N4bp2l1 to elucidate the functions underlying the role of N4bp2l1 in adipocyte differentiation. Our results revealed that N4bp2l1 mRNA expression increased in 3T3-L1 cells in a differentiation-dependent manner. To investigate the transcriptional regulation of N4bp2l1, the 2-kb 5′ region upstream of the mouse N4bp2l1 promoter was cloned into a luciferase vector. Luciferase reporter assays indicated that USF1 induces the N4bp2l1 promoter activity. Electrophoretic mobility shift and chromatin immunoprecipitation assays confirmed that USF1 directly binds to the Ebox in the N4bp2l1 promoter. Furthermore, the expressions of adipocyte differentiation markers significantly decreased in N4bp2l1-knockdown cells compared with those in control cells. Our results demonstrated that N4bp2l1 is a novel USF1 target gene that may be involved in adipogenesis regulation.

•

N4bp2l1 expression is increased in a differentiation-dependent manner in 3T3-L1.

•

N4bp2l1 is a novel USF1 target gene.

•

USF1 directly binds to the Ebox in the N4bp2l1 promoter.

•

Inhibition of 3T3-L1 adipocyte differentiation by N4bp2l1 knockdown.