Comparative FAIRE-seq analysis reveals distinguishing features of the chromatin structure of ground state- and primed-pluripotent cells

Prediction accuracy of regulatory elements from sequence varies by functional sequencing technique

Introduction

Various sequencing based approaches are used to identify and characterize the activities of cis-regulatory elements in a genome-wide fashion. Some of these techniques rely on indirect markers such as histone modifications (ChIP-seq with histone antibodies) or chromatin accessibility (ATAC-seq, DNase-seq, FAIRE-seq), while other techniques use direct measures such as episomal assays measuring the enhancer properties of DNA sequences (STARR-seq) and direct measurement of the binding of transcription factors (ChIP-seq with transcription factor-specific antibodies). The activities of cis-regulatory elements such as enhancers, promoters, and repressors are determined by their sequence and secondary processes such as chromatin accessibility, DNA methylation, and bound histone markers.

Methods

Here, machine learning models are employed to evaluate the accuracy with which cis-regulatory elements identified by various commonly used sequencing techniques can be predicted by their underlying sequence alone to distinguish between cis-regulatory activity that is reflective of sequence content versus secondary processes.

Results and discussion

Models trained and evaluated on D. melanogaster sequences identified through DNase-seq and STARR-seq are significantly more accurate than models trained on sequences identified by H3K4me1, H3K4me3, and H3K27ac ChIP-seq, FAIRE-seq, and ATAC-seq. These results suggest that the activity detected by DNase-seq and STARR-seq can be largely explained by underlying DNA sequence, independent of secondary processes. Experimentally, a subset of DNase-seq and H3K4me1 ChIP-seq sequences were tested for enhancer activity using luciferase assays and compared with previous tests performed on STARR-seq sequences. The experimental data indicated that STARR-seq sequences are substantially enriched for enhancer-specific activity, while the DNase-seq and H3K4me1 ChIP-seq sequences are not. Taken together, these results indicate that the DNase-seq approach identifies a broad class of regulatory elements of which enhancers are a subset and the associated data are appropriate for training models for detecting regulatory activity from sequence alone, STARR-seq data are best for training enhancer-specific sequence models, and H3K4me1 ChIP-seq data are not well suited for training and evaluating sequence-based models for cis-regulatory element prediction.

FAIRE-MS reveals mitotic retention of transcriptional regulators on a proteome-wide scale

Mitosis entails global and dramatic alterations, such as higher-order chromatin organization disruption, concomitant with global transcription downregulation. Cells reliably re-establishing gene expression patterns upon mitotic exit and maintaining cellular identities remain poorly understood. Previous studies indicated that certain transcription factors (TFs) remain associated with individual loci during mitosis and serve as mitotic bookmarkers. However, it is unclear which regulatory factors remain bound to the compacted mitotic chromosomes. We developed formaldehyde-assisted isolation of regulatory elements-coupled mass spectrometry (FAIRE-MS) that combines FAIRE-based open chromatin-associated protein pull-down and mass spectrometry (MS) to quantify the open chromatin-associated proteome during the interphase and mitosis. We identified 189 interphase and mitosis maintained (IM) regulatory factors using FAIRE-MS and found intrinsically disordered proteins and regions (IDP(R)s) are highly enriched, which plays a crucial role in liquid-liquid phase separation (LLPS) and chromatin organization during the cell cycle. Notably, in these IDP(R)s, we identified mitotic bookmarkers, such as CEBPB, HMGB1, and TFAP2A, and several factors, including MAX, HMGB3, hnRNP A2/B1, FUS, hnRNP D, and TIAL1, which are at least partially bound to the mitotic chromosome. Furthermore, it will be essential to study whether these IDP(R)s through LLPS helps cells transit from mitosis to the G1 phase during the cell cycle.

SWI/SNF (BAF) complexes: From framework to a functional role in endothelial mechanotransduction

