BANP opens chromatin and activates CpG-island-regulated genes

Q-rich activation domains: flexible 'rulers' for transcription start site selection?

Recent findings broadened the function of RNA polymerase II (Pol II) proximal promoter motifs from quantitative regulators of transcription to important determinants of transcription start site (TSS) position. These motifs are recognized by transcription factors (TFs) that we propose to term 'ruler' TFs (rTFs), such as NRF1, NF-Y, YY1, ZNF143, BANP, and members of the SP, ETS, and CRE families, sharing as a common feature a glutamine-rich (Q-rich) effector domain also enriched in valine, isoleucine, and threonine (QVIT-rich). We propose that rTFs guide TSS location by constraining the position of the pre-initiation complex (PIC) during its promoter recognition phase through a specialized, and still enigmatic, class of activation domains.

DNA-binding proteins from MBD through ZF to BEN: recognition of cytosine methylation status by one arginine with two conformations

Maintenance methylation, of palindromic CpG dinucleotides at DNA replication forks, is crucial for the faithful mitotic inheritance of genomic 5-methylcytosine (5mC) methylation patterns. MBD proteins use two arginine residues to recognize symmetrically-positioned methyl groups in fully-methylated 5mCpG/5mCpG and 5mCpA/TpG dinucleotides. In contrast, C2H2 zinc finger (ZF) proteins recognize CpG and CpA, whether methylated or not, within longer specific sequences in a site- and strand-specific manner. Unmethylated CpG sites, often within CpG island (CGI) promoters, need protection by protein factors to maintain their hypomethylated status. Members of the BEN domain proteins bind CGCG or CACG elements within CGIs to regulate gene expression. Despite their overall structural diversity, MBD, ZF and BEN proteins all use arginine residues to recognize guanine, adopting either a 'straight-on' or 'oblique' conformation. The straight-on conformation accommodates a methyl group in the (5mC/T)pG dinucleotide, while the oblique conformation can clash with the methyl group of 5mC, leading to preferential binding of unmethylated sequences.

Structural basis of DNA recognition by BEN domain proteins reveals a role for oligomerization in unmethylated DNA selection by BANP

The BEN domain is a newly discovered type of DNA-binding domain that exists in a variety of species. There are nine BEN domain-containing proteins in humans, and most have been shown to have chromatin-related functions. NACC1 preferentially binds to CATG motif-containing sequences and functions primarily as a transcriptional coregulator. BANP and BEND3 preferentially bind DNA bearing unmethylated CpG motifs, and they function as CpG island-binding proteins. To date, the DNA recognition mechanism of quite a few of these proteins remains to be determined. In this study, we solved the crystal structures of the BEN domains of NACC1 and BANP in complex with their cognate DNA substrates. We revealed the details of DNA binding by these BEN domain proteins and unexpectedly revealed that oligomerization is required for BANP to select unmethylated CGCG motif-containing DNA substrates. Our study clarifies the controversies surrounding DNA recognition by BANP and demonstrates a new mechanism by which BANP selects unmethylated CpG motifs and functions as a CpG island-binding protein. This understanding will facilitate further exploration of the physiological functions of the BEN domain proteins in the future.

Evolutionary Innovations in Conserved Regulatory Elements Associate With Developmental Genes in Mammals

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Genetic variation within enhancer sequences is an important contributor to phenotypic variation including evolutionary adaptations and human disease. Certain genes and pathways may be more prone to regulatory evolution than others, with different patterns across diverse organisms, but whether such patterns exist has not been investigated at a sufficient scale. To address this question, we identified signatures of accelerated sequence evolution in conserved enhancer elements throughout the mammalian phylogeny at an unprecedented scale. While different genes and pathways were enriched for regulatory evolution in different parts of the tree, we found a striking overall pattern of pleiotropic genes involved in gene regulatory and developmental processes being enriched for accelerated enhancer evolution. These genes were connected to more enhancers than other genes, which was the basis for having an increased amount of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. Detailed study of one acceleration event in an enhancer of HES1 revealed that sequence evolution led to a new activity domain in the developing limb that emerged concurrently with the evolution of digit reduction in hoofed mammals. Our results provide evidence that enhancer evolution has been a frequent contributor to regulatory innovation at conserved developmental signaling genes in mammals.

