DNA replication and transcription programs respond to the same chromatin cues

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Multi-Scale Imaging of the Dynamic Organization of Chromatin

Chromatin is now regarded as a heterogeneous and dynamic structure occupying a non-random position within the cell nucleus, where it plays a key role in regulating various functions of the genome. This current view of chromatin has emerged thanks to high spatiotemporal resolution imaging, among other new technologies developed in the last decade. In addition to challenging early assumptions of chromatin being regular and static, high spatiotemporal resolution imaging made it possible to visualize and characterize different chromatin structures such as clutches, domains and compartments. More specifically, super-resolution microscopy facilitates the study of different cellular processes at a nucleosome scale, providing a multi-scale view of chromatin behavior within the nucleus in different environments. In this review, we describe recent imaging techniques to study the dynamic organization of chromatin at high spatiotemporal resolution. We also discuss recent findings, elucidated by these techniques, on the chromatin landscape during different cellular processes, with an emphasis on the DNA damage response.

ATXN3 controls DNA replication and transcription by regulating chromatin structure

The deubiquitinating enzyme Ataxin-3 (ATXN3) contains a polyglutamine (PolyQ) region, the expansion of which causes spinocerebellar ataxia type-3 (SCA3). ATXN3 has multiple functions, such as regulating transcription or controlling genomic stability after DNA damage. Here we report the role of ATXN3 in chromatin organization during unperturbed conditions, in a catalytic-independent manner. The lack of ATXN3 leads to abnormalities in nuclear and nucleolar morphology, alters DNA replication timing and increases transcription. Additionally, indicators of more open chromatin, such as increased mobility of histone H1, changes in epigenetic marks and higher sensitivity to micrococcal nuclease digestion were detected in the absence of ATXN3. Interestingly, the effects observed in cells lacking ATXN3 are epistatic to the inhibition or lack of the histone deacetylase 3 (HDAC3), an interaction partner of ATXN3. The absence of ATXN3 decreases the recruitment of endogenous HDAC3 to the chromatin, as well as the HDAC3 nuclear/cytoplasm ratio after HDAC3 overexpression, suggesting that ATXN3 controls the subcellular localization of HDAC3. Importantly, the overexpression of a PolyQ-expanded version of ATXN3 behaves as a null mutant, altering DNA replication parameters, epigenetic marks and the subcellular distribution of HDAC3, giving new insights into the molecular basis of the disease.

The MOZ-BRPF1 acetyltransferase complex in epigenetic crosstalk linked to gene regulation, development, and human diseases

Acetylation of lysine residues on histone tails is an important post-translational modification (PTM) that regulates chromatin dynamics to allow gene transcription as well as DNA replication and repair. Histone acetyltransferases (HATs) are often found in large multi-subunit complexes and can also modify specific lysine residues in non-histone substrates. Interestingly, the presence of various histone PTM recognizing domains (reader domains) in these complexes ensures their specific localization, enabling the epigenetic crosstalk and context-specific activity. In this review, we will cover the biochemical and functional properties of the MOZ-BRPF1 acetyltransferase complex, underlining its role in normal biological processes as well as in disease progression. We will discuss how epigenetic reader domains within the MOZ-BRPF1 complex affect its chromatin localization and the histone acetyltransferase specificity of the complex. We will also summarize how MOZ-BRPF1 is linked to development via controlling cell stemness and how mutations or changes in expression levels of MOZ/BRPF1 can lead to developmental disorders or cancer. As a last touch, we will review the latest drug candidates for these two proteins and discuss the therapeutic possibilities.

Distinct signaling signatures drive compensatory proliferation via S-phase acceleration

Regeneration relies on cell proliferation to restore damaged tissues. Multiple signaling pathways activated by local or paracrine cues have been identified to promote regenerative proliferation. How different types of tissue damage may activate distinct signaling pathways and how these differences converge on regenerative proliferation is less well defined. To better understand how tissue damage and proliferative signals are integrated during regeneration, we investigate models of compensatory proliferation in Drosophila imaginal discs. We find that compensatory proliferation is associated with a unique cell cycle profile, which is characterized by short G1 and G2 phases and, surprisingly, by acceleration of the S-phase. S-phase acceleration can be induced by two distinct signaling signatures, aligning with inflammatory and non-inflammatory tissue damage. Specifically, non-autonomous activation of JAK/STAT and Myc in response to inflammatory damage, or local activation of Ras/ERK and Hippo/Yki in response to elevated cell death, promote accelerated nucleotide incorporation during S-phase. This previously unappreciated convergence of different damaging insults on the same regenerative cell cycle program reconciles previous conflicting observations on proliferative signaling in different tissue regeneration and tumor models.

A Drosophila insulator interacting protein suppresses enhancer-blocking function and modulates replication timing

Insulators play important roles in genome structure and function in eukaryotes. Interactions between a DNA binding insulator protein and its interacting partner proteins define the properties of each insulator site. The different roles of insulator protein partners in the Drosophila genome and how they confer functional specificity remain poorly understood. The Suppressor of Hairy wing [Su(Hw)] insulator is targeted to the nuclear lamina, preferentially localizes at euchromatin/heterochromatin boundaries, and is associated with the gypsy retrotransposon. Insulator activity relies on the ability of the Su(Hw) protein to bind the DNA at specific sites and interact with Mod(mdg4)67.2 and CP190 partner proteins. HP1 and insulator partner protein 1 (HIPP1) is a partner of Su(Hw), but how HIPP1 contributes to the function of Su(Hw) insulator complexes is unclear. Here, we demonstrate that HIPP1 colocalizes with the Su(Hw) insulator complex in polytene chromatin and in stress-induced insulator bodies. We find that the overexpression of either HIPP1 or Su(Hw) or mutation of the HIPP1 crotonase-like domain (CLD) causes defects in cell proliferation by limiting the progression of DNA replication. We also show that HIPP1 overexpression suppresses the Su(Hw) insulator enhancer-blocking function, while mutation of the HIPP1 CLD does not affect Su(Hw) enhancer blocking. These findings demonstrate a functional relationship between HIPP1 and the Su(Hw) insulator complex and suggest that the CLD, while not involved in enhancer blocking, influences cell cycle progression.

