A question of timing: emerging links between transcription and replication

Rif1-Dependent Control of Replication Timing

Successful duplication of the genome requires the accurate replication of billions of base pairs of DNA within a relatively short time frame. Failure to accurately replicate the genome results in genomic instability and a host of diseases. To faithfully and rapidly replicate the genome, DNA replication must be tightly regulated and coordinated with many other nuclear processes. These regulations, however, must also be flexible as replication kinetics can change through development and differentiation. Exactly how DNA replication is regulated and how this regulation changes through development is an active field of research. One aspect of genome duplication where much remains to be discovered is replication timing (RT), which dictates when each segment of the genome is replicated during S phase. All organisms display some level of RT, yet the precise mechanisms that govern RT remain are not fully understood. The study of Rif1, a protein that actively regulates RT from yeast to humans, provides a key to unlock the underlying molecular mechanisms controlling RT. The paradigm for Rif1 function is to delay helicase activation within certain regions of the genome, causing these regions to replicate late in S phase. Many questions, however, remain about the intricacies of Rif1 function. Here, we review the current models for the activity of Rif1 with the goal of trying to understand how Rif1 functions to establish the RT program.

Epigenetic regulation of replication origin assembly: A role for histone H1 and chromatin remodeling factors

During early embryonic development in several metazoans, accurate DNA replication is ensured by high number of replication origins. This guarantees rapid genome duplication coordinated with fast cell divisions. In Xenopus laevis embryos this program switches to one with a lower number of origins at a developmental stage known as mid-blastula transition (MBT) when cell cycle length increases and gene transcription starts. Consistent with this regulation, somatic nuclei replicate poorly when transferred to eggs, suggesting the existence of an epigenetic memory suppressing replication assembly origins at all available sites. Recently, it was shown that histone H1 imposes a non-permissive chromatin configuration preventing replication origin assembly on somatic nuclei. This somatic state can be erased by SSRP1, a subunit of the FACT complex. Here, we further develop the hypothesis that this novel form of epigenetic memory might impact on different areas of vertebrate biology going from nuclear reprogramming to cancer development.

Thoracic underreplication in Drosophila species estimates a minimum genome size and the dynamics of added DNA

Many cells in the thorax of Drosophila were found to stall during replication, a phenomenon known as underreplication. Unlike underreplication in nuclei of salivary and follicle cells, this stall occurs with less than one complete round of replication. This stall point allows precise estimations of early-replicating euchromatin and late-replicating heterochromatin regions, providing a powerful tool to investigate the dynamics of structural change across the genome. We measure underreplication in 132 species across the Drosophila genus and leverage these data to propose a model for estimating the rate at which additional DNA is accumulated as heterochromatin and euchromatin and also predict the minimum genome size for Drosophila. According to comparative phylogenetic approaches, the rates of change of heterochromatin differ strikingly between Drosophila subgenera. Although these subgenera differ in karyotype, there were no differences by chromosome number, suggesting other structural changes may influence accumulation of heterochromatin. Measurements were taken for both sexes, allowing the visualization of genome size and heterochromatin changes for the hypothetical path of XY sex chromosome differentiation. Additionally, the model presented here estimates a minimum genome size in Sophophora remarkably close to the smallest insect genome measured to date, in a species over 200 million years diverged from Drosophila.

Germline DNA replication timing shapes mammalian genome composition

Mammalian DNA replication is a highly organized and regulated process. Large, Mb-sized regions are replicated at defined times along S-phase. Replication Timing (RT) is thought to play a role in shaping the mammalian genome by affecting mutation rates. Previous analyses relied on somatic RT profiles. However, only germline mutations are passed on to offspring and affect genomic composition. Therefore, germ cell RT information is necessary to evaluate the influences of RT on the mammalian genome. We adapted the RT mapping technique for limited amounts of cells, and measured RT from two stages in the mouse germline - primordial germ cells (PGCs) and spermatogonial stem cells (SSCs). RT in germline cells exhibited stronger correlations to both mutation rate and recombination hotspots density than those of RT in somatic tissues, emphasizing the importance of using correct tissues-of-origin for RT profiling. Germline RT maps exhibited stronger correlations to additional genetic features including GC-content, transposable elements (SINEs and LINEs), and gene density. GC content stratification and multiple regression analysis revealed independent contributions of RT to SINE, gene, mutation, and recombination hotspot densities. Together, our results establish a central role for RT in shaping multiple levels of mammalian genome composition.

DNA Replication Timing Enters the Single-Cell Era

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

Differences between the genomes of lymphoblastoid cell lines and blood-derived samples

Lymphoblastoid cell lines (LCLs) represent a convenient research tool for expanding the amount of biologic material available from an individual. LCLs are commonly used as reference materials, most notably from the Genome in a Bottle Consortium. However, the question remains how faithfully LCL-derived genome assemblies represent the germline genome of the donor individual as compared to the genome assemblies derived from peripheral blood mononuclear cells. We present an in-depth comparison of a large collection of LCL- and peripheral blood mononuclear cell-derived genomes in terms of distributions of coverage and copy number alterations. We found significant differences in the depth of coverage and copy number calls, which may be driven by differential replication timing. Importantly, these copy number changes preferentially affect regions closer to genes and with higher GC content. This suggests that genomic studies based on LCLs may display locus-specific biases, and that conclusions based on analysis of depth of coverage and copy number variation may require further scrutiny.

The Genetic Program of Pancreatic β-Cell Replication In Vivo

The molecular program underlying infrequent replication of pancreatic β-cells remains largely inaccessible. Using transgenic mice expressing green fluorescent protein in cycling cells, we sorted live, replicating β-cells and determined their transcriptome. Replicating β-cells upregulate hundreds of proliferation-related genes, along with many novel putative cell cycle components. Strikingly, genes involved in β-cell functions, namely, glucose sensing and insulin secretion, were repressed. Further studies using single-molecule RNA in situ hybridization revealed that in fact, replicating β-cells double the amount of RNA for most genes, but this upregulation excludes genes involved in β-cell function. These data suggest that the quiescence-proliferation transition involves global amplification of gene expression, except for a subset of tissue-specific genes, which are "left behind" and whose relative mRNA amount decreases. Our work provides a unique resource for the study of replicating β-cells in vivo.

Cell therapy from bench to bedside: Hepatocytes from fibroblasts - the truth and myth of transdifferentiation

