Transcription Factor Chromatin Research Articles

Gene expression, transcription factor binding and histone modification predict leaf adaxial–abaxial polarity related genes

Leaf adaxial–abaxial (ad–abaxial) polarity is crucial for leaf morphology and function, but the genetic machinery governing this process remains unclear. To uncover critical genes involved in leaf ad–abaxial patterning, we applied a combination of in silico prediction using machine learning (ML) and experimental analysis. A Random Forest model was trained using genes known to influence ad–abaxial polarity as ground truth. Gene expression data from various tissues and conditions as well as promoter regulation data derived from transcription factor chromatin immunoprecipitation sequencing (ChIP-seq) was used as input, enabling the prediction of novel ad–abaxial polarity-related genes and additional transcription factors. Parallel to this, available and newly-obtained transcriptome data enabled us to identify genes differentially expressed across leaf ad–abaxial sides. Based on these analyses, we obtained a set of 111 novel genes which are involved in leaf ad–abaxial specialization. To explore implications for vegetable crop breeding, we examined the conservation of expression patterns between Arabidopsis and Brassica rapa using single-cell transcriptomics. The results demonstrated the utility of our computational approach for predicting candidate genes in crop species. Our findings expand the understanding of the genetic networks governing leaf ad–abaxial differentiation in agriculturally important vegetables, enhancing comprehension of natural variation impacting leaf morphology and development, with demonstrable breeding applications.

Development of an orally bioavailable mSWI/SNF ATPase degrader and acquired mechanisms of resistance in prostate cancer

Mammalian switch/sucrose nonfermentable (mSWI/SNF) ATPase degraders have been shown to be effective in enhancer-driven cancers by functioning to impede oncogenic transcription factor chromatin accessibility. Here, we developed AU-24118, an orally bioavailable proteolysis-targeting chimera (PROTAC) degrader of mSWI/SNF ATPases (SMARCA2 and SMARCA4) and PBRM1. AU-24118 demonstrated tumor regression in a model of castration-resistant prostate cancer (CRPC) which was further enhanced with combination enzalutamide treatment, a standard of care androgen receptor (AR) antagonist used in CRPC patients. Importantly, AU-24118 exhibited favorable pharmacokinetic profiles in preclinical analyses in mice and rats, and further toxicity testing in mice showed a favorable safety profile. As acquired resistance is common with targeted cancer therapeutics, experiments were designed to explore potential mechanisms of resistance that may arise with long-term mSWI/SNF ATPase PROTAC treatment. Prostate cancer cell lines exposed to long-term treatment with high doses of a mSWI/SNF ATPase degrader developed SMARCA4 bromodomain mutations and ABCB1 (ATP binding cassette subfamily B member 1) overexpression as acquired mechanisms of resistance. Intriguingly, while SMARCA4 mutations provided specific resistance to mSWI/SNF degraders, ABCB1 overexpression provided broader resistance to other potent PROTAC degraders targeting bromodomain-containing protein 4 and AR. The ABCB1 inhibitor, zosuquidar, reversed resistance to all three PROTAC degraders tested. Combined, these findings position mSWI/SNF degraders for clinical translation for patients with enhancer-driven cancers and define strategies to overcome resistance mechanisms that may arise.

Old school, new rules: floral meristem development revealed by 3D gene expression atlases and high-resolution transcription factor-chromatin dynamics.

The intricate morphology of the flower is primarily established within floral meristems in which floral organs will be defined and from where the developing flower will emerge. Floral meristem development involves multiscale-level regulation, including lineage and positional mechanisms for establishing cell-type identity, and transcriptional regulation mediated by changes in the chromatin environment. However, many key aspects of floral meristem development remain to be determined, such as: 1) the exact role of cellular location in connecting transcriptional inputs to morphological outcomes, and 2) the precise interactions between transcription factors and chromatin regulators underlying the transcriptional networks that regulate the transition from cell proliferation to differentiation during floral meristem development. Here, we highlight recent studies addressing these points through newly developed spatial reconstruction techniques and high-resolution transcription factor-chromatin environment interactions in the model plant Arabidopsis thaliana. Specifically, we feature studies that reconstructed 3D gene expression atlases of the floral meristem. We also discuss how the precise timing of floral meristem specification, floral organ patterning, and floral meristem termination is determined through temporally defined epigenetic dynamics for fine-tuning of gene expression. These studies offer fresh insights into the well-established principles of floral meristem development and outline the potential for further advances in this field in an age of integrated, powerful, multiscale resolution approaches.