The Fkh1 Forkhead associated domain promotes ORC binding to a subset of DNA replication origins in budding yeast

The pioneer event in eukaryotic DNA replication is binding of chromosomal DNA by the origin recognitioncomplex (ORC). The ORC-DNA complex directs the formation of origins, the specific chromosomal regions where DNA synthesis initiates. In all eukaryotes, incompletely understood features of chromatin promote ORC-DNA binding. Here, we uncover a role for the Fkh1 (Forkhead homolog) protein and its forkhead associated (FHA) domain in promoting ORC-origin binding and origin activity at a subset of origins in Saccharomyces cerevisiae. Several of the FHA-dependent origins examined required a distinct Fkh1 binding site located 5′ of and proximal to their ORC sites (5′-FKH-T site). Genetic and molecular experiments provided evidence that the Fkh1-FHA domain promoted origin activity directly through Fkh1 binding to this 5′ FKH-T site. Nucleotide substitutions within two relevant origins that enhanced their ORC-DNA affinity bypassed the requirement for their 5′ FKH-T sites and for the Fkh1-FHA domain. Significantly, assessment of ORC-origin binding by ChIPSeq provided evidence that this mechanism was relevant at ∼25% of yeast origins. Thus, the FHA domain of the conserved cell-cycle transcription factor Fkh1 enhanced origin selection in yeast at the level of ORC-origin binding.

Rtt109 slows replication speed by histone N-terminal acetylation

The wrapping of DNA around histone octamers challenges processes that use DNA as their template. In vitro, DNA replication through chromatin depends on histone modifiers, raising the possibility that cells modify histones to optimize fork progression. Rtt109 is an acetyl transferase that acetylates histone H3 before its DNA incorporation on the K56 and N-terminal residues. We previously reported that, in budding yeast, a wave of histone H3 K9 acetylation progresses ∼3–5 kb ahead of the replication fork. Whether this wave contributes to replication dynamics remained unknown. Here, we show that the replication fork velocity increases following deletion of RTT109, the gene encoding the enzyme required for the prereplication H3 acetylation wave. By using histone H3 mutants, we find that Rtt109-dependent N-terminal acetylation regulates fork velocity, whereas K56 acetylation contributes to replication dynamics only when N-terminal acetylation is compromised. We propose that acetylation of newly synthesized histones slows replication by promoting replacement of nucleosomes evicted by the incoming fork, thereby protecting genome integrity.

Comparing DNA replication programs reveals large timing shifts at centromeres of endocycling cells in maize roots

Plant cells undergo two types of cell cycles–the mitotic cycle in which DNA replication is coupled to mitosis, and the endocycle in which DNA replication occurs in the absence of cell division. To investigate DNA replication programs in these two types of cell cycles, we pulse labeled intact root tips of maize (Zea mays) with 5-ethynyl-2’-deoxyuridine (EdU) and used flow sorting of nuclei to examine DNA replication timing (RT) during the transition from a mitotic cycle to an endocycle. Comparison of the sequence-based RT profiles showed that most regions of the maize genome replicate at the same time during S phase in mitotic and endocycling cells, despite the need to replicate twice as much DNA in the endocycle and the fact that endocycling is typically associated with cell differentiation. However, regions collectively corresponding to 2% of the genome displayed significant changes in timing between the two types of cell cycles. The majority of these regions are small with a median size of 135 kb, shift to a later RT in the endocycle, and are enriched for genes expressed in the root tip. We found larger regions that shifted RT in centromeres of seven of the ten maize chromosomes. These regions covered the majority of the previously defined functional centromere, which ranged between 1 and 2 Mb in size in the reference genome. They replicate mainly during mid S phase in mitotic cells but primarily in late S phase of the endocycle. In contrast, the immediately adjacent pericentromere sequences are primarily late replicating in both cell cycles. Analysis of CENH3 enrichment levels in 8C vs 2C nuclei suggested that there is only a partial replacement of CENH3 nucleosomes after endocycle replication is complete. The shift to later replication of centromeres and possible reduction in CENH3 enrichment after endocycle replication is consistent with a hypothesis that centromeres are inactivated when their function is no longer needed.

In traditional cell division, or mitosis, a cell’s genetic material is duplicated and then split between two daughter cells. In contrast, in some specialized cell types, the DNA is duplicated a second time without an intervening division step, resulting in cells that carry twice as much DNA. This phenomenon, which is called the endocycle, is common during plant development. At each step, DNA replication follows an ordered program in which highly compacted DNA is unraveled and replicated in sections at different times during the synthesis (S) phase. In plants, it is unclear whether traditional and endocycle programs are the same, especially since endocycling cells are typically in the process of differentiation. Using root tips of maize, we found that in comparison to replication in the mitotic cell cycle, there is a small portion of the genome whose replication in the endocycle is shifted in time, usually to later in S phase. Some of these regions are scattered around the genome and mostly coincide with active genes. However, the most prominent shifts occur in centromeres. The shift to later replication in centromeres is noteworthy because they orchestrate the process of separating duplicated chromosomes into daughter cells, a function that is not needed in the endocycle.

3D genome organization contributes to genome instability at fragile sites

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

Common fragile sites are regions susceptible to replication stress and are prone to chromosomal instability. Here, the authors, by analyzing the contribution of 3D chromatin organization, identify and characterize a fragility signature and precisely map these fragility regions.

Rif1 Functions in a Tissue-Specific Manner To Control Replication Timing Through Its PP1-Binding Motif

