Topologically associating domains are stable units of replication-timing regulation

Nuclear Envelope Regulation of Oncogenic Processes: Roles in Pancreatic Cancer

Pancreatic cancer is an aggressive and intractable malignancy with high mortality. This is due in part to a high resistance to chemotherapeutics and radiation treatment conferred by diverse regulatory mechanisms. Among these, constituents of the nuclear envelope play a significant role in regulating oncogenesis and pancreatic tumor biology, and this review focuses on three specific components and their roles in cancer. The LINC complex is a nuclear envelope component formed by proteins with SUN and KASH domains that interact in the periplasmic space of the nuclear envelope. These interactions functionally and structurally couple the cytoskeleton to chromatin and facilitates gene regulation informed by cytoplasmic activity. Furthermore, cancer cell invasiveness is impacted by LINC complex biology. The nuclear lamina is adjacent to the inner nuclear membrane of the nuclear envelope and can actively regulate chromatin in addition to providing structural integrity to the nucleus. A disrupted lamina can impart biophysical compromise to nuclear structure and function, as well as form dysfunctional micronuclei that may lead to genomic instability and chromothripsis. In close relationship to the nuclear lamina is the nuclear pore complex, a large megadalton structure that spans both outer and inner membranes of the nuclear envelope. The nuclear pore complex mediates bidirectional nucleocytoplasmic transport and is comprised of specialized proteins called nucleoporins that are overexpressed in many cancers and are diagnostic markers for oncogenesis. Furthermore, recent demonstration of gene regulatory functions for discrete nucleoporins independent of their nuclear trafficking function suggests that these proteins may contribute more to malignant phenotypes beyond serving as biomarkers. The nuclear envelope is thus a complex, intricate regulator of cell signaling, with roles in pancreatic tumorigenesis and general oncogenic transformation.

Conflicts of CpG density and DNA methylation are proximally and distally involved in gene regulation in human and mouse tissues

The relationship between CpG content and DNA methylation has attracted considerable interest in recent years. Direct or indirect methods have been developed to investigate their regulatory functions based on various hypotheses, large cohort studies, and meta-analyses. However, all of these analyses were performed at units of CpG blocks and, thus, the influence of finer genome structure has been neglected. Herein, we present a novel algorithm of base-pair resolution to systematically investigate the relationship between CpG contents and DNA methylation. By introducing the concept of 'complementary index' we examined the methylomes of 34 adult and 7 embryonic tissues and successfully fitted the relationship of DNA methylation and CpG density into a nonlinear mathematical model. A further algorithm was developed to locate the regions where CpG density does not match expectations from the model, termed 'conflict of gap' (COG) regions. Interestingly, COGs are highly concordant in human and mouse and their distributions display a tissue-specific pattern. Based on COG methylation patterns we correctly classified tissues according to their function or origin. We demonstrate that COGs based on our method can reveal more and deeper information than traditional differential methylation region (DMR) approaches. We also found that when COGs are located near to transcription start site (TSS), these regions can determine which promoters will be utilized for initiating gene transcription. Furthermore, COGs located far from the TSS perform as enhancers in terms of histone modification, sequence conservation, transcription factor binding, and DNase I-hypersensitivity.

DNA double-strand breaks as drivers of neural genomic change, function, and disease

Early work from about two decades ago implicated DNA double-strand break (DSB) formation and repair in neuronal development. Findings emerging from recent studies of DSBs in proliferating neural progenitors and in mature, non-dividing neurons suggest important roles of DSBs in brain physiology, aging, cancer, psychiatric and neurodegenerative disorders. We provide an overview of some findings and speculate on what may lie ahead.

Dynamic changes in ORC localization and replication fork progression during tissue differentiation

Background

Genomic regions repressed for DNA replication, resulting in either delayed replication in S phase or underreplication in polyploid cells, are thought to be controlled by inhibition of replication origin activation. Studies in Drosophila polytene cells, however, raised the possibility that impeding replication fork progression also plays a major role.

Results

We exploited genomic regions underreplicated (URs) with tissue specificity in Drosophila polytene cells to analyze mechanisms of replication repression. By localizing the Origin Recognition Complex (ORC) in the genome of the larval fat body and comparing this to ORC binding in the salivary gland, we found that sites of ORC binding show extensive tissue specificity. In contrast, there are common domains nearly devoid of ORC in the salivary gland and fat body that also have reduced density of ORC binding sites in diploid cells. Strikingly, domains lacking ORC can still be replicated in some polytene tissues, showing absence of ORC and origins is insufficient to repress replication. Analysis of the width and location of the URs with respect to ORC position indicates that whether or not a genomic region lacking ORC is replicated is controlled by whether replication forks formed outside the region are inhibited.

Conclusions

These studies demonstrate that inhibition of replication fork progression can block replication across genomic regions that constitutively lack ORC. Replication fork progression can be inhibited in both tissue-specific and genome region-specific ways. Consequently, when evaluating sources of genome instability it is important to consider altered control of replication forks in response to differentiation.

Electronic supplementary material

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

Continuous-Trait Probabilistic Model for Comparing Multi-species Functional Genomic Data

A large amount of multi-species functional genomic data from high-throughput assays are becoming available to help understand the molecular mechanisms for phenotypic diversity across species. However, continuous-trait probabilistic models, which are key to such comparative analysis, remain under-explored. Here we develop a new model, called phylogenetic hidden Markov Gaussian processes (Phylo-HMGP), to simultaneously infer heterogeneous evolutionary states of functional genomic features in a genome-wide manner. Both simulation studies and real data application demonstrate the effectiveness of Phylo-HMGP. Importantly, we applied Phylo-HMGP to analyze a new cross-species DNA replication timing (RT) dataset from the same cell type in five primate species (human, chimpanzee, orangutan, gibbon, and green monkey). We demonstrate that our Phylo-HMGP model enables discovery of genomic regions with distinct evolutionary patterns of RT. Our method provides a generic framework for comparative analysis of multi-species continuous functional genomic signals to help reveal regions with conserved or lineage-specific regulatory roles.

Cellular senescence induces replication stress with almost no affect on DNA replication timing

