Chromatin architecture reorganization during stem cell differentiation

Intrinsic Molecular Processes: Impact on Mutagenesis

Mutations provide resources for genome evolution by generating genetic variability. In addition, mutations act as a driving force leading to disease pathogenesis, and thus have important implications for disease diagnosis, prognosis, and treatment. Understanding the mechanisms underlying how mutations occur is therefore of prime importance for elucidating evolutionary and pathogenic processes. Recent genomics studies have revealed that mutations occur non-randomly across the human genome. In particular, the distribution of mutations is highly associated with intrinsic molecular processes including transcription, chromatin organization, DNA replication timing, and DNA repair. Interplay between intrinsic processes and extrinsic mutagenic exposure may thus imprint a characteristic mutational landscape on tumors. We discuss the impact of intrinsic molecular processes on mutation acquisition in cancer.

Epigenetic landscapes reveal transcription factors that regulate CD8+ T cell differentiation

Dynamic changes in the expression of transcription factors (TFs) can influence the specification of distinct CD8+ T cell fates, but the observation of equivalent expression of TFs among differentially fated precursor cells suggests additional underlying mechanisms. Here we profiled the genome-wide histone modifications, open chromatin and gene expression of naive, terminal-effector, memory-precursor and memory CD8+ T cell populations induced during the in vivo response to bacterial infection. Integration of these data suggested that the expression and binding of TFs contributed to the establishment of subset-specific enhancers during differentiation. We developed a new bioinformatics method using the PageRank algorithm to reveal key TFs that influence the generation of effector and memory populations. The TFs YY1 and Nr3c1, both constitutively expressed during CD8+ T cell differentiation, regulated the formation of terminal-effector cell fates and memory-precursor cell fates, respectively. Our data define the epigenetic landscape of differentiation intermediates and facilitate the identification of TFs with previously unappreciated roles in CD8+ T cell differentiation.

Alterations in Three-Dimensional Organization of the Cancer Genome and Epigenome

The structural and functional basis of the genome is provided by the three-dimensional (3D) chromatin state. To enable accurate gene regulation, enhancer elements and promoter regions are brought into close spatial proximity to ensure proper, cell type-specific gene expression. In cancer, genetic and epigenetic processes can deregulate the transcriptional program. To investigate whether the 3D chromatin state is also disrupted in cancer we performed Hi-C chromosome conformation sequencing in normal and prostate cancer cells and compared the chromatin interaction maps with changes to the genome and epigenome. Notably, we find that additional topologically associated domain (TAD) boundaries are formed in cancer cells resulting in smaller TADs and altered gene expression profiles. The new TAD boundaries are commonly associated with copy-number changes observed in the cancer genome. We also identified new cancer-specific chromatin loops within TADs that are enriched for enhancers and promoters. Finally, we find that many of the long-range epigenetically silenced (LRES) and long-range epigenetically active (LREA) regions in cancer are characterized by differential chromatin interactions. Together our data provide a new insight into charting alterations in higher-order structure and the relationship with genetic, epigenetic, and transcriptional changes across the cancer genome.

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.

A tiling-deletion-based genetic screen for cis-regulatory element identification in mammalian cells

Millions of cis-regulatory elements are predicted in the human genome, but direct evidence for their biological function is still scarce. Here we report a high-throughput method, Cis-Regulatory Element Scan by Tiling-deletion and sequencing (CREST-seq), for unbiased discovery and functional assessment of cis regulatory sequences in the genome. We use it to interrogate the 2Mbp POU5F1 locus in the human embryonic stem cells and identify 45 cis-regulatory elements of POU5F1. A majority of these elements display active chromatin marks, DNase hypersensitivity and occupancy by multiple transcription factors, confirming the utility of chromatin signatures in cis elements mapping. Notably, 17 of them are previously annotated promoters of functionally unrelated genes, and like typical enhancers, they form extensive spatial contacts with the POU5F1 promoter. Taken together, these results support the utility of CREST-seq for large-scale cis regulatory element discovery and point to commonality of enhancer-like promoters in the human genome.

Gene Regulatory Elements, Major Drivers of Human Disease

Gene expression changes, the driving forces for cellular diversity in multicellular organisms, are regulated by a diverse set of gene regulatory elements that direct transcription in specific cells. Mutations in these elements, ranging from chromosomal aberrations to single-nucleotide polymorphisms, are a major cause of human disease. However, we currently have a very limited understanding of how regulatory element genotypes lead to specific phenotypes. In this review, we discuss the various methods of regulatory element identification, the different types of mutations they harbor, and their impact on human disease. We highlight how these variations can affect transcription of multiple genes in gene regulatory networks. In addition, we describe how novel technologies, such as massively parallel reporter assays and CRISPR/Cas9 genome editing, are beginning to provide a better understanding of the functional roles that these elements have and how their alteration can lead to specific phenotypes.

Multi-scale chromatin state annotation using a hierarchical hidden Markov model

Chromatin-state analysis is widely applied in the studies of development and diseases. However, existing methods operate at a single length scale, and therefore cannot distinguish large domains from isolated elements of the same type. To overcome this limitation, we present a hierarchical hidden Markov model, diHMM, to systematically annotate chromatin states at multiple length scales. We apply diHMM to analyse a public ChIP-seq data set. diHMM not only accurately captures nucleosome-level information, but identifies domain-level states that vary in nucleosome-level state composition, spatial distribution and functionality. The domain-level states recapitulate known patterns such as super-enhancers, bivalent promoters and Polycomb repressed regions, and identify additional patterns whose biological functions are not yet characterized. By integrating chromatin-state information with gene expression and Hi-C data, we identify context-dependent functions of nucleosome-level states. Thus, diHMM provides a powerful tool for investigating the role of higher-order chromatin structure in gene regulation.