Endothelial cells (ECs) are constantly subjected to an array of mechanical cues, especially shear stress, due to their luminal placement in the blood vessels. Blood flow can regulate various aspects of endothelial biology and pathophysiology by regulating the endothelial processes at the transcriptomic, proteomic, miRNomic, metabolomics, and epigenomic levels. ECs sense, respond, and adapt to altered blood flow patterns and shear profiles by specialized mechanisms of mechanosensing and mechanotransduction, resulting in qualitative and quantitative differences in their gene expression. Chromatin-regulatory proteins can regulate transcriptional activation by modifying the organization of nucleosomes at promoters, enhancers, silencers, insulators, and locus control regions. Recent research efforts have illustrated that SWI/SNF (SWItch/Sucrose Non-Fermentable) or BRG1/BRM-associated factor (BAF) complex regulates DNA accessibility and chromatin structure. Since the discovery, the gene-regulatory mechanisms of the BAF complex associated with chromatin remodeling have been intensively studied to investigate its role in diverse disease phenotypes. Thus far, it is evident that (1) the SWI/SNF complex broadly regulates the activity of transcriptional enhancers to control lineage-specific differentiation and (2) mutations in the BAF complex proteins lead to developmental disorders and cancers. It is unclear if blood flow can modulate the activity of SWI/SNF complex to regulate EC differentiation and reprogramming. This review emphasizes the integrative role of SWI/SNF complex from a structural and functional standpoint with a special reference to cardiovascular diseases (CVDs). The review also highlights how regulation of this complex by blood flow can lead to the discovery of new therapeutic interventions for the treatment of endothelial dysfunction in vascular diseases.

The Dynamics of Transcriptional Activation by Hepatic Reprogramming Factors

Specific combinations of two transcription factors (Hnf4α plus Foxa1, Foxa2, or Foxa3) can induce direct conversion of mouse fibroblasts into hepatocyte-like cells. However, the molecular mechanisms underlying hepatic reprogramming are largely unknown. Here, we show that the Foxa protein family members and Hnf4α sequentially and cooperatively bind to chromatin to activate liver-specific gene expression. Although all Foxa proteins bind to and open regions of closed chromatin as pioneer factors, Foxa3 has the unique potential of transferring from the distal to proximal regions of the transcription start site of target genes, binding RNA polymerase II, and co-traversing target genes. These distinctive characteristics of Foxa3 are essential for inducing the hepatic fate in fibroblasts. Similar functional coupling of transcription factors to RNA polymerase II may occur in other contexts whereby transcriptional activation can induce cell differentiation.

Genomic methods in profiling DNA accessibility and factor localization

Recent advancements in next-generation sequencing technologies and accompanying reductions in cost have led to an explosion of techniques to examine DNA accessibility and protein localization on chromatin genome-wide. Generally, accessible regions of chromatin are permissive for factor binding and are therefore hotspots for regulation of gene expression; conversely, genomic regions that are highly occupied by histone proteins are not permissive for factor binding and are less likely to be active regulatory regions. Identifying regions of differential accessibility can be useful to uncover putative gene regulatory regions, such as enhancers, promoters, and insulators. In addition, DNA-binding proteins, such as transcription factors that preferentially bind certain DNA sequences and histone proteins that form the core of the nucleosome, play essential roles in all DNA-templated processes. Determining the genomic localization of chromatin-bound proteins is therefore essential in determining functional roles, sequence motifs important for factor binding, and regulatory networks controlling gene expression. In this review, we discuss techniques for determining DNA accessibility and nucleosome positioning (DNase-seq, FAIRE-seq, MNase-seq, and ATAC-seq) and techniques for detecting and functionally characterizing chromatin-bound proteins (ChIP-seq, DamID, and CUT&RUN). These methods have been optimized to varying degrees of resolution, specificity, and ease of use. Here, we outline some advantages and disadvantages of these techniques, their general protocols, and a brief discussion of their development. Together, these complimentary approaches have provided an unparalleled view of chromatin architecture and functional gene regulation.

Detect accessible chromatin using ATAC-sequencing, from principle to applications

Background

Chromatin accessibility is crucial for gene expression regulation in specific cells and in multiple biological processes. Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) is an effective way to reveal chromatin accessibility at a genome-wide level. Through ATAC-seq, produced reads from a small number of cells reflect accessible regions that correspond to nucleosome positioning and transcription factor binding sites, due to probing hyperactive Tn5 transposase to DNA sequence.