Organismal aging entails a gradual decline of normal physiological functions and a major contributor to this decline is withdrawal of the cell cycle, known as senescence. Senescence can result from telomere diminution leading to a finite number of population doublings, known as replicative senescence (RS), or from oncogene overexpression, as a protective mechanism against cancer. Senescence is associated with large-scale chromatin re-organization and changes in gene expression. Replication stress is a complex phenomenon, defined as the slowing or stalling of replication fork progression and/or DNA synthesis, which has serious implications for genome stability, and consequently in human diseases. Aberrant replication fork structures activate the replication stress response leading to the activation of dormant origins, which is thought to be a safeguard mechanism to complete DNA replication on time. However, the relationship between replicative stress and the changes in the spatiotemporal program of DNA replication in senescence progression remains unclear. Here, we studied the DNA replication program during senescence progression in proliferative and pre-senescent cells from donors of various ages by single DNA fiber combing of replicated DNA, origin mapping by sequencing short nascent strands and genome-wide profiling of replication timing (TRT). We demonstrate that, progression into RS leads to reduced replication fork rates and activation of dormant origins, which are the hallmarks of replication stress. However, with the exception of a delay in RT of the CREB5 gene in all pre-senescent cells, RT was globally unaffected by replication stress during entry into either oncogene-induced or RS. Consequently, we conclude that RT alterations associated with physiological and accelerated aging, do not result from senescence progression. Our results clarify the interplay between senescence, aging and replication programs and demonstrate that RT is largely resistant to replication stress.

MCM2 promotes symmetric inheritance of modified histones during DNA replication

During genome replication, parental histones are recycled to newly replicated DNA with their posttranslational modifications (PTMs). Whether sister chromatids inherit modified histones evenly remains unknown. We measured histone PTM partition to sister chromatids in embryonic stem cells. We found that parental histones H3-H4 segregate to both daughter DNA strands with a weak leading-strand bias, skewing partition at topologically associating domain (TAD) borders and enhancers proximal to replication initiation zones. Segregation of parental histones to the leading strand increased markedly in cells with histone-binding mutations in MCM2, part of the replicative helicase, exacerbating histone PTM sister chromatid asymmetry. This work reveals how histones are inherited to sister chromatids and identifies a mechanism by which the replication machinery ensures symmetric cell division.

Decoding topologically associating domains with ultra-low resolution Hi-C data by graph structural entropy

Submegabase-size topologically associating domains (TAD) have been observed in high-throughput chromatin interaction data (Hi-C). However, accurate detection of TADs depends on ultra-deep sequencing and sophisticated normalization procedures. Here we propose a fast and normalization-free method to decode the domains of chromosomes (deDoc) that utilizes structural information theory. By treating Hi-C contact matrix as a representation of a graph, deDoc partitions the graph into segments with minimal structural entropy. We show that structural entropy can also be used to determine the proper bin size of the Hi-C data. By applying deDoc to pooled Hi-C data from 10 single cells, we detect megabase-size TAD-like domains. This result implies that the modular structure of the genome spatial organization may be fundamental to even a small cohort of single cells. Our algorithms may facilitate systematic investigations of chromosomal domains on a larger scale than hitherto have been possible.

Modeling three-dimensional genomic organization in evolution and pathogenesis

The regulation of gene expression is mediated via the complex three-dimensional (3D) conformation of the genetic material and its interactions with various intracellular factors. Various experimental and computational approaches have been developed in recent years for understating the relation between the 3D conformation of the genome and the phenotypes of cells in normal condition and diseases. In this review, we will discuss novel approaches for analyzing and modeling the 3D genomic conformation, focusing on deciphering disease-causing mutations that affect gene expression. We conclude that as this is a very challenging mission, an important direction should involve the comparative analysis of various 3D models from various organisms or cells.

The 3D genome: From fundamental principles to disease and cancer

In higher eukaryotes, the three-dimensional (3D) organization of the genome is intimately related to numerous key biological functions including gene expression, DNA repair and DNA replication regulations. Alteration of this 3D organization is detrimental to the organism and can give rise to a broad range of diseases such as cancers. Here, we review recent advances in the field. We first describe how the genome is packed in 3D to form chromosome territories, compartments and domains. We also give an overview of the recent techniques that allow to map the genome in 3D up to the kilobase resolution. We then discuss potential mechanisms by which genome misfolding can affect proper gene expression by distal enhancers, and how the 3D genome influences the formation of genomic rearrangements.

Pre-established Chromatin Interactions Mediate the Genomic Response to Glucocorticoids

The glucocorticoid receptor (GR) is a hormone-inducible transcription factor involved in metabolic and anti-inflammatory gene expression responses. To investigate what controls interactions between GR binding sites and their target genes, we used in situ Hi-C to generate high-resolution, genome-wide maps of chromatin interactions before and after glucocorticoid treatment. We found that GR binding to the genome typically does not cause new chromatin interactions to target genes but instead acts through chromatin interactions that already exist prior to hormone treatment. Both glucocorticoid-induced and glucocorticoid-repressed genes increased interactions with distal GR binding sites. In addition, while glucocorticoid-induced genes increased interactions with transcriptionally active chromosome compartments, glucocorticoid-repressed genes increased interactions with transcriptionally silent compartments. Lastly, while the architectural DNA-binding proteins CTCF and RAD21 were bound to most chromatin interactions, we found that glucocorticoid-responsive chromatin interactions were depleted for CTCF binding but enriched for RAD21. Together, these findings offer new insights into the mechanisms underlying GC-mediated gene activation and repression.

Genetic predisposition to B-cell acute lymphoblastic leukemia at 14q11.2 is mediated by a CEBPE promoter polymorphism

Acute lymphoblastic leukaemia (ALL) is the most common paediatric malignancy. Genome-wide association studies have shown variation at 14q11.2 influences ALL risk. We sought to decipher causal variant(s) at 14q11.2 and the mechanism of tumorigenesis. We show rs2239630 G>A resides in the promoter of the CCAT enhancer-binding protein epsilon (CEBPE) gene. The rs2239630-A risk allele is associated with increased promotor activity and CEBPE expression. Depletion of CEBPE in ALL cells reduces cell growth, correspondingly CEBPE binds to the promoters of electron transport and energy generation genes. RNA-seq in CEBPE depleted cells demonstrates CEBPE regulates the expression of genes involved in B-cell development (IL7R), apoptosis (BCL2), and methotrexate resistance (RASS4L). CEBPE regulated genes significantly overlapped in CEBPE depleted cells, ALL blasts and IGH-CEBPE translocated ALL. This suggests CEBPE regulates a similar set of genes in each, consistent with a common biological mechanism of leukemogenesis for rs2239630 associated and CEBPE translocated ALL. Finally, we map IGH-CEBPE translocation breakpoints in two cases, implicating RAG recombinase activity in their formation.

Poly(ADP-ribosyl)ation associated changes in CTCF-chromatin binding and gene expression in breast cells