The impact of chromatin structure on gene expression makes it integral to our understanding of developmental and disease processes. Here, the authors introduce a hierarchical hidden Markov model to systematically annotate chromatin states at multiple length scales, and demonstrate its utility for the elucidation of the role of chromatin structure in gene expression.

Pgltools: a genomic arithmetic tool suite for manipulation of Hi-C peak and other chromatin interaction data

Background

Genomic interaction studies use next-generation sequencing (NGS) to examine the interactions between two loci on the genome, with subsequent bioinformatics analyses typically including annotation, intersection, and merging of data from multiple experiments. While many file types and analysis tools exist for storing and manipulating single locus NGS data, there is currently no file standard or analysis tool suite for manipulating and storing paired-genomic-loci: the data type resulting from “genomic interaction” studies. As genomic interaction sequencing data are becoming prevalent, a standard file format and tools for working with these data conveniently and efficiently are needed.

Results

This article details a file standard and novel software tool suite for working with paired-genomic-loci data. We present the paired-genomic-loci (PGL) file standard for genomic-interactions data, and the accompanying analysis tool suite “pgltools”: a cross platform, pypy compatible python package available both as an easy-to-use UNIX package, and as a python module, for integration into pipelines of paired-genomic-loci analyses.

Conclusions

Pgltools is a freely available, open source tool suite for manipulating paired-genomic-loci data. Source code, an in-depth manual, and a tutorial are available publicly at www.github.com/billgreenwald/pgltools, and a python module of the operations can be installed from PyPI via the PyGLtools module.

Single-nucleus Hi-C reveals unique chromatin reorganization at oocyte-to-zygote transition

Chromatin is reprogrammed after fertilization to produce a totipotent zygote with the potential to generate a new organism. The maternal genome inherited from the oocyte and the paternal genome provided by sperm coexist as separate haploid nuclei in the zygote. How these two epigenetically distinct genomes are spatially organized is poorly understood. Existing chromosome conformation capture-based methods are not applicable to oocytes and zygotes owing to a paucity of material. To study three-dimensional chromatin organization in rare cell types, we developed a single-nucleus Hi-C (high-resolution chromosome conformation capture) protocol that provides greater than tenfold more contacts per cell than the previous method. Here we show that chromatin architecture is uniquely reorganized during the oocyte-to-zygote transition in mice and is distinct in paternal and maternal nuclei within single-cell zygotes. Features of genomic organization including compartments, topologically associating domains (TADs) and loops are present in individual oocytes when averaged over the genome, but the presence of each feature at a locus varies between cells. At the sub-megabase level, we observed stochastic clusters of contacts that can occur across TAD boundaries but average into TADs. Notably, we found that TADs and loops, but not compartments, are present in zygotic maternal chromatin, suggesting that these are generated by different mechanisms. Our results demonstrate that the global chromatin organization of zygote nuclei is fundamentally different from that of other interphase cells. An understanding of this zygotic chromatin 'ground state' could potentially provide insights into reprogramming cells to a state of totipotency.

Linkages between changes in the 3D organization of the genome and transcription during myotube differentiation in vitro

Background

The spatial organization of eukaryotic genomes facilitates and reflects the underlying nuclear processes that are occurring in the cell. As such, the spatial organization of a genome represents a window on the genome biology that enables analysis of the nuclear regulatory processes that contribute to mammalian development.

Methods

In this study, Hi-C and RNA-seq were used to capture the genome organization and transcriptome in mouse muscle progenitor cells (C2C12 myoblasts) before and after differentiation to myotubes, in the presence or absence of the cytidine analogue AraC.

Results

We observed significant local and global developmental changes despite high levels of correlation between the myotubes and myoblast genomes. Notably, the genes that exhibited the greatest variation in transcript levels between the different developmental stages were predominately within the euchromatic compartment. There was significant re-structuring and changes in the expression of replication-dependent histone variants within the HIST1 locus. Finally, treating terminally differentiated myotubes with AraC resulted in additional changes to the transcriptome and 3D genome organization of sets of genes that were all involved in pyroptosis.

Conclusions

Collectively, our results provide evidence for muscle cell-specific responses to developmental and environmental stimuli mediated through a chromatin structure mechanism.

Electronic supplementary material

The online version of this article (doi:10.1186/s13395-017-0122-1) contains supplementary material, which is available to authorized users.

Visualization of aging-associated chromatin alterations with an engineered TALE system

Visualization of specific genomic loci in live cells is a prerequisite for the investigation of dynamic changes in chromatin architecture during diverse biological processes, such as cellular aging. However, current precision genomic imaging methods are hampered by the lack of fluorescent probes with high specificity and signal-to-noise contrast. We find that conventional transcription activator-like effectors (TALEs) tend to form protein aggregates, thereby compromising their performance in imaging applications. Through screening, we found that fusing thioredoxin with TALEs prevented aggregate formation, unlocking the full power of TALE-based genomic imaging. Using thioredoxin-fused TALEs (TTALEs), we achieved high-quality imaging at various genomic loci and observed aging-associated (epi) genomic alterations at telomeres and centromeres in human and mouse premature aging models. Importantly, we identified attrition of ribosomal DNA repeats as a molecular marker for human aging. Our study establishes a simple and robust imaging method for precisely monitoring chromatin dynamics in vitro and in vivo.

Potential energy landscapes identify the information-theoretic nature of the epigenome

