Search

Cells respond to many stressors by senescing, acquiring stable growth arrest, morphologic and metabolic changes, and a proinflammatory senescence-associated secretory phenotype. The heterogeneity of senescent cells (SnCs) and senescence-associated secretory phenotype are vast, yet ill characterized. SnCs have diverse roles in health and disease and are therapeutically targetable, making characterization of SnCs and their detection a priority. The Cellular Senescence Network (SenNet), a National Institutes of Health Common Fund initiative, was established to address this need. The goal of SenNet is to map SnCs across the human lifespan to advance diagnostic and therapeutic approaches to improve human health. State-of-the-art methods will be applied to identify, define and map SnCs in 18 human tissues. A common coordinate framework will integrate data to create four-dimensional SnC atlases. Other key SenNet deliverables include innovative tools and technologies to detect SnCs, new SnC biomarkers and extensive public multi-omics datasets. This Perspective lays out the impetus, goals, approaches and products of SenNet.

Dequeker and colleagues performed elegant in vivo, in silico, and in vitro experiments to demonstrate that the MCM complex, an essential DNA replication factor, is an obstacle for the DNA loop formation by cohesin.

Cells respond to a myriad of stressors by senescing, acquiring stable growth arrest, morphologic and metabolic changes, and a senescence-associated-secretory-phenotype (SASP). The heterogeneity of senescent cells (SnCs) and their SASP is vast, yet poorly characterized. SnCs have diverse roles in health and disease and are therapeutically targetable, making characterization of SnCs and harmonization of their nomenclature a priority. The Cellular Senescence Network (SenNet), a NIH Common Fund initiative, will leverage emerging single cell and spatial-omics to identify and map SnCs in numerous organs across the lifespan of humans and mice. A common coordinate framework will integrate the data, using validated, standardized methods, creating public 4-dimensional SnC atlases. Key SenNet deliverables include development of innovative tools/technologies to detect SnCs, biomarker discovery, common annotations to describe SnCs and extensive public data sets. The goal is to comprehensively understand and map SnCs for diagnostic and therapeutic purposes to improve human health.

Currently, there is a significant shortage of transplantable organs for patients in need. Interspecies chimerism and blastocyst complementation are alternatives for generating transplantable human organs in host animals such as pigs to meet this shortage. While successful interspecies chimerism and organ generation have been observed between evolutionarily close species such as rat and mouse, barriers still exist for more distant species pairs such as human–mouse, marmoset–mouse, human–pig, and others. One of the proposed barriers to chimerism is the difference in developmental stages between the donor cells and the host embryo at the time the cells are introduced into the host embryo. Hence, there is a logical effort to stage-match the donor cells with the host embryos for enhancing interspecies chimerism. In this study, we used an in silico approach to simultaneously stage-match the early developing embryos of four species, including human, marmoset, mouse, and pig based on transcriptome similarities. We used an unsupervised clustering algorithm to simultaneously stage-match all four species as well as Spearman’s correlation analyses to stage-match pairs of donor–host species. From our stage-matching analyses, we found that the four stages that best matched with each other are the human blastocyst (E6/E7), the gastrulating mouse embryo (E6–E6.75), the marmoset late inner cell mass, and the pig late blastocyst. We further demonstrated that human pluripotent stem cells best matched with the mouse post-implantation stages. We also performed ontology analysis of the genes upregulated and commonly expressed between donor–host species pairs at their best matched stages. The stage-matching results predicted by this study will inform in vivo and in vitro interspecies chimerism and blastocyst complementation studies and can be used to match donor cells with host embryos between multiple species pairs to enhance chimerism for organogenesis.

Summary The temporal order of DNA replication is regulated during development and is highly correlated with gene expression, histone modifications and 3D genome architecture. We tracked changes in replication timing, gene expression, and chromatin conformation capture (Hi-C) A/B compartments over the first two cell cycles during differentiation of human embryonic stem cells to definitive endoderm. Remarkably, transcriptional programs were irreversibly reprogrammed within the first cell cycle and were largely but not universally coordinated with replication timing changes. Moreover, changes in A/B compartment and several histone modifications that normally correlate strongly with replication timing showed weak correlation during the early cell cycles of differentiation but showed increased alignment in later differentiation stages and in terminally differentiated cell lines. Thus, epigenetic cell fate transitions during early differentiation can occur despite dynamic and discordant changes in otherwise highly correlated genomic properties., Graphical Abstract, Highlights • Stable transcriptional reprogramming toward endoderm within one cell cycle • Replication timing can change in the first S phase of a cell fate change • Early discordance between Replication timing and Hi-C compartments • Alignment of replication timing to compartments increases in later cell cycles, The temporal order of DNA replication is regulated during development, and is highly correlated with gene expression, chromatin structure, and 3D-genome architecture. Gilbert and colleagues find that stable transcriptional changes of human embryonic stem cells to endoderm can occur within a single cell cycle, accompanied by discordance in replication timing and chromatin compartment that align in later differentiation stages.

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes., The authors summarize the data produced by phase III of the Encyclopedia of DNA Elements (ENCODE) project, a resource for better understanding of the human and mouse genomes.