CTCF is an evolutionarily conserved and ubiquitously expressed architectural protein regulating a plethora of cellular functions via different molecular mechanisms. CTCF can undergo a number of post-translational modifications which change its properties and functions. One such modifications linked to cancer is poly(ADP-ribosyl)ation (PARylation). The highly PARylated CTCF form has an apparent molecular mass of 180 kDa (referred to as CTCF180), which can be distinguished from hypo- and non-PARylated CTCF with the apparent molecular mass of 130 kDa (referred to as CTCF130). The existing data accumulated so far have been mainly related to CTCF130. However, the properties of CTCF180 are not well understood despite its abundance in a number of primary tissues. In this study we performed ChIP-seq and RNA-seq analyses in human breast cells 226LDM, which display predominantly CTCF130 when proliferating, but CTCF180 upon cell cycle arrest. We observed that in the arrested cells the majority of sites lost CTCF, whereas fewer sites gained CTCF or remain bound (i.e. common sites). The classical CTCF binding motif was found in the lost and common, but not in the gained sites. The changes in CTCF occupancies in the lost and common sites were associated with increased chromatin densities and altered expression from the neighboring genes. Based on these results we propose a model integrating the CTCF130/180 transition with CTCF-DNA binding and gene expression changes. This study also issues an important cautionary note concerning the design and interpretation of any experiments using cells and tissues where CTCF180 may be present.

Epigenetic maintenance of topological domains in the highly rearranged gibbon genome

The relationship between evolutionary genome remodeling and the three-dimensional structure of the genome remain largely unexplored. Here, we use the heavily rearranged gibbon genome to examine how evolutionary chromosomal rearrangements impact genome-wide chromatin interactions, topologically associating domains (TADs), and their epigenetic landscape. We use high-resolution maps of gibbon-human breaks of synteny (BOS), apply Hi-C in gibbon, measure an array of epigenetic features, and perform cross-species comparisons. We find that gibbon rearrangements occur at TAD boundaries, independent of the parameters used to identify TADs. This overlap is supported by a remarkable genetic and epigenetic similarity between BOS and TAD boundaries, namely presence of CpG islands and SINE elements, and enrichment in CTCF and H3K4me3 binding. Cross-species comparisons reveal that regions orthologous to BOS also correspond with boundaries of large (400-600 kb) TADs in human and other mammalian species. The colocalization of rearrangement breakpoints and TAD boundaries may be due to higher chromatin fragility at these locations and/or increased selective pressure against rearrangements that disrupt TAD integrity. We also examine the small portion of BOS that did not overlap with TAD boundaries and gave rise to novel TADs in the gibbon genome. We postulate that these new TADs generally lack deleterious consequences. Last, we show that limited epigenetic homogenization occurs across breakpoints, irrespective of their time of occurrence in the gibbon lineage. Overall, our findings demonstrate remarkable conservation of chromatin interactions and epigenetic landscape in gibbons, in spite of extensive genomic shuffling.

Promoter Capture Hi-C: High-resolution, Genome-wide Profiling of Promoter Interactions

The three-dimensional organization of the genome is linked to its function. For example, regulatory elements such as transcriptional enhancers control the spatio-temporal expression of their target genes through physical contact, often bridging considerable (in some cases hundreds of kilobases) genomic distances and bypassing nearby genes. The human genome harbors an estimated one million enhancers, the vast majority of which have unknown gene targets. Assigning distal regulatory regions to their target genes is thus crucial to understand gene expression control. We developed Promoter Capture Hi-C (PCHi-C) to enable the genome-wide detection of distal promoter-interacting regions (PIRs), for all promoters in a single experiment. In PCHi-C, highly complex Hi-C libraries are specifically enriched for promoter sequences through in-solution hybrid selection with thousands of biotinylated RNA baits complementary to the ends of all promoter-containing restriction fragments. The aim is to then pull-down promoter sequences and their frequent interaction partners such as enhancers and other potential regulatory elements. After high-throughput paired-end sequencing, a statistical test is applied to each promoter-ligated restriction fragment to identify significant PIRs at the restriction fragment level. We have used PCHi-C to generate an atlas of long-range promoter interactions in dozens of human and mouse cell types. These promoter interactome maps have contributed to a greater understanding of mammalian gene expression control by assigning putative regulatory regions to their target genes and revealing preferential spatial promoter-promoter interaction networks. This information also has high relevance to understanding human genetic disease and the identification of potential disease genes, by linking non-coding disease-associated sequence variants in or near control sequences to their target genes.

The eukaryotic bell-shaped temporal rate of DNA replication origin firing emanates from a balance between origin activation and passivation

The time-dependent rate I(t) of origin firing per length of unreplicated DNA presents a universal bell shape in eukaryotes that has been interpreted as the result of a complex time-evolving interaction between origins and limiting firing factors. Here, we show that a normal diffusion of replication fork components towards localized potential replication origins (p-oris) can more simply account for the I(t) universal bell shape, as a consequence of a competition between the origin firing time and the time needed to replicate DNA separating two neighboring p-oris. We predict the I(t) maximal value to be the product of the replication fork speed with the squared p-ori density. We show that this relation is robustly observed in simulations and in experimental data for several eukaryotes. Our work underlines that fork-component recycling and potential origins localization are sufficient spatial ingredients to explain the universality of DNA replication kinetics.

Before a cell can divide, it must duplicate its DNA. In eukaryotes – organisms such as animals and fungi, which store their DNA in the cell’s nucleus – DNA replication starts at specific sites in the genome called replication origins. At each origin sits a protein complex that will activate when it randomly captures an activating protein that diffuses within the nucleus. Once a replication origin activates or “fires”, the complex then splits into two new complexes that move away from each other as they duplicate the DNA. If an active complex collides with an inactive one at another origin, the latter is inactivated – a phenomenon known as origin passivation. When two active complexes meet, they release the activating proteins, which diffuse away and eventually activate other origins in unreplicated DNA.

The number of origins that activate each minute divided by the length of unreplicated DNA is referred to as the “rate of origin firing”. In all eukaryotes, this rate – also known as I(t) – follows the same pattern. First, it increases until more than half of the DNA is duplicated. Then it decreases until everything is duplicated. This means that, if plotted out, the graph of origin firing rate would always be a bell-shaped curve, even for organisms with genomes of different sizes that have different numbers of origins. The reason for this universal shape remained unclear.

Scientists had tried to create numerical simulations that model the rate of origin firing. However, for these simulations to reproduce the bell-shape curve, a number of untested assumptions had to be made about how DNA replication takes place. In addition, these models ignored the fact that it takes time to replicate the DNA between origins.

To take this time into account, Arbona et al. instead decided to model the replication origins as discrete and distinct entities. This way of building the mathematical model succeeded in reproducing the universal bell curve shape without additional assumptions. With this simulation, the balance between origin activation and passivation is enough to achieve the observed pattern.

The new model also predicts that the maximum rate of origin firing is determined by the speed of DNA replication and the density of origins in the genome. Arbona et al. verified this prediction in yeast, fly, frog and human cells – organisms with different sized genomes that take between 20 minutes and 8 hours to replicate their DNA. Lastly, the prediction also held true in yeast treated with hydroxyurea, an anticancer drug that slows DNA replication.