Hepatocyte transplantation is an alternative to liver transplantation in certain disorders such as inherited liver diseases and liver failure. It is a relatively less complicated surgical procedure, and has the advantage that it can be repeated several times if unsuccessful. Another advantage is that hepatocytes can be isolated from partly damaged livers which are not suitable for liver transplantation. Despite these advantages hepatocyte transplantation is less popular. Important issues are poor engraftment of the transplanted cells and the scarcity of donor hepatocytes. Generation of "hepatocyte like cells"/iHeps from embryonic stem cells (ES) and induced pluripotent stem cells (iPSCs) by directed differentiation is an emerging solution to the latter issue. Direct conversation or trans-differentiation of fibroblasts to "hepatocyte like cells" is another way which is, being explored. However this method has several inherent and technical disadvantages compared to the directed differentiation from ES or iPSC. There are several methods claiming to be "highly efficient" for generating "highly functional" "hepatocyte like cells". Currently different groups are working independently and coming up with differentiation protocols and each group claiming an advantage for their protocol. Directed differentiation protocols need to be designed, compared, analyzed and tweaked systematically and logically than empirically. There is a need for a well-coordinated global initiative comparable to the Human Genome Project to achieve this goal in the near future.

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.

Large transcription units unify copy number variants and common fragile sites arising under replication stress

Copy number variants (CNVs) resulting from genomic deletions and duplications and common fragile sites (CFSs) seen as breaks on metaphase chromosomes are distinct forms of structural chromosome instability precipitated by replication inhibition. Although they share a common induction mechanism, it is not known how CNVs and CFSs are related or why some genomic loci are much more prone to their occurrence. Here we compare large sets of de novo CNVs and CFSs in several experimental cell systems to each other and to overlapping genomic features. We first show that CNV hotpots and CFSs occurred at the same human loci within a given cultured cell line. Bru-seq nascent RNA sequencing further demonstrated that although genomic regions with low CNV frequencies were enriched in transcribed genes, the CNV hotpots that matched CFSs specifically corresponded to the largest active transcription units in both human and mouse cells. Consistently, active transcription units >1 Mb were robust cell-type-specific predictors of induced CNV hotspots and CFS loci. Unlike most transcribed genes, these very large transcription units replicated late and organized deletion and duplication CNVs into their transcribed and flanking regions, respectively, supporting a role for transcription in replication-dependent lesion formation. These results indicate that active large transcription units drive extreme locus- and cell-type-specific genomic instability under replication stress, resulting in both CNVs and CFSs as different manifestations of perturbed replication dynamics.

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.

Replication of alpha-satellite DNA arrays in endogenous human centromeric regions and in human artificial chromosome

In human chromosomes, centromeric regions comprise megabase-size arrays of 171 bp alpha-satellite DNA monomers. The large distances spanned by these arrays preclude their replication from external sites and imply that the repetitive monomers contain replication origins. However, replication within these arrays has not previously been profiled and the role of alpha-satellite DNA in initiation of DNA replication has not yet been demonstrated. Here, replication of alpha-satellite DNA in endogenous human centromeric regions and in de novo formed Human Artificial Chromosome (HAC) was analyzed. We showed that alpha-satellite monomers could function as origins of DNA replication and that replication of alphoid arrays organized into centrochromatin occurred earlier than those organized into heterochromatin. The distribution of inter-origin distances within centromeric alphoid arrays was comparable to the distribution of inter-origin distances on randomly selected non-centromeric chromosomal regions. Depletion of CENP-B, a kinetochore protein that binds directly to a 17 bp CENP-B box motif common to alpha-satellite DNA, resulted in enrichment of alpha-satellite sequences for proteins of the ORC complex, suggesting that CENP-B may have a role in regulating the replication of centromeric regions. Mapping of replication initiation sites in the HAC revealed that replication preferentially initiated in transcriptionally active regions.

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.

Meta-analysis reveals that genes regulated by the Y chromosome in Drosophila melanogaster are preferentially localized to repressive chromatin

The Drosophila Y chromosome is a degenerated, heterochromatic chromosome with few functional genes. Despite this, natural variation on the Y chromosome in D. melanogaster has substantial trans-acting effects on the regulation of X-linked and autosomal genes. It is not clear, however, whether these genes simply represent a random subset of the genome or whether specific functional properties are associated with susceptibility to regulation by Y-linked variation. Here, we present a meta-analysis of four previously published microarray studies of Y-linked regulatory variation (YRV) in D. melanogaster. We show that YRV genes are far from a random subset of the genome: They are more likely to be in repressive chromatin contexts, be expressed tissue specifically, and vary in expression within and between species than non-YRV genes. Furthermore, YRV genes are more likely to be associated with the nuclear lamina than non-YRV genes and are generally more likely to be close to each other in the nucleus (although not along chromosomes). Taken together, these results suggest that variation on the Y chromosome plays a role in modifying how the genome is distributed across chromatin compartments, either via changes in the distribution of DNA-binding proteins or via changes in the spatial arrangement of the genome in the nucleus.

Back to the origin: reconsidering replication, transcription, epigenetics, and cell cycle control

In bacteria, replication is a carefully orchestrated event that unfolds the same way for each bacterium and each cell division. The process of DNA replication in bacteria optimizes cell growth and coordinates high levels of simultaneous replication and transcription. In metazoans, the organization of replication is more enigmatic. The lack of a specific sequence that defines origins of replication has, until recently, severely limited our ability to define the organizing principles of DNA replication. This question is of particular importance as emerging data suggest that replication stress is an important contributor to inherited genetic damage and the genomic instability in tumors. We consider here the replication program in several different organisms including recent genome-wide analyses of replication origins in humans. We review recent studies on the role of cytosine methylation in replication origins, the role of transcriptional looping and gene gating in DNA replication, and the role of chromatin's 3-dimensional structure in DNA replication. We use these new findings to consider several questions surrounding DNA replication in metazoans: How are origins selected? What is the relationship between replication and transcription? How do checkpoints inhibit origin firing? Why are there early and late firing origins? We then discuss whether oncogenes promote cancer through a role in DNA replication and whether errors in DNA replication are important contributors to the genomic alterations and gene fusion events observed in cancer. We conclude with some important areas for future experimentation.

Asynchronous replication, mono-allelic expression, and long range Cis-effects of ASAR6

Mammalian chromosomes initiate DNA replication at multiple sites along their length during each S phase following a temporal replication program. The majority of genes on homologous chromosomes replicate synchronously. However, mono-allelically expressed genes such as imprinted genes, allelically excluded genes, and genes on female X chromosomes replicate asynchronously. We have identified a cis-acting locus on human chromosome 6 that controls this replication-timing program. This locus encodes a large intergenic non-coding RNA gene named Asynchronous replication and Autosomal RNA on chromosome 6, or ASAR6. Disruption of ASAR6 results in delayed replication, delayed mitotic chromosome condensation, and activation of the previously silent alleles of mono-allelic genes on chromosome 6. The ASAR6 gene resides within an ∼1.2 megabase domain of asynchronously replicating DNA that is coordinated with other random asynchronously replicating loci along chromosome 6. In contrast to other nearby mono-allelic genes, ASAR6 RNA is expressed from the later-replicating allele. ASAR6 RNA is synthesized by RNA Polymerase II, is not polyadenlyated, is restricted to the nucleus, and is subject to random mono-allelic expression. Disruption of ASAR6 leads to the formation of bridged chromosomes, micronuclei, and structural instability of chromosome 6. Finally, ectopic integration of cloned genomic DNA containing ASAR6 causes delayed replication of entire mouse chromosomes.

