Eukaryotic chromosomes replicate in a temporal order known as the replication-timing program. In mammals, replication timing is cell-type-specific with at least half the genome switching replication timing during development, primarily in units of 400-800 kilobases ('replication domains'), whose positions are preserved in different cell types, conserved between species, and appear to confine long-range effects of chromosome rearrangements. Early and late replication correlate, respectively, with open and closed three-dimensional chromatin compartments identified by high-resolution chromosome conformation capture (Hi-C), and, to a lesser extent, late replication correlates with lamina-associated domains (LADs). Recent Hi-C mapping has unveiled substructure within chromatin compartments called topologically associating domains (TADs) that are largely conserved in their positions between cell types and are similar in size to replication domains. However, TADs can be further sub-stratified into smaller domains, challenging the significance of structures at any particular scale. Moreover, attempts to reconcile TADs and LADs to replication-timing data have not revealed a common, underlying domain structure. Here we localize boundaries of replication domains to the early-replicating border of replication-timing transitions and map their positions in 18 human and 13 mouse cell types. We demonstrate that, collectively, replication domain boundaries share a near one-to-one correlation with TAD boundaries, whereas within a cell type, adjacent TADs that replicate at similar times obscure replication domain boundaries, largely accounting for the previously reported lack of alignment. Moreover, cell-type-specific replication timing of TADs partitions the genome into two large-scale sub-nuclear compartments revealing that replication-timing transitions are indistinguishable from late-replicating regions in chromatin composition and lamina association and accounting for the reduced correlation of replication timing to LADs and heterochromatin. Our results reconcile cell-type-specific sub-nuclear compartmentalization and replication timing with developmentally stable structural domains and offer a unified model for large-scale chromosome structure and function.

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.

Fluorescence recovery after photobleaching (FRAP) is a widely used imaging technique for measuring protein dynamics in live cells that has provided many important biological insights. Although FRAP presumes that the conversion of a fluorophore from a bright to a dark state is irreversible, GFP as well as other genetically encoded fluorescent proteins now in common use can also exhibit a reversible conversion known as photoswitching. Various studies have shown how photoswitching can cause at least four different artifacts in FRAP, leading to false conclusions about various biological phenomena, including the erroneous identification of anomalous diffusion or the overestimation of the freely diffusible fraction of a cellular protein. Unfortunately, identifying and then correcting these artifacts is difficult. Here we report a new characteristic of an organic fluorophore tetramethylrhodamine bound to the HaloTag protein (TMR-HaloTag), which like GFP can be genetically encoded, but which directly and simply overcomes the artifacts caused by photoswitching in FRAP. We show that TMR exhibits virtually no photoswitching in live cells under typical imaging conditions for FRAP. We also demonstrate that TMR eliminates all of the four reported photoswitching artifacts in FRAP. Finally, we apply this photoswitching-free FRAP with TMR to show that the chromatin decondensation following UV irradiation does not involve loss of nucleosomes from the damaged DNA. In sum, we demonstrate that the TMR Halo label provides a genetically encoded fluorescent tag very well suited for accurate FRAP experiments.

We have examined the distribution of early replicating origins on stretched DNA fibers when nuclei from CHO cells synchronized at different times during G1 phase initiate DNA replication in Xenopus egg extracts. Origins were differentially labeled in vivo versus in vitro to allow a comparison of their relative positions and spacing. With nuclei isolated in the first hour of G1 phase, in vitro origins were distributed throughout a larger number of DNA fibers and did not coincide with in vivo origins. With nuclei isolated 1 h later, a similar total number of in vitro origins were clustered within a smaller number of DNA fibers but still did not coincide with in vivo origins. However, with nuclei isolated later in G1 phase, the positions of many in vitro origins coincided with in vivo origin sites without further change in origin number or density. These results highlight two distinct G1 steps that establish a spatial and temporal program for replication.

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.

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.