A better understanding of DNA replication can help scientists to understand how this process is perturbed in cancers and how drugs that target DNA replication can treat these diseases. Future work will explore how the 3D organization of the genome affects the diffusion of activating proteins within the cell nucleus.

5C-ID: Increased resolution Chromosome-Conformation-Capture-Carbon-Copy with in situ 3C and double alternating primer design

Mammalian genomes are folded in a hierarchy of compartments, topologically associating domains (TADs), subTADs, and looping interactions. Currently, there is a great need to evaluate the link between chromatin topology and genome function across many biological conditions and genetic perturbations. Hi-C can generate genome-wide maps of looping interactions but is intractable for high-throughput comparison of loops across multiple conditions due to the enormous number of reads (>6 Billion) required per library. Here, we describe 5C-ID, a new version of Chromosome-Conformation-Capture-Carbon-Copy (5C) with restriction digest and ligation performed in the nucleus (in situ Chromosome-Conformation-Capture (3C)) and ligation-mediated amplification performed with a double alternating primer design. We demonstrate that 5C-ID produces higher-resolution 3D genome folding maps with reduced spatial noise using markedly lower cell numbers than canonical 5C. 5C-ID enables the creation of high-resolution, high-coverage maps of chromatin loops in up to a 30 Megabase subset of the genome at a fraction of the cost of Hi-C.

Replication timing and nuclear structure

DNA replication proceeds along spatially and temporally coordinated patterns within the nucleus, thus protecting the genome during the synthesis of new genetic material. While we have been able to visualize replication patterns on DNA fibers for 50 years, recent developments and discoveries have provided a greater insight into how DNA replication is controlled. In this review, we highlight many of these discoveries. Of great interest are the physiological role of the replication timing program, cis and trans-acting factors that modulate replication timing and the effects of chromatin structure on the replication timing program. We also discuss future directions in the study of replication timing.

Allele-specific control of replication timing and genome organization during development

DNA replication occurs in a defined temporal order known as the replication-timing (RT) program. RT is regulated during development in discrete chromosomal units, coordinated with transcriptional activity and 3D genome organization. Here, we derived distinct cell types from F1 hybrid musculus × castaneus mouse crosses and exploited the high single-nucleotide polymorphism (SNP) density to characterize allelic differences in RT (Repli-seq), genome organization (Hi-C and promoter-capture Hi-C), gene expression (total nuclear RNA-seq), and chromatin accessibility (ATAC-seq). We also present HARP, a new computational tool for sorting SNPs in phased genomes to efficiently measure allele-specific genome-wide data. Analysis of six different hybrid mESC clones with different genomes (C57BL/6, 129/sv, and CAST/Ei), parental configurations, and gender revealed significant RT asynchrony between alleles across ∼12% of the autosomal genome linked to subspecies genomes but not to parental origin, growth conditions, or gender. RT asynchrony in mESCs strongly correlated with changes in Hi-C compartments between alleles but not as strongly with SNP density, gene expression, imprinting, or chromatin accessibility. We then tracked mESC RT asynchronous regions during development by analyzing differentiated cell types, including extraembryonic endoderm stem (XEN) cells, four male and female primary mouse embryonic fibroblasts (MEFs), and neural precursor cells (NPCs) differentiated in vitro from mESCs with opposite parental configurations. We found that RT asynchrony and allelic discordance in Hi-C compartments seen in mESCs were largely lost in all differentiated cell types, accompanied by novel sites of allelic asynchrony at a considerably smaller proportion of the genome, suggesting that genome organization of homologs converges to similar folding patterns during cell fate commitment.

An APOE-independent cis-eSNP on chromosome 19q13.32 influences tau levels and late-onset Alzheimer's disease risk

Although multiple susceptibility loci for late-onset Alzheimer's disease (LOAD) have been identified, a large portion of the genetic risk for this disease remains unexplained. LOAD risk may be associated with single-nucleotide polymorphisms responsible for changes in gene expression (eSNPs). To detect eSNPs associated with LOAD, we integrated data from LOAD genome-wide association studies and expression quantitative trait loci using Sherlock (a Bayesian statistical method). We identified a cis-regulatory eSNP (rs2927438) located on chromosome 19q13.32, for which subsequent analyses confirmed the association with both LOAD risk and the expression level of several nearby genes. Importantly, rs2927438 may represent an APOE-independent LOAD eSNP according to the weak linkage disequilibrium of rs2927438 with the 2 polymorphisms (rs7412 and rs429358) defining the APOE-ε2, -ε3, and -ε4 alleles. Furthermore, rs2927438 does not influence chromatin interaction events at the APOE locus or cis-regulation of APOE expression. Further exploratory analysis revealed that rs2927438 is significantly associated with tau levels in the cerebrospinal fluid. Our findings suggest that rs2927438 may confer APOE-independent risk for LOAD.

Super-resolution Imaging of Individual Human Subchromosomal Regions in Situ Reveals Nanoscopic Building Blocks of Higher-Order Structure

It is widely recognized that the higher-order spatial organization of the genome, beyond the nucleosome, plays an important role in many biological processes. However, to date, direct information on even such fundamental structural details as the typical sizes and DNA content of these higher-order structures in situ is poorly characterized. Here, we examine the nanoscopic DNA organization within human nuclei using super-resolution direct stochastic optical reconstruction microscopy (dSTORM) imaging and 5-ethynyl-2'-deoxyuridine click chemistry, studying single fully labeled chromosomes within an otherwise unlabeled nuclei to improve the attainable resolution. We find that, regardless of nuclear position, individual subchromosomal regions consist of three different levels of DNA compaction: (i) dispersed chromatin; (ii) nanodomains of sizes ranging tens of nanometers containing a few kilobases (kb) of DNA; and (iii) clusters of nanodomains. Interestingly, the sizes and DNA content of the nanodomains are approximately the same at the nuclear periphery, nucleolar proximity, and nuclear interior, suggesting that these nanodomains share a roughly common higher-order architecture. Overall, these results suggest that DNA compaction within the eukaryote nucleus occurs via the condensation of DNA into few-kb nanodomains of approximately similar structure, with further compaction occurring via the clustering of nanodomains.

Computational prediction of CTCF/cohesin-based intra-TAD loops that insulate chromatin contacts and gene expression in mouse liver