Mammalian chromosomes are duplicated every cell cycle during a precise temporal DNA replication program. Thus, every chromosome contains regions that are replicated early and other regions that are replicated late during each S phase. Most of the genes, present in two copies on homologous chromosomes, replicate synchronously during each S phase. Exceptions to this rule are genes located on X chromosomes, genetically imprinted genes, and genes subject to allelic exclusion. Thus, all mono-allelically expressed genes are subject to asynchronous replication, where one allele replicates before the other. Perhaps the best-studied example of asynchronous replication in mammals occurs during X inactivation in female cells. A large non-coding RNA gene called XIST, located within the X inactivation center, controls the transcriptional silencing and late replication of the inactive X chromosome. We have identified a locus on human chromosome 6 that shares many characteristics with XIST. This chromosome 6 locus encodes a large intergenic non-coding RNA gene, ASAR6, which displays random mono-allelic expression, asynchronous replication, and controls the mono-allelic expression of other genes on chromosome 6. Our work supports a model in which all mammalian chromosomes contain similar cis-acting loci that function to ensure proper chromosome replication, mitotic condensation, mono-allelic expression, and stability of individual chromosomes.

PcG-mediated higher-order chromatin structures modulate replication programs at the Drosophila BX-C

Polycomb group proteins (PcG) exert conserved epigenetic functions that convey maintenance of repressed transcriptional states, via post-translational histone modifications and high order structure formation. During S-phase, in order to preserve cell identity, in addition to DNA information, PcG-chromatin-mediated epigenetic signatures need to be duplicated requiring a tight coordination between PcG proteins and replication programs. However, the interconnection between replication timing control and PcG functions remains unknown. Using Drosophila embryonic cell lines, we find that, while presence of specific PcG complexes and underlying transcription state are not the sole determinants of cellular replication timing, PcG-mediated higher-order structures appear to dictate the timing of replication and maintenance of the silenced state. Using published datasets we show that PRC1, PRC2, and PhoRC complexes differently correlate with replication timing of their targets. In the fully repressed BX-C, loss of function experiments revealed a synergistic role for PcG proteins in the maintenance of replication programs through the mediation of higher-order structures. Accordingly, replication timing analysis performed on two Drosophila cell lines differing for BX-C gene expression states, PcG distribution, and chromatin domain conformation revealed a cell-type-specific replication program that mirrors lineage-specific BX-C higher-order structures. Our work suggests that PcG complexes, by regulating higher-order chromatin structure at their target sites, contribute to the definition and the maintenance of genomic structural domains where genes showing the same epigenetic state replicate at the same time.

DNA replication is a tightly orchestrated process that precisely duplicates the entire genome during cell division to ensure that daughter cells inherit the same genetic information. The genome is replicated following a specific temporal program, where different segments replicate in distinct moments of the S phase correlating with active (early) and repressed (late) transcriptional state of resident genes. Moreover, replicating chromosomal domains are organized in the nuclear space, perhaps to guarantee the conservation of the same topological order in daughter cells. Epigenetic mechanisms, acting via chromatin organization, determine transcriptional states and must be maintained through cell division. Here, we analyzed in detail the link between Polycomb Group (PcG) proteins, higher-order chromatin structure, and replication timing in Drosophila. By using bioinformatic analyses combined with functional experiments, we show that Polycomb Repressive Complex 1 (PRC1), PRC2, and PhoRC differently correlate with replication timing of their targets and that transcription per se does not determine replication timing. Strikingly, by analyzing the PcG-regulated Bithorax Complex, where PRC1, PRC2, and PhoRC complexes are bound to repressed targets, we provide evidence for a synergistic role of PcG proteins in the modulation and maintenance of replication timing through the definition of specific, topologically distinct genomic domains.

Role of androgen receptor in progression of LNCaP prostate cancer cells from G1 to S phase

Background

The androgen receptor (AR) plays a critical role in the proliferation of prostate cancer cells. However, its mechanism of action in proliferation remains unknown. An understanding of the mechanism of AR action in proliferation may lead to the development of effective strategies for the treatment of prostate cancer.

Methodology/Principal Findings

In this study we report that pulse treatment of synchronized LNCaP cells with Casodex, an AR-antagonist, for 4 hours in mid-G₁ phase was sufficient to prevent cells from entering S phase. Since the assembly of pre-replication complex (pre-RC) in G₁ is required for the progression of cells from G₁ to S phase, the effect of Casodex during mid-G₁ suggested that the role of AR in proliferation might be to regulate the assembly of pre-RC. To test this possibility, we investigated the interaction between AR and Cdc6, an essential component of pre-RC in LNCaP cells. AR co-localized and co-immunoprecipitated with Cdc6, and Casodex treatment disrupted this interaction. AR-immunoprecipitate (AR-IP) also contained cyclin E and cyclin A, which play a critical role in pre-RC assembly and cell cycle entry into S phase, and DNA polymerase-α, PCNA, and ribonucleotide reductase, which are essential for the initiation of DNA synthesis. In addition, in cells in S phase, AR co-sedimented with components of the DNA replication machinery of cells that entered S phase.

Conclusions/Significance

Together, these observations suggest a novel role of AR as a component of the pre-RC to exert control over progression of LNCaP cells from G₁ to S phase through a mechanism that is independent of its role as a transcription factor.

DNA replication timing, genome stability and cancer: late and/or delayed DNA replication timing is associated with increased genomic instability

Normal cellular division requires that the genome be faithfully replicated to ensure that unaltered genomic information is passed from one generation to the next. DNA replication initiates from thousands of origins scattered throughout the genome every cell cycle; however, not all origins initiate replication at the same time. A vast amount of work over the years indicates that different origins along each eukaryotic chromosome are activated in early, middle or late S phase. This temporal control of DNA replication is referred to as the replication-timing program. The replication-timing program represents a very stable epigenetic feature of chromosomes. Recent evidence has indicated that the replication-timing program can influence the spatial distribution of mutagenic events such that certain regions of the genome experience increased spontaneous mutagenesis compared to surrounding regions. This influence has helped shape the genomes of humans and other multicellular organisms and can affect the distribution of mutations in somatic cells. It is also becoming clear that the replication-timing program is deregulated in many disease states, including cancer. Aberrant DNA replication timing is associated with changes in gene expression, changes in epigenetic modifications and an increased frequency of structural rearrangements. Furthermore, certain replication timing changes can directly lead to overt genomic instability and may explain unique mutational signatures that are present in cells that have undergone the recently described processes of "chromothripsis" and "kataegis". In this review, we will discuss how the normal replication timing program, as well as how alterations to this program, can contribute to the evolution of the genomic landscape in normal and cancerous cells.

