Topologically associating domains are stable units of replication-timing regulation

Spatio-temporal re-organization of replication foci accompanies replication domain consolidation during human pluripotent stem cell lineage specification

Lineage specification of both mouse and human pluripotent stem cells (PSCs) is accompanied by spatial consolidation of chromosome domains and temporal consolidation of their replication timing. Replication timing and chromatin organization are both established during G1 phase at the timing decision point (TDP). Here, we have developed live cell imaging tools to track spatio-temporal replication domain consolidation during differentiation. First, we demonstrate that the fluorescence ubiquitination cell cycle indicator (Fucci) system is incapable of demarcating G1/S or G2/M cell cycle transitions. Instead, we employ a combination of fluorescent PCNA to monitor S phase progression, cytokinesis to demarcate mitosis, and fluorescent nucleotides to label early and late replication foci and track their 3D organization into sub-nuclear chromatin compartments throughout all cell cycle transitions. We find that, as human PSCs differentiate, the length of S phase devoted to replication of spatially clustered replication foci increases, coincident with global compartmentalization of domains into temporally clustered blocks of chromatin. Importantly, re-localization and anchorage of domains was completed prior to the onset of S phase, even in the context of an abbreviated PSC G1 phase. This approach can also be employed to investigate cell fate transitions in single PSCs, which could be seen to differentiate preferentially from G1 phase. Together, our results establish real-time, live-cell imaging methods for tracking cell cycle transitions during human PSC differentiation that can be applied to study chromosome domain consolidation and other aspects of lineage specification.

Co-regulation of paralog genes in the three-dimensional chromatin architecture

Paralog genes arise from gene duplication events during evolution, which often lead to similar proteins that cooperate in common pathways and in protein complexes. Consequently, paralogs show correlation in gene expression whereby the mechanisms of co-regulation remain unclear. In eukaryotes, genes are regulated in part by distal enhancer elements through looping interactions with gene promoters. These looping interactions can be measured by genome-wide chromatin conformation capture (Hi-C) experiments, which revealed self-interacting regions called topologically associating domains (TADs). We hypothesize that paralogs share common regulatory mechanisms to enable coordinated expression according to TADs. To test this hypothesis, we integrated paralogy annotations with human gene expression data in diverse tissues, genome-wide enhancer-promoter associations and Hi-C experiments in human, mouse and dog genomes. We show that paralog gene pairs are enriched for co-localization in the same TAD, share more often common enhancer elements than expected and have increased contact frequencies over large genomic distances. Combined, our results indicate that paralogs share common regulatory mechanisms and cluster not only in the linear genome but also in the three-dimensional chromatin architecture. This enables concerted expression of paralogs over diverse cell-types and indicate evolutionary constraints in functional genome organization.

CTCF and Cohesin in Genome Folding and Transcriptional Gene Regulation

Genome function, replication, integrity, and propagation rely on the dynamic structural organization of chromosomes during the cell cycle. Genome folding in interphase provides regulatory segmentation for appropriate transcriptional control, facilitates ordered genome replication, and contributes to genome integrity by limiting illegitimate recombination. Here, we review recent high-resolution chromosome conformation capture and functional studies that have informed models of the spatial and regulatory compartmentalization of mammalian genomes, and discuss mechanistic models for how CTCF and cohesin control the functional architecture of mammalian chromosomes.

Controlling gene expression by DNA mechanics: emerging insights and challenges

Transcription initiation is a major control point for the precise regulation of gene expression. Our knowledge of this process has been mainly derived from protein-centric studies wherein cis-regulatory DNA sequences play a passive role, mainly in arranging the protein machinery to coalesce at the transcription start sites of genes in a spatial and temporal-specific manner. However, this is a highly dynamic process in which molecular motors such as RNA polymerase II (RNAPII), helicases, and other transcription factors, alter the level of mechanical force in DNA, rather than simply a set of static DNA-protein interactions. The double helix is a fiber that responds to flexural and torsional stress, which if accumulated, can affect promoter output as well as change DNA and chromatin structure. The relationship between DNA mechanics and the control of early transcription initiation events has been under-investigated. Genomic techniques to display topological stress and conformational variation in DNA across the mammalian genome provide an exciting new insight on the role of DNA mechanics in the early stages of the transcription cycle. Without understanding how torsional and flexural stresses are generated, transmitted, and dissipated, no model of transcription will be complete and accurate.

Genome maintenance in the context of 4D chromatin condensation

The eukaryotic genome is packaged in the three-dimensional nuclear space by forming loops, domains, and compartments in a hierarchical manner. However, when duplicated genomes prepare for segregation, mitotic cells eliminate topologically associating domains and abandon the compartmentalized structure. Alongside chromatin architecture reorganization during the transition from interphase to mitosis, cells halt most DNA-templated processes such as transcription and repair. The intrinsically condensed chromatin serves as a sophisticated signaling module subjected to selective relaxation for programmed genomic activities. To understand the elaborate genome-epigenome interplay during cell cycle progression, the steady three-dimensional genome requires a time scale to form a dynamic four-dimensional and a more comprehensive portrait. In this review, we will dissect the functions of critical chromatin architectural components in constructing and maintaining an orderly packaged chromatin environment. We will also highlight the importance of the spatially and temporally conscious orchestration of chromatin remodeling to ensure high-fidelity genetic transmission.

Selectivity of ORC binding sites and the relation to replication timing, fragile sites, and deletions in cancers