Conclusion

In this review, we summarize both principle and features of ATAC-seq, highlight its applications in basic and clinical research. ATAC-seq has generated comprehensive chromatin accessible maps, and is becoming a powerful tool to understand dynamic gene expression regulation in stem cells, early embryos and tumors.

OCEAN-C: mapping hubs of open chromatin interactions across the genome reveals gene regulatory networks

We develop a method called open chromatin enrichment and network Hi-C (OCEAN-C) for antibody-independent mapping of global open chromatin interactions. By integrating FAIRE-seq and Hi-C, OCEAN-C detects open chromatin interactions enriched by active cis-regulatory elements. We identify more than 10,000 hubs of open chromatin interactions (HOCIs) in human cells, which are mainly active promoters and enhancers bound by many DNA-binding proteins and form interaction networks crucial for gene transcription. In addition to identifying large-scale topological structures, including topologically associated domains and A/B compartments, OCEAN-C can detect HOCI-mediated chromatin interactions that are strongly associated with gene expression, super-enhancers, and broad H3K4me3 domains.

Sirtuin 1 Promotes Deacetylation of Oct4 and Maintenance of Naive Pluripotency

The enhancer landscape is dramatically restructured as naive preimplantation epiblasts transition to the post-implantation state of primed pluripotency. A key factor in this process is Otx2, which is upregulated during the early stages of this transition and ultimately recruits Oct4 to a different set of enhancers. In this study, we discover that the acetylation status of Oct4 regulates the induction of the primed pluripotency gene network. Maintenance of the naive state requires the NAD-dependent deacetylase, SirT1, which deacetylates Oct4. The activity of SirT1 is reduced during the naive-to-primed transition; Oct4 becomes hyper-acetylated and binds to an Otx2 enhancer to induce Otx2 expression. Induction of Otx2 causes the reorganization of acetylated Oct4 and results in the induction of the primed pluripotency gene network. Regulation of Oct4 by SirT1 may link stem cell development to environmental conditions, and it may provide strategies to manipulate epiblast cell state.

Heterarchy of transcription factors driving basal and luminal cell phenotypes in human urothelium

Cell differentiation is affected by complex networks of transcription factors that co-ordinate re-organisation of the chromatin landscape. The hierarchies of these relationships can be difficult to dissect. During in vitro differentiation of normal human uro-epithelial cells, formaldehyde-assisted isolation of regulatory elements (FAIRE-seq) and RNA-seq was used to identify alterations in chromatin accessibility and gene expression changes following activation of the nuclear receptor peroxisome proliferator-activated receptor gamma (PPARγ) as a differentiation-initiating event. Regions of chromatin identified by FAIRE-seq, as having altered accessibility during differentiation, were found to be enriched with sequence-specific binding motifs for transcription factors predicted to be involved in driving basal and differentiated urothelial cell phenotypes, including forkhead box A1 (FOXA1), P63, GRHL2, CTCF and GATA-binding protein 3 (GATA3). In addition, co-occurrence of GATA3 motifs was observed within subsets of differentiation-specific peaks containing P63 or FOXA1. Changes in abundance of GRHL2, GATA3 and P63 were observed in immunoblots of chromatin-enriched extracts. Transient siRNA knockdown of P63 revealed that P63 favoured a basal-like phenotype by inhibiting differentiation and promoting expression of basal marker genes. GATA3 siRNA prevented differentiation-associated downregulation of P63 protein and transcript, and demonstrated positive feedback of GATA3 on PPARG transcript, but showed no effect on FOXA1 transcript or protein expression. This approach indicates that as a transcriptionally regulated programme, urothelial differentiation operates as a heterarchy, wherein GATA3 is able to co-operate with FOXA1 to drive expression of luminal marker genes, but that P63 has potential to transrepress expression of the same genes.

DNA methylation is dispensable for changes in global chromatin architecture but required for chromocentre formation in early stem cell differentiation

