Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position

The histone chaperone CAF-1 safeguards somatic cell identity

Cellular differentiation involves profound remodeling of chromatic landscapes, yet the mechanisms by which somatic cell identity is subsequently maintained remain incompletely understood. To further elucidate regulatory pathways that safeguard the somatic state, we performed two comprehensive RNAi screens targeting chromatin factors during transcription factor-mediated reprogramming of mouse fibroblasts to induced pluripotent stem cells (iPSCs). Remarkably, subunits of the chromatin assembly factor-1 (CAF-1) complex emerged as the most prominent hits from both screens, followed by modulators of lysine sumoylation and heterochromatin maintenance. Optimal modulation of both CAF-1 and transcription factor levels increased reprogramming efficiency by several orders of magnitude and facilitated iPSC formation in as little as 4 days. Mechanistically, CAF-1 suppression led to a more accessible chromatin structure at enhancer elements early during reprogramming. These changes were accompanied by a decrease in somatic heterochromatin domains, increased binding of Sox2 to pluripotency-specific targets and activation of associated genes. Notably, suppression of CAF-1 also enhanced the direct conversion of B cells into macrophages and fibroblasts into neurons. Together, our findings reveal the histone chaperone CAF-1 as a novel regulator of somatic cell identity during transcription factor-induced cell fate transitions and provide a potential strategy to modulate cellular plasticity in a regenerative setting.

Genome-Wide Epigenetic Studies in Human Disease: A Primer on -Omic Technologies

Epigenetic information encoded in covalent modifications of DNA and histone proteins regulates fundamental biological processes through the action of chromatin regulators, transcription factors, and noncoding RNA species. Epigenetic plasticity enables an organism to respond to developmental and environmental signals without genetic changes. However, aberrant epigenetic control plays a key role in pathogenesis of disease. Normal epigenetic states could be disrupted by detrimental mutations and expression alteration of chromatin regulators or by environmental factors. In this primer, we briefly review the epigenetic basis of human disease and discuss how recent discoveries in this field could be translated into clinical diagnosis, prevention, and treatment. We introduce platforms for mapping genome-wide chromatin accessibility, nucleosome occupancy, DNA-binding proteins, and DNA methylation, primarily focusing on the integration of DNA methylation and chromatin immunoprecipitation-sequencing technologies into disease association studies. We highlight practical considerations in applying high-throughput epigenetic assays and formulating analytical strategies. Finally, we summarize current challenges in sample acquisition, experimental procedures, data analysis, and interpretation and make recommendations on further refinement in these areas. Incorporating epigenomic testing into the clinical research arsenal will greatly facilitate our understanding of the epigenetic basis of disease and help identify novel therapeutic targets.

Mapping nucleosome positions using DNase-seq

Although deoxyribonuclease I (DNase I) was used to probe the structure of the nucleosome in the 1960s and 1970s, in the current high-throughput sequencing era, DNase I has mainly been used to study genomic regions devoid of nucleosomes. Here, we reveal for the first time that DNase I can be used to precisely map the (translational) positions of in vivo nucleosomes genome-wide. Specifically, exploiting a distinctive DNase I cleavage profile within nucleosome-associated DNA—including a signature 10.3 base pair oscillation that corresponds to accessibility of the minor groove as DNA winds around the nucleosome—we develop a Bayes-factor–based method that can be used to map nucleosome positions along the genome. Compared to methods that require genetically modified histones, our DNase-based approach is easily applied in any organism, which we demonstrate by producing maps in yeast and human. Compared to micrococcal nuclease (MNase)-based methods that map nucleosomes based on cuts in linker regions, we utilize DNase I cuts both outside and within nucleosomal DNA; the oscillatory nature of the DNase I cleavage profile within nucleosomal DNA enables us to identify translational positioning details not apparent in MNase digestion of linker DNA. Because the oscillatory pattern corresponds to nucleosome rotational positioning, it also reveals the rotational context of transcription factor (TF) binding sites. We show that potential binding sites within nucleosome-associated DNA are often centered preferentially on an exposed major or minor groove. This preferential localization may modulate TF interaction with nucleosome-associated DNA as TFs search for binding sites.

Integration of ATAC-seq and RNA-seq identifies human alpha cell and beta cell signature genes

Objective

Although glucagon-secreting α-cells and insulin-secreting β-cells have opposing functions in regulating plasma glucose levels, the two cell types share a common developmental origin and exhibit overlapping transcriptomes and epigenomes. Notably, destruction of β-cells can stimulate repopulation via transdifferentiation of α-cells, at least in mice, suggesting plasticity between these cell fates. Furthermore, dysfunction of both α- and β-cells contributes to the pathophysiology of type 1 and type 2 diabetes, and β-cell de-differentiation has been proposed to contribute to type 2 diabetes. Our objective was to delineate the molecular properties that maintain islet cell type specification yet allow for cellular plasticity. We hypothesized that correlating cell type-specific transcriptomes with an atlas of open chromatin will identify novel genes and transcriptional regulatory elements such as enhancers involved in α- and β-cell specification and plasticity.

Methods

We sorted human α- and β-cells and performed the “Assay for Transposase-Accessible Chromatin with high throughput sequencing” (ATAC-seq) and mRNA-seq, followed by integrative analysis to identify cell type-selective gene regulatory regions.

Results

