Highly expressed loci are vulnerable to misleading ChIP localization of multiple unrelated proteins

Benchmark of chromatin-protein interaction methods in haploid round spermatids

Introduction

Chromatin-protein interactions are fundamental for regulation of gene transcription. While chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) has long been the gold standard for mapping these interactions, emerging techniques such as CUT&RUN and CUT&Tag, which offer advantages such as low-input requirements and high signal-to-noise ratios, have aroused great attention. However, research addressing the potential biases introduced by enzyme-based tagmentation approaches and comparative assessment with ChIP-seq remain absent.

Methods

This study aims to systematically evaluate and compare the performance of ChIP-seq, CUT&Tag, and CUT&RUN for profiling genome-wide transcription factors and histone modification binding.

Results

Our analysis revealed that all three methods reliably detect histone modifications and transcription factor enrichment, with CUT&Tag standing out for its comparatively higher signal-to-noise ratio. Detailed peak comparison revealed unique and overlapping enrichment among the three techniques. Additionally, CUT&Tag can identify novel CTCF peaks compared with the other two methods. A strong correlation was observed between CUT&Tag signal intensity and chromatin accessibility, highlighting its ability to generate high-resolution signals in accessible regions.

Discussion

The systematic comparison summarizes the differences between CUT&Tag and CUT&RUN in terms of the signal-to-noise ratio and bias toward accessible chromatin. Considering the experimental procedures, signal specificity, and inherent biases, we recommend tailoring the choice of method to the type of chromatin-protein interaction under study. CUT&Tag offers a promising alternative for applications requiring high sensitivity and reduced background noise.

Optimized CUT&RUN protocol for activated primary mouse B cells

ChIP-seq has long been the standard for study of chromatin-protein interactions. However, development of a new technique, CUT&RUN, showed substantial advantages compared to ChIP-seq including higher quality signal while using substantially less sample. While a powerful technique, the original protocol was designed using cell lines and histones as targets. Due to their fragility, this was unsuitable for obtaining high-quality data from activated primary B lymphocytes. To adapt this protocol for B cells, cells were fixed prior to nuclear isolation, and several critical adjustments were introduced to the procedure and reagents. We measured binding of H3K4me3 histone and RNA Polymerase II, detecting robust peaks with as little as 100k nuclei. Additionally, freeze-thaw of B cells prior to processing did not affect results, emphasizing the flexibility of this modified technique. Using the protocol described here will allow one to quantify non-histone proteins bound to DNA from limited numbers of B cells with more efficiency than can be achieved from the current standard, ChIP-seq.

CUT&Tag recovers up to half of ENCODE ChIP-seq histone acetylation peaks

DNA-protein interactions have traditionally been profiled via chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq). Cleavage Under Targets & Tagmentation (CUT&Tag) is a rapidly expanding technique that enables the profiling of such interactions in situ at high sensitivity. However, thorough evaluation and benchmarking against established ChIP-seq datasets are lacking. Here, we comprehensively benchmarked CUT&Tag for H3K27ac and H3K27me3 against published ChIP-seq profiles from ENCODE in K562 cells. Combining multiple new and published CUT&Tag datasets, there was an average recall of 54% known ENCODE peaks for both histone modifications. We tested peak callers MACS2 and SEACR and identified optimal peak calling parameters. Overall, peaks identified by CUT&Tag represent the strongest ENCODE peaks and show the same functional and biological enrichments as ChIP-seq peaks identified by ENCODE. Our workflow systematically evaluates the merits of methodological adjustments, providing a benchmarking framework for the experimental design and analysis of CUT&Tag studies.

SMC translocation is unaffected by an excess of nucleoid associated proteins in vivo

Genome organization is important for DNA replication, gene expression, and chromosome segregation. In bacteria, two large families of proteins, nucleoid-associated proteins (NAPs) and SMC complexes, play important roles in organizing the genome. NAPs are highly abundant DNA-binding proteins that can bend, wrap, bridge, and compact DNA, while SMC complexes load onto the chromosome, translocate on the DNA, and extrude DNA loops. Although SMC complexes are capable of traversing the entire chromosome bound by various NAPs in vivo, it is unclear whether SMC translocation is influenced by NAPs. In this study, using Bacillus subtilis as a model system, we expressed a collection of representative bacterial and archaeal DNA-binding proteins that introduce distinct DNA structures and potentially pose different challenges for SMC movement. By fluorescence microscopy and chromatin immunoprecipitation, we observed that these proteins bound to the genome in characteristic manners. Using genome-wide chromosome conformation capture (Hi-C) assays, we found that the SMC complex traversed these DNA-binding proteins without slowing down. Our findings revealed that the DNA-loop-extruding activity of the SMC complex is unaffected by exogenously expressed DNA-binding proteins, which highlights the robustness of SMC motors on the busy chromatin.

Predicting genome-wide tissue-specific enhancers via combinatorial transcription factor genomic occupancy analysis

Enhancers are non-coding cis-regulatory elements crucial for transcriptional regulation. Mutations in enhancers can disrupt gene regulation, leading to disease phenotypes. Identifying enhancers and their tissue-specific activity is challenging due to their lack of stereotyped sequences. This study presents a sequence-based computational model that uses combinatorial transcription factor (TF) genomic occupancy to predict tissue-specific enhancers. Trained on diverse datasets, including ENCODE and Vista enhancer browser data, the model predicted 25 000 forebrain-specific cis-regulatory modules (CRMs) in the human genome. Validation using biochemical features, disease-associated SNPs, and in vivo zebrafish analysis confirmed its effectiveness. This model aids in predicting enhancers lacking well-characterized chromatin features, complementing experimental approaches in tissue-specific enhancer discovery.

Chromatin endogenous cleavage provides a global view of yeast RNA polymerase II transcription kinetics

