Chromatin signatures of the Drosophila replication program

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.

Mechanisms and regulation of DNA replication initiation in eukaryotes

Cellular DNA replication is initiated through the action of multiprotein complexes that recognize replication start sites in the chromosome (termed origins) and facilitate duplex DNA melting within these regions. In a typical cell cycle, initiation occurs only once per origin and each round of replication is tightly coupled to cell division. To avoid aberrant origin firing and re-replication, eukaryotes tightly regulate two events in the initiation process: loading of the replicative helicase, MCM2-7, onto chromatin by the origin recognition complex (ORC), and subsequent activation of the helicase by its incorporation into a complex known as the CMG. Recent work has begun to reveal the details of an orchestrated and sequential exchange of initiation factors on DNA that give rise to a replication-competent complex, the replisome. Here, we review the molecular mechanisms that underpin eukaryotic DNA replication initiation - from selecting replication start sites to replicative helicase loading and activation - and describe how these events are often distinctly regulated across different eukaryotic model organisms.

Kc167, a widely used Drosophila cell line, contains an active primary piRNA pathway

PIWI family proteins bind to small RNAs known as PIWI-interacting RNAs (piRNAs) and play essential roles in the germline by silencing transposons and by promoting germ cell specification and function. Here we report that the widely used Kc167 cell line, derived from Drosophila melanogaster embryos, expresses piRNAs that are loaded to Aub and Piwi. Kc167 piRNAs are produced by a canonical, primary piRNA biogenesis pathway, from phased processing of precursor transcripts by the Zuc endonuclease, Armi helicase, and dGasz mitochondrial scaffold protein. Kc167 piRNAs derive from cytoplasmic transcripts, notably tRNAs and mRNAs, and their abundance correlates with that of parent transcripts. The expression of Aub is robust in Kc167, that of Piwi is modest, while Ago3 is undetectable, explaining the lack of transposon-related piRNA amplification by the Aub-Ago3, ping-pong mechanism. We propose that the default state of the primary piRNA biogenesis machinery is random transcript sampling to allow generation of piRNAs from any transcript, including newly acquired retrotransposons. This state is unmasked in Kc167, likely because they do not express piRNA cluster transcripts in sufficient amounts and do not amplify transposon piRNAs. We use Kc167 to characterize an inactive isoform of Aub protein. Since most Kc167 piRNAs are genic, they can be mapped uniquely to the genome, facilitating computational analyses. Furthermore, because Kc167 is a widely used and well-characterized cell line that is easily amenable to experimental manipulations, we expect that it will serve as an excellent system to study piRNA biogenesis and piRNA-related factors.

Cryo-EM of dynamic protein complexes in eukaryotic DNA replication

DNA replication in Eukaryotes is a highly dynamic process that involves several dozens of proteins. Some of these proteins form stable complexes that are amenable to high-resolution structure determination by cryo-EM, thanks to the recent advent of the direct electron detector and powerful image analysis algorithm. But many of these proteins associate only transiently and flexibly, precluding traditional biochemical purification. We found that direct mixing of the component proteins followed by 2D and 3D image sorting can capture some very weakly interacting complexes. Even at 2D average level and at low resolution, EM images of these flexible complexes can provide important biological insights. It is often necessary to positively identify the feature-of-interest in a low resolution EM structure. We found that systematically fusing or inserting maltose binding protein (MBP) to selected proteins is highly effective in these situations. In this chapter, we describe the EM studies of several protein complexes involved in the eukaryotic DNA replication over the past decade or so. We suggest that some of the approaches used in these studies may be applicable to structural analysis of other biological systems.

Historical Perspective of Eukaryotic DNA Replication

The replication of the genome of a eukaryotic cell is a complex process requiring the ordered assembly of multiprotein replisomes at many chromosomal sites. The process is strictly controlled during the cell cycle to ensure the complete and faithful transmission of genetic information to progeny cells. Our current understanding of the mechanisms of eukaryotic DNA replication has evolved over a period of more than 30 years through the efforts of many investigators. The aim of this perspective is to provide a brief history of the major advances during this period.

Links of genome replication, transcriptional silencing and chromatin dynamics

Genome replication in multicellular organisms involves duplication of both the genetic material and the epigenetic information stored in DNA and histones. In some cases, the DNA replication process provides a window of opportunity for resetting chromatin marks in the genome of the future daughter cells instead of transferring them identical copies. This crucial step of genome replication depends on the correct function of DNA replication factors and the coordination between replication and transcription in proliferating cells. In fact, the histone composition and modification status appears to be intimately associated with the proliferation potential of cells within developing organs. Here we discuss these topics in the light of recent advances in our understanding of how genome replication, transcriptional silencing and chromatin dynamics are coordinated in proliferating cells.

Accessing the Inaccessible: The Organization, Transcription, Replication, and Repair of Heterochromatin in Plants

Eukaryotic genomes often contain large quantities of potentially deleterious sequences, such as transposons. One strategy for mitigating this risk is to package such sequences into so-called constitutive heterochromatin, where the dense chromatin environment is thought to inhibit transcription by excluding transcription factors and RNA polymerase. This type of model makes it tempting to think of heterochromatin as an inert region that is isolated from the rest of the nucleus. Recent work on heterochromatin, however, reveals that it is a dynamic environment. Despite its dense packaging, heterochromatin must remain accessible for a host of processes, including DNA replication and repair, and, paradoxically, transcription. In plants, transcripts produced by specialized RNA polymerases are used to target regions of the genome for silencing via DNA methylation. Thus, the maintenance of heterochromatin requires a careful balancing act of access and exclusion, which is achieved through the action of a host of interrelated pathways.

Transcriptionally Driven DNA Replication Program of the Human Parasite Leishmania major