We identified numerous transcripts with either α-cell- or β-cell-selective expression and discovered the cell type-selective open chromatin regions that correlate with these gene activation patterns. We confirmed cell type-selective expression on the protein level for two of the top hits from our screen. The “group specific protein” (GC; or vitamin D binding protein) was restricted to α-cells, while CHODL (chondrolectin) immunoreactivity was only present in β-cells. Furthermore, α-cell- and β-cell-selective ATAC-seq peaks were identified to overlap with known binding sites for islet transcription factors, as well as with single nucleotide polymorphisms (SNPs) previously identified as risk loci for type 2 diabetes.

Conclusions

We have determined the genetic landscape of human α- and β-cells based on chromatin accessibility and transcript levels, which allowed for detection of novel α- and β-cell signature genes not previously known to be expressed in islets. Using fine-mapping of open chromatin, we have identified thousands of potential cis-regulatory elements that operate in an endocrine cell type-specific fashion.

•

Defined open chromatin regions in sorted human α- and β-cells using ATAC-seq.

•

Detected type 2 diabetes-associated risk loci in human α- and β-cell open chromatin.

•

Classified human α- and β-cell-specific transcripts using mRNA-seq.

•

Discovered novel human α- and β-cell signature proteins.

•

Identified potential gene regulatory regions by integrating ATAC- and mRNA-seq data.

ANISEED 2015: a digital framework for the comparative developmental biology of ascidians

Ascidians belong to the tunicates, the sister group of vertebrates and are recognized model organisms in the field of embryonic development, regeneration and stem cells. ANISEED is the main information system in the field of ascidian developmental biology. This article reports the development of the system since its initial publication in 2010. Over the past five years, we refactored the system from an initial custom schema to an extended version of the Chado schema and redesigned all user and back end interfaces. This new architecture was used to improve and enrich the description of Ciona intestinalis embryonic development, based on an improved genome assembly and gene model set, refined functional gene annotation, and anatomical ontologies, and a new collection of full ORF cDNAs. The genomes of nine ascidian species have been sequenced since the release of the C. intestinalis genome. In ANISEED 2015, all nine new ascidian species can be explored via dedicated genome browsers, and searched by Blast. In addition, ANISEED provides full functional gene annotation, anatomical ontologies and some gene expression data for the six species with highest quality genomes. ANISEED is publicly available at: http://www.aniseed.cnrs.fr.

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.

A Genome-Wide Perspective on Metabolism

Mammals have at least 210 histologically diverse cell types (Alberts, Molecular biology of the cell. Garland Science, New York, 2008) and the number would be even higher if functional differences are taken into account. The genome in each of these cell types is differentially programmed to express the specific set of genes needed to fulfill the phenotypical requirements of the cell. Furthermore, in each of these cell types, the gene program can be differentially modulated by exposure to external signals such as hormones or nutrients. The basis for the distinct gene programs relies on cell type-selective activation of transcriptional enhancers, which in turn are particularly sensitive to modulation. Until recently we had only fragmented insight into the regulation of a few of these enhancers; however, the recent advances in high-throughput sequencing technologies have enabled the development of a large number of technologies that can be used to obtain genome-wide insight into how genomes are reprogrammed during development and in response to specific external signals. By applying such technologies, we have begun to reveal the cross-talk between metabolism and the genome, i.e., how genomes are reprogrammed in response to metabolites, and how the regulation of metabolic networks is coordinated at the genomic level.

Determining Nucleosome Position at Individual Loci After Biotic Stress Using MNase-qPCR

Nucleosome occupancy in promoter and genic regions can severely influence the transcription levels. Few methods have been established to investigate the nucleosome occupancy along the DNA. In this chapter we describe a detailed protocol to analyze the nucleosome occupancy at a specific locus using MNase-pPCR.

Multiplexed analysis of chromosome conformation at vastly improved sensitivity

Methods for analyzing chromosome conformation in mammalian cells are either low resolution or low throughput and are technically challenging. In next-generation (NG) Capture-C, we have redesigned the Capture-C method to achieve unprecedented levels of sensitivity and reproducibility. NG Capture-C can be used to analyze many genetic loci and samples simultaneously. High-resolution data can be produced with as few as 100,000 cells, and single-nucleotide polymorphisms can be used to generate allele-specific tracks. The method is straightforward to perform and should greatly facilitate the investigation of many questions related to gene regulation as well as the functional dissection of traits examined in genome-wide association studies.

Making sense of GWAS: using epigenomics and genome engineering to understand the functional relevance of SNPs in non-coding regions of the human genome

Considerable progress towards an understanding of complex diseases has been made in recent years due to the development of high-throughput genotyping technologies. Using microarrays that contain millions of single-nucleotide polymorphisms (SNPs), Genome Wide Association Studies (GWASs) have identified SNPs that are associated with many complex diseases or traits. For example, as of February 2015, 2111 association studies have identified 15,396 SNPs for various diseases and traits, with the number of identified SNP-disease/trait associations increasing rapidly in recent years. However, it has been difficult for researchers to understand disease risk from GWAS results. This is because most GWAS-identified SNPs are located in non-coding regions of the genome. It is important to consider that the GWAS-identified SNPs serve only as representatives for all SNPs in the same haplotype block, and it is equally likely that other SNPs in high linkage disequilibrium (LD) with the array-identified SNPs are causal for the disease. Because it was hoped that disease-associated coding variants would be identified if the true casual SNPs were known, investigators have expanded their analyses using LD calculation and fine-mapping. However, such analyses also identified risk-associated SNPs located in non-coding regions. Thus, the GWAS field has been left with the conundrum as to how a single-nucleotide change in a non-coding region could confer increased risk for a specific disease. One possible answer to this puzzle is that the variant SNPs cause changes in gene expression levels rather than causing changes in protein function. This review provides a description of (1) advances in genomic and epigenomic approaches that incorporate functional annotation of regulatory elements to prioritize the disease risk-associated SNPs that are located in non-coding regions of the genome for follow-up studies, (2) various computational tools that aid in identifying gene expression changes caused by the non-coding disease-associated SNPs, and (3) experimental approaches to identify target genes of, and study the biological phenotypes conferred by, non-coding disease-associated SNPs.