Common fragile sites (CFSs) are regions susceptible to replication stress and are hotspots for chromosomal instability in cancer. Several features were suggested to underlie CFS instability, however, these features are prevalent across the genome. Therefore, the molecular mechanisms underlying CFS instability remain unclear. Here, we explore the transcriptional profile and DNA replication timing (RT) under mild replication stress in the context of the 3D genome organization. The results reveal a fragility signature, comprised of a TAD boundary overlapping a highly transcribed large gene with APH-induced RT-delay. This signature enables precise mapping of core fragility regions in known CFSs and identification of novel fragile sites. CFS stability may be compromised by incomplete DNA replication and repair in TAD boundaries core fragility regions leading to genomic instability. The identified fragility signature will allow for a more comprehensive mapping of CFSs and pave the way for investigating mechanisms promoting genomic instability in cancer., Common fragile sites are regions susceptible to replication stress and are prone to chromosomal instability. Here, the authors, by analyzing the contribution of 3D chromatin organization, identify and characterize a fragility signature and precisely map these fragility regions.

To understand the role of the extensive senescence-associated 3D genome reorganization, we generated genome-wide chromatin interaction maps, epigenome, replication-timing, whole-genome bisulfite sequencing, and gene expression profiles from cells entering replicative senescence (RS) or upon oncogene-induced senescence (OIS). We identify senescence-associated heterochromatin domains (SAHDs). Differential intra- versus inter-SAHD interactions lead to the formation of senescence-associated heterochromatin foci (SAHFs) in OIS but not in RS. This OIS-specific configuration brings active genes located in genomic regions adjacent to SAHDs in close spatial proximity and favors their expression. We also identify DNMT1 as a factor that induces SAHFs by promoting HMGA2 expression. Upon DNMT1 depletion, OIS cells transition to a 3D genome conformation akin to that of cells in replicative senescence. These data show how multi-omics and imaging can identify critical features of RS and OIS and discover determinants of acute senescence and SAHF formation. Work at the M.A.M.-R. lab was supported by the European Research Council under the 7th Framework Program FP7/2007-2013 (ERC grant agreement 609989), the European Union’s Horizon 2020 research and innovation programme (grant agreement 676556), the Ministry of Economy and Competitiveness (BFU2017-85926-P), and the Agència de Gestió d’Ajuts Universitaris i de Recerca, AGAUR (SGR468). Work at CRG, BIST, and UPF was in part funded by the Spanish Ministry of Economy and Competitiveness, ‘‘Centro de Excelencia Severo Ochoa 2013-2017’’ (SEV-2012-0208), and ‘‘Centro de Excelencia María de Maeztu 2016-2019.’’ This article/publication is based upon work from COST Action CA18127, supported by COST (European Cooperation in Science and Technology)

Structural variants can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of structural variants in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify structural variants in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of non-coding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel 3D chromatin structural domains. Our results indicate that non-coding structural variations may be underappreciated mutational drivers in cancer genomes.

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.

In metazoans, nuclear DNA is organized during the interphase in negatively supercoiled loops anchored to a compartment or substructure known as the nuclear matrix. The interactions between DNA and the nuclear matrix (NM) are of higher affinity than those between DNA and chromatin proteins since the last ones do not resist the procedures for extracting the NM. The structural interactions DNA-NM constitute a set of topological relationships that define a nuclear higher order structure (NHOS) although there are further higher order levels of organization within the nucleus. So far, the evidence derived from studies with primary hepatocytes and naive B lymphocytes indicates that the NHOS is cell-type specific at the local and at the large-scale level, and so it has been suggested that such NHOS is primary determined by structural and thermodynamic constraints. We carried out a comparative characterization of the NHOS of postmitotic cortical neurons with that of hepatocytes and naive B lymphocytes. Our results indicate that the NHOS of neurons is completely different at the large scale and at the local level from that one observed in hepatocytes or in naive B lymphocytes, confirming on the one hand that the set of structural DNA-NM interactions is cell-type specific and supporting, on the other hand the notion that structural constraints that impinge on chromosomal DNA and the NM are more important for determining this NHOS than functional constraints related to replication and/or transcription. J. Cell. Biochem. 118: 2151-2160, 2017. © 2016 Wiley Periodicals, Inc.