CTCF and cohesin are key drivers of 3D-nuclear organization, anchoring the megabase-scale Topologically Associating Domains (TADs) that segment the genome. Here, we present and validate a computational method to predict cohesin-and-CTCF binding sites that form intra-TAD DNA loops. The intra-TAD loop anchors identified are structurally indistinguishable from TAD anchors regarding binding partners, sequence conservation, and resistance to cohesin knockdown; further, the intra-TAD loops retain key functional features of TADs, including chromatin contact insulation, blockage of repressive histone mark spread, and ubiquity across tissues. We propose that intra-TAD loops form by the same loop extrusion mechanism as the larger TAD loops, and that their shorter length enables finer regulatory control in restricting enhancer-promoter interactions, which enables selective, high-level expression of gene targets of super-enhancers and genes located within repressive nuclear compartments. These findings elucidate the role of intra-TAD cohesin-and-CTCF binding in nuclear organization associated with widespread insulation of distal enhancer activity.

Eukaryotic DNA replication: Orchestrated action of multi-subunit protein complexes

Genome duplication is an essential process to preserve genetic information between generations. The eukaryotic cell cycle is composed of functionally distinct phases: G1, S, G2, and M. One of the key replicative proteins that participate at every stage of DNA replication is the Mcm2-7 complex, a replicative helicase. In the G1 phase, inactive Mcm2-7 complexes are loaded on the replication origins by replication-initiator proteins, ORC and Cdc6. Two kinases, S-CDK and DDK, convert the inactive origin-loaded Mcm2-7 complex to an active helicase, the CMG complex in the S phase. The activated CMG complex begins DNA unwinding and recruits enzymes essential for DNA synthesis to assemble a replisome at the replication fork. After completion of DNA synthesis, the inactive CMG complex on the replicated DNA is removed from chromatin to terminate DNA replication. In this review, we will discuss the structure, function, and regulation of the molecular machines involved in each step of DNA replication.

Profiling DNA Replication Timing Using Zebrafish as an In Vivo Model System

DNA replication timing is an important cellular characteristic, exhibiting significant relationships with chromatin structure, transcription, and DNA mutation rates. Changes in replication timing occur during development and in cancer, but the role replication timing plays in development and disease is not known. Zebrafish were recently established as an in vivo model system to study replication timing. Here is detailed the protocols for using the zebrafish to determine DNA replication timing. After sorting cells from embryos and adult zebrafish, high-resolution genome-wide DNA replication timing patterns can be constructed by determining changes in DNA copy number through analysis of next generation sequencing data. The zebrafish model system allows for evaluation of the replication timing changes that occur in vivo throughout development, and can also be used to assess changes in individual cell types, disease models, or mutant lines. These methods will enable studies investigating the mechanisms and determinants of replication timing establishment and maintenance during development, the role replication timing plays in mutations and tumorigenesis, and the effects of perturbing replication timing on development and disease.

Chromatin conformation regulates the coordination between DNA replication and transcription

Chromatin is the template for the basic processes of replication and transcription, making the maintenance of chromosomal integrity critical for cell viability. To elucidate how dividing cells respond to alterations in chromatin structure, here we analyse the replication programme of primary cells with altered chromatin configuration caused by the genetic ablation of the HMGB1 gene, or three histone H1 genes. We find that loss of chromatin compaction in H1-depleted cells triggers the accumulation of stalled forks and DNA damage as a consequence of transcription–replication conflicts. In contrast, reductions in nucleosome occupancy due to the lack of HMGB1 cause faster fork progression without impacting the initiation landscape or fork stability. Thus, perturbations in chromatin integrity elicit a range of responses in the dynamics of DNA replication and transcription, with different consequences on replicative stress. These findings have broad implications for our understanding of how defects in chromatin structure contribute to genomic instability.

The maintenance of chromatin integrity during replication is critical for cell viability. Here the authors study how dividing cells respond to alterations in chromatin structure and find that these elicit a range of responses in the dynamics of DNA replication and consequences on replicative stress.

Similarity in replication timing between polytene and diploid cells is associated with the organization of the Drosophila genome

Morphologically, polytene chromosomes of Drosophila melanogaster consist of compact “black” bands alternating with less compact “grey” bands and interbands. We developed a comprehensive approach that combines cytological mapping data of FlyBase-annotated genes and novel tools for predicting cytogenetic features of chromosomes on the basis of their protein composition and determined the genomic coordinates for all black bands of polytene chromosome 2R. By a PCNA immunostaining assay, we obtained the replication timetable for all the bands mapped. The results allowed us to compare replication timing between polytene chromosomes in salivary glands and chromosomes from cultured diploid cell lines and to observe a substantial similarity in the global replication patterns at the band resolution level. In both kinds of chromosomes, the intervals between black bands correspond to early replication initiation zones. Black bands are depleted of replication initiation events and are characterized by a gradient of replication timing; therefore, the time of replication completion correlates with the band length. The bands are characterized by low gene density, contain predominantly tissue-specific genes, and are represented by silent chromatin types in various tissues. The borders of black bands correspond well to the borders of topological domains as well as to the borders of the zones showing H3K27me3, SUUR, and LAMIN enrichment. In conclusion, the characteristic pattern of polytene chromosomes reflects partitioning of the Drosophila genome into two global types of domains with contrasting properties. This partitioning is conserved in different tissues and determines replication timing in Drosophila.

A TAD boundary is preserved upon deletion of the CTCF-rich Firre locus

The binding of the transcriptional regulator CTCF to the genome has been implicated in the formation of topologically associated domains (TADs). However, the general mechanisms of folding the genome into TADs are not fully understood. Here we test the effects of deleting a CTCF-rich locus on TAD boundary formation. Using genome-wide chromosome conformation capture (Hi-C), we focus on one TAD boundary on chromosome X harboring ~ 15 CTCF binding sites and located at the long non-coding RNA (lncRNA) locus Firre. Specifically, this TAD boundary is invariant across evolution, tissues, and temporal dynamics of X-chromosome inactivation. We demonstrate that neither the deletion of this locus nor the ectopic insertion of Firre cDNA or its ectopic expression are sufficient to alter TADs in a sex-specific or allele-specific manner. In contrast, Firre's deletion disrupts the chromatin super-loop formation of the inactive X-chromosome. Collectively, our findings suggest that apart from CTCF binding, additional mechanisms may play roles in establishing TAD boundary formation.

Transcriptional Dysregulation of MYC Reveals Common Enhancer-Docking Mechanism

Transcriptional dysregulation of the MYC oncogene is among the most frequent events in aggressive tumor cells, and this is generally accomplished by acquisition of a super-enhancer somewhere within the 2.8 Mb TAD where MYC resides. We find that these diverse cancer-specific super-enhancers, differing in size and location, interact with the MYC gene through a common and conserved CTCF binding site located 2 kb upstream of the MYC promoter. Genetic perturbation of this enhancer-docking site in tumor cells reduces CTCF binding, super-enhancer interaction, MYC gene expression, and cell proliferation. CTCF binding is highly sensitive to DNA methylation, and this enhancer-docking site, which is hypomethylated in diverse cancers, can be inactivated through epigenetic editing with dCas9-DNMT. Similar enhancer-docking sites occur at other genes, including genes with prominent roles in multiple cancers, suggesting a mechanism by which tumor cell oncogenes can generally hijack enhancers. These results provide insights into mechanisms that allow a single target gene to be regulated by diverse enhancer elements in different cell types.