Replication initiation in eukaryotic cells occurs asynchronously throughout S phase, yielding early- and late-replicating regions of the genome, a process known as replication timing (RT). RT changes during development to ensure accurate genome duplication and maintain genome stability. To understand the relative contributions that cell lineage, cell cycle, and replication initiation regulators have on RT, we utilized the powerful developmental systems available in Drosophila melanogaster We generated and compared RT profiles from mitotic cells of different tissues and from mitotic and endocycling cells of the same tissue. Our results demonstrate that cell lineage has the largest effect on RT, whereas switching from a mitotic to an endoreplicative cell cycle has little to no effect on RT. Additionally, we demonstrate that the RT differences we observed in all cases are largely independent of transcriptional differences. We also employed a genetic approach in these same cell types to understand the relative contribution the eukaryotic RT control factor, Rif1, has on RT control. Our results demonstrate that Rif1 can function in a tissue-specific manner to control RT. Importantly, the Protein Phosphatase 1 (PP1) binding motif of Rif1 is essential for Rif1 to regulate RT. Together, our data support a model in which the RT program is primarily driven by cell lineage and is further refined by Rif1/PP1 to ultimately generate tissue-specific RT programs.

Replication timing in Drosophila and its peculiarities in polytene chromosomes

Drosophila melanogaster is one of the popular model organisms in DNA replication studies. Since the 1960s, DNA replication of polytene chromosomes has been extensively studied by cytological methods. In the recent two decades, the progress in our understanding of DNA replication was associated with new techniques. Use of fluorescent dyes increased the resolution of cytological methods significantly. High-throughput methods allowed analysis of DNA replication on a genome scale, as well as its correlation with chromatin structure and gene activi ty. Precise mapping of the cytological structures of polytene chromosomes to the genome assembly allowed comparison of replication between polytene chromosomes and chromosomes of diploid cells. New features of replication characteristic for D. melanogaster were described for both diploid and polytene chromosomes. Comparison of genomic replication profiles revealed a significant similarity between Drosophila and other well-studi ed eukaryotic species, such as human. Early replication is often confined to intensely transcribed gene-dense regions characterized by multiple replication initiation sites. Features of DNA replication in Drosophila might be explained by a compact genome. The organization of replication in polytene chromosomes has much in common with the organization of replication in chromosomes in diploid cells. The most important feature of replication in polytene chromosomes is its low rate and the dependence of S-phase duration on many factors: external and internal, local and global. The speed of replication forks in D. melanogaster polytene chromosomes is affected by SUUR and Rif1 proteins. It is not known yet how universal the mechanisms associated with these factors are, but their study is very promising.

Differences in firing efficiency, chromatin, and transcription underlie the developmental plasticity of the Arabidopsis DNA replication origins

Eukaryotic genome replication depends on thousands of DNA replication origins (ORIs). A major challenge is to learn ORI biology in multicellular organisms in the context of growing organs to understand their developmental plasticity. We have identified a set of ORIs of Arabidopsis thaliana and their chromatin landscape at two stages of post-embryonic development. ORIs associate with multiple chromatin signatures including transcription start sites (TSS) but also proximal and distal regulatory regions and heterochromatin, where ORIs colocalize with retrotransposons. In addition, quantitative analysis of ORI activity led us to conclude that strong ORIs have high GC content and clusters of GGN trinucleotides. Development primarily influences ORI firing strength rather than ORI location. ORIs that preferentially fire at early developmental stages colocalize with GC-rich heterochromatin, but at later stages with transcribed genes, perhaps as a consequence of changes in chromatin features associated with developmental processes. Our study provides the set of ORIs active in an organism at the post-embryo stage that should allow us to study ORI biology in response to development, environment, and mutations with a quantitative approach. In a wider scope, the computational strategies developed here can be transferred to other eukaryotic systems.

Cohesin-Mediated Genome Architecture Does Not Define DNA Replication Timing Domains

3D genome organization is strongly predictive of DNA replication timing in mammalian cells. This work tested the extent to which loop-based genome architecture acts as a regulatory unit of replication timing by using an auxin-inducible system for acute cohesin ablation. Cohesin ablation in a population of cells in asynchronous culture was shown not to disrupt patterns of replication timing as assayed by replication sequencing (RepliSeq) or BrdU-focus microscopy. Furthermore, cohesin ablation prior to S phase entry in synchronized cells was similarly shown to not impact replication timing patterns. These results suggest that cohesin-mediated genome architecture is not required for the execution of replication timing patterns in S phase, nor for the establishment of replication timing domains in G1.

H3K9 Promotes Under-Replication of Pericentromeric Heterochromatin in Drosophila Salivary Gland Polytene Chromosomes

Chromatin structure and its organization contributes to the proper regulation and timing of DNA replication. Yet, the precise mechanism by which chromatin contributes to DNA replication remains incompletely understood. This is particularly true for cell types that rely on polyploidization as a developmental strategy for growth and high biosynthetic capacity. During Drosophila larval development, cells of the salivary gland undergo endoreplication, repetitive rounds of DNA synthesis without intervening cell division, resulting in ploidy values of ~1350C. S phase of these endocycles displays a reproducible pattern of early and late replicating regions of the genome resulting from the activity of the same replication initiation factors that are used in diploid cells. However, unlike diploid cells, the latest replicating regions of polyploid salivary gland genomes, composed primarily of pericentric heterochromatic enriched in H3K9 methylation, are not replicated each endocycle, resulting in under-replicated domains with reduced ploidy. Here, we employ a histone gene replacement strategy in Drosophila to demonstrate that mutation of a histone residue important for heterochromatin organization and function (H3K9) but not mutation of a histone residue important for euchromatin function (H4K16), disrupts proper endoreplication in Drosophila salivary gland polyploid genomes thereby leading to DNA copy gain in pericentric heterochromatin. These findings reveal that H3K9 is necessary for normal levels of under-replication of pericentric heterochromatin and suggest that under-replication at pericentric heterochromatin is mediated through H3K9 methylation.

RNA-on-X 1 and 2 in Drosophila melanogaster fulfill separate functions in dosage compensation