The BAF chromatin remodeler synergizes with RNA polymerase II and transcription factors to evict nucleosomes

Chromatin accessibility is a hallmark of active transcription and entails ATP-dependent nucleosome remodeling, which is carried out by complexes such as Brahma-associated factor (BAF). However, the mechanistic links between transcription, nucleosome remodeling and chromatin accessibility are unclear. Here, we used a chemical–genetic approach coupled with time-resolved chromatin profiling to dissect the interplay between RNA Polymerase II (RNAPII), BAF and DNA-sequence-specific transcription factors in mouse embryonic stem cells. We show that BAF dynamically unwraps and evicts nucleosomes at accessible chromatin regions, while RNAPII promoter-proximal pausing stabilizes BAF chromatin occupancy and enhances ATP-dependent nucleosome eviction by BAF. We find that although RNAPII and BAF dynamically probe both transcriptionally active and Polycomb-repressed genomic regions, pluripotency transcription factor chromatin binding confers locus specificity for productive chromatin remodeling and nucleosome eviction by BAF. Our study suggests a paradigm for how functional synergy between dynamically acting chromatin factors regulates locus-specific nucleosome organization and chromatin accessibility.

Promoter R-Loops Recruit U2AF1 to Modulate Its Phase Separation and RNA Splicing.

R-loops and guanine quadruplexes (G4s) are secondary structures of nucleic acids that are ubiquitously present in cells and are enriched in promoter regions of genes. By employing a bioinformatic approach based on overlap analysis of transcription factor chromatin immunoprecipitation sequencing (ChIP-seq) data sets, we found that many splicing factors, including U2AF1 whose recognition of the 3' splicing site is crucial for pre-mRNA splicing, exhibit pronounced enrichment at endogenous R-loop- and DNA G4-structure loci in promoter regions of human genes. We also revealed that U2AF1 binds directly to R-loops and DNA G4 structures at a low-nM binding affinity. Additionally, we showed the ability of U2AF1 to undergo phase separation, which could be stimulated by binding with R-loops, but not duplex DNA, RNA/DNA hybrid, DNA G4, or single-stranded RNA. We also demonstrated that U2AF1 binds to promoter R-loops in human cells, and this binding competes with U2AF1's interaction with 3' splicing site and leads to augmented distribution of RNA polymerase II (RNAPII) to promoters over gene bodies, thereby modulating cotranscriptional pre-mRNA splicing. Together, we uncovered a group of candidate proteins that can bind to both R-loops and DNA G4s, revealed the direct and strong interactions of U2AF1 with these nucleic acid structures, and established a biochemical rationale for U2AF1's occupancy in gene promoters. We also unveiled that interaction with R-loops promotes U2AF1's phase separation, and our work suggests that U2AF1 modulates pre-mRNA splicing by regulating RNAPII's partition in transcription initiation versus elongation.

Function and regulation of transcription factors during mitosis-to-G1 transition.

