DNA replication timing alterations identify common markers between distinct progeroid diseases

RIF1 regulates early replication timing in murine B cells

The mammalian DNA replication timing (RT) program is crucial for the proper functioning and integrity of the genome. The best-known mechanism for controlling RT is the suppression of late origins of replication in heterochromatin by RIF1. Here, we report that in antigen-activated, hypermutating murine B lymphocytes, RIF1 binds predominantly to early-replicating active chromatin and promotes early replication, but plays a minor role in regulating replication origin activity, gene expression and genome organization in B cells. Furthermore, we find that RIF1 functions in a complementary and non-epistatic manner with minichromosome maintenance (MCM) proteins to establish early RT signatures genome-wide and, specifically, to ensure the early replication of highly transcribed genes. These findings reveal additional layers of regulation within the B cell RT program, driven by the coordinated activity of RIF1 and MCM proteins.

Repli-seq Sample Preparation using Cell Sorting with Cell-Permeant Dyes

Replication timing is significantly correlated with gene expression and chromatin organization, changes dynamically during cell differentiation, and is altered in diseased states. Genome-wide analysis of replication timing is performed in actively replicating cells by Repli-seq. Current methods for Repli-seq require cells to be fixed in large numbers. This is a barrier for sample types that are sensitive to fixation or are in very limited numbers. In this article, we outline optimized methods to process live cells and intact nuclei for Repli-seq. Our protocol enables the processing of a smaller number of cells per sample and reduces processing time and sample loss while obtaining high-quality data. Further, for samples that tend to form clumps and are difficult to dissociate into a single-cell suspension, we also outline methods for isolation, staining, and processing of nuclei for Repli-seq. The Repli-seq data obtained from live cells and intact nuclei are comparable to those obtained from the standard protocols. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Live cell isolation and staining Alternate Protocol: Nuclei isolation and staining.

Clearance of defective muscle stem cells by senolytics restores myogenesis in myotonic dystrophy type 1

Muscle stem cells, the engine of muscle repair, are affected in myotonic dystrophy type 1 (DM1); however, the underlying molecular mechanism and the impact on the disease severity are still elusive. Here, we show using patients' samples that muscle stem cells/myoblasts exhibit signs of cellular senescence in vitro and in situ. Single cell RNAseq uncovers a subset of senescent myoblasts expressing high levels of genes related to the senescence-associated secretory phenotype (SASP). We show that the levels of interleukin-6, a prominent SASP cytokine, in the serum of DM1 patients correlate with muscle weakness and functional capacity limitations. Drug screening revealed that the senolytic BCL-XL inhibitor (A1155463) can specifically remove senescent DM1 myoblasts by inducing their apoptosis. Clearance of senescent cells reduced the expression of SASP, which rescued the proliferation and differentiation capacity of DM1 myoblasts in vitro and enhanced their engraftment following transplantation in vivo. Altogether, this study identifies the pathogenic mechanism associated with muscle stem cell defects in DM1 and opens a therapeutic avenue that targets these defective cells to restore myogenesis.

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.

Cell division drives DNA methylation loss in late-replicating domains in primary human cells

DNA methylation undergoes dramatic age-related changes, first described more than four decades ago. Loss of DNA methylation within partially methylated domains (PMDs), late-replicating regions of the genome attached to the nuclear lamina, advances with age in normal tissues, and is further exacerbated in cancer. We present here experimental evidence that this DNA hypomethylation is directly driven by proliferation-associated DNA replication. Within PMDs, loss of DNA methylation at low-density CpGs in A:T-rich immediate context (PMD solo-WCGWs) tracks cumulative population doublings in primary cell culture. Cell cycle deceleration results in a proportional decrease in the rate of DNA hypomethylation. Blocking DNA replication via Mitomycin C treatment halts methylation loss. Loss of methylation continues unabated after TERT immortalization until finally reaching a severely hypomethylated equilibrium. Ambient oxygen culture conditions increases the rate of methylation loss compared to low-oxygen conditions, suggesting that some methylation loss may occur during unscheduled, oxidative damage repair-associated DNA synthesis. Finally, we present and validate a model to estimate the relative cumulative replicative histories of human cells, which we call "RepliTali" (Replication Times Accumulated in Lifetime).