The origin recognition complex (ORC) binds sites from which DNA replication is initiated. We address ORC binding selectivity in vivo by mapping ∼52,000 ORC2 binding sites throughout the human genome. The ORC binding profile is broader than those of sequence-specific transcription factors, suggesting that ORC is not bound or recruited to specific DNA sequences. Instead, ORC binds nonspecifically to open (DNase I-hypersensitive) regions containing active chromatin marks such as H3 acetylation and H3K4 methylation. ORC sites in early and late replicating regions have similar properties, but there are far more ORC sites in early replicating regions. This suggests that replication timing is due primarily to ORC density and stochastic firing of origins. Computational simulation of stochastic firing from identified ORC sites is in accord with replication timing data. Large genomic regions with a paucity of ORC sites are strongly associated with common fragile sites and recurrent deletions in cancers. We suggest that replication origins, replication timing, and replication-dependent chromosome breaks are determined primarily by the genomic distribution of activator proteins at enhancers and promoters. These activators recruit nucleosome-modifying complexes to create the appropriate chromatin structure that allows ORC binding and subsequent origin firing.

Characterization of the replication timing program of 6 human model cell lines

During the S-phase, the DNA replication process is finely orchestrated and regulated by two programs: the spatial program that determines where replication will start in the genome (Cadoret et al. (2008 Oct 14), Cayrou et al. (2011 Sep), Picard et al. (2014 May 1) [1], [2], [3]), and the temporal program that determines when during the S phase different parts of the genome are replicated and when origins are activated. The temporal program is so well conserved for each cell type from independent individuals [4] that it is possible to identify a cell type from an unknown sample just by determining its replication timing program. Moreover, replicative domains are strongly correlated with the partition of the genome into topological domains (determined by the Hi-C method, Lieberman-Aiden et al. (2009 Oct 9), Pope et al. (2014 Nov 20) [5], [6]). On the one hand, replicative areas are well defined and participate in shaping the spatial organization of the genome for a given cell type. On the other hand, studies on the timing program during cell differentiation showed a certain plasticity of this program according to the stage of cell differentiation Hiratani et al. (2008 Oct 7, 2010 Feb) [7], [8]. Domains where a replication timing change was observed went through a nuclear re-localization. Thus the temporal program of replication can be considered as an epigenetic mark Hiratani and Gilbert (2009 Feb 16) [9]. We present the genomic data of replication timing in 6 human model cell lines: U2OS (GSM2111308), RKO (GSM2111309), HEK 293T (GSM2111310), HeLa (GSM2111311), MRC5-SV (GSM2111312) and K562 (GSM2111313). A short comparative analysis was performed that allowed us to define regions common to the 6 cell lines. These replication timing data can be taken into account when performing studies that use these model cell lines.

Long Neural Genes Harbor Recurrent DNA Break Clusters in Neural Stem/Progenitor Cells

Repair of DNA double-strand breaks (DSBs) by non-homologous end joining is critical for neural development, and brain cells frequently contain somatic genomic variations that might involve DSB intermediates. We now use an unbiased, high-throughput approach to identify genomic regions harboring recurrent DSBs in primary neural stem/progenitor cells (NSPCs). We identify 27 recurrent DSB clusters (RDCs), and remarkably, all occur within gene bodies. Most of these NSPC RDCs were detected only upon mild, aphidicolin-induced replication stress, providing a nucleotide-resolution view of replication-associated genomic fragile sites. The vast majority of RDCs occur in long, transcribed, and late-replicating genes. Moreover, almost 90% of identified RDC-containing genes are involved in synapse function and/or neural cell adhesion, with a substantial fraction also implicated in tumor suppression and/or mental disorders. Our characterization of NSPC RDCs reveals a basis of gene fragility and suggests potential impacts of DNA breaks on neurodevelopment and neural functions.

Gene regulation at a distance: From remote enhancers to 3D regulatory ensembles

Large-scale identification of elements associated with gene expression revealed that many of them are located extremely far from gene transcriptional start sites. We review here the growing evidence that show that distal cis-acting elements provide key instructions to genes, as genetic variation affecting them is growingly identified as an importance source of phenotypic diversity and disease. We discuss the different mechanisms that allow these elements to exert their regulatory functions, in a robust and specific manner, despite the large genomic distances separating them from their target genes. We particularly focus on the role of the structural organization of the genome in guiding such regulatory interactions.

CRISPR/Cas9 genome editing throws descriptive 3-D genome folding studies for a loop

CRISPR/Cas9 genome editing studies have recently shed new light into the causal link between the linear DNA sequence and 3-D chromatin architecture. Here we describe current models for the folding of genomes into a nested hierarchy of higher-order structures and discuss new insights into the organizing principles governing genome folding at each length scale. WIREs Syst Biol Med 2016, 8:286-299. doi: 10.1002/wsbm.1338 For further resources related to this article, please visit the WIREs website.

Mechanism of chromatin segregation to the nuclear periphery in C. elegans embryos

In eukaryotic organisms, gene regulation occurs in the context of chromatin. In the interphase nucleus, euchromatin and heterochromatin occupy distinct space during cell differentiation, with heterochromatin becoming enriched at the nuclear and nucleolar peripheries. This organization is thought to fine-tune gene expression. To elucidate the mechanisms that govern this level of genome organization, screens were carried out in C. elegans which monitored the loss of heterochromatin sequestration at the nuclear periphery. This led to the identification of a novel chromodomain protein, CEC-4 (Caenorhabditis elegans chromodomain protein 4) that mediates the anchoring of H3K9 methylation-bearing chromatin at the nuclear periphery in early to mid-stage embryos. Surprisingly, the loss of CEC-4 does not derepress genes found in heterochromatic domains, nor does it affect differentiation under standard laboratory conditions. On the other hand, CEC-4 contributes to the efficiency with which muscle differentiation is induced following ectopic expression of the master regulator, HLH-1. This is one of the first phenotypes specifically attributed to the ablation of heterochromatin anchoring.

The mutation spectrum in genomic late replication domains shapes mammalian GC content

Genome sequence compositions and epigenetic organizations are correlated extensively across multiple length scales. Replication dynamics, in particular, is highly correlated with GC content. We combine genome-wide time of replication (ToR) data, topological domains maps and detailed functional epigenetic annotations to study the correlations between replication timing and GC content at multiple scales. We find that the decrease in genomic GC content at large scale late replicating regions can be explained by mutation bias favoring A/T nucleotide, without selection or biased gene conversion. Quantification of the free dNTP pool during the cell cycle is consistent with a mechanism involving replication-coupled mutation spectrum that favors AT nucleotides at late S-phase. We suggest that mammalian GC content composition is shaped by independent forces, globally modulating mutation bias and locally selecting on functional element. Deconvoluting these forces and analyzing them on their native scales is important for proper characterization of complex genomic correlations.