Epigenetics is the study of biochemical modifications carrying information independent of DNA sequence, which are heritable through cell division. In 1940, Waddington coined the term "epigenetic landscape" as a metaphor for pluripotency and differentiation, but methylation landscapes have not yet been rigorously computed. Using principles from statistical physics and information theory, we derive epigenetic energy landscapes from whole-genome bisulfite sequencing (WGBS) data that enable us to quantify methylation stochasticity genome-wide using Shannon's entropy, associating it with chromatin structure. Moreover, we consider the Jensen-Shannon distance between sample-specific energy landscapes as a measure of epigenetic dissimilarity and demonstrate its effectiveness for discerning epigenetic differences. By viewing methylation maintenance as a communications system, we introduce methylation channels and show that higher-order chromatin organization can be predicted from their informational properties. Our results provide a fundamental understanding of the information-theoretic nature of the epigenome that leads to a powerful approach for studying its role in disease and aging.

Global reorganisation of cis-regulatory units upon lineage commitment of human embryonic stem cells

Long-range cis-regulatory elements such as enhancers coordinate cell-specific transcriptional programmes by engaging in DNA looping interactions with target promoters. Deciphering the interplay between the promoter connectivity and activity of cis-regulatory elements during lineage commitment is crucial for understanding developmental transcriptional control. Here, we use Promoter Capture Hi-C to generate a high-resolution atlas of chromosomal interactions involving ~22,000 gene promoters in human pluripotent and lineage-committed cells, identifying putative target genes for known and predicted enhancer elements. We reveal extensive dynamics of cis-regulatory contacts upon lineage commitment, including the acquisition and loss of promoter interactions. This spatial rewiring occurs preferentially with predicted changes in the activity of cis-regulatory elements and is associated with changes in target gene expression. Our results provide a global and integrated view of promoter interactome dynamics during lineage commitment of human pluripotent cells.

Diverse Non-genetic, Allele-Specific Expression Effects Shape Genetic Architecture at the Cellular Level in the Mammalian Brain

Interactions between genetic and epigenetic effects shape brain function, behavior, and the risk for mental illness. Random X inactivation and genomic imprinting are epigenetic allelic effects that are well known to influence genetic architecture and disease risk. Less is known about the nature, prevalence, and conservation of other potential epigenetic allelic effects in vivo in the mouse and primate brain. Here we devise genomics, in situ hybridization, and mouse genetics strategies to uncover diverse allelic effects in the brain that are not caused by imprinting or genetic variation. We found allelic effects that are developmental stage and cell type specific, that are prevalent in the neonatal brain, and that cause mosaics of monoallelic brain cells that differentially express wild-type and mutant alleles for heterozygous mutations. Finally, we show that diverse non-genetic allelic effects that impact mental illness risk genes exist in the macaque and human brain. Our findings have potential implications for mammalian brain genetics. VIDEO ABSTRACT.

Erythroid specific activator GATA-1-dependent interactions between CTCF sites around the β-globin locus

CTCF sites (binding motifs for CCCTC-binding factor, an insulator protein) are located considerable distances apart on genomes but are closely positioned in organized chromatin. The close positioning of CTCF sites is often cell type or tissue specific. Here we analyzed chromatin organization in eight CTCF sites around the β-globin locus by 3C assay and explored the roles of erythroid specific transcription activator GATA-1 and KLF1 in it. It was found five CTCF sites convergent to the locus interact with each other in erythroid K562 cells but not in non-erythroid 293 cells. The interaction was decreased by depletion of GATA-1 or KLF1. It accompanied reductions of CTCF and Rad21 occupancies and loss of active chromatin structure at the CTCF sites. Furthermore Rad21 occupancy was reduced in the β-globin locus control region (LCR) hypersensitive sites (HSs) by the depletion of GATA-1 or KLF1. The role of GATA-1 in interaction between CTCF sites was revealed by its ectopic expression in 293 cells and by deletion of a GATA-1 site in the LCR HS2. These findings indicate that erythroid specific activator GATA-1 acts at CTCF sites around the β-globin locus to establish tissue-specific chromatin organization.

Chrom3D: three-dimensional genome modeling from Hi-C and nuclear lamin-genome contacts

Current three-dimensional (3D) genome modeling platforms are limited by their inability to account for radial placement of loci in the nucleus. We present Chrom3D, a user-friendly whole-genome 3D computational modeling framework that simulates positions of topologically-associated domains (TADs) relative to each other and to the nuclear periphery. Chrom3D integrates chromosome conformation capture (Hi-C) and lamin-associated domain (LAD) datasets to generate structure ensembles that recapitulate radial distributions of TADs detected in single cells. Chrom3D reveals unexpected spatial features of LAD regulation in cells from patients with a laminopathy-causing lamin mutation. Chrom3D is freely available on github.

Rfx2 Stabilizes Foxj1 Binding at Chromatin Loops to Enable Multiciliated Cell Gene Expression

Cooperative transcription factor binding at cis-regulatory sites in the genome drives robust eukaryotic gene expression, and many such sites must be coordinated to produce coherent transcriptional programs. The transcriptional program leading to motile cilia formation requires members of the DNA-binding forkhead (Fox) and Rfx transcription factor families and these factors co-localize to cilia gene promoters, but it is not clear how many cilia genes are regulated by these two factors, whether these factors act directly or indirectly, or how these factors act with specificity in the context of a 3-dimensional genome. Here, we use genome-wide approaches to show that cilia genes reside at the boundaries of topological domains and that these areas have low enhancer density. We show that the transcription factors Foxj1 and Rfx2 binding occurs in the promoters of more cilia genes than other known cilia transcription factors and that while Rfx2 binds directly to promoters and enhancers equally, Foxj1 prefers direct binding to enhancers and is stabilized at promoters by Rfx2. Finally, we show that Rfx2 and Foxj1 lie at the anchor endpoints of chromatin loops, suggesting that target genes are activated when Foxj1 bound at distal sites is recruited via a loop created by Rfx2 binding at both sites. We speculate that the primary function of Rfx2 is to stabilize distal enhancers with proximal promoters by operating as a scaffolding factor, bringing key regulatory domains bound by Foxj1 into close physical proximity and enabling coordinated cilia gene expression.