Megabase replication domains along the human genome: relation to chromatin structure and genome organisation

In higher eukaryotes, the absence of specific sequence motifs, marking the origins of replication has been a serious hindrance to the understanding of (i) the mechanisms that regulate the spatio-temporal replication program, and (ii) the links between origins activation, chromatin structure and transcription. In this chapter, we review the partitioning of the human genome into megabased-size replication domains delineated as N-shaped motifs in the strand compositional asymmetry profiles. They collectively span 28.3% of the genome and are bordered by more than 1,000 putative replication origins. We recapitulate the comparison of this partition of the human genome with high-resolution experimental data that confirms that replication domain borders are likely to be preferential replication initiation zones in the germline. In addition, we highlight the specific distribution of experimental and numerical chromatin marks along replication domains. Domain borders correspond to particular open chromatin regions, possibly encoded in the DNA sequence, and around which replication and transcription are highly coordinated. These regions also present a high evolutionary breakpoint density, suggesting that susceptibility to breakage might be linked to local open chromatin fiber state. Altogether, this chapter presents a compartmentalization of the human genome into replication domains that are landmarks of the human genome organization and are likely to play a key role in genome dynamics during evolution and in pathological situations.

Induced transcription results in local changes in chromatin structure, replication timing, and DNA polytenization in a site of intercalary heterochromatin

In salivary gland polytene chromosomes of Drosophila melanogaster, the regions of intercalary heterochromatin are characterized by late replication, under-replication, and genetic silencing. Using Gal4/UAS system, we induced transcription of sequences adjacent to transgene insertions in the band 11A6-9. This activation resulted in a loss of "silent" and appearance of "active" epigenetic marks, recruitment of RNA polymerase II, and formation of a puff. The activated region is now early replicating and shows increased level of DNA polytenization. Notably, all these changes are restricted to the area around the inserts, whereas the rest of the band remains inactive and late replicating. Although only a short area near the insertion site is transcribed, it results in an "open" chromatin conformation in a much broader region. We conclude that regions of intercalary heterochromatin do not form stand-alone units of late replication and under-replication. Every part of such regions can be activated and polytenized independently of other parts.

Mammalian chromosomes contain cis-acting elements that control replication timing, mitotic condensation, and stability of entire chromosomes

Recent studies indicate that mammalian chromosomes contain discrete cis-acting loci that control replication timing, mitotic condensation, and stability of entire chromosomes. Disruption of the large non-coding RNA gene ASAR6 results in late replication, an under-condensed appearance during mitosis, and structural instability of human chromosome 6. Similarly, disruption of the mouse Xist gene in adult somatic cells results in a late replication and instability phenotype on the X chromosome. ASAR6 shares many characteristics with Xist, including random mono-allelic expression and asynchronous replication timing. Additional "chromosome engineering" studies indicate that certain chromosome rearrangements affecting many different chromosomes display this abnormal replication and instability phenotype. These observations suggest that all mammalian chromosomes contain "inactivation/stability centers" that control proper replication, condensation, and stability of individual chromosomes. Therefore, mammalian chromosomes contain four types of cis-acting elements, origins, telomeres, centromeres, and "inactivation/stability centers", all functioning to ensure proper replication, condensation, segregation, and stability of individual chromosomes.

Replication timing and its emergence from stochastic processes

The temporal organization of DNA replication has puzzled cell biologists since before the mechanism of replication was understood. The realization that replication timing correlates with important features, such as transcription, chromatin structure and genome evolution, and is misregulated in cancer and aging has only deepened the fascination. Many ideas about replication timing have been proposed, but most have been short on mechanistic detail. However, recent work has begun to elucidate basic principles of replication timing. In particular, mathematical modeling of replication kinetics in several systems has shown that the reproducible replication timing patterns seen in population studies can be explained by stochastic origin firing at the single-cell level. This work suggests that replication timing need not be controlled by a hierarchical mechanism that imposes replication timing from a central regulator, but instead results from simple rules that affect individual origins.

A machine learning approach for identifying novel cell type-specific transcriptional regulators of myogenesis

Transcriptional enhancers integrate the contributions of multiple classes of transcription factors (TFs) to orchestrate the myriad spatio-temporal gene expression programs that occur during development. A molecular understanding of enhancers with similar activities requires the identification of both their unique and their shared sequence features. To address this problem, we combined phylogenetic profiling with a DNA-based enhancer sequence classifier that analyzes the TF binding sites (TFBSs) governing the transcription of a co-expressed gene set. We first assembled a small number of enhancers that are active in Drosophila melanogaster muscle founder cells (FCs) and other mesodermal cell types. Using phylogenetic profiling, we increased the number of enhancers by incorporating orthologous but divergent sequences from other Drosophila species. Functional assays revealed that the diverged enhancer orthologs were active in largely similar patterns as their D. melanogaster counterparts, although there was extensive evolutionary shuffling of known TFBSs. We then built and trained a classifier using this enhancer set and identified additional related enhancers based on the presence or absence of known and putative TFBSs. Predicted FC enhancers were over-represented in proximity to known FC genes; and many of the TFBSs learned by the classifier were found to be critical for enhancer activity, including POU homeodomain, Myb, Ets, Forkhead, and T-box motifs. Empirical testing also revealed that the T-box TF encoded by org-1 is a previously uncharacterized regulator of muscle cell identity. Finally, we found extensive diversity in the composition of TFBSs within known FC enhancers, suggesting that motif combinatorics plays an essential role in the cellular specificity exhibited by such enhancers. In summary, machine learning combined with evolutionary sequence analysis is useful for recognizing novel TFBSs and for facilitating the identification of cognate TFs that coordinate cell type-specific developmental gene expression patterns.

Late replication domains in polytene and non-polytene cells of Drosophila melanogaster

In D. melanogaster polytene chromosomes, intercalary heterochromatin (IH) appears as large dense bands scattered in euchromatin and comprises clusters of repressed genes. IH displays distinctly low gene density, indicative of their particular regulation. Genes embedded in IH replicate late in the S phase and become underreplicated. We asked whether localization and organization of these late-replicating domains is conserved in a distinct cell type. Using published comprehensive genome-wide chromatin annotation datasets (modENCODE and others), we compared IH organization in salivary gland cells and in a Kc cell line. We first established the borders of 60 IH regions on a molecular map, these regions containing underreplicated material and encompassing ∼12% of Drosophila genome. We showed that in Kc cells repressed chromatin constituted 97% of the sequences that corresponded to IH bands. This chromatin is depleted for ORC-2 binding and largely replicates late. Differences in replication timing between the cell types analyzed are local and affect only sub-regions but never whole IH bands. As a rule such differentially replicating sub-regions display open chromatin organization, which apparently results from cell-type specific gene expression of underlying genes. We conclude that repressed chromatin organization of IH is generally conserved in polytene and non-polytene cells. Yet, IH domains do not function as transcription- and replication-regulatory units, because differences in transcription and replication between cell types are not domain-wide, rather they are restricted to small “islands” embedded in these domains. IH regions can thus be defined as a special class of domains with low gene density, which have narrow temporal expression patterns, and so displaying relatively conserved organization.

