Genome-scale mapping of DNase I hypersensitivity

Genomic context sensitizes regulatory elements to genetic disruption

Enhancer function is frequently investigated piecemeal using truncated reporter assays or single deletion analysis. Thus it remains unclear to what extent enhancer function at native loci relies on surrounding genomic context. Using the Big-IN technology for targeted integration of large DNAs, we analyzed the regulatory architecture of the murine Igf2/H19 locus, a paradigmatic model of enhancer selectivity. We assembled payloads containing a 157-kb functional Igf2/H19 locus and engineered mutations to genetically direct CTCF occupancy at the imprinting control region (ICR) that switches the target gene of the H19 enhancer cluster. Contrasting activity of payloads delivered at the endogenous Igf2/H19 locus or ectopically at Hprt revealed that the Igf2/H19 locus includes additional, previously unknown long-range regulatory elements. Exchanging components of the Igf2/H19 locus with the well-studied Sox2 locus showed that the H19 enhancer cluster functioned poorly out of context, and required its native surroundings to activate Sox2 expression. Conversely, the Sox2 locus control region (LCR) could activate both Igf2 and H19 outside its native context, but its activity was only partially modulated by CTCF occupancy at the ICR. Analysis of regulatory DNA actuation across different cell types revealed that, while the H19 enhancers are tightly coordinated within their native locus, the Sox2 LCR acts more independently. We show that these enhancer clusters typify broader classes of loci genome-wide. Our results show that unexpected dependencies may influence even the most studied functional elements, and our synthetic regulatory genomics approach permits large-scale manipulation of complete loci to investigate the relationship between locus architecture and function.

XL-DNase-Seq: Footprinting Analysis of Dynamic Transcription Factors

We have developed a novel method for genomic footprinting of transcription factors (TFs) that detects potential gene regulatory relationships from DNase-seq data at the nucleotide level. We introduce an assay termed cross-link (XL)-DNase-seq, designed to capture chromatin interactions of dynamic TFs. A mild cross-linking step in XL-DNase-seq improves the detection of DNase-based footprints of dynamic TFs. The footprint strengths and detectability depend on an optimal cross-linking procedure. This method may help extract novel gene regulatory circuits involving previously undetectable TFs. The XL-DNase-seq method is illustrated here for activated mouse macrophage-like cells, which share several features with inflammatory macrophages.

The ENCODE Uniform Analysis Pipelines

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

The ENCODE Uniform Analysis Pipelines

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

Typical Enhancers, Super-Enhancers, and Cancers

Simple Summary

The cancer genome has been exhaustively studied upon the advent of Next-Generation Sequencing technologies. Coding and non-coding sequences have been defined as hotspots of genomic variations that affect the naïve gene expression programs established in normal cells, thus working as endogenous drivers of carcinogenesis. In this review, we comprehensively summarize fundamental aspects of gene expression regulation, with emphasis on the impact of sequence and structural variations mapped across non-coding cis-acting elements of genes encoding for tumor-related transcription factors. Chromatin architecture, epigenome reprogramming, transcriptional enhancers and Super-enhancers, oncogene regulation, cutting-edge technologies, and pharmacological treatment are substantially highlighted.

Abstract

Non-coding segments of the human genome are enriched in cis-regulatory modules that constitute functional elements, such as transcriptional enhancers and Super-enhancers. A hallmark of cancer pathogenesis is the dramatic dysregulation of the “archetype” gene expression profiles of normal human cells. Genomic variations can promote such deficiencies when occurring across enhancers and Super-enhancers, since they affect their mechanistic principles, their functional capacity and specificity, and the epigenomic features of the chromatin microenvironment across which these regulatory elements reside. Here, we comprehensively describe: fundamental mechanisms of gene expression dysregulation in cancers that involve genomic abnormalities within enhancers’ and Super-enhancers’ (SEs) sequences, which alter the expression of oncogenic transcription factors (TFs); cutting-edge technologies applied for the analysis of variation-enriched hotspots of the cancer genome; and pharmacological approaches for the treatment of Super-enhancers’ aberrant function. Finally, we provide an intratumor meta-analysis, which highlights that genomic variations in transcription-factor-driven tumors are accompanied overexpression of genes, a portion of which encodes for additional cancer-related transcription factors.

An effector index to predict target genes at GWAS loci

Drug development and biological discovery require effective strategies to map existing genetic associations to causal genes. To approach this problem, we selected 12 common diseases and quantitative traits for which highly powered genome-wide association studies (GWAS) were available. For each disease or trait, we systematically curated positive control gene sets from Mendelian forms of the disease and from targets of medicines used for disease treatment. We found that these positive control genes were highly enriched in proximity of GWAS-associated single-nucleotide variants (SNVs). We then performed quantitative assessment of the contribution of commonly used genomic features, including open chromatin maps, expression quantitative trait loci (eQTL), and chromatin conformation data. Using these features, we trained and validated an Effector Index (Ei), to map target genes for these 12 common diseases and traits. Ei demonstrated high predictive performance, both with cross-validation on the training set, and an independently derived set for type 2 diabetes. Key predictive features included coding or transcript-altering SNVs, distance to gene, and open chromatin-based metrics. This work outlines a simple, understandable approach to prioritize genes at GWAS loci for functional follow-up and drug development, and provides a systematic strategy for prioritization of GWAS target genes.

Analysis of Enhancers and Transcriptional Networks in Thermogenic Adipocytes

Transcription factor (TF) networks orchestrate the regulation of gene programs in mammalian cells, including white and brown adipocytes. In this protocol, we outline how genomics and transcriptomics data can be integrated to infer causal TFs of a given cellular response or cell type using "Integrated analysis of Motif Activity and Gene Expression changes of transcription factors" (IMAGE). Here, we show how key regulatory TFs controlling white and brown adipocyte gene programs can be predicted from chromatin accessibility and RNA-seq data. Furthermore, we demonstrate how information about target sites and target genes of the predicted key regulators can be integrated to propose testable hypotheses regarding the role and mechanisms of TFs.

Discrete regulatory modules instruct hematopoietic lineage commitment and differentiation