Regulatory T cell memory

Memory for antigen is a defining feature of adaptive immunity. Antigen-specific lymphocyte populations show an increase in number and function after antigen encounter and more rapidly re-expand upon subsequent antigen exposure. Studies of immune memory have primarily focused on effector B cells and T cells with microbial specificity, using prime-challenge models of infection. However, recent work has also identified persistently expanded populations of antigen-specific regulatory T cells that protect against aberrant immune responses. In this Review, we consider the parallels between memory effector T cells and memory regulatory T cells, along with the functional implications of regulatory memory in autoimmunity, antimicrobial host defence and maternal-fetal tolerance. In addition, we discuss emerging evidence for regulatory T cell memory in humans and key unanswered questions in this rapidly evolving field.

cChIP-seq: a robust small-scale method for investigation of histone modifications

BACKGROUND

ChIP-seq is highly utilized for mapping histone modifications that are informative about gene regulation and genome annotations. For example, applying ChIP-seq to histone modifications such as H3K4me1 has facilitated generating epigenomic maps of putative enhancers. This powerful technology, however, is limited in its application by the large number of cells required. ChIP-seq involves extensive manipulation of sample material and multiple reactions with limited quality control at each step, therefore, scaling down the number of cells required has proven challenging. Recently, several methods have been proposed to overcome this limit but most of these methods require extensive optimization to tailor the protocol to the specific antibody used or number of cells being profiled.

RESULTS

Here we describe a robust, yet facile method, which we named carrier ChIP-seq (cChIP-seq), for use on limited cell amounts. cChIP-seq employs a DNA-free histone carrier in order to maintain the working ChIP reaction scale, removing the need to tailor reactions to specific amounts of cells or histone modifications to be assayed. We have applied our method to three different histone modifications, H3K4me3, H3K4me1 and H3K27me3 in the K562 cell line, and H3K4me1 in H1 hESCs. We successfully obtained epigenomic maps for these histone modifications starting with as few as 10,000 cells. We compared cChIP-seq data to data generated as part of the ENCODE project. ENCODE data are the reference standard in the field and have been generated starting from tens of million of cells. Our results show that cChIP-seq successfully recapitulates bulk data. Furthermore, we showed that the differences observed between small-scale ChIP-seq data and ENCODE data are largely to be due to lab-to-lab variability rather than operating on a reduced scale.

CONCLUSIONS

Data generated using cChIP-seq are equivalent to reference epigenomic maps from three orders of magnitude more cells. Our method offers a robust and straightforward approach to scale down ChIP-seq to as low as 10,000 cells. The underlying principle of our strategy makes it suitable for being applied to a vast range of chromatin modifications without requiring expensive optimization. Furthermore, our strategy of a DNA-free carrier can be adapted to most ChIP-seq protocols.

A novel ATAC-seq approach reveals lineage-specific reinforcement of the open chromatin landscape via cooperation between BAF and p63

Background

Open chromatin regions are correlated with active regulatory elements in development and are dysregulated in diseases. The BAF (SWI/SNF) complex is essential for development, and has been demonstrated to remodel reconstituted chromatin in vitro and to control the accessibility of a few individual regions in vivo. However, it remains unclear where and how BAF controls the open chromatin landscape to regulate developmental processes, such as human epidermal differentiation.

Results

Using a novel “on-plate” ATAC-sequencing approach for profiling open chromatin landscapes with a low number of adherent cells, we demonstrate that the BAF complex is essential for maintaining 11.6 % of open chromatin regions in epidermal differentiation. These BAF-dependent open chromatin regions are highly cell-type-specific and are strongly enriched for binding sites for p63, a master epidermal transcription factor. The DNA sequences of p63 binding sites intrinsically favor nucleosome formation and are inaccessible in other cell types without p63 to prevent ectopic activation. In epidermal cells, BAF and p63 mutually recruit each other to maintain 14,853 open chromatin regions. We further demonstrate that BAF and p63 cooperatively position nucleosomes away from p63 binding sites and recruit transcriptional machinery to control tissue differentiation.

Conclusions

BAF displays high specificity in controlling the open chromatin landscape during epidermal differentiation by cooperating with the master transcription factor p63 to maintain lineage-specific open chromatin regions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0840-9) contains supplementary material, which is available to authorized users.

Glioma Stem-like Cells Keep Their H3.3 Variant Levels at Bay

Pediatric glioblastomas (GBMs) commonly harbor mutations in histone variant H3.3, while adult GBMs do not. In this issue of Cancer Cell, Gallo and colleagues demonstrate that adult GBM stem-like cells repress H3.3 expression to maintain self-renewal properties.