The protozoan nucleus

The nucleus is arguably the defining characteristic of eukaryotes, distinguishing their cell organisation from both bacteria and archaea. Though the evolutionary history of the nucleus remains the subject of debate, its emergence differs from several other eukaryotic organelles in that it appears not to have evolved through symbiosis, but by cell membrane elaboration from an archaeal ancestor. Evolution of the nucleus has been accompanied by elaboration of nuclear structures that are intimately linked with most aspects of nuclear genome function, including chromosome organisation, DNA maintenance, replication and segregation, and gene expression controls. Here we discuss the complexity of the nucleus and its substructures in protozoan eukaryotes, with a particular emphasis on divergent aspects in eukaryotic parasites, which shed light on nuclear function throughout eukaryotes and reveal specialisations that underpin pathogen biology.

Spectral identification of topological domains

MOTIVATION

Topological domains have been proposed as the backbone of interphase chromosome structure. They are regions of high local contact frequency separated by sharp boundaries. Genes within a domain often have correlated transcription. In this paper, we present a computational efficient spectral algorithm to identify topological domains from chromosome conformation data (Hi-C data). We consider the genome as a weighted graph with vertices defined by loci on a chromosome and the edge weights given by interaction frequency between two loci. Laplacian-based graph segmentation is then applied iteratively to obtain the domains at the given compactness level. Comparison with algorithms in the literature shows the advantage of the proposed strategy.

RESULTS

An efficient algorithm is presented to identify topological domains from the Hi-C matrix.

AVAILABILITY AND IMPLEMENTATION

The Matlab source code and illustrative examples are available at http://bionetworks.ccmb.med.umich.edu/

CONTACT

: indikar@med.umich.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Replication timing and transcriptional control: beyond cause and effect-part III

DNA replication is essential for faithful transmission of genetic information and is intimately tied to chromosome structure and function. Genome duplication occurs in a defined temporal order known as the replication-timing (RT) program, which is regulated during the cell cycle and development in discrete units referred to as replication domains (RDs). RDs correspond to topologically-associating domains (TADs) and are spatio-temporally compartmentalized in the nucleus. While improvements in experimental tools have begun to reveal glimpses of causality, they have also unveiled complex context-dependent relationships that challenge long recognized correlations of RT to chromatin organization and gene regulation. In particular, RDs/TADs that switch RT during development march to the beat of a different drummer.

Influence of ATM-Mediated DNA Damage Response on Genomic Variation in Human Induced Pluripotent Stem Cells

Genome instability is a potential limitation to the research and therapeutic application of induced pluripotent stem cells (iPSCs). Observed genomic variations reflect the combined activities of DNA damage, cellular DNA damage response (DDR), and selection pressure in culture. To understand the contribution of DDR on the distribution of copy number variations (CNVs) in iPSCs, we mapped CNVs of iPSCs with mutations in the central DDR gene ATM onto genome organization landscapes defined by genome-wide replication timing profiles. We show that following reprogramming the early and late replicating genome is differentially affected by CNVs in ATM-deficient iPSCs relative to wild-type iPSCs. Specifically, the early replicating regions had increased CNV losses during retroviral (RV) reprogramming. This differential CNV distribution was not present after later passage or after episomal reprogramming. Comparison of different reprogramming methods in the setting of defective DDR reveals unique vulnerability of early replicating open chromatin to RV vectors.

A Tale of Two Cities: How Xist and its partners localize to and silence the bicompartmental X

Sex chromosomal dosage compensation in mammals takes the form of X chromosome inactivation (XCI), driven by the non-coding RNA Xist. In contrast to dosage compensation systems of flies and worms, mammalian XCI has to restrict its function to the Xist-producing X chromosome, while leaving autosomes and active X untouched. The mechanisms behind the long-range yet cis-specific localization and silencing activities of Xist have long been enigmatic, but genomics, proteomics, super-resolution microscopy, and innovative genetic approaches have produced significant new insights in recent years. In this review, I summarize and integrate these findings with a particular focus on the redundant yet mutually reinforcing pathways that enable long-term transcriptional repression throughout the soma. This includes an exploration of concurrent epigenetic changes acting in parallel within two distinct compartments of the inactive X. I also examine how Polycomb repressive complexes 1 and 2 and macroH2A may bridge XCI establishment and maintenance. XCI is a remarkable phenomenon that operates across multiple scales, combining changes in nuclear architecture, chromosome topology, chromatin compaction, and nucleosome/nucleotide-level epigenetic cues. Learning how these pathways act in concert likely holds the answer to the riddle posed by Cattanach's and other autosomal translocations: What makes the X especially receptive to XCI?

4D Visualization of replication foci in mammalian cells corresponding to individual replicons

Since the pioneering proposal of the replicon model of DNA replication 50 years ago, the predicted replicons have not been identified and quantified at the cellular level. Here, we combine conventional and super-resolution microscopy of replication sites in live and fixed cells with computational image analysis. We complement these data with genome size measurements, comprehensive analysis of S-phase dynamics and quantification of replication fork speed and replicon size in human and mouse cells. These multidimensional analyses demonstrate that replication foci (RFi) in three-dimensional (3D) preserved somatic mammalian cells can be optically resolved down to single replicons throughout S-phase. This challenges the conventional interpretation of nuclear RFi as replication factories, that is, the complex entities that process multiple clustered replicons. Accordingly, 3D genome organization and duplication can be now followed within the chromatin context at the level of individual replicons.

3D replicon distributions arise from stochastic initiation and domino-like DNA replication progression