DNA replication occurs in a defined temporal order known as the replication timing (RT) program and is regulated during development, coordinated with 3D genome organization and transcriptional activity. However, transcription and RT are not sufficiently coordinated to predict each other, suggesting an indirect relationship. Here, we exploit genome-wide RT profiles from 15 human cell types and intermediate differentiation stages derived from human embryonic stem cells to construct different types of RT regulatory networks. First, we constructed networks based on the coordinated RT changes during cell fate commitment to create highly complex RT networks composed of thousands of interactions that form specific functional subnetwork communities. We also constructed directional regulatory networks based on the order of RT changes within cell lineages, and identified master regulators of differentiation pathways. Finally, we explored relationships between RT networks and transcriptional regulatory networks (TRNs) by combining them into more complex circuitries of composite and bipartite networks. Results identified novel trans interactions linking transcription factors that are core to the regulatory circuitry of each cell type to RT changes occurring in those cell types. These core transcription factors were found to bind cooperatively to sites in the affected replication domains, providing provocative evidence that they constitute biologically significant directional interactions. Our findings suggest a regulatory link between the establishment of cell-type–specific TRNs and RT control during lineage specification.

Common fragile sties (CFSs) are regions susceptible to replication stress and are hotspots for chromosomal instability in cancer. Several features characterizing CFSs have been associated with their instability, however, these features are prevalent across the genome and do not account for all known CFSs. Therefore, the molecular mechanism underlying CFS instability remains unclear. Here, we explored the transcriptional profile and temporal order of DNA replication (replication timing, RT) of cells under replication stress conditions. We show that the RT of only a small portion of the genome is affected by replication stress, and that CFSs are enriched for delayed RT. We identified a signature for chromosomal fragility, comprised of replication stress-induced delay in RT of early/mid S-phase replicating regions within actively transcribed large genes. This fragility signature enabled precise mapping of the core fragility region. Furthermore, the signature enabled the identification of novel fragile sites that were not detected cytogenetically, highlighting the improved sensitivity of our approach for identifying fragile sites. Altogether, this study reveals a link between altered DNA replication and transcription of large genes underlying the mechanism of CFS expression. Thus, investigating the RT and transcriptional changes in cancer may contribute to the understanding of mechanisms promoting genomic instability in cancer.

ENCODE comprises thousands of functional genomics datasets, and the encyclopedia covers hundreds of cell types, providing a universal annotation for genome interpretation. However, for particular applications, it may be advantageous to use a customized annotation. Here, we develop such a custom annotation by leveraging advanced assays, such as eCLIP, Hi-C, and whole-genome STARR-seq on a number of data-rich ENCODE cell types. A key aspect of this annotation is comprehensive and experimentally derived networks of both transcription factors and RNA-binding proteins (TFs and RBPs). Cancer, a disease of system-wide dysregulation, is an ideal application for such a network-based annotation. Specifically, for cancer-associated cell types, we put regulators into hierarchies and measure their network change (rewiring) during oncogenesis. We also extensively survey TF-RBP crosstalk, highlighting how SUB1, a previously uncharacterized RBP, drives aberrant tumor expression and amplifies the effect of MYC, a well-known oncogenic TF. Furthermore, we show how our annotation allows us to place oncogenic transformations in the context of a broad cell space; here, many normal-to-tumor transitions move towards a stem-like state, while oncogene knockdowns show an opposing trend. Finally, we organize the resource into a coherent workflow to prioritize key elements and variants, in addition to regulators. We showcase the application of this prioritization to somatic burdening, cancer differential expression and GWAS. Targeted validations of the prioritized regulators, elements and variants using siRNA knockdowns, CRISPR-based editing, and luciferase assays demonstrate the value of the ENCODE resource., ENCODE is a resource comprising thousands of functional genomic datasets. Here, the authors present custom annotation within ENCODE for cancer, highlighting a workflow that can help prioritise key elements in oncogenesis.

Human B-lineage precursor acute lymphoid leukemias (BCP-ALLs) comprise a group of genetically and clinically distinct disease entities with features of differentiation arrest at known stages of normal B-lineage differentiation. We previously showed BCP-ALL cells display unique and clonally heritable DNA-replication timing (RT) programs; i.e., programs describing the variable order of replication of megabase-scale chromosomal units of DNA in different cell types. To determine the extent to which BCP-ALL RT programs mirror or deviate from specific stages of normal human B-cell differentiation, we transplanted immunodeficient mice with quiescent normal human CD34+ cord blood cells and obtained RT signatures of the regenerating B-lineage populations. We then compared these with RT signatures for leukemic cells from a large cohort of BCP-ALL patients. The results identify BCP-ALL subtype-specific features that resemble specific stages of B-cell differentiation and features that appear associated with relapse. These results suggest the genesis of BCP-ALL involves alterations in RT that reflect clinically relevant leukemia-specific genetic and/or epigenetic changes.SUMMARYGenome-wide DNA replication timing profiles of >100 pediatric leukemic samples and normally differentiating human B-lineage cells isolated from xenografted immunodeficient mice were generated. Comparison of these identified potentially clinically relevant features that both match and deviate from the normal profiles.

