High-throughput localization of functional elements by quantitative chromatin profiling

Enhancer-MDLF: a novel deep learning framework for identifying cell-specific enhancers

Enhancers, noncoding DNA fragments, play a pivotal role in gene regulation, facilitating gene transcription. Identifying enhancers is crucial for understanding genomic regulatory mechanisms, pinpointing key elements and investigating networks governing gene expression and disease-related mechanisms. Existing enhancer identification methods exhibit limitations, prompting the development of our novel multi-input deep learning framework, termed Enhancer-MDLF. Experimental results illustrate that Enhancer-MDLF outperforms the previous method, Enhancer-IF, across eight distinct human cell lines and exhibits superior performance on generic enhancer datasets and enhancer-promoter datasets, affirming the robustness of Enhancer-MDLF. Additionally, we introduce transfer learning to provide an effective and potential solution to address the prediction challenges posed by enhancer specificity. Furthermore, we utilize model interpretation to identify transcription factor binding site motifs that may be associated with enhancer regions, with important implications for facilitating the study of enhancer regulatory mechanisms. The source code is openly accessible at https://github.com/HaoWuLab-Bioinformatics/Enhancer-MDLF.

Genome-Wide Mapping of Active Regulatory Elements Using ATAC-seq

Active cis-regulatory elements (cREs) in eukaryotes are characterized by nucleosomal depletion and, accordingly, higher accessibility. This property has turned out to be immensely useful for identifying cREs genome-wide and tracking their dynamics across different cellular states and is the basis of numerous methods taking advantage of the preferential enzymatic cleavage/labeling of accessible DNA. ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) has emerged as the most versatile and widely adaptable method and has been widely adopted as the standard tool for mapping open chromatin regions. Here, we discuss the current optimal practices and important considerations for carrying out ATAC-seq experiments, primarily in the context of mammalian systems.

SENIES: DNA Shape Enhanced Two-Layer Deep Learning Predictor for the Identification of Enhancers and Their Strength

Identifying enhancers is a critical task in bioinformatics due to their primary role in regulating gene expression. For this reason, various computational algorithms devoted to enhancer identification have been put forward over the years. More features are extracted from the single DNA sequences to boost the performance. Nevertheless, DNA structural information is neglected, which is an essential factor affecting the binding preferences of transcription factors to regulatory elements like enhancers. Here, we propose SENIES, a DNA shape enhanced deep learning predictor, to identify enhancers and their strength. The predictor consists of two layers where the first layer is for enhancer and non-enhancer identification, and the second layer is for predicting the strength of enhancers. Apart from two common sequence-derived features (i.e., one-hot and k-mer), DNA shape is introduced to describe the 3D structures of DNA sequences. Performance comparison with state-of-the-art methods conducted on public datasets demonstrates the effectiveness and robustness of our predictor. The code implementation of SENIES is publicly available at https://github.com/hlju-liye/SENIES.

iEnhancer-DCLA: using the original sequence to identify enhancers and their strength based on a deep learning framework

Enhancers are small regions of DNA that bind to proteins, which enhance the transcription of genes. The enhancer may be located upstream or downstream of the gene. It is not necessarily close to the gene to be acted on, because the entanglement structure of chromatin allows the positions far apart in the sequence to have the opportunity to contact each other. Therefore, identifying enhancers and their strength is a complex and challenging task. In this article, a new prediction method based on deep learning is proposed to identify enhancers and enhancer strength, called iEnhancer-DCLA. Firstly, we use word2vec to convert k-mers into number vectors to construct an input matrix. Secondly, we use convolutional neural network and bidirectional long short-term memory network to extract sequence features, and finally use the attention mechanism to extract relatively important features. In the task of predicting enhancers and their strengths, this method has improved to a certain extent in most evaluation indexes. In summary, we believe that this method provides new ideas in the analysis of enhancers.

Chromatin Dynamics in Digestive System Cancer: Commander and Regulator

Digestive system tumors have a poor prognosis due to complex anatomy, insidious onset, challenges in early diagnosis, and chemoresistance. Epidemiological statistics has verified that digestive system tumors rank first in tumor-related death. Although a great number of studies are devoted to the molecular biological mechanism, early diagnostic markers, and application of new targeted drugs in digestive system tumors, the therapeutic effect is still not satisfactory. Epigenomic alterations including histone modification and chromatin remodeling are present in human cancers and are now known to cooperate with genetic changes to drive the cancer phenotype. Chromatin is the carrier of genetic information and consists of DNA, histones, non-histone proteins, and a small amount of RNA. Chromatin and nucleosomes control the stability of the eukaryotic genome and regulate DNA processes such as transcription, replication, and repair. The dynamic structure of chromatin plays a key role in this regulatory function. Structural fluctuations expose internal DNA and thus provide access to the nuclear machinery. The dynamic changes are affected by various complexes and epigenetic modifications. Variation of chromatin dynamics produces early and superior regulation of the expression of related genes and downstream pathways, thereby controlling tumor development. Intervention at the chromatin level can change the process of cancer earlier and is a feasible option for future tumor diagnosis and treatment. In this review, we introduced chromatin dynamics including chromatin remodeling, histone modifications, and chromatin accessibility, and current research on chromatin regulation in digestive system tumors was also summarized.

Single-Molecule Multikilobase-Scale Profiling of Chromatin Accessibility Using m6A-SMAC-Seq and m6A-CpG-GpC-SMAC-Seq