Lineage commitment and differentiation is driven by the concerted action of master transcriptional regulators at their target chromatin sites. Multiple efforts have characterized the key transcription factors (TFs) that determine the various hematopoietic lineages. However, the temporal interactions between individual TFs and their chromatin targets during differentiation and how these interactions dictate lineage commitment remains poorly understood. Here we perform dense, daily, temporal profiling of chromatin accessibility (DNase I-seq) and gene expression changes (total RNA-seq) along ex vivo human erythropoiesis to comprehensively define developmentally regulated DNase I hypersensitive sites (DHSs) and transcripts. We link both distal DHSs to their target gene promoters and individual TFs to their target DHSs, revealing that the regulatory landscape is organized in distinct sequential regulatory modules that regulate lineage restriction and maturation. Finally, direct comparison of transcriptional dynamics (bulk and single-cell) and lineage potential between erythropoiesis and megakaryopoiesis uncovers differential fate commitment dynamics between the two lineages as they exit the stem and progenitor stage. Collectively, these data provide insights into the temporally regulated synergy of the cis- and the trans-regulatory components underlying hematopoietic lineage commitment and differentiation.

Tissue context determines the penetrance of regulatory DNA variation

Functional assessment of disease-associated sequence variation at non-coding regulatory elements is complicated by their high degree of context sensitivity to both the local chromatin and nuclear environments. Allelic profiling of DNA accessibility across individuals has shown that only a select minority of sequence variation affects transcription factor (TF) occupancy, yet low sequence diversity in human populations means that no experimental assessment is available for the majority of disease-associated variants. Here we describe high-resolution in vivo maps of allelic DNA accessibility in liver, kidney, lung and B cells from 5 increasingly diverged strains of F1 hybrid mice. The high density of heterozygous sites in these hybrids enables precise quantification of effect size and cell-type specificity for hundreds of thousands of variants throughout the mouse genome. We show that chromatin-altering variants delineate characteristic sensitivity profiles for hundreds of TF motifs. We develop a compendium of TF-specific sensitivity profiles accounting for genomic context effects. Finally, we link maps of allelic accessibility to allelic transcript levels in the same samples. This work provides a foundation for quantitative prediction of cell-type specific effects of non-coding variation on TF activity, which will facilitate both fine-mapping and systems-level analyses of common disease-associated variation in human genomes.

Molecular and computational approaches to map regulatory elements in 3D chromatin structure

Epigenetic marks do not change the sequence of DNA but affect gene expression in a cell-type specific manner by altering the activities of regulatory elements. Development of new molecular biology assays, sequencing technologies, and computational approaches enables us to profile the human epigenome in three-dimensional structure genome-wide. Here we describe various molecular biology techniques and bioinformatic tools that have been developed to measure the activities of regulatory elements and their chromatin interactions. Moreover, we list currently available three-dimensional epigenomic data sets that are generated in various human cell types and tissues to assist in the design and analysis of research projects.

MEDEA: analysis of transcription factor binding motifs in accessible chromatin

Deciphering the interplay between chromatin accessibility and transcription factor (TF) binding is fundamental to understanding transcriptional regulation, control of cellular states, and the establishment of new phenotypes. Recent genome-wide chromatin accessibility profiling studies have provided catalogs of putative open regions, where TFs can recognize their motifs and regulate gene expression programs. Here, we present motif enrichment in differential elements of accessibility (MEDEA), a computational tool that analyzes high-throughput chromatin accessibility genomic data to identify cell-type-specific accessible regions and lineage-specific motifs associated with TF binding therein. To benchmark MEDEA, we used a panel of reference cell lines profiled by ENCODE and curated by the ENCODE Project Consortium for the ENCODE-DREAM Challenge. By comparing results with RNA-seq data, ChIP-seq peaks, and DNase-seq footprints, we show that MEDEA improves the detection of motifs associated with known lineage specifiers. We then applied MEDEA to 610 ENCODE DNase-seq data sets, where it revealed significant motifs even when absolute enrichment was low and where it identified novel regulators, such as NRF1 in kidney development. Finally, we show that MEDEA performs well on both bulk and single-cell ATAC-seq data. MEDEA is publicly available as part of our Glossary-GENRE suite for motif enrichment analysis.

FactorNet: A deep learning framework for predicting cell type specific transcription factor binding from nucleotide-resolution sequential data

Due to the large numbers of transcription factors (TFs) and cell types, querying binding profiles of all valid TF/cell type pairs is not experimentally feasible. To address this issue, we developed a convolutional-recurrent neural network model, called FactorNet, to computationally impute the missing binding data. FactorNet trains on binding data from reference cell types to make predictions on testing cell types by leveraging a variety of features, including genomic sequences, genome annotations, gene expression, and signal data, such as DNase I cleavage. FactorNet implements several convenient strategies to reduce runtime and memory consumption. By visualizing the neural network models, we can interpret how the model predicts binding. We also investigate the variables that affect cross-cell type accuracy, and offer suggestions to improve upon this field. Our method ranked among the top teams in the ENCODE-DREAM in vivo Transcription Factor Binding Site Prediction Challenge, achieving first place on six of the 13 final round evaluation TF/cell type pairs, the most of any competing team. The FactorNet source code is publicly available, allowing users to reproduce our methodology from the ENCODE-DREAM Challenge.

Beyond the Exome: The Non-coding Genome and Enhancers in Neurodevelopmental Disorders and Malformations of Cortical Development

The development of the human cerebral cortex is a complex and dynamic process, in which neural stem cell proliferation, neuronal migration, and post-migratory neuronal organization need to occur in a well-organized fashion. Alterations at any of these crucial stages can result in malformations of cortical development (MCDs), a group of genetically heterogeneous neurodevelopmental disorders that present with developmental delay, intellectual disability and epilepsy. Recent progress in genetic technologies, such as next generation sequencing, most often focusing on all protein-coding exons (e.g., whole exome sequencing), allowed the discovery of more than a 100 genes associated with various types of MCDs. Although this has considerably increased the diagnostic yield, most MCD cases remain unexplained. As Whole Exome Sequencing investigates only a minor part of the human genome (1-2%), it is likely that patients, in which no disease-causing mutation has been identified, could harbor mutations in genomic regions beyond the exome. Even though functional annotation of non-coding regions is still lagging behind that of protein-coding genes, tremendous progress has been made in the field of gene regulation. One group of non-coding regulatory regions are enhancers, which can be distantly located upstream or downstream of genes and which can mediate temporal and tissue-specific transcriptional control via long-distance interactions with promoter regions. Although some examples exist in literature that link alterations of enhancers to genetic disorders, a widespread appreciation of the putative roles of these sequences in MCDs is still lacking. Here, we summarize the current state of knowledge on cis-regulatory regions and discuss novel technologies such as massively-parallel reporter assay systems, CRISPR-Cas9-based screens and computational approaches that help to further elucidate the emerging role of the non-coding genome in disease. Moreover, we discuss existing literature on mutations or copy number alterations of regulatory regions involved in brain development. We foresee that the future implementation of the knowledge obtained through ongoing gene regulation studies will benefit patients and will provide an explanation to part of the missing heritability of MCDs and other genetic disorders.