Chromatin immunoprecipitation (ChIP-seq) is the most common approach to observe global binding of proteins to DNA in vivo. The occupancy of transcription factors (TFs) from ChIP-seq agrees well with an alternative method, chromatin endogenous cleavage (ChEC-seq2). However, ChIP-seq and ChEC-seq2 reveal strikingly different patterns of enrichment of yeast RNA polymerase II (RNAPII). We hypothesized that this reflects distinct populations of RNAPII, some of which are captured by ChIP-seq and some of which are captured by ChEC-seq2. RNAPII association with enhancers and promoters - predicted from biochemical studies - is detected well by ChEC-seq2 but not by ChIP-seq. Enhancer/promoter-bound RNAPII correlates with transcription levels and matches predicted occupancy based on published rates of enhancer recruitment, preinitiation assembly, initiation, elongation, and termination. The occupancy from ChEC-seq2 allowed us to develop a stochastic model for global kinetics of RNAPII transcription which captured both the ChEC-seq2 data and changes upon chemical-genetic perturbations to transcription. Finally, RNAPII ChEC-seq2 and kinetic modeling suggests that a mutation in the Gcn4 transcription factor that blocks interaction with the NPC destabilizes promoter-associated RNAPII without altering its recruitment to the enhancer.

Functional characteristics and computational model of abundant hyperactive loci in the human genome

Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

LAP2alpha facilitates myogenic gene expression by preventing nucleoplasmic lamin A/C from spreading to active chromatin regions

A-type lamins form a filamentous meshwork beneath the nuclear membrane that anchors large heterochromatic genomic regions at the nuclear periphery. A-type lamins also exist as a dynamic, non-filamentous pool in the nuclear interior, where they interact with lamin-associated polypeptide 2 alpha (LAP2α). Both proteins associate with largely overlapping euchromatic genomic regions in the nucleoplasm, but the functional significance of this interaction is poorly understood. Here, we report that LAP2α relocates towards regions containing myogenic genes in the early stages of muscle differentiation, possibly facilitating efficient gene regulation, while lamins A and C mostly associate with genomic regions away from these genes. Strikingly, upon depletion of LAP2α, A-type lamins spread across active chromatin and accumulate at regions of active H3K27ac and H3K4me3 histone marks in the vicinity of myogenic genes whose expression is impaired in the absence of LAP2α. Reorganization of A-type lamins on chromatin is accompanied by depletion of the active chromatin mark H3K27ac and a significantly impaired myogenic differentiation. Thus, the interplay of LAP2α and A-type lamins is crucial for proper positioning of intranuclear lamin A/C on chromatin to allow efficient myogenic differentiation.

Chromatin endogenous cleavage provides a global view of yeast RNA polymerase II transcription kinetics

Chromatin immunoprecipitation (ChIP-seq) is the most common approach to observe global binding of proteins to DNA in vivo. The occupancy of transcription factors (TFs) from ChIP-seq agrees well with an alternative method, chromatin endogenous cleavage (ChEC-seq2). However, ChIP-seq and ChEC-seq2 reveal strikingly different patterns of enrichment of yeast RNA polymerase II. We hypothesized that this reflects distinct populations of RNAPII, some of which are captured by ChIP-seq and some of which are captured by ChEC-seq2. RNAPII association with enhancers and promoters - predicted from biochemical studies - is detected well by ChEC-seq2 but not by ChIP-seq. Enhancer/promoter bound RNAPII correlates with transcription levels and matches predicted occupancy based on published rates of enhancer recruitment, preinitiation assembly, initiation, elongation and termination. The occupancy from ChEC-seq2 allowed us to develop a stochastic model for global kinetics of RNAPII transcription which captured both the ChEC-seq2 data and changes upon chemical-genetic perturbations to transcription. Finally, RNAPII ChEC-seq2 and kinetic modeling suggests that a mutation in the Gcn4 transcription factor that blocks interaction with the NPC destabilizes promoter-associated RNAPII without altering its recruitment to the enhancer.

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Functional characteristics and computational model of abundant hyperactive loci in the human genome

Enhancers and promoters are classically considered to be bound by a small set of TFs in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1,003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected 5 distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

Emerging toolkits for decoding the co-occurrence of modified histones and chromatin proteins

In eukaryotes, DNA is packaged into chromatin with the help of highly conserved histone proteins. Together with DNA-binding proteins, posttranslational modifications (PTMs) on these histones play crucial roles in regulating genome function, cell fate determination, inheritance of acquired traits, cellular states, and diseases. While most studies have focused on individual DNA-binding proteins, chromatin proteins, or histone PTMs in bulk cell populations, such chromatin features co-occur and potentially act cooperatively to accomplish specific functions in a given cell. This review discusses state-of-the-art techniques for the simultaneous profiling of multiple chromatin features in low-input samples and single cells, focusing on histone PTMs, DNA-binding, and chromatin proteins. We cover the origins of the currently available toolkits, compare and contrast their characteristic features, and discuss challenges and perspectives for future applications. Studying the co-occurrence of histone PTMs, DNA-binding proteins, and chromatin proteins in single cells will be central for a better understanding of the biological relevance of combinatorial chromatin features, their impact on genomic output, and cellular heterogeneity.

White-opaque switching in Candida albicans: cell biology, regulation, and function

SUMMARYCandida albicans remains a major fungal pathogen colonizing humans and opportunistically invading tissue when conditions are predisposing. Part of the success of C. albicans was attributed to its capacity to form hyphae that facilitate tissue invasion. However, in 1987, a second developmental program was discovered, the "white-opaque transition," a high-frequency reversible switching system that impacted most aspects of the physiology, cell architecture, virulence, and gene expression of C. albicans. For the 15 years following the discovery of white-opaque switching, its role in the biology of C. albicans remained elusive. Then in 2002, it was discovered that in order to mate, C. albicans had to switch from white to opaque, a unique step in a yeast mating program. In 2006, three laboratories simultaneously identified a putative master switch gene, which led to a major quest to elucidate the underlying mechanisms that regulate white-opaque switching. Here, the evolving discoveries related to this complicated phenotypic transition are reviewed in a quasi-chronological order not only to provide a historical perspective but also to highlight several unique characteristics of white-opaque switching, which are fascinating and may be important to the life history and virulence of this persistent pathogen. Many of these characteristics have not been fully investigated, in many cases, leaving intriguing questions unresolved. Some of these include the function of unique channeled pimples on the opaque cell wall, the capacity to form opaque cells in the absence of the master switch gene WOR1, the formation of separate "pathogenic" and "sexual" biofilms, and the possibility that a significant portion of natural strains colonizing the lower gastrointestinal tract may be in the opaque phase. This review addresses many of these characteristics with the intent of engendering interest in resolving questions that remain unanswered.

Growth-regulated co-occupancy of Mediator and Lsm3 at intronic ribosomal protein genes