Epiblast stem cells (EpiSCs), which are pluripotent cells isolated from early post-implantation mouse embryos (E5.5), show both similarities and differences compared to mouse embryonic stem cells (mESCs), isolated earlier from the inner cell mass (ICM) of the E3.5 embryo. Previously, we have observed that while chromatin is very dispersed in E3.5 ICM, compact chromatin domains and chromocentres appear in E5.5 epiblasts after embryo implantation. Given that the observed chromatin re-organization in E5.5 epiblasts coincides with an increase in DNA methylation, in this study, we aimed to examine the role of DNA methylation in chromatin re-organization during the in vitro conversion of ESCs to EpiSCs. The requirement for DNA methylation was determined by converting both wild-type and DNA methylation-deficient ESCs to EpiSCs, followed by structural analysis with electron spectroscopic imaging (ESI). We show that the chromatin re-organization which occurs in vivo can be re-capitulated in vitro during the ESC to EpiSC conversion. Indeed, after 7 days in EpiSC media, compact chromatin domains begin to appear throughout the nuclear volume, creating a chromatin organization similar to E5 epiblasts and embryo-derived EpiSCs. Our data demonstrate that DNA methylation is dispensable for this global chromatin re-organization but required for the compaction of pericentromeric chromatin into chromocentres.

Widespread Chromatin Accessibility at Repetitive Elements Links Stem Cells with Human Cancer

Chromatin regulation is critical for differentiation and disease. However, features linking the chromatin environment of stem cells with disease remain largely unknown. We explored chromatin accessibility in embryonic and multipotent stem cells and unexpectedly identified widespread chromatin accessibility at repetitive elements. Integrating genomic and biochemical approaches, we demonstrate that these sites of increased accessibility are associated with well-positioned nucleosomes marked by distinct histone modifications. Differentiation is accompanied by chromatin remodeling at repetitive elements associated with altered expression of genes in relevant developmental pathways. Remarkably, we found that the chromatin environment of Ewing sarcoma, a mesenchymally derived tumor, is shared with primary mesenchymal stem cells (MSCs). Accessibility at repetitive elements in MSCs offers a permissive environment that is exploited by the critical oncogene responsible for this cancer. Our data demonstrate that stem cells harbor a unique chromatin landscape characterized by accessibility at repetitive elements, a feature associated with differentiation and oncogenesis.

Identification of chromatin accessibility domains in human breast cancer stem cells

Epithelial-to-mesenchymal transition (EMT) is physiological in embryogenesis and wound healing but also associated with the formation of cancer stem cells (CSCs). Many EMT signaling pathways are implicated in CSC formation, but the precise underlying mechanisms of CSC formation remain elusive. We have previously demonstrated that PKC is critical for EMT induction and CSC formation in inducible breast EMT/CSC models. Here, we used formaldehyde-assisted isolation of regulatory elements-sequencing (FAIRE-seq) to investigate DNA accessibility changes after PKC activation and determine how they influence EMT and CSC formation. During EMT, DNA accessibility principally increased in regions distant from transcription start sites, low in CpG content, and enriched with chromatin enhancer marks. ChIP-sequencing revealed that a subset of these regions changed from poised to active enhancers upon stimulation, with some even more acteylated in CSCs. While regions with increased accessibility were enriched for FOX, AP-1, TEAD, and TFAP2 motifs, those containing FOX and AP-1 motif were associated with increased expression of CSC-associated genes, while those with TFAP2 were associated with genes with increased expression in non-CSCs. Silencing of 2 members of the FOX family, FOXN2 and FOXQ1, repressed CSCs and the mesenchymal phenotype and inhibited the CSC gene signature. These novel, PKC-induced DNA accessibility regions help explain how the epigenomic plasticity of cells undergoing EMT leads to CSC gene activation.

Genome-wide targets identification of "core" pluripotency transcription factors with integrated features in human embryonic stem cells