A hallmark feature of active cis-regulatory elements (CREs) in eukaryotes is their nucleosomal depletion and, accordingly, higher accessibility to enzymatic treatment. This property has been the basis of a number of sequencing-based assays for genome-wide identification and tracking the activity of CREs across different biological conditions, such as DNAse-seq, ATAC-seq , NOMeseq, and others. However, the fragmentation of DNA inherent to many of these assays and the limited read length of short-read sequencing platforms have so far not allowed the simultaneous measurement of the chromatin accessibility state of CREs located distally from each other. The combination of labeling accessible DNA with DNA modifications and nanopore sequencing has made it possible to develop such assays. Here, we provide a detailed protocol for carrying out the SMAC-seq assay (Single-Molecule long-read Accessible Chromatin mapping sequencing), in its m6A-SMAC-seq and m6A-CpG-GpC-SMAC-seq variants, together with methods for data processing and analysis, and discuss key experimental and analytical considerations for working with SMAC-seq datasets.

The role of contactin-associated protein-like 2 in neurodevelopmental disease and human cerebral cortex evolution

The contactin-associated protein-like 2 (CNTNAP2) gene is associated with multiple neurodevelopmental disorders, including autism spectrum disorder (ASD), intellectual disability (ID), and specific language impairment (SLI). Experimental work has shown that CNTNAP2 is important for neuronal development and synapse formation. There is also accumulating evidence for the differential use of CNTNAP2 in the human cerebral cortex compared with other primates. Here, we review the current literature on CNTNAP2, including what is known about its expression, disease associations, and molecular/cellular functions. We also review the evidence for its role in human brain evolution, such as the presence of eight human accelerated regions (HARs) within the introns of the gene. While progress has been made in understanding the function(s) of CNTNAP2, more work is needed to clarify the precise mechanisms through which CNTNAP2 acts. Such information will be crucial for developing effective treatments for CNTNAP2 patients. It may also shed light on the longstanding question of what makes us human.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

Simplified Epigenome Profiling Using Antibody-tethered Tagmentation

We previously introduced Cleavage Under Targets & Tagmentation (CUT&Tag), an epigenomic profiling method in which antibody tethering of the Tn5 transposase to a chromatin epitope of interest maps specific chromatin features in small samples and single cells. With CUT&Tag, intact cells or nuclei are permeabilized, followed by successive addition of a primary antibody, a secondary antibody, and a chimeric Protein A-Transposase fusion protein that binds to the antibody. Addition of Mg⁺⁺ activates the transposase and inserts sequencing adapters into adjacent DNA in situ. We have since adapted CUT&Tag to also map chromatin accessibility by simply modifying the transposase activation conditions when using histone H3K4me2, H3K4me3, or Serine-5-phosphorylated RNA Polymerase II antibodies. Using these antibodies, we redirect the tagmentation of accessible DNA sites to produce chromatin accessibility maps with exceptionally high signal-to-noise and resolution. All steps from nuclei to amplified sequencing-ready libraries are performed in single PCR tubes using non-toxic reagents and inexpensive equipment, making our simplified strategy for simultaneous chromatin profiling and accessibility mapping suitable for the lab, home workbench, or classroom.

Highly Parallel Quantification and Compartment Localization of Transcription Factors and Nuclear Proteins

Transcription factors and other chromatin-associated proteins are difficult to quantify comprehensively. Here, we combine facile nuclear sub-fractionation with data-independent acquisition mass spectrometry to achieve rapid, sensitive, and highly parallel quantification of the nuclear proteome in human cells. We apply this approach to quantify the response to acute degradation of BET bromodomains, revealing unexpected chromatin regulatory dynamics. The method is simple and enables system-level study of previously inaccessible chromatin and genome regulators.

Interrogating the Accessible Chromatin Landscape of Eukaryote Genomes Using ATAC-seq

The ATAC-seq assay has emerged as the most useful, versatile, and widely adaptable method for profiling accessible chromatin regions and tracking the activity of cis-regulatory elements (cREs) in eukaryotes. Thanks to its great utility, it is now being applied to map active chromatin in the context of a very wide diversity of biological systems and questions. In the course of these studies, considerable experience working with ATAC-seq data has accumulated and a standard set of computational tasks that need to be carried for most ATAC-seq analyses has emerged. Here, we review and provide examples of common such analytical procedures (including data processing, quality control, peak calling, identifying differentially accessible open chromatin regions, and variable transcription factor (TF) motif accessibility) and discuss recommended optimal practices.

Efficient chromatin accessibility mapping in situ by nucleosome-tethered tagmentation

Chromatin accessibility mapping is a powerful approach to identify potential regulatory elements. A popular example is ATAC-seq, whereby Tn5 transposase inserts sequencing adapters into accessible DNA ('tagmentation'). CUT&Tag is a tagmentation-based epigenomic profiling method in which antibody tethering of Tn5 to a chromatin epitope of interest profiles specific chromatin features in small samples and single cells. Here, we show that by simply modifying the tagmentation conditions for histone H3K4me2 or H3K4me3 CUT&Tag, antibody-tethered tagmentation of accessible DNA sites is redirected to produce chromatin accessibility maps that are indistinguishable from the best ATAC-seq maps. Thus, chromatin accessibility maps can be produced in parallel with CUT&Tag maps of other epitopes with all steps from nuclei to amplified sequencing-ready libraries performed in single PCR tubes in the laboratory or on a home workbench. As H3K4 methylation is produced by transcription at promoters and enhancers, our method identifies transcription-coupled accessible regulatory sites.

Long-range single-molecule mapping of chromatin accessibility in eukaryotes