Mediator is a well-known transcriptional co-regulator and serves as an adaptor between gene-specific regulatory proteins and RNA polymerase II. Studies on the chromatin-bound form of Mediator revealed interactions with additional protein complexes involved in various transcription-related processes, such as the Lsm2-8 complex that is part of the spliceosomal U6 small nuclear ribonucleoprotein complex. Here, we employ Chromatin Immunoprecipitation sequencing (ChIP-seq) of chromatin associated with the Lsm3 protein and the Med1 or Med15 Mediator subunits. We identify 86 genes co-occupied by both Lsm3 and Mediator, of which 73 were intron-containing ribosomal protein genes. In logarithmically growing cells, Mediator primarily binds to their promoter regions but also shows a second, less pronounced occupancy at their 3'-exons. During the late exponential phase, we observe a near-complete transition of Mediator from these promoters to a position in their 3'-ends, overlapping the Lsm3 binding sites ∼250 bp downstream of their last intron-exon boundaries. Using an unbiased RNA sequencing approach, we show that transition of Mediator from promoters to the last exon of these genes correlates to reduction of both their messenger RNA levels and splicing ratios, indicating that the Mediator and Lsm complexes cooperate to control growth-regulated expression of intron-containing ribosomal protein genes at the levels of transcription and splicing.

Regional centromere configuration in the fungal pathogens of the Pneumocystis genus

Centromeres are constricted chromosomal regions that are essential for cell division. In eukaryotes, centromeres display a remarkable architectural and genetic diversity. The basis of centromere-accelerated evolution remains elusive. Here, we focused on Pneumocystis species, a group of mammalian-specific fungal pathogens that form a sister taxon with that of the Schizosaccharomyces pombe, an important genetic model for centromere biology research. Methods allowing reliable continuous culture of Pneumocystis species do not currently exist, precluding genetic manipulation. CENP-A, a variant of histone H3, is the epigenetic marker that defines centromeres in most eukaryotes. Using heterologous complementation, we show that the Pneumocystis CENP-A ortholog is functionally equivalent to CENP-ACnp1 of S. pombe. Using organisms from a short-term in vitro culture or infected animal models and chromatin immunoprecipitation (ChIP)-Seq, we identified CENP-A bound regions in two Pneumocystis species that diverged ~35 million years ago. Each species has a unique short regional centromere (<10 kb) flanked by heterochromatin in 16-17 monocentric chromosomes. They span active genes and lack conserved DNA sequence motifs and repeats. These features suggest an epigenetic specification of centromere function. Analysis of centromeric DNA across multiple Pneumocystis species suggests a vertical transmission at least 100 million years ago. The common ancestry of Pneumocystis and S. pombe centromeres is untraceable at the DNA level, but the overall architectural similarity could be the result of functional constraint for successful chromosomal segregation.IMPORTANCEPneumocystis species offer a suitable genetic system to study centromere evolution in pathogens because of their phylogenetic proximity with the non-pathogenic yeast S. pombe, a popular model for cell biology. We used this system to explore how centromeres have evolved after the divergence of the two clades ~ 460 million years ago. To address this question, we established a protocol combining short-term culture and ChIP-Seq to characterize centromeres in multiple Pneumocystis species. We show that Pneumocystis have short epigenetic centromeres that function differently from those in S. pombe.

Less-is-more: selecting transcription factor binding regions informative for motif inference

Numerous statistical methods have emerged for inferring DNA motifs for transcription factors (TFs) from genomic regions. However, the process of selecting informative regions for motif inference remains understudied. Current approaches select regions with strong ChIP-seq signal for a given TF, assuming that such strong signal primarily results from specific interactions between the TF and its motif. Additionally, these selection approaches do not account for non-target motifs, i.e. motifs of other TFs; they presume the occurrence of these non-target motifs infrequent compared to that of the target motif, and thus assume these have minimal interference with the identification of the target. Leveraging extensive ChIP-seq datasets, we introduced the concept of TF signal 'crowdedness', referred to as C-score, for each genomic region. The C-score helps in highlighting TF signals arising from non-specific interactions. Moreover, by considering the C-score (and adjusting for the length of genomic regions), we can effectively mitigate interference of non-target motifs. Using these tools, we find that in many instances, strong ChIP-seq signal stems mainly from non-specific interactions, and the occurrence of non-target motifs significantly impacts the accurate inference of the target motif. Prioritizing genomic regions with reduced crowdedness and short length markedly improves motif inference. This 'less-is-more' effect suggests that ChIP-seq region selection warrants more attention.

Epigenomic Profiling of B Cell Subsets by CUT&Tag

Epigenetic programs play a key role in regulating the development and function of immune cells. However, conventional methods for profiling epigenetic mechanisms, such as the post-translational modifications to histones, present several technical challenges that prevent a complete understanding of gene regulation. Here, we provide a detailed protocol of the Cleavage Under Targets and Tagmentation (CUT&Tag) chromatin profiling technique for identifying histone modifications in human and mouse lymphocytes.

Integrative Analysis of CUT&Tag and RNA-Seq Data Through Bioinformatics: A Unified Workflow for Enhanced Insights

Cleavage Under Targets and Tagmentation (CUT&Tag) is a recent methodology used for robust epigenomic profiling that, unlike conventional chromatin immunoprecipitation (ChIP-Seq), requires only a limited amount of cells as starting material. RNA sequencing (RNA-Seq) reveals the presence and quantity of RNA in a biological sample, describing the continuously changing cellular transcriptome. The integrated analysis of transcriptional activity, histone modifications, and chromatin accessibility via CUT&Tag is still in its infancy compared to the well-established ChIP-Seq. This chapter describes a robust bioinformatics methodology and workflow to perform an integrative CUT&Tag/RNA-Seq analysis.

Vitamin D status associates with skeletal muscle loss after anterior cruciate ligament reconstruction

