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

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.

Chromatin endogenous cleavage provides a global view of 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.

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.

Sequence characteristics and an accurate 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 seemingly 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.

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.

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.

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.

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.

CUT&RUN Profiling of the Budding Yeast Epigenome

Mapping the epigenome is key to describe the relationship between chromatin landscapes and the control of DNA-based cellular processes such as transcription. Cleavage under targets and release using nuclease (CUT&RUN) is an in situ chromatin profiling strategy in which controlled cleavage by antibody-targeted Micrococcal Nuclease solubilizes specific protein-DNA complexes for paired-end DNA sequencing. When applied to budding yeast, CUT&RUN profiling yields precise genome-wide maps of histone modifications, histone variants, transcription factors, and ATP-dependent chromatin remodelers, while avoiding cross-linking and solubilization issues associated with the most commonly used chromatin profiling technique Chromatin Immunoprecipitation (ChIP). Furthermore, targeted chromatin complexes cleanly released by CUT&RUN can be used as input for a subsequent native immunoprecipitation step (CUT&RUN.ChIP) to simultaneously map two epitopes in single molecules genome-wide. The intrinsically low background and high resolution of CUT&RUN and CUT&RUN.ChIP allows for identification of transient genomic features such as dynamic nucleosome-remodeling intermediates. Starting from cells, one can perform CUT&RUN or CUT&RUN.ChIP and obtain purified DNA for sequencing library preparation in 2 days.

Automated CUT & RUN Using the KingFisher Duo Prime

Elucidating the biological function of histone methyltransferases requires knowledge of the genomic sites at which they act. CUT&RUN represents a valuable alternative to chromatin immunoprecipitation for the mapping of histone methylation patterns, generally producing results of equivalent quality while requiring less sequencing depth, less starting material and less effort. Automated CUT&RUN procedures have been developed to further facilitate chromatin profiling. Here we describe our automated CUT&RUN protocol using the Thermo Fisher KingFisher Duo Prime system.

Master lineage transcription factors anchor trans mega transcriptional complexes at highly accessible enhancer sites to promote long-range chromatin clustering and transcription of distal target genes

The term 'super enhancers' (SE) has been widely used to describe stretches of closely localized enhancers that are occupied collectively by large numbers of transcription factors (TFs) and co-factors, and control the transcription of highly-expressed genes. Through integrated analysis of >600 DNase-seq, ChIP-seq, GRO-seq, STARR-seq, RNA-seq, Hi-C and ChIA-PET data in five human cancer cell lines, we identified a new class of autonomous SEs (aSEs) that are excluded from classic SE calls by the widely used Rank Ordering of Super-Enhancers (ROSE) method. TF footprint analysis revealed that compared to classic SEs and regular enhancers, aSEs are tightly bound by a dense array of master lineage TFs, which serve as anchors to recruit additional TFs and co-factors in trans. In addition, aSEs are preferentially enriched for Cohesins, which likely involve in stabilizing long-distance interactions between aSEs and their distal target genes. Finally, we showed that aSEs can be reliably predicted using a single DNase-seq data or combined with Mediator and/or P300 ChIP-seq. Overall, our study demonstrates that aSEs represent a unique class of functionally important enhancer elements that distally regulate the transcription of highly expressed genes.

Mediator dynamics during heat shock in budding yeast

The Mediator complex is central to transcription by RNA polymerase II (Pol II) in eukaryotes. In budding yeast (Saccharomyces cerevisiae), Mediator is recruited by activators and associates with core promoter regions, where it facilitates preinitiation complex (PIC) assembly, only transiently before Pol II escape. Interruption of the transcription cycle by inactivation or depletion of Kin28 inhibits Pol II escape and stabilizes this association. However, Mediator occupancy and dynamics have not been examined on a genome-wide scale in yeast grown in nonstandard conditions. Here we investigate Mediator occupancy following heat shock or CdCl₂ exposure, with and without depletion of Kin28. We find that Pol II occupancy shows similar dependence on Mediator under normal and heat shock conditions. However, although Mediator association increases at many genes upon Kin28 depletion under standard growth conditions, little or no increase is observed at most genes upon heat shock, indicating a more stable association of Mediator after heat shock. Unexpectedly, Mediator remains associated upstream of the core promoter at genes repressed by heat shock or CdCl₂ exposure whether or not Kin28 is depleted, suggesting that Mediator is recruited by activators but is unable to engage PIC components at these repressed targets. This persistent association is strongest at promoters that bind the HMGB family member Hmo1, and is reduced but not eliminated in hmo1Δ yeast. Finally, we show a reduced dependence on PIC components for Mediator occupancy at promoters after heat shock, further supporting altered dynamics or stronger engagement with activators under these conditions.

Characterizing batch effects and binding site-specific variability in ChIP-seq data