In Drosophila melanogaster, the male-specific lethal (MSL) complex plays a key role in dosage compensation by stimulating expression of male X-chromosome genes. It consists of MSL proteins and two long noncoding RNAs, roX1 and roX2, that are required for spreading of the complex on the chromosome and are redundant in the sense that loss of either does not affect male viability. However, despite rapid evolution, both roX species are present in diverse Drosophilidae species, raising doubts about their full functional redundancy. Thus, we have investigated consequences of deleting roX1 and/or roX2 to probe their specific roles and redundancies in D. melanogaster. We have created a new mutant allele of roX2 and show that roX1 and roX2 have partly separable functions in dosage compensation. In larvae, roX1 is the most abundant variant and the only variant present in the MSL complex when the complex is transmitted (physically associated with the X-chromosome) in mitosis. Loss of roX1 results in reduced expression of the genes on the X-chromosome, while loss of roX2 leads to MSL-independent upregulation of genes with male-biased testis-specific transcription. In roX1 roX2 mutant, gene expression is strongly reduced in a manner that is not related to proximity to high-affinity sites. Our results suggest that high tolerance of mis-expression of the X-chromosome has evolved. We propose that this may be a common property of sex-chromosomes, that dosage compensation is a stochastic process and its precision for each individual gene is regulated by the density of high-affinity sites in the locus.

DNA mismatch repair preferentially safeguards actively transcribed genes

DNA mismatch repair (MMR) is an evolutionally conserved genome maintenance pathway and is well known for its role in maintaining replication fidelity by correcting biosynthetic errors generated during DNA replication. However, recent studies have shown that MMR preferentially protects actively transcribed genes from mutation during both DNA replication and transcription. This review describes the recent discoveries in this area. Potential mechanisms by which MMR safeguards actively transcribed genes are also discussed.

Chromatin conformation and transcriptional activity are permissive regulators of DNA replication initiation in Drosophila

Chromatin structure has emerged as a key contributor to spatial and temporal control over the initiation of DNA replication. However, despite genome-wide correlations between early replication of gene-rich, accessible euchromatin and late replication of gene-poor, inaccessible heterochromatin, a causal relationship between chromatin structure and replication initiation remains elusive. Here, we combined histone gene engineering and whole-genome sequencing in Drosophila to determine how perturbing chromatin structure affects replication initiation. We found that most pericentric heterochromatin remains late replicating in H3K9R mutants, even though H3K9R pericentric heterochromatin is depleted of HP1a, more accessible, and transcriptionally active. These data indicate that HP1a loss, increased chromatin accessibility, and elevated transcription do not result in early replication of heterochromatin. Nevertheless, a small amount of pericentric heterochromatin with increased accessibility replicates earlier in H3K9R mutants. Transcription is de-repressed in these regions of advanced replication but not in those regions of the H3K9R mutant genome that replicate later, suggesting that transcriptional repression may contribute to late replication. We also explored relationships among chromatin, transcription, and replication in euchromatin by analyzing H4K16R mutants. In Drosophila, the X Chromosome gene expression is up-regulated twofold and replicates earlier in XY males than it does in XX females. We found that H4K16R mutation prevents normal male development and abrogates hyperexpression and earlier replication of the male X, consistent with previously established genome-wide correlations between transcription and early replication. In contrast, H4K16R females are viable and fertile, indicating that H4K16 modification is dispensable for genome replication and gene expression.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

H3K36me3-mediated mismatch repair preferentially protects actively transcribed genes from mutation

Histone H3 trimethylation at lysine 36 (H3K36me3) is an important histone mark involved in both transcription elongation and DNA mismatch repair (MMR). It is known that H3K36me3 recruits the mismatch-recognition protein MutSα to replicating chromatin via its physical interaction with MutSα's PWWP domain, but the exact role of H3K36me3 in transcription is undefined. Using ChIP combined with whole-genome DNA sequencing analysis, we demonstrate here that H3K36me3, together with MutSα, is involved in protecting against mutation, preferentially in actively transcribed genomic regions. We found that H3K36me3 and MutSα are much more co-enriched in exons and actively transcribed regions than in introns and nontranscribed regions. The H3K36me3-MutSα co-enrichment correlated with a much lower mutation frequency in exons and actively transcribed regions than in introns and nontranscribed regions. Correspondingly, depleting H3K36me3 or disrupting the H3K36me3-MutSα interaction elevated the spontaneous mutation frequency in actively transcribed genes, but it had little influence on the mutation frequency in nontranscribed or transcriptionally inactive regions. Similarly, H₂O₂-induced mutations, which mainly cause base oxidations, preferentially occurred in actively transcribed genes in MMR-deficient cells. The data presented here suggest that H3K36me3-mediated MMR preferentially safeguards actively transcribed genes not only during replication by efficiently correcting mispairs in early replicating chromatin but also during transcription by directly or indirectly removing DNA lesions associated with a persistently open chromatin structure.

Tethering of CHROMATOR and dCTCF proteins results in decompaction of condensed bands in the Drosophila melanogaster polytene chromosomes but does not affect their transcription and replication timing

Instulator proteins are central to domain organization and gene regulation in the genome. We used ectopic tethering of CHROMATOR (CHRIZ/CHRO) and dCTCF to pre-defined regions of the genome to dissect the influence of these proteins on local chromatin organization, to analyze their interaction with other key chromatin proteins and to evaluate the effects on transcription and replication. Specifically, using UAS-GAL4DBD system, CHRO and dCTCF were artificially recruited into highly compacted polytene chromosome bands that share the features of silent chromatin type known as intercalary heterochromatin (IH). This led to local chromatin decondensation, formation of novel DHSes and recruitment of several "open chromatin" proteins. CHRO tethering resulted in the recruitment of CP190 and Z4 (PZG), whereas dCTCF tethering attracted CHRO, CP190, and Z4. Importantly, formation of a local stretch of open chromatin did not result in the reactivation of silent marker genes yellow and mini-white immediately adjacent to the targeting region (UAS), nor did RNA polII become recruited into this chromatin. The decompacted region retained late replicated, similarly to the wild-type untargeted region.