During cell division, drastic cellular changes characteristic of mitosis result in the inactivation of the transcriptional machinery, and global downregulation of transcription. Sequence-specific transcription factors (TFs) have thus been considered mere bystanders, devoid of any regulatory function during mitosis. This view changed significantly in recent years, upon the conclusion that many TFs associate with condensed chromosomes during cell division, even occupying a fraction of their genomic target sites in mitotic chromatin. This finding was at the origin of the concept of mitotic bookmarking by TFs, proposed as a mechanism to propagate gene regulatory information across cell divisions, by facilitating the reactivation of specific bookmarked genes. While the underlying mechanisms and biological significance of this model remain elusive, recent developments in this fast-moving field have cast new light into TF activity during mitosis, beyond a bookmarking role. Here, we start by reviewing the most recent findings on the complex nature of TF-chromatin interactions during mitosis, and on mechanisms that may regulate them. Next, and in light of recent reports describing how transcription is reinitiated in temporally distinct waves during mitosis-to-G1 transition, we explore how TFs may contribute to defining this hierarchical gene expression process. Finally, we discuss how TF activity during mitotic exit may impact the acquisition of cell identity upon cell division, and propose a model that integrates dynamic changes in TF-chromatin interactions during this cell-cycle period, with the execution of cell-fate decisions.

Mechanistic insights into genetic susceptibility to prostate cancer.

Prostate cancer (PCa) is the second most common cancer in men and is a highly heritable disease that affects millions of individuals worldwide. Genome-wide association studies have to date discovered nearly 270 genetic loci harboring hundreds of single nucleotide polymorphisms (SNPs) that are associated with PCa susceptibility. In contrast, the functional characterization of the mechanisms underlying PCa risk association is still growing. Given that PCa risk-associated SNPs are highly enriched in noncoding cis-regulatory genomic regions, accumulating evidence suggests a widespread modulation of transcription factor chromatin binding and allelic enhancer activity by these noncoding SNPs, thereby dysregulating gene expression. Emerging studies have shown that a proportion of noncoding variants can modulate the formation of transcription factor complexes at enhancers and CTCF-mediated 3D genome architecture. Interestingly, DNA methylation-regulated CTCF binding could orchestrate a long-range chromatin interaction between PCa risk enhancer and causative genes. Additionally, one-causal-variant-two-risk genes or multiple-risk-variant-multiple-genes are prevalent in some PCa risk-associated loci. In this review, we will discuss the current understanding of the general principles of SNP-mediated gene regulation, experimental advances, and functional evidence supporting the mechanistic roles of several PCa genetic loci with potential clinical impact on disease prevention and treatment.

Construction of Condition-Specific Gene Regulatory Network Using Kernel Canonical Correlation Analysis.

Gene expression profile or transcriptome can represent cellular states, thus understanding gene regulation mechanisms can help understand how cells respond to external stress. Interaction between transcription factor (TF) and target gene (TG) is one of the representative regulatory mechanisms in cells. In this paper, we present a novel computational method to construct condition-specific transcriptional networks from transcriptome data. Regulatory interaction between TFs and TGs is very complex, specifically multiple-to-multiple relations. Experimental data from TF Chromatin Immunoprecipitation sequencing is useful but produces one-to-multiple relations between TF and TGs. On the other hand, co-expression networks of genes can be useful for constructing condition transcriptional networks, but there are many false positive relations in co-expression networks. In this paper, we propose a novel method to construct a condition-specific and combinatorial transcriptional network, applying kernel canonical correlation analysis (kernel CCA) to identify multiple-to-multiple TF–TG relations in certain biological condition. Kernel CCA is a well-established statistical method for computing the correlation of a group of features vs. another group of features. We, therefore, employed kernel CCA to embed TFs and TGs into a new space where the correlation of TFs and TGs are reflected. To demonstrate the usefulness of our network construction method, we used the blood transcriptome data for the investigation on the response to high fat diet in a human and an arabidopsis data set for the investigation on the response to cold/heat stress. Our method detected not only important regulatory interactions reported in previous studies but also novel TF–TG relations where a module of TF is regulating a module of TGs upon specific stress.

Transcription factor chromatin profiling genome-wide using uliCUT&RUN in single cells and individual blastocysts.