Complete duplication of large metazoan chromosomes requires thousands of potential initiation sites, only a small fraction of which are selected in each cell cycle. Assembly of the replication machinery is highly conserved and tightly regulated during the cell cycle, but the sites of initiation are highly flexible, and their temporal order of firing is regulated at the level of large-scale multi-replicon domains. Importantly, the number of replication forks must be quickly adjusted in response to replication stress to prevent genome instability. Here we argue that large genomes are divided into domains for exactly this reason. Once established, domain structure abrogates the need for precise initiation sites and creates a scaffold for the evolution of other chromosome functions.

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 corresponding to topologically-associating domains (TADs) that are spatially compartmentalized in the nucleus. Correlations of RT to chromatin organization and gene regulation have been known for decades but causal and mechanistic links remain unknown. The complete elucidation of these intriguing liaisons is critical to understand the connection between the three-dimensional organization of the nucleus and cellular function. Here, we discuss emerging evidence providing new insights into regulation of chromosome architecture, transcription and its connections with RT.

The replication timing (RT) program has been linked to many key biological processes including cell fate commitment, 3D chromatin organization and transcription regulation. Significant technology progress now allows to characterize the RT program in the entire human genome in a high-throughput and high-resolution fashion. These experiments suggest that RT changes dynamically during development in coordination with gene activity. Since RT is such a fundamental biological process, we believe that an effective quantitative profile of the local RT program from a diverse set of cell types in various developmental stages and lineages can provide crucial biological insights for a genomic locus. In the present study, we explored recurrent and spatially coherent combinatorial profiles from 42 RT programs collected from multiple lineages at diverse differentiation states. We found that a Hidden Markov Model with 15 hidden states provide a good model to describe these genome-wide RT profiling data. Each of the hidden state represents a unique combination of RT profiles across different cell types which we refer to as “RT states”. To understand the biological properties of these RT states, we inspected their relationship with chromatin states, gene expression, functional annotation and 3D chromosomal organization. We found that the newly defined RT states possess interesting genome-wide functional properties that add complementary information to the existing annotation of the human genome.AUTHOR SUMMARYThe replication timing (RT) program is an important cellular mechanism and has been linked to many key biological processes including cell fate commitment, 3D chromatin organization and transcription regulation. Significant technology progress now allows us to characterize the RT program in the entire human genome. Results from these experiments suggest that RT changes dynamically across different developmental stages. Since RT is such a fundamental biological process, we believe that the local RT program from a diverse set of cell types in various developmental stages can provide crucial biological insights for a genomic locus. In the present study, we explored combinatorial profiles from 42 RT programs collected from multiple lineages at diverse differentiation states. We developed a statistical model consist of 15 “RT states” to describe these genome-wide RT profiling data. To understand the biological properties of these RT states, we inspected the relationship between RT states and other types of functional annotations of the genome. We found that the newly defined RT states possess interesting genome-wide functional properties that add complementary information to the existing annotation of the human genome.

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.

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.

SUMMARYThe temporal order of DNA replication (replication timing, RT) is highly coupled with genome architecture, butcis-elements regulating spatio-temporal control of replication have remained elusive. We performed an extensive series of CRISPR mediated deletions and inversions and high-resolution capture Hi-C of a pluripotency associated domain (DppA2/4) in mouse embryonic stem cells. Whereas CTCF mediated loops and chromatin domain boundaries were dispensable, deletion of three intra-domain prominent CTCF-independent 3D contact sites caused a domain-wide delay in RT, shift in sub-nuclear chromatin compartment and loss of transcriptional activity, These “early replication control elements” (ERCEs) display prominent chromatin features resembling enhancers/promoters and individual and pair-wise deletions of the ERCEs confirmed their partial redundancy and interdependency in controlling domain-wide RT and transcription. Our results demonstrate that discretecis-regulatory elements mediate domain-wide RT, chromatin compartmentalization, and transcription, representing a major advance in dissecting the relationship between genome structure and function.Highlightscis-elements (ERCEs) regulate large scale chromosome structure and functionMultiple ERCEs cooperatively control domain-wide replicationERCEs harbor prominent active chromatin features and form CTCF-independent loopsERCEs enable genetic dissection of large-scale chromosome structure-function.

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and thenecessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