Replication origins and timing of temporal replication in budding yeast: how to solve the conundrum?

Similarly to metazoans, the budding yeast Saccharomyces cereviasiae replicates its genome with a defined timing. In this organism, well-defined, site-specific origins, are efficient and fire in almost every round of DNA replication. However, this strategy is neither conserved in the fission yeast Saccharomyces pombe, nor in Xenopus or Drosophila embryos, nor in higher eukaryotes, in which DNA replication initiates asynchronously throughout S phase at random sites. Temporal and spatial controls can contribute to the timing of replication such as Cdk activity, origin localization, epigenetic status or gene expression. However, a debate is going on to answer the question how individual origins are selected to fire in budding yeast. Two opposing theories were proposed: the "replicon paradigm" or "temporal program" vs. the "stochastic firing". Recent data support the temporal regulation of origin activation, clustering origins into temporal blocks of early and late replication. Contrarily, strong evidences suggest that stochastic processes acting on origins can generate the observed kinetics of replication without requiring a temporal order. In mammalian cells, a spatiotemporal model that accounts for a partially deterministic and partially stochastic order of DNA replication has been proposed. Is this strategy the solution to reconcile the conundrum of having both organized replication timing and stochastic origin firing also for budding yeast? In this review we discuss this possibility in the light of our recent study on the origin activation, suggesting that there might be a stochastic component in the temporal activation of the replication origins, especially under perturbed conditions.

Determinants and dynamics of genome accessibility

In eukaryotes, all DNA-templated reactions occur in the context of chromatin. Nucleosome packaging inherently restricts DNA accessibility for regulatory proteins but also provides an opportunity to regulate DNA-based processes through modulating nucleosome positions and local chromatin structure. Recent advances in genome-scale methods are yielding increasingly detailed profiles of the genomic distribution of nucleosomes, their modifications and their modifiers. The picture now emerging is one in which the dynamic control of genome accessibility is governed by contributions from DNA sequence, ATP-dependent chromatin remodelling and nucleosome modifications. Here we discuss the interplay of these processes by reviewing our current understanding of how chromatin access contributes to the regulation of transcription, replication and repair.

Paucity and preferential suppression of transgenes in late replication domains of the D. melanogaster genome

Background

Eukaryotic genomes are organized in extended domains with distinct features intimately linking genome structure, replication pattern and chromatin state. Recently we identified a set of long late replicating euchromatic regions that are underreplicated in salivary gland polytene chromosomes of D. melanogaster.

Results

Here we demonstrate that these underreplicated regions (URs) have a low density of P-element and piggyBac insertions compared to the genome average or neighboring regions. In contrast, Minos-based transposons show no paucity in URs but have a strong bias to testis-specific genes. We estimated the suppression level in 2,852 stocks carrying a single P-element by analysis of eye color determined by the mini-white marker gene and demonstrate that the proportion of suppressed transgenes in URs is more than three times higher than in the flanking regions or the genomic average. The suppressed transgenes reside in intergenic, genic or promoter regions of the annotated genes. We speculate that the low insertion frequency of P-elements and piggyBacs in URs partially results from suppression of transgenes that potentially could prevent identification of transgenes due to complete suppression of the marker gene. In a similar manner, the proportion of suppressed transgenes is higher in loci replicating late or very late in Kc cells and these loci have a lower density of P-elements and piggyBac insertions. In transgenes with two marker genes suppression of mini-white gene in eye coincides with suppression of yellow gene in bristles.

Conclusions

Our results suggest that the late replication domains have a high inactivation potential apparently linked to the silenced or closed chromatin state in these regions, and that such inactivation potential is largely maintained in different tissues.

Mutational bias shaping fly copy number variation: implications for genome evolution

Copy number variants (CNVs) underlie several genomic disorders and are a major source of genetic innovation. Consequently, any bias affecting their placement in the genome will impact our understanding of human disease and genome evolution. Here we report a mutational bias affecting CNVs that generates different probabilities of duplication and deletion across the genome in association with DNA replication time. We show that this mutational bias has important consequences for genome evolution by leading to different probabilities of gene duplication for different classes of genes and by linking the probability of gene duplication with the transcriptional activity of genes.

Genome-wide dynamics of replication timing revealed by in vitro models of mouse embryogenesis

Differentiation of mouse embryonic stem cells (mESCs) is accompanied by changes in replication timing. To explore the relationship between replication timing and cell fate transitions, we constructed genome-wide replication-timing profiles of 22 independent mouse cell lines representing 10 stages of early mouse development, and transcription profiles for seven of these stages. Replication profiles were cell-type specific, with 45% of the genome exhibiting significant changes at some point during development that were generally coordinated with changes in transcription. Comparison of early and late epiblast cell culture models revealed a set of early-to-late replication switches completed at a stage equivalent to the post-implantation epiblast, prior to germ layer specification and down-regulation of key pluripotency transcription factors [POU5F1 (also known as OCT4)/NANOG/SOX2] and coinciding with the emergence of compact chromatin near the nuclear periphery. These changes were maintained in all subsequent lineages (lineage-independent) and involved a group of irreversibly down-regulated genes, at least some of which were repositioned closer to the nuclear periphery. Importantly, many genomic regions of partially reprogrammed induced pluripotent stem cells (piPSCs) failed to re-establish ESC-specific replication-timing and transcription programs. These regions were enriched for lineage-independent early-to-late changes, which in female cells included the inactive X chromosome. Together, these results constitute a comprehensive "fate map" of replication-timing changes during early mouse development. Moreover, they support a model in which a distinct set of replication domains undergoes a form of "autosomal Lyonization" in the epiblast that is difficult to reprogram and coincides with an epigenetic commitment to differentiation prior to germ layer specification.

Increased rate of human mutations where DNA and RNA polymerases collide

Gene density and orientation of genes in eukaryotes seem to be correlated with the replication origin and the mutation rate is greater in late replicating regions; however, the reason for these patterns is unknown. Here, we investigate predicted replication origins in the human genome and find that levels of polymorphism as well as divergence from the chimpanzee genome are greater in genes transcribed on the lagging strand than those on the leading strand. This might be caused by interference between RNA and DNA polymerases, and avoidance of collisions between these enzymes might be an evolutionary force shaping gene orientation and density surrounding replication start sites. Physical constraints might have a larger influence on genome evolution than previously thought.