DNA replication dynamics in cells from higher eukaryotes follows very complex but highly efficient mechanisms. However, the principles behind initiation of potential replication origins and emergence of typical patterns of nuclear replication sites remain unclear. Here, we propose a comprehensive model of DNA replication in human cells that is based on stochastic, proximity-induced replication initiation. Critical model features are: spontaneous stochastic firing of individual origins in euchromatin and facultative heterochromatin, inhibition of firing at distances below the size of chromatin loops and a domino-like effect by which replication forks induce firing of nearby origins. The model reproduces the empirical temporal and chromatin-related properties of DNA replication in human cells. We advance the one-dimensional DNA replication model to a spatial model by taking into account chromatin folding in the nucleus, and we are able to reproduce the spatial and temporal characteristics of the replication foci distribution throughout S-phase.

Replication origins: determinants or consequences of nuclear organization?

Chromosome replication, gene expression and chromatin assembly all occur on the same template, necessitating a tight spatial and temporal coordination to maintain genomic stability. The distribution of replication initiation events is responsive to local and global changes in chromatin structure and is affected by transcriptional activity. Concomitantly, replication origin sequences, which determine the locations of replication initiation events, can affect chromatin structure and modulate transcriptional efficiency. The flexibility observed in the replication initiation landscape might help achieve complete and accurate genome duplication while coordinating the DNA replication program with transcription and other nuclear processes in a cell-type specific manner. This review discusses the relationships among replication origin distribution, local and global chromatin structures and concomitant nuclear metabolic processes.

Atypical Teratoid/Rhabdoid Tumors Are Comprised of Three Epigenetic Subgroups with Distinct Enhancer Landscapes

Atypical teratoid/rhabdoid tumor (ATRT) is one of the most common brain tumors in infants. Although the prognosis of ATRT patients is poor, some patients respond favorably to current treatments, suggesting molecular inter-tumor heterogeneity. To investigate this further, we genetically and epigenetically analyzed 192 ATRTs. Three distinct molecular subgroups of ATRTs, associated with differences in demographics, tumor location, and type of SMARCB1 alterations, were identified. Whole-genome DNA and RNA sequencing found no recurrent mutations in addition to SMARCB1 that would explain the differences between subgroups. Whole-genome bisulfite sequencing and H3K27Ac chromatin-immunoprecipitation sequencing of primary tumors, however, revealed clear differences, leading to the identification of subgroup-specific regulatory networks and potential therapeutic targets.

How to rule the nucleus: divide et impera

Genome-wide molecular studies have provided new insights into the organization of nuclear chromatin by revealing the presence of chromatin domains of differing transcriptional activity, frequency of cis-interactions, proximity to scaffolding structures and replication timing. These studies have not only brought our understanding of genome function to a new level, but also offered functional insight for many phenomena observed in microscopic studies. In this review, we discuss the major principles of nuclear organization based on the spatial segregation of euchromatin and heterochromatin, as well as the dynamic genome rearrangements occurring during cell differentiation and development. We hope to unite the existing molecular and microscopic data on genome organization to get a holistic view of the nucleus, and propose a model, in which repeat repertoire together with scaffolding structures blueprint the functional nuclear architecture.

Spatially coordinated replication and minimization of expression noise constrain three-dimensional organization of yeast genome

Despite recent advances, the underlying functional constraints that shape the three-dimensional organization of eukaryotic genome are not entirely clear. Through comprehensive multivariate analyses of genome-wide datasets, we show that cis and trans interactions in yeast genome have significantly distinct functional associations. In particular, (i) the trans interactions are constrained by coordinated replication and co-varying mutation rates of early replicating domains through interactions among early origins, while cis interactions are constrained by coordination of late replication through interactions among late origins; (ii)cis and trans interactions exhibit differential preference for nucleosome occupancy; (iii)cis interactions are also constrained by the essentiality and co-fitness of interacting genes. Essential gene clusters associate with high average interaction frequency, relatively short-range interactions of low variance, and exhibit less fluctuations in chromatin conformation, marking a physically restrained state of engaged loci that, we suggest, is important to mitigate the epigenetic errors by restricting the spatial mobility of loci. Indeed, the genes with lower expression noise associate with relatively short-range interactions of lower variance and exhibit relatively higher average interaction frequency, a property that is conserved across Escherichia coli,yeast, and mESCs. Altogether, our observations highlight the coordination of replication and the minimization of expression noise, not necessarily co-expression of genes, as potent evolutionary constraints shaping the spatial organization of yeast genome.

Reorganization of chromosome architecture in replicative cellular senescence

Replicative cellular senescence is a fundamental biological process characterized by an irreversible arrest of proliferation. Senescent cells accumulate a variety of epigenetic changes, but the three-dimensional (3D) organization of their chromatin is not known. We applied a combination of whole-genome chromosome conformation capture (Hi-C), fluorescence in situ hybridization, and in silico modeling methods to characterize the 3D architecture of interphase chromosomes in proliferating, quiescent, and senescent cells. Although the overall organization of the chromatin into active (A) and repressive (B) compartments and topologically associated domains (TADs) is conserved between the three conditions, a subset of TADs switches between compartments. On a global level, the Hi-C interaction matrices of senescent cells are characterized by a relative loss of long-range and gain of short-range interactions within chromosomes. Direct measurements of distances between genetic loci, chromosome volumes, and chromatin accessibility suggest that the Hi-C interaction changes are caused by a significant reduction of the volumes occupied by individual chromosome arms. In contrast, centromeres oppose this overall compaction trend and increase in volume. The structural model arising from our study provides a unique high-resolution view of the complex chromosomal architecture in senescent cells.

Human and murine erythropoiesis

PURPOSE OF REVIEW

Research into the fundamental mechanisms of erythropoiesis has provided critical insights into inherited and acquired disorders of the erythrocyte. Studies of human erythropoiesis have primarily utilized in-vitro systems, whereas murine models have provided insights from in-vivo studies. This report reviews recent insights into human and murine erythropoiesis gained from transcriptome-based analyses.