BACKGROUNDAlthough 25-hydroxyvitamin D [25(OH)D] concentrations of 30 ng/mL or higher are known to reduce injury risk and boost strength, the influence on anterior cruciate ligament reconstruction (ACLR) outcomes remains unexamined. This study aimed to define the vitamin D signaling response to ACLR, assess the relationship between vitamin D status and muscle fiber cross-sectional area (CSA) and bone density outcomes, and discover vitamin D receptor (VDR) targets after ACLR.METHODSTwenty-one young, healthy, physically active participants with recent ACL tears were enrolled (17.8 ± 3.2 years, BMI 26.0 ± 3.5 kg/m2). Data were collected through blood samples, vastus lateralis biopsies, dual energy x-ray bone density measurements, and isokinetic dynamometer measures at baseline, 1 week, 4 months, and 6 months after ACLR. The biopsies facilitated CSA, Western blotting, RNA-seq, and VDR ChIP-seq analyses.RESULTSACLR surgery led to decreased circulating bioactive vitamin D and increased VDR and activating enzyme expression in skeletal muscle 1 week after ACLR. Participants with less than 30 ng/mL 25(OH)D levels (n = 13) displayed more significant quadriceps fiber CSA loss 1 week and 4 months after ACLR than those with 30 ng/mL or higher (n = 8; P < 0.01 for post hoc comparisons; P = 0.041 for time × vitamin D status interaction). RNA-seq and ChIP-seq data integration revealed genes associated with energy metabolism and skeletal muscle recovery, potentially mediating the impact of vitamin D status on ACLR recovery. No difference in bone mineral density losses between groups was observed.CONCLUSIONCorrecting vitamin D status prior to ACLR may aid in preserving skeletal muscle during recovery.FUNDINGNIH grants R01AR072061, R01AR071398-04S1, and K99AR081367.

Locus-specific proteome decoding reveals Fpt1 as a chromatin-associated negative regulator of RNA polymerase III assembly

Transcription of tRNA genes by RNA polymerase III (RNAPIII) is tuned by signaling cascades. The emerging notion of differential tRNA gene regulation implies the existence of additional regulatory mechanisms. However, tRNA gene-specific regulators have not been described. Decoding the local chromatin proteome of a native tRNA gene in yeast revealed reprogramming of the RNAPIII transcription machinery upon nutrient perturbation. Among the dynamic proteins, we identified Fpt1, a protein of unknown function that uniquely occupied RNAPIII-regulated genes. Fpt1 binding at tRNA genes correlated with the efficiency of RNAPIII eviction upon nutrient perturbation and required the transcription factors TFIIIB and TFIIIC but not RNAPIII. In the absence of Fpt1, eviction of RNAPIII was reduced, and the shutdown of ribosome biogenesis genes was impaired upon nutrient perturbation. Our findings provide support for a chromatin-associated mechanism required for RNAPIII eviction from tRNA genes and tuning the physiological response to changing metabolic demands.

Adaptation of CUT&RUN for use in African trypanosomes

This Cleavage Under Targets and Release Using Nuclease (CUT&RUN) protocol produces genomic occupancy data for a protein of interest in the protozoan parasite Trypanosoma brucei. The data produced is analyzed in a similar way as that produced by ChIP-seq. While we describe the protocol for parasites carrying an epitope tag for the protein of interest, antibodies against the native protein could be used for the same purpose.

Systematic Dissection of Sequence Features Affecting the Binding Specificity of a Pioneer Factor Reveals Binding Synergy Between FOXA1 and AP-1

Despite the unique ability of pioneer transcription factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called ChIP-ISO to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1. Combining ChIP-ISO with in vitro and neural network analyses, we find that 1) FOXA1 binding is strongly affected by co-binding TFs AP-1 and CEBPB, 2) FOXA1 and AP-1 show binding cooperativity in vitro, 3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin, and 4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

The selfish yeast plasmid exploits a SWI/SNF-type chromatin remodeling complex for hitchhiking on chromosomes and ensuring high-fidelity propagation

Extra-chromosomal selfish DNA elements can evade the risk of being lost at every generation by behaving as chromosome appendages, thereby ensuring high fidelity segregation and stable persistence in host cell populations. The yeast 2-micron plasmid and episomes of the mammalian gammaherpes and papilloma viruses that tether to chromosomes and segregate by hitchhiking on them exemplify this strategy. We document for the first time the utilization of a SWI/SNF-type chromatin remodeling complex as a conduit for chromosome association by a selfish element. One principal mechanism for chromosome tethering by the 2-micron plasmid is the bridging interaction of the plasmid partitioning proteins (Rep1 and Rep2) with the yeast RSC2 complex and the plasmid partitioning locus STB. We substantiate this model by multiple lines of evidence derived from genomics, cell biology and interaction analyses. We describe a Rep-STB bypass system in which a plasmid engineered to non-covalently associate with the RSC complex mimics segregation by chromosome hitchhiking. Given the ubiquitous prevalence of SWI/SNF family chromatin remodeling complexes among eukaryotes, it is likely that the 2-micron plasmid paradigm or analogous ones will be encountered among other eukaryotic selfish elements.

A R-loop sensing pathway mediates the relocation of transcribed genes to nuclear pore complexes

Nuclear pore complexes (NPCs) have increasingly recognized interactions with the genome, as exemplified in yeast, where they bind transcribed or damaged chromatin. By combining genome-wide approaches with live imaging of model loci, we uncover a correlation between NPC association and the accumulation of R-loops, which are genotoxic structures formed through hybridization of nascent RNAs with their DNA templates. Manipulating hybrid formation demonstrates that R-loop accumulation per se, rather than transcription or R-loop-dependent damages, is the primary trigger for relocation to NPCs. Mechanistically, R-loop-dependent repositioning involves their recognition by the ssDNA-binding protein RPA, and SUMO-dependent interactions with NPC-associated factors. Preventing R-loop-dependent relocation leads to lethality in hybrid-accumulating conditions, while NPC tethering of a model hybrid-prone locus attenuates R-loop-dependent genetic instability. Remarkably, this relocation pathway involves molecular factors similar to those required for the association of stalled replication forks with NPCs, supporting the existence of convergent mechanisms for sensing transcriptional and genotoxic stresses.

The impact of the chromatin binding DEK protein in hematopoiesis and acute myeloid leukemia

Hematopoiesis is an exquisitely regulated process of cellular differentiation to create diverse cell types of the blood. Genetic mutations, or aberrant regulation of gene transcription, can interrupt normal hematopoiesis. This can have dire pathological consequences, including acute myeloid leukemia (AML), in which generation of the myeloid lineage of differentiated cells is interrupted. In this literature review, we discuss how the chromatin remodeling DEK protein can control hematopoietic stem cell quiescence, hematopoietic progenitor cell proliferation, and myelopoiesis. We further discuss the oncogenic consequences of the t(6;9) chromosomal translocation, which creates the DEK-NUP214 (aka: DEK-CAN) fusion gene, during the pathogenesis of AML. Combined, the literature indicates that DEK is crucial for maintaining homeostasis of hematopoietic stem and progenitor cells, including myeloid progenitors.

