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

ASAR lncRNAs control DNA replication timing through interactions with multiple hnRNP/RNA binding proteins

ASARs are a family of very-long noncoding RNAs that control replication timing on individual human autosomes, and are essential for chromosome stability. The eight known ASAR lncRNAs remain closely associated with their parent chromosomes. Analysis of RNA-protein interaction data (from ENCODE) revealed numerous RBPs with significant interactions with multiple ASAR lncRNAs, with several hnRNPs as abundant interactors. An ~7 kb domain within the ASAR6-141 lncRNA shows a striking density of RBP interaction sites. Genetic deletion and ectopic integration assays indicate that this ~7 kb RNA binding protein domain contains functional sequences for controlling replication timing of entire chromosomes in cis. shRNA-mediated depletion of 10 different RNA binding proteins, including HNRNPA1, HNRNPC, HNRNPL, HNRNPM, HNRNPU, or HNRNPUL1, results in dissociation of ASAR lncRNAs from their chromosome territories, and disrupts the synchronous replication that occurs on all autosome pairs, recapitulating the effect of individual ASAR knockouts on a genome-wide scale. Our results further demonstrate the role that ASARs play during the temporal order of genome-wide replication, and we propose that ASARs function as essential RNA scaffolds for the assembly of hnRNP complexes that help maintain the structural integrity of each mammalian chromosome.

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, mediated by Polycomb repression. We show that histone genes are expressed with allelic imbalance in mESCs, are involved in haplotype-specific chromatin contact marked by H3K27me3, and are targets of Polycomb repression through conditional knockouts of Ezh2 or Ring1b. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

A rewiring of DNA replication mediated by MRE11 exonuclease underlies primed-to-naive cell de-differentiation

Mouse embryonic stem cells (mESCs) in the primed pluripotency state, which resembles the post-implantation epiblast, can be de-differentiated in culture to a naive state that resembles the pre-implantation inner cell mass. We report that primed-to-naive mESC transition entails a significant slowdown of DNA replication forks and the compensatory activation of dormant origins. Using isolation of proteins on nascent DNA coupled to mass spectrometry, we identify key changes in replisome composition that are responsible for these effects. Naive mESC forks are enriched in MRE11 nuclease and other DNA repair proteins. MRE11 is recruited to newly synthesized DNA in response to transcription-replication conflicts, and its inhibition or genetic downregulation in naive mESCs is sufficient to restore the fork rate of primed cells. Transcriptomic analyses indicate that MRE11 exonuclease activity is required for the complete primed-to-naive mESC transition, demonstrating a direct link between DNA replication dynamics and the mESC de-differentiation process.

Stem cell transcriptional profiles from mouse subspecies reveal cis-regulatory evolution at translation genes

A key goal of evolutionary genomics is to harness molecular data to draw inferences about selective forces that have acted on genomes. The field progresses in large part through the development of advanced molecular-evolution analysis methods. Here we explored the intersection between classical sequence-based tests for selection and an empirical expression-based approach, using stem cells from Mus musculus subspecies as a model. Using a test of directional, cis-regulatory evolution across genes in pathways, we discovered a unique program of induction of translation genes in stem cells of the Southeast Asian mouse M. m. castaneus relative to its sister taxa. As a complement, we used sequence analyses to find population-genomic signatures of selection in M. m. castaneus, at the upstream regions of the translation genes, including at transcription factor binding sites. We interpret our data under a model of changes in lineage-specific pressures across Mus musculus in stem cells with high translational capacity. Together, our findings underscore the rigor of integrating expression and sequence-based methods to generate hypotheses about evolutionary events from long ago.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Emergence of replication timing during early mammalian development

DNA replication enables genetic inheritance across the kingdoms of life. Replication occurs with a defined temporal order known as the replication timing (RT) programme, leading to organization of the genome into early- or late-replicating regions. RT is cell-type specific, is tightly linked to the three-dimensional nuclear organization of the genome1,2 and is considered an epigenetic fingerprint3. In spite of its importance in maintaining the epigenome4, the developmental regulation of RT in mammals in vivo has not been explored. Here, using single-cell Repli-seq5, we generated genome-wide RT maps of mouse embryos from the zygote to the blastocyst stage. Our data show that RT is initially not well defined but becomes defined progressively from the 4-cell stage, coinciding with strengthening of the A and B compartments. We show that transcription contributes to the precision of the RT programme and that the difference in RT between the A and B compartments depends on RNA polymerase II at zygotic genome activation. Our data indicate that the establishment of nuclear organization precedes the acquisition of defined RT features and primes the partitioning of the genome into early- and late-replicating domains. Our work sheds light on the establishment of the epigenome at the beginning of mammalian development and reveals the organizing principles of genome organization.