Faithful inheritance of eukaryotic genomes requires the orchestrated activation of multiple DNA replication origins (ORIs). Although origin firing is mechanistically conserved, how origins are specified and selected for activation varies across different model systems. Here, we provide a complete analysis of the nucleosomal landscape and replication program of the human parasite Leishmania major, building on a better evolutionary understanding of replication organization in Eukarya. We found that active transcription is a driving force for the nucleosomal organization of the L. major genome and that both the spatial and the temporal program of DNA replication can be explained as associated to RNA polymerase kinetics. This simple scenario likely provides flexibility and robustness to deal with the environmental changes that impose alterations in the genetic programs during parasitic life cycle stages. Our findings also suggest that coupling replication initiation to transcription elongation could be an ancient solution used by eukaryotic cells for origin maintenance.

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.

The Histone Variant H3.3 Is Enriched at Drosophila Amplicon Origins but Does Not Mark Them for Activation

Eukaryotic DNA replication begins from multiple origins. The origin recognition complex (ORC) binds origin DNA and scaffolds assembly of a prereplicative complex (pre-RC), which is subsequently activated to initiate DNA replication. In multicellular eukaryotes, origins do not share a strict DNA consensus sequence, and their activity changes in concert with chromatin status during development, but mechanisms are ill-defined. Previous genome-wide analyses in Drosophila and other organisms have revealed a correlation between ORC binding sites and the histone variant H3.3. This correlation suggests that H3.3 may designate origin sites, but this idea has remained untested. To address this question, we examined the enrichment and function of H3.3 at the origins responsible for developmental gene amplification in the somatic follicle cells of the Drosophila ovary. We found that H3.3 is abundant at these amplicon origins. H3.3 levels remained high when replication initiation was blocked, indicating that H3.3 is abundant at the origins before activation of the pre-RC. H3.3 was also enriched at the origins during early oogenesis, raising the possibility that H3.3 bookmarks sites for later amplification. However, flies null mutant for both of the H3.3 genes in Drosophila did not have overt defects in developmental gene amplification or genomic replication, suggesting that H3.3 is not essential for the assembly or activation of the pre-RC at origins. Instead, our results imply that the correlation between H3.3 and ORC sites reflects other chromatin attributes that are important for origin function.

CHANGE POINT ANALYSIS OF HISTONE MODIFICATIONS REVEALS EPIGENETIC BLOCKS LINKING TO PHYSICAL DOMAINS

Histone modification is a vital epigenetic mechanism for transcriptional control in eukaryotes. High-throughput techniques have enabled whole-genome analysis of histone modifications in recent years. However, most studies assume one combination of histone modification invariantly translates to one transcriptional output regardless of local chromatin environment. In this study we hypothesize that, the genome is organized into local domains that manifest similar enrichment pattern of histone modification, which leads to orchestrated regulation of expression of genes with relevant biological functions. We propose a multivariate Bayesian Change Point (BCP) model to segment the Drosophila melanogaster genome into consecutive blocks on the basis of combinatorial patterns of histone marks. By modeling the sparse distribution of histone marks with a zero-inflated Gaussian mixture, our partitions capture local BLOCKs that manifest relatively homogeneous enrichment pattern of histone marks. We further characterized BLOCKs by their transcription levels, distribution of genes, degree of co-regulation and GO enrichment. Our results demonstrate that these BLOCKs, although inferred merely from histone modifications, reveal strong relevance with physical domains, which suggests their important roles in chromatin organization and coordinated gene regulation.

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.

Tools for Targeted Genome Engineering of Established Drosophila Cell Lines

We describe an adaptation of φC31 integrase-mediated targeted cassette exchange for use in Drosophila cell lines. Single copies of an attP-bounded docking platform carrying a GFP-expression marker, with or without insulator elements flanking the attP sites, were inserted by P-element transformation into the Kc167 and Sg4 cell lines; each of the resulting docking-site lines carries a single mapped copy of one of the docking platforms. Vectors for targeted substitution contain a cloning cassette flanked by attB sites. Targeted substitution occurs by integrase-mediated substitution between the attP sites (integrated) and the attB sites (vector). We describe procedures for isolating cells carrying the substitutions and for eliminating the products of secondary off-target events. We demonstrate the technology by integrating a cassette containing a Cu(2+)-inducible mCherry marker, and we report the expression properties of those lines. When compared with clonal lines made by traditional transformation methods, which lead to the illegitimate insertion of tandem arrays, targeted insertion lines give more uniform expression, lower basal expression, and higher induction ratios. Targeted substitution, though intricate, affords results that should greatly improve comparative expression assays-a major emphasis of cell-based studies.

Choosing a suitable method for the identification of replication origins in microbial genomes

As the replication of genomic DNA is arguably the most important task performed by a cell and given that it is controlled at the initiation stage, the events that occur at the replication origin play a central role in the cell cycle. Making sense of DNA replication origins is important for improving our capacity to study cellular processes and functions in the regulation of gene expression, genome integrity in much finer detail. Thus, clearly comprehending the positions and sequences of replication origins which are fundamental to chromosome organization and duplication is the first priority of all. In view of such important roles of replication origins, tremendous work has been aimed at identifying and testing the specificity of replication origins. A number of computational tools based on various skew types have been developed to predict replication origins. Using various in silico approaches such as Ori-Finder, and databases such as DoriC, researchers have predicted the locations of replication origins sites for thousands of bacterial chromosomes and archaeal genomes. Based on the predicted results, we should choose an effective method for identifying and confirming the interactions at origins of replication. Here we describe the main existing experimental methods that aimed to determine the replication origin regions and list some of the many the practical applications of these methods.

DNA replication origin activation in space and time

DNA replication begins with the assembly of pre-replication complexes (pre-RCs) at thousands of DNA replication origins during the G1 phase of the cell cycle. At the G1-S-phase transition, pre-RCs are converted into pre-initiation complexes, in which the replicative helicase is activated, leading to DNA unwinding and initiation of DNA synthesis. However, only a subset of origins are activated during any S phase. Recent insights into the mechanisms underlying this choice reveal how flexibility in origin usage and temporal activation are linked to chromosome structure and organization, cell growth and differentiation, and replication stress.

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.

Links between genome replication and chromatin landscapes