A variegated model of transcription factor function in the immune system

Specific combinations of transcription factors (TFs) control the gene expression programs that underlie specialized immune responses. Previous models of TF function in immunocytes had restricted each TF to a single functional categorization [e.g., lineage-defining (LDTFs) vs. signal-dependent TFs (SDTFs)] within one cell type. Synthesizing recent results, we instead propose a variegated model of immunological TF function, whereby many TFs have flexible and different roles across distinct cell states, contributing to cell phenotypic diversity. We discuss evidence in support of this variegated model, describe contextual inputs that enable TF diversification, and look to the future to imagine warranted experimental and computational tools to build quantitative and predictive models of immunocyte gene regulatory networks.

Broadly Applicable Control Approaches Improve Accuracy of ChIP-Seq Data

Chromatin ImmunoPrecipitation (ChIP) is a widely used method for the analysis of protein-DNA interactions in vivo; however, ChIP has pitfalls, particularly false-positive signal enrichment that permeates the data. We have developed a new approach to control for non-specific enrichment in ChIP that involves the expression of a non-genome-binding protein targeted in the IP alongside the experimental target protein due to the sharing of epitope tags. ChIP of the protein provides a "sensor" for non-specific enrichment that can be used for the normalization of the experimental data, thereby correcting for non-specific signals and improving data quality as validated against known binding sites for several proteins that we tested, including Fkh1, Orc1, Mcm4, and Sir2. We also tested a DNA-binding mutant approach and showed that, when feasible, ChIP of a site-specific DNA-binding mutant of the target protein is likely an ideal control. These methods vastly improve our ChIP-seq results in S. cerevisiae and should be applicable in other systems.

Context dependent function of the transcriptional regulator Rap1 in gene silencing and activation in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, heterochromatin is formed through interactions between site-specific DNA-binding factors, including the transcriptional activator Rap1, and Sir proteins. Despite a vast understanding of the establishment and maintenance of Sir-silenced chromatin, the mechanism of gene silencing by Sir proteins has remained a mystery. Utilizing high resolution chromatin immunoprecipitation, we found that Rap1, the native activator of the bi-directional HML α promoter, bound its recognition sequence in silenced chromatin and its binding was enhanced by the presence of Sir proteins. In contrast to prior results, various components of transcription machinery were not able to access HML α in the silenced state. These findings disproved the long-standing model of indiscriminate steric occlusion by Sir proteins and led to investigation of the transcriptional activator Rap1 in Sir-silenced chromatin. Using a highly sensitive assay that monitors loss-of-silencing events, we identified a novel role for promoter-bound Rap1 in the maintenance of silent chromatin through interactions with the Sir complex. We also found that promoter-bound Rap1 activated HML α when in an expressed state, and aided in the transition from transcription initiation to elongation. Highlighting the importance of epigenetic context in transcription factor function, these results point toward a model in which the duality of Rap1 function was mediated by local chromatin environment rather than binding-site availability.

Significance Statement

The coarse partitioning of the genome into regions of active euchromatin and repressed heterochromatin is an important, and conserved, level gene expression regulation in eukaryotes. Repressor Activator Protein (Rap1) is a transcription factor that promotes the activation of genes when recruited to promoters, and aids in the establishment of heterochromatin through interactions with silencer elements. Here, we investigate the role of Rap1 when bound to a promoter in silent chromatin and dissect the context-specific epigenetic cues that regulate the dual properties of this transcription factor. Together, our data highlight the importance of protein-protein interactions and local chromatin state on transcription factor function.

Characterizing the targets of transcription regulators by aggregating ChIP-seq and perturbation expression data sets

Mapping the gene targets of chromatin-associated transcription regulators (TRs) is a major goal of genomics research. ChIP-seq of TRs and experiments that perturb a TR and measure the differential abundance of gene transcripts are a primary means by which direct relationships are tested on a genomic scale. It has been reported that there is a poor overlap in the evidence across gene regulation strategies, emphasizing the need for integrating results from multiple experiments. Although research consortia interested in gene regulation have produced a valuable trove of high-quality data, there is an even greater volume of TR-specific data throughout the literature. In this study, we show a workflow for the identification, uniform processing, and aggregation of ChIP-seq and TR perturbation experiments for the ultimate purpose of ranking human and mouse TR-target interactions. Focusing on an initial set of eight regulators (ASCL1, HES1, MECP2, MEF2C, NEUROD1, PAX6, RUNX1, and TCF4), we identified 497 experiments suitable for analysis. We used this corpus to examine data concordance, to identify systematic patterns of the two data types, and to identify putative orthologous interactions between human and mouse. We build upon commonly used strategies to forward a procedure for aggregating and combining these two genomic methodologies, assessing these rankings against independent literature-curated evidence. Beyond a framework extensible to other TRs, our work also provides empirically ranked TR-target listings, as well as transparent experiment-level gene summaries for community use.

A dual-activity topoisomerase complex promotes both transcriptional activation and repression in response to starvation

Topoisomerases are required to release topological stress generated by RNA polymerase II (RNAPII) during transcription. Here, we show that in response to starvation, the complex of topoisomerase 3b (TOP3B) and TDRD3 can enhance not only transcriptional activation, but also repression, which mimics other topoisomerases that can also alter transcription in both directions. The genes enhanced by TOP3B-TDRD3 are enriched with long and highly-expressed ones, which are also preferentially stimulated by other topoisomerases, suggesting that different topoisomerases may recognize their targets through a similar mechanism. Specifically, human HCT116 cells individually inactivated for TOP3B, TDRD3 or TOP3B topoisomerase activity, exhibit similarly disrupted transcription for both starvation-activated genes (SAGs) and starvation-repressed genes (SRGs). Responding to starvation, both TOP3B-TDRD3 and the elongating form of RNAPII exhibit concomitantly increased binding to TOP3B-dependent SAGs, at binding sites that overlap. Notably, TOP3B inactivation decreases the binding of elongating RNAPII to TOP3B-dependent SAGs while increased it to SRGs. Furthermore, TOP3B-ablated cells display reduced transcription of several autophagy-associated genes and autophagy per se. Our data suggest that TOP3B-TDRD3 can promote both transcriptional activation and repression by regulating RNAPII distribution. In addition, the findings that it can facilitate autophagy may account for the shortened lifespan of Top3b-KO mice.