Quantification of absolute transcription factor binding affinities in the native chromatin context using BANC-seq

Transcription factor binding across the genome is regulated by DNA sequence and chromatin features. However, it is not yet possible to quantify the impact of chromatin context on transcription factor binding affinities. Here, we report a method called binding affinities to native chromatin by sequencing (BANC-seq) to determine absolute apparent binding affinities of transcription factors to native DNA across the genome. In BANC-seq, a concentration range of a tagged transcription factor is added to isolated nuclei. Concentration-dependent binding is then measured per sample to quantify apparent binding affinities across the genome. BANC-seq adds a quantitative dimension to transcription factor biology, which enables stratification of genomic targets based on transcription factor concentration and prediction of transcription factor binding sites under non-physiological conditions, such as disease-associated overexpression of (onco)genes. Notably, whereas consensus DNA binding motifs for transcription factors are important to establish high-affinity binding sites, these motifs are not always strictly required to generate nanomolar-affinity interactions in the genome.

Single-cell DNA replication dynamics in genomically unstable cancers

Dysregulated DNA replication is both a cause and a consequence of aneuploidy, yet the dynamics of DNA replication in aneuploid cell populations remains understudied. We developed a new method, PERT, for inferring cell-specific DNA replication states from single-cell whole genome sequencing, and investigated clone-specific DNA replication dynamics in >50,000 cells obtained from a collection of aneuploid and clonally heterogeneous cell lines, xenografts and primary cancer tissues. Clone replication timing (RT) profiles correlated with future copy number changes in serially passaged cell lines. Cell type was the strongest determinant of RT heterogeneity, while whole genome doubling and mutational process were associated with accumulation of late S-phase cells and weaker RT associations. Copy number changes affecting chromosome X had striking impact on RT, with loss of the inactive X allele shifting replication earlier, and loss of inactive Xq resulting in reactivation of Xp. Finally, analysis of time series xenografts illustrate how cell cycle distributions approximate clone proliferation, recapitulating expected relationships between proliferation and fitness in treatment-naive and chemotherapeutic contexts.

Allele pairing at Sun1-enriched domains at the nuclear periphery via T1A3 tandem DNA repeats

Spatiotemporal gene regulation is fundamental to the biology of diploid cells. Therefore, effective communication between two alleles and their geometry in the nucleus is important. However, the mechanism that fine-tunes the expression from each of the two alleles of an autosome is enigmatic. Establishing an allele-specific gene expression visualization system in living cells, we show that alleles of biallelically expressed Cth and Ttc4 genes are paired prior to acquiring monoallelic expression. We found that active alleles of monoallelic genes are preferentially localized at Sun1-enriched domains at the nuclear periphery. These peripherally localized active DNA loci are enriched with adenine-thymidine-rich tandem repeats that interact with Hnrnpd and reside in a Hi-C-defined A compartment within the B compartment. Our results demonstrate the biological significance of T ₁ A ₃ tandem repeat sequences in genome organization and how the regulation of gene expression, at the level of individual alleles, relates to their spatial arrangement.

Conserved noncoding sequences correlate with distant gene contacts in Arabidopsis and Brassica

Physical contact between genes distant on chromosomes is a potentially important way for genes to coordinate their expressions. To investigate the potential importance of distant contacts, we performed high-throughput chromatin conformation capture (Hi-C) experiments on leaf nuclei isolated from Brassica rapa and Brassica oleracea. We then combined our results with published Hi-C data from Arabidopsis thaliana. We found that distant genes come into physical contact and do so preferentially between the proximal promoter of one gene and the downstream region of another gene. Genes with higher numbers of conserved noncoding sequences (CNSs) nearby were more likely to have contact with distant genes. With more CNSs came higher numbers of transcription factor binding sites and more histone modifications associated with the activity. In addition, for the genes we studied, distant contacting genes with CNSs were more likely to be transcriptionally coordinated. These observations suggest that CNSs may enrich active histone modifications and recruit transcription factors, correlating with distant contacts to ensure coordinated expression. This study advances our knowledge of gene contacts and provides insights into the relationship between CNSs and distant gene contacts in plants.

Optimized Repli-seq: improved DNA replication timing analysis by next-generation sequencing