DNA replication occurs in a defined temporal order during S phase, known as the replication timing programme, which is regulated not only during the cell cycle but also during the process of development and differentiation. The units of replication timing regulation, known as replication domains (RDs), frequently comprise several nearly synchronously firing replication origins. Replication domains correspond to topologically associating domains (TADs) mapped by chromatin conformation capture methods and are likely to be the molecular equivalents of replication foci observed using cytogenetic methods. Both TAD and replication foci are considered to be stable structural units of chromosomes, conserved through the cell cycle and development, and accordingly, the boundaries of RDs also appear to be stable in different cell types. During both normal development and progression of disease, distinct cell states are characterized by unique replication timing signatures, with approximately half of genomic RDs switching replication timing between these cell states. Advances in functional genomics provide hope that we can soon gain an understanding of the cause and consequence of the replication timing programme and its myriad correlations with chromatin context and gene regulation.

The temporal order of DNA replication (replication timing, RT) is highly coupled with genome architecture, but cis-elements regulating spatio-temporal control of replicatio have remained elusive. We performed an extensive series of CRISPR mediated deletions and inversions and high-resolution capture Hi-C of a pluripotency associated domain (DppA2/4) in mouse embryonic stem cells. Whereas CTCF mediated loops and chromatin domain boundaries were dispensable, deletion of three intra-domain prominent CTCF-independent 3D contact sites caused a domain-wide delay in RT, shift in sub-nuclear chromatin compartment and loss of transcriptional activity, These “early replication control elements” (ERCEs) display prominent chromatin features resembling enhancers/promoters and individual and pair-wise deletions of the ERCEs confirmed their partial redundancy and interdependency in controlling domain-wide RT and transcription. Our results demonstrate that discrete cis-regulatory elements mediate domain-wide RT, chromatin compartmentalization, and transcription, representing a major advance in dissecting the relationship between genome structure and function.

DNA replication occurs in a defined temporal order known as the replication-timing (RT) program and is regulated during development, coordinated with 3D genome organization and transcriptional activity. Here, we exploit genome-wide RT profiles from 15 human cell types and intermediate differentiation stages derived from human embryonic stem cells to construct different types of RT regulatory networks. First, we constructed networks based on the coordinated RT changes during cell fate commitment to create RT networks composed of specific functional sub-network communities. We also constructed directional regulatory networks based on the order of RT changes within cell lineages and identified master regulators of differentiation pathways. Finally, we explored relationships between RT networks and transcriptional regulatory networks (TRNs), by combining them into more complex circuitries of composite and bipartite networks. Our findings show that RT networks can be exploited to dissect the cellular mechanisms that regulate lineage specification and cellular identity maintenance.

International audience; Progeroid syndromes are rare genetic disorders that phenotypically resemble natural aging. Different causal mutations have been identified, but no molecular alterations have been identified that are in common to these diseases. DNA replication timing (RT) is a robust cell type-specific epigenetic feature highly conserved in the same cell types from different individuals but altered in disease. Here, we characterized DNA RT program alterations in Hutchinson–Gilford progeria syndrome (HGPS) and Rothmund–Thomson syndrome (RTS) patients compared with natural aging and cellular senescence. Our results identified a progeroid-specific RT signature that is common to cells from three HGPS and three RTS patients and distinguishes them from healthy individuals across a wide range of ages. Among the RT abnormalities, we identified the tumor protein p63 gene (TP63) as a gene marker for progeroid syndromes. By using the redifferentiation of four patient-derived induced pluripotent stem cells as a model for the onset of progeroid syndromes, we tracked the progression of RT abnormalities during development, revealing altered RT of the TP63 gene as an early event in disease progression of both HGPS and RTS. Moreover, the RT abnormalities in progeroid patients were associated with altered isoform expression of TP63. Our findings demonstrate the value of RT studies to identify biomarkers not detected by other methods, reveal abnormal TP63 RT as an early event in progeroid disease progression, and suggest TP63 gene regulation as a potential therapeutic target.

DNA replication occurs in a defined temporal order known as the replication-timing (RT) program and is regulated during development, coordinated with 3D genome organization and transcriptional activity. However, transcription and RT are not sufficiently coordinated to predict each other, suggesting an indirect relationship. Here, we exploit genome-wide RT profiles from 15 human cell types and intermediate differentiation stages derived from human embryonic stem cells to construct different types of RT regulatory networks. First, we constructed networks based on the coordinated RT changes during cell fate commitment to create highly complex RT networks composed of thousands of interactions that form specific functional sub-network communities. We also constructed directional regulatory networks based on the order of RT changes within cell lineages and identified master regulators of differentiation pathways. Finally, we explored relationships between RT networks and transcriptional regulatory networks (TRNs), by combining them into more complex circuitries of composite and bipartite networks, revealing novel trans interactions between transcription factors and downstream RT changes that were validated with ChIP-seq data. Our findings suggest a regulatory link between the establishment of cell type specific TRNs and RT control during lineage specification.