Transcription shapes 3D chromatin organization by interacting with loop extrusion

Cohesin folds mammalian interphase chromosomes by extruding the chromatin fiber into numerous loops. "Loop extrusion" can be impeded by chromatin-bound factors, such as CTCF, which generates characteristic and functional chromatin organization patterns. It has been proposed that transcription relocalizes or interferes with cohesin and that active promoters are cohesin loading sites. However, the effects of transcription on cohesin have not been reconciled with observations of active extrusion by cohesin. To determine how transcription modulates extrusion, we studied mouse cells in which we could alter cohesin abundance, dynamics, and localization by genetic "knockouts" of the cohesin regulators CTCF and Wapl. Through Hi-C experiments, we discovered intricate, cohesin-dependent contact patterns near active genes. Chromatin organization around active genes exhibited hallmarks of interactions between transcribing RNA polymerases (RNAPs) and extruding cohesins. These observations could be reproduced by polymer simulations in which RNAPs were moving barriers to extrusion that obstructed, slowed, and pushed cohesins. The simulations predicted that preferential loading of cohesin at promoters is inconsistent with our experimental data. Additional ChIP-seq experiments showed that the putative cohesin loader Nipbl is not predominantly enriched at promoters. Therefore, we propose that cohesin is not preferentially loaded at promoters and that the barrier function of RNAP accounts for cohesin accumulation at active promoters. Altogether, we find that RNAP is an extrusion barrier that is not stationary, but rather, translocates and relocalizes cohesin. Loop extrusion and transcription might interact to dynamically generate and maintain gene interactions with regulatory elements and shape functional genomic organization.

NetProphet 3: a machine learning framework for transcription factor network mapping and multi-omics integration

MOTIVATION

Many methods have been proposed for mapping the targets of transcription factors (TFs) from gene expression data. It is known that combining outputs from multiple methods can improve performance. To date, outputs have been combined by using either simplistic formulae, such as geometric mean, or carefully hand-tuned formulae that may not generalize well to new inputs. Finally, the evaluation of accuracy has been challenging due to the lack of genome-scale, ground-truth networks.

RESULTS

We developed NetProphet3, which combines scores from multiple analyses automatically, using a tree boosting algorithm trained on TF binding location data. We also developed three independent, genome-scale evaluation metrics. By these metrics, NetProphet3 is more accurate than other commonly used packages, including NetProphet 2.0, when gene expression data from direct TF perturbations are available. Furthermore, its integration mode can forge a consensus network from gene expression data and TF binding location data.

AVAILABILITY AND IMPLEMENTATION

All data and code are available at https://zenodo.org/record/7504131#.Y7Wu3i-B2x8.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Advances in Off-Target Detection for CRISPR-Based Genome Editing

The CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-based genome editing system exhibits marked potential for both gene editing and gene therapy, and its continuous improvement contributes to its great clinical potential. However, the largest hindrance to its application in clinical practice is the presence of off-target effects (OTEs). Thus, in addition to continuous optimization of the CRISPR system to reduce and eventually eliminate OTEs, further development of unbiased genome-wide detection of OTEs is key for its successful clinical application. This article summarizes detection strategies for OTEs of different CRISPR systems, to provide detailed guidance for the detection of OTEs in CRISPR-based genome editing.

Prevalent and dynamic binding of the cell cycle checkpoint kinase Rad53 to gene promoters

Replication of the genome must be coordinated with gene transcription and cellular metabolism, especially following replication stress in the presence of limiting deoxyribonucleotides. The Saccharomyces cerevisiae Rad53 (CHEK2 in mammals) checkpoint kinase plays a major role in cellular responses to DNA replication stress. Cell cycle regulated, genome-wide binding of Rad53 to chromatin was examined. Under replication stress, the kinase bound to sites of active DNA replication initiation and fork progression, but unexpectedly to the promoters of about 20% of genes encoding proteins involved in multiple cellular functions. Rad53 promoter binding correlated with changes in expression of a subset of genes. Rad53 promoter binding to certain genes was influenced by sequence-specific transcription factors and less by checkpoint signaling. However, in checkpoint mutants, untimely activation of late-replicating origins reduces the transcription of nearby genes, with concomitant localization of Rad53 to their gene bodies. We suggest that the Rad53 checkpoint kinase coordinates genome-wide replication and transcription under replication stress conditions.

SESAME-catalyzed H3T11 phosphorylation inhibits Dot1-catalyzed H3K79me3 to regulate autophagy and telomere silencing

The glycolytic enzyme, pyruvate kinase Pyk1 maintains telomere heterochromatin by phosphorylating histone H3T11 (H3pT11), which promotes SIR (silent information regulator) complex binding at telomeres and prevents autophagy-mediated Sir2 degradation. However, the exact mechanism of action for H3pT11 is poorly understood. Here, we report that H3pT11 directly inhibits Dot1-catalyzed H3K79 tri-methylation (H3K79me3) and uncover how this histone crosstalk regulates autophagy and telomere silencing. Mechanistically, Pyk1-catalyzed H3pT11 directly reduces the binding of Dot1 to chromatin and inhibits Dot1-catalyzed H3K79me3, which leads to transcriptional repression of autophagy genes and reduced autophagy. Despite the antagonism between H3pT11 and H3K79me3, they work together to promote the binding of SIR complex at telomeres to maintain telomere silencing. Furthermore, we identify Reb1 as a telomere-associated factor that recruits Pyk1-containing SESAME (Serine-responsive SAM-containing Metabolic Enzyme) complex to telomere regions to phosphorylate H3T11 and prevent the invasion of H3K79me3 from euchromatin into heterochromatin to maintain telomere silencing. Together, these results uncover a histone crosstalk and provide insights into dynamic regulation of silent heterochromatin and autophagy in response to cell metabolism.

Control of ribosomal RNA synthesis by hematopoietic transcription factors