Integrated Functional Genomic Analysis Enables Annotation of Kidney Genome-Wide Association Study Loci

BACKGROUND

Linking genetic risk loci identified by genome-wide association studies (GWAS) to their causal genes remains a major challenge. Disease-associated genetic variants are concentrated in regions containing regulatory DNA elements, such as promoters and enhancers. Although researchers have previously published DNA maps of these regulatory regions for kidney tubule cells and glomerular endothelial cells, maps for podocytes and mesangial cells have not been available.

METHODS

We generated regulatory DNA maps (DNase-seq) and paired gene expression profiles (RNA-seq) from primary outgrowth cultures of human glomeruli that were composed mainly of podocytes and mesangial cells. We generated similar datasets from renal cortex cultures, to compare with those of the glomerular cultures. Because regulatory DNA elements can act on target genes across large genomic distances, we also generated a chromatin conformation map from freshly isolated human glomeruli.

RESULTS

We identified thousands of unique regulatory DNA elements, many located close to transcription factor genes, which the glomerular and cortex samples expressed at different levels. We found that genetic variants associated with kidney diseases (GWAS) and kidney expression quantitative trait loci were enriched in regulatory DNA regions. By combining GWAS, epigenomic, and chromatin conformation data, we functionally annotated 46 kidney disease genes.

CONCLUSIONS

We demonstrate a powerful approach to functionally connect kidney disease-/trait-associated loci to their target genes by leveraging unique regulatory DNA maps and integrated epigenomic and genetic analysis. This process can be applied to other kidney cell types and will enhance our understanding of genome regulation and its effects on gene expression in kidney disease.

Anchor: trans-cell type prediction of transcription factor binding sites

The ENCyclopedia of DNA Elements (ENCODE) consortium has generated transcription factor (TF) binding ChIP-seq data covering hundreds of TF proteins and cell types; however, due to limits on time and resources, only a small fraction of all possible TF-cell type pairs have been profiled. One solution is to build machine learning models trained on currently available epigenomic data sets that can be applied to the remaining missing pairs. A major challenge is that TF binding sites are cell-type-specific, which can be attributed to cellular contexts such as chromatin accessibility. Meanwhile, indirect TF-DNA binding and interactions between TFs complicate this regulatory process. Technical issues such as sequencing biases and batch effects render the prediction task even more challenging. Many pioneering efforts have been made to predict TF binding profiles based on DNA sequence and DNase-seq footprints, but to what extent a model can be generalized to completely untested cell conditions remains unknown. In this study, we describe our first place solution to the 2017 ENCODE-DREAM in vivo TF binding site prediction challenge. By carefully addressing multisource biases and information imbalance across cell types, we created a pipeline that significantly outperforms the current state-of-the-art methods. The proposed method is sufficiently complex enough to model nonlinear interactions between TF binding motifs and chromatin accessibility information up to 1500 bp from the genomic region of interest.

Mapping and Dynamics of Regulatory DNA in Maturing Arabidopsis thaliana Siliques

The genome is reprogrammed during development to produce diverse cell types, largely through altered expression and activity of key transcription factors. The accessibility and critical functions of epidermal cells have made them a model for connecting transcriptional events to development in a range of model systems. In Arabidopsis thaliana and many other plants, fertilization triggers differentiation of specialized epidermal seed coat cells that have a unique morphology caused by large extracellular deposits of polysaccharides. Here, we used DNase I-seq to generate regulatory landscapes of A. thaliana seeds at two critical time points in seed coat maturation (4 and 7 DPA), enriching for seed coat cells with the INTACT method. We found over 3,000 developmentally dynamic regulatory DNA elements and explored their relationship with nearby gene expression. The dynamic regulatory elements were enriched for motifs for several transcription factors families; most notably the TCP family at the earlier time point and the MYB family at the later one. To assess the extent to which the observed regulatory sites in seeds added to previously known regulatory sites in A. thaliana, we compared our data to 11 other data sets generated with 7-day-old seedlings for diverse tissues and conditions. Surprisingly, over a quarter of the regulatory, i.e. accessible, bases observed in seeds were novel. Notably, plant regulatory landscapes from different tissues, cell types, or developmental stages were more dynamic than those generated from bulk tissue in response to environmental perturbations, highlighting the importance of extending studies of regulatory DNA to single tissues and cell types during development.

Genome-wide analysis of chromatin accessibility using ATAC-seq

Programs of gene transcription are controlled by cis-acting DNA elements, including enhancers, silencers, and promoters. Local accessibility of chromatin has proven to be a highly informative structural feature for identifying such regulatory elements, which tend to be relatively open due to their interactions with proteins. Recently, ATAC-seq (assay for transposase-accessible chromatin using sequencing) has emerged as one of the most powerful approaches for genome-wide chromatin accessibility profiling. This method assesses DNA accessibility using hyperactive Tn5 transposase, which simultaneously cuts DNA and inserts sequencing adaptors, preferentially in regions of open chromatin. ATAC-seq is a relatively simple procedure which can be applied to only a few thousand cells. It is well-suited to developing embryos of sea urchins and other echinoderms, which are a prominent experimental model for understanding the genomic control of animal development. In this chapter, we present a protocol for applying ATAC-seq to embryonic cells of sea urchins.

Changes in chromatin accessibility between Arabidopsis stem cells and mesophyll cells illuminate cell type-specific transcription factor networks