Fine-mapping cellular QTLs with RASQUAL and ATAC-seq

When cellular traits are measured using high-throughput DNA sequencing, quantitative trait loci (QTLs) manifest as fragment count differences between individuals and allelic differences within individuals. We present RASQUAL (Robust Allele-Specific Quantitation and Quality Control), a new statistical approach for association mapping that models genetic effects and accounts for biases in sequencing data using a single, probabilistic framework. RASQUAL substantially improves fine-mapping accuracy and sensitivity relative to existing methods in RNA-seq, DNase-seq and ChIP-seq data. We illustrate how RASQUAL can be used to maximize association detection by generating the first map of chromatin accessibility QTLs (caQTLs) in a European population using ATAC-seq. Despite a modest sample size, we identified 2,707 independent caQTLs (at a false discovery rate of 10%) and demonstrated how RASQUAL and ATAC-seq can provide powerful information for fine-mapping gene-regulatory variants and for linking distal regulatory elements with gene promoters. Our results highlight how combining between-individual and allele-specific genetic signals improves the functional interpretation of noncoding variation.

MLL5 Orchestrates a Cancer Self-Renewal State by Repressing the Histone Variant H3.3 and Globally Reorganizing Chromatin

Mutations in the histone 3 variant H3.3 have been identified in one-third of pediatric glioblastomas (GBMs), but not in adult tumors. Here we show that H3.3 is a dynamic determinant of functional properties in adult GBM. H3.3 is repressed by mixed lineage leukemia 5 (MLL5) in self-renewing GBM cells. MLL5 is a global epigenetic repressor that orchestrates reorganization of chromatin structure by punctuating chromosomes with foci of compacted chromatin, favoring tumorigenic and self-renewing properties. Conversely, H3.3 antagonizes self-renewal and promotes differentiation. We exploited these epigenetic states to rationally identify two small molecules that effectively curb cancer stem cell properties in a preclinical model. Our work uncovers a role for MLL5 and H3.3 in maintaining self-renewal hierarchies in adult GBM.

Diversity and Divergence of Dinoflagellate Histone Proteins

Histone proteins and the nucleosomal organization of chromatin are near-universal eukaroytic features, with the exception of dinoflagellates. Previous studies have suggested that histones do not play a major role in the packaging of dinoflagellate genomes, although several genomic and transcriptomic surveys have detected a full set of core histone genes. Here, transcriptomic and genomic sequence data from multiple dinoflagellate lineages are analyzed, and the diversity of histone proteins and their variants characterized, with particular focus on their potential post-translational modifications and the conservation of the histone code. In addition, the set of putative epigenetic mark readers and writers, chromatin remodelers and histone chaperones are examined. Dinoflagellates clearly express the most derived set of histones among all autonomous eukaryote nuclei, consistent with a combination of relaxation of sequence constraints imposed by the histone code and the presence of numerous specialized histone variants. The histone code itself appears to have diverged significantly in some of its components, yet others are conserved, implying conservation of the associated biochemical processes. Specifically, and with major implications for the function of histones in dinoflagellates, the results presented here strongly suggest that transcription through nucleosomal arrays happens in dinoflagellates. Finally, the plausible roles of histones in dinoflagellate nuclei are discussed.

Favorable genomic environments for cis-regulatory evolution: A novel theoretical framework

Cis-regulatory changes are arguably the primary evolutionary source of animal morphological diversity. With the recent explosion of genome-wide comparisons of the cis-regulatory content in different animal species is now possible to infer general principles underlying enhancer evolution. However, these studies have also revealed numerous discrepancies and paradoxes, suggesting that the mechanistic causes and modes of cis-regulatory evolution are still not well understood and are probably much more complex than generally appreciated. Here, we argue that the mutational mechanisms and genomic regions generating new regulatory activities must comply with the constraints imposed by the molecular properties of cis-regulatory elements (CREs) and the organizational features of long-range chromatin interactions. Accordingly, we propose a new integrative evolutionary framework for cis-regulatory evolution based on two major premises for the origin of novel enhancer activity: (i) an accessible chromatin environment and (ii) compatibility with the 3D structure and interactions of pre-existing CREs. Mechanisms and DNA sequences not fulfilling these premises, will be less likely to have a measurable impact on gene expression and as such, will have a minor contribution to the evolution of gene regulation. Finally, we discuss current comparative cis-regulatory data under the light of this new evolutionary model, and propose that the two most prominent mechanisms for the evolution of cis-regulatory changes are the overprinting of ancestral CREs and the exaptation of transposable elements.

Leukemia-Associated Cohesin Mutants Dominantly Enforce Stem Cell Programs and Impair Human Hematopoietic Progenitor Differentiation

Recurrent mutations in cohesin complex proteins have been identified in pre-leukemic hematopoietic stem cells and during the early development of acute myeloid leukemia and other myeloid malignancies. Although cohesins are involved in chromosome separation and DNA damage repair, cohesin complex functions during hematopoiesis and leukemic development are unclear. Here, we show that mutant cohesin proteins block differentiation of human hematopoietic stem and progenitor cells (HSPCs) in vitro and in vivo and enforce stem cell programs. These effects are restricted to immature HSPC populations, where cohesin mutants show increased chromatin accessibility and likelihood of transcription factor binding site occupancy by HSPC regulators including ERG, GATA2, and RUNX1, as measured by ATAC-seq and ChIP-seq. Epistasis experiments show that silencing these transcription factors rescues the differentiation block caused by cohesin mutants. Together, these results show that mutant cohesins impair HSPC differentiation by controlling chromatin accessibility and transcription factor activity, possibly contributing to leukemic disease.