Ribosomal RNAs (rRNAs) are the most abundant cellular RNAs, and their synthesis from rDNA repeats by RNA polymerase I accounts for the bulk of all transcription. Despite substantial variation in rRNA transcription rates across cell types, little is known about cell-type-specific factors that bind rDNA and regulate rRNA transcription to meet tissue-specific needs. Using hematopoiesis as a model system, we mapped about 2,200 ChIP-seq datasets for 250 transcription factors (TFs) and chromatin proteins to human and mouse rDNA and identified robust binding of multiple TF families to canonical TF motifs on rDNA. Using a 47S-FISH-Flow assay developed for nascent rRNA quantification, we demonstrated that targeted degradation of C/EBP alpha (CEBPA), a critical hematopoietic TF with conserved rDNA binding, caused rapid reduction in rRNA transcription due to reduced RNA Pol I occupancy. Our work identifies numerous potential rRNA regulators and provides a template for dissection of TF roles in rRNA transcription.

Building integrative functional maps of gene regulation

Every cell in the human body inherits a copy of the same genetic information. The three billion base pairs of DNA in the human genome, and the roughly 50 000 coding and non-coding genes they contain, must thus encode all the complexity of human development and cell and tissue type diversity. Differences in gene regulation, or the modulation of gene expression, enable individual cells to interpret the genome differently to carry out their specific functions. Here we discuss recent and ongoing efforts to build gene regulatory maps, which aim to characterize the regulatory roles of all sequences in a genome. Many researchers and consortia have identified such regulatory elements using functional assays and evolutionary analyses; we discuss the results, strengths and shortcomings of their approaches. We also discuss new techniques the field can leverage and emerging challenges it will face while striving to build gene regulatory maps of ever-increasing resolution and comprehensiveness.

Epigenetics of neural differentiation: Spotlight on enhancers

Neural induction, both in vivo and in vitro, includes cellular and molecular changes that result in phenotypic specialization related to specific transcriptional patterns. These changes are achieved through the implementation of complex gene regulatory networks. Furthermore, these regulatory networks are influenced by epigenetic mechanisms that drive cell heterogeneity and cell-type specificity, in a controlled and complex manner. Epigenetic marks, such as DNA methylation and histone residue modifications, are highly dynamic and stage-specific during neurogenesis. Genome-wide assessment of these modifications has allowed the identification of distinct non-coding regulatory regions involved in neural cell differentiation, maturation, and plasticity. Enhancers are short DNA regulatory regions that bind transcription factors (TFs) and interact with gene promoters to increase transcriptional activity. They are of special interest in neuroscience because they are enriched in neurons and underlie the cell-type-specificity and dynamic gene expression profiles. Classification of the full epigenomic landscape of neural subtypes is important to better understand gene regulation in brain health and during diseases. Advances in novel next-generation high-throughput sequencing technologies, genome editing, Genome-wide association studies (GWAS), stem cell differentiation, and brain organoids are allowing researchers to study brain development and neurodegenerative diseases with an unprecedented resolution. Herein, we describe important epigenetic mechanisms related to neurogenesis in mammals. We focus on the potential roles of neural enhancers in neurogenesis, cell-fate commitment, and neuronal plasticity. We review recent findings on epigenetic regulatory mechanisms involved in neurogenesis and discuss how sequence variations within enhancers may be associated with genetic risk for neurological and psychiatric disorders.

Transcription factors perform a 2-step search of the nucleus

Transcription factors regulate gene expression by binding to regulatory DNA and recruiting regulatory protein complexes. The DNA-binding and protein-binding functions of transcription factors are traditionally described as independent functions performed by modular protein domains. Here, I argue that genome binding can be a 2-part process with both DNA-binding and protein-binding steps, enabling transcription factors to perform a 2-step search of the nucleus to find their appropriate binding sites in a eukaryotic genome. I support this hypothesis with new and old results in the literature, discuss how this hypothesis parsimoniously resolves outstanding problems, and present testable predictions.

Capturing Protein-Nucleic Acid Interactions by High-Intensity Laser-Induced Covalent Crosslinking

Interactions of DNA with structural proteins such as histones, regulatory proteins, and enzymes play a crucial role in major cellular processes such as transcription, replication and repair. The in vivo mapping and characterization of the binding sites of the involved biomolecules are of primary importance for a better understanding of genomic deployment that is implicated in tissue and developmental stage-specific gene expression regulation. The most powerful and commonly used approach to date is immunoprecipitation of chemically cross-linked chromatin (XChIP) coupled with sequencing analysis (ChIP-seq). While the resolution and the sensitivity of the high-throughput sequencing techniques have been constantly improved little progress has been achieved in the crosslinking step. Because of its low efficiency the use of the conventional UVC lamps remains very limited while the formaldehyde method was established as the "gold standard" crosslinking agent. Efficient biphotonic crosslinking of directly interacting nucleic acid-protein complexes by a single short UV laser pulse has been introduced as an innovative technique for overcoming limitations of conventionally used chemical and photochemical approaches. In this survey, the main available methods including the laser approach are critically reviewed for their ability to generate DNA-protein crosslinks in vitro model systems and cells.

GoPeaks: histone modification peak calling for CUT&Tag

Genome-wide mapping of histone modifications is critical to understanding transcriptional regulation. CUT&Tag is a new method for profiling histone modifications, offering improved sensitivity and decreased cost compared with ChIP-seq. Here, we present GoPeaks, a peak calling method specifically designed for histone modification CUT&Tag data. We compare the performance of GoPeaks against commonly used peak calling algorithms to detect histone modifications that display a range of peak profiles and are frequently used in epigenetic studies. We find that GoPeaks robustly detects genome-wide histone modifications and, notably, identifies a substantial number of H3K27ac peaks with improved sensitivity compared to other standard algorithms.

Genomic Occupancy of the Bromodomain Protein Bdf3 Is Dynamic during Differentiation of African Trypanosomes from Bloodstream to Procyclic Forms