DNA replication through a chromatin environment

Compaction of the genome into the nuclear space is achieved by wrapping DNA around octameric assemblies of histone proteins to form nucleosomes, the fundamental repeating unit of chromatin. Aside from providing a means by which to fit larger genomes into the cell, chromatinization of DNA is a crucial means by which the cell regulates access to the genome. While the complex role that chromatin plays in gene transcription has been appreciated for a long time, it is now also apparent that crucial aspects of DNA replication are linked to the biology of chromatin. This review will focus on recent advances in our understanding of how the chromatin environment influences key aspects of DNA replication.This article is part of the themed issue 'Chromatin modifiers and remodellers in DNA repair and signalling'.

Epigenetic differences between naïve and primed pluripotent stem cells

It has been 8 years since the concept of naïve and primed pluripotent stem cell states was first proposed. Both are states of pluripotency, but exhibit slightly different properties. The naïve state represents the cellular state of the preimplantation mouse blastocyst inner cell mass, while the primed state is representative of the post-implantation epiblast cells. These two cell types exhibit clearly distinct developmental potential, as evidenced by the fact that naïve cells are able to contribute to blastocyst chimeras, while primed cells cannot. However, the epigenetic differences that underlie the distinct developmental potential of these cell types remain unclear, which is rather surprising given the large amount of active investigation over the years. Elucidating such epigenetic differences should lead to a better understanding of the fundamental properties of these states of pluripotency and the means by which the naïve-to-primed transition occurs, which may provide insights into the essence of stem cell commitment.

Genome-wide analysis of replication timing by next-generation sequencing with E/L Repli-seq

This protocol is an extension to: Nat. Protoc. 6, 870-895 (2014); doi:10.1038/nprot.2011.328; published online 02 June 2011Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early- and late-replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and subnuclear position. Moreover, RT is regulated during development and is altered in diseases. Here, we describe E/L Repli-seq, an extension of our Repli-chip protocol. E/L Repli-seq is a rapid, robust and relatively inexpensive protocol for analyzing RT by next-generation sequencing (NGS), allowing genome-wide assessment of how cellular processes are linked to RT. Briefly, cells are pulse-labeled with BrdU, and early and late S-phase fractions are sorted by flow cytometry. Labeled nascent DNA is immunoprecipitated from both fractions and sequenced. Data processing leads to a single bedGraph file containing the ratio of nascent DNA from early versus late S-phase fractions. The results are comparable to those of Repli-chip, with the additional benefits of genome-wide sequence information and an increased dynamic range. We also provide computational pipelines for downstream analyses, for parsing phased genomes using single-nucleotide polymorphisms (SNPs) to analyze RT allelic asynchrony, and for direct comparison to Repli-chip data. This protocol can be performed in up to 3 d before sequencing, and requires basic cellular and molecular biology skills, as well as a basic understanding of Unix and R.

Correlative live and super-resolution imaging reveals the dynamic structure of replication domains

Combining fluorescence labeling with live-cell confocal and correlative super-resolution microscopy, Xiang et al. characterize biophysical parameters defining the internal organization, spacing, and mechanical coupling of replication domains.

Chromosome organization in higher eukaryotes controls gene expression, DNA replication, and DNA repair. Genome mapping has revealed the functional units of chromatin at the submegabase scale as self-interacting regions called topologically associating domains (TADs) and showed they correspond to replication domains (RDs). A quantitative structural and dynamic description of RD behavior in the nucleus is, however, missing because visualization of dynamic subdiffraction-sized RDs remains challenging. Using fluorescence labeling of RDs combined with correlative live and super-resolution microscopy in situ, we determined biophysical parameters to characterize the internal organization, spacing, and mechanical coupling of RDs. We found that RDs are typically 150 nm in size and contain four co-replicating regions spaced 60 nm apart. Spatially neighboring RDs are spaced 300 nm apart and connected by highly flexible linker regions that couple their motion only <550 nm. Our pipeline allows a robust quantitative characterization of chromosome structure in situ and provides important biophysical parameters to understand general principles of chromatin organization.

Developing in 3D: the role of CTCF in cell differentiation

CTCF is a highly conserved zinc-finger DNA-binding protein that mediates interactions between distant sequences in the genome. As a consequence, CTCF regulates enhancer-promoter interactions and contributes to the three-dimensional organization of the genome. Recent studies indicate that CTCF is developmentally regulated, suggesting that it plays a role in cell type-specific genome organization. Here, we review these studies and discuss how CTCF functions during the development of various cell and tissue types, ranging from embryonic stem cells and gametes, to neural, muscle and cardiac cells. We propose that the lineage-specific control of CTCF levels, and its partnership with lineage-specific transcription factors, allows for the control of cell type-specific gene expression via chromatin looping.

Dissecting super-enhancer hierarchy based on chromatin interactions

Recent studies have highlighted super-enhancers (SEs) as important regulatory elements for gene expression, but their intrinsic properties remain incompletely characterized. Through an integrative analysis of Hi-C and ChIP-seq data, here we find that a significant fraction of SEs are hierarchically organized, containing both hub and non-hub enhancers. Hub enhancers share similar histone marks with non-hub enhancers, but are distinctly associated with cohesin and CTCF binding sites and disease-associated genetic variants. Genetic ablation of hub enhancers results in profound defects in gene activation and local chromatin landscape. As such, hub enhancers are the major constituents responsible for SE functional and structural organization.

Genome-Wide Analysis of the Arabidopsis Replication Timing Program

Eukaryotes use a temporally regulated process, known as the replication timing program, to ensure that their genomes are fully and accurately duplicated during S phase. Replication timing programs are predictive of genomic features and activity and are considered to be functional readouts of chromatin organization. Although replication timing programs have been described for yeast and animal systems, much less is known about the temporal regulation of plant DNA replication or its relationship to genome sequence and chromatin structure. We used the thymidine analog, 5-ethynyl-2'-deoxyuridine, in combination with flow sorting and Repli-Seq to describe, at high-resolution, the genome-wide replication timing program for Arabidopsis (Arabidopsis thaliana) Col-0 suspension cells. We identified genomic regions that replicate predominantly during early, mid, and late S phase, and correlated these regions with genomic features and with data for chromatin state, accessibility, and long-distance interaction. Arabidopsis chromosome arms tend to replicate early while pericentromeric regions replicate late. Early and mid-replicating regions are gene-rich and predominantly euchromatic, while late regions are rich in transposable elements and primarily heterochromatic. However, the distribution of chromatin states across the different times is complex, with each replication time corresponding to a mixture of states. Early and mid-replicating sequences interact with each other and not with late sequences, but early regions are more accessible than mid regions. The replication timing program in Arabidopsis reflects a bipartite genomic organization with early/mid-replicating regions and late regions forming separate, noninteracting compartments. The temporal order of DNA replication within the early/mid compartment may be modulated largely by chromatin accessibility.