The multiciliated cell extends hundreds of motile cilia to produce fluid flow in the airways and other organ systems. The formation of this specialized cell type requires the coordinated expression of hundreds of genes in order to produce all the protein parts motile cilia require. While a relatively small number of transcription factors has been identified that promote gene expression during multiciliate cell differentiation, it is not clear how they work together to coordinate the expression of genes required for multiple motile ciliation. Here, we show that two transcription factors known to drive cilia formation, Foxj1 and Rfx2, play complementary roles wherein Foxj1 activates target genes but tends not to bind near them in the genome, whereas Rfx2 can’t activate target genes by itself but instead acts as a scaffold by localizing Foxj1 to the proper targets. These results suggest not only a mechanism by which complex gene expression is coordinated in multiciliated cells, but also how transcriptional programs in general could be modular and deployed across different cellular contexts with the same basic promoter configuration.

Comprehensive characterization of neutrophil genome topology

Neutrophils are responsible for the first line of defense against invading pathogens. Their nuclei are uniquely structured as multiple lobes that establish a highly constrained nuclear environment. Here we found that neutrophil differentiation was not associated with large-scale changes in the number and sizes of topologically associating domains (TADs). However, neutrophil genomes were enriched for long-range genomic interactions that spanned multiple TADs. Population-based simulation of spherical and toroid genomes revealed declining radii of gyration for neutrophil chromosomes. We found that neutrophil genomes were highly enriched for heterochromatic genomic interactions across vast genomic distances, a process named supercontraction. Supercontraction involved genomic regions located in the heterochromatic compartment in both progenitors and neutrophils or genomic regions that switched from the euchromatic to the heterochromatic compartment during neutrophil differentiation. Supercontraction was accompanied by the repositioning of centromeres, pericentromeres, and long interspersed nuclear elements (LINEs) to the neutrophil nuclear lamina. We found that Lamin B receptor expression was required to attach centromeric and pericentromeric repeats but not LINE-1 elements to the lamina. Differentiating neutrophils also repositioned ribosomal DNA and mininucleoli to the lamina-a process that was closely associated with sharply reduced ribosomal RNA expression. We propose that large-scale chromatin reorganization involving supercontraction and recruitment of heterochromatin and nucleoli to the nuclear lamina facilitates the folding of the neutrophil genome into a confined geometry imposed by a multilobed nuclear architecture.

Perspectives: using polymer modeling to understand the formation and function of nuclear compartments

Compartmentalization is a ubiquitous feature of cellular function. In the nucleus, early observations revealed a non-random spatial organization of the genome with a large-scale segregation between transcriptionally active—euchromatin—and silenced—heterochromatin—parts of the genome. Recent advances in genome-wide mapping and imaging techniques have strikingly improved the resolution at which nuclear genome folding can be analyzed and have revealed a multiscale spatial compartmentalization with increasing evidences that such compartment may indeed result from and participate to genome function. Understanding the underlying mechanisms of genome folding and in particular the link to gene regulation requires a cross-disciplinary approach that combines the new high-resolution techniques with computational modeling of chromatin and chromosomes. In this perspective article, we first present how the copolymer theoretical framework can account for the genome compartmentalization. We then suggest, in a second part, that compartments may act as a “nanoreactor,” increasing the robustness of either activation or repression by enhancing the local concentration of regulators. We conclude with the need to develop a new framework, namely the “living chromatin” model that will allow to explicitly investigate the coupling between spatial compartmentalization and gene regulation.

Dynamic Reorganization of Extremely Long-Range Promoter-Promoter Interactions between Two States of Pluripotency

Serum-to-2i interconversion of mouse embryonic stem cells (mESCs) is a valuable in vitro model for early embryonic development. To assess whether 3D chromatin organization changes during this transition, we established Capture Hi-C with target-sequence enrichment of DNase I hypersensitive sites. We detected extremely long-range intra- and inter-chromosomal interactions between a small subset of H3K27me3 marked bivalent promoters involving the Hox clusters in serum-grown cells. Notably, these promoter-mediated interactions are not present in 2i ground-state pluripotent mESCs but appear upon their further development into primed-like serum mESCs. Reverting serum mESCs to ground-state 2i mESCs removes these promoter-promoter interactions in a spatiotemporal manner. H3K27me3, which is largely absent at bivalent promoters in ground-state 2i mESCs, is necessary, but not sufficient, to establish these interactions, as confirmed by Capture Hi-C on Eed(-/-) serum mESCs. Our results implicate H3K27me3 and PRC2 as critical players in chromatin alteration during priming of ESCs for differentiation.

HiC-bench: comprehensive and reproducible Hi-C data analysis designed for parameter exploration and benchmarking

BACKGROUND

Chromatin conformation capture techniques have evolved rapidly over the last few years and have provided new insights into genome organization at an unprecedented resolution. Analysis of Hi-C data is complex and computationally intensive involving multiple tasks and requiring robust quality assessment. This has led to the development of several tools and methods for processing Hi-C data. However, most of the existing tools do not cover all aspects of the analysis and only offer few quality assessment options. Additionally, availability of a multitude of tools makes scientists wonder how these tools and associated parameters can be optimally used, and how potential discrepancies can be interpreted and resolved. Most importantly, investigators need to be ensured that slight changes in parameters and/or methods do not affect the conclusions of their studies.

RESULTS