Determining chromatin-associated protein localization across the genome has provided insight into the functions of DNA-binding proteins and their connections to disease. However, established protocols requiring large quantities of cell or tissue samples currently limit applications for clinical and biomedical research in this field. Furthermore, most technologies have been optimized to assess abundant histone protein localization, prohibiting the investigation of nonhistone protein localization in low cell numbers. We recently described a protocol to profile chromatin-associated protein localization in as low as one cell: ultra-low-input cleavage under targets and release using nuclease (uliCUT&RUN). Optimized from chromatin immunocleavage and CUT&RUN, uliCUT&RUN is a tethered enzyme-based protocol that utilizes a combination of recombinant protein, antibody recognition and stringent purification to selectively target proteins of interest and isolate the associated DNA. Performed in native conditions, uliCUT&RUN profiles protein localization to chromatin with low input and high precision. Compared with other profiling technologies, uliCUT&RUN can determine nonhistone protein chromatin occupancies in low cell numbers, permitting the investigation into the molecular functions of a range of DNA-binding proteins within rare samples. From sample preparation to sequencing library submission, the uliCUT&RUN protocol takes <2 d to perform, with the accompanying data analysis timeline dependent on experience level.

Chromatin information content landscapes inform transcription factor and DNA interactions

Interactions between transcription factors and chromatin are fundamental to genome organization and regulation and, ultimately, cell state. Here, we use information theory to measure signatures of organized chromatin resulting from transcription factor-chromatin interactions encoded in the patterns of the accessible genome, which we term chromatin information enrichment (CIE). We calculate CIE for hundreds of transcription factor motifs across human samples and identify two classes: low and high CIE. The 10–20% of common and tissue-specific high CIE transcription factor motifs, associate with higher protein–DNA residence time, including different binding site subclasses of the same transcription factor, increased nucleosome phasing, specific protein domains, and the genetic control of both chromatin accessibility and gene expression. These results show that variations in the information encoded in chromatin architecture reflect functional biological variation, with implications for cell state dynamics and memory.

HAND1 and BARX1 Act as Transcriptional and Anatomic Determinants of Malignancy in Gastrointestinal Stromal Tumor

Gastrointestinal stromal tumor (GIST) arises from interstitial cells of Cajal (ICC) or their precursors, which are present throughout the gastrointestinal tract. Although gastric GIST is commonly indolent and small intestine GIST more aggressive, a molecular understanding of disease behavior would inform therapy decisions in GIST. Although a core transcription factor (TF) network is conserved across GIST, accessory TFs HAND1 and BARX1 are expressed in a disease state-specific pattern. Here, we characterize two divergent transcriptional programs maintained by HAND1 and BARX1, and evaluate their association with clinical outcomes. We evaluated RNA sequencing and TF chromatin immunoprecipitation with sequencing in GIST samples and cultured cells for transcriptional programs associated with HAND1 and BARX1. Multiplexed tissue-based cyclic immunofluorescence and IHC evaluated tissue- and cell-level expression of TFs and their association with clinical factors. We show that HAND1 is expressed in aggressive GIST, modulating KIT and core TF expression and supporting proliferative cellular programs. In contrast, BARX1 is expressed in indolent and micro-GISTs. HAND1 and BARX1 expression were superior predictors of relapse-free survival, as compared with standard risk stratification, and they predict progression-free survival on imatinib. Reflecting the developmental origins of accessory TF programs, HAND1 was expressed solely in small intestine ICCs, whereas BARX1 expression was restricted to gastric ICCs. Our results define anatomic and transcriptional determinants of GIST and molecular origins of clinical phenotypes. Assessment of HAND1 and BARX1 expression in GIST may provide prognostic information and improve clinical decisions on the administration of adjuvant therapy.

One-Day TALEN Assembly Protocol and a Dual-Tagging System for Genome Editing.