The human genome is divided into functional units that replicate at specific times during S-phase. This temporal program is known as replication timing (RT) and is coordinated with the spatial organization of the genome and transcriptional activity. RT is also cell type-specific, dynamically regulated during development, and alterations in RT are observed in multiple diseases. Thus, the precise measure of RT is critical to understand the role of RT in gene function regulation. Distinct methods for assaying the RT program exist; however, conventional methods require thousands of cells as input, prohibiting its applicability to samples with limited cell numbers such as those from disease patients or from early developing embryos. Although single-cell RT analyses have been developed, these methods are low throughput, require generation of numerous libraries, increased sequencing costs, and produce low resolution data. Here, we developed an improved method to measure RT genome-wide that enables high-resolution analysis of low input samples. This method incorporates direct cell sorting into lysis buffer, as well as DNA fragmentation and library preparation in a single tube, resulting in higher yields, increased quality, and reproducibility with decreased costs. We also performed a systematic data processing analysis to provide standardized parameters for RT measurement. This optimized method facilitates RT analysis and will enable its application to a broad range of studies investigating the role of RT in gene expression, nuclear architecture, and disease.

Epigenetic control of chromosome-associated lncRNA genes essential for replication and stability

ASARs are long noncoding RNA genes that control replication timing of entire human chromosomes in cis. The three known ASAR genes are located on human chromosomes 6 and 15, and are essential for chromosome integrity. To identify ASARs on all human chromosomes we utilize a set of distinctive ASAR characteristics that allow for the identification of hundreds of autosomal loci with epigenetically controlled, allele-restricted behavior in expression and replication timing of coding and noncoding genes, and is distinct from genomic imprinting. Disruption of noncoding RNA genes at five of five tested loci result in chromosome-wide delayed replication and chromosomal instability, validating their ASAR activity. In addition to the three known essential cis-acting chromosomal loci, origins, centromeres, and telomeres, we propose that all mammalian chromosomes also contain "Inactivation/Stability Centers" that display allele-restricted epigenetic regulation of protein coding and noncoding ASAR genes that are essential for replication and stability of each chromosome.

Allele-specific differential regulation of monoallelically expressed autosomal genes in the cardiac lineage

Each mammalian autosomal gene is represented by two alleles in diploid cells. To our knowledge, no insights have been made in regard to allele-specific regulatory mechanisms of autosomes. Here we use allele-specific single cell transcriptomic analysis to elucidate the establishment of monoallelic gene expression in the cardiac lineage. We find that monoallelically expressed autosomal genes in mESCs and mouse blastocyst cells are differentially regulated based on the genetic background of the parental alleles. However, the genetic background of the allele does not affect the establishment of monoallelic genes in differentiated cardiomyocytes. Additionally, we observe epigenetic differences between deterministic and random autosomal monoallelic genes. Moreover, we also find a greater contribution of the maternal versus paternal allele to the development and homeostasis of cardiac tissue and in cardiac health, highlighting the importance of maternal influence in male cardiac tissue homeostasis. Our findings emphasize the significance of allele-specific insights into gene regulation in development, homeostasis and disease.

Telomere-to-telomere human DNA replication timing profiles

The spatiotemporal organization of DNA replication produces a highly robust and reproducible replication timing profile. Sequencing-based methods for assaying replication timing genome-wide have become commonplace, but regions of high repeat content in the human genome have remained refractory to analysis. Here, we report the first nearly-gapless telomere-to-telomere replication timing profiles in human, using the T2T-CHM13 genome assembly and sequencing data for five cell lines. We find that replication timing can be successfully assayed in centromeres and large blocks of heterochromatin. Centromeric regions replicate in mid-to-late S-phase and contain replication-timing peaks at a similar density to other genomic regions, while distinct families of heterochromatic satellite DNA differ in their bias for replicating in late S-phase. The high degree of consistency in centromeric replication timing across chromosomes within each cell line prompts further investigation into the mechanisms dictating that some cell lines replicate their centromeres earlier than others, and what the consequences of this variation are.

Screening thousands of transcribed coding and non-coding regions reveals sequence determinants of RNA polymerase II elongation potential