Bacterial artificial chromosomes establish replication timing and sub-nuclear compartment de novo as extra-chromosomal vectors

The role of DNA sequence in determining replication timing (RT) and chromatin higher order organization remains elusive. To address this question, we have developed an extra-chromosomal replication system (E-BACs) consisting of ∼200 kb human bacterial artificial chromosomes (BACs) modified with Epstein-Barr virus (EBV) stable segregation elements. E-BACs were stably maintained as autonomous mini-chromosomes in EBNA1-expressing HeLa or human induced pluripotent stem cells (hiPSCs) and established distinct RT patterns. An E-BAC harboring an early replicating chromosomal region replicated early during S phase, while E-BACs derived from RT transition regions (TTRs) and late replicating regions replicated in mid to late S phase. Analysis of E-BAC interactions with cellular chromatin (4C-seq) revealed that the early replicating E-BAC interacted broadly throughout the genome and preferentially with the early replicating compartment of the nucleus. In contrast, mid- to late-replicating E-BACs interacted with more specific late replicating chromosomal segments, some of which were shared between different E-BACs. Together, we describe a versatile system in which to study the structure and function of chromosomal segments that are stably maintained separately from the influence of cellular chromosome context.

CDK activity provides temporal and quantitative cues for organizing genome duplication

In eukaryotes, the spatial and temporal organization of genome duplication gives rise to distinctive profiles of replication origin usage along the chromosomes. While it has become increasingly clear that these programs are important for cellular physiology, the mechanisms by which they are determined and modulated remain elusive. Replication initiation requires the function of cyclin-dependent kinases (CDKs), which associate with various cyclin partners to drive cell proliferation. Surprisingly, although we possess detailed knowledge of the CDK regulators and targets that are crucial for origin activation, little is known about whether CDKs play a critical role in establishing the genome-wide pattern of origin selection. We have addressed this question in the fission yeast, taking advantage of a simplified cell cycle network in which cell proliferation is driven by a single cyclin-CDK module. This system allows us to precisely control CDK activity in vivo using chemical genetics. First, in contrast to previous reports, our results clearly show that distinct cyclin-CDK pairs are not essential for regulating specific subsets of origins and for establishing a normal replication program. Importantly, we then demonstrate that the timing at which CDK activity reaches the S phase threshold is critical for the organization of replication in distinct efficiency domains, while the level of CDK activity at the onset of S phase is a dose-dependent modulator of overall origin efficiencies. Our study therefore implicates these different aspects of CDK regulation as versatile mechanisms for shaping the architecture of DNA replication across the genome.

The duplication of the genetic material is a highly conserved and tightly regulated process that is essential for cell proliferation. DNA synthesis initiates at sites called origins that are distributed throughout the genome. Replication origins are not all used equivalently, and their patterns of activation along the chromosomes give rise to specific profiles, or programs, of DNA replication. These programs change during development and in response to external stimuli, and we have previously shown that they have important consequences for cellular function. However, we still do not understand the mechanisms by which cells establish different replication patterns. Here we investigate the role of the cyclin-dependent kinase (CDK) family of proteins, whose activities are critical for cell cycle progression, in regulating the organization of genome duplication. Taking advantage of a system that allows us to precisely modulate CDK activity levels in living cells, we demonstrate that both the temporal and quantitative controls of CDK function are crucial for determining distinct programs of DNA replication. Our work therefore uncovers a fundamental link between CDK activity, a central input in a variety of cellular and developmental processes, and the architecture of genome duplication.

PREP1 tumor suppressor protects the late-replicating DNA by controlling its replication timing and symmetry

The synthesis of middle-to-late-replicating DNA can be affected independently of the rest of the genome by down-regulating the tumor suppressor PREP1 (PKNOX1). Indeed, DNA combing shows that PREP1 down-regulation affects DNA replication rate, increases the number of simultaneously firing origins and the asymmetry of DNA replication, leading to DNA damage. Genome-wide analysis of replication timing by Repli-seq shows that, upon PREP1 down-regulation, 25% of the genome is replicated earlier in the S-phase. The targeted DNA sequences correspond to Lamin-Associated Domains (LADs), and include late-replicating (LRRs) and temporal transition regions (TTRs). Notably, the distribution of PREP1 DNA binding sites and of its target genes indicates that DNA replication defects are independent of the overall PREP1 transcriptional activity. Finally, PREP1 down-regulation causes a substantial decrease in Lamin B1 levels. This suggests that DNA is released from the nuclear lamina earlier than in the control cells and is available for replication, thus explaining timing defects and DNA damage.This is the first evidence that the replication timing of a specific fraction of the human genome is affected by PREP1 tumor suppressor. This previously unknown function might significantly contribute to the genomic instability observed in human tumors.

Sequence-based multiscale modeling for high-throughput chromosome conformation capture (Hi-C) data analysis

In this paper, we introduce sequence-based multiscale modeling for biomolecular data analysis. We employ spectral clustering method in our modeling and reveal the difference between sequence-based global scale clustering and local scale clustering. Essentially, two types of distances, i.e., Euclidean (or spatial) distance and genomic (or sequential) distance, can be used in data clustering. Clusters from sequence-based global scale models optimize spatial distances, meaning spatially adjacent loci are more likely to be assigned into the same cluster. Sequence-based local scale models, on the other hand, result in clusters that optimize genomic distances. That is to say, in these models, sequentially adjoining loci tend to be cluster together. We propose two sequence-based multiscale models (SeqMMs) for the study of chromosome hierarchical structures, including genomic compartments and topological associated domains (TADs). We find that genomic compartments are determined only by global scale information in the Hi-C data. The removal of all the local interactions within a band region as large as 10 Mb in genomic distance has almost no significant influence on the final compartment results. Further, in TAD analysis, we find that when the sequential scale is small, a tiny variation of diagonal band region in a contact map will result in a great change in the predicted TAD boundaries. When the scale value is larger than a threshold value, the TAD boundaries become very consistent. This threshold value is highly related to TAD sizes. By the comparison of our results with those previously obtained using a spectral clustering model, we find that our method is more robust and reliable. Finally, we demonstrate that almost all TAD boundaries from both clustering methods are local minimum of a TAD summation function.