To address these issues (compare, explore and reproduce), we introduce HiC-bench, a configurable computational platform for comprehensive and reproducible analysis of Hi-C sequencing data. HiC-bench performs all common Hi-C analysis tasks, such as alignment, filtering, contact matrix generation and normalization, identification of topological domains, scoring and annotation of specific interactions using both published tools and our own. We have also embedded various tasks that perform quality assessment and visualization. HiC-bench is implemented as a data flow platform with an emphasis on analysis reproducibility. Additionally, the user can readily perform parameter exploration and comparison of different tools in a combinatorial manner that takes into account all desired parameter settings in each pipeline task. This unique feature facilitates the design and execution of complex benchmark studies that may involve combinations of multiple tool/parameter choices in each step of the analysis. To demonstrate the usefulness of our platform, we performed a comprehensive benchmark of existing and new TAD callers exploring different matrix correction methods, parameter settings and sequencing depths. Users can extend our pipeline by adding more tools as they become available.

CONCLUSIONS

HiC-bench consists an easy-to-use and extensible platform for comprehensive analysis of Hi-C datasets. We expect that it will facilitate current analyses and help scientists formulate and test new hypotheses in the field of three-dimensional genome organization.

Most brain disease-associated and eQTL haplotypes are not located within transcription factor DNase-seq footprints in brain

Dense genotyping approaches have revealed much about the genetic architecture both of gene expression and disease susceptibility. However, assigning causality to genetic variants associated with a transcriptomic or phenotypic trait presents a far greater challenge. The development of epigenomic resources by ENCODE, the Epigenomic Roadmap and others has led to strategies that seek to infer the likely functional variants underlying these genome-wide association signals. It is known, for example, that such variants tend to be located within areas of open chromatin, as detected by techniques such as DNase-seq and FAIRE-seq. We aimed to assess what proportion of variants associated with phenotypic or transcriptomic traits in the human brain are located within transcription factor binding sites. The bioinformatic tools, Wellington and HINT, were used to infer transcription factor footprints from existing DNase-seq data derived from central nervous system tissues with high spatial resolution. This dataset was then employed to assess the likely contribution of altered transcription factor binding to both expression quantitative trait loci (eQTL) and genome-wide association study (GWAS) signals. Surprisingly, we show that most haplotypes associated with GWAS or eQTL phenotypes are located outside of DNase-seq footprints. This could imply that DNase-seq footprinting is too insensitive an approach to identify a large proportion of true transcription factor binding sites. Importantly, this suggests that prioritising variants for genome engineering studies to establish causality will continue to be frustrated by an inability of footprinting to identify the causative variant within a haplotype.

Disruption of the 3D cancer genome blueprint

Recent advances in chromosome conformation capture technologies are improving the current appreciation of how 3D genome architecture affects its function in different cell types and disease. Long-range chromatin interactions are organized into topologically associated domains, which are known to play a role in constraining gene expression patterns. However, in cancer cells there are alterations in the 3D genome structure, which impacts on gene regulation. Disruption of topologically associated domains architecture can result in alterations in chromatin interactions that bring new regulatory elements and genes together, leading to altered expression of oncogenes and tumor suppressor genes. Here, we discuss the impact of genetic and epigenetic changes in cancer and how this affects the spatial organization of chromatin. Understanding how disruptions to the 3D architecture contribute to the cancer genome will provide novel insights into the principles of epigenetic gene regulation in cancer and mechanisms responsible for cancer associated mutations and rearrangements.

Coupling 1D modifications and 3D nuclear organization: data, models and function

Over the past decade, advances in molecular methods have strikingly improved the resolution at which nuclear genome folding can be analyzed. This revealed a wealth of conserved features organizing the one dimensional DNA molecule into tridimensional nuclear domains. In this review, we briefly summarize the main findings and highlight how models based on polymer physics shed light on the principles underlying the formation of these domains. Finally, we discuss the mechanistic similarities allowing self-organization of these structures and the functional importance of these in the maintenance of transcriptional programs.

Expression of Terminal Effector Genes in Mammalian Neurons Is Maintained by a Dynamic Relay of Transient Enhancers

Generic spinal motor neuron identity is established by cooperative binding of programming transcription factors (TFs), Isl1 and Lhx3, to motor-neuron-specific enhancers. How expression of effector genes is maintained following downregulation of programming TFs in maturing neurons remains unknown. High-resolution exonuclease (ChIP-exo) mapping revealed that the majority of enhancers established by programming TFs are rapidly deactivated following Lhx3 downregulation in stem-cell-derived hypaxial motor neurons. Isl1 is released from nascent motor neuron enhancers and recruited to new enhancers bound by clusters of Onecut1 in maturing neurons. Synthetic enhancer reporter assays revealed that Isl1 operates as an integrator factor, translating the density of Lhx3 or Onecut1 binding sites into transient enhancer activity. Importantly, independent Isl1/Lhx3- and Isl1/Onecut1-bound enhancers contribute to sustained expression of motor neuron effector genes, demonstrating that outwardly stable expression of terminal effector genes in postmitotic neurons is controlled by a dynamic relay of stage-specific enhancers.

A computational strategy to adjust for copy number in tumor Hi-C data

MOTIVATION

The Hi-C technology was designed to decode the three-dimensional conformation of the genome. Despite progress towards more and more accurate contact maps, several systematic biases have been demonstrated to affect the resulting data matrix. Here we report a new source of bias that can arise in tumor Hi-C data, which is related to the copy number of genomic DNA. To address this bias, we designed a chromosome-adjusted iterative correction method called caICB. Our caICB correction method leads to significant improvements when compared with the original iterative correction in terms of eliminating copy number bias.

AVAILABILITY AND IMPLEMENTATION

The method is available at https://bitbucket.org/mthjwu/hicapp CONTACT: michor@jimmy.harvard.eduSupplementary information: Supplementary data are available at Bioinformatics online.