Embryonic stem cells (ESCs) play an important role in developmental biology which is still lacking clear molecular mechanisms. The "core" transcription factors (TFs) including OCT4, SOX2 and NANOG are essential for maintaining the stemness of ESCs. But the downstream targets of these "core" TFs are still ambiguous. Based on support vector machine (SVM) technology, this study develops a label method algorithm (LMA) for genome-wide target identification of "core" TFs in humans, which eliminates the need for negative training samples. This method integrates histone modifications and TF binding motifs as identification features. Compared with a previous mapping-convergence (M-C) algorithm, the LMA can provide more stable and reliable predictions. 4796, 3166 and 4384 target genes of OCT4, SOX2 and NANOG, respectively, were identified with the LMA model. Then verifications of the predicted targets were carried out based on their functional consistency and their connection degree in networks from a computational system biology perspective. The results showed that the targets of "core" TFs present higher gene functional similarity and closer connection distance than background levels.

Chromatin features, RNA polymerase II and the comparative expression of lens genes encoding crystallins, transcription factors, and autophagy mediators

PURPOSE

Gene expression correlates with local chromatin structure. Our studies have mapped histone post-translational modifications, RNA polymerase II (pol II), and transcription factor Pax6 in lens chromatin. These data represent the first genome-wide insights into the relationship between lens chromatin structure and lens transcriptomes and serve as an excellent source for additional data analysis and refinement. The principal lens proteins, the crystallins, are encoded by predominantly expressed mRNAs; however, the regulatory mechanisms underlying their high expression in the lens remain poorly understood.

METHODS

The formaldehyde-assisted identification of regulatory regions (FAIRE-Seq) was employed to analyze newborn lens chromatin. ChIP-seq and RNA-seq data published earlier (GSE66961) have been used to assist in FAIRE-seq data interpretation. RNA transcriptomes from murine lens epithelium, lens fibers, erythrocytes, forebrain, liver, neurons, and pancreas were compared to establish the gene expression levels of the most abundant mRNAs versus median gene expression across other differentiated cells.

RESULTS

Normalized RNA expression data from multiple tissues show that crystallins rank among the most highly expressed genes in mammalian cells. These findings correlate with the extremely high abundance of pol II all across the crystallin loci, including crystallin genes clustered on chromosomes 1 and 5, as well as within regions of "open" chromatin, as identified by FAIRE-seq. The expression levels of mRNAs encoding DNA-binding transcription factors (e.g., Foxe3, Hsf4, Maf, Pax6, Prox1, Sox1, and Tfap2a) revealed that their transcripts form "clusters" of abundant mRNAs in either lens fibers or lens epithelium. The expression of three autophagy regulatory mRNAs, encoding Tfeb, FoxO1, and Hif1α, was found within a group of lens preferentially expressed transcription factors compared to the E12.5 forebrain.

CONCLUSIONS

This study reveals novel features of lens chromatin, including the remarkably high abundance of pol II at the crystallin loci that exhibit features of "open" chromatin. Hsf4 ranks among the most abundant fiber cell-preferred DNA-binding transcription factors. Notable transcripts, including Atf4, Ctcf, E2F4, Hey1, Hmgb1, Mycn, RXRβ, Smad4, Sp1, and Taf1 (transcription factors) and Ctsd, Gabarapl1, and Park7 (autophagy regulators) have been identified with high levels of expression in lens fibers, which suggests specific roles in lens fiber cell terminal differentiation.

High throughput technologies for the functional discovery of mammalian enhancers: new approaches for understanding transcriptional regulatory network dynamics

Completion of the human and mouse genomes has inspired new initiatives to obtain a global understanding of the functional regulatory networks governing gene expression. Enhancers are primary regulatory DNA elements determining precise spatio- and temporal gene expression patterns, but the observation that they can function at any distance from the gene(s) they regulate has made their genome-wide characterization challenging. Since traditional, single reporter approaches would be unable to accomplish this enormous task, high throughput technologies for mapping chromatin features associated with enhancers have emerged as an effective surrogate for enhancer discovery. However, the last few years have witnessed the development of several new innovative approaches that can effectively screen for and discover enhancers based on their functional activation of transcription using massively parallel reporter systems. In addition to their application for genome annotation, these new high throughput functional approaches open new and exciting avenues for modeling gene regulatory networks.