Cell differentiation is driven by changes in the activity of transcription factors (TFs) and subsequent alterations in transcription. To study this process, differences in TF binding between cell types can be deduced by probing chromatin accessibility. We used cell type-specific nuclear purification followed by the assay for transposase-accessible chromatin (ATAC-seq) to delineate differences in chromatin accessibility and TF regulatory networks between stem cells of the shoot apical meristem (SAM) and differentiated leaf mesophyll cells in Arabidopsis thaliana. Chromatin accessibility profiles of SAM stem cells and leaf mesophyll cells were very similar at a qualitative level, yet thousands of regions having quantitatively different chromatin accessibility were also identified. Analysis of the genomic regions preferentially accessible in each cell type identified hundreds of overrepresented TF-binding motifs, highlighting sets of TFs that are probably important for each cell type. Within these sets, we found evidence for extensive co-regulation of target genes by multiple TFs that are preferentially expressed in each cell type. Interestingly, the TFs within each of these cell type-enriched sets also showed evidence of extensively co-regulating each other. We further found that preferentially accessible chromatin regions in mesophyll cells tended to also be substantially accessible in the stem cells, whereas the converse was not true. This observation suggests that the generally higher accessibility of regulatory elements in stem cells might contribute to their developmental plasticity. This work demonstrates the utility of cell type-specific chromatin accessibility profiling for the rapid development of testable models of regulatory control differences between cell types.

Global analysis of primary mesenchyme cell cis-regulatory modules by chromatin accessibility profiling

Background

The developmental gene regulatory network (GRN) that underlies skeletogenesis in sea urchins and other echinoderms is a paradigm of GRN structure, function, and evolution. This transcriptional network is deployed selectively in skeleton-forming primary mesenchyme cells (PMCs) of the early embryo. To advance our understanding of this model developmental GRN, we used genome-wide chromatin accessibility profiling to identify and characterize PMC cis-regulatory modules (CRMs).

Results

ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) analysis of purified PMCs provided a global picture of chromatin accessibility in these cells. We used both ATAC-seq and DNase-seq (DNase I hypersensitive site sequencing) to identify > 3000 sites that exhibited increased accessibility in PMCs relative to other embryonic cell lineages, and provide both computational and experimental evidence that a large fraction of these sites represent bona fide skeletogenic CRMs. Putative PMC CRMs were preferentially located near genes differentially expressed by PMCs and consensus binding sites for two key transcription factors in the PMC GRN, Alx1 and Ets1, were enriched in these CRMs. Moreover, a high proportion of candidate CRMs drove reporter gene expression specifically in PMCs in transgenic embryos. Surprisingly, we found that PMC CRMs were partially open in other embryonic lineages and exhibited hyperaccessibility as early as the 128-cell stage.

Conclusions

Our work provides a comprehensive picture of chromatin accessibility in an early embryonic cell lineage. By identifying thousands of candidate PMC CRMs, we significantly enhance the utility of the sea urchin skeletogenic network as a general model of GRN architecture and evolution. Our work also shows that differential chromatin accessibility, which has been used for the high-throughput identification of enhancers in differentiated cell types, is a powerful approach for the identification of CRMs in early embryonic cells. Lastly, we conclude that in the sea urchin embryo, CRMs that control the cell type-specific expression of effector genes are hyperaccessible several hours in advance of gene activation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4542-z) contains supplementary material, which is available to authorized users.

Rapid chromatin repression by Aire provides precise control of immune tolerance

Aire mediates the expression of tissue-specific antigens in thymic epithelial cells to promote tolerance against self-reactive T lymphocytes. However, the mechanism that allows expression of tissue-specific genes at levels that prevent harm is unknown. Here we show that Brg1 generates accessibility at tissue-specific loci to impose central tolerance. We found that Aire has an intrinsic repressive function that restricts chromatin accessibility and opposes Brg1 across the genome. Aire exerted this repressive influence within minutes after recruitment to chromatin and restrained the amplitude of active transcription. Disease-causing mutations that impair Aire-induced activation also impair the protein's repressive function, which indicates dual roles for Aire. Together, Brg1 and Aire fine-tune the expression of tissue-specific genes at levels that prevent toxicity yet promote immune tolerance.

Application of CRISPR/Cas to Understand Cis- and Trans-Regulatory Elements in Plants

The recent emergence of the CRISPR/Cas system as a genome editing tool enables simple, fast, and efficient induction of DNA double-strand breaks at precise positions in the genome. This has proven extremely useful for analysis and modification of protein-coding sequences. Regulatory sequences have received much less attention, but can now be quickly and easily disrupted as well. Editing of cis-regulatory elements (CRE) offers considerable potential for crop improvement via fine-tuning of gene expression that cannot be achieved by simple KO mutations, but its widespread application is still hampered by a lack of precise knowledge about functional motifs in CRE. As demonstrated for mammalian cells, CRISPR/Cas is also extremely useful for the identification and analysis of CRE in their native environment on a large scale using tiling screens. Transcriptional complexes are another promising target for crop genome editing, as demonstrated for pathogen resistance and regulation of flowering. The development of more diverse and sophisticated CRISPR/Cas tools for genome editing will allow even more efficient and powerful approaches for editing of regulatory sequences in the future.

Single-cell whole-genome analyses by Linear Amplification via Transposon Insertion (LIANTI)

Single-cell genomics is important for biology and medicine. However, current whole-genome amplification (WGA) methods are limited by low accuracy of copy-number variation (CNV) detection and low amplification fidelity. Here we report an improved single-cell WGA method, Linear Amplification via Transposon Insertion (LIANTI), which outperforms existing methods, enabling micro-CNV detection with kilobase resolution. This allowed direct observation of stochastic firing of DNA replication origins, which differs from cell to cell. We also show that the predominant cytosine-to-thymine mutations observed in single-cell genomics often arise from the artifact of cytosine deamination upon cell lysis. However, identifying single-nucleotide variations (SNVs) can be accomplished by sequencing kindred cells. We determined the spectrum of SNVs in a single human cell after ultraviolet radiation, revealing their nonrandom genome-wide distribution.

Transcriptional Dynamics at Brain Enhancers: from Functional Specialization to Neurodegeneration

Over the last decade, the noncoding part of the genome has been shown to harbour thousands of cis-regulatory elements, such as enhancers, that activate well-defined gene expression programs. Driven by the development of numerous techniques, many of these elements are now identified in multiple tissues and cell types, and their characteristics as well as importance in development and disease are becoming increasingly clear. Here, we provide an overview of the insights that were gained from the analysis of noncoding gene regulatory elements in the brain and describe their potential contribution to cell type specialization, brain function and neurodegenerative disease.

Inference of cell type specific regulatory networks on mammalian lineages