Precise regulation of transcription by RNA polymerase II (RNAPII) is critical for organismal growth and development. However, what determines whether an engaged RNAPII will synthesize a full-length transcript or terminate prematurely is poorly understood. Notably, RNAPII is far more susceptible to termination when transcribing non-coding RNAs than when synthesizing protein-coding mRNAs, but the mechanisms underlying this are unclear. To investigate the impact of transcribed sequence on elongation potential, we developed a method to screen the effects of thousands of INtegrated Sequences on Expression of RNA and Translation using high-throughput sequencing (INSERT-seq). We found that higher AT content in non-coding RNAs, rather than specific sequence motifs, drives RNAPII termination. Further, we demonstrate that 5' splice sites autonomously stimulate processive transcription, even in the absence of polyadenylation signals. Our results reveal a potent role for the transcribed sequence in dictating gene output and demonstrate the power of INSERT-seq toward illuminating these contributions.

Cell cycle dynamics and developmental dynamics of the 3D genome: toward linking the two timescales

In the mammalian cell nucleus, chromosomes are folded differently in interphase and mitosis. Interphase chromosomes are relatively decondensed and display at least two unique layers of higher-order organization: topologically associating domains (TADs) and cell-type-specific A/B compartments, which correlate well with early/late DNA replication timing (RT). In mitosis, these structures rapidly disappear but are gradually reconstructed during G1 phase, coincident with the establishment of the RT program. However, these structures also change dynamically during cell differentiation and reprogramming, and yet we are surprisingly ignorant about the relationship between their cell cycle dynamics and developmental dynamics. In this review, we summarize the recent findings on this topic, discuss how these two processes might be coordinated with each other and its potential significance.

Asynchronous Replication Timing: A Mechanism for Monoallelic Choice During Development

Developmental programming is carried out by a sequence of molecular choices that epigenetically mark the genome to generate the stable cell types which make up the total organism. A number of important processes, such as genomic imprinting, selection of immune or olfactory receptors, and X-chromosome inactivation in females are dependent on the ability to stably choose one single allele in each cell. In this perspective, we propose that asynchronous replication timing (ASRT) serves as the basis for a sophisticated universal mechanism for mediating and maintaining these decisions.

Whole chromosome loss and genomic instability in mouse embryos after CRISPR-Cas9 genome editing

Karyotype alterations have emerged as on-target complications from CRISPR-Cas9 genome editing. However, the events that lead to these karyotypic changes in embryos after Cas9-treatment remain unknown. Here, using imaging and single-cell genome sequencing of 8-cell stage embryos, we track both spontaneous and Cas9-induced karyotype aberrations through the first three divisions of embryonic development. We observe the generation of abnormal structures of the nucleus that arise as a consequence of errors in mitosis, including micronuclei and chromosome bridges, and determine their contribution to common karyotype aberrations including whole chromosome loss that has been recently reported after editing in embryos. Together, these data demonstrate that Cas9-mediated germline genome editing can lead to unwanted on-target side effects, including major chromosome structural alterations that can be propagated over several divisions of embryonic development.

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

Uncovering Divergence in Gene Expression Regulation in the Adaptation of Yeast to Nitrogen Scarcity

Saccharomyces cerevisiae rewires its transcriptional output to survive stressful environments, such as nitrogen scarcity under fermentative conditions. Although divergence in nitrogen metabolism among natural yeast populations has been reported, the impact of regulatory genetic variants modulating gene expression and nitrogen consumption remains to be investigated. Here, we employed an F1 hybrid from two contrasting S. cerevisiae strains, providing a controlled genetic environment to map cis factors involved in the divergence of gene expression regulation in response to nitrogen scarcity. We used a dual approach to obtain genome-wide allele-specific profiles of chromatin accessibility, transcription factor binding, and gene expression through ATAC-seq (assay for transposase accessible chromatin) and RNA-seq (transcriptome sequencing). We observed large variability in allele-specific expression and accessibility between the two genetic backgrounds, with a third of these differences specific to a deficient nitrogen environment. Furthermore, we discovered events of allelic bias in gene expression correlating with allelic bias in transcription factor binding solely under nitrogen scarcity, where the majority of these transcription factors orchestrates the nitrogen catabolite repression regulatory pathway and demonstrates a cis × environment-specific response. Our approach allowed us to find cis variants modulating gene expression, chromatin accessibility, and allelic differences in transcription factor binding in response to low nitrogen culture conditions. IMPORTANCE Historically, coding variants were prioritized when searching for causal mechanisms driving adaptation of natural populations to stressful environments. However, the recent focus on noncoding variants demonstrated their ubiquitous role in adaptation. Here, we performed genome-wide regulatory variation profiles between two divergent yeast strains when facing nitrogen nutritional stress. The open chromatin availability of several regulatory regions changes in response to nitrogen scarcity. Importantly, we describe regulatory events that deviate between strains. Our results demonstrate a widespread variation in gene expression regulation between naturally occurring populations in response to stressful environments.