Predictable dynamic program of timing of DNA replication in human cells

The organization of mammalian DNA replication is poorly understood. We have produced high-resolution dynamic maps of the timing of replication in human erythroid, mesenchymal, and embryonic stem (ES) cells using TimEX, a method that relies on gaussian convolution of massive, highly redundant determinations of DNA copy-number variations during S phase to produce replication timing profiles. We first obtained timing maps of 3% of the genome using high-density oligonucleotide tiling arrays and then extended the TimEX method genome-wide using massively parallel sequencing. We show that in untransformed human cells, timing of replication is highly regulated and highly synchronous, and that many genomic segments are replicated in temporal transition regions devoid of initiation, where replication forks progress unidirectionally from origins that can be hundreds of kilobases away. Absence of initiation in one transition region is shown at the molecular level by single molecule analysis of replicated DNA (SMARD). Comparison of ES and erythroid cells replication patterns revealed that these cells replicate about 20% of their genome in different quarters of S phase. Importantly, we detected a strong inverse relationship between timing of replication and distance to the closest expressed gene. This relationship can be used to predict tissue-specific timing of replication profiles from expression data and genomic annotations. We also provide evidence that early origins of replication are preferentially located near highly expressed genes, that mid-firing origins are located near moderately expressed genes, and that late-firing origins are located far from genes.

Open chromatin encoded in DNA sequence is the signature of 'master' replication origins in human cells

For years, progress in elucidating the mechanisms underlying replication initiation and its coupling to transcriptional activities and to local chromatin structure has been hampered by the small number (approximately 30) of well-established origins in the human genome and more generally in mammalian genomes. Recent in silico studies of compositional strand asymmetries revealed a high level of organization of human genes around 1000 putative replication origins. Here, by comparing with recently experimentally identified replication origins, we provide further support that these putative origins are active in vivo. We show that regions approximately 300-kb wide surrounding most of these putative replication origins that replicate early in the S phase are hypersensitive to DNase I cleavage, hypomethylated and present a significant enrichment in genomic energy barriers that impair nucleosome formation (nucleosome-free regions). This suggests that these putative replication origins are specified by an open chromatin structure favored by the DNA sequence. We discuss how this distinctive attribute makes these origins, further qualified as 'master' replication origins, priviledged loci for future research to decipher the human spatio-temporal replication program. Finally, we argue that these 'master' origins are likely to play a key role in genome dynamics during evolution and in pathological situations.

SUMOylation enhances DNA methyltransferase 1 activity

DNA methylation regulates gene expression through a complex network of protein-protein and protein-DNA interactions in chromatin. The maintenance methylase, DNMT1 (DNA methyltransferase 1), is a prominent enzyme in the process that is linked to DNA replication and drives the heritable nature of epigenetic modifications. The mechanistic details that explain how DNMT1 catalytic action is directed and regulated in chromatin are important in our overall understanding of gene control. In this work, we show that DNMT1 is modified by SUMOylation and we have mapped these SUMOylation sites by defined mutations. SUMOylated DNMT1 is catalytically active on genomic DNA in vivo and we find that SUMOylation significantly enhances the methylase activity of DNMT1 both in vitro and in chromatin. These data suggest that SUMOylation modulates the endogenous activity of a prominent epigenetic maintenance pathway in somatic cells.

Universal temporal profile of replication origin activation in eukaryotes

Although replication proteins are conserved among eukaryotes, the sequence requirements for replication initiation differ between species. In all species, however, replication origins fire asynchronously throughout S phase. The temporal program of origin firing is reproducible in cell populations but largely probabilistic at the single-cell level. The mechanisms and the significance of this program are unclear. Replication timing has been correlated with gene activity in metazoans but not in yeast. One potential role for a temporal regulation of origin firing is to minimize fluctuations in replication end time and avoid persistence of unreplicated DNA in mitosis. Here, we have extracted the population-averaged temporal profiles of replication initiation rates for S. cerevisiae, S. pombe, D. melanogaster, X. laevis and H. sapiens from genome-wide replication timing and DNA combing data. All the profiles have a strikingly similar shape, increasing during the first half of S phase then decreasing before its end. A previously proposed minimal model of stochastic initiation modulated by accumulation of a recyclable, limiting replication-fork factor and fork-promoted initiation of new origins, quantitatively described the observed profiles without requiring new implementations.

The selective pressure for timely completion of genome replication and optimal usage of replication proteins that must be imported into the cell nucleus can explain the generic shape of the profiles. We have identified a universal behavior of eukaryotic replication initiation that transcends the mechanisms of origin specification. The population-averaged efficiency of replication origin usage changes during S phase in a strikingly similar manner in a highly diverse set of eukaryotes. The quantitative model previously proposed for origin activation in X. laevis can be generalized to explain this evolutionary conservation.

Chromatin state marks cell-type- and gender-specific replication of the Drosophila genome

Duplication of eukaryotic genomes during S phase is coordinated in space and time. In order to identify zones of initiation and cell-type- as well as gender-specific plasticity of DNA replication, we profiled replication timing, histone acetylation, and transcription throughout the Drosophila genome. We observed two waves of replication initiation with many distinct zones firing in early-S phase and multiple, less defined peaks at the end of S phase, suggesting that initiation becomes more promiscuous in late-S phase. A comparison of different cell types revealed widespread plasticity of replication timing on autosomes. Most occur in large regions, but only half coincide with local differences in transcription. In contrast to confined autosomal differences, a global shift in replication timing occurs throughout the single male X chromosome. Unlike in females, the dosage-compensated X chromosome replicates almost exclusively early. This difference occurs at sites that are not transcriptionally hyperactivated, but show increased acetylation of Lys 16 of histone H4 (H4K16ac). This suggests a transcription-independent, yet chromosome-wide process related to chromatin. Importantly, H4K16ac is also enriched at initiation zones as well as early replicating regions on autosomes during S phase. Together, our study reveals novel organizational principles of DNA replication of the Drosophila genome and suggests that H4K16ac is more closely correlated with replication timing than is transcription.

Introduction to molecular combing: genomics, DNA replication, and cancer

The sequencing of the human genome inaugurated a new era in both fundamental and applied genetics. At the same time, the emergence of new technologies for probing the genome has transformed the field of pharmaco-genetics and made personalized genomic profiling and high-throughput screening of new therapeutic agents all but a matter of routine. One of these technologies, molecular combing, has served to bridge the technical gap between the examination of gross chromosomal abnormalities and sequence-specific alterations. Molecular combing provides a new perspective on the structure and dynamics of the human genome at the whole genome and sub-chromosomal levels with a resolution ranging from a few kilobases up to a megabase and more. Originally developed to study genetic rearrangements and to map genes for positional cloning, recent advances have extended the spectrum of its applications to studying the real-time dynamics of the replication of the genome. Understanding how the genome is replicated is essential for elucidating the mechanisms that both maintain genome integrity and result in the instabilities leading to human genetic disease and cancer. In the following, we will examine recent discoveries and advances due to the application of molecular combing to new areas of research in the fields of molecular cytogenetics and cancer genomics.