Transcriptional regulatory networks are at the core of establishing cell type specific gene expression programs. In mammalian systems, such regulatory networks are determined by multiple levels of regulation, including by transcription factors, chromatin environment, and three-dimensional organization of the genome. Recent efforts to measure diverse regulatory genomic datasets across multiple cell types and tissues offer unprecedented opportunities to examine the context-specificity and dynamics of regulatory networks at a greater resolution and scale than before. In parallel, numerous computational approaches to analyze these data have emerged that serve as important tools for understanding mammalian cell type specific regulation. In this article, we review recent computational approaches to predict the expression and sequence-based regulators of a gene's expression level and examine long-range gene regulation. We highlight promising approaches, insights gained, and open challenges that need to be overcome to build a comprehensive picture of cell type specific transcriptional regulatory networks.

The birth of a human-specific neural gene by incomplete duplication and gene fusion

BACKGROUND

Gene innovation by duplication is a fundamental evolutionary process but is difficult to study in humans due to the large size, high sequence identity, and mosaic nature of segmental duplication blocks. The human-specific gene hydrocephalus-inducing 2, HYDIN2, was generated by a 364 kbp duplication of 79 internal exons of the large ciliary gene HYDIN from chromosome 16q22.2 to chromosome 1q21.1. Because the HYDIN2 locus lacks the ancestral promoter and seven terminal exons of the progenitor gene, we sought to characterize transcription at this locus by coupling reverse transcription polymerase chain reaction and long-read sequencing.

RESULTS

5' RACE indicates a transcription start site for HYDIN2 outside of the duplication and we observe fusion transcripts spanning both the 5' and 3' breakpoints. We observe extensive splicing diversity leading to the formation of altered open reading frames (ORFs) that appear to be under relaxed selection. We show that HYDIN2 adopted a new promoter that drives an altered pattern of expression, with highest levels in neural tissues. We estimate that the HYDIN duplication occurred ~3.2 million years ago and find that it is nearly fixed (99.9%) for diploid copy number in contemporary humans. Examination of 73 chromosome 1q21 rearrangement patients reveals that HYDIN2 is deleted or duplicated in most cases.

CONCLUSIONS

Together, these data support a model of rapid gene innovation by fusion of incomplete segmental duplications, altered tissue expression, and potential subfunctionalization or neofunctionalization of HYDIN2 early in the evolution of the Homo lineage.

ALTRE: workflow for defining ALTered Regulatory Elements using chromatin accessibility data

Summary

Regulatory elements regulate gene transcription, and their location and accessibility is cell-type specific, particularly for enhancers. Mapping and comparing chromatin accessibility between different cell types may identify mechanisms involved in cellular development and disease progression. To streamline and simplify differential analysis of regulatory elements genome-wide using chromatin accessibility data, such as DNase-seq, ATAC-seq, we developed ALTRE ( ALT ered R egulatory E lements), an R package and associated R Shiny web app. ALTRE makes such analysis accessible to a wide range of users-from novice to practiced computational biologists.

Availability and Implementation

https://github.com/Mathelab/ALTRE.

Contact

ewy.mathe@osumc.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Deep homology in the age of next-generation sequencing

The principle of homology is central to conceptualizing the comparative aspects of morphological evolution. The distinctions between homologous or non-homologous structures have become blurred, however, as modern evolutionary developmental biology (evo-devo) has shown that novel features often result from modification of pre-existing developmental modules, rather than arising completely de novo. With this realization in mind, the term 'deep homology' was coined, in recognition of the remarkably conserved gene expression during the development of certain animal structures that would not be considered homologous by previous strict definitions. At its core, it can help to formulate an understanding of deeper layers of ontogenetic conservation for anatomical features that lack any clear phylogenetic continuity. Here, we review deep homology and related concepts in the context of a gene expression-based homology discussion. We then focus on how these conceptual frameworks have profited from the recent rise of high-throughput next-generation sequencing. These techniques have greatly expanded the range of organisms amenable to such studies. Moreover, they helped to elevate the traditional gene-by-gene comparison to a transcriptome-wide level. We will end with an outlook on the next challenges in the field and how technological advances might provide exciting new strategies to tackle these questions.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

The Cartography of UV-induced DNA Damage Formation and DNA Repair

DNA damage presents a barrier to DNA-templated biochemical processes, including gene expression and faithful DNA replication. Compromised DNA repair leads to mutations, enhancing the risk for genetic diseases and cancer development. Conventional experimental approaches to study DNA damage required a researcher to choose between measuring bulk damage over the entire genome, with little or no resolution regarding a specific location, and obtaining data specific to a locus of interest, without a global perspective. Recent advances in high-throughput genomic tools overcame these limitations and provide high-resolution measurements simultaneously across the genome. In this review, we discuss the available methods for measuring DNA damage and their repair, focusing on genomewide assays for pyrimidine photodimers, the major types of damage induced by ultraviolet irradiation. These new genomic assays will be a powerful tool in identifying key components of genome stability and carcinogenesis.

Identification of Regulatory DNA Elements Using Genome-wide Mapping of DNase I Hypersensitive Sites during Tomato Fruit Development

Development and ripening of tomato fruit are precisely controlled by transcriptional regulation, which depends on the orchestrated accessibility of regulatory proteins to promoters and other cis-regulatory DNA elements. This accessibility and its effect on gene expression play a major role in defining the developmental process. To understand the regulatory mechanism and functional elements modulating morphological and anatomical changes during fruit development, we generated genome-wide high-resolution maps of DNase I hypersensitive sites (DHSs) from the fruit tissues of the tomato cultivar "Moneymaker" at 20 days post anthesis as well as break stage. By exploring variation of DHSs across fruit development stages, we pinpointed the most likely hypersensitive sites related to development-specific genes. By detecting binding motifs on DHSs of these development-specific genes or genes in the ascorbic acid biosynthetic pathway, we revealed the common regulatory elements contributing to coordinating gene transcription of plant ripening and specialized metabolic pathways. Our results contribute to a better understanding of the regulatory dynamics of genes involved in tomato fruit development and ripening.

Enhancer regions show high histone H3.3 turnover that changes during differentiation