Genome-Wide Analysis of the Arabidopsis Replication Timing Program

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

A mechanism for epigenetic control of DNA replication

DNA replication origins in hyperacetylated euchromatin fire preferentially during early S phase. However, how acetylation controls DNA replication timing is unknown. TICRR/TRESLIN is an essential protein required for the initiation of DNA replication. Here, we report that TICRR physically interacts with the acetyl-histone binding bromodomain (BRD) and extraterminal (BET) proteins BRD2 and BRD4. Abrogation of this interaction impairs TICRR binding to acetylated chromatin and disrupts normal S-phase progression. Our data reveal a novel function for BET proteins and establish the TICRR-BET interaction as a potential mechanism for epigenetic control of DNA replication.

Repressive Chromatin in Caenorhabditis elegans: Establishment, Composition, and Function

Chromatin is organized and compacted in the nucleus through the association of histones and other proteins, which together control genomic activity. Two broad types of chromatin can be distinguished: euchromatin, which is generally transcriptionally active, and heterochromatin, which is repressed. Here we examine the current state of our understanding of repressed chromatin in Caenorhabditis elegans, focusing on roles of histone modifications associated with repression, such as methylation of histone H3 lysine 9 (H3K9me2/3) or the Polycomb Repressive Complex 2 (MES-2/3/6)-deposited modification H3K27me3, and on proteins that recognize these modifications. Proteins involved in chromatin repression are important for development, and have demonstrated roles in nuclear organization, repetitive element silencing, genome integrity, and the regulation of euchromatin. Additionally, chromatin factors participate in repression with small RNA pathways. Recent findings shed light on heterochromatin function and regulation in C. elegans, and should inform our understanding of repressed chromatin in other animals.

DNA Replication Control During Drosophila Development: Insights into the Onset of S Phase, Replication Initiation, and Fork Progression

Proper control of DNA replication is critical to ensure genomic integrity during cell proliferation. In addition, differential regulation of the DNA replication program during development can change gene copy number to influence cell size and gene expression. Drosophila melanogaster serves as a powerful organism to study the developmental control of DNA replication in various cell cycle contexts in a variety of differentiated cell and tissue types. Additionally, Drosophila has provided several developmentally regulated replication models to dissect the molecular mechanisms that underlie replication-based copy number changes in the genome, which include differential underreplication and gene amplification. Here, we review key findings and our current understanding of the developmental control of DNA replication in the contexts of the archetypal replication program as well as of underreplication and differential gene amplification. We focus on the use of these latter two replication systems to delineate many of the molecular mechanisms that underlie the developmental control of replication initiation and fork elongation.

Genomic Studies Reveal New Aspects of the Biology of DNA Damaging Agents

A flurry of papers has appeared recently to force a rethinking of our understanding of how chemicals, light, and metal complexes damage our genomes. Conventional wisdom was that damaging agents were indiscriminate and it was statistical bad luck, coupled with evolutionary selection, that drove mutational signatures after exposure of DNA to damaging agents. Recent data, however, suggests that primary DNA damage itself does not drive mutational signatures; instead, it is the selectivity of repair pathways on different regions of the genome that is decisive. In particular, genomic regions shielded by transcription factors or packed densely in nucleosomes are poorly repaired by nucleotide excision repair and are far more susceptible to mutation. There are plenty of approved therapies, the mode-of-action of which is to alkylate DNA, and although historically efforts have been focused on understanding how chemicals modify DNA, these new findings suggest that focus should be shifted to understanding genome-wide repair specificities when different types of alkylation damage occur.

Retrotransposons are specified as DNA replication origins in the gene-poor regions of Arabidopsis heterochromatin

Genomic stability depends on faithful genome replication. This is achieved by the concerted activity of thousands of DNA replication origins (ORIs) scattered throughout the genome. The DNA and chromatin features determining ORI specification are not presently known. We have generated a high-resolution genome-wide map of 3230 ORIs in cultured Arabidopsis thaliana cells. Here, we focused on defining the features associated with ORIs in heterochromatin. In pericentromeric gene-poor domains ORIs associate almost exclusively with the retrotransposon class of transposable elements (TEs), in particular of the Gypsy family. ORI activity in retrotransposons occurs independently of TE expression and while maintaining high levels of H3K9me2 and H3K27me1, typical marks of repressed heterochromatin. ORI-TEs largely colocalize with chromatin signatures defining GC-rich heterochromatin. Importantly, TEs with active ORIs contain a local GC content higher than the TEs lacking them. Our results lead us to conclude that ORI colocalization with retrotransposons is determined by their transposition mechanism based on transcription, and a specific chromatin landscape. Our detailed analysis of ORIs responsible for heterochromatin replication has implications on the mechanisms of ORI specification in other multicellular organisms in which retrotransposons are major components of heterochromatin and of the entire genome.

Replicating Large Genomes: Divide and Conquer

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

Repliscan: a tool for classifying replication timing regions

BACKGROUND

Replication timing experiments that use label incorporation and high throughput sequencing produce peaked data similar to ChIP-Seq experiments. However, the differences in experimental design, coverage density, and possible results make traditional ChIP-Seq analysis methods inappropriate for use with replication timing.

RESULTS

To accurately detect and classify regions of replication across the genome, we present Repliscan. Repliscan robustly normalizes, automatically removes outlying and uninformative data points, and classifies Repli-seq signals into discrete combinations of replication signatures. The quality control steps and self-fitting methods make Repliscan generally applicable and more robust than previous methods that classify regions based on thresholds.

CONCLUSIONS

Repliscan is simple and effective to use on organisms with different genome sizes. Even with analysis window sizes as small as 1 kilobase, reliable profiles can be generated with as little as 2.4x coverage.

PEG3 control on the mammalian MSL complex