Chromatin profiling in single cells has been extremely challenging and almost exclusively limited to histone proteins. In cases where single-cell methods have shown promise, many require highly specialized equipment or cell type-specific protocols and are relatively low throughput. Here, we combine the advantages of tagmentation, linear amplification, and combinatorial indexing to produce a high-throughput single-cell DNA binding site mapping method that is simple, inexpensive, and capable of multiplexing several independent samples per experiment. Targeted insertion of promoters sequencing (TIP-seq) uses Tn5 fused to proteinA to insert a T7 RNA polymerase promoter adjacent to a chromatin protein of interest. Linear amplification of flanking DNA with T7 polymerase before sequencing library preparation provides ∼10-fold higher unique reads per single cell compared with other methods. We applied TIP-seq to map histone modifications, RNA polymerase II (RNAPII), and transcription factor CTCF binding sites in single human and mouse cells.

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.

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.

The laboratory mouse shares the majority of its protein-coding genes with humans, making it the premier model organism in biomedical research, yet the two mammals differ in significant ways. To gain greater insights into both shared and species-specific transcriptional and cellular regulatory programs in the mouse, the Mouse ENCODE Consortium has mapped transcription, DNase I hypersensitivity, transcription factor binding, chromatin modifications and replication domains throughout the mouse genome in diverse cell and tissue types. By comparing with the human genome, we not only confirm substantial conservation in the newly annotated potential functional sequences, but also find a large degree of divergence of sequences involved in transcriptional regulation, chromatin state and higher order chromatin organization. Our results illuminate the wide range of evolutionary forces acting on genes and their regulatory regions, and provide a general resource for research into mammalian biology and mechanisms of human diseases.

The genomic neighborhood of a gene influences its activity, a behavior that is attributable in part to domain-scale regulation. Previous genomic studies have identified many types of regulatory domains. However, due to the difficulty of integrating genomics data sets, the relationships among these domain types are poorly understood. Semi-automated genome annotation (SAGA) algorithms facilitate human interpretation of heterogeneous collections of genomics data by simultaneously partitioning the human genome and assigning labels to the resulting genomic segments. However, existing SAGA methods cannot integrate inherently pairwise chromatin conformation data. We developed a new computational method, called graph-based regularization (GBR), for expressing a pairwise prior that encourages certain pairs of genomic loci to receive the same label in a genome annotation. We used GBR to exploit chromatin conformation information during genome annotation by encouraging positions that are close in 3D to occupy the same type of domain. Using this approach, we produced a model of chromatin domains in eight human cell types, thereby revealing the relationships among known domain types. Through this model, we identified clusters of tightly regulated genes expressed in only a small number of cell types, which we term "specific expression domains." We found that domain boundaries marked by promoters and CTCF motifs are consistent between cell types even when domain activity changes. Finally, we showed that GBR can be used to transfer information from well-studied cell types to less well-characterized cell types during genome annotation, making it possible to produce high-quality annotations of the hundreds of cell types with limited available data.

How subunit dosage contributes to the assembly and function of multimeric complexes is an important question with implications in understanding biochemical, evolutionary, and disease mechanisms. Toward identifying pathways that are susceptible to decreased gene dosage, we performed a genome-wide screen for haploinsufficient (HI) genes that guard against genome instability in Saccharomyces cerevisiae. This led to the identification of all three genes (SPC97, SPC98, and TUB4) encoding the evolutionarily conserved γ-tubulin small complex (γ-TuSC), which nucleates microtubule assembly. We found that hemizygous γ-TuSC mutants exhibit higher rates of chromosome loss and increases in anaphase spindle length and elongation velocities. Fluorescence microscopy, fluorescence recovery after photobleaching, electron tomography, and model convolution simulation of spc98/+ mutants revealed improper regulation of interpolar (iMT) and kinetochore (kMT) microtubules in anaphase. The underlying cause is likely due to reduced levels of Tub4, as overexpression of TUB4 suppressed the spindle and chromosome segregation defects in spc98/+ mutants. We propose that γ-TuSC is crucial for balanced assembly between iMTs and kMTs for spindle organization and accurate chromosome segregation. Taken together, the results show how gene dosage studies provide critical insights into the assembly and function of multisubunit complexes that may not be revealed by using traditional studies with haploid gene deletion or conditional alleles.