Comprehensive analysis of DNA replication timing across 184 cell lines suggests a role for MCM10 in replication timing regulation

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. Cells from a single patient carrying MCM10 mutations demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in the mutated MCM10 cells were predominantly comprised of replication delays and initiation site gains and losses. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel candidate modulator of DNA replication timing.

Mapping Replication Timing in Single Mammalian Cells

Replication timing (RT) is the temporal order in which genomic DNA is replicated during S phase. Early and late replication correlate with transcriptionally active and inactive chromatin compartments, but mechanistic links between large-scale chromosome structure, transcription, and replication are still enigmatic. A proper RT program is necessary to maintain the global epigenome that defines cell identity, suggesting that RT is critical for epigenome integrity by facilitating the assembly of different types of chromatin at different times during S phase. RT is regulated during development and has been found to be altered in disease. Thus, RT can identify stable epigenetic differences distinguishing cell types, and can be used to help stratify patient outcomes and identify markers of disease. Most methods to profile RT require thousands of S-phase cells. In cases where cells are rare (e.g., early-stage embryos or rare primary cell types) or consist of a heterogeneous mixture of cell states (e.g., differentiation intermediates), or when the interest is in determining the degree of stable epigenetic heterogeneity within a population of cells, single-cell measurements of RT are necessary. We have previously developed single cell Repli-seq, a method to measure replication timing in single cells using DNA copy number quantification. To date, however, single-cell Repli-seq suffers from relatively low throughput and high costs. Here, we describe an improved single-cell Repli-seq protocol that uses degenerate oligonucleotide-primed PCR (DOP-PCR) for uniform whole-genome amplification and uniquely barcoded primers that permit early pooling of single-cell samples into a single library preparation. We also provide a bioinformatics platform for analysis of the data. The improved throughput and decreased costs of this method relative to previously published single-cell Repli-seq protocols should make it considerably more accessible to a broad range of investigators. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Whole Genome Amplification (WGA) of single cells and sequence library construction. Basic Protocol 2: Deriving and displaying single-cell replication timing data from whole genome sequencing.

Collective regulation of chromatin modifications predicts replication timing during cell cycle

Replication timing (RT) associates with genome architecture, while having a mixed relationship to histone marks. By profiling replication at high resolution and assessing broad histone marks across the cell cycle at the resolution of RT with and without genetic perturbation, we address the causal relationship between histone marks and RT. Four primary chromatin states, including an uncharacterized H3K36me2 state, emerge and define 97% of the mappable genome. RT and local replication patterns (e.g., initiation zones) quantitatively associate with chromatin states, histone mark dynamics, and spatial chromatin structure. Manipulation of broad histone marks and enhancer elements by overexpressing the histone H3 lysine 9/36 tri-demethylase KDM4A impacts RT across 11% of the genome. Broad histone modification changes were strong predictors of the observed RT alterations. Lastly, replication within H3K36me2-enriched neighborhoods is sensitive to KDM4A overexpression and is controlled at a megabase scale. These studies establish a role for collective chromatin mark regulation in modulating RT.

Mammalian DNA Replication Timing

Immediately following the discovery of the structure of DNA and the semi-conservative replication of the parental DNA sequence into two new DNA strands, it became apparent that DNA replication is organized in a temporal and spatial fashion during the S phase of the cell cycle, correlated with the large-scale organization of chromatin in the nucleus. After many decades of limited progress, technological advances in genomics, genome engineering, and imaging have finally positioned the field to tackle mechanisms underpinning the temporal and spatial regulation of DNA replication and the causal relationships between DNA replication and other features of large-scale chromosome structure and function. In this review, we discuss these major recent discoveries as well as expectations for the coming decade.