Understanding three-dimensional chromatin organization in diploid genomes

The three-dimensional (3D) organization of chromatin in the nucleus of diploid eukaryotic organisms has fascinated biologists for many years. Despite major progress in chromatin conformation studies, current knowledge regarding the spatial organization of diploid (maternal and paternal) genomes is still limited. Recent advances in Hi-C technology and data processing approaches have enabled construction of diploid Hi-C contact maps. These maps greatly accelerated the pace of novel discoveries in haplotype-resolved 3D genome studies, revealing the role of allele biased chromatin conformation in transcriptional regulation. Here, we review emerging concepts and haplotype phasing strategies of Hi-C data in 3D diploid genome studies. We discuss new insights on homologous chromosomal organization and the interplay between allelic biased chromatin architecture and several nuclear functions, explaining how haplotype-resolved Hi-C technologies have been used to resolve important biological questions.

GAMIBHEAR: whole-genome haplotype reconstruction from Genome Architecture Mapping data

Motivation

Genome Architecture Mapping (GAM) was recently introduced as a digestion- and ligation-free method to detect chromatin conformation. Orthogonal to existing approaches based on chromatin conformation capture (3C), GAM’s ability to capture both inter- and intra-chromosomal contacts from low amounts of input data makes it particularly well suited for allele-specific analyses in a clinical setting. Allele-specific analyses are powerful tools to investigate the effects of genetic variants on many cellular phenotypes including chromatin conformation, but require the haplotypes of the individuals under study to be known a priori. So far, however, no algorithm exists for haplotype reconstruction and phasing of genetic variants from GAM data, hindering the allele-specific analysis of chromatin contact points in non-model organisms or individuals with unknown haplotypes.

Results

We present GAMIBHEAR, a tool for accurate haplotype reconstruction from GAM data. GAMIBHEAR aggregates allelic co-observation frequencies from GAM data and employs a GAM-specific probabilistic model of haplotype capture to optimize phasing accuracy. Using a hybrid mouse embryonic stem cell line with known haplotype structure as a benchmark dataset, we assess correctness and completeness of the reconstructed haplotypes, and demonstrate the power of GAMIBHEAR to infer accurate genome-wide haplotypes from GAM data.

Availability and implementation

GAMIBHEAR is available as an R package under the open-source GPL-2 license at https://bitbucket.org/schwarzlab/gamibhear.

Supplementary information

Supplementary data are available at Bioinformatics online.

Chromosomal coordination and differential structure of asynchronous replicating regions

Stochastic asynchronous replication timing (AS-RT) is a phenomenon in which the time of replication of each allele is different, and the identity of the early allele varies between cells. By taking advantage of stable clonal pre-B cell populations derived from C57BL6/Castaneous mice, we have mapped the genome-wide AS-RT loci, independently of genetic differences. These regions are characterized by differential chromatin accessibility, mono-allelic expression and include new gene families involved in specifying cell identity. By combining population level mapping with single cell FISH, our data reveal the existence of a novel regulatory program that coordinates a fixed relationship between AS-RT regions on any given chromosome, with some loci set to replicate in a parallel and others set in the anti-parallel orientation. Our results show that AS-RT is a highly regulated epigenetic mark established during early embryogenesis that may be used for facilitating the programming of mono-allelic choice throughout development.

Most regions of the mammalian genome replicate both alleles in a synchronous manner, but some loci have been found to replicate asynchronously and the time of replication of each allele is different. Here the authors, by employing clonal mouse cells from a hybrid strain chart replication timing over the entire genome, using polymorphisms to distinguish between the paternal and maternal alleles.

Formation of a multi-layered 3-dimensional structure of the heterochromatin compartment during early mammalian development

During embryogenesis in mammals, the 3-dimensional (3D) genome organization changes globally in parallel with transcription changes in a cell-type specific manner. This involves the progressive formation of heterochromatin, the best example of which is the inactive X chromosome (Xi) in females, originally discovered as a compact 3D structure at the nuclear periphery known as the Barr body. The heterochromatin formation on the autosomes and the Xi is tightly associated with the differentiation state and the developmental potential of cells, making it an ideal readout of the cellular epigenetic state. At a glance, the heterochromatin appears to be uniform. However, recent studies are beginning to reveal a more complex picture, with multiple hierarchical levels co-existing within the heterochromatin compartment. Such hierarchical levels appear to exist in the heterochromatin compartment on autosomes as well as on the Xi. Here, we review recent progress in our understanding of the 3D genome organization changes during the period of differentiation surrounding pluripotency in vivo and in vitro, with a focus on the heterochromatin compartment. We first look at the whole genome, then focus on the Xi, and discuss their differences and similarities. Finally, we present a unified view of how the heterochromatin compartment is formed and regulated during early development. In particular, we emphasize that there are multiple layers within the heterochromatic compartment on both the autosomes and the Xi, with regulatory mechanisms common and specific to each layer.