The Genetics of Transcription Factor DNA Binding Variation

Most complex trait-associated variants are located in non-coding regulatory regions of the genome, where they have been shown to disrupt transcription factor (TF)-DNA binding motifs. Variable TF-DNA interactions are therefore increasingly considered as key drivers of phenotypic variation. However, recent genome-wide studies revealed that the majority of variable TF-DNA binding events are not driven by sequence alterations in the motif of the studied TF. This observation implies that the molecular mechanisms underlying TF-DNA binding variation and, by extrapolation, inter-individual phenotypic variation are more complex than originally anticipated. Here, we summarize the findings that led to this important paradigm shift and review proposed mechanisms for local, proximal, or distal genetic variation-driven variable TF-DNA binding. In addition, we discuss the biomedical implications of these findings for our ability to dissect the molecular role(s) of non-coding genetic variants in complex traits, including disease susceptibility.

Decoding the role of regulatory element polymorphisms in complex disease

Genetic variation in gene regulatory elements contributes to diverse human diseases, ranging from rare and severe developmental defects to common and complex diseases such as obesity and diabetes. Early examples of regulatory mechanisms of human diseases involve large chromosomal rearrangements that change the regulatory connections within the genome. Single nucleotide variants in regulatory elements can also contribute to disease, potentially via demonstrated associations with changes in transcription factor binding, enhancer activity, post-translational histone modifications, long-range enhancer-promoter interactions, or RNA polymerase recruitment. Establishing causality between non-coding genetic variants, gene regulation, and disease has recently become more feasible with advances in genome-editing and epigenome-editing technologies. As establishing causal regulatory mechanisms of diseases becomes routine, functional annotation of target genes is likely to emerge as a major bottleneck for translation into patient benefits. In this review, we discuss the history and recent advances in understanding the regulatory mechanisms of human disease, and new challenges likely to be encountered once establishing those mechanisms becomes rote.

Activin/Smad2-induced Histone H3 Lys-27 Trimethylation (H3K27me3) Reduction Is Crucial to Initiate Mesendoderm Differentiation of Human Embryonic Stem Cells

Differentiation of human embryonic stem cells into mesendoderm (ME) is directed by extrinsic signals and intrinsic epigenetic modifications. However, the dynamics of these epigenetic modifications and the mechanisms by which extrinsic signals regulate the epigenetic modifications during the initiation of ME differentiation remain elusive. In this study, we report that levels of histone H3 Lys-27 trimethylation (H3K27me3) decrease during ME initiation, which is essential for subsequent differentiation induced by the combined effects of activin and Wnt signaling. Furthermore, we demonstrate that activin mediates the H3K27me3 decrease via the Smad2-mediated reduction of EZH2 protein level. Our results suggest a two-step process of ME initiation: first, epigenetic priming via removal of H3K27me3 marks and, second, transcription activation. Our findings demonstrate a critical role of H3K27me3 priming and a direct interaction between extrinsic signals and epigenetic modifications during ME initiation.

Characterization of hundreds of regulatory landscapes in developing limbs reveals two regimes of chromatin folding

Complex regulatory landscapes control the pleiotropic transcriptional activities of developmental genes. For most genes, the number, location, and dynamics of their associated regulatory elements are unknown. In this work, we characterized the three-dimensional chromatin microarchitecture and regulatory landscape of 446 limb-associated gene loci in mouse using Capture-C, ChIP-seq, and RNA-seq in forelimb, hindlimb at three developmental stages, and midbrain. The fine mapping of chromatin interactions revealed a strong preference for functional genomic regions such as repressed or active domains. By combining chromatin marks and interaction peaks, we annotated more than 1000 putative limb enhancers and their associated genes. Moreover, the analysis of chromatin interactions revealed two regimes of chromatin folding, one producing interactions stable across tissues and stages and another one associated with tissue and/or stage-specific interactions. Whereas stable interactions associate strongly with CTCF/RAD21 binding, the intensity of variable interactions correlates with changes in underlying chromatin modifications, specifically at the viewpoint and at the interaction site. In conclusion, this comprehensive data set provides a resource for the characterization of hundreds of limb-associated regulatory landscapes and a framework to interpret the chromatin folding dynamics observed during embryogenesis.

Integrating Epigenomics into the Understanding of Biomedical Insight

Epigenetics is one of the most rapidly expanding fields in biomedical research, and the popularity of the high-throughput next-generation sequencing (NGS) highlights the accelerating speed of epigenomics discovery over the past decade. Epigenetics studies the heritable phenotypes resulting from chromatin changes but without alteration on DNA sequence. Epigenetic factors and their interactive network regulate almost all of the fundamental biological procedures, and incorrect epigenetic information may lead to complex diseases. A comprehensive understanding of epigenetic mechanisms, their interactions, and alterations in health and diseases genome widely has become a priority in biological research. Bioinformatics is expected to make a remarkable contribution for this purpose, especially in processing and interpreting the large-scale NGS datasets. In this review, we introduce the epigenetics pioneering achievements in health status and complex diseases; next, we give a systematic review of the epigenomics data generation, summarize public resources and integrative analysis approaches, and finally outline the challenges and future directions in computational epigenomics.

Genome-wide mapping and analysis of chromosome architecture

Chromosomes of eukaryotes adopt highly dynamic and complex hierarchical structures in the nucleus. The three-dimensional (3D) organization of chromosomes profoundly affects DNA replication, transcription and the repair of DNA damage. Thus, a thorough understanding of nuclear architecture is fundamental to the study of nuclear processes in eukaryotic cells. Recent years have seen rapid proliferation of technologies to investigate genome organization and function. Here, we review experimental and computational methodologies for 3D genome analysis, with special focus on recent advances in high-throughput chromatin conformation capture (3C) techniques and data analysis.