During meiotic prophase in male mammals, the X and Y chromosomes condense to form a macrochromatin body, termed the sex, or XY, body, within which X- and Y-linked genes are transcriptionally repressed. The molecular basis and biological function of both sex body formation and meiotic sex chromosome inactivation (MSCI) are unknown. A phosphorylated form of H2AX, a histone H2A variant implicated in DNA repair, accumulates in the sex body in a manner independent of meiotic recombination-associated double-strand breaks. Here we show that the X and Y chromosomes of histone H2AX-deficient spermatocytes fail to condense to form a sex body, do not initiate MSCI, and exhibit severe defects in meiotic pairing. Moreover, other sex body proteins, including macroH2A1.2 and XMR, do not preferentially localize with the sex chromosomes in the absence of H2AX. Thus, H2AX is required for the chromatin remodeling and associated silencing in male meiosis.

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,000 N 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.

Primary human hepatocytes isolated from patient biopsies represent the most physiologically relevant cell culture model for hepatitis C virus (HCV) infection, but these primary cells are not readily accessible, display individual variability, and are largely refractory to genetic manipulation. Hepatocyte-like cells differentiated from pluripotent stem cells provide an attractive alternative as they not only overcome these shortcomings but can also provide an unlimited source of noncancer cells for both research and cell therapy. Despite its promise, the permissiveness to HCV infection of differentiated human hepatocyte-like cells (DHHs) has not been explored. Here we report a novel infection model based on DHHs derived from human embryonic (hESCs) and induced pluripotent stem cells (iPSCs). DHHs generated in chemically defined media under feeder-free conditions were subjected to infection by both HCV derived in cell culture (HCVcc) and patient-derived virus (HCVser). Pluripotent stem cells and definitive endoderm were not permissive for HCV infection whereas hepatic progenitor cells were persistently infected and secreted infectious particles into culture medium. Permissiveness to infection was correlated with induction of the liver-specific microRNA-122 and modulation of cellular factors that affect HCV replication. RNA interference directed toward essential cellular cofactors in stem cells resulted in HCV-resistant hepatocyte-like cells after differentiation. The ability to infect cultured cells directly with HCV patient serum, to study defined stages of viral permissiveness, and to produce genetically modified cells with desired phenotypes all have broad significance for host-pathogen interactions and cell therapy.

Dynamic access to genetic information is central to organismal development and environmental response. Consequently, genomic processes must be regulated by mechanisms that alter genome function relatively rapidly. Conventional chromatin immunoprecipitation (ChIP) experiments measure transcription factor occupancy, but give no indication of kinetics and are poor predictors of transcription factor function at a given locus. To measure transcription-factor-binding dynamics across the genome, we performed competition ChIP (refs 6, 7) with a sequence-specific Saccharomyces cerevisiae transcription factor, Rap1 (ref. 8). Rap1-binding dynamics and Rap1 occupancy were only weakly correlated (R(2) = 0.14), but binding dynamics were more strongly linked to function than occupancy. Long Rap1 residence was coupled to transcriptional activation, whereas fast binding turnover, which we refer to as 'treadmilling', was linked to low transcriptional output. Thus, DNA-binding events that seem identical by conventional ChIP may have different underlying modes of interaction that lead to opposing functional outcomes. We propose that transcription factor binding turnover is a major point of regulation in determining the functional consequences of transcription factor binding, and is mediated mainly by control of competition between transcription factors and nucleosomes. Our model predicts a clutch-like mechanism that rapidly engages a treadmilling transcription factor into a stable binding state, or vice versa, to modulate transcription factor function.