Epstein-Barr virus episome stability is coupled to a delay in replication timing

The temporal regulation of DNA replication is thought to be important for chromosome organization and genome stability. We show here that Epstein-Barr virus (EBV) genomes replicate in mid- to late S phase and that agents that accelerate replication timing of EBV reduce viral genome stability. Hydroxyurea (HU) treatment, which is known to eliminate EBV episomes, shifted EBV replication to earlier times in the cell cycle. HU treatment correlated with hyperacetylation of histone H3 and loss of telomere repeat factor 2 (TRF2) binding at the EBV origin of plasmid replication (OriP). Deletion of TRF2 binding sites within OriP or short hairpin RNA depletion of TRF2 advanced the replication timing of OriP-containing plasmids. Inhibitors of class I histone deacetylases (HDACs) increased histone acetylation at OriP, advanced the replication timing of EBV, and reduced EBV genome copy number. We also show that HDAC1 and -2 form a stable complex with TRF2 at OriP and that HU treatment inhibits HDAC activity. We propose that the TRF2-HDAC complex enhances EBV episome stability by providing a checkpoint that delays replication initiation at OriP.

ReplicationDomain: a visualization tool and comparative database for genome-wide replication timing data

BACKGROUND

Eukaryotic DNA replication is regulated at the level of large chromosomal domains (0.5-5 megabases in mammals) within which replicons are activated relatively synchronously. These domains replicate in a specific temporal order during S-phase and our genome-wide analyses of replication timing have demonstrated that this temporal order of domain replication is a stable property of specific cell types.

RESULTS

We have developed ReplicationDomain http://www.replicationdomain.org as a web-based database for analysis of genome-wide replication timing maps (replication profiles) from various cell lines and species. This database also provides comparative information of transcriptional expression and is configured to display any genome-wide property (for instance, ChIP-Chip or ChIP-Seq data) via an interactive web interface. Our published microarray data sets are publicly available. Users may graphically display these data sets for a selected genomic region and download the data displayed as text files, or alternatively, download complete genome-wide data sets. Furthermore, we have implemented a user registration system that allows registered users to upload their own data sets. Upon uploading, registered users may choose to: (1) view their data sets privately without sharing; (2) share with other registered users; or (3) make their published or "in press" data sets publicly available, which can fulfill journal and funding agencies' requirements for data sharing.

CONCLUSION

ReplicationDomain is a novel and powerful tool to facilitate the comparative visualization of replication timing in various cell types as well as other genome-wide chromatin features and is considerably faster and more convenient than existing browsers when viewing multi-megabase segments of chromosomes. Furthermore, the data upload function with the option of private viewing or sharing of data sets between registered users should be a valuable resource for the scientific community.

Spatio-temporal dynamics of replication and transcription sites in the mammalian cell nucleus

To study when and where active genes replicated in early S phase are transcribed, a series of pulse-chase experiments are performed to label replicating chromatin domains (RS) in early S phase and subsequently transcription sites (TS) after chase periods of 0 to 24 h. Surprisingly, transcription activity throughout these chase periods did not show significant colocalization with early RS chromatin domains. Application of novel image segmentation and proximity algorithms, however, revealed close proximity of TS with the labeled chromatin domains independent of chase time. In addition, RNA polymerase II was highly proximal and showed significant colocalization with both TS and the chromatin domains. Based on these findings, we propose that chromatin activated for transcription dynamically unfolds or "loops out" of early RS chromatin domains where it can interact with RNA polymerase II and other components of the transcriptional machinery. Our results further suggest that the early RS chromatin domains are transcribing genes throughout the cell cycle and that multiple chromatin domains are organized around the same transcription factory.

Global reorganization of replication domains during embryonic stem cell differentiation

DNA replication in mammals is regulated via the coordinate firing of clusters of replicons that duplicate megabase-sized chromosome segments at specific times during S-phase. Cytogenetic studies show that these “replicon clusters” coalesce as subchromosomal units that persist through multiple cell generations, but the molecular boundaries of such units have remained elusive. Moreover, the extent to which changes in replication timing occur during differentiation and their relationship to transcription changes has not been rigorously investigated. We have constructed high-resolution replication-timing profiles in mouse embryonic stem cells (mESCs) before and after differentiation to neural precursor cells. We demonstrate that chromosomes can be segmented into multimegabase domains of coordinate replication, which we call “replication domains,” separated by transition regions whose replication kinetics are consistent with large originless segments. The molecular boundaries of replication domains are remarkably well conserved between distantly related ESC lines and induced pluripotent stem cells. Unexpectedly, ESC differentiation was accompanied by the consolidation of smaller differentially replicating domains into larger coordinately replicated units whose replication time was more aligned to isochore GC content and the density of LINE-1 transposable elements, but not gene density. Replication-timing changes were coordinated with transcription changes for weak promoters more than strong promoters, and were accompanied by rearrangements in subnuclear position. We conclude that replication profiles are cell-type specific, and changes in these profiles reveal chromosome segments that undergo large changes in organization during differentiation. Moreover, smaller replication domains and a higher density of timing transition regions that interrupt isochore replication timing define a novel characteristic of the pluripotent state.

Microscopy studies have suggested that chromosomal DNA is composed of multiple, megabase-sized segments, each replicated at different times during S-phase of the cell cycle. However, a molecular definition of these coordinately replicated sequences and the stability of the boundaries between them has not been established. We constructed genome-wide replication-timing maps in mouse embryonic stem cells, identifying multimegabase coordinately replicated chromosome segments—“replication domains”—separated by remarkably distinct temporal boundaries. These domain boundaries were shared between several unrelated embryonic stem cell lines, including somatic cells reprogrammed to pluripotency (so-called induced pluripotent stem cells). However, upon differentiation to neural precursor cells, domains encompassing approximately 20% of the genome changed their replication timing, temporally consolidating into fewer, larger replication domains that were conserved between different neural precursor cell lines. Domains that changed replication timing showed a unique sequence composition, a strongly biased directionality for changes in resident gene expression, and altered radial positioning within the three-dimensional space in the cell nucleus, suggesting that changes in replication timing are related to the reorganization of higher-order chromosome structure and function during differentiation. Moreover, the property of smaller discordantly replicating domains may define a novel characteristic of pluripotency.

Analyzing the temporal order of DNA replication across the genome during embryonic stem cell differentiation reveals stable boundaries between coordinately replicated regions that consolidate into fewer, larger domains during differentiation.