Multiple sources of variability can bias ChIP-seq data toward inferring transcription factor (TF) binding profiles. As ChIP-seq datasets increase in public repositories, it is now possible and necessary to account for complex sources of variability in ChIP-seq data analysis. We find that two types of variability, the batch effects by sequencing laboratories and differences between biological replicates, not associated with changes in condition or state, vary across genomic sites. This implies that observed differences between samples from different conditions or states, such as cell-type, must be assessed statistically, with an understanding of the distribution of obscuring noise. We present a statistical approach that characterizes both differences of interests and these source of variability through the parameters of a mixed effects model. We demonstrate the utility of our approach on a CTCF binding dataset composed of 211 samples representing 90 different cell-types measured across three different laboratories. The results revealed that sites exhibiting large variability were associated with sequence characteristics such as GC-content and low complexity. Finally, we identified TFs associated with high-variance CTCF sites using TF motifs documented in public databases, pointing the possibility of these being false positives if the sources of variability are not properly accounted for.

Biologically relevant transfer learning improves transcription factor binding prediction

BACKGROUND

Deep learning has proven to be a powerful technique for transcription factor (TF) binding prediction but requires large training datasets. Transfer learning can reduce the amount of data required for deep learning, while improving overall model performance, compared to training a separate model for each new task.

RESULTS

We assess a transfer learning strategy for TF binding prediction consisting of a pre-training step, wherein we train a multi-task model with multiple TFs, and a fine-tuning step, wherein we initialize single-task models for individual TFs with the weights learned by the multi-task model, after which the single-task models are trained at a lower learning rate. We corroborate that transfer learning improves model performance, especially if in the pre-training step the multi-task model is trained with biologically relevant TFs. We show the effectiveness of transfer learning for TFs with ~ 500 ChIP-seq peak regions. Using model interpretation techniques, we demonstrate that the features learned in the pre-training step are refined in the fine-tuning step to resemble the binding motif of the target TF (i.e., the recipient of transfer learning in the fine-tuning step). Moreover, pre-training with biologically relevant TFs allows single-task models in the fine-tuning step to learn useful features other than the motif of the target TF.

CONCLUSIONS

Our results confirm that transfer learning is a powerful technique for TF binding prediction.

Anomaly detection in genomic catalogues using unsupervised multi-view autoencoders

Background

Accurate identification of Transcriptional Regulator binding locations is essential for analysis of genomic regions, including Cis Regulatory Elements. The customary NGS approaches, predominantly ChIP-Seq, can be obscured by data anomalies and biases which are difficult to detect without supervision.

Results

Here, we develop a method to leverage the usual combinations between many experimental series to mark such atypical peaks. We use deep learning to perform a lossy compression of the genomic regions’ representations with multiview convolutions. Using artificial data, we show that our method correctly identifies groups of correlating series and evaluates CRE according to group completeness. It is then applied to the ReMap database’s large volume of curated ChIP-seq data. We show that peaks lacking known biological correlators are singled out and less confirmed in real data. We propose normalization approaches useful in interpreting black-box models.

Conclusion

Our approach detects peaks that are less corroborated than average. It can be extended to other similar problems, and can be interpreted to identify correlation groups. It is implemented in an open-source tool called atyPeak.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04359-2.

Defining Proximity Proteomics of Histone Modifications by Antibody-mediated Protein A-APEX2 Labeling

Proximity labeling catalyzed by promiscuous enzymes, such as APEX2, has emerged as a powerful approach to characterize multiprotein complexes and protein–protein interactions. However, current methods depend on the expression of exogenous fusion proteins and cannot be applied to identify proteins surrounding post-translationally modified proteins. To address this limitation, we developed a new method to label proximal proteins of interest by antibody-mediated protein A-ascorbate peroxidase 2 (pA-APEX2) labeling (AMAPEX). In this method, a modified protein is bound in situ by a specific antibody, which then tethers a pA-APEX2 fusion protein. Activation of APEX2 labels the nearby proteins with biotin; the biotinylated proteins are then purified using streptavidin beads and identified by mass spectrometry. We demonstrated the utility of this approach by profiling the proximal proteins of histone modifications including H3K27me3, H3K9me3, H3K4me3, H4K5ac, and H4K12ac, as well as verifying the co-localization of these identified proteins with bait proteins by published ChIP-seq analysis and nucleosome immunoprecipitation. Overall, AMAPEX is an efficient method to identify proteins that are proximal to modified histones.

Genome-wide variability in recombination activity is associated with meiotic chromatin organization

Recombination enables reciprocal exchange of genomic information between parental chromosomes and successful segregation of homologous chromosomes during meiosis. Errors in this process lead to negative health outcomes, whereas variability in recombination rate affects genome evolution. In mammals, most crossovers occur in hotspots defined by PRDM9 motifs, although PRDM9 binding peaks are not all equally hot. We hypothesize that dynamic patterns of meiotic genome folding are linked to recombination activity. We apply an integrative bioinformatics approach to analyze how three-dimensional (3D) chromosomal organization during meiosis relates to rates of double-strand-break (DSB) and crossover (CO) formation at PRDM9 binding peaks. We show that active, spatially accessible genomic regions during meiotic prophase are associated with DSB-favored loci, which further adopt a transient locally active configuration in early prophase. Conversely, crossover formation is depleted among DSBs in spatially accessible regions during meiotic prophase, particularly within gene bodies. We also find evidence that active chromatin regions have smaller average loop sizes in mammalian meiosis. Collectively, these findings establish that differences in chromatin architecture along chromosomal axes are associated with variable recombination activity. We propose an updated framework describing how 3D organization of brush-loop chromosomes during meiosis may modulate recombination.