Mapping open chromatin regions has emerged as a widely used tool for identifying active regulatory elements in eukaryotes. However, existing approaches, limited by reliance on DNA fragmentation and short-read sequencing, cannot provide information about large-scale chromatin states or reveal coordination between the states of distal regulatory elements. We have developed a method for profiling the accessibility of individual chromatin fibers, a single-molecule long-read accessible chromatin mapping sequencing assay (SMAC-seq), enabling the simultaneous, high-resolution, single-molecule assessment of chromatin states at multikilobase length scales. Our strategy is based on combining the preferential methylation of open chromatin regions by DNA methyltransferases with low sequence specificity, in this case EcoGII, an N6-methyladenosine (m6A) methyltransferase, and the ability of nanopore sequencing to directly read DNA modifications. We demonstrate that aggregate SMAC-seq signals match bulk-level accessibility measurements, observe single-molecule nucleosome and transcription factor protection footprints, and quantify the correlation between chromatin states of distal genomic elements.

The transcription factor PU.1 mediates enhancer-promoter looping that is required for IL-1β eRNA and mRNA transcription in mouse melanoma and macrophage cell lines

The DNA-binding protein PU.1 is a myeloid lineage-determining and pioneering transcription factor due to its ability to bind "closed" genomic sites and maintain "open" chromatin state for myeloid lineage-specific genes. The precise mechanism of PU.1 in cell type-specific programming is yet to be elucidated. The melanoma cell line B16BL6, although it is nonmyeloid lineage, expressed Toll-like receptors and activated the transcription factor NF-κB upon stimulation by the bacterial cell wall component lipopolysaccharide. However, it did not produce cytokines, such as IL-1β mRNA. Ectopic PU.1 expression induced remodeling of a novel distal enhancer (located ∼10 kbp upstream of the IL-1β transcription start site), marked by nucleosome depletion, enhancer-promoter looping, and histone H3 lysine 27 acetylation (H3K27ac). PU.1 induced enhancer-promoter looping and H3K27ac through two distinct PU.1 regions. These PU.1-dependent events were independently required for subsequent signal-dependent and co-dependent events: NF-κB recruitment and further H3K27ac, both of which were required for enhancer RNA (eRNA) transcription. In murine macrophage RAW264.7 cells, these PU.1-dependent events were constitutively established and readily expressed eRNA and subsequently IL-1β mRNA by lipopolysaccharide stimulation. In summary, this study showed a sequence of epigenetic events in programming IL-1β transcription by the distal enhancer priming and eRNA production mediated by PU.1 and the signal-dependent transcription factor NF-κB.

Evolution of DNAase I Hypersensitive Sites in MHC Regulatory Regions of Primates

It has been challenging to determine the disease-causing variant(s) for most major histocompatibility complex (MHC)-associated diseases. However, it is becoming increasingly clear that regulatory variation is pervasive and a fundamentally important mechanism governing phenotypic diversity and disease susceptibility. We gathered DNase I data from 136 human cells to characterize the regulatory landscape of the MHC region, including 4867 DNase I hypersensitive sites (DHSs). We identified thousands of regulatory elements that have been gained or lost in the human or chimpanzee genomes since their evolutionary divergence. We compared alignments of the DHS across six primates and found 149 DHSs with convincing evidence of positive and/or purifying selection. Of these DHSs, compared to neutral sequences, 24 evolved rapidly in the human lineage. We identified 15 instances of transcription-factor-binding motif gains, such as USF, MYC, MAX, MAFK, STAT1, PBX3, etc, and observed 16 GWAS (genome-wide association study) SNPs associated with diseases within these 24 DHSs using FIMO (Find Individual Motif Occurrences) and UCSC (University of California, Santa Cruz) ChIP-seq data. Combining eQTL and Hi-C data, our results indicated that there were five SNPs located in human gains motifs affecting the corresponding gene's expression, two of which closely matched DHS target genes. In addition, a significant SNP, rs7756521, at genome-wide significant level likely affects DDR expression and represents a causal genetic variant for HIV-1 control. These results indicated that species-specific motif gains or losses of rapidly evolving DHSs in the primate genomes might play a role during adaptation evolution and provided some new evidence for a potentially causal role for these GWAS SNPs.

DHSpred: support-vector-machine-based human DNase I hypersensitive sites prediction using the optimal features selected by random forest

DNase I hypersensitive sites (DHSs) are genomic regions that provide important information regarding the presence of transcriptional regulatory elements and the state of chromatin. Therefore, identifying DHSs in uncharacterized DNA sequences is crucial for understanding their biological functions and mechanisms. Although many experimental methods have been proposed to identify DHSs, they have proven to be expensive for genome-wide application. Therefore, it is necessary to develop computational methods for DHS prediction. In this study, we proposed a support vector machine (SVM)-based method for predicting DHSs, called DHSpred (DNase I Hypersensitive Site predictor in human DNA sequences), which was trained with 174 optimal features. The optimal combination of features was identified from a large set that included nucleotide composition and di- and trinucleotide physicochemical properties, using a random forest algorithm. DHSpred achieved a Matthews correlation coefficient and accuracy of 0.660 and 0.871, respectively, which were 3% higher than those of control SVM predictors trained with non-optimized features, indicating the efficiency of the feature selection method. Furthermore, the performance of DHSpred was superior to that of state-of-the-art predictors. An online prediction server has been developed to assist the scientific community, and is freely available at: http://www.thegleelab.org/DHSpred.html.

Human Accelerated Regions and Other Human-Specific Sequence Variations in the Context of Evolution and Their Relevance for Brain Development