Structural variants can contribute to oncogenesis through a variety of mechanisms, yet, despite their importance, the identification of structural variants in cancer genomes remains challenging. Here, we present an integrative framework for comprehensively identifying structural variation in cancer genomes. For the first time, we apply next-generation optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole genome sequencing to systematically detect SVs in a variety of cancer cells.Using this approach, we identify and characterize structural variants in up to 29 commonly used normal and cancer cell lines. We find that each method has unique strengths in identifying different classes of structural variants and at different scales, suggesting that integrative approaches are likely the only way to comprehensively identify structural variants in the genome. Studying the impact of the structural variants in cancer cell lines, we identify widespread structural variation events affecting the functions of non-coding sequences in the genome, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel 3D chromatin structural domains.These results underscore the importance of comprehensive structural variant identification and indicate that non-coding structural variation may be an underappreciated mutational process in cancer genomes.

Cycling 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 sub-nuclear position. Moreover, RT is regulated during development and is altered in disease. Exploring mechanisms linking RT to other cellular processes in normal and diseased cells will be facilitated by rapid and robust methods with which to measure RT genome wide. Here, we describe a rapid, robust and relatively inexpensive protocol to analyze genome-wide RT by next-generation sequencing (NGS). This protocol yields highly reproducible results across laboratories and platforms. We also provide computational pipelines for analysis, parsing phased genomes using single nucleotide polymorphisms (SNP) for analyzing RT allelic asynchrony, and for direct comparison to Repli-chip data obtained by analyzing nascent DNA by microarrays.

DNA replication occurs in a defined temporal order during S phase, known as the replication timing programme, which is regulated not only during the cell cycle but also during the process of development and differentiation. The units of replication timing regulation, known as replication domains (RDs), frequently comprise several nearly synchronously firing replication origins. Replication domains correspond to topologically associating domains (TADs) mapped by chromatin conformation capture methods and are likely to be the molecular equivalents of replication foci observed using cytogenetic methods. Both TAD and replication foci are considered to be stable structural units of chromosomes, conserved through the cell cycle and development, and accordingly, the boundaries of RDs also appear to be stable in different cell types. During both normal development and progression of disease, distinct cell states are characterized by unique replication timing signatures, with approximately half of genomic RDs switching replication timing between these cell states. Advances in functional genomics provide hope that we can soon gain an understanding of the cause and consequence of the replication timing programme and its myriad correlations with chromatin context and gene regulation.

In metazoans, nuclear DNA is organized during the interphase in negatively supercoiled loops anchored to a compartment or substructure known as the nuclear matrix. The interactions between DNA and the nuclear matrix (NM) are of higher affinity than those between DNA and chromatin proteins since the last ones do not resist the procedures for extracting the NM. The structural interactions DNA-NM constitute a set of topological relationships that define a nuclear higher order structure (NHOS) although there are further higher order levels of organization within the nucleus. So far, the evidence derived from studies with primary hepatocytes and naïve B lymphocytes indicates that the NHOS is cell-type specific at the local and at the large-scale level, and so it has been suggested that such NHOS is primary determined by structural and thermodynamic constraints. We carried out a comparative characterization of the NHOS of postmitotic cortical neurons with that of hepatocytes and naïve B lymphocytes. Our results indicate that the NHOS of neurons is completely different at the large scale and at the local level from that one observed in hepatocytes or in naïve B lymphocytes, confirming on the one hand that the set of structural DNA-NM interactions is cell-type specific and supporting, on the other hand the notion that structural constraints that impinge on chromosomal DNA and the NM are more important for determining this NHOS than functional constraints related to replication and/or transcription. J. Cell. Biochem. 118: 2151-2160, 2017. © 2016 Wiley Periodicals, Inc.

Nuclear DNA of metazoans is organized in supercoiled loops anchored to a proteinaceous substructure known as the nuclear matrix (NM). DNA is anchored to the NM by non-coding sequences known as matrix attachment regions (MARs). There are no consensus sequences for identification of MARs and not all potential MARs are actually bound to the NM constituting loop attachment regions (LARs). Fundamental processes of nuclear physiology occur at macromolecular complexes organized on the NM; thus, the topological organization of DNA loops must be important. Here, we describe a general method for determining the structural DNA loop organization in any large genomic region with a known sequence. The method exploits the topological properties of loop DNA attached to the NM and elementary topological principles such as that points in a deformable string (DNA) can be positionally mapped relative to a position-reference invariant (NM), and from such mapping, the configuration of the string in third dimension can be deduced. Therefore, it is possible to determine the specific DNA loop configuration without previous characterization of the LARs involved. We determined in hepatocytes and B-lymphocytes of the rat the DNA loop organization of a genomic region that contains four members of the albumin gene family.