Impact of the gut microbiota on enhancer accessibility in gut intraepithelial lymphocytes

The gut microbiota impacts many aspects of host biology including immune function. One hypothesis is that microbial communities induce epigenetic changes with accompanying alterations in chromatin accessibility, providing a mechanism that allows a community to have sustained host effects even in the face of its structural or functional variation. We used Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) to define chromatin accessibility in predicted enhancer regions of intestinal αβ⁺ and γδ⁺ intraepithelial lymphocytes purified from germ-free mice, their conventionally raised (CONV-R) counterparts, and mice reared germ free and then colonized with CONV-R gut microbiota at the end of the suckling-weaning transition. Characterizing genes adjacent to traditional enhancers and super-enhancers revealed signaling networks, metabolic pathways, and enhancer-associated transcription factors affected by the microbiota. Our results support the notion that epigenetic modifications help define microbial community-affiliated functional features of host immune cell lineages.

The p53 Family Coordinates Wnt and Nodal Inputs in Mesendodermal Differentiation of Embryonic Stem Cells

In this study, we outline a regulatory network that involves the p53 tumor suppressor family and the Wnt pathway acting together with the TGF-β pathway in mesendodermal differentiation of mouse and human embryonic stem cells. Knockout of all three members, p53, p63, and p73, shows that the p53 family is essential for mesendoderm specification during exit from pluripotency in embryos and in culture. Wnt3 and its receptor Fzd1 are direct p53 family target genes in this context, and induction of Wnt signaling by p53 is critical for activation of mesendodermal differentiation genes. Globally, Wnt3-activated Tcf3 and nodal-activated Smad2/3 transcription factors depend on each other for co-occupancy of target enhancers associated with key differentiation loci. Our results therefore highlight an unanticipated role for p53 family proteins in a regulatory network that integrates essential Wnt-Tcf and nodal-Smad inputs in a selective and interdependent way to drive mesendodermal differentiation of pluripotent cells.

Insulated Neighborhoods: Structural and Functional Units of Mammalian Gene Control

Understanding how transcriptional enhancers control over 20,000 protein-coding genes to maintain cell-type-specific gene expression programs in all human cells is a fundamental challenge in regulatory biology. Recent studies suggest that gene regulatory elements and their target genes generally occur within insulated neighborhoods, which are chromosomal loop structures formed by the interaction of two DNA sites bound by the CTCF protein and occupied by the cohesin complex. Here, we review evidence that insulated neighborhoods provide for specific enhancer-gene interactions, are essential for both normal gene activation and repression, form a chromosome scaffold that is largely preserved throughout development, and are perturbed by genetic and epigenetic factors in disease. Insulated neighborhoods are a powerful paradigm for gene control that provides new insights into development and disease.

The Shh Topological Domain Facilitates the Action of Remote Enhancers by Reducing the Effects of Genomic Distances

Gene expression often requires interaction between promoters and distant enhancers, which occur within the context of highly organized topologically associating domains (TADs). Using a series of engineered chromosomal rearrangements at the Shh locus, we carried out an extensive fine-scale characterization of the factors that govern the long-range regulatory interactions controlling Shh expression. We show that Shh enhancers act pervasively, yet not uniformly, throughout the TAD. Importantly, changing intra-TAD distances had no impact on Shh expression. In contrast, inversions disrupting the TAD altered global folding of the region and prevented regulatory contacts in a distance-dependent manner. Our data indicate that the Shh TAD promotes distance-independent contacts between distant regions that would otherwise interact only sporadically, enabling functional communication between them. In large genomes where genomic distances per se can limit regulatory interactions, this function of TADs could be as essential for gene expression as the formation of insulated neighborhoods.

A Compendium of Chromatin Contact Maps Reveals Spatially Active Regions in the Human Genome

The three-dimensional configuration of DNA is integral to all nuclear processes in eukaryotes, yet our knowledge of the chromosome architecture is still limited. Genome-wide chromosome conformation capture studies have uncovered features of chromatin organization in cultured cells, but genome architecture in human tissues has yet to be explored. Here, we report the most comprehensive survey to date of chromatin organization in human tissues. Through integrative analysis of chromatin contact maps in 21 primary human tissues and cell types, we find topologically associating domains highly conserved in different tissues. We also discover genomic regions that exhibit unusually high levels of local chromatin interactions. These frequently interacting regions (FIREs) are enriched for super-enhancers and are near tissue-specifically expressed genes. They display strong tissue-specificity in local chromatin interactions. Additionally, FIRE formation is partially dependent on CTCF and the Cohesin complex. We further show that FIREs can help annotate the function of non-coding sequence variants.

Regulation of disease-associated gene expression in the 3D genome

Genetic variation associated with disease often appears in non-coding parts of the genome. Understanding the mechanisms by which this phenomenon leads to disease is necessary to translate results from genetic association studies to the clinic. Assigning function to this type of variation is notoriously difficult because the human genome harbours a complex regulatory landscape with a dizzying array of transcriptional regulatory sequences, such as enhancers that have unpredictable, promiscuous and context-dependent behaviour. In this Review, we discuss how technological advances have provided increasingly detailed information on genome folding; for example, genome folding forms loops that bring enhancers and target genes into close proximity. We also now know that enhancers function within topologically associated domains, which are structural and functional units of chromosomes. Studying disease-associated mutations and chromosomal rearrangements in the context of the 3D genome will enable the identification of dysregulated target genes and aid the progression from descriptive genetic association results to discovering molecular mechanisms underlying disease.

Mapping 3D genome architecture through in situ DNase Hi-C