Post-embryonic organogenesis in plants requires the continuous production of cells in the organ primordia, their expansion and a coordinated exit to differentiation. Genome replication is one of the most important processes that occur during the cell cycle, as the maintenance of genomic integrity is of primary relevance for development. As it is chromatin that must be duplicated, a strict coordination occurs between DNA replication, the deposition of new histones, and the introduction of histone modifications and variants. In turn, the chromatin landscape affects several stages during genome replication. Thus, chromatin accessibility is crucial for the initial stages and to specify the location of DNA replication origins with different chromatin signatures. The chromatin landscape also determines the timing of activation during the S phase. Genome replication must occur fully, but only once during each cell cycle. The re-replication avoidance mechanisms rely primarily on restricting the availability of certain replication factors; however, the presence of specific histone modifications are also revealed as contributing to the mechanisms that avoid re-replication, in particular for heterochromatin replication. We provide here an update of genome replication mostly focused on data from Arabidopsis, and the advances that genomic approaches are likely to provide in the coming years. The data available, both in plants and animals, point to the relevance of the chromatin landscape in genome replication, and require a critical evaluation of the existing views about the nature of replication origins, the mechanisms of origin specification and the relevance of epigenetic modifications for genome replication.

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.

Mapping of Replication Origins in the X Inactivation Center of Vole Microtus levis Reveals Extended Replication Initiation Zone

DNA replication initiates at specific positions termed replication origins. Genome-wide studies of human replication origins have shown that origins are organized into replication initiation zones. However, only few replication initiation zones have been described so far. Moreover, few origins were mapped in other mammalian species besides human and mouse. Here we analyzed pattern of short nascent strands in the X inactivation center (XIC) of vole Microtus levis in fibroblasts, trophoblast stem cells, and extraembryonic endoderm stem cells and confirmed origins locations by ChIP approach. We found that replication could be initiated in a significant part of XIC. We also analyzed state of XIC chromatin in these cell types. We compared origin localization in the mouse and vole XIC. Interestingly, origins associated with gene promoters are conserved in these species. The data obtained allow us to suggest that the X inactivation center of M. levis is one extended replication initiation zone.

Drosophila Rif1 is an essential gene and controls late developmental events by direct interaction with PP1-87B

Rif1, identified as a regulator of telomerase in yeast, is an evolutionarily conserved protein and functions in diverse processes including telomere length regulation, epigenetic gene regulation, anti-checkpoint activity, DNA repair and establishing timing of firing at replication origins. Previously we had identified that all Rif1 homologues have PP1 interacting SILK-RVxF motif. In the present study, we show that Drosophila Rif1 is essential for normal fly development and loss of dRif1 impairs temporal regulation of initiation of DNA replication. In multiple tissues dRif1 colocalizes with HP1, a protein known to orchestrate timing of replication in fly. dRif1 associates with chromosomes in a mitotic stage-dependent manner coinciding with dephosphorylation of histones. Ectopic expression of dRif1 causes enlarged larval imaginal discs and early pupal lethality which is completely reversed by co-expression of PP1 87B, the major protein phosphatase in Drosophila, indicating genetic and functional interaction. These findings suggest that dRif1 is an adaptor for PP1 and functions by recruiting PP1 to multiple sites on the chromosome.

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.

Effect of replication on epigenetic memory and consequences on gene transcription

Gene activity in eukaryotes is in part regulated at the level of chromatin through the assembly of local chromatin states that are more or less permissive to transcription. How do these chromatin states achieve their functions and whether or not they contribute to the epigenetic inheritance of the transcriptional program remain to be elucidated. In cycling cells, stability is indeed strongly challenged by the periodic occurrence of replication and cell division. To address this question, we perform simulations of the stochastic dynamics of chromatin states when driven out-of-equilibrium by periodic perturbations. We show how epigenetic memory is significantly affected by the cell cycle length. In addition, we develop a simple model to connect the epigenetic state to the transcriptional state and gene activity. In particular, it suggests that replication may induce transcriptional bursting at repressive loci. Finally, we discuss how our findings-effect of replication and link to gene transcription-have original and deep implications to various biological contexts of epigenetic memory.

Genome-wide chromatin footprinting reveals changes in replication origin architecture induced by pre-RC assembly

Start sites of DNA replication are marked by the origin recognition complex (ORC), which coordinates Mcm2–7 helicase loading to form the prereplicative complex (pre-RC). Belsky et al. “footprinted” nucleosomes, transcription factors, and replication proteins at multiple points during the Saccharomyces cerevisiae cell cycle. This revealed a precise ORC-dependent footprint at 269 origins in G2. A separate class of inefficient origins exhibited protein occupancy only in G1. The local chromatin environment restricts the loading of the Mcm2–7 double hexamer either upstream of or downstream from the ACS.

Start sites of DNA replication are marked by the origin recognition complex (ORC), which coordinates Mcm2–7 helicase loading to form the prereplicative complex (pre-RC). Although pre-RC assembly is well characterized in vitro, the process is poorly understood within the local chromatin environment surrounding replication origins. To reveal how the chromatin architecture modulates origin selection and activation, we “footprinted” nucleosomes, transcription factors, and replication proteins at multiple points during the Saccharomyces cerevisiae cell cycle. Our nucleotide-resolution protein occupancy profiles resolved a precise ORC-dependent footprint at 269 origins in G2. A separate class of inefficient origins exhibited protein occupancy only in G1, suggesting that stable ORC chromatin association in G2 is a determinant of origin efficiency. G1 nucleosome remodeling concomitant with pre-RC assembly expanded the origin nucleosome-free region and enhanced activation efficiency. Finally, the local chromatin environment restricts the loading of the Mcm2–7 double hexamer either upstream of or downstream from the ARS consensus sequence (ACS).

DNA replication and transcription programs respond to the same chromatin cues