Control of plant cell differentiation by histone modification and DNA methylation

How cells differentiate and acquire diverse arrays of determined states in multicellular organisms is a fundamental and yet unanswered question in biology. Molecular genetic studies over the last few decades have identified many transcriptional regulators that activate or repress gene expression to promote cell differentiation in plant development. What has recently emerged as an additional important regulatory layer is the control at the epigenetic level by which locus-specific DNA methylation and histone modification alter the chromatin state and limit the expression of key developmental regulators to specific windows of time and space. Accumulating evidence suggests that histone acetylation is commonly linked with active transcription and this mechanism is adopted to control sequential progression of cell differentiation. Histone H3 trimethylation at lysine 27 and DNA methylation are both associated with gene repression, and these mechanisms are often utilised to promote and/or maintain the differentiated status of plant cells.

Intrauterine growth restriction perturbs nucleosome depletion at a growth hormone-responsive element in the mouse IGF-1 gene

Intrauterine growth restriction (IUGR) is a common human pregnancy complication. IUGR offspring carry significant postnatal risk for early-onset metabolic syndrome, which is associated with persistent reduction in IGF-1 protein expression. We have previously shown that preadolescent IUGR male mice have decreased hepatic IGF-1 mRNA and circulating IGF-1 protein at postnatal day 21, the age when growth hormone (GH) normally upregulates hepatic IGF-1 expression. Here we studied nucleosome occupancy and CpG methylation at a putative growth hormone-responsive element in intron 2 (in2GHRE) of the hepatic IGF-1 gene in normal, sham-operated, and IUGR mice. Nucleosome occupancy and CpG methylation were determined in embryonic stem cells (ESCs) and in liver at postnatal days 14, 21, and 42. For CpG methylation, additional time points out to 2 yr were analyzed. We confirmed the putative mouse in2GHRE was GH-responsive, and in normal mice, a single nucleosome was displaced from the hepatic in2GHRE by postnatal day 21, which exposed two STAT5b DNA binding sites. Nucleosome displacement correlated with developmentally programmed CpG demethylation. Finally, IUGR significantly altered the nucleosome-depleted region (NDR) at the in2GHRE of IGF-1 on postnatal day 21, with either complete absence of the NDR or with a shifted NDR exposing only one of two STAT5b DNA binding sites. An NDR shift was also seen in offspring of sham-operated mothers. We conclude that prenatal insult such as IUGR or anesthesia/surgery could perturb the proper formation of a well-positioned NDR at the mouse hepatic IGF-1 in2GHRE necessary for transitioning to an open chromatin state.

The Genotype and Phenotype (GaP) registry: a living biobank for the analysis of quantitative traits

We describe the development of the Genotype and Phenotype (GaP) Registry, a living biobank of normal volunteers who are genotyped for genetic markers related to human disease. Participants in the GaP can be recalled for hypothesis driven study of disease associated genetic variants. The GaP has facilitated functional studies of several autoimmune disease associated loci including Csk, Blk, PDRM1 (Blimp-1) and PTPN22. It is likely that expansion of such living biobank registries will play an important role in studying and understanding the function of disease associated alleles in complex disease.

Architectural and Functional Commonalities between Enhancers and Promoters

With the explosion of genome-wide studies of regulated transcription, it has become clear that traditional definitions of enhancers and promoters need to be revisited. These control elements can now be characterized in terms of their local and regional architecture, their regulatory components, including histone modifications and associated binding factors, and their functional contribution to transcription. This Review discusses unifying themes between promoters and enhancers in transcriptional regulatory mechanisms.

Transcriptional Enhancers: Bridging the Genome and Phenome

Enhancers play a major role in animal development by modulating spatiotemporal expression of genes. They interact with sequence-specific transcriptional regulators in response to internal and external cues to bring about transcriptional changes, thus serving as the critical link between an organism's genome and its phenotypic traits. Deciphering the biology of enhancers is a key to understanding the genetic basis of common human diseases. Although a large number of candidate enhancers have been annotated through genome-wide analyses of chromatin accessibility, transcription factor binding, and histone modification in diverse cell types, efforts to characterize their biological roles in human diseases have only begun. Recent experiments have suggested a role for the three-dimensional chromatin architecture in regulation of gene expression by enhancers.

Massively parallel cis-regulatory analysis in the mammalian central nervous system

Cis-regulatory elements (CREs, e.g., promoters and enhancers) regulate gene expression, and variants within CREs can modulate disease risk. Next-generation sequencing has enabled the rapid generation of genomic data that predict the locations of CREs, but a bottleneck lies in functionally interpreting these data. To address this issue, massively parallel reporter assays (MPRAs) have emerged, in which barcoded reporter libraries are introduced into cells, and the resulting barcoded transcripts are quantified by next-generation sequencing. Thus far, MPRAs have been largely restricted to assaying short CREs in a limited repertoire of cultured cell types. Here, we present two advances that extend the biological relevance and applicability of MPRAs. First, we adapt exome capture technology to instead capture candidate CREs, thereby tiling across the targeted regions and markedly increasing the length of CREs that can be readily assayed. Second, we package the library into adeno-associated virus (AAV), thereby allowing delivery to target organs in vivo. As a proof of concept, we introduce a capture library of about 46,000 constructs, corresponding to roughly 3500 DNase I hypersensitive (DHS) sites, into the mouse retina by ex vivo plasmid electroporation and into the mouse cerebral cortex by in vivo AAV injection. We demonstrate tissue-specific cis-regulatory activity of DHSs and provide examples of high-resolution truncation mutation analysis for multiplex parsing of CREs. Our approach should enable massively parallel functional analysis of a wide range of CREs in any organ or species that can be infected by AAV, such as nonhuman primates and human stem cell–derived organoids.