RECENT FINDINGS

The availability of high-throughput genomic methodologies has allowed attainment of detailed gene expression data from cells at varying developmental and differentiation stages of erythropoiesis. Transcriptome analyses of human and murine reveal both stage and species-specific similarities and differences across terminal erythroid differentiation. Erythroid-specific long noncoding RNAs exhibit poor sequence conservation between human and mouse. Genome-wide analyses of alternative splicing reveal that complex, dynamic, stage-specific programs of alternative splicing program are utilized during terminal erythroid differentiation. Transcriptome data provide a significant resource for understanding mechanisms of normal and perturbed erythropoiesis. Understanding these processes will provide innovative strategies to detect, diagnose, prevent, and treat hematologic disease.

SUMMARY

Understanding the shared and different mechanisms controlling human and murine erythropoiesis will allow investigators to leverage the best model system to provide insights in normal and perturbed erythropoiesis.

Active medulloblastoma enhancers reveal subgroup-specific cellular origins

Medulloblastoma is a highly malignant paediatric brain tumour, often inflicting devastating consequences on the developing child. Genomic studies have revealed four distinct molecular subgroups with divergent biology and clinical behaviour. An understanding of the regulatory circuitry governing the transcriptional landscapes of medulloblastoma subgroups, and how this relates to their respective developmental origins, is lacking. Here, using H3K27ac and BRD4 chromatin immunoprecipitation followed by sequencing (ChIP-seq) coupled with tissue-matched DNA methylation and transcriptome data, we describe the active cis-regulatory landscape across 28 primary medulloblastoma specimens. Analysis of differentially regulated enhancers and super-enhancers reinforced inter-subgroup heterogeneity and revealed novel, clinically relevant insights into medulloblastoma biology. Computational reconstruction of core regulatory circuitry identified a master set of transcription factors, validated by ChIP-seq, that is responsible for subgroup divergence, and implicates candidate cells of origin for Group 4. Our integrated analysis of enhancer elements in a large series of primary tumour samples reveals insights into cis-regulatory architecture, unrecognized dependencies, and cellular origins.

Nuclear Architecture Organized by Rif1 Underpins the Replication-Timing Program

DNA replication is temporally and spatially organized in all eukaryotes, yet the molecular control and biological function of the replication-timing program are unclear. Rif1 is required for normal genome-wide regulation of replication timing, but its molecular function is poorly understood. Here we show that in mouse embryonic stem cells, Rif1 coats late-replicating domains and, with Lamin B1, identifies most of the late-replicating genome. Rif1 is an essential determinant of replication timing of non-Lamin B1-bound late domains. We further demonstrate that Rif1 defines and restricts the interactions between replication-timing domains during the G1 phase, thereby revealing a function of Rif1 as organizer of nuclear architecture. Rif1 loss affects both number and replication-timing specificity of the interactions between replication-timing domains. In addition, during the S phase, Rif1 ensures that replication of interacting domains is temporally coordinated. In summary, our study identifies Rif1 as the molecular link between nuclear architecture and replication-timing establishment in mammals.

Histones and histone modifications in perinuclear chromatin anchoring: from yeast to man

It is striking that within a eukaryotic nucleus, the genome can assume specific spatiotemporal distributions that correlate with the cell's functional states. Cell identity itself is determined by distinct sets of genes that are expressed at a given time. On the level of the individual gene, there is a strong correlation between transcriptional activity and associated histone modifications. Histone modifications act by influencing the recruitment of non-histone proteins and by determining the level of chromatin compaction, transcription factor binding, and transcription elongation. Accumulating evidence also shows that the subnuclear position of a gene or domain correlates with its expression status. Thus, the question arises whether this spatial organization results from or determines a gene's chromatin status. Although the association of a promoter with the inner nuclear membrane (INM) is neither necessary nor sufficient for repression, the perinuclear sequestration of heterochromatin is nonetheless conserved from yeast to man. How does subnuclear localization influence gene expression? Recent work argues that the common denominator between genome organization and gene expression is the modification of histones and in some cases of histone variants. This provides an important link between local chromatin structure and long-range genome organization in interphase cells. In this review, we will evaluate how histones contribute to the latter, and discuss how this might help to regulate genes crucial for cell differentiation.

Replication landscape of the human genome

Despite intense investigation, human replication origins and termini remain elusive. Existing data have shown strong discrepancies. Here we sequenced highly purified Okazaki fragments from two cell types and, for the first time, quantitated replication fork directionality and delineated initiation and termination zones genome-wide. Replication initiates stochastically, primarily within non-transcribed, broad (up to 150 kb) zones that often abut transcribed genes, and terminates dispersively between them. Replication fork progression is significantly co-oriented with the transcription. Initiation and termination zones are frequently contiguous, sometimes separated by regions of unidirectional replication. Initiation zones are enriched in open chromatin and enhancer marks, even when not flanked by genes, and often border ‘topologically associating domains' (TADs). Initiation zones are enriched in origin recognition complex (ORC)-binding sites and better align to origins previously mapped using bubble-trap than λ-exonuclease. This novel panorama of replication reveals how chromatin and transcription modulate the initiation process to create cell-type-specific replication programs.

The physical origin and termination sites of DNA replication in human cells have remained elusive. Here the authors use Okazaki fragment sequencing to reveal global replication patterns and show how chromatin and transcription modulate the process.

APOBEC-induced mutations in human cancers are strongly enriched on the lagging DNA strand during replication