DNA replication is a dynamic process that occurs in a temporal order along each of the chromosomes. A consequence of the temporally coordinated activation of replication origins is the establishment of broad domains (>100 kb) that replicate either early or late in S phase. This partitioning of the genome into early and late replication domains is important for maintaining genome stability, gene dosage, and epigenetic inheritance; however, the molecular mechanisms that define and establish these domains are poorly understood. The modENCODE Project provided an opportunity to investigate the chromatin features that define the Drosophila replication timing program in multiple cell lines. The majority of early and late replicating domains in the Drosophila genome were static across all cell lines; however, a small subset of domains was dynamic and exhibited differences in replication timing between the cell lines. Both origin selection and activation contribute to defining the DNA replication program. Our results suggest that static early and late replicating domains were defined at the level of origin selection (ORC binding) and likely mediated by chromatin accessibility. In contrast, dynamic domains exhibited low ORC densities in both cell types, suggesting that origin activation and not origin selection governs the plasticity of the DNA replication program. Finally, we show that the male-specific early replication of the X chromosome is dependent on the dosage compensation complex (DCC), suggesting that the transcription and replication programs respond to the same chromatin cues. Specifically, MOF-mediated hyperacetylation of H4K16 on the X chromosome promotes both the up-regulation of male-specific transcription and origin activation.

Extensive duplication of the Wolbachia DNA in chromosome four of Drosophila ananassae

Background

Lateral gene transfer (LGT) from bacterial Wolbachia endosymbionts has been detected in ~20% of arthropod and nematode genome sequencing projects. Many of these transfers are large and contain a substantial part of the Wolbachia genome.

Results

Here, we re-sequenced three D. ananassae genomes from Asia and the Pacific that contain large LGTs from Wolbachia. We find that multiple copies of the Wolbachia genome are transferred to the Drosophila nuclear genome in all three lines. In the D. ananassae line from Indonesia, the copies of Wolbachia DNA in the nuclear genome are nearly identical in size and sequence yielding an even coverage of mapped reads over the Wolbachia genome. In contrast, the D. ananassae lines from Hawaii and India show an uneven coverage of mapped reads over the Wolbachia genome suggesting that different parts of these LGTs are present in different copy numbers. In the Hawaii line, we find that this LGT is underrepresented in third instar larvae indicative of being heterochromatic. Fluorescence in situ hybridization of mitotic chromosomes confirms that the LGT in the Hawaii line is heterochromatic and represents ~20% of the sequence on chromosome 4 (dot chromosome, Muller element F).

Conclusions

This collection of related lines contain large lateral gene transfers composed of multiple Wolbachia genomes that constitute >2% of the D. ananassae genome (~5 Mbp) and partially explain the abnormally large size of chromosome 4 in D. ananassae.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1097) contains supplementary material, which is available to authorized users.

Tethering of SUUR and HP1 proteins results in delayed replication of euchromatic regions in Drosophila melanogaster polytene chromosomes

We analyze how artificial targeting of Suppressor of Under-Replication (SUUR) and HP1 proteins affects DNA replication in the "open," euchromatic regions. Normally these regions replicate early in the S phase and display no binding of either SUUR or HP1. These proteins were expressed as fusions with DNA-binding domain of GAL4 and recruited to multimerized UAS integrated in three euchromatic sites of the polytene X chromosome: 3B, 8D, and 18B. Using PCNA staining as a marker of ongoing replication, we showed that targeting of SUUR(GAL4DBD) and HP1(GAL4DBD) results in delayed replication of appropriate euchromatic regions. Specifically, replication at these regions starts early, much like in the absence of the fusion proteins; however, replication completion is significantly delayed. Notably, delayed replication was insufficient to induce underreplication. Recruitment of SUUR(GAL4DBD) and HP1(GAL4DBD) had distinct effects on expression of a mini-white reporter, found near UAS. Whereas SUUR(GAL4DBD) had no measurable influence on mini-white expression, HP1(GAL4DBD) targeting silenced mini-white, even in the absence of functional SU(VAR)3-9. Furthermore, recruitment of SUUR(GAL4DBD) and HP1(GAL4DBD) had distinct effects on the protein composition of target regions. HP1(GAL4DBD) but not SUUR(GAL4DBD) could displace an open chromatin marker, CHRIZ, from the tethering sites.

CNV instability associated with DNA replication dynamics: evidence for replicative mechanisms in CNV mutagenesis

Copy number variation (CNV) in the human genome is of vital importance to human health and evolution of our species. However, much of the molecular basis of CNV mutagenesis remains to be elucidated. Considering the DNA replication model of 'fork stalling and template switching' for CNV formation, we hypothesized that replication fork progression could be important for CNV mutagenesis. However, molecular assays of replication fork progression at the genome level are technically challenging. Instead, we conducted an estimation of DNA replication dynamics, as the statistic R, using the readily available data of replication timing. Small R-values can reflect 'slowed' replication, which could result from less fork initiation, reduced fork speed or fork barriers. We generated genome-wide profiles of R in the genomes of human, mouse and Drosophila. Intriguingly, the CNV breakpoints in all three genomes showed significantly biased distributions toward the genomic regions with small R-values, suggesting potential replication stress-induced CNV instability. Notably, among the human CNVs with distinct breakpoint junction characteristics, the homology-mediated and VNTR-mediated CNVs contribute the most to the correlation between CNV instability and the statistic R, consistent with the recent findings in the C. elegans and yeast genomes of repeat-induced DNA replication error and consequent CNV formation. The statistic R may reflect both replication stress and the effect of local genome architecture on fork progression. Our concordant observations suggest an important role for DNA replicative mechanisms in CNV mutagenesis and genome instability.

DNA copy-number control through inhibition of replication fork progression

Proper control of DNA replication is essential to ensure faithful transmission of genetic material and prevent chromosomal aberrations that can drive cancer progression and developmental disorders. DNA replication is regulated primarily at the level of initiation and is under strict cell-cycle regulation. Importantly, DNA replication is highly influenced by developmental cues. In Drosophila, specific regions of the genome are repressed for DNA replication during differentiation by the SNF2 domain-containing protein SUUR through an unknown mechanism. We demonstrate that SUUR is recruited to active replication forks and mediates the repression of DNA replication by directly inhibiting replication fork progression instead of functioning as a replication fork barrier. Mass spectrometry identification of SUUR-associated proteins identified the replicative helicase member CDC45 as a SUUR-associated protein, supporting a role for SUUR directly at replication forks. Our results reveal that control of eukaryotic DNA copy number can occur through the inhibition of replication fork progression.