Glioblastoma Stem Cells Respond to Differentiation Cues but Fail to Undergo Commitment and Terminal Cell-Cycle Arrest

Glioblastoma (GBM) is an aggressive brain tumor whose growth is driven by stemcell-like cells. BMP signaling triggers cell-cycle exit and differentiation of GBM stemcells (GSCs) and, therefore, might have therapeutic value. However, the epigenetic mechanisms that accompany differentiation remain poorly defined. It is also unclear whether cell-cycle arrest is terminal. Herewe find only a subset ofGSCcultures exhibit astrocyte differentiation in response to BMP. Although overtly differentiated non-cycling astrocytes are generated, they remain vulnerable to cell-cycle re-entry and fail to appropriately reconfigure DNA methylation patterns. Chromatin accessibility mapping identified loci that failed to alter in response to BMP and these were enriched in SOX transcription factor-binding motifs. SOX transcription factors, therefore, may limit differentiation commitment. A similar propensity for cell-cycle re-entry and de-differentiation was observed in GSC-derived oligodendrocyte-like cells. These findings highlight significant obstacles to BMP-induced differentiation as therapy forGBM.

Cellular and chromatin dynamics of antibody-secreting plasma cells

Plasma cells are terminally differentiated B cells responsible for maintaining protective serum antibody titers. Despite their clinical importance, our understanding of the linear genomic features and chromatin structure of plasma cells is incomplete. The plasma cell differentiation program can be triggered by different signals and in multiple, diverse peripheral B cell subsets. This heterogeneity raises questions about the gene regulatory circuits required for plasma cell specification. Recently, new regulators of plasma cell differentiation have been identified and the enhancer landscapes of naïve B cells have been described. Other studies have revealed that the bone marrow niche harbors heterogeneous plasma cell subsets. Still undefined are the minimal requirements to become a plasma cell and what molecular features make peripheral B cell subsets competent to become antibody-secreting plasma cells. New technologies promise to reveal underlying chromatin configurations that promote efficient antibody secretion. For further resources related to this article, please visit the WIREs website.

Dose-dependent role of the cohesin complex in normal and malignant hematopoiesis

Cohesin complex members have recently been identified as putative tumor suppressors in hematologic and epithelial malignancies. The cohesin complex guides chromosome segregation; however, cohesin mutant leukemias do not show genomic instability. We hypothesized that reduced cohesin function alters chromatin structure and disrupts cis-regulatory architecture of hematopoietic progenitors. We investigated the consequences of Smc3 deletion in normal and malignant hematopoiesis. Biallelic Smc3 loss induced bone marrow aplasia with premature sister chromatid separation and revealed an absolute requirement for cohesin in hematopoietic stem cell (HSC) function. In contrast, Smc3 haploinsufficiency increased self-renewal in vitro and in vivo, including competitive transplantation. Smc3 haploinsufficiency reduced coordinated transcriptional output, including reduced expression of transcription factors and other genes associated with lineage commitment. Smc3 haploinsufficiency cooperated with Flt3-ITD to induce acute leukemia in vivo, with potentiated Stat5 signaling and altered nucleolar topology. These data establish a dose dependency for cohesin in regulating chromatin structure and HSC function.

Prioritization Of Nonsynonymous Single Nucleotide Variants For Exome Sequencing Studies Via Integrative Learning On Multiple Genomic Data

The rapid advancement of next generation sequencing technology has greatly accelerated the progress for understanding human inherited diseases via such innovations as exome sequencing. Nevertheless, the identification of causative variants from sequencing data remains a great challenge. Traditional statistical genetics approaches such as linkage analysis and association studies have limited power in analyzing exome sequencing data, while relying on simply filtration strategies and predicted functional implications of mutations to pinpoint pathogenic variants are prone to produce false positives. To overcome these limitations, we herein propose a supervised learning approach, termed snvForest, to prioritize candidate nonsynonymous single nucleotide variants for a specific type of disease by integrating 11 functional scores at the variant level and 8 association scores at the gene level. We conduct a series of large-scale in silico validation experiments, demonstrating the effectiveness of snvForest across 2,511 diseases of different inheritance styles and the superiority of our approach over two state-of-the-art methods. We further apply snvForest to three real exome sequencing data sets of epileptic encephalophathies and intellectual disability to show the ability of our approach to identify causative de novo mutations for these complex diseases. The online service and standalone software of snvForest are found at http://bioinfo.au.tsinghua.edu.cn/jianglab/snvforest.

Dynamics of chromatin accessibility and epigenetic state in response to UV damage