APOBEC3A and APOBEC3B, cytidine deaminases of the APOBEC family, are among the main factors causing mutations in human cancers. APOBEC deaminates cytosines in single-stranded DNA (ssDNA). A fraction of the APOBEC-induced mutations occur as clusters ("kataegis") in single-stranded DNA produced during repair of double-stranded breaks (DSBs). However, the properties of the remaining 87% of nonclustered APOBEC-induced mutations, the source and the genomic distribution of the ssDNA where they occur, are largely unknown. By analyzing genomic and exomic cancer databases, we show that >33% of dispersed APOBEC-induced mutations occur on the lagging strand during DNA replication, thus unraveling the major source of ssDNA targeted by APOBEC in cancer. Although methylated cytosine is generally more mutation-prone than nonmethylated cytosine, we report that methylation reduces the rate of APOBEC-induced mutations by a factor of roughly two. Finally, we show that in cancers with extensive APOBEC-induced mutagenesis, there is almost no increase in mutation rates in late replicating regions (contrary to other cancers). Because late-replicating regions are depleted in exons, this results in a 1.3-fold higher fraction of mutations residing within exons in such cancers. This study provides novel insight into the APOBEC-induced mutagenesis and describes the peculiarity of the mutational processes in cancers with the signature of APOBEC-induced mutations.

A diffusion model for the coordination of DNA replication in Schizosaccharomyces pombe

The locations of proteins and epigenetic marks on the chromosomal DNA sequence are believed to demarcate the eukaryotic genome into distinct structural and functional domains that contribute to gene regulation and genome organization. However, how these proteins and epigenetic marks are organized in three dimensions remains unknown. Recent advances in proximity-ligation methodologies and high resolution microscopy have begun to expand our understanding of these spatial relationships. Here we use polymer models to examine the spatial organization of epigenetic marks, euchromatin and heterochromatin, and origins of replication within the Schizosaccharomyces pombe genome. These models incorporate data from microscopy and proximity-ligation experiments that inform on the positions of certain elements and contacts within and between chromosomes. Our results show a striking degree of compartmentalization of epigenetic and genomic features and lead to the proposal of a diffusion based mechanism, centred on the spindle pole body, for the coordination of DNA replication in S. pombe.

Hierarchical folding and reorganization of chromosomes are linked to transcriptional changes in cellular differentiation

Mammalian chromosomes fold into arrays of megabase‐sized topologically associating domains (TADs), which are arranged into compartments spanning multiple megabases of genomic DNA. TADs have internal substructures that are often cell type specific, but their higher‐order organization remains elusive. Here, we investigate TAD higher‐order interactions with Hi‐C through neuronal differentiation and show that they form a hierarchy of domains‐within‐domains (metaTADs) extending across genomic scales up to the range of entire chromosomes. We find that TAD interactions are well captured by tree‐like, hierarchical structures irrespective of cell type. metaTAD tree structures correlate with genetic, epigenomic and expression features, and structural tree rearrangements during differentiation are linked to transcriptional state changes. Using polymer modelling, we demonstrate that hierarchical folding promotes efficient chromatin packaging without the loss of contact specificity, highlighting a role far beyond the simple need for packing efficiency.

Implications of Epigenetic Variability within a Cell Population for "Cell Type" Classification

Here, we propose a new approach to defining nerve “cell types” in reaction to recent advances in single cell analysis. Among cells previously thought to be equivalent, considerable differences in global gene expression and biased tendencies among differing developmental fates have been demonstrated within multiple lineages. The model of classifying cells into distinct types thus has to be revised to account for this intrinsic variability. A “cell type” could be a group of cells that possess similar, but not necessarily identical properties, variable within a spectrum of epigenetic adjustments that permit its developmental path toward a specific function to be achieved. Thus, the definition of a cell type is becoming more similar to the definition of a species: sharing essential properties with other members of its group, but permitting a certain amount of deviation in aspects that do not seriously impact function. This approach accommodates, even embraces the spectrum of natural variation found in various cell populations and consequently avoids the fallacy of false equivalence. For example, developing neurons will react to their microenvironments with epigenetic changes resulting in slight changes in gene expression and morphology. Addressing the new questions implied here will have significant implications for developmental neurobiology.

Minor Loops in Major Folds: Enhancer-Promoter Looping, Chromatin Restructuring, and Their Association with Transcriptional Regulation and Disease

The organization and folding of chromatin within the nucleus can determine the outcome of gene expression. Recent technological advancements have enabled us to study chromatin interactions in a genome-wide manner at high resolution. These studies have increased our understanding of the hierarchy and dynamics of chromatin domains that facilitate cognate enhancer–promoter looping, defining the transcriptional program of different cell types. In this review, we focus on vertebrate chromatin long-range interactions as they relate to transcriptional regulation. In addition, we describe how the alteration of boundaries that mark discrete regions in the genome with high interaction frequencies within them, called topological associated domains (TADs), could lead to various phenotypes, including human diseases, which we term as “TADopathies.”

Glucocorticoid receptor and nuclear factor kappa-b affect three-dimensional chromatin organization

Background

The impact of signal-dependent transcription factors, such as glucocorticoid receptor and nuclear factor kappa-b, on the three-dimensional organization of chromatin remains a topic of discussion. The possible scenarios range from remodeling of higher order chromatin architecture by activated transcription factors to recruitment of activated transcription factors to pre-established long-range interactions.

Results

Using circular chromosome conformation capture coupled with next generation sequencing and high-resolution chromatin interaction analysis by paired-end tag sequencing of P300, we observed agonist-induced changes in long-range chromatin interactions, and uncovered interconnected enhancer–enhancer hubs spanning up to one megabase. The vast majority of activated glucocorticoid receptor and nuclear factor kappa-b appeared to join pre-existing P300 enhancer hubs without affecting the chromatin conformation. In contrast, binding of the activated transcription factors to loci with their consensus response elements led to the increased formation of an active epigenetic state of enhancers and a significant increase in long-range interactions within pre-existing enhancer networks. De novo enhancers or ligand-responsive enhancer hubs preferentially interacted with ligand-induced genes.

Conclusions

We demonstrate that, at a subset of genomic loci, ligand-mediated induction leads to active enhancer formation and an increase in long-range interactions, facilitating efficient regulation of target genes. Therefore, our data suggest an active role of signal-dependent transcription factors in chromatin and long-range interaction remodeling.

Electronic supplementary material

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