La replicación del dna es indispensable para la transmisión de la información genética y permite copiar el genoma con gran exactitud. Desde el siglo pasado se propuso el modelo del replicón para explicar el mecanismo general de duplicación del genoma en bacterias. Estudios posteriores en la levadura permitieron identificar proteínas y secuencias de dna que participan en el inicio de la replicación en forma similar a lo descrito en procariontes, esto condujo a intentar generalizar el modelo del replicón a los eucariontes. Se han descrito algunos factores clave en el proceso de replicación que están conservados desde la levadura hasta el humano. Sin embargo, todavía no se comprende cómo se determinan los sitios de inicio de la replicación y cuál es la estructura del replicón en los metazoarios. En este artículo se sugiere que la organización topológica del dna en el núcleo celular determina la estructura y función de los replicones en los eucariontes superiores.

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.

Duplication of the genome in mammalian cells occurs in a defined temporal order referred to as its replication-timing (RT) program. RT changes dynamically during development, regulated in units of 400–800 kb referred to as replication domains (RDs). Changes in RT are generally coordinated with transcriptional competence and changes in subnuclear position. We generated genome-wide RT profiles for 26 distinct human cell types, including embryonic stem cell (hESC)-derived, primary cells and established cell lines representing intermediate stages of endoderm, mesoderm, ectoderm, and neural crest (NC) development. We identified clusters of RDs that replicate at unique times in each stage (RT signatures) and confirmed global consolidation of the genome into larger synchronously replicating segments during differentiation. Surprisingly, transcriptome data revealed that the well-accepted correlation between early replication and transcriptional activity was restricted to RT-constitutive genes, whereas two-thirds of the genes that switched RT during differentiation were strongly expressed when late replicating in one or more cell types. Closer inspection revealed that transcription of this class of genes was frequently restricted to the lineage in which the RT switch occurred, but was induced prior to a late-to-early RT switch and/or down-regulated after an early-to-late RT switch. Analysis of transcriptional regulatory networks showed that this class of genes contains strong regulators of genes that were only expressed when early replicating. These results provide intriguing new insight into the complex relationship between transcription and RT regulation during human development.

Genome-wide DNA replication timing (RT) profiles reflect the global 3D chromosome architecture of cells. They also provide a comprehensive and unique megabase-scale picture of the cellular epigenetic state. Thus normal differentiation involves reproducible changes in RT and transformation generally perturbs these, although the potential effects of altered RT on the properties of transformed cells remain largely unknown. A major challenge to interrogating these issues in human acute lymphoid leukemia (ALL) is the low proliferative activity of most of the cells, which may be further reduced in cryopreserved samples and difficult to overcome in vitro. In contrast, the ability of many human ALL cell populations to expand when transplanted in highly immunodeficient mice is well documented. To examine the stability of DNA RT profiles of serially passaged xenografts of primary human B- and T-ALL cells, we first devised a method that circumvents the need for BrdU incorporation to distinguish early versus late S-phase cells. Using this and more standard protocols, we found consistent strong retention in xenografts of the original patient-specific RT features, for all 8 primary human ALL cases surveyed (7 B-ALLs and one T-ALL). Moreover, in a case where genomic analyses indicated changing subclonal dynamics in serial passages, the RT profiles tracked concordantly. These results show that DNA RT is a relatively stable feature of human ALLs propagated in immunodeficient mice. In addition, they suggest the power of this approach for future interrogation of the origin and consequences of altered DNA RT in these diseases.

La replicación del dna es indispensable para la transmisión de la información genética y permite copiar el genoma con gran exactitud. Desde el siglo pasado se propuso el modelo del replicón para explicar el mecanismo general de duplicación del genoma en bacterias. Estudios posteriores en la levadura permitieron identificar proteínas y secuencias de dna que participan en el inicio de la replicación en forma similar a lo descrito en procariontes, esto condujo a intentar generalizar el modelo del replicón a los eucariontes. Se han descrito algunos factores clave en el proceso de replicación que están conservados desde la levadura hasta el humano. Sin embargo, todavía no se comprende cómo se determinan los sitios de inicio de la replicación y cuál es la estructura del replicón en los metazoarios. En este artículo se sugiere que la organización topológica del dna en el núcleo celular determina la estructura y función de los replicones en los eucariontes superiores.