The organization of DNA into chromatin is dynamic; nucleosomes are frequently displaced to facilitate the ability of regulatory proteins to access specific DNA elements. To gain insight into nucleosome dynamics, and to follow how dynamics change during differentiation, we used a technique called time-ChIP to quantitatively assess histone H3.3 turnover genome-wide during differentiation of mouse ESCs. We found that, without prior assumptions, high turnover could be used to identify regions involved in gene regulation. High turnover was seen at enhancers, as observed previously, with particularly high turnover at super-enhancers. In contrast, regions associated with the repressive Polycomb-Group showed low turnover in ESCs. Turnover correlated with DNA accessibility. Upon differentiation, numerous changes in H3.3 turnover rates were observed, the majority of which occurred at enhancers. Thus, time-ChIP measurement of histone turnover shows that active enhancers are unusually dynamic in ESCs and changes in highly dynamic nucleosomes predominate at enhancers during differentiation.

Next generation sequencing technology and genomewide data analysis: Perspectives for retinal research

The advent of high throughput next generation sequencing (NGS) has accelerated the pace of discovery of disease-associated genetic variants and genomewide profiling of expressed sequences and epigenetic marks, thereby permitting systems-based analyses of ocular development and disease. Rapid evolution of NGS and associated methodologies presents significant challenges in acquisition, management, and analysis of large data sets and for extracting biologically or clinically relevant information. Here we illustrate the basic design of commonly used NGS-based methods, specifically whole exome sequencing, transcriptome, and epigenome profiling, and provide recommendations for data analyses. We briefly discuss systems biology approaches for integrating multiple data sets to elucidate gene regulatory or disease networks. While we provide examples from the retina, the NGS guidelines reviewed here are applicable to other tissues/cell types as well.

Decoding transcriptional enhancers: Evolving from annotation to functional interpretation

Deciphering the intricate molecular processes that orchestrate the spatial and temporal regulation of genes has become an increasingly major focus of biological research. The differential expression of genes by diverse cell types with a common genome is a hallmark of complex cellular functions, as well as the basis for multicellular life. Importantly, a more coherent understanding of gene regulation is critical for defining developmental processes, evolutionary principles and disease etiologies. Here we present our current understanding of gene regulation by focusing on the role of enhancer elements in these complex processes. Although functional genomic methods have provided considerable advances to our understanding of gene regulation, these assays, which are usually performed on a genome-wide scale, typically provide correlative observations that lack functional interpretation. Recent innovations in genome editing technologies have placed gene regulatory studies at an exciting crossroads, as systematic, functional evaluation of enhancers and other transcriptional regulatory elements can now be performed in a coordinated, high-throughput manner across the entire genome. This review provides insights on transcriptional enhancer function, their role in development and disease, and catalogues experimental tools commonly used to study these elements. Additionally, we discuss the crucial role of novel techniques in deciphering the complex gene regulatory landscape and how these studies will shape future research.

MNase titration reveals differences between nucleosome occupancy and chromatin accessibility

Chromatin accessibility plays a fundamental role in gene regulation. Nucleosome placement, usually measured by quantifying protection of DNA from enzymatic digestion, can regulate accessibility. We introduce a metric that uses micrococcal nuclease (MNase) digestion in a novel manner to measure chromatin accessibility by combining information from several digests of increasing depths. This metric, MACC (MNase accessibility), quantifies the inherent heterogeneity of nucleosome accessibility in which some nucleosomes are seen preferentially at high MNase and some at low MNase. MACC interrogates each genomic locus, measuring both nucleosome location and accessibility in the same assay. MACC can be performed either with or without a histone immunoprecipitation step, and thereby compares histone and non-histone protection. We find that changes in accessibility at enhancers, promoters and other regulatory regions do not correlate with changes in nucleosome occupancy. Moreover, high nucleosome occupancy does not necessarily preclude high accessibility, which reveals novel principles of chromatin regulation.

Assay for transposase-accessible chromatin and circularized chromosome conformation capture, two methods to explore the regulatory landscapes of genes in zebrafish

Accurate transcriptional control of genes is fundamental for the correct functioning of organs and developmental processes. This control depends on the interplay between the promoter of genes and other noncoding sequences, whose interaction is mediated by 3D chromatin arrangements. Thus, the detailed description of transcriptional regulatory landscapes is essential to understand the mechanisms of transcriptional regulation. However, to achieve that, two important challenges have to be faced: (1) the identification of the noncoding sequences that contribute to gene transcription and (2) the association of these sequences to the respective genes they control. In this chapter, we describe two protocols that allow overcoming these important challenges: the assay for transposase-accessible chromatin using sequencing (ATAC-seq) and circularized chromosome conformation capture (4C-seq). ATAC-seq is a very efficient technique that, using a very low number of cells as starting material, allows the identification of active chromatin regions genome wide, whereas 4C-seq detects the subset of sequences that interact specifically with the promoter of a given gene. When combined, both techniques provide a comprehensive snapshot of the regulatory landscapes of developmental genes. The protocols we present here have been optimized for teleost fish samples, zebrafish and medaka, allowing the in-depth study of transcriptional regulation in these two emerging animal models. Given the amenability and easy genetic manipulation of these two experimental systems, we anticipate that they will be important in revealing general principles of the vertebrate regulatory genome.

Genome-Wide Chromatin Landscape Transitions Identify Novel Pathways in Early Commitment to Osteoblast Differentiation

Bone continuously undergoes remodeling by a tightly regulated process that involves osteoblast differentiation from Mesenchymal Stem Cells (MSC). However, commitment of MSC to osteoblastic lineage is a poorly understood process. Chromatin organization functions as a molecular gatekeeper of DNA functions. Detection of sites that are hypersensitive to Dnase I has been used for detailed examination of changes in response to hormones and differentiation cues. To investigate the early steps in commitment of MSC to osteoblasts, we used a model human temperature-sensitive cell line, hFOB. When shifted to non-permissive temperature, these cells undergo "spontaneous" differentiation that takes several weeks, a process that is greatly accelerated by osteogenic induction media. We performed Dnase I hypersensitivity assays combined with deep sequencing to identify genome-wide potential regulatory events in cells undergoing early steps of commitment to osteoblasts. Massive reorganization of chromatin occurred within hours of differentiation. Whereas ~30% of unique DHS sites were located in the promoters, the majority was outside of the promoters, designated as enhancers. Many of them were at novel genomic sites and need to be confirmed experimentally. We developed a novel method for identification of cellular networks based solely on DHS enhancers signature correlated to gene expression. The analysis of enhancers that were unique to differentiating cells led to identification of bone developmental program encompassing 147 genes that directly or indirectly participate in osteogenesis. Identification of these pathways provided an unprecedented view of genomic regulation during early steps of differentiation and changes related to WNT, AP-1 and other pathways may have therapeutic implications.