CTCF-Mediated Genome Architecture Regulates the Dosage of Mitotically Stable Mono-allelic Expression of Autosomal Genes

The mechanisms that guide the clonally stable random mono-allelic expression of autosomal genes remain enigmatic. We show that (1) mono-allelically expressed (MAE) genes are assorted and insulated from bi-allelically expressed (BAE) genes through CTCF-mediated chromatin loops; (2) the cell-type-specific dynamics of mono-allelic expression coincides with the gain and loss of chromatin insulator sites; (3) dosage of MAE genes is more sensitive to the loss of chromatin insulation than that of BAE genes; and (4) inactive alleles of MAE genes are significantly more insulated than active alleles and are de-repressed upon CTCF depletion. This alludes to a topology wherein the inactive alleles of MAE genes are insulated from the spatial interference of transcriptional states from the neighboring bi-allelic domains via CTCF-mediated loops. We propose that CTCF functions as a typical insulator on inactive alleles, but facilitates transcription through enhancer-linking on active allele of MAE genes, indicating widespread allele-specific regulatory roles of CTCF.

Distinct features of nucleolus-associated domains in mouse embryonic stem cells

Heterochromatin in eukaryotic interphase cells frequently localizes to the nucleolar periphery (nucleolus-associated domains (NADs)) and the nuclear lamina (lamina-associated domains (LADs)). Gene expression in somatic cell NADs is generally low, but NADs have not been characterized in mammalian stem cells. Here, we generated the first genome-wide map of NADs in mouse embryonic stem cells (mESCs) via deep sequencing of chromatin associated with biochemically purified nucleoli. As we had observed in mouse embryonic fibroblasts (MEFs), the large type I subset of NADs overlaps with constitutive LADs and is enriched for features of constitutive heterochromatin, including late replication timing and low gene density and expression levels. Conversely, the type II NAD subset overlaps with loci that are not lamina-associated, but in mESCs, type II NADs are much less abundant than in MEFs. mESC NADs are also much less enriched in H3K27me3 modified regions than are NADs in MEFs. Additionally, comparision of MEF and mESC NADs revealed enrichment of developmentally regulated genes in cell-type-specific NADs. Together, these data indicate that NADs are a developmentally dynamic component of heterochromatin. These studies implicate association with the nucleolar periphery as a mechanism for developmentally regulated gene expression and will facilitate future studies of NADs during mESC differentiation.

High-resolution Repli-Seq defines the temporal choreography of initiation, elongation and termination of replication in mammalian cells

BACKGROUND

DNA replication in mammalian cells occurs in a defined temporal order during S phase, known as the replication timing (RT) programme. Replication timing is developmentally regulated and correlated with chromatin conformation and local transcriptional potential. Here, we present RT profiles of unprecedented temporal resolution in two human embryonic stem cell lines, human colon carcinoma line HCT116, and mouse embryonic stem cells and their neural progenitor derivatives.

RESULTS

Fine temporal windows revealed a remarkable degree of cell-to-cell conservation in RT, particularly at the very beginning and ends of S phase, and identified 5 temporal patterns of replication in all cell types, consistent with varying degrees of initiation efficiency. Zones of replication initiation (IZs) were detected throughout S phase and interacted in 3D space preferentially with other IZs of similar firing time. Temporal transition regions were resolved into segments of uni-directional replication punctuated at specific sites by small, inefficient IZs. Sites of convergent replication were divided into sites of termination or large constant timing regions consisting of many synchronous IZs in tandem. Developmental transitions in RT occured mainly by activating or inactivating individual IZs or occasionally by altering IZ firing time, demonstrating that IZs, rather than individual origins, are the units of developmental regulation. Finally, haplotype phasing revealed numerous regions of allele-specific and allele-independent asynchronous replication. Allele-independent asynchronous replication was correlated with the presence of previously mapped common fragile sites.

CONCLUSIONS

Altogether, these data provide a detailed temporal choreography of DNA replication in mammalian cells.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

The Histone Code of Senescence