Discovery of lineage-specific somatic copy number variation (CNV) in mammals has led to debate over whether CNVs are mutations that propagate disease or whether they are a normal, and even essential, aspect of cell biology. We show that 1,000N polyploid trophoblast giant cells (TGCs) of the mouse placenta contain 47 regions, totaling 138 Megabases, where genomic copies are underrepresented (UR). UR domains originate from a subset of late-replicating heterochromatic regions containing gene deserts and genes involved in cell adhesion and neurogenesis. While lineage-specific CNVs have been identified in mammalian cells, classically in the immune system where V(D)J recombination occurs, we demonstrate that CNVs form during gestation in the placenta by an underreplication mechanism, not by recombination nor deletion. Our results reveal that large scale CNVs are a normal feature of the mammalian placental genome, which are regulated systematically during embryogenesis and are propagated by a mechanism of underreplication., Author Summary Generally, every mammalian cell has the same complement of each part of its genome. However, copy number variation (CNV) can occur, where, compared to the rest of its genome, a cell has either more or less of a specific genomic region. It is unknown whether CNVs cause disease, or whether they are a normal aspect of cell biology. We investigated CNVs in polyploid trophoblast giant cells (TGCs) of the mouse placenta, which have up to 1,000 copies of the genome in each cell. We found that there are 47 regions with decreased copy number in TGCs, which we call underrepresented (UR) domains. These domains are marked in the TGC progenitor cells and we suggest that they gradually form during gestation due to slow replication versus fast replication of the rest of the genome. While UR domains contain cell adhesion and neuronal genes, they also contain significantly fewer genes than other genomic regions. Our results demonstrate that CNVs are a normal feature of the mammalian placental genome, which are regulated systematically during pregnancy.

Background In the interphase nucleus of metazoan cells DNA is organized in supercoiled loops anchored to a nuclear matrix (NM). There is varied evidence indicating that DNA replication occurs in replication factories organized upon the NM and that DNA loops may correspond to the actual replicons in vivo. In normal rat liver the hepatocytes are arrested in G0 but they synchronously re-enter the cell cycle after partial-hepatectomy leading to liver regeneration in vivo. We have previously determined in quiescent rat hepatocytes that a 162 kbp genomic region containing members of the albumin gene family is organized into five structural DNA loops. Results In the present work we tracked down the movement relative to the NM of DNA sequences located at different points within such five structural DNA loops during the S phase and after the return to cellular quiescence during liver regeneration. Our results indicate that looped DNA moves sequentially towards the NM during replication and then returns to its original position in newly quiescent cells, once the liver regeneration has been achieved. Conclusions Looped DNA moves in a sequential fashion, as if reeled in, towards the NM during DNA replication in vivo thus supporting the notion that the DNA template is pulled progressively towards the replication factories on the NM so as to be replicated. These results provide further evidence that the structural DNA loops correspond to the actual replicons in vivo.

s / Placenta 34 (2013) A1–A99 A8 migrants to high altitude. Our global hypothesis is that evolution has favored genes that may permit greater fetal growth even in the presence of hypoxia. Methods: Using microarrays, we compared placental RNA from 45 Native Americans versus Europeans. Our criteria for differentially expressed genes (DEGs) were a fold change > 1.25 with an adjusted p-value < 0.01. Results: There were only 10 DEGs in Andeans versus Europeans at 400 m. The most pronounced was localized to chromosome 19, is a member of the pregnancy specific glycoprotein family, and its expression is increased in Andeans. In contrast, there were 54 altitude-induced DEGs in the Europeans and 40 in the Andeans. These are ancestry dependent DEGs because there were only two shared between ancestry groups (Figure 1). There were 103 DEGs between the high and low altitude individuals irrespective of ancestry, including NOS3, a nitric oxide synthase, associated with the regulation of oxygen in hypoxic environments. Conclusion: These findings support the concept that the ancestry-environment interaction differs between Native American and European groups. Andeans at 3600 m showed enrichment of genes related to the HIF1 pathway, while there is over-representation of DEGs related to inflammation and coagulation in Europeans. Our current analyses target differences in placental morphometrics in relation to the DEGs reported here. Support: 1R211HD068954-01, NIHHD042737, NSFBCS0309142

Biological networks inherently have uncertain topologies. This arises from many factors. For instance, interactions between molecules may or may not take place under varying conditions. Genetic or epigenetic mutations may also alter biological processes like transcription or translation. This uncertainty is often modeled by associating each interaction with a probability value. Studying biological networks under this probabilistic model has already been shown to yield accurate and insightful analysis of interaction data. However, the problem of assigning accurate probability values to interactions remains unresolved. In this paper, we present a novel method for computing interaction probabilities in signaling networks based on transcription levels of genes. The transcription levels define the signal reachability probability between membrane receptors and transcription factors. Our method computes the interaction probabilities that minimize the gap between the observed and the computed signal reachability probabilities. We evaluate our method on four signaling networks from the Kyoto Encyclopedia of Genes and Genomes (KEGG). For each network, we compute its edge probabilities using the gene expression profiles for seven major leukemia subtypes. We use these values to analyze how the stress induced by different leukemia subtypes affects signaling interactions.