This study developed a new rapid transcription activator-like effector nuclease (TALEN) preparation protocol by thoroughly redesigning the widely used Golden Gate TALEN and TAL Effector Kit 2.0. The new protocol can be used to prepare any custom 18-bp binding TALENs in just one day (about 12 h), more rapidly than CRISPR. This protocol used a set of linear monomers, a final TALE-FokI backbone plasmid, and a pipeline to assemble the ready-to-use TALEN expression plasmid, which were all newly developed for this study. The set of linear monomers can be easily produced and reproduced by high-fidelity polymerase chain reaction (PCR) amplification in a 96-well plate using a pair of universal primers. Most important of all, our rapid TALEN construction pipeline can easily obtain many positive colonies with high efficiency (over 80%). By preparing five pairs of TALENs targeting five NF-κB genes (RELA, RELB, CREL,NFKB1, and NFKB2) and editing these genes in different cell lines (293T, HepG2, and PANC1), this study demonstrated that the new protocol has high efficiency, reproducibility, reliability, and applicability. Moreover, this study showed that the fabricated TALEN has much higher editing efficiency than CRISPR. Finally, this study developed a dual-tagging system for simultaneously tagging target proteins and successfully edited cells with a streptavidin-binding peptide (SBP) or AviTag via homology-directed repair, which could have wide applications in protein (antigen) preparation, immunoprecipitation, and a transcription factor chromatin immunoprecipitation assay.

Determinants of transcription factor regulatory range

Characterization of the genomic distances over which transcription factor (TF) binding influences gene expression is important for inferring target genes from TF chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Here we systematically examine the relationship between thousands of TF and histone modification ChIP-seq data sets with thousands of gene expression profiles. We develop a model for integrating these data, which reveals two classes of TFs with distinct ranges of regulatory influence, chromatin-binding preferences, and auto-regulatory properties. We find that the regulatory range of the same TF bound within different topologically associating domains (TADs) depend on intrinsic TAD properties such as local gene density and G/C content, but also on the TAD chromatin states. Our results suggest that considering TF type, binding distance to gene locus, as well as chromatin context is important in identifying implicated TFs from GWAS SNPs.

C/EBPα and GATA-2 Mutations Induce Bilineage Acute Erythroid Leukemia through Transformation of a Neomorphic Neutrophil-Erythroid Progenitor

SummaryAcute erythroid leukemia (AEL) commonly involves both myeloid and erythroid lineage transformation. However, the mutations that cause AEL and the cell(s) that sustain the bilineage leukemia phenotype remain unknown. We here show that combined biallelic Cebpa and Gata2 zinc finger-1 (ZnF1) mutations cooperatively induce bilineage AEL, and that the major leukemia-initiating cell (LIC) population has a neutrophil-monocyte progenitor (NMP) phenotype. In pre-leukemic NMPs Cebpa and Gata2 mutations synergize by increasing erythroid transcription factor (TF) expression and erythroid TF chromatin access, respectively, thereby installing ectopic erythroid potential. This erythroid-permissive chromatin conformation is retained in bilineage LICs. These results demonstrate that synergistic transcriptional and epigenetic reprogramming by leukemia-initiating mutations can generate neomorphic pre-leukemic progenitors, defining the lineage identity of the resulting leukemia.

Genetic and Epigenetic Control of Heart Development.

A transcriptional program implemented by transcription factors and epigenetic regulators governs cardiac development and disease. Mutations in these factors are important causes of congenital heart disease. Here, we review selected recent advances in our understanding of the transcriptional and epigenetic control of heart development, including determinants of cardiac transcription factor chromatin occupancy, the gene regulatory network that regulates atrial septation, the chromatin landscape and cardiac gene regulation, and the role of Brg/Brahma-associated factor (BAF), nucleosome remodeling and histone deacetylation (NuRD), and Polycomb epigenetic regulatory complexes in heart development.

A reference map of murine cardiac transcription factor chromatin occupancy identifies dynamic and conserved enhancers

Mapping the chromatin occupancy of transcription factors (TFs) is a key step in deciphering developmental transcriptional programs. Here we use biotinylated knockin alleles of seven key cardiac TFs (GATA4, NKX2-5, MEF2A, MEF2C, SRF, TBX5, TEAD1) to sensitively and reproducibly map their genome-wide occupancy in the fetal and adult mouse heart. These maps show that TF occupancy is dynamic between developmental stages and that multiple TFs often collaboratively occupy the same chromatin region through indirect cooperativity. Multi-TF regions exhibit features of functional regulatory elements, including evolutionary conservation, chromatin accessibility, and activity in transcriptional enhancer assays. H3K27ac, a feature of many enhancers, incompletely overlaps multi-TF regions, and multi-TF regions lacking H3K27ac retain conservation and enhancer activity. TEAD1 is a core component of the cardiac transcriptional network, co-occupying cardiac regulatory regions and controlling cardiomyocyte-specific gene functions. Our study provides a resource for deciphering the cardiac transcriptional regulatory network and gaining insights into the molecular mechanisms governing heart development.