Senescence is the end point of a complex cellular response that proceeds through a set of highly regulated steps. Initially, the permanent cell-cycle arrest that characterizes senescence is a pro-survival response to irreparable DNA damage. The maintenance of this prolonged condition requires the adaptation of the cells to an unfavorable, demanding and stressful microenvironment. This adaptation is orchestrated through a deep epigenetic resetting. A first wave of epigenetic changes builds a dam on irreparable DNA damage and sustains the pro-survival response and the cell-cycle arrest. Later on, a second wave of epigenetic modifications allows the genomic reorganization to sustain the transcription of pro-inflammatory genes. The balanced epigenetic dynamism of senescent cells influences physiological processes, such as differentiation, embryogenesis and aging, while its alteration leads to cancer, neurodegeneration and premature aging. Here we provide an overview of the most relevant histone modifications, which characterize senescence, aging and the activation of a prolonged DNA damage response.

Integrative characterization of G-Quadruplexes in the three-dimensional chromatin structure

DNA molecules are highly compacted in the eukaryotic nucleus where distal regulatory elements reach their targets through three-dimensional chromosomal interactions. G-quadruplexes, stable four-stranded non-canonical DNA structures, can change local chromatin organization through the exclusion of nucleosomes. However, the relationship between G-quadruplexes and higher-order genome organization remains unknown. Here, we found that G-quadruplexes are significantly enriched at boundaries of topological associated domains (TADs). Architectural protein occupancy, which plays critical roles in the formation of TADs, was highly correlated with the content of G-quadruplexes at TAD boundaries. Moreover, adjacent boundaries containing G-quadruplexes frequently interacted with each other because of the high enrichment of architectural protein binding sites. Similar to CCCTC-binding factor (CTCF) binding sites, G-quadruplexes also showed strong insulation ability in the separation of adjacent regions. Additionally, the insulation ability of CTCF binding sites and TAD boundaries was significantly reinforced by G-quadruplexes. Furthermore, G-quadruplex motifs on different strands were associated with the orientation of CTCF binding sites. These findings suggest a potential role for G-quadruplexes in loop extrusion. The enrichment of transcription factor binding sites (TFBSs) around regulatory elements containing G-quadruplexes led to frequent interactions between regulatory elements containing G-quadruplexes. Intriguingly, more than 99% of G-quadruplexes overlapped with TFBSs. The binding sites of CTCF and cohesin proteins were preferentially located surrounding G-quadruplexes. Accordingly, we proposed a new mechanism of long-distance gene regulation in which G-quadruplexes are involved in distal interactions between enhancers and promoters.

Replication timing networks reveal a link between transcription regulatory circuits and replication timing control

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.

Single-cell DNA replication profiling identifies spatiotemporal developmental dynamics of chromosome organization

In mammalian cells, chromosomes are partitioned into megabase-sized topologically associating domains (TADs). TADs can be in either A (active) or B (inactive) subnuclear compartments, which exhibit early and late replication timing (RT), respectively. Here, we show that A/B compartments change coordinately with RT changes genome wide during mouse embryonic stem cell (mESC) differentiation. While A to B compartment changes and early to late RT changes were temporally inseparable, B to A changes clearly preceded late to early RT changes and transcriptional activation. Compartments changed primarily by boundary shifting, altering the compartmentalization of TADs facing the A/B compartment interface, which was conserved during reprogramming and confirmed in individual cells by single-cell Repli-seq. Differentiating mESCs altered single-cell Repli-seq profiles gradually but uniformly, transiently resembling RT profiles of epiblast-derived stem cells (EpiSCs), suggesting that A/B compartments might also change gradually but uniformly toward a primed pluripotent state. These results provide insights into how megabase-scale chromosome organization changes in individual cells during differentiation.

Rapid Irreversible Transcriptional Reprogramming in Human Stem Cells Accompanied by Discordance between Replication Timing and Chromatin Compartment

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.

Replication dynamics of individual loci in single living cells reveal changes in the degree of replication stochasticity through S phase

Eukaryotic genomes are replicated under the control of a highly sophisticated program during the restricted time period corresponding to S phase. The most widely used replication timing assays, which are performed on populations of millions of cells, suggest that most of the genome is synchronously replicated on homologous chromosomes. We investigated the stochastic nature of this temporal program, by comparing the precise replication times of allelic loci within single vertebrate cells progressing through S phase at six loci replicated from very early to very late. We show that replication timing is strictly controlled for the three loci replicated in the first half of S phase. Out of the three loci replicated in the second part of S phase, two present a significantly more stochastic pattern. Surprisingly, we find that the locus replicated at the very end of S phase, presents stochasticity similar to those replicated in early S phase. We suggest that the richness of loci in efficient origins of replication, which decreases from early- to late-replicating regions, and the strength of interaction with the nuclear lamina may underlie the variation of timing control during S phase.