Peg3 (paternally expressed gene 3) encodes a DNA-binding protein that functions as a transcriptional repressor. Recent studies revealed that PEG3 binds to Msl1 (male-specific lethal 1) and Msl3, the two main components of the MSL complex. In the current study, we investigated potential roles of Peg3 in controlling its downstream genes through H4K16ac, the histone modification by the MSL complex. According to the results, complete removal of PEG3 resulted in up-regulation of Msl1 and Msl3, and subsequently an increase in the global levels of H4K16ac, confirming PEG3 as a transcriptional repressor for MSL during mammalian development. Genome-wide analyses further revealed that about 10% of the entire gene catalogue was affected in the MEF cells lacking PEG3, displaying the increased levels of H4K16ac in their promoter regions. The expression levels of a small subset of the affected genes were up-regulated in the MEF cells lacking PEG3. Interestingly, three Hox clusters also exhibited changes in the levels of H4K16ac, suggesting potential roles of PEG3 and MSL in the regulation of Hox clusters. Overall, the current study reports that Peg3 may control its downstream genes through mammalian MSL.

Order from clutter: selective interactions at mammalian replication origins

Mammalian chromosome duplication progresses in a precise order and is subject to constraints that are often relaxed in developmental disorders and malignancies. Molecular information about the regulation of DNA replication at the chromatin level is lacking because protein complexes that initiate replication seem to bind chromatin indiscriminately. High-throughput sequencing and mathematical modelling have yielded detailed genome-wide replication initiation maps. Combining these maps and models with functional genetic analyses suggests that distinct DNA-protein interactions at subgroups of replication initiation sites (replication origins) modulate the ubiquitous replication machinery and supports an emerging model that delineates how indiscriminate DNA-binding patterns translate into a consistent, organized replication programme.

Variable Autosomal and X Divergence Near and Far from Genes Affects Estimates of Male Mutation Bias in Great Apes

Male mutation bias, when more mutations are passed on via the male germline than via the female germline, is observed across mammals. One common way to infer the magnitude of male mutation bias, α, is to compare levels of neutral sequence divergence between genomic regions that spend different amounts of time in the male and female germline. For great apes, including human, we show that estimates of divergence are reduced in putatively unconstrained regions near genes relative to unconstrained regions far from genes. Divergence increases with increasing distance from genes on both the X chromosome and autosomes, but increases faster on the X chromosome than autosomes. As a result, ratios of X/A divergence increase with increasing distance from genes and corresponding estimates of male mutation bias are significantly higher in intergenic regions near genes versus far from genes. Future studies in other species will need to carefully consider the effect that genomic location will have on estimates of male mutation bias.

Spatio-temporal re-organization of replication foci accompanies replication domain consolidation during human pluripotent stem cell lineage specification

Lineage specification of both mouse and human pluripotent stem cells (PSCs) is accompanied by spatial consolidation of chromosome domains and temporal consolidation of their replication timing. Replication timing and chromatin organization are both established during G1 phase at the timing decision point (TDP). Here, we have developed live cell imaging tools to track spatio-temporal replication domain consolidation during differentiation. First, we demonstrate that the fluorescence ubiquitination cell cycle indicator (Fucci) system is incapable of demarcating G1/S or G2/M cell cycle transitions. Instead, we employ a combination of fluorescent PCNA to monitor S phase progression, cytokinesis to demarcate mitosis, and fluorescent nucleotides to label early and late replication foci and track their 3D organization into sub-nuclear chromatin compartments throughout all cell cycle transitions. We find that, as human PSCs differentiate, the length of S phase devoted to replication of spatially clustered replication foci increases, coincident with global compartmentalization of domains into temporally clustered blocks of chromatin. Importantly, re-localization and anchorage of domains was completed prior to the onset of S phase, even in the context of an abbreviated PSC G1 phase. This approach can also be employed to investigate cell fate transitions in single PSCs, which could be seen to differentiate preferentially from G1 phase. Together, our results establish real-time, live-cell imaging methods for tracking cell cycle transitions during human PSC differentiation that can be applied to study chromosome domain consolidation and other aspects of lineage specification.

Selectivity of ORC binding sites and the relation to replication timing, fragile sites, and deletions in cancers

The origin recognition complex (ORC) binds sites from which DNA replication is initiated. We address ORC binding selectivity in vivo by mapping ∼52,000 ORC2 binding sites throughout the human genome. The ORC binding profile is broader than those of sequence-specific transcription factors, suggesting that ORC is not bound or recruited to specific DNA sequences. Instead, ORC binds nonspecifically to open (DNase I-hypersensitive) regions containing active chromatin marks such as H3 acetylation and H3K4 methylation. ORC sites in early and late replicating regions have similar properties, but there are far more ORC sites in early replicating regions. This suggests that replication timing is due primarily to ORC density and stochastic firing of origins. Computational simulation of stochastic firing from identified ORC sites is in accord with replication timing data. Large genomic regions with a paucity of ORC sites are strongly associated with common fragile sites and recurrent deletions in cancers. We suggest that replication origins, replication timing, and replication-dependent chromosome breaks are determined primarily by the genomic distribution of activator proteins at enhancers and promoters. These activators recruit nucleosome-modifying complexes to create the appropriate chromatin structure that allows ORC binding and subsequent origin firing.

Timing the Drosophila Mid-Blastula Transition: A Cell Cycle-Centered View

At the mid-blastula transition (MBT), externally developing embryos refocus from increasing cell number to elaboration of the body plan. Studies in Drosophila reveal a sequence of changes in regulators of Cyclin:Cdk1 that increasingly restricts the activity of this cell cycle kinase to slow cell cycles during early embryogenesis. By reviewing these events, we provide an outline of the mechanisms slowing the cell cycle at and around the time of MBT. The perspectives developed should provide a guiding paradigm for the study of other MBT changes as the embryo transits from maternal control to a regulatory program centered on the expression of zygotic genes.

Methylation of histone H4 lysine 20 by PR-Set7 ensures the integrity of late replicating sequence domains in Drosophila