Incomplete replication generates somatic DNA alterations within Drosophila polytene salivary gland cells

DNA replication remains unfinished in many Drosophila polyploid cells, which harbor disproportionately fewer copies of late-replicating chromosomal regions. Using NextGen sequencing of DNA from giant polytene cells of the larval salivary gland, Yarosh and Spradling show that sporadic, incomplete replication during the endocycle S phase alters the Drosophila genome at thousands of sites that differ in every cell; similar events occur in the ovary. The authors propose that the extensive somatic DNA instability described here underlies position effect variegation and molds the structure of polytene chromosomes.

DNA replication remains unfinished in many Drosophila polyploid cells, which harbor disproportionately fewer copies of late-replicating chromosomal regions. By analyzing paired-end high-throughput sequence data from polytene larval salivary gland cells, we define 112 underreplicated (UR) euchromatic regions 60–480 kb in size. To determine the effects of underreplication on genome integrity, we analyzed anomalous read pairs and breakpoint reads throughout the euchromatic genome. Each UR euchromatic region contains many different deletions 10–500 kb in size, while very few deletions are present in fully replicated chromosome regions or UR zones from embryo DNA. Thus, during endocycles, stalled forks within UR regions break and undergo local repair instead of remaining stable and generating nested forks. As a result, each salivary gland cell contains hundreds of unique deletions that account for their copy number reductions. Similar UR regions and deletions were observed in ovarian DNA, suggesting that incomplete replication, fork breakage, and repair occur widely in polytene cells. UR regions are enriched in genes encoding immunoglobulin superfamily proteins and contain many neurally expressed and homeotic genes. We suggest that the extensive somatic DNA instability described here underlies position effect variegation, molds the structure of polytene chromosomes, and should be investigated for possible functions.

DNA replication in nurse cell polytene chromosomes of Drosophila melanogaster otu mutants

Drosophila cell lines are used extensively to study replication timing, yet data about DNA replication in larval and adult tissues are extremely limited. To address this gap, we traced DNA replication in polytene chromosomes from nurse cells of Drosophila melanogaster otu mutants using bromodeoxyuridine incorporation. Importantly, nurse cells are of female germline origin, unlike the classical larval salivary glands, that are somatic. In contrast to salivary gland polytene chromosomes, where replication begins simultaneously across all puffs and interbands, replication in nurse cells is first observed at several specific chromosomal regions. For instance, in the chromosome 2L, these include the regions 31B-E and 37E and proximal parts of 34B and 35B, with the rest of the decondensed chromosomal regions joining replication process a little later. We observed that replication timing of pericentric heterochromatin in nurse cells was shifted from late S phase to early and mid stages. Curiously, chromosome 4 may represent a special domain of the genome, as it replicates on its own schedule which is uncoupled from the rest of the chromosomes. Finally, we report that SUUR protein, an established marker of late replication in salivary gland polytene chromosomes, does not always colocalize with late-replicating regions in nurse cells.

Origin replication complex binding, nucleosome depletion patterns, and a primary sequence motif can predict origins of replication in a genome with epigenetic centromeres

Origins of DNA replication are key genetic elements, yet their identification remains elusive in most organisms. In previous work, we found that centromeres contain origins of replication (ORIs) that are determined epigenetically in the pathogenic yeast Candida albicans. In this study, we used origin recognition complex (ORC) binding and nucleosome occupancy patterns in Saccharomyces cerevisiae and Kluyveromyces lactis to train a machine learning algorithm to predict the position of active arm (noncentromeric) origins in the C. albicans genome. The model identified bona fide active origins as determined by the presence of replication intermediates on nondenaturing two-dimensional (2D) gels. Importantly, these origins function at their native chromosomal loci and also as autonomously replicating sequences (ARSs) on a linear plasmid. A “mini-ARS screen” identified at least one and often two ARS regions of ≥100 bp within each bona fide origin. Furthermore, a 15-bp AC-rich consensus motif was associated with the predicted origins and conferred autonomous replicating activity to the mini-ARSs. Thus, while centromeres and the origins associated with them are epigenetic, arm origins are dependent upon critical DNA features, such as a binding site for ORC and a propensity for nucleosome exclusion.

DNA replication machinery is highly conserved, yet the definition of exactly what specifies a replication origin differs in different species. Here, we utilized computational genomics to predict origin locations in Candida albicans by combining locations of binding sites for the conserved origin replication complex, necessary for replication initiation, together with chromatin organization patterns. We identified predicted sequences that exhibited bona fide origin function and developed a linear plasmid assay to delimit the DNA fragments necessary for origin function. Additionally, we found that a short AC-rich motif, which is enriched in predicted origins, is required for origin function. Thus, we demonstrated a new machine learning paradigm for identification of potential origins from a genome with no prior information. Furthermore, this work suggests that C. albicans has two different types of origins: “hard-wired” arm origins that rely upon specific sequence motifs and “epigenetic” centromeric origins that are recruited to kinetochores in a sequence-independent manner.

pENCODE: a plant encyclopedia of DNA elements

ENCODE projects exist for many eukaryotes, including humans, but as of yet no defined project exists for plants. A plant ENCODE would be invaluable to the research community and could be more readily produced than its metazoan equivalents by capitalizing on the preexisting infrastructure provided from similar projects. Collecting and normalizing plant epigenomic data for a range of species will facilitate hypothesis generation, cross-species comparisons, annotation of genomes, and an understanding of epigenomic functions throughout plant evolution. Here, we discuss the need for such a project, outline the challenges it faces, and suggest ways forward to build a plant ENCODE.

Genetic organization of interphase chromosome bands and interbands in Drosophila melanogaster