Progerin impairs 3D genome organization and induces fragile telomeres by limiting the dNTP pools

Chromatin organization within the nuclear volume is essential to regulate many aspects of its function and to safeguard its integrity. A key player in this spatial scattering of chromosomes is the nuclear envelope (NE). The NE tethers large chromatin domains through interaction with the nuclear lamina and other associated proteins. This organization is perturbed in cells from Hutchinson–Gilford progeria syndrome (HGPS), a genetic disorder characterized by premature aging features. Here, we show that HGPS-related lamina defects trigger an altered 3D telomere organization with increased contact sites between telomeres and the nuclear lamina, and an altered telomeric chromatin state. The genome-wide replication timing signature of these cells is perturbed, with a shift to earlier replication for regions that normally replicate late. As a consequence, we detected a higher density of replication forks traveling simultaneously on DNA fibers, which relies on limiting cellular dNTP pools to support processive DNA synthesis. Remarkably, increasing dNTP levels in HGPS cells rescued fragile telomeres, and improved the replicative capacity of the cells. Our work highlights a functional connection between NE dysfunction and telomere homeostasis in the context of premature aging.

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

Artificial Intelligence in Epigenetic Studies: Shedding Light on Rare Diseases

More than 7,000 rare diseases (RDs) exist worldwide, affecting approximately 350 million people, out of which only 5% have treatment. The development of novel genome sequencing techniques has accelerated the discovery and diagnosis in RDs. However, most patients remain undiagnosed. Epigenetics has emerged as a promise for diagnosis and therapies in common disorders (e.g., cancer) with several epimarkers and epidrugs already approved and used in clinical practice. Hence, it may also become an opportunity to uncover new disease mechanisms and therapeutic targets in RDs. In this “big data” age, the amount of information generated, collected, and managed in (bio)medicine is increasing, leading to the need for its rapid and efficient collection, analysis, and characterization. Artificial intelligence (AI), particularly deep learning, is already being successfully applied to analyze genomic information in basic research, diagnosis, and drug discovery and is gaining momentum in the epigenetic field. The application of deep learning to epigenomic studies in RDs could significantly boost discovery and therapy development. This review aims to collect and summarize the application of AI tools in the epigenomic field of RDs. The lower number of studies found, specific for RDs, indicate that this is a field open to expansion, following the results obtained for other more common disorders.

A Link between Replicative Stress, Lamin Proteins, and Inflammation

Double-stranded breaks (DSB), the most toxic DNA lesions, are either a consequence of cellular metabolism, programmed as in during V(D)J recombination, or induced by anti-tumoral therapies or accidental genotoxic exposure. One origin of DSB sources is replicative stress, a major source of genome instability, especially when the integrity of the replication forks is not properly guaranteed. To complete stalled replication, restarting the fork requires complex molecular mechanisms, such as protection, remodeling, and processing. Recently, a link has been made between DNA damage accumulation and inflammation. Indeed, defects in DNA repair or in replication can lead to the release of DNA fragments in the cytosol. The recognition of this self-DNA by DNA sensors leads to the production of inflammatory factors. This beneficial response activating an innate immune response and destruction of cells bearing DNA damage may be considered as a novel part of DNA damage response. However, upon accumulation of DNA damage, a chronic inflammatory cellular microenvironment may lead to inflammatory pathologies, aging, and progression of tumor cells. Progress in understanding the molecular mechanisms of DNA damage repair, replication stress, and cytosolic DNA production would allow to propose new therapeutical strategies against cancer or inflammatory diseases associated with aging. In this review, we describe the mechanisms involved in DSB repair, the replicative stress management, and its consequences. We also focus on new emerging links between key components of the nuclear envelope, the lamins, and DNA repair, management of replicative stress, and inflammation.

Replication timing alterations in leukemia affect clinically relevant chromosome domains