On the identification of differentially-active transcription factors from ATAC-seq data

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (a chromVAR-limma workflow with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs-in our case NFkB-at the protein level.

Genome access is transcription factor-specific and defined by nucleosome position

Mammalian gene expression is controlled by transcription factors (TFs) that engage sequence motifs in a chromatinized genome, where nucleosomes can restrict DNA access. Yet, how nucleosomes affect individual TFs remains unclear. Here, we measure the ability of over one hundred TF motifs to recruit TFs in a defined chromosomal locus in mouse embryonic stem cells. This identifies a set sufficient to enable the binding of TFs with diverse tissue specificities, functions, and DNA-binding domains. These chromatin-competent factors are further classified when challenged to engage motifs within a highly phased nucleosome. The pluripotency factors OCT4-SOX2 preferentially engage non-nucleosomal and entry-exit motifs, but not nucleosome-internal sites, a preference that also guides binding genome wide. By contrast, factors such as BANP, REST, or CTCF engage throughout, causing nucleosomal displacement. This supports that TFs vary widely in their sensitivity to nucleosomes and that genome access is TF specific and influenced by nucleosome position in the cell.

CpG Methylation of Receptor Activator NF-κB (RANK) Gene Promoter Region Delineates Senescence-Related Decrease of RANK Gene Expression

While the rapid decrease in estrogen is well known as the main cause of postmenopausal osteoporosis in women, the precise pathogenesis of senile osteoporosis in the elderly regardless of gender is largely unknown. The age-related epigenetic regulation of receptor activator NF-κB (RANK) gene expression was investigated with the use of a high-passaged mouse osteoclast progenitor cell line, RAW264.7, as an in vitro model of aging. In the RAW264.7 cells after repeated passages, receptor RANK expression was downregulated, resulting in decreased soluble RANK ligand (sRANKL)-induced osteoclastogenesis, expression of tartrate-resistant acid phosphatase-5b (TRAcP) and cathepsin K (CTSK). Methylation-specific PCR and bisulfite mapping revealed hypermethylation of CpG-loci located in the RANK gene promoter in multiple-passaged cells. ICON probe-mediated in situ assessment of methylated-cytosine at the CpG loci revealed an increase in the percentage of methylated RAW264.7 cells in the RANK gene in a passage-dependent manner. Conversely, upon treatment with demethylating agent 5-aza-2-deoxycytidine (5-aza-dC), high-passaged RAW264.7 cells displayed restored expression of the RANK gene, osteoclastogenesis, TRAcP and CTSK. Ex vivo cultures of splenic macrophages from young (10.5 W) and aged (12 M) mice also showed that CpG methylation was predominant in the aged animals, resulting in reduced RANK expression and osteoclastogenesis. Reduced RANK expression by age-related accumulation of DNA methylation, albeit in a limited population of osteoclast precursor cells, might be, at least in part, indicative of low-turnover bone characteristic of senile osteoporosis.

On the identification of differentially-active transcription factors from ATAC-seq data

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (a chromVAR-limma workflow with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs - in our case NFkB - at the protein level.

guidedNOMe-seq quantifies chromatin states at single allele resolution for hundreds of custom regions in parallel

Since the introduction of next generation sequencing technologies, the field of epigenomics has evolved rapidly. However, most commonly used assays are enrichment-based methods and thus only semi-quantitative. Nucleosome occupancy and methylome sequencing (NOMe-seq) allows for quantitative inference of chromatin states with single locus resolution, but this requires high sequencing depth and is therefore prohibitively expensive to routinely apply to organisms with large genomes. To overcome this limitation, we introduce guidedNOMe-seq, where we combine NOMe profiling with large scale sgRNA synthesis and Cas9-mediated region-of-interest (ROI) liberation. To facilitate quantitative comparisons between multiple samples, we additionally develop an R package to standardize differential analysis of any type of NOMe-seq data. We extensively benchmark guidedNOMe-seq in a proof-of-concept study, dissecting the interplay of ChAHP and CTCF on chromatin. In summary we present a cost-effective, scalable, and customizable target enrichment extension to the existing NOMe-seq protocol allowing genome-scale quantification of nucleosome occupancy and transcription factor binding at single allele resolution.

Deciphering lineage specification during early embryogenesis in mouse gastruloids using multilayered proteomics

Gastrulation is a critical stage in embryonic development during which the germ layers are established. Advances in sequencing technologies led to the identification of gene regulatory programs that control the emergence of the germ layers and their derivatives. However, proteome-based studies of early mammalian development are scarce. To overcome this, we utilized gastruloids and a multilayered mass spectrometry-based proteomics approach to investigate the global dynamics of (phospho) protein expression during gastruloid differentiation. Our findings revealed many proteins with temporal expression and unique expression profiles for each germ layer, which we also validated using single-cell proteomics technology. Additionally, we profiled enhancer interaction landscapes using P300 proximity labeling, which revealed numerous gastruloid-specific transcription factors and chromatin remodelers. Subsequent degron-based perturbations combined with single-cell RNA sequencing (scRNA-seq) identified a critical role for ZEB2 in mouse and human somitogenesis. Overall, this study provides a rich resource for developmental and synthetic biology communities endeavoring to understand mammalian embryogenesis.

Position-dependent function of human sequence-specific transcription factors

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

Structural homology-based identification of BEN domain proteins in Poxviruses

BEND family transcription factors directly interact with DNA through BEN domains and have been found across metazoan species. Interestingly, certain insect and mammalian viruses have also hijacked Bend genes into their genome. However, the phylogenetic classification and evolution of these viral BEN domains remain unclear. Building on our previous finding that in silico method accurately determine the 3D model of BEN domains, we used AlphaFold2 to predict the tertiary structures of poxviral BEN domains for comprehensive homologous comparison. We revealed that the majority of poxviral BEN modules exhibit characteristics of type II BEN. Additionally, electrostatic surface potential analysis found various poxviral BEN domains, including the first BEN of OPG067 in Orthopoxvirus, the third BEN of OPG067 in Yatapoxvirus and the third BEN of MC036R in MCV, have positively charged protein surfaces, indicating a structural basis for DNA loading. Notably, MC036R shares structural resemblance with human BEND3, as they both contain four BEN domains and an intrinsically disordered region. In summary, our discoveries provide deeper insights into the functional roles of BEN proteins within poxviruses.

Systematic assessment of ISWI subunits shows that NURF creates local accessibility for CTCF

Catalytic activity of the imitation switch (ISWI) family of remodelers is critical for nucleosomal organization and DNA binding of certain transcription factors, including the insulator protein CTCF. Here we define the contribution of individual subcomplexes by deriving a panel of isogenic mouse stem cell lines, each lacking one of six ISWI accessory subunits. Individual deletions of subunits of either CERF, RSF, ACF, WICH or NoRC subcomplexes only moderately affect the chromatin landscape, while removal of the NURF-specific subunit BPTF leads to a strong reduction in chromatin accessibility and SNF2H ATPase localization around CTCF sites. This affects adjacent nucleosome occupancy and CTCF binding. At a group of sites with reduced chromatin accessibility, CTCF binding persists but cohesin occupancy is reduced, resulting in decreased insulation. These results suggest that CTCF binding can be separated from its function as an insulator in nuclear organization and identify a specific role for NURF in mediating SNF2H localization and chromatin opening at bound CTCF sites.

Exploring affinity chromatography in proteomics: A comprehensive review

Over the past decades, the proteomics field has undergone rapid growth. Progress in mass spectrometry and bioinformatics, together with separation methods, has brought many innovative approaches to the study of the molecular biology of the cell. The potential of affinity chromatography was recognized immediately after its first application in proteomics, and since that time, it has become one of the cornerstones of many proteomic protocols. Indeed, this chromatographic technique exploiting the specific binding between two molecules has been employed for numerous purposes, from selective removal of interfering (over)abundant proteins or enrichment of scarce biomarkers in complex biological samples to mapping the post-translational modifications and protein interactions with other proteins, nucleic acids or biologically active small molecules. This review presents a comprehensive survey of this versatile analytical tool in current proteomics. To navigate the reader, the haphazard space of affinity separations is classified according to the experiment's aims and the separated molecule's nature. Different types of available ligands and experimental strategies are discussed in further detail for each of the mentioned procedures.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

Integrative cross-omics and cross-context analysis elucidates molecular links underlying genetic effects on complex traits

Genetic effects on functionally related 'omic' traits often co-occur in relevant cellular contexts, such as tissues. Motivated by the multi-tissue methylation quantitative trait loci (mQTLs) and expression QTLs (eQTLs) analysis, we propose X-ING (Cross-INtegrative Genomics) for cross-omics and cross-context integrative analysis. X-ING takes as input multiple matrices of association statistics, each obtained from different omics data types across multiple cellular contexts. It models the latent binary association status of each statistic, captures the major association patterns among omics data types and contexts, and outputs the posterior mean and probability for each input statistic. X-ING enables the integration of effects from different omics data with varying effect distributions. In the multi-tissue cis-association analysis, X-ING shows improved detection and replication of mQTLs by integrating eQTL maps. In the trans-association analysis, X-ING reveals an enrichment of trans-associations in many disease/trait-relevant tissues.

Evolutionary innovation in conserved regulatory elements across the mammalian tree of life

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Evolution of enhancer sequences can alter target gene expression without causing detrimental misexpression in other contexts. It has long been thought that this modularity allows evolutionary changes in enhancers to escape pleiotropic constraints, which is especially important for evolutionary constrained developmental patterning genes. However, there is still little data supporting this hypothesis. Here we identified signatures of accelerated evolution in conserved enhancer elements across the mammalian phylogeny. We found that pleiotropic genes involved in gene regulatory and developmental processes were enriched for accelerated sequence evolution within their enhancer elements. These genes were associated with an excess number of enhancers compared to other genes, and due to this they exhibit a substantial degree of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. We studied one acceleration event in depth and found that its sequence evolution led to the emergence of a new enhancer activity domain that may be involved in the evolution of digit reduction in hoofed mammals. Our results provide tangible evidence that enhancer evolution has been a frequent contributor to modifications involving constrained developmental signaling genes in mammals.

Loss of the DYRK1A Protein Kinase Results in the Reduction in Ribosomal Protein Gene Expression, Ribosome Mass and Reduced Translation

Ribosomal proteins (RPs) are evolutionary conserved proteins that are essential for protein translation. RP expression must be tightly regulated to ensure the appropriate assembly of ribosomes and to respond to the growth demands of cells. The elements regulating the transcription of RP genes (RPGs) have been characterized in yeast and Drosophila, yet how cells regulate the production of RPs in mammals is less well understood. Here, we show that a subset of RPG promoters is characterized by the presence of the palindromic TCTCGCGAGA motif and marked by the recruitment of the protein kinase DYRK1A. The presence of DYRK1A at these promoters is associated with the enhanced binding of the TATA-binding protein, TBP, and it is negatively correlated with the binding of the GABP transcription factor, establishing at least two clusters of RPGs that could be coordinately regulated. However, DYRK1A silencing leads to a global reduction in RPGs mRNAs, pointing at DYRK1A activities beyond those dependent on its chromatin association. Significantly, cells in which DYRK1A is depleted have reduced RP levels, fewer ribosomes, reduced global protein synthesis and a smaller size. We therefore propose a novel role for DYRK1A in coordinating the expression of genes encoding RPs, thereby controlling cell growth in mammals.

A genetic screen identifies BEND3 as a regulator of bivalent gene expression and global DNA methylation

Epigenetic mechanisms are essential to establish and safeguard cellular identities in mammals. They dynamically regulate the expression of genes, transposable elements and higher-order chromatin structures. Consequently, these chromatin marks are indispensable for mammalian development and alterations often lead to disease, such as cancer. Bivalent promoters are especially important during differentiation and development. Here we used a genetic screen to identify new regulators of a bivalent repressed gene. We identify BEND3 as a regulator of hundreds of bivalent promoters, some of which it represses, and some of which it activates. We show that BEND3 is recruited to a CpG-containg consensus site that is present in multiple copies in many bivalent promoters. Besides having direct effect on the promoters it binds, the loss of BEND3 leads to genome-wide gains of DNA methylation, which are especially marked at regions normally protected by the TET enzymes. DNA hydroxymethylation is reduced in Bend3 mutant cells, possibly as consequence of altered gene expression leading to diminished alpha-ketoglutarate production, thus lowering TET activity. Our results clarify the direct and indirect roles of an important chromatin regulator, BEND3, and, more broadly, they shed light on the regulation of bivalent promoters.

Differential responses to maternal diabetes in embryo and visceral yolk sac

Introduction: Maternal diabetes during pregnancy is well known to be associated with a higher risk for structural birth defects in the offspring. Recent searches for underlying mechanisms have largely focused on aberrant processes in the embryo itself, although prior research in rodent models implicated dysfunction also of the visceral yolk sac. The objective of our research was to investigate both tissues within the conceptus simultaneously. Methods: We conducted unbiased transcriptome profiling by RNA sequencing on pairs of individual yolk sacs and their cognate embryos, using the non-obese diabetic (NOD) mouse model. The analysis was performed at gestational day 8.5 on morphologically normal specimen to circumvent confounding by defective development. Results: Even with large sample numbers (n = 33 in each group), we observed considerable variability of gene expression, primarily driven by exposure to maternal diabetes, and secondarily by developmental stage of the embryo. Only a moderate number of genes changed expression in the yolk sac, while in the embryo, the exposure distinctly influenced the relationship of gene expression levels to developmental progression, revealing a possible role for altered cell cycle regulation in the response. Also affected in embryos under diabetic conditions were genes involved in cholesterol biosynthesis and NAD metabolism pathways. Discussion: Exposure to maternal diabetes during gastrulation changes transcriptomic profiles in embryos to a substantially greater effect than in the corresponding yolk sacs, indicating that despite yolk sac being of embryonic origin, different mechanisms control transcriptional activity in these tissues. The effects of maternal diabetes on expression of many genes that are correlated with developmental progression (i.e. somite stage) highlight the importance of considering developmental maturity in the interpretation of transcriptomic data. Our analyses identified cholesterol biosynthesis and NAD metabolism as novel pathways not previously implicated in diabetic pregnancies. Both NAD and cholesterol availability affect a wide variety of cellular signaling processes, and can be modulated by diet, implying that prevention of adverse outcomes from diabetic pregnancies may require broad interventions, particularly in the early stages of pregnancy.

Transposable elements as tissue-specific enhancers in cancers of endodermal lineage

Transposable elements (TE) are repetitive genomic elements that harbor binding sites for human transcription factors (TF). A regulatory role for TEs has been suggested in embryonal development and diseases such as cancer but systematic investigation of their functions has been limited by their widespread silencing in the genome. Here, we utilize unbiased massively parallel reporter assay data using a whole human genome library to identify TEs with functional enhancer activity in two human cancer types of endodermal lineage, colorectal and liver cancers. We show that the identified TE enhancers are characterized by genomic features associated with active enhancers, such as epigenetic marks and TF binding. Importantly, we identify distinct TE subfamilies that function as tissue-specific enhancers, namely MER11- and LTR12-elements in colon and liver cancers, respectively. These elements are bound by distinct TFs in each cell type, and they have predicted associations to differentially expressed genes. In conclusion, these data demonstrate how different cancer types can utilize distinct TEs as tissue-specific enhancers, paving the way for comprehensive understanding of the role of TEs as bona fide enhancers in the cancer genomes.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Predicting the impact of sequence motifs on gene regulation using single-cell data

The binding of transcription factors at proximal promoters and distal enhancers is central to gene regulation. Identifying regulatory motifs and quantifying their impact on expression remains challenging. Using a convolutional neural network trained on single-cell data, we infer putative regulatory motifs and cell type-specific importance. Our model, scover, explains 29% of the variance in gene expression in multiple mouse tissues. Applying scover to distal enhancers identified using scATAC-seq from the developing human brain, we identify cell type-specific motif activities in distal enhancers. Scover can identify regulatory motifs and their importance from single-cell data where all parameters and outputs are easily interpretable.

Determining DNA-Protein Binding Affinities and Specificities from Crude Lysates Using a Combined SILAC/TMT Labeling Strategy

In recent years, quantitative mass spectrometry-based interaction proteomics technology has proven very useful in identifying specific DNA-protein interactions using single pull-downs from crude lysates. Here, we applied a SILAC/TMT-based higher-order multiplexing approach to develop an interaction proteomics workflow called Protein-nucleic acid Affinity and Specificity quantification by MAss spectrometry in Nuclear extracts or PASMAN. In PASMAN, DNA pull-downs using a concentration range of specific and control DNA baits are performed in SILAC-labeled nuclear extracts. MS1-based quantification to determine specific DNA-protein interactions is then combined with sequential TMT-based quantification of fragmented SILAC peptides, allowing the generation of Hill-like curves and determination of apparent binding affinities. We benchmarked PASMAN using the SP/KLF motif and further applied it to gain insights into two CGCG-containing consensus DNA motifs. These motifs are recognized by two BEN domain-containing proteins, BANP and BEND3, which we find to interact with these motifs with distinct affinities. Finally, we profiled the BEND3 proximal proteome, revealing the NuRD complex as the major BEND3 proximal protein complex in vivo. In summary, PASMAN represents, to our knowledge, the first higher-order multiplexing-based interaction proteomics method that can be used to decipher specific DNA-protein interactions and their apparent affinities in various biological and pathological contexts.

A genome-wide screen reveals new regulators of the 2-cell-like cell state

In mammals, only the zygote and blastomeres of the early embryo are totipotent. This totipotency is mirrored in vitro by mouse '2-cell-like cells' (2CLCs), which appear at low frequency in cultures of embryonic stem cells (ESCs). Because totipotency is not completely understood, we carried out a genome-wide CRISPR knockout screen in mouse ESCs, searching for mutants that reactivate the expression of Dazl, a gene expressed in 2CLCs. Here we report the identification of four mutants that reactivate Dazl and a broader 2-cell-like signature: the E3 ubiquitin ligase adaptor SPOP, the Zinc-Finger transcription factor ZBTB14, MCM3AP, a component of the RNA processing complex TREX-2, and the lysine demethylase KDM5C. All four factors function upstream of DPPA2 and DUX, but not via p53. In addition, SPOP binds DPPA2, and KDM5C interacts with ncPRC1.6 and inhibits 2CLC gene expression in a catalytic-independent manner. These results extend our knowledge of totipotency, a key phase of organismal life.

Remodeling oncogenic transcriptomes by small molecules targeting NONO

Much of the human proteome is involved in mRNA homeostasis, but most RNA-binding proteins lack chemical probes. Here we identify electrophilic small molecules that rapidly and stereoselectively decrease the expression of transcripts encoding the androgen receptor and its splice variants in prostate cancer cells. We show by chemical proteomics that the compounds engage C145 of the RNA-binding protein NONO. Broader profiling revealed that covalent NONO ligands suppress an array of cancer-relevant genes and impair cancer cell proliferation. Surprisingly, these effects were not observed in cells genetically disrupted for NONO, which were instead resistant to NONO ligands. Reintroduction of wild-type NONO, but not a C145S mutant, restored ligand sensitivity in NONO-disrupted cells. The ligands promoted NONO accumulation in nuclear foci and stabilized NONO-RNA interactions, supporting a trapping mechanism that may prevent compensatory action of paralog proteins PSPC1 and SFPQ. These findings show that NONO can be co-opted by covalent small molecules to suppress protumorigenic transcriptional networks.

Readout of histone methylation by Trim24 locally restricts chromatin opening by p53

The genomic binding sites of the transcription factor (TF) and tumor suppressor p53 are unusually diverse with regard to their chromatin features, including histone modifications, raising the possibility that the local chromatin environment can contextualize p53 regulation. Here, we show that epigenetic characteristics of closed chromatin, such as DNA methylation, do not influence the binding of p53 across the genome. Instead, the ability of p53 to open chromatin and activate its target genes is locally restricted by its cofactor Trim24. Trim24 binds to both p53 and unmethylated histone 3 lysine 4 (H3K4), thereby preferentially localizing to those p53 sites that reside in closed chromatin, whereas it is deterred from accessible chromatin by H3K4 methylation. The presence of Trim24 increases cell viability upon stress and enables p53 to affect gene expression as a function of the local chromatin state. These findings link H3K4 methylation to p53 function and illustrate how specificity in chromatin can be achieved, not by TF-intrinsic sensitivity to histone modifications, but by employing chromatin-sensitive cofactors that locally modulate TF function.

grandR: a comprehensive package for nucleotide conversion RNA-seq data analysis

Metabolic labeling of RNA is a powerful technique for studying the temporal dynamics of gene expression. Nucleotide conversion approaches greatly facilitate the generation of data but introduce challenges for their analysis. Here we present grandR, a comprehensive package for quality control, differential gene expression analysis, kinetic modeling, and visualization of such data. We compare several existing methods for inference of RNA synthesis rates and half-lives using progressive labeling time courses. We demonstrate the need for recalibration of effective labeling times and introduce a Bayesian approach to study the temporal dynamics of RNA using snapshot experiments.

Unanticipated broad phylogeny of BEN DNA-binding domains revealed by structural homology searches

Organization of protein sequences into domain families is a foundation for cataloging and investigating protein functions. However, long-standing strategies based on primary amino acid sequences are blind to the possibility that proteins with dissimilar sequences could have comparable tertiary structures. Building on our recent findings that in silico structural predictions of BEN family DNA-binding domains closely resemble their experimentally determined crystal structures, we exploited the AlphaFold2 database for comprehensive identification of BEN domains. Indeed, we identified numerous novel BEN domains, including members of new subfamilies. For example, while no BEN domain factors had previously been annotated in C. elegans, this species actually encodes multiple BEN proteins. These include key developmental timing genes of orphan domain status, sel-7 and lin-14, the latter being the central target of the founding miRNA lin-4. We also reveal that the domain of unknown function 4806 (DUF4806), which is widely distributed across metazoans, is structurally similar to BEN and comprises a new subtype. Surprisingly, we find that BEN domains resemble both metazoan and non-metazoan homeodomains in 3D conformation and preserve characteristic residues, indicating that despite their inability to be aligned by conventional methods, these DNA-binding modules are probably evolutionarily related. Finally, we broaden the application of structural homology searches by revealing novel human members of DUF3504, which exists on diverse proteins with presumed or known nuclear functions. Overall, our work strongly expands this recently identified family of transcription factors and illustrates the value of 3D structural predictions to annotate protein domains and interpret their functions.

Targeting tumor endothelial cells with methyltransferase inhibitors: Mechanisms of action and the potential of combination therapy

Tumor endothelial cells (TECs) reside in the inner lining of blood vessels and represent a promising target for targeted cancer therapy. DNA methylation is a chemical process that involves the transfer of a methyl group to a specific base in the DNA strand, catalyzed by DNA methyltransferase (DNMT). DNMT inhibitors (DNMTis) can inhibit the activity of DNMTs, thereby preventing the transfer of methyl groups from s-adenosyl methionine (SAM) to cytosine. Currently, the most viable therapy for TECs is the development of DNMTis to release cancer suppressor genes from their repressed state. In this review, we first outline the characteristics of TECs and describe the development of tumor blood vessels and TECs. Abnormal DNA methylation is closely linked to tumor initiation, progression, and cell carcinogenesis, as evidenced by numerous studies. Therefore, we summarize the role of DNA methylation and DNA methyltransferase and the therapeutic potential of four types of DNMTi in targeting TECs. Finally, we discuss the accomplishments, challenges, and opportunities associated with combination therapy with DNMTis for TECs.

Exploring the Genetic Predisposition to Epigenetic Changes in Alzheimer's Disease

Alzheimer's disease (AD) is a prevalent type of dementia in elderly populations with a significant genetic component. The accumulating evidence suggests that AD involves a reconfiguration of the epigenetic landscape, including DNA methylation, post-translational modification of histone proteins, and chromatin remodeling. Along with environmental factors, individual specific genetic features play a considerable role in the formation of epigenetic architecture. In this study, we attempt to identify the non-coding regulatory SNPs (rSNPs) able to affect the epigenetic mechanisms in AD. To this end, the multi-omics approach is used. The GEO (Gene Expression Omnibus) available data (GSE153875) for AD patients and controls are integrated to reveal the rSNPs that display allele-specific features in both ChIP-seq profiles of four histone modifications and RNA-seq. Furthermore, we analyze the presence of rSNPs in the promoters of genes reported to be differentially expressed between AD and the normal brain (AD-related genes) and involved in epigenetic regulation according to the EpiFactors database. We also searched for the rSNPs in the promoters of the genes coding for transcription regulators of the identified AD-related genes. These regulators were selected based on the corresponding ChIP-seq peaks (ENCODE) in the promoter regions of these genes. Finally, we formed a panel of rSNPs localized to the promoters of genes that contribute to the epigenetic landscape in AD and, thus, to the genetic predisposition for this disease.

Structural insights into DNA recognition by the BEN domain of the transcription factor BANP

The BEN domain-containing transcription factors regulate transcription by recruiting chromatin-modifying factors to specific chromatin regions via their DNA-binding BEN domains. The BEN domain of BANP has been shown to bind to a CGCG DNA sequence or an AAA-containing MARs (matrix attachment regions) DNA sequence. Consistent with these in vivo observations, we identified an optimal DNA binding sequence of AAATCTCG by PBM (protein binding microarray), which was also confirmed by our ITC (Isothermal Titration Calorimetry) and mutagenesis results to uncover additional mechanistic details about DNA binding by the BEN domain of BANP. We then determined crystal structures of the BANP BEN domain in apo form and in complex with a CGCG-containing DNA, respectively, which revealed that the BANP BEN domain mainly used the electrostatic interactions to bind DNA with some base-specific interactions with the TC motifs. Our ITC results also showed that BANP bound to unmethylated and methylated DNAs with comparable binding affinities. Our complex structure of BANP-mCGCG revealed that the BANP BEN domain bound to the unmethylated and methylated DNAs in a similar mode and cytosine methylation did not get involved in binding, which is also consistent with our observations from the complex structures of the BEND6 BEN domain with the CGCG or CGmCG DNAs. Taken together, our results further elucidate the elements important for DNA recognition and transcriptional regulation by the BANP BEN domain-containing transcription factor.

Single-molecule footprinting identifies context-dependent regulation of enhancers by DNA methylation

Enhancers are cis-regulatory elements that control the establishment of cell identities during development. In mammals, enhancer activation is tightly coupled with DNA demethylation. However, whether this epigenetic remodeling is necessary for enhancer activation is unknown. Here, we adapted single-molecule footprinting to measure chromatin accessibility and transcription factor binding as a function of the presence of methylation on the same DNA molecules. We leveraged natural epigenetic heterogeneity at active enhancers to test the impact of DNA methylation on their chromatin accessibility in multiple cell lineages. Although reduction of DNA methylation appears dispensable for the activity of most enhancers, we identify a class of cell-type-specific enhancers where DNA methylation antagonizes the binding of transcription factors. Genetic perturbations reveal that chromatin accessibility and transcription factor binding require active demethylation at these loci. Thus, in addition to safeguarding the genome from spurious activation, DNA methylation directly controls transcription factor occupancy at active enhancers.

Timed exercise stabilizes behavioral rhythms but not molecular programs in the brain's suprachiasmatic clock

Timed daily access to a running-wheel (scheduled voluntary exercise; SVE) synchronizes rodent circadian rhythms and promotes stable, 24h rhythms in animals with genetically targeted impairment of neuropeptide signaling (Vipr2 ^-/- mice). Here we used RNA-seq and/or qRT-PCR to assess how this neuropeptide signaling impairment as well as SVE shapes molecular programs in the brain clock (suprachiasmatic nuclei; SCN) and peripheral tissues (liver and lung). Compared to Vipr2 ^+/+ animals, the SCN transcriptome of Vipr2 ^-/- mice showed extensive dysregulation which included core clock components, transcription factors, and neurochemicals. Furthermore, although SVE stabilized behavioral rhythms in these animals, the SCN transcriptome remained dysregulated. The molecular programs in the lung and liver of Vipr2 ^-/- mice were partially intact, although their response to SVE differed to that of these peripheral tissues in the Vipr2 ^+/+ mice. These findings highlight that SVE can correct behavioral abnormalities in circadian rhythms without causing large scale alterations to the SCN transcriptome.

The homeodomain of Oct4 is a dimeric binder of methylated CpG elements

Oct4 is essential to maintain pluripotency and has a pivotal role in establishing the germline. Its DNA-binding POU domain was recently found to bind motifs with methylated CpG elements normally associated with epigenetic silencing. However, the mode of binding and the consequences of this capability has remained unclear. Here, we show that Oct4 binds to a compact palindromic DNA element with a methylated CpG core (CpGpal) in alternative states of pluripotency and during cellular reprogramming towards induced pluripotent stem cells (iPSCs). During cellular reprogramming, typical Oct4 bound enhancers are uniformly demethylated, with the prominent exception of the CpGpal sites where DNA methylation is often maintained. We demonstrate that Oct4 cooperatively binds the CpGpal element as a homodimer, which contrasts with the ectoderm-expressed POU factor Brn2. Indeed, binding to CpGpal is Oct4-specific as other POU factors expressed in somatic cells avoid this element. Binding assays combined with structural analyses and molecular dynamic simulations show that dimeric Oct4-binding to CpGpal is driven by the POU-homeodomain whilst the POU-specific domain is detached from DNA. Collectively, we report that Oct4 exerts parts of its regulatory function in the context of methylated DNA through a DNA recognition mechanism that solely relies on its homeodomain.

Functional impact of multi-omic interactions in breast cancer subtypes

Multi-omic approaches are expected to deliver a broader molecular view of cancer. However, the promised mechanistic explanations have not quite settled yet. Here, we propose a theoretical and computational analysis framework to semi-automatically produce network models of the regulatory constraints influencing a biological function. This way, we identified functions significantly enriched on the analyzed omics and described associated features, for each of the four breast cancer molecular subtypes. For instance, we identified functions sustaining over-representation of invasion-related processes in the basal subtype and DNA modification processes in the normal tissue. We found limited overlap on the omics-associated functions between subtypes; however, a startling feature intersection within subtype functions also emerged. The examples presented highlight new, potentially regulatory features, with sound biological reasons to expect a connection with the functions. Multi-omic regulatory networks thus constitute reliable models of the way omics are connected, demonstrating a capability for systematic generation of mechanistic hypothesis.

Genome-wide association study for the primary feather color trait in a native Chinese duck

Background: To reveal candidate genes and the molecular genetic mechanism underlying primary feather color trait in ducks, a genome-wide association study (GWAS) for the primary feather color trait was performed based on the genotyping-by-sequencing (GBS) technology for a native Chinese female duck, Longyan Shan-ma ducks. Methods: Blood genomic DNA from 314 female Longyan Shan-ma duck were genotyped using GBS technology. A GWAS for the primary feather color trait with genome variations was performed using an univariate linear mixed model based on all SNPs in autosomes. Results: Seven genome-wide significant single nucleotide polymorphisms (SNPs, Bonferroni-adjusted p-value <8.03 × 10^-7) within the introns of the genes STARD9, ZNF106, SLC7A5, and BANP genes were associated with the primary feather color trait. Twenty-two genome-wide suggestive SNPs (Bonferroni-adjusted p-value <1.61 × 10^-5) of 17 genes (besides ZNF106 and SLC7A5) were also identified. Seven SNPs were located at one 0.22 Mb region (38.65-38.87 Mb) on chromosome 5, and six SNPs were located at one 0.31 Mb region (19.53-19.84 Mb) on chromosome 11. The functions of STARD9, SLC7A5, BANP, LOC101798015, and IPMK were involved pigmentation and follicle development, especially, STARD9 upregulated expression in black feather (haplotype-CCCC) bulb tissue compared with in pockmarked feather (haplotype-TGTT) bulb tissue, implicating these genes as candidate genes for primary feather color trait. Conclusion: The preliminarily findings suggested candidate genes and regions, and the genetic basis of primary feather color trait in a female duck.

Identification of copy number variations in the genome of Dairy Gir cattle

Studying structural variants that can control complex traits is relevant for dairy cattle production, especially for animals that are tolerant to breeding conditions in the tropics, such as the Dairy Gir cattle. This study identified and characterized high confidence copy number variation regions (CNVR) in the Gir breed genome. A total of 38 animals were whole-genome sequenced, and 566 individuals were genotyped with a high-density SNP panel, among which 36 animals had both sequencing and SNP genotyping data available. Two sets of high confidence CNVR were established: one based on common CNV identified in the studied population (CNVR_POP), and another with CNV identified in sires with both sequence and SNP genotyping data available (CNVR_ANI). We found 10 CNVR_POP and 45 CNVR_ANI, which covered 1.05 Mb and 4.4 Mb of the bovine genome, respectively. Merging these CNV sets for functional analysis resulted in 48 unique high confidence CNVR. The overlapping genes were previously related to embryonic mortality, environmental adaptation, evolutionary process, immune response, longevity, mammary gland, resistance to gastrointestinal parasites, and stimuli recognition, among others. Our results contribute to a better understanding of the Gir breed genome. Moreover, the CNV identified in this study can potentially affect genes related to complex traits, such as production, health, and reproduction.

Evidence that direct inhibition of transcription factor binding is the prevailing mode of gene and repeat repression by DNA methylation

Cytosine methylation efficiently silences CpG-rich regulatory regions of genes and repeats in mammalian genomes. To what extent this entails direct inhibition of transcription factor (TF) binding versus indirect inhibition via recruitment of methyl-CpG-binding domain (MBD) proteins is unclear. Here we show that combinatorial genetic deletions of all four proteins with functional MBDs in mouse embryonic stem cells, derived neurons or a human cell line do not reactivate genes or repeats with methylated promoters. These do, however, become activated by methylation-restricted TFs if DNA methylation is removed. We identify several causal TFs in neurons, including ONECUT1, which is methylation sensitive only at a motif variant. Rampantly upregulated retrotransposons in methylation-free neurons feature a CRE motif, which activates them in the absence of DNA methylation via methylation-sensitive binding of CREB1. Our study reveals methylation-sensitive TFs in vivo and argues that direct inhibition, rather than indirect repression by the tested MBD proteins, is the prevailing mechanism of methylation-mediated repression at regulatory regions and repeats.

The transcriptomic landscape of elderly acute myeloid leukemia identifies B7H3 and BANP as a favorable signature in high-risk patients

Acute myeloid leukemia (AML) in the elderly remains a clinical challenge, with a five-year overall survival rate below 10%. The current ELN 2017 genetic risk classification considers cytogenetic and mutational characteristics to stratify fit AML patients into different prognostic groups. However, this classification is not validated for elderly patients treated with a non-intensive approach, and its performance may be suboptimal in this context. Indeed, the transcriptomic landscape of AML in the elderly has been less explored and it might help stratify this group of patients. In the current study, we analyzed the transcriptome of 224 AML patients > 65 years-old at diagnosis treated in the Spanish PETHEMA-FLUGAZA clinical trial in order to identify new prognostic biomarkers in this population. We identified a specific transcriptomic signature for high-risk patients with mutated TP53 or complex karyotype, revealing that low expression of B7H3 gene with high expression of BANP gene identifies a subset of high-risk AML patients surviving more than 12 months. This result was further validated in the BEAT AML cohort. This unique signature highlights the potential of transcriptomics to identify prognostic biomarkers in in elderly AML.

Combining NGN2 programming and dopaminergic patterning for a rapid and efficient generation of hiPSC-derived midbrain neurons

The use of human derived induced pluripotent stem cells (hiPSCs) differentiated to dopaminergic (DA) neurons offers a valuable experimental model to decorticate the cellular and molecular mechanisms of Parkinson's disease (PD) pathogenesis. However, the existing approaches present with several limitations, notably the lengthy time course of the protocols and the high variability in the yield of DA neurons. Here we report on the development of an improved approach that combines neurogenin-2 programming with the use of commercially available midbrain differentiation kits for a rapid, efficient, and reproducible directed differentiation of hiPSCs to mature and functional induced DA (iDA) neurons, with minimum contamination by other brain cell types. Gene expression analysis, associated with functional characterization examining neurotransmitter release and electrical recordings, support the functional identity of the iDA neurons to A9 midbrain neurons. iDA neurons showed selective vulnerability when exposed to 6-hydroxydopamine, thus providing a viable in vitro approach for modeling PD and for the screening of small molecules with neuroprotective proprieties.

Banp regulates DNA damage response and chromosome segregation during the cell cycle in zebrafish retina

Btg3-associated nuclear protein (Banp) was originally identified as a nuclear matrix-associated region (MAR)-binding protein and it functions as a tumor suppressor. At the molecular level, Banp regulates transcription of metabolic genes via a CGCG-containing motif called the Banp motif. However, its physiological roles in embryonic development are unknown. Here, we report that Banp is indispensable for the DNA damage response and chromosome segregation during mitosis. Zebrafish banp mutants show mitotic cell accumulation and apoptosis in developing retina. We found that DNA replication stress and tp53-dependent DNA damage responses were activated to induce apoptosis in banp mutants, suggesting that Banp is required for regulation of DNA replication and DNA damage repair. Furthermore, consistent with mitotic cell accumulation, chromosome segregation was not smoothly processed from prometaphase to anaphase in banp morphants, leading to a prolonged M-phase. Our RNA- and ATAC-sequencing identified 31 candidates for direct Banp target genes that carry the Banp motif. Interestingly, a DNA replication fork regulator, wrnip1, and two chromosome segregation regulators, cenpt and ncapg, are included in this list. Thus, Banp directly regulates transcription of wrnip1 for recovery from DNA replication stress, and cenpt and ncapg for chromosome segregation during mitosis. Our findings provide the first in vivo evidence that Banp is required for cell-cycle progression and cell survival by regulating DNA damage responses and chromosome segregation during mitosis.

In order for a cell to divide, it must progress through a series of carefully controlled steps known as the cell cycle. First, the cell replicates its DNA and both copies get segregated to opposite ends. The cell then splits into two and each new cell receives a copy of the duplicated genetic material. If any of the stages in the cell cycle become disrupted or mis-regulated this can lead to uncontrolled divisions that may result in cancer.

Researchers have often used a structure within the eye known as the retina to study the cell cycle in zebrafish and other animals as cells in the retina rapidly divide in a highly controlled manner. A protein called Banp is known to help stop tumors from growing in humans and mice, but its normal role in the body, particularly the cell cycle, has remained unclear. To investigate, Babu et al. studied the retina of mutant zebrafish that were unable to make the Banp protein.

The experiments revealed that two stress responses indicating DNA damage or defects in copying DNA were active in the retinal cells of the mutant zebrafish. This suggested that Banp allows cell to progress through the cell cycle by repairing any DNA damage that may arise during replication. Banp does this by activating the gene for another protein called Wrnip1. Babu et al. also found that Banp helps segregate the two copies of DNA during cell division by promoting the activation of two other proteins called Cenpt and Ncapg. Further experiments identified 31 genes that were directly regulated by Banp.

These findings demonstrate that Banp is required for zebrafish cells to be able to accurately copy their DNA and divide in to two new cells. In the future, the work of Babu et al. will provide a useful resource to investigate how tumors grow and spread around the body, and may contribute to the development of new treatments for cancer.

Generating specificity in genome regulation through transcription factor sensitivity to chromatin

Cell type-specific gene expression relies on transcription factors (TFs) binding DNA sequence motifs embedded in chromatin. Understanding how motifs are accessed in chromatin is crucial to comprehend differential transcriptional responses and the phenotypic impact of sequence variation. Chromatin obstacles to TF binding range from DNA methylation to restriction of DNA access by nucleosomes depending on their position, composition and modification. In vivo and in vitro approaches now enable the study of TF binding in chromatin at unprecedented resolution. Emerging insights suggest that TFs vary in their ability to navigate chromatin states. However, it remains challenging to link binding and transcriptional outcomes to molecular characteristics of TFs or the local chromatin substrate. Here, we discuss our current understanding of how TFs access DNA in chromatin and novel techniques and directions towards a better understanding of this critical step in genome regulation.

Highly enriched BEND3 prevents the premature activation of bivalent genes during differentiation

Bivalent genes are ready for activation upon the arrival of developmental cues. Here, we report that BEND3 is a CpG island (CGI)-binding protein that is enriched at regulatory elements. The cocrystal structure of BEND3 in complex with its target DNA reveals the structural basis for its DNA methylation-sensitive binding property. Mouse embryos ablated of Bend3 died at the pregastrulation stage. Bend3 null embryonic stem cells (ESCs) exhibited severe defects in differentiation, during which hundreds of CGI-containing bivalent genes were prematurely activated. BEND3 is required for the stable association of polycomb repressive complex 2 (PRC2) at bivalent genes that are highly occupied by BEND3, which suggests a reining function of BEND3 in maintaining high levels of H3K27me3 at these bivalent genes in ESCs to prevent their premature activation in the forthcoming developmental stage.

Histochemistry, Cytochemistry and Epigenetics

Over the past few decades, many researchers have individually identified tumor-related genes, and have accumulated information on their basic research in a database. With the development of technology that can comprehensively test the expression status within a short time, oncogene panel testing has become attainable. On the other hand, changes in gene expression that do not depend on changes in base sequences, that is, epigenetics, or more comprehensively, epigenomes, are also highly involved in the development and progression of disease. Oncogene panel tests tend to focus on DNA base mutations such as point mutations, deletions, duplications, and chimera formation. Elucidation leads to correct interpretation of diseases and treatment choices, and we are in an era where integrated understanding of the genome and epigenome is indispensable. In this review, we make every effort to cover a wide range of knowledge, including data on histone protein modification, non-coding (nc)RNA and DNA methylation, and recent application trials for demonstrating epigenetic alterations in histologic and cytologic specimens. We hope this review will help marshal the knowledge accumulated by researchers involved in genomic and epigenomic studies.

Sequence determinants of human gene regulatory elements

DNA can determine where and when genes are expressed, but the full set of sequence determinants that control gene expression is unknown. Here, we measured the transcriptional activity of DNA sequences that represent an ~100 times larger sequence space than the human genome using massively parallel reporter assays (MPRAs). Machine learning models revealed that transcription factors (TFs) generally act in an additive manner with weak grammar and that most enhancers increase expression from a promoter by a mechanism that does not appear to involve specific TF–TF interactions. The enhancers themselves can be classified into three types: classical, closed chromatin and chromatin dependent. We also show that few TFs are strongly active in a cell, with most activities being similar between cell types. Individual TFs can have multiple gene regulatory activities, including chromatin opening and enhancing, promoting and determining transcription start site (TSS) activity, consistent with the view that the TF binding motif is the key atomic unit of gene expression.

Analysis of massively parallel reporter assays measuring the transcriptional activity of DNA sequences indicates that most transcription factor (TF) activity is additive and does not rely on specific TF–TF interactions. Individual TFs can have different gene regulatory activities.

Distinct structural bases for sequence-specific DNA binding by mammalian BEN domain proteins

The BEN domain is a recently recognized DNA binding module that is present in diverse metazoans and certain viruses. Several BEN domain factors are known as transcriptional repressors, but, overall, relatively little is known of how BEN factors identify their targets in humans. In particular, X-ray structures of BEN domain:DNA complexes are only known for Drosophila factors bearing a single BEN domain, which lack direct vertebrate orthologs. Here, we characterize several mammalian BEN domain (BD) factors, including from two NACC family BTB-BEN proteins and from BEND3, which has four BDs. In vitro selection data revealed sequence-specific binding activities of isolated BEN domains from all of these factors. We conducted detailed functional, genomic, and structural studies of BEND3. We show that BD4 is a major determinant for in vivo association and repression of endogenous BEND3 targets. We obtained a high-resolution structure of BEND3-BD4 bound to its preferred binding site, which reveals how BEND3 identifies cognate DNA targets and shows differences with one of its non-DNA-binding BEN domains (BD1). Finally, comparison with our previous invertebrate BEN structures, along with additional structural predictions using AlphaFold2 and RoseTTAFold, reveal distinct strategies for target DNA recognition by different types of BEN domain proteins. Together, these studies expand the DNA recognition activities of BEN factors and provide structural insights into sequence-specific DNA binding by mammalian BEN proteins.