Genome-Scale Analysis of Cell-Specific Regulatory Codes Using Nuclear Enzymes

High-throughput sequencing technologies have made it possible for biologists to generate genome-wide profiles of chromatin features at the nucleotide resolution. Enzymes such as nucleases or transposes have been instrumental as a chromatin-probing agent due to their ability to target accessible chromatin for cleavage or insertion. On the scale of a few hundred base pairs, preferential action of the nuclear enzymes on accessible chromatin allows mapping of cell state-specific accessibility in vivo. Such accessible regions contain functionally important regulatory sites, including promoters and enhancers, which undergo active remodeling for cells adapting in a dynamic environment. DNase-seq and the more recent ATAC-seq are two assays that are gaining popularity. Deep sequencing of DNA libraries from these assays, termed genomic footprinting, has been proposed to enable the comprehensive construction of protein occupancy profiles over the genome at the nucleotide level. Recent studies have discovered limitations of genomic footprinting which reduce the scope of detectable proteins. In addition, the identification of putative factors that bind to the observed footprints remains challenging. Despite these caveats, the methodology still presents significant advantages over alternative techniques such as ChIP-seq or FAIRE-seq. Here we describe computational approaches and tools for analysis of chromatin accessibility and genomic footprinting. Proper experimental design and assay-specific data analysis ensure the detection sensitivity and maximize retrievable information. The enzyme-based chromatin profiling approaches represent a powerful and evolving methodology which facilitates our understanding of how the genome is regulated.

Competition between DNA methylation and transcription factors determines binding of NRF1

Eukaryotic transcription factors (TFs) are key determinants of gene activity, yet they bind only a fraction of their corresponding DNA sequence motifs in any given cell type. Chromatin has the potential to restrict accessibility of binding sites; however, in which context chromatin states are instructive for TF binding remains mainly unknown. To explore the contribution of DNA methylation to constrained TF binding, we mapped DNase-I-hypersensitive sites in murine stem cells in the presence and absence of DNA methylation. Methylation-restricted sites are enriched for TF motifs containing CpGs, especially for those of NRF1. In fact, the TF NRF1 occupies several thousand additional sites in the unmethylated genome, resulting in increased transcription. Restoring de novo methyltransferase activity initiates remethylation at these sites and outcompetes NRF1 binding. This suggests that binding of DNA-methylation-sensitive TFs relies on additional determinants to induce local hypomethylation. In support of this model, removal of neighbouring motifs in cis or of a TF in trans causes local hypermethylation and subsequent loss of NRF1 binding. This competition between DNA methylation and TFs in vivo reveals a case of cooperativity between TFs that acts indirectly via DNA methylation. Methylation removal by methylation-insensitive factors enables occupancy of methylation-sensitive factors, a principle that rationalizes hypomethylation of regulatory regions.

Large-scale identification of sequence variants influencing human transcription factor occupancy in vivo

The function of human regulatory regions depends exquisitely on their local genomic environment and on cellular context, complicating experimental analysis of common disease- and trait-associated variants that localize within regulatory DNA. We use allelically resolved genomic DNase I footprinting data encompassing 166 individuals and 114 cell types to identify >60,000 common variants that directly influence transcription factor occupancy and regulatory DNA accessibility in vivo. The unprecedented scale of these data enables systematic analysis of the impact of sequence variation on transcription factor occupancy in vivo. We leverage this analysis to develop accurate models of variation affecting the recognition sites for diverse transcription factors and apply these models to discriminate nearly 500,000 common regulatory variants likely to affect transcription factor occupancy across the human genome. The approach and results provide a new foundation for the analysis and interpretation of noncoding variation in complete human genomes and for systems-level investigation of disease-associated variants.

Functional footprinting of regulatory DNA

Regulatory regions harbor multiple transcription factor (TF) recognition sites; however, the contribution of individual sites to regulatory function remains challenging to define. We describe an approach that exploits the error-prone nature of genome editing-induced double-strand break repair to map functional elements within regulatory DNA at nucleotide resolution. We demonstrate the approach on a human erythroid enhancer, revealing single TF recognition sites that gate the majority of downstream regulatory function.

Comparative analysis of gene regulatory networks: from network reconstruction to evolution

Regulation of gene expression is central to many biological processes. Although reconstruction of regulatory circuits from genomic data alone is therefore desirable, this remains a major computational challenge. Comparative approaches that examine the conservation and divergence of circuits and their components across strains and species can help reconstruct circuits as well as provide insights into the evolution of gene regulatory processes and their adaptive contribution. In recent years, advances in genomic and computational tools have led to a wealth of methods for such analysis at the sequence, expression, pathway, module, and entire network level. Here, we review computational methods developed to study transcriptional regulatory networks using comparative genomics, from sequence to functional data. We highlight how these methods use evolutionary conservation and divergence to reliably detect regulatory components as well as estimate the extent and rate of divergence. Finally, we discuss the promise and open challenges in linking regulatory divergence to phenotypic divergence and adaptation.

Cardiac gene expression data and in silico analysis provide novel insights into human and mouse taste receptor gene regulation

G protein-coupled receptors are the principal mediators of the sweet, umami, bitter, and fat taste qualities in mammals. Intriguingly, the taste receptors are also expressed outside of the oral cavity, including in the gut, airways, brain, and heart, where they have additional functions and contribute to disease. However, there is little known about the mechanisms governing the transcriptional regulation of taste receptor genes. Following our recent delineation of taste receptors in the heart, we investigated the genomic loci encoding for taste receptors to gain insight into the regulatory mechanisms that drive their expression in the heart. Gene expression analyses of healthy and diseased human and mouse hearts showed coordinated expression for a subset of chromosomally clustered taste receptors. This chromosomal clustering mirrored the cardiac expression profile, suggesting that a common gene regulatory block may control the taste receptor locus. We identified unique domains with strong regulatory potential in the vicinity of taste receptor genes. We also performed de novo motif enrichment in the proximal promoter regions and found several overrepresented DNA motifs in cardiac taste receptor gene promoters corresponding to ubiquitous and cardiac-specific transcription factor binding sites. Thus, combining cardiac gene expression data with bioinformatic analyses, this study has provided insights into the noncoding regulatory landscape for taste GPCRs. These findings also have broader relevance for the study of taste GPCRs outside of the classical gustatory system, where understanding the mechanisms controlling the expression of these receptors may have implications for future therapeutic development.