A reference map of murine cardiac transcription factor chromatin occupancy identifies dynamic and conserved enhancers

A reference map of murine cardiac transcription factor chromatin occupancy identifies dynamic and conserved enhancers

Inherent DNA‐binding specificities of the HIF‐1α and HIF‐2α transcription factors in chromatin

Hypoxia‐inducible factor (HIF) is the major transcriptional regulator of cellular responses to hypoxia. The two principal HIF‐α isoforms, HIF‐1α and HIF‐2α, are progressively stabilized in response to hypoxia and form heterodimers with HIF‐1β to activate a broad range of transcriptional responses. Here, we report on the pan‐genomic distribution of isoform‐specific HIF binding in response to hypoxia of varying severity and duration, and in response to genetic ablation of each HIF‐α isoform. Our findings reveal that, despite an identical consensus recognition sequence in DNA, each HIF heterodimer loads progressively at a distinct repertoire of cell‐type‐specific sites across the genome, with little evidence of redistribution under any of the conditions examined. Marked biases towards promoter‐proximal binding of HIF‐1 and promoter‐distant binding of HIF‐2 were observed under all conditions and were consistent in multiple cell type. The findings imply that each HIF isoform has an inherent property that determines its binding distribution across the genome, which might be exploited to therapeutically target the specific transcriptional output of each isoform independently.

Targeted DamID reveals differential binding of mammalian pluripotency factors.

The precise control of gene expression by transcription factor networks is crucial to organismal development. The predominant approach for mapping transcription factor-chromatin interactions has been chromatin immunoprecipitation (ChIP). However, ChIP requires a large number of homogeneous cells and antisera with high specificity. A second approach, DamID, has the drawback that high levels of Dam methylase are toxic. Here, we modify our targeted DamID approach (TaDa) to enable cell type-specific expression in mammalian systems, generating an inducible system (mammalian TaDa or MaTaDa) to identify genome-wide protein/DNA interactions in 100 to 1000 times fewer cells than ChIP-based approaches. We mapped the binding sites of two key pluripotency factors, OCT4 and PRDM14, in mouse embryonic stem cells, epiblast-like cells and primordial germ cell-like cells (PGCLCs). PGCLCs are an important system for elucidating primordial germ cell development in mice. We monitored PRDM14 binding during the specification of PGCLCs, identifying direct targets of PRDM14 that are key to understanding its crucial role in PGCLC development. We show that MaTaDa is a sensitive and accurate method for assessing cell type-specific transcription factor binding in limited numbers of cells.

DNA methylation at enhancers identifies distinct breast cancer lineages

Breast cancers exhibit genome-wide aberrant DNA methylation patterns. To investigate how these affect the transcriptome and which changes are linked to transformation or progression, we apply genome-wide expression–methylation quantitative trait loci (emQTL) analysis between DNA methylation and gene expression. On a whole genome scale, in cis and in trans, DNA methylation and gene expression have remarkably and reproducibly conserved patterns of association in three breast cancer cohorts (n = 104, n = 253 and n = 277). The expression–methylation quantitative trait loci associations form two main clusters; one relates to tumor infiltrating immune cell signatures and the other to estrogen receptor signaling. In the estrogen related cluster, using ChromHMM segmentation and transcription factor chromatin immunoprecipitation sequencing data, we identify transcriptional networks regulated in a cell lineage-specific manner by DNA methylation at enhancers. These networks are strongly dominated by ERα, FOXA1 or GATA3 and their targets were functionally validated using knockdown by small interfering RNA or GRO-seq analysis after transcriptional stimulation with estrogen.