Drosophila melanogaster polytene chromosomes display specific banding pattern; the underlying genetic organization of this pattern has remained elusive for many years. In the present paper, we analyze 32 cytology-mapped polytene chromosome interbands. We estimated molecular locations of these interbands, described their molecular and genetic organization and demonstrate that polytene chromosome interbands contain the 5' ends of housekeeping genes. As a rule, interbands display preferential "head-to-head" orientation of genes. They are enriched for "broad" class promoters characteristic of housekeeping genes and associate with open chromatin proteins and Origin Recognition Complex (ORC) components. In two regions, 10A and 100B, coding sequences of genes whose 5'-ends reside in interbands map to constantly loosely compacted, early-replicating, so-called "grey" bands. Comparison of expression patterns of genes mapping to late-replicating dense bands vs genes whose promoter regions map to interbands shows that the former are generally tissue-specific, whereas the latter are represented by ubiquitously active genes. Analysis of RNA-seq data (modENCODE-FlyBase) indicates that transcripts from interband-mapping genes are present in most tissues and cell lines studied, across most developmental stages and upon various treatment conditions. We developed a special algorithm to computationally process protein localization data generated by the modENCODE project and show that Drosophila genome has about 5700 sites that demonstrate all the features shared by the interbands cytologically mapped to date.

Underreplicated regions in Drosophila melanogaster are enriched with fast-evolving genes and highly conserved noncoding sequences

Many late replicating regions are underreplicated in polytene chromosomes of Drosophila melanogaster. These regions contain silenced chromatin and overlap long syntenic blocks of conserved gene order in drosophilids. In this report we show that in D. melanogaster the underreplicated regions are enriched with fast-evolving genes lacking homologs in distant species such as mosquito or human, indicating that the phylogenetic conservation of genes correlates with replication timing and chromatin status. Drosophila genes without human homologs located in the underreplicated regions have higher nonsynonymous substitution rate and tend to encode shorter proteins when compared with those in the adjacent regions. At the same time, the underreplicated regions are enriched with ultraconserved elements and highly conserved noncoding sequences, especially in introns of very long genes indicating the presence of an extensive regulatory network that may be responsible for the conservation of gene order in these regions. The regions have a modest preference for long noncoding RNAs but are depleted for small nucleolar RNAs, microRNAs, and transfer RNAs. Our results demonstrate that the underreplicated regions have a specific genic composition and distinct pattern of evolution.

Cell lines

We review the properties and uses of cell lines in Drosophila research, emphasizing the variety of lines, the large body of genomic and transcriptional data available for many of the lines, and the variety of ways the lines have been used to provide tools for and insights into the developmental, molecular, and cell biology of Drosophila and mammals.

Chromatin structure and replication origins: determinants of chromosome replication and nuclear organization

The DNA replication program is, in part, determined by the epigenetic landscape that governs local chromosome architecture and directs chromosome duplication. Replication must coordinate with other biochemical processes occurring concomitantly on chromatin, such as transcription and remodeling, to insure accurate duplication of both genetic and epigenetic features and to preserve genomic stability. The importance of genome architecture and chromatin looping in coordinating cellular processes on chromatin is illustrated by two recent sets of discoveries. First, chromatin-associated proteins that are not part of the core replication machinery were shown to affect the timing of DNA replication. These chromatin-associated proteins could be working in concert, or perhaps in competition, with the transcriptional machinery and with chromatin modifiers to determine the spatial and temporal organization of replication initiation events. Second, epigenetic interactions are mediated by DNA sequences that determine chromosomal replication. In this review, we summarize recent findings and current models linking spatial and temporal regulation of the replication program with epigenetic signaling. We discuss these issues in the context of the genome's three-dimensional structure with an emphasis on events occurring during the initiation of DNA replication.

Chromatin meets the cell cycle

The cell cycle is one of the most comprehensively studied biological processes, due primarily to its significance in growth and development, and its deregulation in many human disorders. Studies using a diverse set of model organisms, including yeast, worms, flies, frogs, mammals, and plants, have greatly expanded our knowledge of the cell cycle and have contributed to the universally accepted view of how the basic cell cycle machinery is regulated. In addition to the oscillating activity of various cyclin-dependent kinase (CDK)-cyclin complexes, a plethora of proteins affecting various aspects of chromatin dynamics has been shown to be essential for cell proliferation during plant development. Furthermore, it was reported recently that core cell cycle regulators control gene expression by modifying histone patterns. This review focuses on the intimate relationship between the cell cycle and chromatin. It describes the dynamics and functions of chromatin structures throughout cell cycle progression and discusses the role of heterochromatin as a barrier against re-replication and endoreduplication. It also proposes that core plant cell cycle regulators control gene expression in a manner similar to that described in mammals. At present, our challenge in plants is to define the complete set of effectors and actors that coordinate cell cycle progression and chromatin structure and to understand better the functional interplay between these two processes.

The spatiotemporal program of DNA replication is associated with specific combinations of chromatin marks in human cells

The duplication of mammalian genomes is under the control of a spatiotemporal program that orchestrates the positioning and the timing of firing of replication origins. The molecular mechanisms coordinating the activation of about predicted origins remain poorly understood, partly due to the intrinsic rarity of replication bubbles, making it difficult to purify short nascent strands (SNS). The precise identification of origins based on the high-throughput sequencing of SNS constitutes a new methodological challenge. We propose a new statistical method with a controlled resolution, adapted to the detection of replication origins from SNS data. We detected an average of 80,000 replication origins in different cell lines. To evaluate the consistency between different protocols, we compared SNS detections with bubble trapping detections. This comparison demonstrated a good agreement between genome-wide methods, with 65% of SNS-detected origins validated by bubble trapping, and 44% of bubble trapping origins validated by SNS origins, when compared at the same resolution. We investigated the interplay between the spatial and the temporal programs of replication at fine scales. We show that most of the origins detected in regions replicated in early S phase are shared by all the cell lines investigated whereas cell-type-specific origins tend to be replicated in late S phase. We shed a new light on the key role of CpG islands, by showing that 80% of the origins associated with CGIs are constitutive. Our results further show that at least 76% of CGIs are origins of replication. The analysis of associations with chromatin marks at different timing of cell division revealed new potential epigenetic regulators driving the spatiotemporal activity of replication origins. We highlight the potential role of H4K20me1 and H3K27me3, the coupling of which is correlated with increased efficiency of replication origins, clearly identifying those marks as potential key regulators of replication origins.