DNA Replication Timing Enters the Single-Cell Era

In mammalian cells, DNA replication timing is controlled at the level of megabase (Mb)-sized chromosomal domains and correlates well with transcription, chromatin structure, and three-dimensional (3D) genome organization. Because of these properties, DNA replication timing is an excellent entry point to explore genome regulation at various levels and a variety of studies have been carried out over the years. However, DNA replication timing studies traditionally required at least tens of thousands of cells, and it was unclear whether the replication domains detected by cell population analyses were preserved at the single-cell level. Recently, single-cell DNA replication profiling methods became available, which revealed that the Mb-sized replication domains detected by cell population analyses were actually well preserved in individual cells. In this article, we provide a brief overview of our current knowledge on DNA replication timing regulation in mammals based on cell population studies, outline the findings from single-cell DNA replication profiling, and discuss future directions and challenges.

Genome-wide stability of the DNA replication program in single mammalian cells

Here, we report a single-cell DNA replication sequencing method, scRepli-seq, a genome-wide methodology that measures copy number differences between replicated and unreplicated DNA. Using scRepli-seq, we demonstrate that replication-domain organization is conserved among individual mouse embryonic stem cells (mESCs). Differentiated mESCs exhibited distinct profiles, which were also conserved among cells. Haplotype-resolved scRepli-seq revealed similar replication profiles of homologous autosomes, while the inactive X chromosome was clearly replicated later than its active counterpart. However, a small degree of cell-to-cell replication-timing heterogeneity was present, which was smallest at the beginning and the end of S phase. In addition, developmentally regulated domains were found to deviate from others and showed a higher degree of heterogeneity, thus suggesting a link to developmental plasticity. Moreover, allelic expression imbalance was found to strongly associate with replication-timing asynchrony. Our results form a foundation for single-cell-level understanding of DNA replication regulation and provide insights into three-dimensional genome organization.

Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication

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 the necessity 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.

Linking the organization of DNA replication with genome maintenance

The spatial and temporal organization of genome duplication, also referred to as the replication program, is defined by the distribution and the activities of the sites of replication initiation across the genome. Alterations to the replication profile are associated with cell fate changes during development and in pathologies, but the importance of undergoing S phase with distinct and specific programs remains largely unexplored. We have recently addressed this question, focusing on the interplay between the replication program and genome maintenance. In particular, we demonstrated that when cells encounter challenges to DNA synthesis, the organization of DNA replication drives the response to replication stress that is mediated by the ATR/Rad3 checkpoint pathway, thus shaping the pattern of genome instability along the chromosomes. In this review, we present the major findings of our study and discuss how they may bring new perspectives to our understanding of the biological importance of the replication program.

Parallels between Mammalian Mechanisms of Monoallelic Gene Expression

Different types of monoallelic gene expression are present in mammals, some of which are highly flexible, whereas others are more rigid. These include allelic exclusion at antigen receptor loci, the expression of olfactory receptor genes, genomic imprinting, X-chromosome inactivation, and random monoallelic expression (MAE). Although these processes play diverse biological roles, and arose through different selective pressures, the underlying epigenetic mechanisms show striking resemblances. Regulatory transcriptional events are important in all systems, particularly in the specification of MAE. Combined with comparative studies between species, this suggests that the different MAE systems found in mammals may have evolved from analogous ancestral processes.

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.

Regulation of the program of DNA replication by CDK: new findings and perspectives

Progression through the cell cycle is driven by the activities of the cyclin-dependent kinase (CDK) family of enzymes, which establish an ordered passage through the cell cycle phases. CDK activity is crucial for the cellular transitions from G1 to S and G2 to M, which are highly controlled to promote the faithful duplication of the genetic material and the transmission of the genome into daughter cells, respectively. While oscillations in CDK activity are essential for cell division, how its specific dynamics may shape cellular processes remains an open question. Recently, we have investigated the potential role of CDK in establishing the profile of replication initiation along the chromosomes, also referred to as the replication program. Our results demonstrated that the timing and level of CDK activity at G1/S provide two critical and independent inputs that modulate the pattern of origin usage. In this review, we will present the conclusions of our study and discuss the implications of our findings for cellular function and physiology.

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.