Human B-cell 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 that BCP-ALL cells display unique and clonally heritable, stable DNA replication timing (RT) programs (ie, programs describing the variable order of replication and subnuclear 3D architecture 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 with varied genetic subtypes and outcomes. The results identify BCP-ALL subtype-specific features that resemble specific stages of B-cell differentiation and features that seem to be associated with relapse. These results suggest that the genesis of BCP-ALL involves alterations in RT that reflect biologically significant and potentially clinically relevant leukemia-specific epigenetic changes.

Chromatin and Nuclear Architecture: Shaping DNA Replication in 3D

In eukaryotes, DNA replication progresses through a finely orchestrated temporal and spatial program. The 3D genome structure and nuclear architecture have recently emerged as fundamental determinants of the replication program. Factors with established roles in replication have been recognized as genome organization regulators. Exploiting paradigms from yeasts and mammals, we discuss how DNA replication is regulated in time and space through DNA-associated trans-acting factors, diffusible limiting replication initiation factors, higher-order chromatin folding, dynamic origin localization, and specific nuclear microenvironments. We present an integrated model for the regulation of DNA replication in 3D and highlight the importance of accurate spatio-temporal regulation of DNA replication in physiology and disease.

Chromosome Conformation Capture and Beyond: Toward an Integrative View of Chromosome Structure and Function

Rapidly developing technologies have recently fueled an exciting era of discovery in the field of chromosome structure and nuclear organization. In addition to chromosome conformation capture (3C) methods, new alternative techniques have emerged to study genome architecture and biological processes in the nucleus, often in single or living cells. This sets an unprecedented stage for exploring the mechanisms that link chromosome structure and biological function. Here we review popular as well as emerging approaches to study chromosome organization, focusing on the contribution of complementary methodologies to our understanding of structures revealed by 3C methods and their biological implications, and discuss the next technical and conceptual frontiers.

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.

RT States: systematic annotation of the human genome using cell type-specific replication timing programs

MOTIVATION

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.

RESULTS

In this 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.

AVAILABILITY AND IMPLEMENTATION

R scripts for inferring HMM models and Perl scripts for further analysis are available https://github.com/PouletAxel/script_HMM_Replication_timing.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

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.

Epigenetic Regulation of Skin Cells in Natural Aging and Premature Aging Diseases

Skin undergoes continuous renewal throughout an individual’s lifetime relying on stem cell functionality. However, a decline of the skin regenerative potential occurs with age. The accumulation of senescent cells over time probably reduces tissue regeneration and contributes to skin aging. Keratinocytes and dermal fibroblasts undergo senescence in response to several intrinsic or extrinsic stresses, including telomere shortening, overproduction of reactive oxygen species, diet, and sunlight exposure. Epigenetic mechanisms directly regulate skin homeostasis and regeneration, but they also mark cell senescence and the natural and pathological aging processes. Progeroid syndromes represent a group of clinical and genetically heterogeneous pathologies characterized by the accelerated aging of various tissues and organs, including skin. Skin cells from progeroid patients display molecular hallmarks that mimic those associated with naturally occurring aging. Thus, investigations on progeroid syndromes strongly contribute to disclose the causal mechanisms that underlie the aging process. In the present review, we discuss the role of epigenetic pathways in skin cell regulation during physiologic and premature aging.

Recent advances in understanding DNA replication: cell type-specific adaptation of the DNA replication program

DNA replication is an essential process occurring prior to cell division. Cell division coupled to proliferation ensures the growth and renewal of a large variety of specialized cell types generated during embryonic development. Changes in the DNA replication program occur during development. Embryonic undifferentiated cells show a high replication rate and fast proliferation, whereas more differentiated cells are characterized by reduced DNA synthesis and a low proliferation rate. Hence, the DNA replication program must adapt to the specific features of cells committed to different fates. Recent findings on DNA synthesis regulation in different cell types open new perspectives for developing efficient and more adapted therapies to treat various diseases such as genetic diseases and cancer. This review will put the emphasis on recent progress made in this field.

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

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

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

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