Replication is the mechanism by which genomes are duplicated into two exact copies. Genomic stability is under the control of a spatiotemporal program that orchestrates both the positioning and the timing of firing of about 50,000 replication starting points, also called replication origins. Replication bubbles found at origins have been very difficult to map due to their short lifespan. Moreover, with the flood of data characterizing new sequencing technologies, the precise statistical analysis of replication data has become an additional challenge. We propose a new method to map replication origins on the human genome, and we assess the reliability of our finding using experimental validation and comparison with origins maps obtained by bubble trapping. This fine mapping then allowed us to identify potential regulators of the replication dynamics. Our study highlights the key role of CpG Islands and identifies new potential epigenetic regulators (methylation of lysine 4 on histone H4, and tri-methylation of lysine 27 on histone H3) whose coupling is correlated with an increase in the efficiency of replication origins, suggesting those marks as potential key regulators of replication. Overall, our study defines new potentially important pathways that might regulate the sequential firing of origins during genome duplication.

A new light on DNA replication from the inactive X chromosome

While large portions of the mammalian genome are known to replicate sequentially in a distinct, tissue-specific order, recent studies suggest that the inactive X chromosome is duplicated rapidly via random, synchronous DNA synthesis at numerous adjacent regions. The rapid duplication of the inactive X chromosome was observed in high-resolution studies visualizing DNA replication patterns in the nucleus, and by allele-specific DNA sequencing studies measuring the extent of DNA synthesis. These studies conclude that inactive X chromosomes complete replication earlier than previously thought and suggest that the strict order of DNA replication detected in the majority of genomic regions is not preserved in non-transcribed, "silent" chromatin. These observations alter current concepts about the regulation of DNA replication in non-transcribed portions of the genome in general and in the inactive X-chromosome in particular.

GC-rich DNA elements enable replication origin activity in the methylotrophic yeast Pichia pastoris

The well-studied DNA replication origins of the model budding and fission yeasts are A/T-rich elements. However, unlike their yeast counterparts, both plant and metazoan origins are G/C-rich and are associated with transcription start sites. Here we show that an industrially important methylotrophic budding yeast, Pichia pastoris, simultaneously employs at least two types of replication origins--a G/C-rich type associated with transcription start sites and an A/T-rich type more reminiscent of typical budding and fission yeast origins. We used a suite of massively parallel sequencing tools to map and dissect P. pastoris origins comprehensively, to measure their replication dynamics, and to assay the global positioning of nucleosomes across the genome. Our results suggest that some functional overlap exists between promoter sequences and G/C-rich replication origins in P. pastoris and imply an evolutionary bifurcation of the modes of replication initiation.

Complex correlations: replication timing and mutational landscapes during cancer and genome evolution

A recent flurry of reports correlates replication timing (RT) with mutation rates during both evolution and cancer. Specifically, point mutations and copy number losses correlate with late replication, while copy number gains and other rearrangements correlate with early replication. In some cases, plausible mechanisms have been proposed. Point mutation rates may reflect temporal variation in repair mechanisms. Transcription-induced double-strand breaks are expected to occur in transcriptionally active early replicating chromatin. Fusion partners are generally in close proximity, and chromatin in close proximity replicates at similar times. However, temporal enrichment of copy number gains and losses remains an enigma. Moreover, many conclusions are compromised by a lack of matched RT and sequence datasets, the filtering out of developmental variation in RT, and the use of somatic cell lines to make inferences about germline evolution.

Combinatorial modeling of chromatin features quantitatively predicts DNA replication timing in Drosophila

In metazoans, each cell type follows a characteristic, spatio-temporally regulated DNA replication program. Histone modifications (HMs) and chromatin binding proteins (CBPs) are fundamental for a faithful progression and completion of this process. However, no individual HM is strictly indispensable for origin function, suggesting that HMs may act combinatorially in analogy to the histone code hypothesis for transcriptional regulation. In contrast to gene expression however, the relationship between combinations of chromatin features and DNA replication timing has not yet been demonstrated. Here, by exploiting a comprehensive data collection consisting of 95 CBPs and HMs we investigated their combinatorial potential for the prediction of DNA replication timing in Drosophila using quantitative statistical models. We found that while combinations of CBPs exhibit moderate predictive power for replication timing, pairwise interactions between HMs lead to accurate predictions genome-wide that can be locally further improved by CBPs. Independent feature importance and model analyses led us to derive a simplified, biologically interpretable model of the relationship between chromatin landscape and replication timing reaching 80% of the full model accuracy using six model terms. Finally, we show that pairwise combinations of HMs are able to predict differential DNA replication timing across different cell types. All in all, our work provides support to the existence of combinatorial HM patterns for DNA replication and reveal cell-type independent key elements thereof, whose experimental investigation might contribute to elucidate the regulatory mode of this fundamental cellular process.

Fungal genes in context: genome architecture reflects regulatory complexity and function

Gene context determines gene expression, with local chromosomal environment most influential. Comparative genomic analysis is often limited in scope to conserved or divergent gene and protein families, and fungi are well suited to this approach with low functional redundancy and relatively streamlined genomes. We show here that one aspect of gene context, the amount of potential upstream regulatory sequence maintained through evolution, is highly predictive of both molecular function and biological process in diverse fungi. Orthologs with large upstream intergenic regions (UIRs) are strongly enriched in information processing functions, such as signal transduction and sequence-specific DNA binding, and, in the genus Aspergillus, include the majority of experimentally studied, high-level developmental and metabolic transcriptional regulators. Many uncharacterized genes are also present in this class and, by implication, may be of similar importance. Large intergenic regions also share two novel sequence characteristics, currently of unknown significance: they are enriched for plus-strand polypyrimidine tracts and an information-rich, putative regulatory motif that was present in the last common ancestor of the Pezizomycotina. Systematic consideration of gene UIR in comparative genomics, particularly for poorly characterized species, could help reveal organisms’ regulatory priorities.