Cell cycle control in the early embryonic development of aquatic animal species

The cell cycle is integrated with many aspects of embryonic development. Not only is proper control over the pace of cell proliferation important, but also the timing of cell cycle progression is coordinated with transcription, cell migration, and cell differentiation. Due to the ease with which the embryos of aquatic organisms can be observed and manipulated, they have been a popular choice for embryologists throughout history. In the cell cycle field, aquatic organisms have been extremely important because they have played a major role in the discovery and analysis of key regulators of the cell cycle. In particular, the frog Xenopus laevis has been instrumental for understanding how the basic embryonic cell cycle is regulated. More recently, the zebrafish has been used to understand how the cell cycle is remodeled during vertebrate development and how it is regulated during morphogenesis. This review describes how some of the unique strengths of aquatic species have been leveraged for cell cycle research and suggests how species such as Xenopus and zebrafish will continue to reveal the roles of the cell cycle in human biology and disease.

BRPF3-HBO1 regulates replication origin activation and histone H3K14 acetylation

During DNA replication, thousands of replication origins are activated across the genome. Chromatin architecture contributes to origin specification and usage, yet it remains unclear which chromatin features impact on DNA replication. Here, we perform a RNAi screen for chromatin regulators implicated in replication control by measuring RPA accumulation upon replication stress. We identify six factors required for normal rates of DNA replication and characterize a function of the bromodomain and PHD finger-containing protein 3 (BRPF3) in replication initiation. BRPF3 forms a complex with HBO1 that specifically acetylates histone H3K14, and genomewide analysis shows high enrichment of BRPF3, HBO1 and H3K14ac at ORC1-binding sites and replication origins found in the vicinity of TSSs. Consistent with this, BRPF3 is necessary for H3K14ac at selected origins and efficient origin activation. CDC45 recruitment, but not MCM2-7 loading, is impaired in BRPF3-depleted cells, identifying a BRPF3-dependent function of HBO1 in origin activation that is complementary to its role in licencing. We thus propose that BRPF3-HBO1 acetylation of histone H3K14 around TSS facilitates efficient activation of nearby replication origins.

Chromosome Architecture and Genome Organization

How the same DNA sequences can function in the three-dimensional architecture of interphase nucleus, fold in the very compact structure of metaphase chromosomes and go precisely back to the original interphase architecture in the following cell cycle remains an unresolved question to this day. The strategy used to address this issue was to analyze the correlations between chromosome architecture and the compositional patterns of DNA sequences spanning a size range from a few hundreds to a few thousands Kilobases. This is a critical range that encompasses isochores, interphase chromatin domains and boundaries, and chromosomal bands. The solution rests on the following key points: 1) the transition from the looped domains and sub-domains of interphase chromatin to the 30-nm fiber loops of early prophase chromosomes goes through the unfolding into an extended chromatin structure (probably a 10-nm "beads-on-a-string" structure); 2) the architectural proteins of interphase chromatin, such as CTCF and cohesin sub-units, are retained in mitosis and are part of the discontinuous protein scaffold of mitotic chromosomes; 3) the conservation of the link between architectural proteins and their binding sites on DNA through the cell cycle explains the "mitotic memory" of interphase architecture and the reversibility of the interphase to mitosis process. The results presented here also lead to a general conclusion which concerns the existence of correlations between the isochore organization of the genome and the architecture of chromosomes from interphase to metaphase.

Large-Scale Chromatin Structure-Function Relationships during the Cell Cycle and Development: Insights from Replication Timing

Chromosome architecture has received a lot of attention since the recent development of genome-scale methods to measure chromatin interactions (Hi-C), enabling the first sequence-based models of chromosome tertiary structure. A view has emerged of chromosomes as a string of structural units (topologically associating domains; TADs) whose boundaries persist through the cell cycle and development. TADs with similar chromatin states tend to aggregate, forming spatially segregated chromatin compartments. However, high-resolution Hi-C has revealed substructure within TADs (subTADs) that poses a challenge for models that attribute significance to structural units at any given scale. More than 20 years ago, the DNA replication field independently identified stable structural (and functional) units of chromosomes (replication foci) as well as spatially segregated chromatin compartments (early and late foci), but lacked the means to link these units to genomic map units. Genome-wide studies of replication timing (RT) have now merged these two disciplines by identifying individual units of replication regulation (replication domains; RDs) that correspond to TADs and are arranged in 3D to form spatiotemporally segregated subnuclear compartments. Furthermore, classifying RDs/TADs by their constitutive versus developmentally regulated RT has revealed distinct classes of chromatin organization, providing unexpected insight into the relationship between large-scale chromosome structure and function.

Gene regulation during development in the light of topologically associating domains

During embryonic development, complex transcriptional programs govern the precision of gene expression. Many key developmental genes are regulated via cis-regulatory elements that are located far away in the linear genome. How sequences located hundreds of kilobases away from a promoter can influence its activity has been the subject of numerous speculations, which all underline the importance of the 3D-organization of the genome. The recent advent of chromosome conformation capture techniques has put into focus the subdivision of the genome into topologically associating domains (TADs). TADs may influence regulatory activities on multiple levels. The relative invariance of TAD limits across cell types suggests that they may form fixed structural domains that could facilitate and/or confine long-range regulatory interactions. However, most recent studies suggest that interactions within TADs are more variable and dynamic than initially described. Hence, different models are emerging regarding how TADs shape the complex 3D conformations, and thereafter influence the networks of cis-interactions that govern gene expression during development. For further resources related to this article, please visit the WIREs website.

De novo identification of replication-timing domains in the human genome by deep learning

Motivation: The de novo identification of the initiation and termination zones—regions that replicate earlier or later than their upstream and downstream neighbours, respectively—remains a key challenge in DNA replication.