The DNA end-binding protein Ku regulates silencing at the internal HML and HMR loci in Saccharomyces cerevisiae

Heterochromatin resides near yeast telomeres and at the cryptic mating-type loci, HML and HMR, where it silences transcription of the alpha- and a-mating-type genes, respectively. Ku is a conserved DNA end-binding protein that binds telomeres and regulates silencing in yeast. The role of Ku in silencing is thought to be limited to telomeric silencing. Here, we tested whether Ku contributes to silencing at HML or HMR. Mutant analysis revealed that yKu70 and Sir1 act collectively to silence the mating-type genes at HML and HMR. In addition, loss of yKu70 function leads to expression of different reporter genes inserted at HMR. Quantitative chromatin-immunoprecipitation experiments revealed that yKu70 binds to HML and HMR and that binding of Ku to these internal loci is dependent on Sir4. The interaction between yKu70 and Sir4 was characterized further and found to be dependent on Sir2 but not on Sir1, Sir3, or yKu80. These observations reveal that, in addition to its ability to bind telomeric DNA ends and aid in the silencing of genes at telomeres, Ku binds to internal silent loci via protein-protein interactions and contributes to the efficient silencing of these loci.

Replication-timing-correlated spatial chromatin arrangements in cancer and in primate interphase nuclei

Using published high-resolution data on S-phase replication timing, we determined the three-dimensional (3D) nuclear arrangement of 33 very-early-replicating and 31 very-late-replicating loci. We analyzed diploid human, non-human primate and rearranged tumor cells by 3D fluorescence in situ hybridization with the aim of investigating the impact of chromosomal structural changes on the nuclear organization of these loci. Overall, their topology was found to be largely conserved between cell types, species and in tumor cells. Early-replicating loci were localized in the nuclear interior, whereas late-replicating loci showed a broader distribution with a higher preference for the periphery than for late-BrdU-incorporation foci. However, differences in the spatial arrangement of early and late loci of chromosome 2, as compared with those from chromosome 5, 7 and 17, argue against replication timing as a major driving force for the 3D radial genome organization in human lymphoblastoid cell nuclei. Instead, genomic properties, and local gene density in particular, were identified as the decisive parameters. Further detailed comparisons of chromosome 7 loci in primate and tumor cells suggest that the inversions analyzed influence nuclear topology to a greater extent than the translocations, thus pointing to geometrical constraints in the 3D conformation of a chromosome territory.

Dynamics and interplay of nuclear architecture, genome organization, and gene expression

The organization of the genome in the nucleus of a eukaryotic cell is fairly complex and dynamic. Various features of the nuclear architecture, including compartmentalization of molecular machines and the spatial arrangement of genomic sequences, help to carry out and regulate nuclear processes, such as DNA replication, DNA repair, gene transcription, RNA processing, and mRNA transport. Compartmentalized multiprotein complexes undergo extensive modifications or exchange of protein subunits, allowing for an exquisite dynamics of structural components and functional processes of the nucleus. The architecture of the interphase nucleus is linked to the spatial arrangement of genes and gene clusters, the structure of chromatin, and the accessibility of regulatory DNA elements. In this review, we discuss recent studies that have provided exciting insight into the interplay between nuclear architecture, genome organization, and gene expression.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Transcription factor Runx1 recruits the polyomavirus replication origin to replication factories

Eukaryotic DNA replication takes place in the replication factories, where replication proteins are properly assembled to form replication forks. Thus, recruitment of DNA replication origins to the replication factories must be the key step for the regulation of DNA replication. The transcription factor Runx1 associates with the nuclear matrix, the putative substructure of DNA replication factories. An earlier report from our laboratory showed that Runx1 activates polyomavirus DNA replication, and that this requires its nuclear matrix-binding activity. Here, we show that Runx1 activates polyomavirus DNA replication by stimulating the binding of the viral-encoded replication initiator/helicase, large T antigen, to its replication origin. We found that newly replicated polyomavirus DNA is associated with the nuclear matrix and that large T antigen is targeted to replication factories, suggesting that polyomavirus is replicated in replication factories on the nuclear matrix. Although Runx1 did not co-localize with large T antigen-containing foci by itself, it co-localized with large T antigen-containing replication factories during Runx1-dependent polyomavirus DNA replication. These observations together suggest that Runx1 recruits the polyomavirus replication origin to the replication factory on the nuclear matrix, and that this requires the nuclear matrix-binding activity of Runx1.

ATM-mediated transcriptional and developmental responses to gamma-rays in Arabidopsis

ATM (Ataxia Telangiectasia Mutated) is an essential checkpoint kinase that signals DNA double-strand breaks in eukaryotes. Its depletion causes meiotic and somatic defects in Arabidopsis and progressive motor impairment accompanied by several cell deficiencies in patients with ataxia telangiectasia (AT). To obtain a comprehensive view of the ATM pathway in plants, we performed a time-course analysis of seedling responses by combining confocal laser scanning microscopy studies of root development and genome-wide expression profiling of wild-type (WT) and homozygous ATM-deficient mutants challenged with a dose of γ-rays (IR) that is sublethal for WT plants. Early morphologic defects in meristematic stem cells indicated that AtATM, an Arabidopsis homolog of the human ATM gene, is essential for maintaining the quiescent center and controlling the differentiation of initial cells after exposure to IR. Results of several microarray experiments performed with whole seedlings and roots up to 5 h post-IR were compiled in a single table, which was used to import gene information and extract gene sets. Sequence and function homology searches; import of spatio-temporal, cell cycling, and mutant-constitutive expression characteristics; and a simplified functional classification system were used to identify novel genes in all functional classes. The hundreds of radiomodulated genes identified were not a random collection, but belonged to functional pathways such as those of the cell cycle; cell death and repair; DNA replication, repair, and recombination; and transcription; translation; and signaling, indicating the strong cell reprogramming and double-strand break abrogation functions of ATM checkpoints. Accordingly, genes in all functional classes were either down or up-regulated concomitantly with downregulation of chromatin deacetylases or upregulation of acetylases and methylases, respectively. Determining the early transcriptional indicators of prolonged S-G2 phases that coincided with cell proliferation delay, or an anticipated subsequent auxin increase, accelerated cell differentiation or death, was used to link IR-regulated hallmark functions and tissue phenotypes after IR. The transcription burst was almost exclusively AtATM-dependent or weakly AtATR-dependent, and followed two major trends of expression in atm: (i)-loss or severe attenuation and delay, and (ii)-inverse and/or stochastic, as well as specific, enabling one to distinguish IR/ATM pathway constituents. Our data provide a large resource for studies on the interaction between plant checkpoints of the cell cycle, development, hormone response, and DNA repair functions, because IR-induced transcriptional changes partially overlap with the response to environmental stress. Putative connections of ATM to stem cell maintenance pathways after IR are also discussed.