Similar Evolutionary Trajectories for Retrotransposon Accumulation in Mammals

The factors guiding retrotransposon insertion site preference are not well understood. Different types of retrotransposons share common replication machinery and yet occupy distinct genomic domains. Autonomous long interspersed elements accumulate in gene-poor domains and their nonautonomous short interspersed elements accumulate in gene-rich domains. To determine genomic factors that contribute to this discrepancy we analyzed the distribution of retrotransposons within the framework of chromosomal domains and regulatory elements. Using comparative genomics, we identified large-scale conserved patterns of retrotransposon accumulation across several mammalian genomes. Importantly, retrotransposons that were active after our sample-species diverged accumulated in orthologous regions. This suggested a similar evolutionary interaction between retrotransposon activity and conserved genome architecture across our species. In addition, we found that retrotransposons accumulated at regulatory element boundaries in open chromatin, where accumulation of particular retrotransposon types depended on insertion size and local regulatory element density. From our results, we propose a model where density and distribution of genes and regulatory elements canalize retrotransposon accumulation. Through conservation of synteny, gene regulation and nuclear organization, mammalian genomes with dissimilar retrotransposons follow similar evolutionary trajectories.

Transcriptional decomposition reveals active chromatin architectures and cell specific regulatory interactions

Transcriptional regulation is tightly coupled with chromosomal positioning and three-dimensional chromatin architecture. However, it is unclear what proportion of transcriptional activity is reflecting such organisation, how much can be informed by RNA expression alone and how this impacts disease. Here, we develop a computational transcriptional decomposition approach separating the proportion of expression associated with genome organisation from independent effects not directly related to genomic positioning. We show that positionally attributable expression accounts for a considerable proportion of total levels and is highly informative of topological associating domain activities and organisation, revealing boundaries and chromatin compartments. Furthermore, expression data alone accurately predict individual enhancer-promoter interactions, drawing features from expression strength, stabilities, insulation and distance. We characterise predictions in 76 human cell types, observing extensive sharing of domains, yet highly cell-type-specific enhancer-promoter interactions and strong enrichments in relevant trait-associated variants. Overall, our work demonstrates a close relationship between transcription and chromatin architecture.

Transcription factors orchestrate dynamic interplay between genome topology and gene regulation during cell reprogramming

Chromosomal architecture is known to influence gene expression, yet its role in controlling cell fate remains poorly understood. Reprogramming of somatic cells into pluripotent stem cells (PSCs) by the transcription factors (TFs) OCT4, SOX2, KLF4 and MYC offers an opportunity to address this question but is severely limited by the low proportion of responding cells. We have recently developed a highly efficient reprogramming protocol that synchronously converts somatic into pluripotent stem cells. Here, we used this system to integrate time-resolved changes in genome topology with gene expression, TF binding and chromatin-state dynamics. The results showed that TFs drive topological genome reorganization at multiple architectural levels, often before changes in gene expression. Removal of locus-specific topological barriers can explain why pluripotency genes are activated sequentially, instead of simultaneously, during reprogramming. Together, our results implicate genome topology as an instructive force for implementing transcriptional programs and cell fate in mammals.

Single-cell replication profiling to measure stochastic variation in mammalian replication timing

Mammalian DNA replication is regulated via multi-replicon segments that replicate in a defined temporal order during S-phase. Further, early/late replication of RDs corresponds to active/inactive chromatin interaction compartments. Although replication origins are selected stochastically, variation in replication timing is poorly understood. Here we devise a strategy to measure variation in replication timing using DNA copy number in single mouse embryonic stem cells. We find that borders between replicated and unreplicated DNA are highly conserved between cells, demarcating active and inactive compartments of the nucleus. Fifty percent of replication events deviated from their average replication time by ± 15% of S phase. This degree of variation is similar between cells, between homologs within cells and between all domains genomewide, regardless of their replication timing. These results demonstrate that stochastic variation in replication timing is independent of elements that dictate timing or extrinsic environmental variation.

The 4D nucleome project

The 4D Nucleome Network aims to develop and apply approaches to map the structure and dynamics of the human and mouse genomes in space and time with the goal of gaining deeper mechanistic insights into how the nucleus is organized and functions. The project will develop and benchmark experimental and computational approaches for measuring genome conformation and nuclear organization, and investigate how these contribute to gene regulation and other genome functions. Validated experimental technologies will be combined with biophysical approaches to generate quantitative models of spatial genome organization in different biological states, both in cell populations and in single cells.

Signaling by Steroid Hormones in the 3D Nuclear Space

Initial studies showed that ligand-activated hormone receptors act by binding to the proximal promoters of individual target genes. Genome-wide studies have now revealed that regulation of transcription by steroid hormones mainly depends on binding of the receptors to distal regulatory elements. Those distal elements, either enhancers or silencers, act on the regulation of target genes by chromatin looping to the gene promoters. In the nucleus, this level of chromatin folding is integrated within dynamic higher orders of genome structures, which are organized in a non-random fashion. Terminally differentiated cells exhibit a tissue-specific three-dimensional (3D) organization of the genome that favors or restrains the activity of transcription factors and modulates the function of steroid hormone receptors, which are transiently activated upon hormone exposure. Conversely, integration of the hormones signal may require modifications of the 3D organization to allow appropriate transcriptional outcomes. In this review, we summarize the main levels of organization of the genome, review how they can modulate the response to steroids in a cell specific manner and discuss the role of receptors in shaping and rewiring the structure in response to hormone. Taking into account the dynamics of 3D genome organization will contribute to a better understanding of the pleiotropic effects of steroid hormones in normal and cancer cells.

Sub-kb Hi-C in D. melanogaster reveals conserved characteristics of TADs between insect and mammalian cells

Topologically associating domains (TADs) are fundamental elements of the eukaryotic genomic structure. However, recent studies suggest that the insulating complexes, CTCF/cohesin, present at TAD borders in mammals are absent from those in Drosophila melanogaster, raising the possibility that border elements are not conserved among metazoans. Using in situ Hi-C with sub-kb resolution, here we show that the D. melanogaster genome is almost completely partitioned into >4000 TADs, nearly sevenfold more than previously identified. The overwhelming majority of these TADs are demarcated by the insulator complexes, BEAF-32/CP190, or BEAF-32/Chromator, indicating that these proteins may play an analogous role in flies as that of CTCF/cohesin in mammals. Moreover, extended regions previously thought to be unstructured are shown to consist of small contiguous TADs, a property also observed in mammals upon re-examination. Altogether, our work demonstrates that fundamental features associated with the higher-order folding of the genome are conserved from insects to mammals.