The methylation state of lysine 20 on histone H4 (H4K20) has been linked to chromatin compaction, transcription, DNA repair and DNA replication. Monomethylation of H4K20 (H4K20me1) is mediated by the cell cycle-regulated histone methyltransferase PR-Set7. PR-Set7 depletion in mammalian cells results in defective S phase progression and the accumulation of DNA damage, which has been partially attributed to defects in origin selection and activation. However, these studies were limited to only a handful of mammalian origins, and it remains unclear how PR-Set7 and H4K20 methylation impact the replication program on a genomic scale. We employed genetic, cytological, and genomic approaches to better understand the role of PR-Set7 and H4K20 methylation in regulating DNA replication and genome stability in Drosophila cells. We find that deregulation of H4K20 methylation had no impact on origin activation throughout the genome. Instead, depletion of PR-Set7 and loss of H4K20me1 results in the accumulation of DNA damage and an ATR-dependent cell cycle arrest. Coincident with the ATR-dependent cell cycle arrest, we find increased DNA damage that is specifically limited to late replicating regions of the Drosophila genome, suggesting that PR-Set7-mediated monomethylation of H4K20 is critical for maintaining the genomic integrity of late replicating domains.

Replication timing and transcriptional control: beyond cause and effect-part III

DNA replication is essential for faithful transmission of genetic information and is intimately tied to chromosome structure and function. Genome duplication occurs in a defined temporal order known as the replication-timing (RT) program, which is regulated during the cell cycle and development in discrete units referred to as replication domains (RDs). RDs correspond to topologically-associating domains (TADs) and are spatio-temporally compartmentalized in the nucleus. While improvements in experimental tools have begun to reveal glimpses of causality, they have also unveiled complex context-dependent relationships that challenge long recognized correlations of RT to chromatin organization and gene regulation. In particular, RDs/TADs that switch RT during development march to the beat of a different drummer.

Transcriptional Regulators Compete with Nucleosomes Post-replication

Every nucleosome across the genome must be disrupted and reformed when the replication fork passes, but how chromatin organization is re-established following replication is unknown. To address this problem, we have developed Mapping In vivo Nascent Chromatin with EdU and sequencing (MINCE-seq) to characterize the genome-wide location of nucleosomes and other chromatin proteins behind replication forks at high temporal and spatial resolution. We find that the characteristic chromatin landscape at Drosophila promoters and enhancers is lost upon replication. The most conspicuous changes are at promoters that have high levels of RNA polymerase II (RNAPII) stalling and DNA accessibility and show specific enrichment for the BRM remodeler. Enhancer chromatin is also disrupted during replication, suggesting a role for transcription factor (TF) competition in nucleosome re-establishment. Thus, the characteristic nucleosome landscape emerges from a uniformly packaged genome by the action of TFs, RNAPII, and remodelers minutes after replication fork passage.

Mechanisms of nuclear lamina growth in interphase

The nuclear lamina represents a multifunctional platform involved in such diverse yet interconnected processes as spatial organization of the genome, maintenance of mechanical stability of the nucleus, regulation of transcription and replication. Most of lamina activities are exerted through tethering of lamina-associated chromatin domains (LADs) to the nuclear periphery. Yet, the lamina is a dynamic structure demonstrating considerable expansion during the cell cycle to accommodate increased number of LADs formed during DNA replication. We analyzed dynamics of nuclear growth during interphase and changes in lamina structure as a function of cell cycle progression. The nuclear lamina demonstrates steady growth from G1 till G2, while quantitative analysis of lamina meshwork by super-resolution microscopy revealed that microdomain organization of the lamina is maintained, with lamin A and lamin B microdomain periodicity and interdomain gap sizes unchanged. FRAP analysis, in contrast, demonstrated differences in lamin A and B1 exchange rates; the latter showing higher recovery rate in S-phase cells. In order to further analyze the mechanism of lamina growth in interphase, we generated a lamina-free nuclear envelope in living interphase cells by reversible hypotonic shock. The nuclear envelope in nuclear buds formed after such a treatment initially lacked lamins, and analysis of lamina formation revealed striking difference in lamin A and B1 assembly: lamin A reassembled within 30 min post-treatment, whereas lamin B1 did not incorporate into the newly formed lamina at all. We suggest that in somatic cells lamin B1 meshwork growth is coordinated with replication of LADs, and lamin A meshwork assembly seems to be chromatin-independent process.

BRPF3-HBO1 regulates replication origin activation and histone H3K14 acetylation

During DNA replication, thousands of replication origins are activated across the genome. Chromatin architecture contributes to origin specification and usage, yet it remains unclear which chromatin features impact on DNA replication. Here, we perform a RNAi screen for chromatin regulators implicated in replication control by measuring RPA accumulation upon replication stress. We identify six factors required for normal rates of DNA replication and characterize a function of the bromodomain and PHD finger-containing protein 3 (BRPF3) in replication initiation. BRPF3 forms a complex with HBO1 that specifically acetylates histone H3K14, and genomewide analysis shows high enrichment of BRPF3, HBO1 and H3K14ac at ORC1-binding sites and replication origins found in the vicinity of TSSs. Consistent with this, BRPF3 is necessary for H3K14ac at selected origins and efficient origin activation. CDC45 recruitment, but not MCM2-7 loading, is impaired in BRPF3-depleted cells, identifying a BRPF3-dependent function of HBO1 in origin activation that is complementary to its role in licencing. We thus propose that BRPF3-HBO1 acetylation of histone H3K14 around TSS facilitates efficient activation of nearby replication origins.

Large-Scale Chromatin Structure-Function Relationships during the Cell Cycle and Development: Insights from Replication Timing