According to the transcription factory model, localized transcription sites composed of immobilized polymerase molecules transcribe chromatin by reeling it through the transcription site and extruding it to form a surrounding domain of recently transcribed decondensed chromatin. Although transcription sites have been identified in various cells, surrounding domains of recently transcribed decondensed chromatin have not. We report evidence that transcription sites associated with a tandem gene array in mouse cells are indeed surrounded by or adjacent to a domain of decondensed chromatin composed of sequences from the gene array. Formation of this decondensed domain requires transcription and topoisomerase IIalpha activity. The decondensed domain is enriched for the trimethyl H3K36 mark that is associated with recently transcribed chromatin in yeast and several mammalian systems. Consistent with this, chromatin immunoprecipitation demonstrates a comparable enrichment of this mark in transcribed sequences at the tandem gene array. These results provide new support for the pol II factory model, in which an immobilized polymerase molecule extrudes decondensed, transcribed sequences into its surroundings.

Three-dimensional genome organisation and replication timing are known to be correlated, however, it remains unknown whether nuclear architecture overall plays an instructive role in the replication-timing programme and, if so, how. Here we demonstrate that RIF1 is a molecular hub that co-regulates both processes. Both nuclear organisation and replication timing depend upon the interaction between RIF1 and PP1. However, whereas nuclear architecture requires the full complement of RIF1 and its interaction with PP1, replication timing is not sensitive to RIF1 dosage. The role of RIF1 in replication timing also extends beyond its interaction with PP1. Availing of this separation-of-function approach, we have therefore identified in RIF1 dual function the molecular bases of the co-dependency of the replication-timing programme and nuclear architecture.

Chromatin boundary elements contribute to the partitioning of mammalian genomes into topological domains to regulate gene expression. Certain boundary elements are adopted as DNA insulators for safe and stable transgene expression in mammalian cells. These elements, however, are ill-defined and less characterized in the non-coding genome, partially due to the lack of a platform to readily evaluate boundary-associated activities of putative DNA sequences. Here we report SHIELD (Site-specific Heterochromatin Insertion of Elements at Lamina-associated Domains), a platform tailored for the high-throughput screening of barrier-type DNA elements in human cells. SHIELD takes advantage of the high specificity of serine integrase at heterochromatin, and exploits the natural heterochromatin spreading inside lamina-associated domains (LADs) for the discovery of potent barrier elements. We adopt SHIELD to evaluate the barrier activity of 1000 DNA elements in a high-throughput manner and identify 8 candidates with barrier activities comparable to the core region of cHS4 element in human HCT116 cells. We anticipate SHIELD could facilitate the discovery of novel barrier DNA elements from the non-coding genome in human cells.

Background: It has been known for many years that in metazoan cells, replication origins are organised into clusters where origins within each cluster fire near-synchronously. Despite clusters being a fundamental organising principle of metazoan DNA replication, the genomic location of origin clusters has not been documented. Methods: We synchronised human U2OS by thymidine block and release followed by L-mimosine block and release to create a population of cells progressing into S phase with a high degree of synchrony. At different times after release into S phase, cells were pulsed with EdU; the EdU-labelled DNA was then pulled down, sequenced and mapped onto the human genome. Results: The early replicating DNA showed features at a range of scales. Wavelet analysis showed that the major feature of the early replicating DNA was at a size of 500 kb, consistent with clusters of replication origins. Over the first two hours of S phase, these Replicon Cluster Domains broadened in width, consistent with their being enlarged by the progression of replication forks at their outer boundaries. The total replication signal associated with each Replicon Cluster Domain varied considerably, and this variation was reproducible and conserved over time. We provide evidence that this variability in replication signal was at least in part caused by Replicon Cluster Domains being activated at different times in different cells in the population. We also provide evidence that adjacent clusters had a statistical preference for being activated in sequence across a group, consistent with the 'domino' model of replication focus activation order observed by microscopy. Conclusions: We show that early replicating DNA is organised into Replicon Cluster Domains that behave as expected of replicon clusters observed by DNA fibre analysis. The coordinated activation of different Replicon Cluster Domains can generate the replication timing programme by which the genome is duplicated.