Epigenetic landscape for initiation of DNA replication

The key genetic process of DNA replication is initiated at specific sites referred to as replication origins. In eukaryotes, origins of DNA replication are not specified by a defined nucleotide sequence. Recent studies have shown that the structural context and topology of DNA sequence, chromatin features, and its transcriptional activity play an important role in origin choice. During differentiation and development, significant changes in chromatin organization and transcription occur, influencing origin activity and choice. In the last few years, a number of different genome-wide studies have broadened the understanding of replication origin regulation. In this review, we discuss the epigenetic factors and mechanisms that modulate origin choice and firing.

A unique epigenetic signature is associated with active DNA replication loci in human embryonic stem cells

The cellular epigenetic landscape changes as pluripotent stem cells differentiate to somatic cells or when differentiated cells transform to a cancerous state. These epigenetic changes are commonly correlated with differences in gene expression. Whether active DNA replication is also associated with distinct chromatin environments in these developmentally and phenotypically diverse cell types has not been known. Here, we used BrdU-seq to map active DNA replication loci in human embryonic stem cells (hESCs), normal primary fibroblasts and a cancer cell line, and correlated these maps to the epigenome. In all cell lines, the majority of BrdU peaks were enriched in euchromatin and at DNA repetitive elements, especially at microsatellite repeats, and coincided with previously determined replication origins. The most prominent BrdU peaks were shared between all cells but a sizable fraction of the peaks were specific to each cell type and associated with cell type-specific genes. Surprisingly, the BrdU peaks that were common to all cell lines were associated with H3K18ac, H3K56ac, and H4K20me1 histone marks only in hESCs but not in normal fibroblasts or cancer cells. Depletion of the histone acetyltransferases for H3K18 and H3K56 dramatically decreased the number and intensity of BrdU peaks in hESCs. Our data reveal a unique epigenetic signature that distinguishes active replication loci in hESCs from normal somatic or malignant cells.

DNA replication origins

The onset of genomic DNA synthesis requires precise interactions of specialized initiator proteins with DNA at sites where the replication machinery can be loaded. These sites, defined as replication origins, are found at a few unique locations in all of the prokaryotic chromosomes examined so far. However, replication origins are dispersed among tens of thousands of loci in metazoan chromosomes, thereby raising questions regarding the role of specific nucleotide sequences and chromatin environment in origin selection and the mechanisms used by initiators to recognize replication origins. Close examination of bacterial and archaeal replication origins reveals an array of DNA sequence motifs that position individual initiator protein molecules and promote initiator oligomerization on origin DNA. Conversely, the need for specific recognition sequences in eukaryotic replication origins is relaxed. In fact, the primary rule for origin selection appears to be flexibility, a feature that is modulated either by structural elements or by epigenetic mechanisms at least partly linked to the organization of the genome for gene expression.

DNA replication timing

Patterns of replication within eukaryotic genomes correlate with gene expression, chromatin structure, and genome evolution. Recent advances in genome-scale mapping of replication kinetics have allowed these correlations to be explored in many species, cell types, and growth conditions, and these large data sets have allowed quantitative and computational analyses. One striking new correlation to emerge from these analyses is between replication timing and the three-dimensional structure of chromosomes. This correlation, which is significantly stronger than with any single histone modification or chromosome-binding protein, suggests that replication timing is controlled at the level of chromosomal domains. This conclusion dovetails with parallel work on the heterogeneity of origin firing and the competition between origins for limiting activators to suggest a model in which the stochastic probability of individual origin firing is modulated by chromosomal domain structure to produce patterns of replication. Whether these patterns have inherent biological functions or simply reflect higher-order genome structure is an open question.

[Regulation of DNA replication timing]

Although distinct chromatin types have been long known to replicate at different timepoints of S phase, fine replication control has only recently become considered as an epigenetic phenomenon. It is now clear that in course of differentiation significant changes in genome replication timing occur, and these changes are intimately linked with the changes in transcriptional activity and nuclear architecture. Temporally coordinate replication is organized spatially into discrete units having specific chromosomal organization and function. Even though the functional aspects of such tight control of replication timing remain to be explored, one can confidently consider the replication program as yet another fundamental feature characteristic of the given differentiation state. The present review touches upon the molecular mechanisms of spatial and temporal control of replication timing, involving individual replication origins as well as large chromatin domains.

Single-molecule analysis of combinatorial epigenomic states in normal and tumor cells

Proper placement of epigenetic marks on DNA and histones is fundamental to normal development, and perturbations contribute to a variety of disease states. Combinations of marks act together to control gene expression; therefore, detecting their colocalization is important, but because of technical challenges, such measurements are rarely reported. Instead, measurements of epigenetic marks are typically performed one at a time in a population of cells, and their colocalization is inferred by association. Here, we describe a single-molecule analytical approach that can perform direct detection of multiple epigenetic marks simultaneously and use it to identify mechanisms coordinating placement of three gene silencing marks, trimethylated histone H3 lysine 9, lysine 27 (H3K9me3, H3K27me3), and cytosine methylation (mC), in the normal and cancer genome. We show that H3K9me3 and mC are present together on individual chromatin fragments in mouse embryonic stem cells and that half of the H3K9me3 marks require mC for their placement. In contrast, mC and H3K27me3 coincidence is rare, and in fact, mC antagonizes H3K27me3 in both embryonic stem cells and primary mouse fibroblasts, indicating this antagonism is shared among primary cells. However, upon immortalization or tumorigenic transformation of mouse fibroblasts, mC is required for complete H3K27me3 placement. Importantly, in human promyelocytic cells, H3K27me3 is also dependent on mC. Because aberrant placement of gene silencing marks at tumor suppressor genes contributes to tumor progression, the improper dependency of H3K27me3 by mC in immortalized cells is likely to be fundamental to cancer. Our platform can enable other studies involving coordination of epigenetic marks and leverage efforts to discover disease biomarkers and epigenome-modifying drugs.