Chromosome architecture has received a lot of attention since the recent development of genome-scale methods to measure chromatin interactions (Hi-C), enabling the first sequence-based models of chromosome tertiary structure. A view has emerged of chromosomes as a string of structural units (topologically associating domains; TADs) whose boundaries persist through the cell cycle and development. TADs with similar chromatin states tend to aggregate, forming spatially segregated chromatin compartments. However, high-resolution Hi-C has revealed substructure within TADs (subTADs) that poses a challenge for models that attribute significance to structural units at any given scale. More than 20 years ago, the DNA replication field independently identified stable structural (and functional) units of chromosomes (replication foci) as well as spatially segregated chromatin compartments (early and late foci), but lacked the means to link these units to genomic map units. Genome-wide studies of replication timing (RT) have now merged these two disciplines by identifying individual units of replication regulation (replication domains; RDs) that correspond to TADs and are arranged in 3D to form spatiotemporally segregated subnuclear compartments. Furthermore, classifying RDs/TADs by their constitutive versus developmentally regulated RT has revealed distinct classes of chromatin organization, providing unexpected insight into the relationship between large-scale chromosome structure and function.

Topologically associating domains and their long-range contacts are established during early G1 coincident with the establishment of the replication-timing program

Mammalian genomes are partitioned into domains that replicate in a defined temporal order. These domains can replicate at similar times in all cell types (constitutive) or at cell type-specific times (developmental). Genome-wide chromatin conformation capture (Hi-C) has revealed sub-megabase topologically associating domains (TADs), which are the structural counterparts of replication domains. Hi-C also segregates inter-TAD contacts into defined 3D spatial compartments that align precisely to genome-wide replication timing profiles. Determinants of the replication-timing program are re-established during early G1 phase of each cell cycle and lost in G2 phase, but it is not known when TAD structure and inter-TAD contacts are re-established after their elimination during mitosis. Here, we use multiplexed 4C-seq to study dynamic changes in chromatin organization during early G1. We find that both establishment of TADs and their compartmentalization occur during early G1, within the same time frame as establishment of the replication-timing program. Once established, this 3D organization is preserved either after withdrawal into quiescence or for the remainder of interphase including G2 phase, implying 3D structure is not sufficient to maintain replication timing. Finally, we find that developmental domains are less well compartmentalized than constitutive domains and display chromatin properties that distinguish them from early and late constitutive domains. Overall, this study uncovers a strong connection between chromatin re-organization during G1, establishment of replication timing, and its developmental control.

DNA sequence templates adjacent nucleosome and ORC sites at gene amplification origins in Drosophila

Eukaryotic origins of DNA replication are bound by the origin recognition complex (ORC), which scaffolds assembly of a pre-replicative complex (pre-RC) that is then activated to initiate replication. Both pre-RC assembly and activation are strongly influenced by developmental changes to the epigenome, but molecular mechanisms remain incompletely defined. We have been examining the activation of origins responsible for developmental gene amplification in Drosophila. At a specific time in oogenesis, somatic follicle cells transition from genomic replication to a locus-specific replication from six amplicon origins. Previous evidence indicated that these amplicon origins are activated by nucleosome acetylation, but how this affects origin chromatin is unknown. Here, we examine nucleosome position in follicle cells using micrococcal nuclease digestion with Ilumina sequencing. The results indicate that ORC binding sites and other essential origin sequences are nucleosome-depleted regions (NDRs). Nucleosome position at the amplicons was highly similar among developmental stages during which ORC is or is not bound, indicating that being an NDR is not sufficient to specify ORC binding. Importantly, the data suggest that nucleosomes and ORC have opposite preferences for DNA sequence and structure. We propose that nucleosome hyperacetylation promotes pre-RC assembly onto adjacent DNA sequences that are disfavored by nucleosomes but favored by ORC.

Lessons from modENCODE

The modENCODE (Model Organism Encyclopedia of DNA Elements) Consortium aimed to map functional elements-including transcripts, chromatin marks, regulatory factor binding sites, and origins of DNA replication-in the model organisms Drosophila melanogaster and Caenorhabditis elegans. During its five-year span, the consortium conducted more than 2,000 genome-wide assays in developmentally staged animals, dissected tissues, and homogeneous cell lines. Analysis of these data sets provided foundational insights into genome, epigenome, and transcriptome structure and the evolutionary turnover of regulatory pathways. These studies facilitated a comparative analysis with similar data types produced by the ENCODE Consortium for human cells. Genome organization differs drastically in these distant species, and yet quantitative relationships among chromatin state, transcription, and cotranscriptional RNA processing are deeply conserved. Of the many biological discoveries of the modENCODE Consortium, we highlight insights that emerged from integrative studies. We focus on operational and scientific lessons that may aid future projects of similar scale or aims in other, emerging model systems.

Dynamic changes in replication timing and gene expression during lineage specification of human pluripotent stem cells

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.

High-resolution profiling of Drosophila replication start sites reveals a DNA shape and chromatin signature of metazoan origins

At every cell cycle, faithful inheritance of metazoan genomes requires the concerted activation of thousands of DNA replication origins. However, the genetic and chromatin features defining metazoan replication start sites remain largely unknown. Here, we delineate the origin repertoire of the Drosophila genome at high resolution. We address the role of origin-proximal G-quadruplexes and suggest that they transiently stall replication forks in vivo. We dissect the chromatin configuration of replication origins and identify a rich spatial organization of chromatin features at initiation sites. DNA shape and chromatin configurations, not strict sequence motifs, mark and predict origins in higher eukaryotes. We further examine the link between transcription and origin firing and reveal that modulation of origin activity across cell types is intimately linked to cell-type-specific transcriptional programs. Our study unravels conserved origin features and provides unique insights into the relationship among DNA topology, chromatin, transcription, and replication initiation across metazoa.

The hunt for origins of DNA replication in multicellular eukaryotes

Origins of DNA replication (ORIs) occur at defined regions in the genome. Although DNA sequence defines the position of ORIs in budding yeast, the factors for ORI specification remain elusive in metazoa. Several methods have been used recently to map ORIs in metazoan genomes with the hope that features for ORI specification might emerge. These methods are reviewed here with analysis of their advantages and shortcomings. The various factors that may influence ORI selection for initiation of DNA replication are discussed.