Trypanosoma brucei, the causative agent of human and animal African trypanosomiasis, cycles between a mammalian host and a tsetse fly vector. The parasite undergoes huge changes in morphology and metabolism during adaptation to each host environment. These changes are reflected in the different transcriptomes of parasites living in each host. However, it remains unclear whether chromatin-interacting proteins help mediate these changes. Bromodomain proteins localize to transcription start sites in bloodstream parasites, but whether the localization of bromodomain proteins changes as parasites differentiate from bloodstream to insect stages remains unknown. To address this question, we performed cleavage under target and release using nuclease (CUT&RUN) against bromodomain protein 3 (Bdf3) in parasites differentiating from bloodstream to insect forms. We found that Bdf3 occupancy at most loci increased at 3 h following onset of differentiation and decreased thereafter. A number of sites with increased bromodomain protein occupancy lie proximal to genes with altered transcript levels during differentiation, such as procyclins, procyclin-associated genes, and invariant surface glycoproteins. Most Bdf3-occupied sites are observed throughout differentiation. However, one site appears de novo during differentiation and lies proximal to the procyclin gene locus housing genes essential for remodeling surface proteins following transition to the insect stage. These studies indicate that occupancy of chromatin-interacting proteins is dynamic during life cycle stage transitions and provide the groundwork for future studies on the effects of changes in bromodomain protein occupancy. Additionally, the adaptation of CUT&RUN for Trypanosoma brucei provides other researchers with an alternative to chromatin immunoprecipitation (ChIP). IMPORTANCE The parasite Trypanosoma brucei is the causative agent of human and animal African trypanosomiasis (sleeping sickness). Trypanosomiasis, which affects humans and cattle, is fatal if untreated. Existing drugs have significant side effects. Thus, these parasites impose a significant human and economic burden in sub-Saharan Africa, where trypanosomiasis is endemic. T. brucei cycles between the mammalian host and a tsetse fly vector, and parasites undergo huge changes in morphology and metabolism to adapt to different hosts. Here, we show that DNA-interacting bromodomain protein 3 (Bdf3) shows changes in occupancy at its binding sites as parasites transition from the bloodstream to the insect stage. Additionally, a new binding site appears near the locus responsible for remodeling of parasite surface proteins during transition to the insect stage. Understanding the mechanisms behind host adaptation is important for understanding the life cycle of the parasite.

Function and dynamics of the Mediator complex: novel insights and new frontiers

The Mediator complex was discovered in the early 1990s as a biochemically fractionated factor from yeast extracts that was necessary for activator-stimulated transcriptional activation to be observed in in vitro transcription assays. The structure of this large, multi-protein complex is now understood in great detail, and novel genetic approaches have provided rich insights into its dynamics during transcriptional activation and the mechanism by which it facilitates activated transcription. Here I review recent findings and unanswered questions regarding Mediator dynamics, the roles of individual subunits, and differences between its function in yeast and metazoan cells.

Virtual ChIP-seq: predicting transcription factor binding by learning from the transcriptome

Existing methods for computational prediction of transcription factor (TF) binding sites evaluate genomic regions with similarity to known TF sequence preferences. Most TF binding sites, however, do not resemble known TF sequence motifs, and many TFs are not sequence-specific. We developed Virtual ChIP-seq, which predicts binding of individual TFs in new cell types, integrating learned associations with gene expression and binding, TF binding sites from other cell types, and chromatin accessibility data in the new cell type. This approach outperforms methods that predict TF binding solely based on sequence preference, predicting binding for 36 TFs (MCC>0.3).

nMOWChIP-seq: low-input genome-wide mapping of non-histone targets

Genome-wide profiling of interactions between genome and various functional proteins is critical for understanding regulatory processes involved in development and diseases. Conventional assays require a large number of cells and high-quality data on tissue samples are scarce. Here we optimized a low-input chromatin immunoprecipitation followed by sequencing (ChIP-seq) technology for profiling RNA polymerase II (Pol II), transcription factor (TF), and enzyme binding at the genome scale. The new approach produces high-quality binding profiles using 1,000-50,000 cells. We used the approach to examine the binding of Pol II and two TFs (EGR1 and MEF2C) in cerebellum and prefrontal cortex of mouse brain and found that their binding profiles are highly reflective of the functional differences between the two brain regions. Our analysis reveals the potential for linking genome-wide TF or Pol II profiles with neuroanatomical origins of brain cells.

Biomimetic mimicry of formaldehyde-induced DNA-protein crosslinks in the confined space of a metal-organic framework

DNA-protein crosslinks (DPCs) are highly toxic DNA lesions induced by crosslinking agents such as formaldehyde (HCHO). Building artificial models to simulate the crosslinking process would advance our understanding of the underlying mechanisms and therefore develop coping strategies accordingly. Herein we report the design and synthesis of a Zn-based metal-organic framework with mixed ligands of 2,6-diaminopurine and amine-functionalized dicarboxylate, representing DNA and protein residues, respectively. Combined characterization techniques allow us to demonstrate the unusual efficiency of HCHO-crosslinking within the confined space of the titled MOF. Particularly, in situ single-crystal X-ray diffraction studies reveal a sequential methylene-knitting process upon HCHO addition, along with strong fluorescence that was not interfered with by other metabolites, glycine, and Tris. This work has successfully constructed a purine-based metal-organic framework with unoccupied Watson-Crick sites, serving as a crystalline model for HCHO-induced DPCs by mimicking the confinement effect of protein/DNA interactions.

MafB, WDR77, and ß-catenin interact with each other and have similar genome association profiles

MafB (a bZIP transcription factor), ß-catenin (the ultimate target of the Wnt signal transduction pathway that acts as a transcriptional co-activator of LEF/TCF proteins), and WDR77 (a transcriptional co-activator of multiple hormone receptors) are important for breast cellular transformation. Unexpectedly, these proteins interact directly with each other, and they have similar genomic binding profiles. Furthermore, while some of these common target sites coincide with those bound by LEF/TCF, the majority are located just downstream of transcription initiation sites at a position near paused RNA polymerase (Pol II) and the +1 nucleosome. Occupancy levels of these factors at these promoter-proximal sites are strongly correlated with the level of paused Pol II and transcriptional activity.

Regulation of the Alternative Neural Transcriptome by ELAV/Hu RNA Binding Proteins

The process of alternative polyadenylation (APA) generates multiple 3' UTR isoforms for a given locus, which can alter regulatory capacity and on occasion change coding potential. APA was initially characterized for a few genes, but in the past decade, has been found to be the rule for metazoan genes. While numerous differences in APA profiles have been catalogued across genetic conditions, perturbations, and diseases, our knowledge of APA mechanisms and biology is far from complete. In this review, we highlight recent findings regarding the role of the conserved ELAV/Hu family of RNA binding proteins (RBPs) in generating the broad landscape of lengthened 3' UTRs that is characteristic of neurons. We relate this to their established roles in alternative splicing, and summarize ongoing directions that will further elucidate the molecular strategies for neural APA, the in vivo functions of ELAV/Hu RBPs, and the phenotypic consequences of these regulatory paradigms in neurons.