DNase I hypersensitivity analysis of the mouse brain and retina identifies region-specific regulatory elements

BACKGROUND

The brain, spinal cord, and neural retina comprise the central nervous system (CNS) of vertebrates. Understanding the regulatory mechanisms that underlie the enormous cell-type diversity of the CNS is a significant challenge. Whole-genome mapping of DNase I-hypersensitive sites (DHSs) has been used to identify cis-regulatory elements in many tissues. We have applied this approach to the mouse CNS, including developing and mature neural retina, whole brain, and two well-characterized brain regions, the cerebellum and the cerebral cortex.

RESULTS

For the various regions and developmental stages of the CNS that we analyzed, there were approximately the same number of DHSs; however, there were many DHSs unique to each CNS region and developmental stage. Many of the DHSs are likely to mark enhancers that are specific to the specific CNS region and developmental stage. We validated the DNase I mapping approach for identification of CNS enhancers using the existing VISTA Browser database and with in vivo and in vitro electroporation of the retina. Analysis of transcription factor consensus sites within the DHSs shows distinct region-specific profiles of transcriptional regulators particular to each region. Clustering developmentally dynamic DHSs in the retina revealed enrichment of developmental stage-specific transcriptional regulators. Additionally, we found reporter gene activity in the retina driven from several previously uncharacterized regulatory elements surrounding the neurodevelopmental gene Otx2. Identification of DHSs shared between mouse and human showed region-specific differences in the evolution of cis-regulatory elements.

CONCLUSIONS

Overall, our results demonstrate the potential of genome-wide DNase I mapping to cis-regulatory questions regarding the regional diversity within the CNS. These data represent an extensive catalogue of potential cis-regulatory elements within the CNS that display region and temporal specificity, as well as a set of DHSs common to CNS tissues. Further examination of evolutionary conservation of DHSs between CNS regions and different species may reveal important cis-regulatory elements in the evolution of the mammalian CNS.

Integrative analysis of 111 reference human epigenomes

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

CENP-A nucleosomes localize to transcription factor hotspots and subtelomeric sites in human cancer cells

Background

The histone H3 variant CENP-A is normally tightly regulated to ensure only one centromere exists per chromosome. Native CENP-A is often found overexpressed in human cancer cells and a range of human tumors. Consequently, CENP-A misregulation is thought to contribute to genome instability in human cancers. However, the consequences of such overexpression have not been directly elucidated in human cancer cells.

Results

To investigate native CENP-A overexpression, we sought to uncover CENP-A-associated defects in human cells. We confirm that CENP-A is innately overexpressed in several colorectal cancer cell lines. In such cells, we report that a subset of structurally distinct CENP-A-containing nucleosomes associate with canonical histone H3, and with the transcription-coupled chaperones ATRX and DAXX. Furthermore, such hybrid CENP-A nucleosomes localize to DNase I hypersensitive and transcription factor binding sites, including at promoters of genes across the human genome. A distinct class of CENP-A hotspots also accumulates at subtelomeric chromosomal locations, including at the 8q24/Myc region long-associated with genomic instability. We show this 8q24 accumulation of CENP-A can also be seen in early stage primary colorectal tumors.

Conclusions

Our data demonstrate that excess CENP-A accumulates at noncentromeric locations in the human cancer genome. These findings suggest that ectopic CENP-A nucleosomes could alter the state of the chromatin fiber, potentially impacting gene regulation and chromosome fragility.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-8935-8-2) contains supplementary material, which is available to authorized users.

Improved DNase-seq protocol facilitates high resolution mapping of DNase I hypersensitive sites in roots in Arabidopsis thaliana

BACKGROUND

Identifying cis-regulatory elements is critical in understanding the direct and indirect regulatory mechanisms of gene expression. Current approaches include DNase-seq, a technique that combines sensitivity to the nonspecific endonuclease DNase I with high throughput sequencing to identify regions of regulatory DNA on a genome-wide scale. While this method was originally developed for human cell lines, later adaptations made the processing of plant tissues possible. Challenges still remain in processing recalcitrant tissues that have low DNA content.

RESULTS

By removing steps requiring the use of gel agarose plugs in DNase-seq, we were able to significantly reduce the time required to perform the protocol by at least 2 days, while also making possible the processing of difficult plant tissues. We refer to this simplified protocol as DNase I SIM (for simplified in-nucleus method). We were able to successfully create DNase-seq libraries for both leaf and root tissues in Arabidopsis using DNase I SIM.

CONCLUSION

This protocol simplifies and facilitates generation of DNase-seq libraries from plant tissues for high resolution mapping of DNase I hypersensitive sites.

Mouse regulatory DNA landscapes reveal global principles of cis-regulatory evolution

To study the evolutionary dynamics of regulatory DNA, we mapped >1.3 million deoxyribonuclease I-hypersensitive sites (DHSs) in 45 mouse cell and tissue types, and systematically compared these with human DHS maps from orthologous compartments. We found that the mouse and human genomes have undergone extensive cis-regulatory rewiring that combines branch-specific evolutionary innovation and loss with widespread repurposing of conserved DHSs to alternative cell fates, and that this process is mediated by turnover of transcription factor (TF) recognition elements. Despite pervasive evolutionary remodeling of the location and content of individual cis-regulatory regions, within orthologous mouse and human cell types the global fraction of regulatory DNA bases encoding recognition sites for each TF has been strictly conserved. Our findings provide new insights into the evolutionary forces shaping mammalian regulatory DNA landscapes.

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.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

Assaying the epigenome in limited numbers of cells

Spectacular advances in the throughput of DNA sequencing have allowed genome-wide analysis of epigenetic features such as methylation, nucleosome position and post-translational modification, chromatin accessibility and connectivity, and transcription factor binding. However, for rare or precious biological samples, input requirements of many of these methods limit their application. In this review we discuss recent advances for low-input genome-wide analysis of chromatin immunoprecipitation, methylation, DNA accessibility, and chromatin conformation.