Epigenetic mechanisms determine the access of regulatory factors to DNA during events such as transcription and the DNA damage response. However, the global response of histone modifications and chromatin accessibility to UV exposure remains poorly understood. Here, we report that UV exposure results in a genome-wide reduction in chromatin accessibility, while the distribution of the active regulatory mark H3K27ac undergoes massive reorganization. Genomic loci subjected to epigenetic reprogramming upon UV exposure represent target sites for sequence-specific transcription factors. Most of these are distal regulatory regions, highlighting their importance in the cellular response to UV exposure. Furthermore, UV exposure results in an extensive reorganization of super-enhancers, accompanied by expression changes of associated genes, which may in part contribute to the stress response. Taken together, our study provides the first comprehensive resource for genome-wide chromatin changes upon UV irradiation in relation to gene expression and elucidates new aspects of this relationship.

Transcriptional enhancers: functional insights and role in human disease

In recent years, studies of cis-regulatory mechanisms have evolved from a predominant focus on promoter regions to the realization that spatial and temporal gene regulation is frequently driven by long-range enhancer clusters that operate within chromosomal compartments. This increased understanding of genome function, together with the emergence of technologies that enable whole-genome sequencing of patients’ DNAs, open the prospect of dissecting the role of cis-regulatory defects in human disease. In this review we discuss how recent epigenomic studies have provided insights into the function of transcriptional enhancers. We then present examples that illustrate how integrative genomics can help uncover enhancer sequence variants underlying Mendelian and common polygenic human disease.

msCentipede: Modeling Heterogeneity across Genomic Sites and Replicates Improves Accuracy in the Inference of Transcription Factor Binding

Understanding global gene regulation depends critically on accurate annotation of regulatory elements that are functional in a given cell type. CENTIPEDE, a powerful, probabilistic framework for identifying transcription factor binding sites from tissue-specific DNase I cleavage patterns and genomic sequence content, leverages the hypersensitivity of factor-bound chromatin and the information in the DNase I spatial cleavage profile characteristic of each DNA binding protein to accurately infer functional factor binding sites. However, the model for the spatial profile in this framework fails to account for the substantial variation in the DNase I cleavage profiles across different binding sites. Neither does it account for variation in the profiles at the same binding site across multiple replicate DNase I experiments, which are increasingly available. In this work, we introduce new methods, based on multi-scale models for inhomogeneous Poisson processes, to account for such variation in DNase I cleavage patterns both within and across binding sites. These models account for the spatial structure in the heterogeneity in DNase I cleavage patterns for each factor. Using DNase-seq measurements assayed in a lymphoblastoid cell line, we demonstrate the improved performance of this model for several transcription factors by comparing against the Chip-seq peaks for those factors. Finally, we explore the effects of DNase I sequence bias on inference of factor binding using a simple extension to our framework that allows for a more flexible background model. The proposed model can also be easily applied to paired-end ATAC-seq and DNase-seq data. msCentipede, a Python implementation of our algorithm, is available at http://rajanil.github.io/msCentipede.

Hematopoietic Differentiation Is Required for Initiation of Acute Myeloid Leukemia

Mutations in acute myeloid leukemia (AML)-associated oncogenes often arise in hematopoietic stem cells (HSCs) and promote acquisition of leukemia stem cell (LSC) phenotypes. However, as LSCs often share features of lineage-restricted progenitors, the relative contribution of differentiation status to LSC transformation is unclear. Using murine MLL-AF9 and MOZ-TIF2 AML models, we show that myeloid differentiation to granulocyte macrophage progenitors (GMPs) is critical for LSC generation. Disrupting GMP formation by deleting the lineage-restricted transcription factor C/EBPa blocked normal granulocyte formation and prevented initiation of AML. However, restoring myeloid differentiation in C/EBPa mutants with inflammatory cytokines reestablished AML transformation capacity. Genomic analyses of GMPs, including gene expression and H3K79me2 profiling in conjunction with ATAC-seq, revealed a permissive genomic environment for activation of a minimal transcription program shared by GMPs and LSCs. Together, these findings show that myeloid differentiation is a prerequisite for LSC formation and AML development, providing insights for therapeutic development.

Early enhancer establishment and regulatory locus complexity shape transcriptional programs in hematopoietic differentiation

We carried out an integrative analysis of enhancer landscape and gene expression dynamics during hematopoietic differentiation using DNase-seq, histone mark ChIP-seq and RNA sequencing to model how the early establishment of enhancers and regulatory locus complexity govern gene expression changes at cell state transitions. We found that high-complexity genes-those with a large total number of DNase-mapped enhancers across the lineage-differ architecturally and functionally from low-complexity genes, achieve larger expression changes and are enriched for both cell type-specific and transition enhancers, which are established in hematopoietic stem and progenitor cells and maintained in one differentiated cell fate but lost in others. We then developed a quantitative model to accurately predict gene expression changes from the DNA sequence content and lineage history of active enhancers. Our method suggests a new mechanistic role for PU.1 at transition peaks during B cell specification and can be used to correct assignments of enhancers to genes.

Bacterial infection remodels the DNA methylation landscape of human dendritic cells