With the advent of massively parallel sequencing, considerable work has gone into adapting chromosome conformation capture (3C) techniques to study chromosomal architecture at a genome-wide scale. We recently demonstrated that the inactive murine X chromosome adopts a bipartite structure using a novel 3C protocol, termed in situ DNase Hi-C. Like traditional Hi-C protocols, in situ DNase Hi-C requires that chromatin be chemically cross-linked, digested, end-repaired, and proximity-ligated with a biotinylated bridge adaptor. The resulting ligation products are optionally sheared, affinity-purified via streptavidin bead immobilization, and subjected to traditional next-generation library preparation for Illumina paired-end sequencing. Importantly, in situ DNase Hi-C obviates the dependence on a restriction enzyme to digest chromatin, instead relying on the endonuclease DNase I. Libraries generated by in situ DNase Hi-C have a higher effective resolution than traditional Hi-C libraries, which makes them valuable in cases in which high sequencing depth is allowed for, or when hybrid capture technologies are expected to be used. The protocol described here, which involves ∼4 d of bench work, is optimized for the study of mammalian cells, but it can be broadly applicable to any cell or tissue of interest, given experimental parameter optimization.

RUNX1 contributes to higher-order chromatin organization and gene regulation in breast cancer cells

RUNX1 is a transcription factor functioning both as an oncogene and a tumor suppressor in breast cancer. RUNX1 alters chromatin structure in cooperation with chromatin modifier and remodeling enzymes. In this study, we examined the relationship between RUNX1-mediated transcription and genome organization. We characterized genome-wide RUNX1 localization and performed RNA-seq and Hi-C in RUNX1-depleted and control MCF-7 breast cancer cells. RNA-seq analysis showed that RUNX1 depletion led to up-regulation of genes associated with chromatin structure and down-regulation of genes related to extracellular matrix biology, as well as NEAT1 and MALAT1 lncRNAs. Our ChIP-Seq analysis supports a prominent role for RUNX1 in transcriptional activation. About 30% of all RUNX1 binding sites were intergenic, indicating diverse roles in promoter and enhancer regulation and suggesting additional functions for RUNX1. Hi-C analysis of RUNX1-depleted cells demonstrated that overall three-dimensional genome organization is largely intact, but indicated enhanced association of RUNX1 near Topologically Associating Domain (TAD) boundaries and alterations in long-range interactions. These results suggest an architectural role for RUNX1 in fine-tuning local interactions rather than in global organization. Our results provide novel insight into RUNX1-mediated perturbations of higher-order genome organization that are functionally linked with RUNX1-dependent compromised gene expression in breast cancer cells.

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

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

Genomic Alterations of Non-Coding Regions Underlie Human Cancer: Lessons from T-ALL

It has been appreciated for decades that somatic genomic alterations that change coding sequences of proto-oncogenes, translocate enhancers/promoters near proto-oncogenes, or create fusion oncogenes can drive cancer by inducing oncogenic activities. An explosion of genome-wide technologies over the past decade has fueled discoveries of the roles of three-dimensional chromosome structure and powerful cis-acting elements (super-enhancers) in regulating gene transcription. In recent years, studies of human T cell acute lymphoblastic leukemia (T-ALL) using genome-wide technologies have provided paradigms for how non-coding genomic region alterations can disrupt 3D chromosome architecture or establish super-enhancers to activate oncogenic transcription of proto-oncogenes. These studies raise important issues to consider with the objective of leveraging basic knowledge into new diagnostic and therapeutic opportunities for cancer patients.

Chromatin accessibility contributes to simultaneous mutations of cancer genes

Somatic mutations of many cancer genes tend to co-occur (termed co-mutations) in certain patterns during tumor initiation and progression. However, the genetic and epigenetic mechanisms that contribute to the co-mutations of these cancer genes have yet to be explored. Here, we systematically investigated the association between the somatic co-mutations of cancer genes and high-order chromatin conformation. Significantly, somatic point co-mutations in protein-coding genes were closely associated with high-order spatial chromatin folding. We propose that these regions be termed Spatial Co-mutation Hotspots (SCHs) and report their occurrence in different cancer types. The conserved mutational signatures and DNA sequences flanking these point co-mutations, as well as CTCF-binding sites, are also enriched within the SCH regions. The genetic alterations that are harboured in the same SCHs tend to disrupt cancer driver genes involved in multiple signalling pathways. The present work demonstrates that high-order spatial chromatin organisation may contribute to the somatic co-mutations of certain cancer genes during tumor development.

Chromosome conformation elucidates regulatory relationships in developing human brain

Three-dimensional physical interactions within chromosomes dynamically regulate gene expression in a tissue-specific manner. However, the 3D organization of chromosomes during human brain development and its role in regulating gene networks dysregulated in neurodevelopmental disorders, such as autism or schizophrenia, are unknown. Here we generate high-resolution 3D maps of chromatin contacts during human corticogenesis, permitting large-scale annotation of previously uncharacterized regulatory relationships relevant to the evolution of human cognition and disease. Our analyses identify hundreds of genes that physically interact with enhancers gained on the human lineage, many of which are under purifying selection and associated with human cognitive function. We integrate chromatin contacts with non-coding variants identified in schizophrenia genome-wide association studies (GWAS), highlighting multiple candidate schizophrenia risk genes and pathways, including transcription factors involved in neurogenesis, and cholinergic signalling molecules, several of which are supported by independent expression quantitative trait loci and gene expression analyses. Genome editing in human neural progenitors suggests that one of these distal schizophrenia GWAS loci regulates FOXG1 expression, supporting its potential role as a schizophrenia risk gene. This work provides a framework for understanding the effect of non-coding regulatory elements on human brain development and the evolution of cognition, and highlights novel mechanisms underlying neuropsychiatric disorders.