Results: Building on advances in deep learning, we developed a novel hybrid architecture combining a pre-trained, deep neural network and a hidden Markov model (DNN-HMM) for the de novo identification of replication domains using replication timing profiles. Our results demonstrate that DNN-HMM can significantly outperform strong, discriminatively trained Gaussian mixture model–HMM (GMM-HMM) systems and other six reported methods that can be applied to this challenge. We applied our trained DNN-HMM to identify distinct replication domain types, namely the early replication domain (ERD), the down transition zone (DTZ), the late replication domain (LRD) and the up transition zone (UTZ), using newly replicated DNA sequencing (Repli-Seq) data across 15 human cells. A subsequent integrative analysis revealed that these replication domains harbour unique genomic and epigenetic patterns, transcriptional activity and higher-order chromosomal structure. Our findings support the ‘replication-domain’ model, which states (1) that ERDs and LRDs, connected by UTZs and DTZs, are spatially compartmentalized structural and functional units of higher-order chromosomal structure, (2) that the adjacent DTZ-UTZ pairs form chromatin loops and (3) that intra-interactions within ERDs and LRDs tend to be short-range and long-range, respectively. Our model reveals an important chromatin organizational principle of the human genome and represents a critical step towards understanding the mechanisms regulating replication timing.

Availability and implementation: Our DNN-HMM method and three additional algorithms can be freely accessed at https://github.com/wenjiegroup/DNN-HMM. The replication domain regions identified in this study are available in GEO under the accession ID GSE53984.

Contact:shuwj@bmi.ac.cn or boxc@bmi.ac.cn

Supplementary information:Supplementary data are available at Bioinformatics online.

Systematic Integration of Brain eQTL and GWAS Identifies ZNF323 as a Novel Schizophrenia Risk Gene and Suggests Recent Positive Selection Based on Compensatory Advantage on Pulmonary Function

Genome-wide association studies have identified multiple risk variants and loci that show robust association with schizophrenia. Nevertheless, it remains unclear how these variants confer risk to schizophrenia. In addition, the driving force that maintains the schizophrenia risk variants in human gene pool is poorly understood. To investigate whether expression-associated genetic variants contribute to schizophrenia susceptibility, we systematically integrated brain expression quantitative trait loci and genome-wide association data of schizophrenia using Sherlock, a Bayesian statistical framework. Our analyses identified ZNF323 as a schizophrenia risk gene (P = 2.22×10(-6)). Subsequent analyses confirmed the association of the ZNF323 and its expression-associated single nucleotide polymorphism rs1150711 in independent samples (gene-expression: P = 1.40×10(-6); single-marker meta-analysis in the combined discovery and replication sample comprising 44123 individuals: P = 6.85×10(-10)). We found that the ZNF323 was significantly downregulated in hippocampus and frontal cortex of schizophrenia patients (P = .0038 and P = .0233, respectively). Evidence for pleiotropic effects was detected (association of rs1150711 with lung function and gene expression of ZNF323 in lung: P = 6.62×10(-5) and P = 9.00×10(-5), respectively) with the risk allele (T allele) for schizophrenia acting as protective allele for lung function. Subsequent population genetics analyses suggest that the risk allele (T) of rs1150711 might have undergone recent positive selection in human population. Our findings suggest that the ZNF323 is a schizophrenia susceptibility gene whose expression may influence schizophrenia risk. Our study also illustrates a possible mechanism for maintaining schizophrenia risk variants in the human gene pool.

Active chromatin and transcription play a key role in chromosome partitioning into topologically associating domains

Recent advances enabled by the Hi-C technique have unraveled many principles of chromosomal folding that were subsequently linked to disease and gene regulation. In particular, Hi-C revealed that chromosomes of animals are organized into topologically associating domains (TADs), evolutionary conserved compact chromatin domains that influence gene expression. Mechanisms that underlie partitioning of the genome into TADs remain poorly understood. To explore principles of TAD folding in Drosophila melanogaster, we performed Hi-C and poly(A)⁺ RNA-seq in four cell lines of various origins (S2, Kc167, DmBG3-c2, and OSC). Contrary to previous studies, we find that regions between TADs (i.e., the inter-TADs and TAD boundaries) in Drosophila are only weakly enriched with the insulator protein dCTCF, while another insulator protein Su(Hw) is preferentially present within TADs. However, Drosophila inter-TADs harbor active chromatin and constitutively transcribed (housekeeping) genes. Accordingly, we find that binding of insulator proteins dCTCF and Su(Hw) predicts TAD boundaries much worse than active chromatin marks do. Interestingly, inter-TADs correspond to decompacted inter-bands of polytene chromosomes, whereas TADs mostly correspond to densely packed bands. Collectively, our results suggest that TADs are condensed chromatin domains depleted in active chromatin marks, separated by regions of active chromatin. We propose the mechanism of TAD self-assembly based on the ability of nucleosomes from inactive chromatin to aggregate, and lack of this ability in acetylated nucleosomal arrays. Finally, we test this hypothesis by polymer simulations and find that TAD partitioning may be explained by different modes of inter-nucleosomal interactions for active and inactive chromatin.

Chromosomal Loop Domains Direct the Recombination of Antigen Receptor Genes

RAG initiates antibody V(D)J recombination in developing lymphocytes by generating "on-target" DNA breaks at matched pairs of bona fide recombination signal sequences (RSSs). We employ bait RAG-generated breaks in endogenous or ectopically inserted RSS pairs to identify huge numbers of RAG "off-target" breaks. Such breaks occur at the simple CAC motif that defines the RSS cleavage site and are largely confined within convergent CTCF-binding element (CBE)-flanked loop domains containing bait RSS pairs. Marked orientation dependence of RAG off-target activity within loops spanning up to 2 megabases implies involvement of linear tracking. In this regard, major RAG off-targets in chromosomal translocations occur as convergent RSS pairs at enhancers within a loop. Finally, deletion of a CBE-based IgH locus element disrupts V(D)J recombination domains and, correspondingly, alters RAG on- and off-target distributions within IgH. Our findings reveal how RAG activity is developmentally focused and implicate mechanisms by which chromatin domains harness biological processes within them.