DNA methylation is an epigenetic mark thought to be robust to environmental perturbations on a short time scale. Here, we challenge that view by demonstrating that the infection of human dendritic cells (DCs) with a live pathogenic bacteria is associated with rapid and active demethylation at thousands of loci, independent of cell division. We performed an integrated analysis of data on genome-wide DNA methylation, histone mark patterns, chromatin accessibility, and gene expression, before and after infection. We found that infection-induced demethylation rarely occurs at promoter regions and instead localizes to distal enhancer elements, including those that regulate the activation of key immune transcription factors. Active demethylation is associated with extensive epigenetic remodeling, including the gain of histone activation marks and increased chromatin accessibility, and is strongly predictive of changes in the expression levels of nearby genes. Collectively, our observations show that active, rapid changes in DNA methylation in enhancers play a previously unappreciated role in regulating the transcriptional response to infection, even in nonproliferating cells.

Genomic approaches to studying human-specific developmental traits

Changes in developmental regulatory programs drive both disease and phenotypic differences among species. Linking human-specific traits to alterations in development is challenging, because we have lacked the tools to assay and manipulate regulatory networks in human and primate embryonic cells. This field was transformed by the sequencing of hundreds of genomes – human and non-human – that can be compared to discover the regulatory machinery of genes involved in human development. This approach has identified thousands of human-specific genome alterations in developmental genes and their regulatory regions. With recent advances in stem cell techniques, genome engineering, and genomics, we can now test these sequences for effects on developmental gene regulation and downstream phenotypes in human cells and tissues.

Enhancer divergence and cis-regulatory evolution in the human and chimp neural crest

cis-regulatory changes play a central role in morphological divergence, yet the regulatory principles underlying emergence of human traits remain poorly understood. Here, we use epigenomic profiling from human and chimpanzee cranial neural crest cells to systematically and quantitatively annotate divergence of craniofacial cis-regulatory landscapes. Epigenomic divergence is often attributable to genetic variation within TF motifs at orthologous enhancers, with a novel motif being most predictive of activity biases. We explore properties of this cis-regulatory change, revealing the role of particular retroelements, uncovering broad clusters of species-biased enhancers near genes associated with human facial variation, and demonstrating that cis-regulatory divergence is linked to quantitative expression differences of crucial neural crest regulators. Our work provides a wealth of candidates for future evolutionary studies and demonstrates the value of "cellular anthropology," a strategy of using in-vitro-derived embryonic cell types to elucidate both fundamental and evolving mechanisms underlying morphological variation in higher primates.

Transcriptional Enhancers in the Regulation of T Cell Differentiation

The changes in phenotype and function that characterize the differentiation of naïve T cells to effector and memory states are underscored by large-scale, coordinated, and stable changes in gene expression. In turn, these changes are choreographed by the interplay between transcription factors and epigenetic regulators that act to restructure the genome, ultimately ensuring lineage-appropriate gene expression. Here, we focus on the mechanisms that control T cell differentiation, with a particular focus on the role of regulatory elements encoded within the genome, known as transcriptional enhancers (TEs). We discuss the central role of TEs in regulating T cell differentiation, both in health and disease.

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.

Epigenomic Landscapes Reflect Neuronal Diversity

Epigenomic profiling of complex tissues obscures regulatory elements that distinguish one cell type from another. In this issue of Neuron, Mo et al. (2015) apply cell-type-specific profiling to mouse neuronal subtypes and discover an unprecedented level of neuronal diversity.

Single-cell transcriptome analysis reveals coordinated ectopic gene-expression patterns in medullary thymic epithelial cells

Expression of tissue-restricted self antigens (TRAs) in medullary thymic epithelial cells (mTECs) is essential for the induction of self-tolerance and prevents autoimmunity, with each TRA being expressed in only a few mTECs. How this process is regulated in single mTECs and is coordinated at the population level, such that the varied single-cell patterns add up to faithfully represent TRAs, is poorly understood. Here we used single-cell RNA sequencing and obtained evidence of numerous recurring TRA-co-expression patterns, each present in only a subset of mTECs. Co-expressed genes clustered in the genome and showed enhanced chromatin accessibility. Our findings characterize TRA expression in mTECs as a coordinated process that might involve local remodeling of chromatin and thus ensures a comprehensive representation of the immunological self.

Exploiting genomics and natural genetic variation to decode macrophage enhancers

The mammalian genome contains on the order of a million enhancer-like regions that are required to establish the identities and functions of specific cell types. Here, we review recent studies in immune cells that have provided insight into the mechanisms that selectively activate certain enhancers in response to cell lineage and environmental signals. We describe a working model wherein distinct classes of transcription factors define the repertoire of active enhancers in macrophages through collaborative and hierarchical interactions, and discuss important challenges to this model, specifically providing examples from T cells. We conclude by discussing the use of natural genetic variation as a powerful approach for decoding transcription factor combinations that play dominant roles in establishing the enhancer landscapes, and the potential that these insights have for advancing our understanding of the molecular causes of human disease.

Genomic-wide transcriptional profiling in primary myoblasts reveals Runx1-regulated genes in muscle regeneration

In response to muscle damage the muscle adult stem cells are activated and differentiate into myoblasts that regenerate the damaged tissue. We have recently showed that following myopathic damage the level of the Runx1 transcription factor (TF) is elevated and that during muscle regeneration this TF regulates the balance between myoblast proliferation and differentiation (Umansky et al.). We employed Runx1-dependent gene expression, Chromatin Immunoprecipitation sequencing (ChIP-seq), Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) and histone H3K4me1/H3K27ac modification analyses to identify a subset of Runx1-regulated genes that are co-occupied by the TFs MyoD and c-Jun and are involved in muscle regeneration (Umansky et al.). The data is available at the GEO database under the superseries accession number GSE56131.