The review discusses, in a format of a timeline, the studies of different types of genetic variants, present in Homo sapiens, but absent in all other primate, mammalian, or vertebrate species, tested so far. The main characteristic of these variants is that they are found in regions of high evolutionary conservation. These sequence variations include single nucleotide substitutions (called human accelerated regions), deletions, and segmental duplications. The rationale for finding such variations in the human genome is that they could be responsible for traits, specific to our species, of which the human brain is the most remarkable. As became obvious, the vast majority of human-specific single nucleotide substitutions are found in noncoding, likely regulatory regions. A number of genes, associated with these human-specific alleles, often through novel enhancer activity, were in fact shown to be implicated in human-specific development of certain brain areas, including the prefrontal cortex. Human-specific deletions may remove regulatory sequences, such as enhancers. Segmental duplications, because of their large size, create new coding sequences, like new functional paralogs. Further functional study of these variants will shed light on evolution of our species, as well as on the etiology of neurodevelopmental disorders.

eFORGE: A Tool for Identifying Cell Type-Specific Signal in Epigenomic Data

Epigenome-wide association studies (EWAS) provide an alternative approach for studying human disease through consideration of non-genetic variants such as altered DNA methylation. To advance the complex interpretation of EWAS, we developed eFORGE (http://eforge.cs.ucl.ac.uk/), a new standalone and web-based tool for the analysis and interpretation of EWAS data. eFORGE determines the cell type-specific regulatory component of a set of EWAS-identified differentially methylated positions. This is achieved by detecting enrichment of overlap with DNase I hypersensitive sites across 454 samples (tissues, primary cell types, and cell lines) from the ENCODE, Roadmap Epigenomics, and BLUEPRINT projects. Application of eFORGE to 20 publicly available EWAS datasets identified disease-relevant cell types for several common diseases, a stem cell-like signature in cancer, and demonstrated the ability to detect cell-composition effects for EWAS performed on heterogeneous tissues. Our approach bridges the gap between large-scale epigenomics data and EWAS-derived target selection to yield insight into disease etiology.

Selecting optimal combinations of transcription factors to promote axon regeneration: Why mechanisms matter

Recovery from injuries to the central nervous system, including spinal cord injury, is constrained in part by the intrinsically low ability of many CNS neurons to mount an effective regenerative growth response. To improve outcomes, it is essential to understand and ultimately reverse these neuron-intrinsic constraints. Genetic manipulation of key transcription factors (TFs), which act to orchestrate production of multiple regeneration-associated genes, has emerged as a promising strategy. It is likely that no single TF will be sufficient to fully restore neuron-intrinsic growth potential, and that multiple, functionally interacting factors will be needed. An extensive literature, mostly from non-neural cell types, has identified potential mechanisms by which TFs can functionally synergize. Here we examine four potential mechanisms of TF/TF interaction; physical interaction, transcriptional cross-regulation, signaling-based cross regulation, and co-occupancy of regulatory DNA. For each mechanism, we consider how existing knowledge can be used to guide the discovery and effective use of TF combinations in the context of regenerative neuroscience. This mechanistic insight into TF interactions is needed to accelerate the design of effective TF-based interventions to relieve neuron-intrinsic constraints to regeneration and to foster recovery from CNS injury.

Repression of the Heat Shock Response Is a Programmed Event at the Onset of Reproduction

The heat shock response (HSR) is essential for proteostasis and cellular health. In metazoans, aging is associated with a decline in quality control, thus increasing the risk for protein conformational disease. Here, we show that in C. elegans, the HSR declines precipitously over a 4 hr period in early adulthood coincident with the onset of reproductive maturity. Repression of the HSR occurs due to an increase in H3K27me3 marks at stress gene loci, the timing of which is determined by reduced expression of the H3K27 demethylase jmjd-3.1. This results in a repressed chromatin state that interferes with HSF-1 binding and suppresses transcription initiation in response to stress. The removal of germline stem cells preserves jmjd-3.1 expression, suppresses the accumulation of H3K27me3 at stress gene loci, and maintains the HSR. These findings suggest that competing requirements of the germline and soma dictate organismal stress resistance as animals begin reproduction.

Comprehensive identification and analysis of human accelerated regulatory DNA

It has long been hypothesized that changes in gene regulation have played an important role in human evolution, but regulatory DNA has been much more difficult to study compared with protein-coding regions. Recent large-scale studies have created genome-scale catalogs of DNase I hypersensitive sites (DHSs), which demark potentially functional regulatory DNA. To better define regulatory DNA that has been subject to human-specific adaptive evolution, we performed comprehensive evolutionary and population genetics analyses on over 18 million DHSs discovered in 130 cell types. We identified 524 DHSs that are conserved in nonhuman primates but accelerated in the human lineage (haDHS), and estimate that 70% of substitutions in haDHSs are attributable to positive selection. Through extensive computational and experimental analyses, we demonstrate that haDHSs are often active in brain or neuronal cell types; play an important role in regulating the expression of developmentally important genes, including many transcription factors such as SOX6, POU3F2, and HOX genes; and identify striking examples of adaptive regulatory evolution that may have contributed to human-specific phenotypes. More generally, our results reveal new insights into conserved and adaptive regulatory DNA in humans and refine the set of genomic substrates that distinguish humans from their closest living primate relatives.

Chromatin accessibility: a window into the genome

Transcriptional activation throughout the eukaryotic lineage has been tightly linked with disruption of nucleosome organization at promoters, enhancers, silencers, insulators and locus control regions due to transcription factor binding. Regulatory DNA thus coincides with open or accessible genomic sites of remodeled chromatin. Current chromatin accessibility assays are used to separate the genome by enzymatic or chemical means and isolate either the accessible or protected locations. The isolated DNA is then quantified using a next-generation sequencing platform. Wide application of these assays has recently focused on the identification of the instrumental epigenetic changes responsible for differential gene expression, cell proliferation, functional diversification and disease development. Here we discuss the limitations and advantages of current genome-wide chromatin accessibility assays with especial attention on experimental precautions and sequence data analysis. We conclude with our perspective on future improvements necessary for moving the field of chromatin profiling forward.

The role of BAF (mSWI/SNF) complexes in mammalian neural development

The BAF (mammalian SWI/SNF) complexes are a family of multi-subunit ATP-dependent chromatin remodelers that use ATP hydrolysis to alter chromatin structure. Distinct BAF complex compositions are possible through combinatorial assembly of homologous subunit families and can serve non-redundant functions. In mammalian neural development, developmental stage-specific BAF assemblies are found in embryonic stem cells, neural progenitors and postmitotic neurons. In particular, the neural progenitor-specific BAF complexes are essential for controlling the kinetics and mode of neural progenitor cell division, while neuronal BAF function is necessary for the maturation of postmitotic neuronal phenotypes as well as long-term memory formation. The microRNA-mediated mechanism for transitioning from npBAF to nBAF complexes is instructive for the neuronal fate and can even convert fibroblasts into neurons. The high frequency of BAF subunit mutations in neurological disorders underscores the rate-determining role of BAF complexes in neural development, homeostasis, and plasticity.

Genome-scale mapping of DNase I hypersensitivity

DNase I-seq is a global and high-resolution method that uses the nonspecific endonuclease DNase I to map chromatin accessibility. These accessible regions, designated as DNase I hypersensitive sites (DHSs), define the regulatory features, (e.g., promoters, enhancers, insulators, and locus control regions) of complex genomes. In this unit, methods are described for nuclei isolation, digestion of nuclei with limiting concentrations of DNase I, and the biochemical fractionation of DNase I hypersensitive sites in preparation for high-throughput sequencing. DNase I-seq is an unbiased and robust method that is not predicated on an a priori understanding of regulatory patterns or chromatin features.

Regulation of gene expression: probing DNA–protein interactionsin vivoandin vitro

Tremendous efforts have been put together over the last several years to complete the entire sequencing of the human genome. As we enter the proteomic era, when the major aim is understanding which gene encodes which protein, the time has also come to identify their precise function inside the astonishing signaling network required to accomplish all cellular functions. Understanding when, why and how a gene is expressed has now become a necessity toward identifying all the regulatory pathways that mediate cellular processes such as differentiation, migration, replication, DNA repair and apoptosis. Regulation of gene transcription is a process that is primarily under the influence of nuclear-located transcription factors. Consequently, identifying which protein activates or represses a specific gene is a prerequisite for understanding cell fate and function. The current state of, and recent advances in, transcriptional regulation approaches are reviewed here, with special emphasis on new technologies required when probing for DNA-protein interactions. This review explores different strategies aimed at identifying both the regulatory sequences of any given gene and the trans-acting regulatory factors that recognize these elements as their target sites in the nucleus. Ongoing developments in the fields of nanotechnology, RNA silencing and protein modeling toward the investigation of DNA-protein interactions and their relevance in the battle against cancer are discussed.

A common set of DNA regulatory elements shapes Drosophila appendages

Animals have body parts made of similar cell types located at different axial positions, such as limbs. The identity and distinct morphology of each structure is often specified by the activity of different "master regulator" transcription factors. Although similarities in gene expression have been observed between body parts made of similar cell types, how regulatory information in the genome is differentially utilized to create morphologically diverse structures in development is not known. Here, we use genome-wide open chromatin profiling to show that among the Drosophila appendages, the same DNA regulatory modules are accessible throughout the genome at a given stage of development, except at the loci encoding the master regulators themselves. In addition, open chromatin profiles change over developmental time, and these changes are coordinated between different appendages. We propose that master regulators create morphologically distinct structures by differentially influencing the function of the same set of DNA regulatory modules.

Enhancing myogenic differentiation of pluripotent stem cells with small molecule inducers

Pluripotent stem cells are able to differentiate into many types of cell lineages in response to differentiation cues. However, a pure population of lineage-specific cells is desirable for any potential clinical application. Therefore, induction of the pluripotent stem cells with lineage-specific regulatory signals, or small molecule inducers, is a prerequisite for effectively directing lineage specification for cell-based therapeutics. In this article, we provide in-depth analysis of recent research findings on small molecule inducers of the skeletal muscle lineage. We also provide perspectives on how different signaling pathways and chromatin dynamics converge to direct the differentiation of skeletal myocytes.

Enhancers: five essential questions

It is estimated that the human genome contains hundreds of thousands of enhancers, so understanding these gene-regulatory elements is a crucial goal. Several fundamental questions need to be addressed about enhancers, such as how do we identify them all, how do they work, and how do they contribute to disease and evolution? Five prominent researchers in this field look at how much we know already and what needs to be done to answer these questions.

Probing DNA shape and methylation state on a genomic scale with DNase I

DNA binding proteins find their cognate sequences within genomic DNA through recognition of specific chemical and structural features. Here we demonstrate that high-resolution DNase I cleavage profiles can provide detailed information about the shape and chemical modification status of genomic DNA. Analyzing millions of DNA backbone hydrolysis events on naked genomic DNA, we show that the intrinsic rate of cleavage by DNase I closely tracks the width of the minor groove. Integration of these DNase I cleavage data with bisulfite sequencing data for the same cell type's genome reveals that cleavage directly adjacent to cytosine-phosphate-guanine (CpG) dinucleotides is enhanced at least eightfold by cytosine methylation. This phenomenon we show to be attributable to methylation-induced narrowing of the minor groove. Furthermore, we demonstrate that it enables simultaneous mapping of DNase I hypersensitivity and regional DNA methylation levels using dense in vivo cleavage data. Taken together, our results suggest a general mechanism by which CpG methylation can modulate protein-DNA interaction strength via the remodeling of DNA shape.

Technical considerations for functional sequencing assays

Next-generation sequencing technologies need careful design of experiments and evaluation of results to meet field requirements. Here we discuss technical considerations for these high-throughput assays, together with criteria to assess the quality of the results and the necessary validation.

Genomic approaches for the discovery of CFTR regulatory elements

Non-coding regions of the human genome contain vast regulatory potential that contributes to the coordination of gene expression. Indeed, regulatory elements can reside large genomic distances from the promoters of genes they control. Here we describe approaches recently used to identify functional elements within the complex CFTR locus.

Evidence for an intronic cis-regulatory element within CD69 gene

CD69 is one of the earliest proteins expressed after leukocyte activation and its engagement is essential in the control of innate and adaptive immune responses. Inducible CD69 expression is strongly controlled at the transcriptional level. The molecular basis for developmental- and stage-specific regulation in T cells is beginning to be elucidated while it remains largely unknown in the rest of immune cells. DNase I hypersensitivity experiments in lymphocytes identified a novel hypersensitive region within mouse and human intron I, which was inducible upon stimulation. In silico analysis of CD69 gene revealed that this open chromatin region was present in different cell types and was associated with positioned nucleosomes. Analysis of histone post-translational modifications of intron I indicated that acetylation and lysine 4 dimethylation of histone H3 were dynamically regulated during thymocyte development and were constitutively high in resting and stimulated mature T lymphocytes. Thus, we provide evidence for the existence of a cis-acting element in intron I that is more accessible to DNase I digestion and that it is developmentally regulated at the chromatin level.

Detecting DNaseI-hypersensitivity sites with MLPA

DNaseI-hypersensitive sites within chromatin are indicative of genomic loci with regulatory function. Several techniques have been described for analyzing these regions, but are either laborious, offer low-throughput possibilities, or are expensive. We have developed a new approach based on a modified version of multiplex ligation-dependent probe amplification (MLPA). Using this method, it is possible to analyse up to 50 defined genomic regions for DNaseI-hypersensitivity in a single PCR-based reaction. This chapter outlines the approach and discusses the critical features of each step of the procedure.

Regulation of Myf5 Early Enhancer by Histone Acetyltransferase p300 during Stem Cell Differentiation

Skeletal myogenesis is an intricate process coordinated temporally by multiple myogenic regulatory factors (MRF) including Myf5, which is the first MRF expressed and marks the commitment of skeletal muscle lineage. The expression of Myf5 gene during early embryogenesis is controlled by a set of enhancer elements, and requires the histone acetyltransferase (HAT) activity of transcriptional coactivator p300. However, it is unclear as to how different regulatory signals converge at enhancer elements to regulate early Myf5 gene expression, and if p300 is directly involved. We show here that p300 associates with the Myf5 early enhancer at the early stage of stem cell differentiation, and its HAT activity is important for the recruitment of β-catenin to this early enhancer. In addition, histone H3-K27 acetylation, but not H3-K9/14, is intimately connected to the p300 HAT activity. Thus, p300 is directly involved in the regulation of the Myf5 early enhancer, and is important for specific histone acetylation and transcription factor recruitment. This connection of p300 HAT activity with H3-K27 acetylation and β-catenin signalling during myogenic differentiation in vitro offers a molecular insight into the enhancer-elements participation observed in embryonic development. In addition, pluripotent stem cell differentiation is a valuable system to dissect the signal-dependent regulation of specific enhancer element during cell fate determinations.

Rapid empirical discovery of optimal peptides for targeted proteomics

We report a method for high-throughput, cost-efficient empirical discovery of optimal proteotypic peptides and fragment ions for targeted proteomics applications using in vitro-synthesized proteins. We demonstrate the approach using human transcription factors – which are typically difficult, low-abundance – targets with an overall success rate of 98%. We show further that targeted proteomic assays developed using our approach facilitate robust in vivo quantification of human transcription factors.

A two-stage evolutionary approach for effective classification of hypersensitive DNA sequences

Hypersensitive (HS) sites in genomic sequences are reliable markers of DNA regulatory regions that control gene expression. Annotation of regulatory regions is important in understanding phenotypical differences among cells and diseases linked to pathologies in protein expression. Several computational techniques are devoted to mapping out regulatory regions in DNA by initially identifying HS sequences. Statistical learning techniques like Support Vector Machines (SVM), for instance, are employed to classify DNA sequences as HS or non-HS. This paper proposes a method to automate the basic steps in designing an SVM that improves the accuracy of such classification. The method proceeds in two stages and makes use of evolutionary algorithms. An evolutionary algorithm first designs optimal sequence motifs to associate explicit discriminating feature vectors with input DNA sequences. A second evolutionary algorithm then designs SVM kernel functions and parameters that optimally separate the HS and non-HS classes. Results show that this two-stage method significantly improves SVM classification accuracy. The method promises to be generally useful in automating the analysis of biological sequences, and we post its source code on our website.

Assessing the effects of symmetry on motif discovery and modeling

Background

Identifying the DNA binding sites for transcription factors is a key task in modeling the gene regulatory network of a cell. Predicting DNA binding sites computationally suffers from high false positives and false negatives due to various contributing factors, including the inaccurate models for transcription factor specificity. One source of inaccuracy in the specificity models is the assumption of asymmetry for symmetric models.

Methodology/Principal Findings

Using simulation studies, so that the correct binding site model is known and various parameters of the process can be systematically controlled, we test different motif finding algorithms on both symmetric and asymmetric binding site data. We show that if the true binding site is asymmetric the results are unambiguous and the asymmetric model is clearly superior to the symmetric model. But if the true binding specificity is symmetric commonly used methods can infer, incorrectly, that the motif is asymmetric. The resulting inaccurate motifs lead to lower sensitivity and specificity than would the correct, symmetric models. We also show how the correct model can be obtained by the use of appropriate measures of statistical significance.

Conclusions/Significance

This study demonstrates that the most commonly used motif-finding approaches usually model symmetric motifs incorrectly, which leads to higher than necessary false prediction errors. It also demonstrates how alternative motif-finding methods can correct the problem, providing more accurate motif models and reducing the errors. Furthermore, it provides criteria for determining whether a symmetric or asymmetric model is the most appropriate for any experimental dataset.

The complex transcription regulatory landscape of our genome: control in three dimensions

The non-coding part of our genome contains sequence motifs that can control gene transcription over distance. Here, we discuss functional genomics studies that uncover and characterize these sequences across the mammalian genome. The picture emerging is of a genome being a complex regulatory landscape. We explore the principles that underlie the wiring of regulatory DNA sequences and genes. We argue transcriptional control over distance can be understood when considering action in the context of the folded genome. Genome topology is expected to differ between individual cells, and this may cause variegated expression. High-resolution three-dimensional genome topology maps, ultimately of single cells, are required to understand the cis-regulatory networks that underlie cellular transcriptomes.

esBAF facilitates pluripotency by conditioning the genome for LIF/STAT3 signalling and by regulating polycomb function

Signaling by the cytokine LIF and its downstream transcription factor, STAT3, prevents differentiation of pluripotent embryonic stem cells (ESCs) by opposing MAP kinase signaling. This contrasts with most cell types where STAT3signaling induces differentiation. We find that STAT3binding across the pluripotent genome is dependent upon Brg, the ATPase subunit of a specialized chromatin remodeling complex (esBAF) found in ESCs. Brg is required to establish chromatin accessibility at STAT3 binding targets, in essence preparing these sites to respond to LIF signaling. Moreover, Brg deletion leads to rapid Polycomb (PcG) binding and H3K27me3-mediated silencing of many Brg-activated targets genome-wide, including the target genes of the LIF signaling pathway. Hence, one crucial role of Brg in ESCs involves its ability to potentiate LIF signaling by opposing PcG. Contrary to expectations, Brg also facilitates PcG function at classical PcG target including all four Hox loci, reinforcing their repression in ESCs. These findings reveal that esBAF does not simply antagonize PcG, but rather, the two chromatin regulators act both antagonistically and synergistically with the common goal of supporting pluripotency.

Extensive chromatin remodelling and establishment of transcription factor 'hotspots' during early adipogenesis

Adipogenesis is tightly controlled by a complex network of transcription factors acting at different stages of differentiation. Peroxisome proliferator-activated receptor γ (PPARγ) and CCAAT/enhancer-binding protein (C/EBP) family members are key regulators of this process. We have employed DNase I hypersensitive site analysis to investigate the genome-wide changes in chromatin structure that accompany the binding of adipogenic transcription factors. These analyses revealed a dramatic and dynamic modulation of the chromatin landscape during the first hours of adipocyte differentiation that coincides with cooperative binding of multiple early transcription factors (including glucocorticoid receptor, retinoid X receptor, Stat5a, C/EBPβ and -δ) to transcription factor 'hotspots'. Our results demonstrate that C/EBPβ marks a large number of these transcription factor 'hotspots' before induction of differentiation and chromatin remodelling and is required for their establishment. Furthermore, a subset of early remodelled C/EBP-binding sites persists throughout differentiation and is later occupied by PPARγ, indicating that early C/EBP family members, in addition to their well-established role in activation of PPARγ transcription, may act as pioneering factors for PPARγ binding.

Chromatin remodeling complex NURF regulates thymocyte maturation

The maturation of T cells requires signaling from both cytokine and T-cell receptors to gene targets in chromatin, but how chromatin architecture influences this process is largely unknown. Here we show that thymocyte maturation post-positive selection is dependent on the nucleosome remodeling factor (NURF). Depletion of Bptf (bromodomain PHD finger transcription factor), the largest NURF subunit, in conditional mouse mutants results in developmental arrest beyond the CD4(+) CD8(int) stage without affecting cellular proliferation, cellular apoptosis, or coreceptor gene expression. In the Bptf mutant, specific subsets of genes important for thymocyte development show aberrant expression. We also observed defects in DNase I-hypersensitive chromatin structures at Egr1, a prototypical Bptf-dependent gene that is required for efficient thymocyte development. Moreover, chromatin binding of the sequence-specific factor Srf (serum response factor) to Egr1 regulatory sites is dependent on Bptf function. Physical interactions between NURF and Srf suggest a model in which Srf recruits NURF to facilitate transcription factor binding at Bptf-dependent genes. These findings provide evidence for causal connections between NURF, transcription factor occupancy, and gene regulation during thymocyte development.

A dynamic Bayesian network for identifying protein-binding footprints from single molecule-based sequencing data

Motivation: A global map of transcription factor binding sites (TFBSs) is critical to understanding gene regulation and genome function. DNaseI digestion of chromatin coupled with massively parallel sequencing (digital genomic footprinting) enables the identification of protein-binding footprints with high resolution on a genome-wide scale. However, accurately inferring the locations of these footprints remains a challenging computational problem.

Results: We present a dynamic Bayesian network-based approach for the identification and assignment of statistical confidence estimates to protein-binding footprints from digital genomic footprinting data. The method, DBFP, allows footprints to be identified in a probabilistic framework and outperforms our previously described algorithm in terms of precision at a fixed recall. Applied to a digital footprinting data set from Saccharomyces cerevisiae, DBFP identifies 4679 statistically significant footprints within intergenic regions. These footprints are mainly located near transcription start sites and are strongly enriched for known TFBSs. Footprints containing no known motif are preferentially located proximal to other footprints, consistent with cooperative binding of these footprints. DBFP also identifies a set of statistically significant footprints in the yeast coding regions. Many of these footprints coincide with the boundaries of antisense transcripts, and the most significant footprints are enriched for binding sites of the chromatin-associated factors Abf1 and Rap1.

Contact:jay.hesselberth@ucdenver.edu; william-noble@u.washington.edu

Supplementary information:Supplementary material is available at Bioinformatics online.

Next-generation genomics: an integrative approach

Integrating results from diverse experiments is an essential process in our effort to understand the logic of complex systems, such as development, homeostasis and responses to the environment. With the advent of high-throughput methods--including genome-wide association (GWA) studies, chromatin immunoprecipitation followed by sequencing (ChIP-seq) and RNA sequencing (RNA-seq)--acquisition of genome-scale data has never been easier. Epigenomics, transcriptomics, proteomics and genomics each provide an insightful, and yet one-dimensional, view of genome function; integrative analysis promises a unified, global view. However, the large amount of information and diverse technology platforms pose multiple challenges for data access and processing. This Review discusses emerging issues and strategies related to data integration in the era of next-generation genomics.

Cell-type-specific long-range looping interactions identify distant regulatory elements of the CFTR gene

Identification of regulatory elements and their target genes is complicated by the fact that regulatory elements can act over large genomic distances. Identification of long-range acting elements is particularly important in the case of disease genes as mutations in these elements can result in human disease. It is becoming increasingly clear that long-range control of gene expression is facilitated by chromatin looping interactions. These interactions can be detected by chromosome conformation capture (3C). Here, we employed 3C as a discovery tool for identification of long-range regulatory elements that control the cystic fibrosis transmembrane conductance regulator gene, CFTR. We identified four elements in a 460-kb region around the locus that loop specifically to the CFTR promoter exclusively in CFTR expressing cells. The elements are located 20 and 80 kb upstream; and 109 and 203 kb downstream of the CFTR promoter. These elements contain DNase I hypersensitive sites and histone modification patterns characteristic of enhancers. The elements also interact with each other and the latter two activate the CFTR promoter synergistically in reporter assays. Our results reveal novel long-range acting elements that control expression of CFTR and suggest that 3C-based approaches can be used for discovery of novel regulatory elements.

Fine mapping of chromatin structure in Drosophila melanogaster embryos using micrococcal nuclease

The structure of chromatin in eukaryotes exerts significant influences on many DNA related processes, including transcription, replication, recombination and repair. A useful tool for mapping chromatin structure is micrococcal nuclease (MNase), which induces double-strand breaks within nucleosome linker regions, and with more extensive digestion, single-strand nicks within the nucleosome itself. Many studies, carried out largely with microbes and cell cultures, have used MNase to determine the positions of nucleosomes within a region of DNA to identify dynamic changes induced during gene regulation. To measure similar processes in a developmental context, we turned to a tractable model system, the Drosophila embryo. Here we describe a protocol that enables MNase mapping of the enhancer chromatin structure in the embryo, and show how it can be used to identify structural changes on a cis-regulatory element targeted by the Knirps repressor.

Integration of distinct intracellular signaling pathways at distal regulatory elements directs T-bet expression in human CD4+ T cells

T-bet is a key regulator controlling Th1 cell development. This factor is not expressed in naive CD4(+) T cells, and the mechanisms controlling expression of T-bet are incompletely understood. In this study, we defined regulatory elements at the human T-bet locus and determined how signals originating at the TCR and at cytokine receptors are integrated to induce chromatin modifications and expression of this gene during human Th1 cell differentiation. We found that T cell activation induced two strong DNase I-hypersensitive sites (HS) and rapid histone acetylation at these elements in CD4(+) T cells. Histone acetylation and T-bet expression were strongly inhibited by cyclosporine A, and we detected binding of NF-AT to a HS in vivo. IL-12 and IFN-gamma signaling alone were not sufficient to induce T-bet expression in naive CD4(+) T cells, but enhanced T-bet expression in TCR/CD28-stimulated cells. We detected a third HS 12 kb upstream of the mRNA start site only in developing Th1 cells, which was bound by IL-12-induced STAT4. Our data suggest that T-bet locus remodeling and gene expression are initiated by TCR-induced NF-AT recruitment and amplified by IL-12-mediated STAT4 binding to distinct distal regulatory elements during human Th1 cell differentiation.